• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    擇伐對思茅松自然種群結(jié)構(gòu)和空間分布格局的影響

    2015-02-13 01:17:13黃小波李帥鋒蘇建榮
    生態(tài)學(xué)報 2015年24期
    關(guān)鍵詞:思茅松齡級天然林

    黃小波, 李帥鋒, 蘇建榮,*, 蘇 磊

    1中國林業(yè)科學(xué)研究院資源昆蟲研究所, 昆明 650224 2國家林業(yè)局普洱森林生態(tài)系統(tǒng)定位研究站, 昆明 650224 3中國科學(xué)院植物研究所植被與環(huán)境變化國家重點(diǎn)實(shí)驗(yàn)室, 北京 100093

    擇伐對思茅松自然種群結(jié)構(gòu)和空間分布格局的影響

    黃小波1,2, 李帥鋒1,2, 蘇建榮1,2,*, 蘇 磊3

    1中國林業(yè)科學(xué)研究院資源昆蟲研究所, 昆明 650224 2國家林業(yè)局普洱森林生態(tài)系統(tǒng)定位研究站, 昆明 650224 3中國科學(xué)院植物研究所植被與環(huán)境變化國家重點(diǎn)實(shí)驗(yàn)室, 北京 100093

    思茅松是云南省重要的材脂兼用樹種,而擇伐是思茅松天然林管理的有效手段,也是獲取木材的重要途徑?;趯υ颇鲜【肮瓤h思茅松天然林的樣地調(diào)查數(shù)據(jù),采用種群徑級結(jié)構(gòu)代替年齡結(jié)構(gòu)、Ripley的L函數(shù)點(diǎn)格局等方法,對思茅松天然林在不同強(qiáng)度擇伐10a后的種群結(jié)構(gòu)和空間分布格局及其不同發(fā)育階段在不同尺度下的相互關(guān)系進(jìn)行了系統(tǒng)研究和分析。結(jié)果表明:擇伐10a后3種類型的思茅松自然種群的胸徑、樹高之間呈顯著的正相關(guān)且可用冪指數(shù)模擬預(yù)測模型。未擇伐的思茅松種群呈衰退型,輕度擇伐(12.9%)后種群呈穩(wěn)定型,中度擇伐(29.3%)后的種群屬增長型,適當(dāng)?shù)闹卸葥穹ビ欣谧匀环N群的更新;未擇伐與中度擇伐后的思茅松種群分布格局整體呈隨機(jī)型,但后者較前者更趨向于均勻分布,輕度擇伐后的思茅松種群分布格局隨著尺度的增加由隨機(jī)轉(zhuǎn)向聚集分布。擇伐對不同徑級思茅松天然林種群空間分布格局有重要影響,未擇伐的思茅松自然種群隨著發(fā)育階段的發(fā)展,多數(shù)尺度上由隨機(jī)分布逐漸轉(zhuǎn)變?yōu)榫奂植?,至大樹則又以隨機(jī)分布為主;輕度擇伐和中度擇伐的思茅松自然種群則隨著發(fā)育階段的發(fā)展而由完全的隨機(jī)分布逐漸轉(zhuǎn)變?yōu)榫奂植?,中度擇伐的聚集分布則更為明顯。思茅松不同發(fā)育階段間在所有尺度上基本都呈顯著正相關(guān),表明思茅松種群種內(nèi)競爭弱,有利于種群的維持。

    思茅松;胸徑-樹高模型;種群結(jié)構(gòu);擇伐強(qiáng)度;點(diǎn)格局分析;空間關(guān)聯(lián)性

    種群結(jié)構(gòu)和空間格局分析是了解種群的生物學(xué)特性、環(huán)境因子相互作用的基本規(guī)律及其與空間分布相聯(lián)系的生態(tài)過程的基礎(chǔ)[1],對揭示群落結(jié)構(gòu)的形成、動態(tài)變化和維持機(jī)制有重要意義[2]。種群年齡結(jié)構(gòu)是種群內(nèi)不同年齡個體數(shù)量的分布情況,既反映了不同齡級個體的組配情況,也反映了種群的動態(tài)變化和發(fā)展趨勢[3]。種群空間分布格局是森林群落結(jié)構(gòu)的重要特征和生態(tài)學(xué)中研究的熱點(diǎn)[4]。植物種群的分布格局與尺度有很大的關(guān)系[5-6];此外還受種群分布區(qū)生境異質(zhì)性和種子擴(kuò)散限制等的影響[7]。

    Ripley提出的點(diǎn)格局分析[8]因克服了傳統(tǒng)格局分析方法上種群的分布類型和空間尺度的關(guān)系都得不到全面反映的缺點(diǎn)而被廣泛應(yīng)用。目前,國外在這方面的研究多集中在熱帶雨林[9-11]。20世紀(jì)90年代張金屯教授將點(diǎn)格局分析方法引入我國,近年圍繞甘肅大野口青海云杉種群[12]、古田山常綠闊葉林優(yōu)勢樹種[13]、小興安嶺闊葉紅松林[14]、長白山闊葉紅松林[15]、滇西北云南紅豆杉種群[16]等物種的分布格局及空間關(guān)聯(lián)性的研究,很好地探究了多尺度上植物種群的生態(tài)學(xué)特性及過程。擇伐是森林中常見的干擾因素,擇伐對種群分布格局影響的研究日益受到重視[17-19],有學(xué)者認(rèn)為適度的擇伐可加快林分生長[20],改變林分空間分布格局[21];胡云云等[18]和沈林等[22]的研究指出,輕度擇伐更加有利于林分空間結(jié)構(gòu)的優(yōu)化,而陳輝榮等[23]的研究則表明,輕度和中度擇伐后林分空間結(jié)構(gòu)趨于優(yōu)化。

    思茅松(Pinuskesiyavar.langbianensis)是以三針一束為主的松樹,主要分布于云南省西南部哀牢山西坡以西的亞熱帶南部,在西藏東南部低海拔地區(qū)也有零星分布[24]。是我國亞熱帶西南部山地的代表種[25]。思茅松具有生長迅速、材質(zhì)優(yōu)良、松脂產(chǎn)量高等特點(diǎn),是云南省重要的材脂兼用樹種和主要造林樹種,據(jù)云南省2007年森林資源連續(xù)清查結(jié)果顯示,思茅松林面積占云南省總林地面積的3.71%[24],具有重要的經(jīng)濟(jì)價值、森林生態(tài)服務(wù)功能和碳匯效益[26-27]。目前,針對思茅松的研究,主要涉及思茅松天然林群落生態(tài)學(xué)[28]、光合生理生態(tài)[29]、化學(xué)成分分析[30]、種質(zhì)資源[31]等方面,還未見關(guān)于思茅松天然林種群空間分布格局的報道。目前,擇伐是思茅松林區(qū)的主要采伐方式,本文利用點(diǎn)格局法分析了自然狀態(tài)下的思茅松種群以及不同強(qiáng)度擇伐10年后思茅松種群結(jié)構(gòu)、空間分布格局及其空間關(guān)聯(lián)性,找出其內(nèi)在的空間分布規(guī)律及產(chǎn)生的原因,以期為以后的思茅松天然林保護(hù)與合理利用提供理論依據(jù)。

    1 研究地概況

    研究地選擇云南省西南部景谷縣境內(nèi),景谷縣位于滇中高原南緣,瀾滄江以東,北回歸線橫穿縣境中部,地理位置為N22°49′—N23°52′,E100°02′—101°07′,海拔600—2920m。境內(nèi)地形分割破碎起伏較大,山高坡陡,為典型的中心切割窄谷地形。地帶性氣候?qū)倌蟻啛釒Ъ撅L(fēng)氣候,四季不分明,干濕季分明,日溫差較大,當(dāng)?shù)啬昶骄鶜鉁?0.1℃,≥10℃年活動積溫7360.9℃。極端最高氣溫38.6℃,極端最低氣溫為0.2℃,年均降水量1026—1402mm,年均蒸發(fā)量1601—2022mm,無霜期353d。地帶性土壤以紅棕壤為主,呈明顯的垂直分布。

    該區(qū)域的自然植被主要為思茅松天然林,同時季風(fēng)常綠闊葉林也有少量分布。思茅松天然林的群落結(jié)構(gòu)分喬木、灌木、草本3層。喬木層主要分兩層,喬木上層主要以思茅松為主,其它物種數(shù)量較少,明顯區(qū)別于喬木下層;喬木下層主要組成物種有高山栲(Castanopsisdelavayi)、華南石櫟(Lithocarpusfenestratus)、小果栲(C.fleuryi)、紅木荷(Schimawallichii)、黃毛青岡(Cyclobalanopsisdelavayi)、滇青岡(C.glaucoides)、毛葉青岡(C.kerrii)和槲櫟(Quercusaliena)等。灌木層常見種類有母豬果(Helicianilagirica)、臭莢蒾(Viburnumfoetidum)、余甘子(Phyllanthusemblica)、米飯花(Lyoniaovalifolia)、大葉算盤子(Glochidionlanceolarium)、鹽膚木(Rhuschinensis)、野漆(Toxicodendronsuccedaneum)、黑面神(Breyniafruticosa)、多花野牡丹(Melastomapolyanthum)、艾膠算盤子(G.lanceolarium)、傘形紫金牛(Ardisiacorymbifera)、大葉千斤拔(Flemingiamacrophylla)等。草本層常見種類有三點(diǎn)金(Desmodiumtriflorum)、長尖莎草(Cyperuscuspidatus)、斑茅(Saccharumarundinaceum)、無芒竹葉草(Oplismenuscompositus)、紫莖澤蘭(Eupatoriumadenphorum)和擬鱗毛蕨(Kuniwatsukiacuspidata)等物種。

    2 研究方法

    2.1 樣地設(shè)置與調(diào)查

    2003年6月在全面勘查的基礎(chǔ)上,在景谷縣思茅松自然種群分布的集中區(qū)域,選擇立地條件和群落結(jié)構(gòu)基本一致的地段,設(shè)置60m×45m的樣地3塊。樣地1為未擇伐標(biāo)準(zhǔn)樣地,作為對照,樣地2伐前總株數(shù)為85株,擇伐株數(shù)為11株,擇伐強(qiáng)度為12.9%,樣地3伐前總株數(shù)為92株,擇伐株數(shù)為27株,擇伐強(qiáng)度為29.3%。參照鄭麗鳳等[32]的研究,將樣地2和樣地3的擇伐強(qiáng)度定義為輕度和中度。出于短期經(jīng)濟(jì)效益的考慮,當(dāng)?shù)亓謭霾扇 翱炒罅粜。?yōu)留劣”的方式進(jìn)行擇伐,伐后天然更新。2013年6月,采用“相鄰格子法”進(jìn)行每木調(diào)查,以5m×5m的小樣方為基本單元,記錄樣方內(nèi)全部思茅松的胸徑、樹高及空間坐標(biāo),并測定樣地的經(jīng)緯度、海拔等立地因子。表1為擇伐10a后樣地基本情況介紹。

    表1 擇伐10a后研究樣地基本概況Table 1 General status of study sites in different selective cutting after ten years

    2.2 胸徑-樹高模型

    在林業(yè)調(diào)查中,胸徑、樹高是重要的測量因子,??捎脕碛嬎銟涓刹姆e、立地質(zhì)量評價和在生長與收獲模型中的應(yīng)用等[33],本文利用樹高、胸徑調(diào)查數(shù)據(jù)建立思茅松的胸徑-樹高預(yù)測模型,模型建立使用SPSS軟件17.0軟件。

    2.3 種群結(jié)構(gòu)劃分

    在相同的環(huán)境條件下,同一樹種的徑級和齡級對環(huán)境的反應(yīng)規(guī)律具有一致性[34],故本文采用徑級代替齡級的方法進(jìn)行年齡結(jié)構(gòu)的分析。根據(jù)云南思茅松生活史,按高度(H)和胸徑(DBH)大小,將云南思茅松劃分為5個不同生長階段,即Ⅰ:幼苗,H<1.3m;Ⅱ:小樹,2.5cm≤DBH<7.5cm;Ⅲ:中樹,7.5cm≤DBH<22.5cm;Ⅳ:大樹,22.5cm≤DBH<47.5cm;Ⅴ:老樹,47.5cm≤DBH,在此基礎(chǔ)上進(jìn)行空間分布格局分析。

    思茅松的徑級結(jié)構(gòu)通過徑級大小分為11級,將DBH<2.5cm作為Ⅰ級,以5cm作為級距進(jìn)行等級劃分,47.5cm≤DBH<52.5cm為Ⅺ級。

    思茅松的高度級結(jié)構(gòu)分為7級,H<1.3m為Ⅰ級,Ⅱ級1.3—5m,5m以上每隔5m為一個高度級,≥25m為Ⅶ級。統(tǒng)計以上各徑級和高度級的個體數(shù),繪制種群結(jié)構(gòu)圖。

    2.4 點(diǎn)格局分析

    本文采用Ripley的L函數(shù)進(jìn)行空間點(diǎn)格局分析。Ripley的L函數(shù)是由Ripley的K函數(shù)發(fā)展而來。Ripley的K函數(shù)是以植物個體坐標(biāo)為基礎(chǔ)分析種群空間格局的工具,該分析方法可分析不同尺度下種群分布格局和各齡級個體的空間關(guān)聯(lián),是目前廣泛使用的分析種群空間分布格局的重要方法[35]。用Monte-Carlo擬合檢驗(yàn)計算上下包跡線,即置信區(qū)間,并以空間尺度r為橫坐標(biāo),上下包跡線為縱坐標(biāo)作圖。單變量分析中,種群實(shí)際分布數(shù)據(jù)計算所得的不同尺度下的函數(shù)值在包跡線內(nèi),則符合隨機(jī)分布;若在包跡線上,則為集群分布;若在包跡線下,則為均勻分布。在雙變量分析中,種群實(shí)際分布數(shù)據(jù)計算得到的不同尺度下的函數(shù)值高于上包跡線說明兩類格局為顯著正相關(guān),在上下包跡線之間為無關(guān)聯(lián),低于下包跡線則為顯著負(fù)相關(guān)。

    點(diǎn)格局和不同齡級間空間關(guān)聯(lián)分析通過生態(tài)學(xué)軟件ADE-4完成,采用的空間尺度為樣地最短邊長的一半即22.5m,步長取1.0m,Monte-Carlo隨機(jī)模擬的數(shù)目為1000次,得到由上下包跡線圍成的99%置信區(qū)間[16]。

    3 結(jié)果與分析

    3.1 思茅松胸徑-樹高模型

    未擇伐的思茅松自然種群和不同強(qiáng)度擇伐10a后思茅松胸徑-樹高模型如圖1所示,從圖中可以看出,未擇伐和不同擇伐強(qiáng)度的思茅松胸徑樹高均有顯著關(guān)系,呈冪函數(shù)關(guān)系。分析胸徑-樹高關(guān)系可為今后的思茅松胸徑分布的預(yù)測提供理論依據(jù)。

    圖1 思茅松胸徑-樹高關(guān)系Fig.1 Relationships between DBH and tree height of Pinus kesiya var.langbianensis population

    3.2 種群結(jié)構(gòu)

    思茅松不同類型種群的徑級結(jié)構(gòu)如圖2所示。不同類型的思茅松種群徑級分布存在明顯差異。未擇伐的思茅松自然種群徑級結(jié)構(gòu)的最高峰出現(xiàn)在Ⅶ級(大樹)(圖2),另一高峰出現(xiàn)在Ⅰ級(幼苗),Ⅶ—Ⅺ(大樹及老樹)的個體數(shù)為41株,占整個種群的53.2%,即使林下具有一定的幼苗更新,但整個未擇伐思茅松種群呈衰退趨勢。輕度擇伐的思茅松自然種群的徑級結(jié)構(gòu)呈紡錘型,說明輕度擇伐后思茅松種群屬于穩(wěn)定型結(jié)構(gòu)。中度擇伐的思茅松種群Ⅰ級幼苗約占整個種群的85.3%,雖然缺失Ⅲ、Ⅹ、Ⅺ級,但種群呈倒“J”型(圖2),中度擇伐種群屬增長型種群,林下存在著豐富的幼苗和幼樹,種群更新良好,但缺乏中等徑級的個體。

    圖2 思茅松種群的徑級結(jié)構(gòu)Fig.2 DBH class structure of Pinus kesiya var.langbianensis populationⅠ:DBH<2.5cm;Ⅱ:2.5cm≤DBH<7.5cm;Ⅲ:7.5cm≤DBH<12.5cm;Ⅳ:12.5cm≤DBH<17.5cm;Ⅴ:17.5cm≤DBH<22.5cm;Ⅵ:22.5cm≤DBH<27.5cm;Ⅶ:27.5cm≤DBH<32.5cm;Ⅷ:32.5cm≤DBH<37.5cm;Ⅸ:37.5cm≤DBH<42.5cm;Ⅹ:42.5cm≤DBH<47.5cm;Ⅺ:47.5cm≤DBH<52.5cm

    圖3 思茅松種群高度級結(jié)構(gòu)Fig.3 Height class structure of Pinus kesiya var.langbianensis populationⅠ:H<1.3m;Ⅱ:1.3m≤H<5m;Ⅲ:5m≤H<10m;Ⅳ:10m≤H<15m;Ⅴ:15m≤DBH<20m;Ⅵ:20m≤H<25m;Ⅶ:25m≤H

    思茅松種群的個體數(shù)隨高度級的分布如圖3所示,高度級結(jié)構(gòu)可作為齡級估測很好的一個參數(shù),未擇伐思茅松種群個體數(shù)大部分集中在Ⅵ和Ⅶ級,總數(shù)為48株,約占總個體的62.34%,幼苗和較小的思茅松個體數(shù)量嚴(yán)重不足,影響思茅松種群的更新。輕度擇伐后,思茅松個體數(shù)多集中在中等高度級上,呈現(xiàn)穩(wěn)定的發(fā)展態(tài)勢。中度擇伐后,Ⅰ級個體最多,為591株,占總數(shù)的85.3%,雖有高度級的缺失,但種群更新能力強(qiáng),有利于種群發(fā)展。

    3.3 思茅松自然種群的空間分布格局

    3個思茅松種群個體分布如圖4所示,它直觀地表現(xiàn)出了云南思茅松在樣地內(nèi)的空間分布狀態(tài)。未擇伐、輕度擇伐和中度擇伐的思茅松個體分別有77株、123株、693株。其中未擇伐的思茅松種群中有老樹1株,大樹45株,中樹14株,小樹5株,幼苗12株。輕度擇伐強(qiáng)度的思茅松種群有大樹29株,中樹61株,小樹24株,幼苗9株。中度擇伐強(qiáng)度的思茅松種群中有大樹54株,中樹14株,小樹34株,幼苗591株??傮w來看,隨著擇伐強(qiáng)度的增加,思茅松的種群數(shù)量越來越高,且中度擇伐影響的最高。

    圖4 不同類型思茅松種群的個體分布位點(diǎn)圖Fig.4 Mapped points pattern of different Pinus kesiya var.langbianensis population

    3個不同類型的思茅松種群的點(diǎn)格局分析結(jié)果見圖5,未擇伐和中度擇伐的思茅松種群的全部個體在所研究的空間尺度(0—22.5m)范圍內(nèi),分布格局呈隨機(jī)分布,且中度擇伐的思茅松自然種群在所有尺度上比未擇伐更傾向于均勻分布。而輕度擇伐種群在0—16m尺度上呈隨機(jī)分布,隨后轉(zhuǎn)向集群分布,且聚集度有緩慢增長趨勢(圖5)。

    圖5 思茅松種群的空間分布格局 Fig.5 Spatial patterns of Pinus kesiya var.langbianensis populationsL(r)表示Ripley′s K函數(shù)值;圖中實(shí)線為L(r)值,表示實(shí)際數(shù)據(jù)計算所得的值;虛線為包跡線,表示所模擬的99%置信區(qū)間

    3.4 思茅松自然種群各生長階段空間分布

    樹種在不同發(fā)育階段會表現(xiàn)不同的空間分布格局,這與森林群落的自然稀疏過程、干擾格局以及環(huán)境的變化有著密切的關(guān)系[2]。不同類型的思茅松種群不同發(fā)育階段個體的空間分布格局及其與尺度的關(guān)系如圖6所示,未擇伐的思茅松自然種群大樹和小樹在小尺度上呈集群分布,隨著尺度的增加則呈現(xiàn)隨機(jī)分布,不同的是中樹在13—22.5m的尺度上又表現(xiàn)為集群分布,小樹在0—5m的尺度上呈隨機(jī)分布,在>5m的尺度上呈集群分布,幼苗在0—20m的尺度呈隨機(jī)分布,在>20m尺度上表現(xiàn)為集群分布;輕度擇伐和中度擇伐的思茅松自然種群在大部分尺度上隨著發(fā)育階段的增加由隨機(jī)分布向聚集分布轉(zhuǎn)換,輕度擇伐的思茅松自然種群的幼苗在所有尺度上呈隨機(jī)分布,小樹在小尺度上(0—2m)呈聚集分布,2—15m尺度上呈隨機(jī)分布,到15—22.5m尺度上呈集群分布,中樹在0—19m尺度上呈隨機(jī)分布,>19m尺度上呈集群分布,大樹在0—3m尺度上分布格局呈隨機(jī)分布,>3m尺度上呈集群分布;中度擇伐的思茅松自然種群中幼苗在不同尺度下均呈隨機(jī)分布,小樹在0—1m的小尺度上表現(xiàn)為集群分布,在2—22.5m的尺度上呈隨機(jī)分布。

    圖6 思茅松種群不同生長階段的空間分布格局Fig.6 Spatial patterns of Pinus kesiya var.langbianensis population among different growth stagesA: 未擇伐幼苗;B: 未擇伐小樹;C: 未擇伐中樹;D: 未擇伐大樹;1: 輕度擇伐;2: 中度擇伐

    3.5 思茅松種群各生長階段空間關(guān)聯(lián)

    3種思茅松類型種群不同徑級間的關(guān)聯(lián)性總體情況見表2。未擇伐的思茅松種群除中樹與幼樹之間在0—8m尺度上無關(guān)聯(lián)或相互獨(dú)立外,其它生長階段的思茅松相互之間在所有尺度上均表現(xiàn)為顯著正相關(guān);與此類似,輕度擇伐的思茅松種群,除中樹與小樹之間在2—4m尺度上無關(guān)聯(lián)外,其它生長階段的思茅松相互之間在所有尺度上也表現(xiàn)為顯著正相關(guān),而中度擇伐的思茅松種群除大樹和中樹在10—13m尺度上無關(guān)聯(lián),其它生長階段的思茅松相互之間在所有尺度上亦表現(xiàn)為顯著正相關(guān)。3個種群不同生長階段之間均無負(fù)相關(guān)空間關(guān)聯(lián)。

    4 結(jié)論和討論

    4.1 擇伐對思茅松自然種群結(jié)構(gòu)與空間分布的影響

    經(jīng)過擇伐10a后的更新,思茅松的齡級結(jié)構(gòu)由未擇伐的衰退型,隨著擇伐強(qiáng)度的增加變成穩(wěn)定型,繼而成為增長型的齡級結(jié)構(gòu)。長期以來,景谷縣作為云南省思茅松天然林分布最為集中的區(qū)域之一[31],通過皆伐與擇伐獲取木材是思茅松森林經(jīng)營管理的重要模式,而擇伐可以明顯改變林分的空間結(jié)構(gòu)[22]。思茅松是當(dāng)?shù)丶撅L(fēng)常綠闊葉林正向演替的先鋒樹種[28],擇伐后思茅松天然林形成大小不同的“林窗”斑塊,利于思茅松種子得以成功定居并能更新成幼苗,由于其生長迅速,大量的思茅松中小徑級個體可以較快的占領(lǐng)喬木層中,而林窗越大,則其個體越多,形成充足的幼苗更新層。而未擇伐思茅松齡級結(jié)構(gòu)呈衰退型的主要特征為幼齡期個體缺乏,原因可能有兩個:一是思茅松幼苗生長不耐蔭涼,在林下存在更新障礙[28],不適應(yīng)在林冠層郁閉度較高的林內(nèi),未擇伐的思茅松天然林喬木下層多是季風(fēng)常綠闊葉林的優(yōu)勢物種,森林郁閉度較大,發(fā)芽率較低;二是林下灌木和草本植物豐富,種間競爭激烈,導(dǎo)致思茅松幼苗的死亡率較高。林下更新能力差是思茅松天然林的普遍現(xiàn)象,使得未擇伐思茅松天然林年齡結(jié)構(gòu)呈現(xiàn)衰退型。

    表2 思茅松種群各生長階段空間關(guān)聯(lián)Table 2 Spatial associations of Pinus kesiya var.langbianensis population among different growth stages

    +: 顯著正相關(guān);0: 無相關(guān);BIG: big;ME: middle;SM: small;SE: seedling

    本研究發(fā)現(xiàn),中度擇伐和未擇伐的思茅松種群的空間分布格局在總體上都呈隨機(jī)分布,而輕度擇伐思茅松種群則隨著尺度的增加由隨機(jī)分布轉(zhuǎn)向聚集分布。研究表明,種群的空間分布格局通常會受多種因素的影響,如種子擴(kuò)散限制、環(huán)境異質(zhì)性、種內(nèi)種間競爭等[16,36]。未擇伐思茅松種群總體呈隨機(jī)分布原因與彭少麟等[37]的研究結(jié)果較為一致,即衰退型年齡結(jié)構(gòu)的空間分布呈隨機(jī)分布。中度擇伐的思茅松種群空間分布呈隨機(jī)分布的原因與其繁殖特性和環(huán)境因子有重要關(guān)系。思茅松的種子較小,數(shù)量充沛,且種子具長翅,傳播距離較遠(yuǎn),可隨機(jī)分布在中度擇伐后形成的較大“林窗”內(nèi);思茅松天然林林木自然稀疏開始早,且較強(qiáng)烈,也是其呈隨機(jī)分布的原因之一[38],同時,中度擇伐相對于皆伐等屬于中等強(qiáng)度的人為干擾,中等強(qiáng)度的人為干擾的林分易呈隨機(jī)分布[22]。而中度擇伐較未擇伐思茅松自然種群更傾向于均勻分布,研究表明思茅松易在林內(nèi)大的林窗更新,且可迅速形成單一優(yōu)勢林分[31],其種子的繁殖特性是其偏向均勻分布的主要原因。而輕度擇伐思茅松種群在<16m尺度上呈隨機(jī)分布的主要原因可能與銀杉分布格局的成因相似[39],即是中小尺度的環(huán)境條件相似、種子散布隨機(jī)所致,而在較大尺度上呈集群分布,主要是由于輕度擇伐后產(chǎn)生極其有限的“林窗”,呈集群分布可抵御外部競爭,利用有限的空間資源。

    擇伐強(qiáng)度使不同發(fā)育階段的思茅松種群的空間分布格局出現(xiàn)顯著差異,這與擇伐改變了思茅松種群的齡級結(jié)構(gòu)有重要關(guān)系。其中未擇伐過的思茅松幼苗在0—20m范圍內(nèi)呈隨機(jī)分布,隨著尺度的增加,趨向于聚集分布;對于小樹和中樹,都是先在一定的尺度上呈隨機(jī)分布,然后隨著尺度的增加呈聚集分布;而大樹在小尺度上呈聚集分布,在大尺度上屬于隨機(jī)分布。究其原因,可能是由于較小齡級的個體所需的資源較少,遇到的競爭相對較弱,所以最終表現(xiàn)為聚集分布。隨著齡級的增大,比如大樹,種群個體遇到的生存壓力導(dǎo)致種內(nèi)競爭加劇,這種密度制約因素最終導(dǎo)致自疏效應(yīng),致使個體減少,表現(xiàn)為最后的隨機(jī)分布[7]。輕度擇伐和中度擇伐后幼苗階段在所有尺度上表現(xiàn)為隨機(jī)分布,主要因?yàn)槠渥陨淼姆敝程匦砸约吧L階段的自疏和它疏效應(yīng)所致,輕度擇伐后,隨著齡級的增長,分布格局有隨機(jī)分布向聚集分布變化的趨勢,原因主要是對于競爭有限的資源,各齡級思茅松個體競爭能力的強(qiáng)弱不同所致;中度擇伐后,大樹和中樹呈明顯的聚集分布,且中樹的聚集程度大于大樹,這是由于中度擇伐的作業(yè)方式以采大留小為原則,且強(qiáng)度較大,故呈聚集分布,這反映了森林管理者的利用方式對種群空間分布干擾的影響。

    4.2 擇伐對思茅松自然種群間空間關(guān)聯(lián)的影響

    未擇伐的思茅松種群和不同強(qiáng)度擇伐10a后的思茅松種群各齡級間的基本上都呈顯著的正關(guān)聯(lián)。通常來說,同一物種不同齡級在空間上的相互關(guān)聯(lián)性是對一定時期內(nèi)植物種群組成個體之間相互關(guān)系的靜態(tài)描述,這種關(guān)系不僅包括空間分布關(guān)系,也包含著個體間的功能關(guān)系,反映出種群的現(xiàn)狀和種群內(nèi)的相互關(guān)系[40]。本研究得出的正關(guān)聯(lián)的結(jié)果是由同一樹種的個體生物學(xué)特性的一致性決定的,也就是說對環(huán)境的要求和適應(yīng)性是一致的。體現(xiàn)了植物利用資源的相似性和生態(tài)位的重疊[7],也說明不同齡級個體在空間上是交錯分布,這樣有利于對各種資源的充分利用,對整個種群的生存和發(fā)展是有利的。

    4.3 思茅松天然林森林經(jīng)營對策與建議

    自20世紀(jì)50年代以來,思茅松林就是森林采伐的重要對象,加之近年采脂產(chǎn)業(yè)需求量大、桉樹和一些經(jīng)濟(jì)作物的推廣種植,使得思茅松天然林面積數(shù)量銳減;而20世紀(jì)80年代以來發(fā)展的思茅松人工林存在生產(chǎn)率低、受病蟲害威脅、林地地力衰退等問題。為充分利用現(xiàn)有思茅松天然林,促進(jìn)思茅松天然林的可持續(xù)發(fā)展,根據(jù)本研究提出如下建議:

    首先,未擇伐的思茅松天然林種群呈衰退型,思茅松林下存在幼苗更新障礙[28],擇伐后有效的打破了該障礙,使種群結(jié)構(gòu)呈穩(wěn)定型和增長型,因此,在思茅松天然林經(jīng)營管理過程中,應(yīng)遵循種群結(jié)構(gòu)和空間格局分布規(guī)律,通過適度的擇伐,如擇伐老樹,清理林下活地被物等撫育措施,創(chuàng)造林窗和林隙,增加林內(nèi)光照強(qiáng)度,降低幼苗所承受的競爭壓力,促進(jìn)幼苗的生長,促進(jìn)自然更新。其次,應(yīng)改善和促進(jìn)擇伐方式,在思茅松聚集分布的區(qū)域應(yīng)是重點(diǎn)采伐的地方,可以多采用群團(tuán)狀擇伐,形成較大林窗,促進(jìn)林內(nèi)小徑級思茅松生長。最后,在擇伐強(qiáng)度上,可選擇中等強(qiáng)度擇伐,促進(jìn)種群更新強(qiáng)度,在林下儲備更多的思茅松幼樹,思茅松林在水平結(jié)構(gòu)和垂直結(jié)構(gòu)上達(dá)到合理配置以保持種群的持續(xù)穩(wěn)定發(fā)展。

    [1] 侯向陽, 韓進(jìn)軒.長白山紅松林主要樹種空間格局的模擬分析.植物生態(tài)學(xué)報, 1997, 21(3): 242-249.

    [2] Greig-Smith P.Quantitative Plant Ecology.California: University of California Press, 1983.

    [3] 蘇志堯, 吳大榮, 陳光北.粵北天然林優(yōu)勢種群結(jié)構(gòu)與空間格局動態(tài).應(yīng)用生態(tài)學(xué)報, 2000, 11(3): 337-341.

    [4] Getzin S, Wiegand T, Wiegand K, He F.Heterogeneity influences spatial patterns and demographics in forest stands.Journal of Ecology, 2008, 96(4): 807-820.

    [5] 張金屯.植物種群空間分布的點(diǎn)格局分析.植物生態(tài)學(xué)報, 1998, 22(4): 344-349.

    [6] Fangliang H, Legendre P, LaFrankie J V.Distribution patterns of tree species in a Malaysian tropical rain forest.Journal of Vegetation Science, 1997, 8(1): 105-114.

    [7] 張俊艷, 成克武, 臧潤國.海南島熱帶自然針葉林主要樹種的空間格局及關(guān)聯(lián)性.生物多樣性, 2014, 22(2): 129-140.

    [8] Ripley B D.Spatial Statistics.New York: John Wiley &Sons, 2005.

    [9] Condit R, Ashton P S, Baker P, Bunyavejchewin S, Gunatilleke S, Gunatilleke N, Hubbell S P, Forester R B, Itoh A, LaFrankie J V, Lee H S, Losos E, Manokaran N, Sukumar R, Yamakura T.Spatial patterns in the distribution of tropical tree species.Science, 2000, 288(5470): 1414-1418.

    [10] Engelbrecht B M J, Comita L S, Condit R, Kursar T A, Tyree M T, Turner B L, Hubbell S P.Drought sensitivity shapes species distribution patterns in tropical forests.Nature, 2007, 447(7140): 80-82.

    [11] Nieder J, Engwald S, Klawun M, Barthlott W.Spatial distribution of vascular epiphytes (including Hemiepiphytes) in a lowland amazonian rain forest (Surumoni Crane Plot) of Southern Venezuela1.Biotropica, 2000, 32(3): 385-396.

    [12] 趙麗瓊, 黃華國, 梁大雙, 張曉麗.甘肅大野口青海云杉種群的空間分布格局.北京林業(yè)大學(xué)學(xué)報, 2010, 32(4): 59-64.

    [13] 李立, 陳建華, 任海保, 米湘成, 于明堅, 楊波.古田山常綠闊葉林優(yōu)勢樹種甜櫧和木荷的空間格局分析.植物生態(tài)學(xué)報, 2010, 34(3): 241-252.

    [14] 侯紅亞, 王立海.小興安嶺闊葉紅松林物種組成及主要種群的空間分布格局.應(yīng)用生態(tài)學(xué)報, 2013, 24(11): 3043-3049.

    [15] 王蕾, 張春雨, 趙秀海.長白山闊葉紅松林的空間分布格局.林業(yè)科學(xué), 2009, 45(5): 54-59.

    [16] 李帥鋒, 劉萬德, 蘇建榮, 郎學(xué)東, 張志鈞.云南蘭坪云南紅豆杉種群年齡結(jié)構(gòu)與空間分布格局分析.西北植物學(xué)報, 2013, 33(4): 792-799.

    [17] 周蔚, 楊華, 亢新剛, 沈林.擇伐強(qiáng)度對長白山區(qū)天然云冷杉針闊混交林空間結(jié)構(gòu)的影響.西北林學(xué)院學(xué)報, 2012, 27(4): 7-12.

    [18] 胡云云, 閔志強(qiáng), 高延, 馮啟祥.擇伐對天然云冷杉林林分生長和結(jié)構(gòu)的影響.林業(yè)科學(xué), 2011, 47(2): 15-24.

    [19] 周建云, 李榮, 張文輝, 何景峰.不同間伐強(qiáng)度下遼東櫟種群結(jié)構(gòu)特征與空間分布格局.林業(yè)科學(xué), 2012, 48(4): 149-155.

    [20] M?kinen H, Isom?ki A.Thinning intensity and long-term changes in increment and stem form of Scots pine trees.Forest Ecology and Management, 2004, 203(1-3): 21-34.

    [21] 董希斌.森林擇伐對林分的影響.東北林業(yè)大學(xué)學(xué)報, 2002, 30(5): 15-18.

    [22] 沈林, 楊華, 亢新剛, 岳剛.擇伐強(qiáng)度對天然云冷杉林空間分布格局的影響.中南林業(yè)科技大學(xué)學(xué)報, 2013, 33(1): 68-74.

    [23] 陳輝榮, 周新年, 蔡瑞添, 胡喜生, 巫志龍, 周成軍.天然林不同強(qiáng)度擇伐后林分空間結(jié)構(gòu)變化動態(tài).植物科學(xué)學(xué)報, 2012, 30(3): 230-237.

    [24] 李帥鋒, 蘇建榮, 劉萬德, 郎學(xué)東, 張志鈞, 蘇磊, 賈呈鑫卓, 楊華景.思茅松天然群體種實(shí)表型變異.植物生態(tài)學(xué)報, 2013, 37(11): 998-1009.

    [25] 李福秀, 張曉明, 周體林, 張德國.思茅松針葉束水培生根的研究.西南林學(xué)院學(xué)報, 2005, 25(1): 24-27.

    [26] 溫慶忠, 趙元藩, 陳曉鳴, 楊子祥, 艾建林, 楊曉松.中國思茅松林生態(tài)服務(wù)功能價值動態(tài)研究.林業(yè)科學(xué)研究, 2010, 23(5): 671-677.

    [27] 李江, 孟夢, 邱瓊, 朱宏濤, 翟明普, 陳宏偉, 郭永清, 馮弦, 劉永剛.思茅松中幼齡人工林生物量及生產(chǎn)力動態(tài).東北林業(yè)大學(xué)學(xué)報, 2010, 38(8): 36-38.

    [28] 李帥鋒, 劉萬德, 蘇建榮, 張志鈞, 劉慶云.普洱季風(fēng)常綠闊葉林次生演替中木本植物幼苗更新特征.生態(tài)學(xué)報, 2012, 32(18): 5653-5662.

    [29] 劉娟, 梁軍生, 王健敏, 陳曉鳴, 楊子祥, 陳航.思茅松干季光合生理日動態(tài)及光響應(yīng)特征分析.林業(yè)科學(xué)研究, 2009, 22(5): 677-682.

    [30] 鄭華, 張弘, 李文良, 趙虹, 張汝國, 于連松.景東地區(qū)云南松和思茅松揮發(fā)物組分含量的比較.福建農(nóng)林大學(xué)學(xué)報: 自然科學(xué)版, 2010, 39(4): 366-369.

    [31] 李帥鋒, 蘇建榮, 劉萬德, 郎學(xué)東, 張志鈞, 蘇磊, 楊華景.云南省思茅松林群落數(shù)量分類及物種多樣性與自然環(huán)境的關(guān)系.生態(tài)學(xué)雜志, 2013, 32(12): 3152-3159.

    [32] 鄭麗鳳, 周新年.擇伐強(qiáng)度對中亞熱帶天然針闊混交林林分空間結(jié)構(gòu)的影響.武漢植物學(xué)研究, 2010, 27(5): 515-521.

    [33] Dorado F C, Anta M B, Parresol B R,lvarez González J G.A stochastic height-diameter model for maritime pine ecoregions in Galicia (northwestern Spain).Annals of Forest Science, 2005, 62(5): 455-465.

    [34] Frost I, Rydin H.Spatial pattern and size distribution of the animal-dispersed treeQuercusroburin two spruce-dominated forests.Ecoscience, 2000, 7(1): 38-44.

    [35] Gray L, He F.Spatial point-pattern analysis for detecting density-dependent competition in a boreal chronosequence of Alberta.Forest Ecology and Management, 2009, 259(1): 98-106.

    [36] 張興旺, 張小平, 郭傳友, 張強(qiáng).皖北石灰?guī)r山地青檀種群不同發(fā)育階段的點(diǎn)格局分析.生態(tài)學(xué)雜志, 2013, 32(3): 542-550.

    [37] 彭少麟, 王伯蓀.鼎湖山森林群落分析——Ⅲ 種群分布格局.熱帶亞熱帶森林生態(tài)系統(tǒng)研究, 1984, 2(1): 24-37.

    [38] 黃三祥, 張赟, 趙秀海.山西太岳山油松種群的空間分布格局.福建林學(xué)院學(xué)報, 2009, 29(3): 269-273.

    [39] 謝宗強(qiáng), 陳偉烈, 劉正宇, 江明喜, 黃漢東.銀杉種群的空間分布格局.植物學(xué)報, 1999, 41(1): 95-101.

    [40] 王磊, 孫啟武, 郝朝運(yùn), 田勝尼, 張姍姍, 陳一錕, 張小平.皖南山區(qū)南方紅豆杉種群不同齡級立木的點(diǎn)格局分析.應(yīng)用生態(tài)學(xué)報, 2010, 21(2): 272-278.

    Influences of selective cutting on the structure and spatial distribution patterns of a natural population ofPinuskesiyavar.langbianensis

    HUANG Xiaobo1,2,LI Shuaifeng1,2,SU Jianrong1,2,*, SU Lei3

    1ResearchInstituteofResourceInsects,ChineseAcademyofForestry,Kunming650224,China2ThePu′erForestEcosystemResearchStation,StateForestryAdministration,Kunming650224,China3StateKeyLaboratoryofVegetationandEnvironmentChange,InstituteofBotany,ChineseAcademyofSciences,Beijing100093,China

    Pinuskesiyavar.langbianensisis an important source of timber and resin production in Yunnan Province, China.In forest management, selective cutting is an effective means ofP.kesiyaforest management, and is also the primary way to obtain timber.This study is based on data from sampling plots of aP.kesiyanatural forest in Jinggu County, Yunnan Province.The method of substituting space for time using the L function (based on the derivatives of Ripley′s K function) was used to analyze the influence of 10years of different selective cutting intensities on the structure and spatial distribution patterns of aP.kesiyanatural forest population.This study also investigated the relationship betweenP.kesiya′s development stages and spatial scales.The results showed a significantly positive correlation between DBH and tree height after 10years of selective cutting, and are used as a power model to simulate a prediction model.TheP.kesiyapopulation without selective cutting showed a decling population type.The population under low intensity selective cutting (12.9%) showed a stable type, and the diameter-class structure of the population under medium intensity selective cutting (29.3%) presented an invert “J” shape, was an expanded type, and had more individuals with the small diameter class.This indicates that medium intensity selective cutting is useful for promoting regeneration.BothP.kesiyapopulations without selective cutting and with selective cutting under medium intensity were randomly distributed, but the latter tended to be more regularly distributed than the former.The spatial pattern of low intensity cutting changed from random to aggregated with the increase of scale.Selective cutting intensities changed the spatial distribution patterns of different DBH classes, and also influenced the overall distribution patterns ofP.kesiya.With the progression of developmental stages, the distribution of trees under no selective cutting changed from random to aggregate.For large trees, the distribution pattern was most often random;the distribution patterns under low and medium intensity cutting changed from completely random to aggregated with the development, and an aggregated distribution at medium intensity selective cutting is more apparent than low intensity selective cutting.The different growth stages ofP.kesiyapopulations were significantly positively correlated at all scales.Intraspecific competition was weak, which had the advantage of population maintenance.Low and medium intensity selective cutting may enhance the stability of the forest community spatial structure.Medium selective cutting is more effective at forest structure maintenance than low selective cutting.Our results have important theoretical and practical implications for establishing a development plan for the sustainable management ofP.kesiyaforests.

    Pinuskesiyavar.langbianensis;DBH-Height equation;population structure;selective cutting intensity;point pattern analysis;spatial association

    云南省科技計劃項(xiàng)目(2013RA004);中國林科院中央級公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(riricaf2012001Z);林業(yè)公益性行業(yè)科研專項(xiàng)項(xiàng)目(201404211)

    2014-06-05; < class="emphasis_bold">網(wǎng)絡(luò)出版日期:

    日期:2015-05-20

    10.5846/stxb201406051159

    *通訊作者Corresponding author.E-mail: jianrongsu@vip.sina.com

    黃小波, 李帥鋒, 蘇建榮, 蘇磊.擇伐對思茅松自然種群結(jié)構(gòu)和空間分布格局的影響.生態(tài)學(xué)報,2015,35(24):8241-8250.

    Huang X B,Li S F,Su J R, Su L.Influences of selective cutting on the structure and spatial distribution patterns of a natural population ofPinuskesiyavar.langbianensis.Acta Ecologica Sinica,2015,35(24):8241-8250.

    猜你喜歡
    思茅松齡級天然林
    堵河源自然保護(hù)區(qū)四照花種群結(jié)構(gòu)及動態(tài)分析
    氣象因子和坡向?qū)λ济┧僧a(chǎn)脂量的影響1)
    山西省天然林保護(hù)工程二期建設(shè)成效與經(jīng)驗(yàn)探討
    氣象因子和坡向?qū)λ济┧僧a(chǎn)脂量的影響1)
    天津:全面停止天然林商業(yè)性采伐
    太行山南麓不同齡級荊條光合特性及光響應(yīng)研究
    天保工程:把天然林都保護(hù)起來
    綠色中國(2019年19期)2019-11-26 07:13:16
    天然林保護(hù)任重而道遠(yuǎn)
    綠色中國(2019年19期)2019-11-26 07:13:16
    山西龐泉溝自然保護(hù)區(qū)遼東櫟種群生命表研究
    不同基質(zhì)配比對思茅松多穴盤幼苗生長的影響
    午夜免费男女啪啪视频观看 | 尤物成人国产欧美一区二区三区| 深爱激情五月婷婷| 亚洲性夜色夜夜综合| 亚洲欧美精品自产自拍| 你懂的网址亚洲精品在线观看 | 久久精品影院6| 久久亚洲精品不卡| 亚洲美女黄片视频| 精品人妻偷拍中文字幕| 久久精品91蜜桃| 色综合色国产| 亚洲无线观看免费| 日韩一区二区视频免费看| 国产探花在线观看一区二区| 十八禁国产超污无遮挡网站| 丝袜喷水一区| 国产综合懂色| 成年女人毛片免费观看观看9| 秋霞在线观看毛片| 欧美xxxx黑人xx丫x性爽| 三级男女做爰猛烈吃奶摸视频| 国产黄色视频一区二区在线观看 | 黄色欧美视频在线观看| 中出人妻视频一区二区| 日本撒尿小便嘘嘘汇集6| 亚洲成a人片在线一区二区| 99久国产av精品国产电影| 97热精品久久久久久| 一区二区三区免费毛片| 亚洲国产精品成人综合色| 亚洲欧美清纯卡通| 一边摸一边抽搐一进一小说| 卡戴珊不雅视频在线播放| 最近中文字幕高清免费大全6| 寂寞人妻少妇视频99o| 女同久久另类99精品国产91| 免费看日本二区| 亚洲精品国产成人久久av| 精品一区二区三区视频在线| 国产三级在线视频| 日本-黄色视频高清免费观看| h日本视频在线播放| 亚洲精华国产精华液的使用体验 | 男插女下体视频免费在线播放| 精品免费久久久久久久清纯| 久久精品综合一区二区三区| 岛国在线免费视频观看| 性色avwww在线观看| 久久久久国内视频| 欧美zozozo另类| 美女内射精品一级片tv| 亚洲高清免费不卡视频| 毛片女人毛片| 一区二区三区免费毛片| 综合色丁香网| 国产黄a三级三级三级人| 色在线成人网| 国产精品综合久久久久久久免费| av中文乱码字幕在线| 国产精品av视频在线免费观看| 午夜福利高清视频| av天堂在线播放| 波多野结衣高清无吗| 久久久国产成人精品二区| 五月玫瑰六月丁香| 搡老岳熟女国产| 国产三级中文精品| 国产成人福利小说| 日日撸夜夜添| 搡老岳熟女国产| av福利片在线观看| 97超碰精品成人国产| 综合色av麻豆| 1024手机看黄色片| 我要搜黄色片| 国产私拍福利视频在线观看| 免费在线观看影片大全网站| 午夜福利高清视频| 国产精品久久久久久亚洲av鲁大| 国产av在哪里看| 国产真实伦视频高清在线观看| 国产精品免费一区二区三区在线| 亚洲美女视频黄频| 最近视频中文字幕2019在线8| 成人性生交大片免费视频hd| 一本精品99久久精品77| 色哟哟哟哟哟哟| 免费av毛片视频| 日韩av在线大香蕉| 九九爱精品视频在线观看| 亚洲色图av天堂| 日本色播在线视频| 超碰av人人做人人爽久久| 欧美中文日本在线观看视频| 精品少妇黑人巨大在线播放 | 国产精品人妻久久久影院| 日韩欧美精品免费久久| 国产综合懂色| 精品一区二区三区人妻视频| 亚洲精品粉嫩美女一区| 亚洲婷婷狠狠爱综合网| 99久久成人亚洲精品观看| 少妇熟女aⅴ在线视频| 人人妻人人澡欧美一区二区| 91久久精品国产一区二区三区| 男人和女人高潮做爰伦理| 亚洲天堂国产精品一区在线| 精品一区二区三区视频在线观看免费| 五月玫瑰六月丁香| 乱系列少妇在线播放| 男女那种视频在线观看| 免费在线观看影片大全网站| 91久久精品国产一区二区三区| 国产精品美女特级片免费视频播放器| 青春草视频在线免费观看| 国产午夜福利久久久久久| 久久久久久久久大av| 男女视频在线观看网站免费| 午夜久久久久精精品| 亚洲欧美精品综合久久99| 伦精品一区二区三区| 亚洲真实伦在线观看| 噜噜噜噜噜久久久久久91| 国产美女午夜福利| 日日摸夜夜添夜夜添av毛片| 国产亚洲精品综合一区在线观看| 欧美zozozo另类| 少妇人妻一区二区三区视频| 日日啪夜夜撸| h日本视频在线播放| 亚洲电影在线观看av| 免费在线观看影片大全网站| 亚洲成人精品中文字幕电影| 热99在线观看视频| 变态另类丝袜制服| 变态另类丝袜制服| 免费看美女性在线毛片视频| 18禁在线无遮挡免费观看视频 | 国产人妻一区二区三区在| 美女高潮的动态| 国产精品伦人一区二区| 亚洲av中文字字幕乱码综合| 国产免费一级a男人的天堂| 最新中文字幕久久久久| 国产一区亚洲一区在线观看| 亚洲国产色片| 国内精品久久久久精免费| 免费av不卡在线播放| 丰满的人妻完整版| 欧美性感艳星| 观看免费一级毛片| 不卡视频在线观看欧美| 少妇被粗大猛烈的视频| 69人妻影院| 午夜激情欧美在线| 搡女人真爽免费视频火全软件 | 久久这里只有精品中国| 国产 一区精品| 91午夜精品亚洲一区二区三区| 日韩人妻高清精品专区| 精华霜和精华液先用哪个| 一区福利在线观看| 国产国拍精品亚洲av在线观看| 成人特级av手机在线观看| 一个人免费在线观看电影| 毛片女人毛片| 99久国产av精品| 国产又黄又爽又无遮挡在线| 日韩欧美在线乱码| 国产白丝娇喘喷水9色精品| 国产女主播在线喷水免费视频网站 | 欧美国产日韩亚洲一区| 91久久精品电影网| a级毛片a级免费在线| 亚洲欧美精品综合久久99| 伦理电影大哥的女人| 九九久久精品国产亚洲av麻豆| 1024手机看黄色片| 赤兔流量卡办理| 神马国产精品三级电影在线观看| 岛国在线免费视频观看| 亚洲国产精品国产精品| 看片在线看免费视频| 亚洲在线自拍视频| 国产午夜精品论理片| 日韩,欧美,国产一区二区三区 | 美女 人体艺术 gogo| 国产精品久久久久久久久免| 亚洲国产精品国产精品| 欧美+日韩+精品| 不卡一级毛片| 国产精品国产三级国产av玫瑰| 亚洲五月天丁香| 日本欧美国产在线视频| 老司机福利观看| 91在线观看av| 男人的好看免费观看在线视频| 国产午夜精品久久久久久一区二区三区 | 久久国产乱子免费精品| 国产av一区在线观看免费| 久久精品夜色国产| 国模一区二区三区四区视频| av女优亚洲男人天堂| 国产国拍精品亚洲av在线观看| 九九久久精品国产亚洲av麻豆| 最新中文字幕久久久久| 可以在线观看毛片的网站| 老司机午夜福利在线观看视频| 久久久久久久久久黄片| 免费看美女性在线毛片视频| 1024手机看黄色片| 精品久久久久久成人av| 日韩欧美国产在线观看| 内射极品少妇av片p| 麻豆久久精品国产亚洲av| 日本色播在线视频| 国产国拍精品亚洲av在线观看| 毛片一级片免费看久久久久| 少妇熟女欧美另类| 观看免费一级毛片| 国产精品乱码一区二三区的特点| 色哟哟哟哟哟哟| www日本黄色视频网| 欧美一区二区国产精品久久精品| videossex国产| 欧美丝袜亚洲另类| 国产精品久久电影中文字幕| 亚洲欧美成人精品一区二区| 久久欧美精品欧美久久欧美| 十八禁网站免费在线| 欧美性感艳星| 国产高清视频在线播放一区| 91久久精品电影网| 蜜臀久久99精品久久宅男| 在线观看一区二区三区| 日韩欧美三级三区| 岛国在线免费视频观看| 亚洲欧美中文字幕日韩二区| 九九爱精品视频在线观看| 91狼人影院| 国产视频一区二区在线看| 亚洲国产精品sss在线观看| 国产一区二区激情短视频| 美女内射精品一级片tv| 久久99热6这里只有精品| av在线老鸭窝| 国语自产精品视频在线第100页| 国产成人精品久久久久久| 日本 av在线| 俄罗斯特黄特色一大片| 国产欧美日韩一区二区精品| 老女人水多毛片| 亚洲最大成人中文| 2021天堂中文幕一二区在线观| 欧美性猛交黑人性爽| 在线免费十八禁| 丝袜喷水一区| 亚洲经典国产精华液单| 亚洲成人精品中文字幕电影| 99在线人妻在线中文字幕| 深爱激情五月婷婷| 久久九九热精品免费| 亚洲第一电影网av| 波野结衣二区三区在线| 婷婷色综合大香蕉| 国产精品爽爽va在线观看网站| 特大巨黑吊av在线直播| 日韩一本色道免费dvd| 六月丁香七月| 真人做人爱边吃奶动态| 97碰自拍视频| 国产精品嫩草影院av在线观看| 亚洲在线自拍视频| 精品一区二区三区视频在线观看免费| 一区福利在线观看| 亚洲熟妇中文字幕五十中出| 少妇猛男粗大的猛烈进出视频 | 美女被艹到高潮喷水动态| 熟女电影av网| 久久99热这里只有精品18| 久久久国产成人免费| 欧美+日韩+精品| 女生性感内裤真人,穿戴方法视频| 一级av片app| 欧美三级亚洲精品| 久久久久久国产a免费观看| 欧美极品一区二区三区四区| 国产精品久久久久久亚洲av鲁大| 亚洲四区av| 不卡视频在线观看欧美| 欧美一区二区国产精品久久精品| 欧美国产日韩亚洲一区| 亚洲成人中文字幕在线播放| 久久精品夜色国产| 三级国产精品欧美在线观看| 欧美激情在线99| 成人欧美大片| 波多野结衣高清作品| 亚洲av电影不卡..在线观看| 亚洲精品在线观看二区| 丝袜美腿在线中文| av在线亚洲专区| 波多野结衣高清作品| 欧美高清性xxxxhd video| 国内少妇人妻偷人精品xxx网站| 五月玫瑰六月丁香| 国产在视频线在精品| 久久精品人妻少妇| 禁无遮挡网站| 色5月婷婷丁香| 亚洲内射少妇av| 免费黄网站久久成人精品| 精品一区二区三区视频在线观看免费| 国产女主播在线喷水免费视频网站 | 最好的美女福利视频网| 亚洲精品色激情综合| 午夜老司机福利剧场| 在线播放国产精品三级| 亚洲国产日韩欧美精品在线观看| 最好的美女福利视频网| 国产av不卡久久| 精品一区二区三区视频在线| 一级毛片aaaaaa免费看小| 日本黄大片高清| 91久久精品国产一区二区三区| 久久国内精品自在自线图片| 久久鲁丝午夜福利片| 亚洲av免费高清在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产成人精品久久久久久| 欧美激情在线99| 午夜日韩欧美国产| 长腿黑丝高跟| 国产成人a区在线观看| 久久精品久久久久久噜噜老黄 | 日本免费a在线| 国产精品人妻久久久久久| 精品一区二区免费观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品亚洲一区二区| 亚洲国产色片| 日韩一本色道免费dvd| 99国产精品一区二区蜜桃av| 免费看光身美女| 有码 亚洲区| 一级毛片aaaaaa免费看小| 亚洲av.av天堂| 精品久久久久久成人av| 少妇被粗大猛烈的视频| 午夜激情欧美在线| 欧美又色又爽又黄视频| 99在线视频只有这里精品首页| 在线天堂最新版资源| 我要搜黄色片| 日本熟妇午夜| 中文字幕熟女人妻在线| 国产精品不卡视频一区二区| 精品久久国产蜜桃| 久久久欧美国产精品| 国产成人91sexporn| 亚洲av不卡在线观看| 搡老岳熟女国产| 少妇人妻精品综合一区二区 | 少妇熟女aⅴ在线视频| 日韩欧美三级三区| 国产 一区 欧美 日韩| eeuss影院久久| 国产精品一二三区在线看| 在线观看av片永久免费下载| 久久午夜亚洲精品久久| 黄色一级大片看看| 我要搜黄色片| 色5月婷婷丁香| 免费看美女性在线毛片视频| 97超碰精品成人国产| 国产精品精品国产色婷婷| 欧美最黄视频在线播放免费| 国产成人福利小说| 高清毛片免费看| 一级黄色大片毛片| 在线观看美女被高潮喷水网站| 欧洲精品卡2卡3卡4卡5卡区| 赤兔流量卡办理| 亚洲精品日韩av片在线观看| 国产一区二区激情短视频| 国产精品一区www在线观看| 91在线观看av| 欧美日韩精品成人综合77777| 精品免费久久久久久久清纯| av黄色大香蕉| ponron亚洲| 日韩成人av中文字幕在线观看 | 日本免费一区二区三区高清不卡| 久久久久久久午夜电影| 亚洲精品456在线播放app| 亚洲成人av在线免费| 久久精品综合一区二区三区| 国产精品伦人一区二区| 精品不卡国产一区二区三区| 我要看日韩黄色一级片| 国产精品一区二区免费欧美| 亚洲精品一区av在线观看| 日韩亚洲欧美综合| 亚洲中文字幕日韩| 麻豆国产av国片精品| 欧美三级亚洲精品| 精品一区二区三区视频在线| 免费在线观看影片大全网站| 欧美成人一区二区免费高清观看| 亚洲经典国产精华液单| 国产一区二区三区在线臀色熟女| 欧美潮喷喷水| 一级毛片电影观看 | 精品一区二区免费观看| 又爽又黄a免费视频| 成人国产麻豆网| 高清毛片免费看| 午夜a级毛片| 特大巨黑吊av在线直播| 搡女人真爽免费视频火全软件 | 久久久久久久久久成人| 成人漫画全彩无遮挡| 欧美+日韩+精品| av在线天堂中文字幕| 在线国产一区二区在线| 久久精品影院6| 一夜夜www| 国产精品日韩av在线免费观看| 午夜爱爱视频在线播放| 人妻制服诱惑在线中文字幕| 久久婷婷人人爽人人干人人爱| 成人av在线播放网站| 少妇熟女aⅴ在线视频| 国产精品人妻久久久影院| 久久精品国产99精品国产亚洲性色| 国产高清不卡午夜福利| 我要看日韩黄色一级片| 国产 一区 欧美 日韩| 精品人妻偷拍中文字幕| 欧美日韩乱码在线| 亚洲在线自拍视频| 国产真实伦视频高清在线观看| 男插女下体视频免费在线播放| 男人的好看免费观看在线视频| 亚洲无线观看免费| 丰满人妻一区二区三区视频av| 能在线免费观看的黄片| 99视频精品全部免费 在线| 中文字幕熟女人妻在线| 精品人妻一区二区三区麻豆 | 十八禁网站免费在线| 久久久久久久久久黄片| 欧美色视频一区免费| 尾随美女入室| 99久久精品热视频| 熟女人妻精品中文字幕| 亚洲,欧美,日韩| 亚洲美女视频黄频| 欧美色视频一区免费| 一个人看视频在线观看www免费| 国产午夜精品久久久久久一区二区三区 | 少妇高潮的动态图| 两个人视频免费观看高清| 黄色配什么色好看| 欧美在线一区亚洲| h日本视频在线播放| 成人特级av手机在线观看| 国产激情偷乱视频一区二区| 国产欧美日韩精品亚洲av| 免费观看精品视频网站| 3wmmmm亚洲av在线观看| 国产欧美日韩一区二区精品| 熟妇人妻久久中文字幕3abv| 菩萨蛮人人尽说江南好唐韦庄 | 中国美白少妇内射xxxbb| 波多野结衣高清作品| 中文字幕人妻熟人妻熟丝袜美| 久久亚洲精品不卡| 午夜精品国产一区二区电影 | 国产又黄又爽又无遮挡在线| 亚洲av.av天堂| 免费人成在线观看视频色| 在线观看美女被高潮喷水网站| 亚洲专区国产一区二区| 亚洲人与动物交配视频| www日本黄色视频网| 精品99又大又爽又粗少妇毛片| 日韩欧美三级三区| 欧美最黄视频在线播放免费| а√天堂www在线а√下载| 亚洲欧美日韩东京热| 欧美在线一区亚洲| 日韩,欧美,国产一区二区三区 | 久久精品人妻少妇| 久久人人爽人人片av| 欧美日韩精品成人综合77777| 亚洲不卡免费看| 国产成人91sexporn| 别揉我奶头 嗯啊视频| 在线免费观看的www视频| 亚洲人成网站高清观看| 成年女人永久免费观看视频| 日日摸夜夜添夜夜添av毛片| or卡值多少钱| 午夜视频国产福利| 18禁在线播放成人免费| 国产精品国产高清国产av| 中国美女看黄片| 欧美日韩一区二区视频在线观看视频在线 | 内射极品少妇av片p| 69av精品久久久久久| 看十八女毛片水多多多| 我要搜黄色片| 亚洲精品一卡2卡三卡4卡5卡| 18禁裸乳无遮挡免费网站照片| 1000部很黄的大片| 中文字幕精品亚洲无线码一区| 无遮挡黄片免费观看| 精品午夜福利视频在线观看一区| www.色视频.com| 丝袜喷水一区| 久久久成人免费电影| 联通29元200g的流量卡| 亚洲性夜色夜夜综合| 国产精品国产高清国产av| 六月丁香七月| 国产熟女欧美一区二区| 欧美色视频一区免费| 日韩一本色道免费dvd| 成人鲁丝片一二三区免费| 精品国产三级普通话版| 99热这里只有是精品50| 国产国拍精品亚洲av在线观看| 大又大粗又爽又黄少妇毛片口| 精品午夜福利在线看| 天堂av国产一区二区熟女人妻| 97超视频在线观看视频| 日本一本二区三区精品| 国产成人91sexporn| 晚上一个人看的免费电影| 村上凉子中文字幕在线| 成人av在线播放网站| 一本久久中文字幕| 国产伦一二天堂av在线观看| 美女xxoo啪啪120秒动态图| 欧美+日韩+精品| 久久6这里有精品| 免费看光身美女| 嫩草影院入口| 国产视频一区二区在线看| 亚洲四区av| 嫩草影视91久久| 久久天躁狠狠躁夜夜2o2o| 亚洲久久久久久中文字幕| 男女下面进入的视频免费午夜| 亚洲国产高清在线一区二区三| 一边摸一边抽搐一进一小说| av福利片在线观看| 男人的好看免费观看在线视频| 日韩一本色道免费dvd| 最近2019中文字幕mv第一页| 日本色播在线视频| 中文字幕人妻熟人妻熟丝袜美| 国产精品三级大全| 久久精品国产亚洲av涩爱 | 国产精品爽爽va在线观看网站| a级一级毛片免费在线观看| 丰满乱子伦码专区| 伊人久久精品亚洲午夜| 97热精品久久久久久| 日本-黄色视频高清免费观看| 一进一出抽搐动态| 91av网一区二区| 久久人妻av系列| 国产精品三级大全| 国内精品一区二区在线观看| 男插女下体视频免费在线播放| 久久热精品热| 欧美成人免费av一区二区三区| 色吧在线观看| 亚洲中文日韩欧美视频| 日日干狠狠操夜夜爽| 婷婷精品国产亚洲av| 最近在线观看免费完整版| 国产午夜精品久久久久久一区二区三区 | 麻豆乱淫一区二区| 色在线成人网| 精品午夜福利视频在线观看一区| 久久久久久伊人网av| 少妇的逼水好多| av卡一久久| 日韩一本色道免费dvd| 欧美成人免费av一区二区三区| 国产蜜桃级精品一区二区三区| 精品国内亚洲2022精品成人| 国产精品爽爽va在线观看网站| 国产亚洲精品av在线| 狠狠狠狠99中文字幕| 99久久九九国产精品国产免费| 国产精华一区二区三区| 男女啪啪激烈高潮av片| 亚洲av电影不卡..在线观看| 97超视频在线观看视频| 国产色爽女视频免费观看| 九色成人免费人妻av| 欧美人与善性xxx| 久久精品国产清高在天天线| 亚洲人成网站在线观看播放| 最近中文字幕高清免费大全6| 欧美最新免费一区二区三区| 亚洲国产色片| 久久精品人妻少妇| 麻豆乱淫一区二区| 国产高清激情床上av|