• <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ì)配比對思茅松多穴盤幼苗生長的影響
    国产精品久久视频播放| 婷婷丁香在线五月| 麻豆一二三区av精品| 九色亚洲精品在线播放| 日本五十路高清| 变态另类丝袜制服| 一本久久中文字幕| 正在播放国产对白刺激| 欧美不卡视频在线免费观看 | 精品卡一卡二卡四卡免费| 黄网站色视频无遮挡免费观看| 97人妻精品一区二区三区麻豆 | 日韩av在线大香蕉| 亚洲欧美激情综合另类| 91精品国产国语对白视频| 国产99白浆流出| 精品一品国产午夜福利视频| 亚洲aⅴ乱码一区二区在线播放 | 性少妇av在线| 久久精品国产综合久久久| 午夜视频精品福利| 黄色视频不卡| 国产欧美日韩综合在线一区二区| 中文字幕av电影在线播放| 日本黄色视频三级网站网址| 男女床上黄色一级片免费看| 午夜精品在线福利| 夜夜躁狠狠躁天天躁| 亚洲色图av天堂| 老司机在亚洲福利影院| 一级a爱片免费观看的视频| 啦啦啦免费观看视频1| videosex国产| 热99re8久久精品国产| 免费在线观看日本一区| 悠悠久久av| 99热只有精品国产| 午夜精品在线福利| 免费一级毛片在线播放高清视频 | 国产成人欧美在线观看| 老司机深夜福利视频在线观看| 人人妻,人人澡人人爽秒播| 亚洲精品粉嫩美女一区| 国产麻豆成人av免费视频| 精品国内亚洲2022精品成人| 午夜视频精品福利| 亚洲性夜色夜夜综合| 人人妻人人澡人人看| 成人欧美大片| 久久久久亚洲av毛片大全| 国产黄a三级三级三级人| 后天国语完整版免费观看| 成人18禁高潮啪啪吃奶动态图| 在线观看66精品国产| 国产欧美日韩一区二区三| 久久精品91无色码中文字幕| 九色国产91popny在线| 黑人巨大精品欧美一区二区mp4| 亚洲精品久久成人aⅴ小说| 啪啪无遮挡十八禁网站| 亚洲精品国产色婷婷电影| av免费在线观看网站| 亚洲片人在线观看| 久久久久国内视频| 我的亚洲天堂| 亚洲 国产 在线| 90打野战视频偷拍视频| 两人在一起打扑克的视频| 午夜日韩欧美国产| 亚洲国产精品999在线| 亚洲精品中文字幕一二三四区| 中文字幕色久视频| 手机成人av网站| 男女下面进入的视频免费午夜 | 香蕉国产在线看| 精品一区二区三区av网在线观看| 午夜视频精品福利| 久久中文字幕一级| 国产精品久久电影中文字幕| 精品一区二区三区视频在线观看免费| 19禁男女啪啪无遮挡网站| 99香蕉大伊视频| 一级黄色大片毛片| 乱人伦中国视频| 嫩草影院精品99| 日本撒尿小便嘘嘘汇集6| 亚洲午夜理论影院| 精品久久久久久,| 最近最新中文字幕大全免费视频| 丰满人妻熟妇乱又伦精品不卡| 国产精品1区2区在线观看.| 男女做爰动态图高潮gif福利片 | 国产精品99久久99久久久不卡| 久久亚洲精品不卡| 在线观看免费日韩欧美大片| 精品久久久久久久人妻蜜臀av | 国产av又大| 免费观看精品视频网站| 女人高潮潮喷娇喘18禁视频| av视频在线观看入口| 亚洲avbb在线观看| 成人特级黄色片久久久久久久| 久久精品影院6| 在线观看免费午夜福利视频| 无限看片的www在线观看| 精品久久久久久,| 99国产精品一区二区蜜桃av| 国产日韩一区二区三区精品不卡| 91字幕亚洲| 99在线人妻在线中文字幕| 亚洲人成电影免费在线| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕人成人乱码亚洲影| 这个男人来自地球电影免费观看| 欧美一级a爱片免费观看看 | 午夜福利高清视频| 免费观看人在逋| 免费久久久久久久精品成人欧美视频| 亚洲成人久久性| av有码第一页| 91麻豆av在线| 久久久久久免费高清国产稀缺| 久久久水蜜桃国产精品网| 无人区码免费观看不卡| 757午夜福利合集在线观看| 国产高清激情床上av| 日韩大码丰满熟妇| 一级毛片高清免费大全| 午夜福利欧美成人| 一区在线观看完整版| 波多野结衣巨乳人妻| 亚洲少妇的诱惑av| 久久精品91无色码中文字幕| 精品国内亚洲2022精品成人| 久久狼人影院| 给我免费播放毛片高清在线观看| svipshipincom国产片| 啦啦啦观看免费观看视频高清 | 国产av精品麻豆| 免费看美女性在线毛片视频| 大香蕉久久成人网| 亚洲欧美精品综合久久99| 免费观看精品视频网站| 9热在线视频观看99| 国产欧美日韩一区二区三区在线| 精品一品国产午夜福利视频| 大陆偷拍与自拍| 日韩大码丰满熟妇| 欧美成人免费av一区二区三区| 精品久久久久久久人妻蜜臀av | 国产成人av教育| 女人爽到高潮嗷嗷叫在线视频| 97人妻精品一区二区三区麻豆 | 免费在线观看亚洲国产| 亚洲最大成人中文| 丰满的人妻完整版| 国产麻豆69| 国产精品一区二区免费欧美| 久久精品亚洲熟妇少妇任你| 99热只有精品国产| 国产单亲对白刺激| 亚洲中文字幕一区二区三区有码在线看 | 亚洲精华国产精华精| aaaaa片日本免费| 久久青草综合色| 国产乱人伦免费视频| 一进一出抽搐gif免费好疼| 宅男免费午夜| 欧美成人性av电影在线观看| 在线观看日韩欧美| www.熟女人妻精品国产| 99久久99久久久精品蜜桃| 久久久久久人人人人人| 久热这里只有精品99| 国产亚洲精品久久久久5区| 国产熟女午夜一区二区三区| 国内久久婷婷六月综合欲色啪| 国产乱人伦免费视频| 99在线人妻在线中文字幕| 久久精品国产99精品国产亚洲性色 | 国产xxxxx性猛交| 一边摸一边抽搐一进一小说| 黑人操中国人逼视频| 国产又色又爽无遮挡免费看| 18禁黄网站禁片午夜丰满| 美女国产高潮福利片在线看| av有码第一页| 色尼玛亚洲综合影院| 悠悠久久av| 狠狠狠狠99中文字幕| 国产日韩一区二区三区精品不卡| 久久国产精品男人的天堂亚洲| 一本久久中文字幕| 少妇熟女aⅴ在线视频| 宅男免费午夜| 久久草成人影院| 中文亚洲av片在线观看爽| 天天添夜夜摸| 国产成人av激情在线播放| 最好的美女福利视频网| 18禁裸乳无遮挡免费网站照片 | 美女国产高潮福利片在线看| 国产成人精品在线电影| 最新美女视频免费是黄的| 亚洲少妇的诱惑av| 久9热在线精品视频| 欧美绝顶高潮抽搐喷水| 久久久水蜜桃国产精品网| 99久久综合精品五月天人人| 老熟妇仑乱视频hdxx| а√天堂www在线а√下载| 国内精品久久久久精免费| 日韩成人在线观看一区二区三区| av超薄肉色丝袜交足视频| 欧美日韩亚洲国产一区二区在线观看| 天天躁夜夜躁狠狠躁躁| 久久久久久久久久久久大奶| 美国免费a级毛片| 午夜福利,免费看| 两人在一起打扑克的视频| 成人特级黄色片久久久久久久| 香蕉丝袜av| 最新美女视频免费是黄的| 涩涩av久久男人的天堂| 麻豆国产av国片精品| 老司机深夜福利视频在线观看| 亚洲久久久国产精品| 亚洲一区高清亚洲精品| 亚洲精品国产区一区二| 电影成人av| 欧美一级a爱片免费观看看 | 村上凉子中文字幕在线| 真人一进一出gif抽搐免费| 免费在线观看视频国产中文字幕亚洲| 亚洲精品美女久久av网站| 国产亚洲欧美98| 亚洲欧美精品综合一区二区三区| 国产精品,欧美在线| 免费在线观看黄色视频的| 亚洲国产欧美网| 亚洲av成人不卡在线观看播放网| 十八禁人妻一区二区| 搡老熟女国产l中国老女人| 国产精品免费一区二区三区在线| 可以免费在线观看a视频的电影网站| 亚洲精品中文字幕在线视频| 久久精品成人免费网站| 激情在线观看视频在线高清| 国产精品精品国产色婷婷| 成人18禁高潮啪啪吃奶动态图| 一级a爱片免费观看的视频| 日韩有码中文字幕| 久久精品91蜜桃| av在线天堂中文字幕| 黄色a级毛片大全视频| 久久国产精品男人的天堂亚洲| 国产欧美日韩综合在线一区二区| 在线观看免费视频日本深夜| 一级片免费观看大全| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品美女久久久久99蜜臀| 国产激情久久老熟女| 精品人妻1区二区| 真人一进一出gif抽搐免费| 精品福利观看| 狠狠狠狠99中文字幕| 精品久久久久久久人妻蜜臀av | 免费久久久久久久精品成人欧美视频| 亚洲国产毛片av蜜桃av| 非洲黑人性xxxx精品又粗又长| 亚洲国产高清在线一区二区三 | 久久午夜亚洲精品久久| 久热爱精品视频在线9| 欧美日韩亚洲综合一区二区三区_| 欧美一级a爱片免费观看看 | 欧美黑人精品巨大| 99久久综合精品五月天人人| 涩涩av久久男人的天堂| 国产麻豆69| 麻豆久久精品国产亚洲av| 久久亚洲真实| 在线十欧美十亚洲十日本专区| 女生性感内裤真人,穿戴方法视频| svipshipincom国产片| 亚洲国产精品sss在线观看| 国产麻豆成人av免费视频| 一区二区三区激情视频| 老鸭窝网址在线观看| 欧美乱码精品一区二区三区| 久久久久久久午夜电影| 免费观看精品视频网站| 一二三四社区在线视频社区8| www.精华液| 久久国产乱子伦精品免费另类| 久久香蕉激情| 妹子高潮喷水视频| 一本久久中文字幕| 中文亚洲av片在线观看爽| 色播亚洲综合网| 国产免费男女视频| 国产av精品麻豆| 久久国产亚洲av麻豆专区| 亚洲第一青青草原| 天堂影院成人在线观看| 一区福利在线观看| 精品日产1卡2卡| 亚洲自拍偷在线| 很黄的视频免费| 久久婷婷成人综合色麻豆| 国产成人影院久久av| 欧美激情久久久久久爽电影 | 侵犯人妻中文字幕一二三四区| 午夜亚洲福利在线播放| 国产免费男女视频| 亚洲精品中文字幕一二三四区| 99久久99久久久精品蜜桃| 国产一区二区在线av高清观看| 久久午夜综合久久蜜桃| 免费av毛片视频| 成熟少妇高潮喷水视频| 99re在线观看精品视频| 欧美另类亚洲清纯唯美| 国产一区在线观看成人免费| 久久久精品国产亚洲av高清涩受| 天天躁夜夜躁狠狠躁躁| 三级毛片av免费| 午夜视频精品福利| 制服诱惑二区| 国产亚洲精品久久久久久毛片| 老鸭窝网址在线观看| 亚洲一区中文字幕在线| 久久久久久久精品吃奶| 午夜激情av网站| cao死你这个sao货| 美女国产高潮福利片在线看| 母亲3免费完整高清在线观看| 美国免费a级毛片| 国产一卡二卡三卡精品| 亚洲中文av在线| 久久久久久久午夜电影| 国产精品美女特级片免费视频播放器 | 国产国语露脸激情在线看| 欧美绝顶高潮抽搐喷水| 欧美老熟妇乱子伦牲交| 亚洲人成电影免费在线| 久久精品亚洲熟妇少妇任你| 99热只有精品国产| 在线av久久热| 国产99白浆流出| 国产精品av久久久久免费| 欧美国产日韩亚洲一区| 亚洲 国产 在线| 亚洲五月天丁香| 欧美 亚洲 国产 日韩一| 丝袜美腿诱惑在线| 一本大道久久a久久精品| 天天添夜夜摸| 亚洲自偷自拍图片 自拍| av在线天堂中文字幕| 在线观看舔阴道视频| 91成年电影在线观看| 一级毛片精品| 亚洲av日韩精品久久久久久密| 黄色a级毛片大全视频| 亚洲成av片中文字幕在线观看| 久久人人爽av亚洲精品天堂| 亚洲熟妇熟女久久| 色综合婷婷激情| 国产一区二区激情短视频| 免费在线观看日本一区| 久久久久久大精品| 成人永久免费在线观看视频| 性色av乱码一区二区三区2| 不卡av一区二区三区| 日本免费a在线| 欧美日本视频| 日日干狠狠操夜夜爽| 99国产精品一区二区三区| а√天堂www在线а√下载| 少妇熟女aⅴ在线视频| 国产成人精品无人区| 亚洲成av片中文字幕在线观看| 淫妇啪啪啪对白视频| 十八禁网站免费在线| 99精品在免费线老司机午夜| 国产精品爽爽va在线观看网站 | 国产精品国产高清国产av| 亚洲国产精品成人综合色| 久久久久国产精品人妻aⅴ院| 日本a在线网址| 国产精品乱码一区二三区的特点 | 精品国产超薄肉色丝袜足j| 午夜精品在线福利| 欧美+亚洲+日韩+国产| 极品人妻少妇av视频| av中文乱码字幕在线| 性欧美人与动物交配| 久久国产精品男人的天堂亚洲| 一卡2卡三卡四卡精品乱码亚洲| 久久狼人影院| 国产aⅴ精品一区二区三区波| 亚洲片人在线观看| 制服人妻中文乱码| 国产欧美日韩精品亚洲av| 日日摸夜夜添夜夜添小说| 露出奶头的视频| 亚洲一区二区三区不卡视频| 免费人成视频x8x8入口观看| 日韩中文字幕欧美一区二区| 在线观看免费午夜福利视频| 午夜福利18| 一卡2卡三卡四卡精品乱码亚洲| 婷婷六月久久综合丁香| 久久人妻av系列| 岛国视频午夜一区免费看| 精品一区二区三区av网在线观看| 成人亚洲精品一区在线观看| 国产精品乱码一区二三区的特点 | 99久久久亚洲精品蜜臀av| 99久久国产精品久久久| 亚洲美女黄片视频| 日本在线视频免费播放| 深夜精品福利| 丝袜美足系列| 搡老妇女老女人老熟妇| 国产精品98久久久久久宅男小说| 嫁个100分男人电影在线观看| 国产一区二区三区视频了| 久久香蕉激情| 久久久国产欧美日韩av| 亚洲精品国产一区二区精华液| 两个人免费观看高清视频| 黑人欧美特级aaaaaa片| 91成年电影在线观看| 国产免费男女视频| 久久精品成人免费网站| 欧美一区二区精品小视频在线| 国产麻豆成人av免费视频| 亚洲美女黄片视频| 高清在线国产一区| 中文字幕人妻丝袜一区二区| 黑人操中国人逼视频| 欧美成人性av电影在线观看| 亚洲中文字幕日韩| 人人妻人人澡人人看| 精品久久久久久,| √禁漫天堂资源中文www| 国产私拍福利视频在线观看| 国产午夜福利久久久久久| 久久香蕉激情| 亚洲欧美日韩无卡精品| 色尼玛亚洲综合影院| 国产精品一区二区三区四区久久 | 国产精品99久久99久久久不卡| 这个男人来自地球电影免费观看| 日韩免费av在线播放| 午夜精品在线福利| 亚洲熟女毛片儿| 国产精品精品国产色婷婷| 日本 欧美在线| 一边摸一边抽搐一进一出视频| 满18在线观看网站| 色婷婷久久久亚洲欧美| 久久人妻熟女aⅴ| 国产亚洲精品一区二区www| 黄色视频,在线免费观看| 色播在线永久视频| 亚洲av日韩精品久久久久久密| 日日夜夜操网爽| 岛国视频午夜一区免费看| 两个人看的免费小视频| 亚洲第一青青草原| 人妻丰满熟妇av一区二区三区| 国产欧美日韩精品亚洲av| 亚洲欧美日韩高清在线视频| 欧美黄色片欧美黄色片| 久久国产精品人妻蜜桃| 精品人妻在线不人妻| 亚洲男人的天堂狠狠| 精品一品国产午夜福利视频| 两人在一起打扑克的视频| 午夜福利高清视频| 亚洲五月天丁香| 99久久99久久久精品蜜桃| 国产一区二区在线av高清观看| 亚洲色图av天堂| 嫩草影院精品99| 黄色毛片三级朝国网站| 日本精品一区二区三区蜜桃| www.熟女人妻精品国产| 午夜福利一区二区在线看| 午夜福利成人在线免费观看| 亚洲成人精品中文字幕电影| 桃色一区二区三区在线观看| 国产精品久久久久久亚洲av鲁大| bbb黄色大片| 精品电影一区二区在线| av福利片在线| 久久久久久人人人人人| 久9热在线精品视频| 久久婷婷成人综合色麻豆| 国产欧美日韩一区二区精品| 亚洲国产欧美日韩在线播放| 两性夫妻黄色片| 色哟哟哟哟哟哟| 精品国产一区二区久久| 国产亚洲精品第一综合不卡| 99久久精品国产亚洲精品| av中文乱码字幕在线| aaaaa片日本免费| 精品人妻1区二区| aaaaa片日本免费| 波多野结衣高清无吗| 人人妻人人澡人人看| ponron亚洲| 99久久国产精品久久久| 久久人人97超碰香蕉20202| 久久香蕉激情| 丝袜在线中文字幕| 国产精华一区二区三区| av免费在线观看网站| 1024视频免费在线观看| 手机成人av网站| 啪啪无遮挡十八禁网站| av免费在线观看网站| 丝袜人妻中文字幕| 变态另类成人亚洲欧美熟女 | 精品久久久久久久久久免费视频| 日韩欧美一区二区三区在线观看| www.精华液| av视频免费观看在线观看| 国内久久婷婷六月综合欲色啪| 国产av一区二区精品久久| 高清在线国产一区| 久久中文看片网| 成人18禁高潮啪啪吃奶动态图| 99精品欧美一区二区三区四区| 50天的宝宝边吃奶边哭怎么回事| 国产欧美日韩精品亚洲av| 国产亚洲精品综合一区在线观看 | 欧美成人性av电影在线观看| 老汉色∧v一级毛片| 欧美成人免费av一区二区三区| 天堂√8在线中文| 99国产精品一区二区三区| 欧美色视频一区免费| 日韩成人在线观看一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 最近最新中文字幕大全免费视频| 一区二区三区高清视频在线| 日韩欧美在线二视频| 国产成人精品在线电影| 亚洲欧美激情在线| 国产单亲对白刺激| 欧美日本中文国产一区发布| 亚洲色图av天堂| 精品无人区乱码1区二区| 久久精品国产清高在天天线| 国产成+人综合+亚洲专区| 禁无遮挡网站| 国产一区在线观看成人免费| 精品国内亚洲2022精品成人| 老鸭窝网址在线观看| 国产精品影院久久| 九色国产91popny在线| 国产伦人伦偷精品视频| 色哟哟哟哟哟哟| 熟女少妇亚洲综合色aaa.| 国产一区二区在线av高清观看| 国产精品98久久久久久宅男小说| 欧美色欧美亚洲另类二区 | 中文字幕人妻丝袜一区二区| 女性被躁到高潮视频| 18禁黄网站禁片午夜丰满| 国产欧美日韩一区二区三| 美国免费a级毛片| 欧美色视频一区免费| 亚洲精品久久成人aⅴ小说| 精品卡一卡二卡四卡免费| 亚洲五月婷婷丁香| 免费少妇av软件| 看片在线看免费视频| 又大又爽又粗| 美女大奶头视频| 黄色毛片三级朝国网站| 丝袜人妻中文字幕| 一级a爱片免费观看的视频| 久久 成人 亚洲| 18美女黄网站色大片免费观看| av天堂在线播放| 欧美一级a爱片免费观看看 | 精品国产乱码久久久久久男人| 亚洲精品在线美女| 婷婷丁香在线五月| 91成年电影在线观看| 亚洲精品在线美女| 免费久久久久久久精品成人欧美视频| 精品久久久久久久毛片微露脸| 国产精品98久久久久久宅男小说| 国产区一区二久久| 变态另类丝袜制服| 精品不卡国产一区二区三区| 精品欧美一区二区三区在线| 成人亚洲精品一区在线观看| 精品一区二区三区视频在线观看免费| 国产亚洲欧美精品永久| 亚洲人成77777在线视频| 亚洲一区二区三区色噜噜| √禁漫天堂资源中文www| 久久人妻福利社区极品人妻图片| 999精品在线视频| 亚洲男人的天堂狠狠| 久久中文字幕一级| ponron亚洲| 精品第一国产精品|