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

    黃龍山天然次生林遼東櫟種群空間格局

    2016-07-21 01:28:23王道亮李衛(wèi)忠王廣儒
    生態(tài)學(xué)報 2016年9期

    王道亮,李衛(wèi)忠,*,曹 鑄,王 慶,王廣儒

    1 西北農(nóng)林科技大學(xué), 林學(xué)院,楊凌 712100 2 延安市黃龍山林業(yè)局,黃龍 715700

    ?

    黃龍山天然次生林遼東櫟種群空間格局

    王道亮1,李衛(wèi)忠1,*,曹鑄1,王慶1,王廣儒2

    1 西北農(nóng)林科技大學(xué), 林學(xué)院,楊凌712100 2 延安市黃龍山林業(yè)局,黃龍715700

    摘要:種群的空間分布格局是研究種群特征、種間相互作用以及種群與環(huán)境關(guān)系的重要手段。為更好地了解黃龍山天然次生林內(nèi)遼東櫟(Quercus wutaishanica)的種群特征,在該林區(qū)中油松(Pinus tabulaeformis)-遼東櫟-白樺(Betula platyphylla)天然次生混交林和遼東櫟天然次生純林內(nèi)分別選取并設(shè)立1塊1 hm2的樣地。根據(jù)樣地調(diào)查規(guī)范對兩塊樣地中的喬木進(jìn)行空間定位并進(jìn)行每木調(diào)查。研究參照靜態(tài)生命表,采用徑階替代年齡的方法,將遼東櫟種群不同生長階段的分段問題轉(zhuǎn)換為根據(jù)胸徑(DBH)的大小劃分生長階段。研究中將遼東櫟劃分為4個不同生長階段:幼齡階段A, 3 cm ≤DBH<10 cm;中齡階段B, 10 cm ≤DBH<20 cm;近熟齡階段C, 20 cm ≤DBH<30 cm;成過熟齡階段D, DBH≥30 cm。利用O-ring統(tǒng)計分析方法的手段,分析2 塊樣地中遼東櫟的種群空間分布格局。本次研究分別采用單變量、雙變量分析方法分析各生長階段的空間格局、不同生長階段間相互關(guān)聯(lián)性。結(jié)果表明:(1) 遼東櫟種群在兩塊樣地內(nèi)的徑級結(jié)構(gòu)均符合倒“J”型,從徑級分布角度說明兩塊樣地的天然更新狀況都良好;(2) 遼東櫟的空間分布格局與尺度有密切關(guān)系, 兩塊樣地在大、中尺度水平上基本表現(xiàn)為隨機(jī)分布,但是在較小尺度上均表現(xiàn)為聚集分布,在樣地Ⅰ中0—6 m尺度上尤為突出;(3) 同一樣地內(nèi)各生長階段的空間分布格局不同,幼齡個體在小尺度均表現(xiàn)為聚集分布,在樣地Ⅱ中0—11 m尺度上尤為突出;(4)樣地中不同生長階段間相互關(guān)聯(lián)性也不盡相同,樣地Ⅰ中各生長階段間關(guān)聯(lián)性均不明顯;樣地Ⅱ中,幼齡的個體與中齡的個體間在小尺度上表現(xiàn)出顯著的正相關(guān)性, 而與近熟個體和成、過熟個體在小尺度上表現(xiàn)為明顯的負(fù)相關(guān)性。在2—5 m尺度上,中齡個體與成過熟個體也表現(xiàn)出一定的負(fù)相關(guān)性。遼東櫟種群是黃龍山天然次生林重要的建群種之一,研究其不同生長階段個體的空間分布規(guī)律有助于掌握黃龍山區(qū)遼東櫟生長現(xiàn)狀;研究不同生長階段的空間關(guān)聯(lián)性有利于深入了解其群落的發(fā)展趨勢,并為對種群聚集性較大的區(qū)域開展結(jié)構(gòu)化森林經(jīng)營提供重要依據(jù)。

    關(guān)鍵詞:天然次生林;空間分布格局;空間關(guān)聯(lián)性;O-ring統(tǒng)計分析;遼東櫟

    物種是森林結(jié)構(gòu)的重要組成,其發(fā)育的形態(tài)不僅與所包含的物種相關(guān),更與當(dāng)?shù)氐淖匀画h(huán)境因子有重要關(guān)系[1]。森林的物種組成特征一方面表征著森林的發(fā)育程度,另一方面也與森林生態(tài)功能的發(fā)揮有著不可分割的聯(lián)系[2- 4]。種群空間分布格局是在種內(nèi)競爭、種間競爭及自然環(huán)境等各種因素共同影響下形成的結(jié)果[5- 7],對其展開深入分析,不僅可以對物種分布規(guī)律有一定的了解,還可以分析推斷其發(fā)育過程[4,8- 9],是研究種群特征、種間相互作用以及種群與環(huán)境關(guān)系的重要手段[2,5],有利于更好地理解種群格局形成的潛在生態(tài)過程[5]。

    黃龍山林區(qū)地處黃土高原溝壑區(qū),為陜西渭北旱塬與陜北黃土高原的過渡帶, 被譽(yù)為“黃土高原的綠色明珠”,是我國水土保持的重要防護(hù)林區(qū)。遼東櫟是黃龍山林區(qū)重要的優(yōu)勢種之一,屬暖溫帶落葉闊葉櫟林的重要樹種[10]。

    目前,在黃龍山林區(qū),對遼東櫟的研究內(nèi)容主要集中于生態(tài)位[11]、生物多樣性[12- 13]、群落穩(wěn)定性[14]等方面,對于其不同生長階段的空間分布及其空間關(guān)聯(lián)性尚未見報道。以往文獻(xiàn)主要以胸徑為依據(jù),通過叢生法[15]、方差/均值比法[16]、g(r)函數(shù)法[17- 18]等方法對物種空間格局分布進(jìn)行研究。O-ring函數(shù)在優(yōu)勢種群的空間分布及種間關(guān)聯(lián)性研究方面日漸多見[5,19- 20],而將O-ring函數(shù)運用于單一樹種不同生長階段之間關(guān)聯(lián)性的文章卻不多見[21]。本次研究選取黃龍山林區(qū)內(nèi)遼東櫟林和油松-遼東櫟-白樺混交林,分別選設(shè)一塊1 hm2級樣地,并對符合調(diào)查條件的喬木進(jìn)行每木檢尺。本文根據(jù)遼東櫟徑階,結(jié)合靜態(tài)生命表[22],按照齡組的劃分方法,將其劃分為不同生長階段,并運用O-ring點格局分析法進(jìn)行空間格局分析及空間關(guān)聯(lián)性分析,以期對種群聚集性較大的區(qū)域?qū)嵤┙Y(jié)構(gòu)化森林經(jīng)營。

    1研究區(qū)概況

    研究區(qū)位于陜西省黃龍山天然林保護(hù)區(qū)內(nèi),地理坐標(biāo)為東經(jīng)109°38′49″—110°12′47″ E,北緯35°28′46″—36°02′01″ N,位于溫帶半濕潤與半干旱氣候的過渡地帶,屬黃土高原丘陵溝壑區(qū),海拔1100—1300 m,年平均降水量為611.8 mm,年均氣溫為8.6 ℃,地帶性土壤為褐土和灰褐色森林土[23],森林植被屬暖溫帶落葉混交林地帶,主要優(yōu)勢種包括油松、遼東櫟、白樺等,天然植被覆蓋率高。

    2研究方法

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

    根據(jù)空間格局研究樣地的要求,樣地必須有一定的規(guī)模才能保證物種的各種格局均能出現(xiàn)[4],且樣地應(yīng)以長方形或正方形最佳[24]。2013 年8 月,根據(jù)樣地設(shè)置技術(shù)要求及研究目標(biāo)[25],通過對虎溝門林場進(jìn)行全面踏查,在該林場設(shè)立了2 塊遼東櫟生長狀況良好并且具有典型特征的樣地:樣地Ⅰ油松-遼東櫟-白樺混交林樣地和樣地Ⅱ遼東櫟純林樣地,其樣地面積均為1 hm2。采用相鄰網(wǎng)格法,用高精度的拓普康全站儀將樣地劃分為25 塊20×20 m2的小樣地,對兩塊樣地內(nèi)喬木層展開調(diào)查,主要調(diào)查因子包含樹種、胸徑(DBH) (cm)、樹高(m)、冠幅(m)、生長狀態(tài)等,同時調(diào)查記錄各林木的空間位置和坐標(biāo)以及環(huán)境因子如土壤、枯落物等,并進(jìn)行逐株掛牌。樣地統(tǒng)計結(jié)果如表1。

    根據(jù)調(diào)查結(jié)果可知:樣地Ⅰ中,喬木層主要由油松、遼東櫟、白樺等組成,灌木層主要由華北紫丁香(Syringaoblata)、茶條槭(Acerginnala)等組成,草本主要由大披針苔草(Carexlanceolata)、高羊茅(Festucaarundinacea)等組成。樣地Ⅱ中,喬木層主要由遼東櫟組成,灌木層主要由野山楂(Crataeguscuneata)、土莊繡線菊(Spiraeapubescens)等組成,草本層主要由大披針苔草、高羊茅、紫菀(Astertataricus)等組成。

    表1 樣地調(diào)查概況表

    遼東櫟樹體高大,木質(zhì)堅硬,因而不適合利用生長錐法等傳統(tǒng)方法進(jìn)行年齡的測定[26],而在相似立地條件和自然環(huán)境下的相同物種其齡級和徑級對環(huán)境的反映規(guī)律具有一致性[27],樣地調(diào)查開始之前,對樣地林分平均胸徑進(jìn)行估測。經(jīng)初步預(yù)測,本研究選取的兩塊樣地平均胸徑均不大于12 cm,結(jié)合林分結(jié)構(gòu)規(guī)律[28],以林分平均胸徑的0.4 倍所在徑階作為起測徑階。因此,文本起測徑階為4徑階,從而確定起測胸徑為3.0 cm。根據(jù)樣地調(diào)查結(jié)果,結(jié)合遼東櫟種群靜態(tài)生命表[22],將其劃分為A—D 4個階段:幼齡階段A, 3 cm ≤DBH<10 cm;中齡階段B, 10 cm ≤DBH<20 cm;近熟齡階段C, 20 cm ≤DBH<30 cm;成過熟齡階段D, DBH≥30 cm。其分布結(jié)果如圖1。

    圖1 遼東櫟空間分布圖Fig.1 Spatial distribution of Quercus wutaishanica

    2.2數(shù)據(jù)處理

    點格局分析法(point pattern analysis)是以植物種的個體在空間的坐標(biāo)為基本數(shù)據(jù),以其空間分布點圖為基礎(chǔ)進(jìn)行格局分析。該方法廣泛應(yīng)用于植物種群多尺度空間分布格局及兩個物種之間多尺度空間關(guān)聯(lián)的研究。

    O-ring函數(shù)是在RipleyK(d)函數(shù)的基礎(chǔ)上,結(jié)合Mark相關(guān)函數(shù)的一種新的點格局分析方法[5]。O-ring采用環(huán)形尺度分析方法,以寬度為1 m的環(huán)帶替代圓形尺度分析,消除了在RipleyK(d)函數(shù)以圓形1 m步長為尺度的研究分析所帶來的大小尺度的累計效應(yīng)[29]。因此,本研究中采用單變量的O-ring函數(shù)[O(r)]分析不同生長階段個體空間分布格局,用雙變量的O-ring函數(shù)[O12(r)]分析兩個階段間的空間關(guān)聯(lián)性。

    零假設(shè)模型的篩選是保證數(shù)據(jù)分析正確與否的關(guān)鍵[5],因此在數(shù)據(jù)預(yù)處理階段應(yīng)仔細(xì)篩選。在進(jìn)行單變量O-ring分析時,根據(jù)蒙特卡羅方法求得各生長階段空間分布的置信區(qū)間,如果函數(shù)值處在置信區(qū)間,則接受零假設(shè),空間分布為隨機(jī)分布;如果在置信區(qū)間以上,則為聚集分布;如果函數(shù)值在置信區(qū)間以下,則為均勻分布。在群落的發(fā)育過程中,高齡個體先于低齡個體出現(xiàn),而低齡個體空間的分布受高齡個體分布的影響。假設(shè)種群在發(fā)育過程中高齡個體對低齡個體有一定的庇護(hù)作用[30],采用先決條件零假設(shè)模型進(jìn)行雙變量O-ring空間關(guān)聯(lián)分析。根據(jù)蒙特卡羅模擬結(jié)果,如果函數(shù)值在置信區(qū)間之間,表明接受零假設(shè),二者之間相互獨立;如果在置信區(qū)間以上,則二者空間上正關(guān)聯(lián);如果在置信區(qū)間以下,則二者空間上負(fù)關(guān)聯(lián)。

    數(shù)據(jù)處理過程主要通過Programita進(jìn)行,而不同生長階段的點格局分布則在GS+ 9.0軟件中實現(xiàn)。采用空間尺度為0—50 m,根據(jù)相應(yīng)的零假設(shè)檢驗?zāi)P停M(jìn)行99次MentoKalo(蒙特卡羅)模擬,得到上下包跡線構(gòu)成的99%置信區(qū)間[19,31- 34]。

    3結(jié)果與分析

    3.1遼東櫟種群特征及生長階段劃分

    經(jīng)調(diào)查(圖2),兩塊樣地中遼東櫟各徑級個體樹木統(tǒng)計如圖2顯示,在樣地Ⅰ、Ⅱ中,遼東櫟在幼齡、中齡、近熟齡和成過熟齡階段個體數(shù)目分別為:349、161、19、20 株; 241、242、139、46 株。由生長趨勢分布可以看出,樣地Ⅰ內(nèi)遼東櫟種群徑級結(jié)構(gòu)分布更加接近于倒“J”型,以小徑級的幼齡樹居多,更新狀態(tài)良好;樣地Ⅱ中幼齡林木略少于樣地Ⅰ,且從其空間分布狀況對比可知(圖1),樣地Ⅰ中幼樹簇生現(xiàn)象更加明顯。

    圖2 樣地Ⅰ和樣地Ⅱ遼東櫟徑級分布圖 Fig.2 Distribution of growth stages of Quercus wutaishanica in PlotⅠ and Polt Ⅱ

    樣地Ⅰ內(nèi)林木達(dá)到檢尺胸徑3 cm的活立木2432 株,其中,油松808 株,遼東櫟549 株,白樺481 株。樣地內(nèi)喬木平均胸徑為9.86 cm,胸高總斷面積為18.57 m2/hm2。遼東櫟的平均胸徑為10.22 cm,平均樹高為8.21 m,胸高斷面積為4.50 m2/hm2。樣地Ⅱ內(nèi)林木達(dá)到檢尺胸徑3 cm的活立木1018 株,其中遼東櫟有668 株。樣地內(nèi)喬木平均胸徑為12.98 cm,胸高總斷面積為13.47 m2/hm2,遼東櫟的平均胸徑為15.36 cm,平均樹高為8.03 m,胸高斷面積為12.38 m2/hm2。通過對比我們可以得出,遼東櫟的平均胸徑在樣地Ⅰ要小于樣地Ⅱ,而樹高正好相反,樣地Ⅰ內(nèi)遼東櫟要比樣地Ⅱ內(nèi)遼東櫟更高。分析其原因:在相似立地條件下,樣地Ⅰ內(nèi)林木株數(shù)密度要大于樣地Ⅱ,樣地Ⅰ內(nèi)林木競爭環(huán)境要比樣地Ⅱ中苛刻,對其林木生長和分布也產(chǎn)生一定影響。且樣地Ⅰ、Ⅱ內(nèi)郁閉度分別為0.79、0.48,在相同光照條件下,下層林木對光照的獲得量差異較為明顯;在競爭過程中,林木個體向高處生長以獲得更多的光資源。而在樣地Ⅱ中在光照充足的情況下這種競爭相對不明顯,且林木發(fā)育狀況比樣地Ⅰ中要好。因此,混交林中遼東櫟與競爭相對較弱的遼東櫟純林相比,樹木表現(xiàn)出高而細(xì)的特點。

    3.2遼東櫟種群空間分布格局

    根據(jù)調(diào)查結(jié)果(圖1)明顯可以看出:樣地Ⅰ幼齡樹明顯多于樣地Ⅱ,且有多處為萌生叢生分布,十分符合次生林的特點。通過對兩塊樣地進(jìn)行單變量O-ring統(tǒng)計分析發(fā)現(xiàn)(圖3),隨尺度變化,其空間分布格局具有不同的特點:樣地Ⅰ內(nèi)遼東櫟在尺度6 m以內(nèi)隨著研究尺度的增大聚集強(qiáng)度逐漸減弱,在11 m、16 m、35—37 m處表現(xiàn)為輕度聚集,其他尺度則均以隨機(jī)分布為主。樣地Ⅱ內(nèi)遼東櫟主要表現(xiàn)為隨機(jī)分布,僅在1—1.8 m之間,3—5.7 m之間,7—9 m之間表現(xiàn)為聚集分布。

    圖3 樣地Ⅰ和樣地Ⅱ空間分布格局Fig.3 Spatial pattern of Quercus wutaishanica in PlotⅠ and Polt Ⅱ

    3.3遼東櫟種群不同齡級的空間分布格局

    兩個樣地的遼東櫟在4 個不同的生長階段的點格局分析表明(圖4,圖5),兩個樣地中遼東櫟胸徑小于20 cm的幼齡和中齡個體空間分布稍有不同:樣地Ⅰ中,遼東櫟個體基本表現(xiàn)為隨機(jī)分布;樣地Ⅱ中胸徑位于3—10 cm之間的個體在11 m以內(nèi)為聚集分布,11 m以外為隨機(jī)分布,10—20 cm之間的個體在2 m以內(nèi)及3—4 m、5.4—7 m之間表現(xiàn)為聚集分布,其他為隨機(jī)分布。在胸徑大于20 cm的近熟和成過熟階段個體,均表現(xiàn)為隨機(jī)分布,并有呈均勻分布的趨勢。

    圖4 樣地Ⅰ各齡級空間分布格局Fig.4 Spatial pattern of each stage of Quercus wutaishanica in plotⅠA: 幼齡階段A; B: 中齡階段B; C: 近熟齡階段C; D: 成過熟齡D

    圖5 樣地Ⅱ各齡級空間分布格局Fig.5 Spatial pattern of each stage of Quercus wutaishanica in plot ⅡA: 幼齡階段A; B: 中齡階段B; C: 近熟齡階段C; D: 成過熟齡D

    3.4遼東櫟種群不同齡級間的空間關(guān)聯(lián)性

    通過對兩塊樣地進(jìn)行雙變量分析結(jié)果顯示(圖6,圖7),樣地Ⅰ中遼東櫟各生長階段之間關(guān)聯(lián)性均不明顯,而樣地Ⅱ中各生長階段之間表現(xiàn)不同:樣地Ⅱ中幼齡個體與中齡個體間在小尺度上表現(xiàn)出顯著的正相關(guān)性,而與近熟齡個體和成過熟齡個體在小尺度上表現(xiàn)為明顯的負(fù)相關(guān)性。此外,在樣地Ⅱ中2—5 m尺度上,中齡個體與成過熟齡個體也表現(xiàn)出一定的負(fù)相關(guān)性。

    圖6 樣地Ⅰ各齡級間關(guān)聯(lián)性Fig.6 Associations of different growth stages of Quercus wutaishanica in PlotⅠA: 幼齡階段A; B: 中齡階段B; C: 近熟齡階段C; D: 成過熟齡D

    圖7 樣地Ⅱ各齡級間關(guān)聯(lián)性Fig.7 Associations of different growth stages of Quercus wutaishanica in Plot ⅡA: 幼齡階段A; B: 中齡階段B; C: 近熟齡階段C; D: 成過熟齡D

    4結(jié)論與討論

    遼東櫟是黃龍山林區(qū)植被中重要的建群種之一。在自然環(huán)境條件下,不同生長階段的種群處于不同的林層,受到周圍環(huán)境的影響程度不盡相同,因此其表現(xiàn)出的分布格局也大不相同[35]。根據(jù)林木生長規(guī)律[36]可知:生長條件相似、年齡階段相同的同種個體之間存在激烈的種內(nèi)競爭,包括光照、養(yǎng)分等;生長條件類似、生長階段不同的同種個體,由于競爭原因,高齡個體會對低齡個體產(chǎn)生一定的影響[30]。因此,本文根據(jù)遼東櫟靜態(tài)生命表將遼東櫟劃分為4個不同生長階段,在研究各生長階段空間分布格局的同時,分析其不同生長階段的相互關(guān)聯(lián)性,對了解該物種生長狀況有著重要意義[37]。

    (1)遼東櫟的空間分布格局與尺度有密切關(guān)系:在2塊典型的天然次生林樣地中,遼東櫟在小尺度上均表現(xiàn)為聚集分布,隨尺度增加,其聚集特點逐漸減弱,在中、大尺度上均表現(xiàn)為隨機(jī)分布。兩塊樣地中遼東櫟在小尺度上表現(xiàn)稍有不同:樣地Ⅰ中,遼東櫟個體在0—11 m以內(nèi)為明顯的聚集分布并隨尺度的增大逐漸減弱;而樣地Ⅱ中,在2—6 m尺度上,遼東櫟表現(xiàn)為聚集分布,且隨尺度變化特征不明顯。

    (2)不同樣地相同生長階段其空間分布特征不同:樣地Ⅰ的分布主要表現(xiàn)為隨機(jī)分布,而在樣地Ⅱ中,幼、中齡個體在小尺度上表現(xiàn)為聚集分布,樣地Ⅱ中幼樹發(fā)育要比樣地Ⅰ中有環(huán)境優(yōu)勢;在同一樣地不同生長階段林木的相互關(guān)聯(lián)分析發(fā)現(xiàn),樣地Ⅰ中,低齡個體大多分布于高齡階段樹木樹冠和根系范圍之外,高齡階段樹木對低齡樹木發(fā)育影響程度不顯著,且從圖6中各不同生長階段間相互關(guān)系可以看出,各階段間的相關(guān)性基本表現(xiàn)為不明顯。在樣地Ⅱ中,成、過熟個體對中幼齡個體生長具有負(fù)作用。這與伊力塔[27]、胡秀娟[30]等關(guān)于遼東櫟整體空間分布規(guī)律具有一致性。

    (3)相同樣地不同生長階段的關(guān)聯(lián)性也不盡相同:樣地Ⅰ中各階段之間并沒有太大的聯(lián)系,而樣地Ⅱ中,在小尺度上,幼、中齡個體分別與近、成過熟個體有著一定的負(fù)相關(guān)。從圖1分析可以得出,樣地Ⅱ中高齡階段個體和低齡個體的分布重合較少,因此可以推斷,在高齡個體生長的樹冠和根系覆蓋范圍內(nèi),幼齡個體因受空間、養(yǎng)分、陽光等環(huán)境因子的限制而存活率相對較低,處于高齡個體樹冠及根系覆蓋范圍外的低齡個體,由于空間、養(yǎng)分、陽光等環(huán)境條件較優(yōu),其發(fā)育空間競爭性較低,其個體的發(fā)育狀況相對要好,從而表現(xiàn)為負(fù)相關(guān)性;而在樣地Ⅰ中,由于郁閉度較高(0.79),物種之間的激烈競爭導(dǎo)致幼樹更新地點的選擇有很大的局限性。樣地Ⅰ、Ⅱ中的林分郁閉度分別為0.79、0.48,通過比較可以得出,樣地Ⅰ內(nèi)上層木對陽光的遮擋作用對下層更新影響很大,樣地Ⅱ的條件相對要比樣地Ⅰ更有利于林下更新。針對這種情況,在進(jìn)行森林撫育的時,可以對樣地Ⅰ中上層木進(jìn)行適當(dāng)疏伐,以保證林下植物的生長及林分更新。

    對遼東櫟次生林而言,種群成簇萌生是影響其空間分布的重要因素[38],但是環(huán)境的異質(zhì)性是物種后期發(fā)展的主導(dǎo)因素。種群在環(huán)境影響及其自身生長特性的共同影響下,形成本群落生存發(fā)展的有利結(jié)構(gòu)與其他物種競爭環(huán)境資源,確保自身的正常發(fā)展[18]。研究結(jié)果表明,在兩塊樣地中,低齡個體樹木株數(shù)比高齡個體樹木株數(shù)整體要多,其徑級結(jié)構(gòu)基本呈現(xiàn)倒“J”型,兩塊樣地中遼東櫟種群徑級結(jié)構(gòu)分布較為合理[25]。樣地Ⅰ中喬木株數(shù)密度為2432 株/ hm2,樣地Ⅱ中喬木株數(shù)密度為1018 株/ hm2,在相似立地條件下,樣地Ⅰ競爭環(huán)境要比樣地Ⅱ中苛刻,對其林木生長和分布也產(chǎn)生一定影響。

    種群的空間格局對尺度有著很大的依賴性[27],尺度不同,其表現(xiàn)出來的特征也不盡相同。本研究采用公頃級樣地對黃龍山林區(qū)遼東櫟空間分布格局開展研究。對于遼東櫟次生林而言,在較小尺度上,主要由其種群起源、萌生以及種子的傳播等生物特性決定。而隨著尺度的增大,則由地形地貌、光照、水分等自然條件產(chǎn)生的影響逐漸取代物種自身產(chǎn)生的影響,從而隨自然條件形成不同空間分布格局[19]。研究不同生長階段的空間關(guān)聯(lián)性有利于深入了解其群落的發(fā)展趨勢[39],并對其空間格局的形成過程有一定的了解。而研究種群分布的空間格局有助于充分了解林分的空間分布格局[17],并有針對性的對種群聚集性較大的區(qū)域開展結(jié)構(gòu)化森林經(jīng)營。

    研究結(jié)果表明(圖1),黃龍山林區(qū)遼東櫟純林和油松-遼東櫟-白樺林中,松櫟樺混交林種群結(jié)構(gòu)相對穩(wěn)定,其徑級結(jié)構(gòu)要比遼東櫟純林中更加符合倒“J”型特點。按照群落演替的一般規(guī)律,通過一段時間的競爭、淘汰,遼東櫟幼齡個體在發(fā)育過程中逐漸完成自然稀疏,在此過程中,進(jìn)行適當(dāng)?shù)娜斯じ蓴_如疏伐、撫育等既可以加速該過程,又能改善其內(nèi)部結(jié)構(gòu),使其更加有利于群落的發(fā)育。群落的形成影響因素有許多[4,33],不是簡單的空間格局分析就能完成的。因此,遼東櫟與其他物種(如油松、白樺等)、物種與環(huán)境因子(如土壤、坡度等)之間的相互關(guān)系及影響還有待進(jìn)一步研究,以期更深入地認(rèn)識遼東櫟空間格局形成的起因與過程,更好地為當(dāng)?shù)厣纸?jīng)營提供依據(jù)。

    參考文獻(xiàn)(References):

    [1]方堅, 王孝安, 郭華, 朱志紅. 黃土高原馬欄林區(qū)遼東櫟林種內(nèi)、種間競爭研究. 西北植物學(xué)報, 2007, 27(2): 334- 339.

    [2]Nakashizuka T. Species coexistence in temperate, mixed deciduous forests. Trends in Ecology & Evolution, 2001, 16(4): 205- 210.

    [3]Moloney K A, Wiegand T. Rings, circles, and null-models for point pattern analysis in ecology. Oikos, 2004, 104(2): 209- 229.

    [4]白小軍. 大興安嶺典型林區(qū)主要植物點格局分布研究. 現(xiàn)代農(nóng)業(yè)科技, 2012, (9): 194- 196, 198- 198.

    [5]趙中華, 惠剛盈, 胡艷波, 周紅敏, 張連金. 2種類型闊葉紅松林優(yōu)勢種群空間分布格局及其關(guān)聯(lián)性. 林業(yè)科學(xué)研究, 2011, 24(5): 554- 562.

    [6]Schurr F M, Bossdorf O, Milton S J, Schumacher J. Spatial pattern formation in semi-arid shrubland: a priori predicted versus observed pattern characteristics. Plant Ecology, 2004, 173(2): 271- 282.

    [7]楊曉東, 呂光輝, 張雪梅, 孫麗君, 何書林. 艾比湖濕地自然保護(hù)區(qū)8個喬灌木種群空間分布格局分析. 植物資源與環(huán)境學(xué)報, 2010, 19(4): 37- 42.

    [8]任珩, 趙成章. 高寒退化草地狼毒與賴草種群空間格局及競爭關(guān)系. 生態(tài)學(xué)報, 2013, 33(2): 435- 442.

    [9]繆寧, 劉世榮, 史作民, 喻泓, 劉興良. 川西亞高山紅樺-岷江冷杉林優(yōu)勢種群的空間格局分析. 應(yīng)用生態(tài)學(xué)報, 2009, 20(6): 1263- 1270.

    [10]侯繼華, 黃建輝, 馬克平. 東靈山遼東櫟林主要樹種種群11年動態(tài)變化. 植物生態(tài)學(xué)報, 2004, 28(5): 609- 615.

    [11]康永祥, 岳軍偉, 雷瑞德, 梁宗鎖, 陳永劍. 陜北黃龍山遼東櫟群落優(yōu)勢種群生態(tài)位研究. 西北植物學(xué)報, 2008, 28(3): 574- 581.

    [12]康永祥, 岳軍偉, 張巧明. 黃龍山林區(qū)遼東櫟群落類型劃分及其生物多樣性研究. 西北林學(xué)院學(xué)報, 2007, 22(3): 7- 10.

    [13]張曉輝, 周建云, 何景峰, 盧彥昌, 張文輝. 黃龍山林區(qū)主要森林群落下物種多樣性研究. 西北林學(xué)院學(xué)報, 2007, 22(4): 39- 44.

    [14]郭其強(qiáng), 張文輝, 曹旭平. 基于模糊綜合評判的森林群落穩(wěn)定性評價體系模型構(gòu)建——以黃龍山主要森林群落為例. 林業(yè)科學(xué), 2009, 45(10): 19- 24.

    [15]馮云, 馬克明, 張育新, 祁建. 北京東靈山地區(qū)遼東櫟(Quercusliaotungensis)種群生活史特征與空間分布. 生態(tài)學(xué)雜志, 2009, 28(8): 1443- 1448.

    [16]張育新, 馬克明, 祁建, 馮云, 張潔瑜. 北京東靈山海拔梯度上遼東櫟種群結(jié)構(gòu)和空間分布. 生態(tài)學(xué)報, 2009, 29(6): 2789- 2796.

    [17]袁志良, 王婷, 朱學(xué)靈, 沙迎迎, 葉永忠. 寶天曼落葉闊葉林樣地栓皮櫟種群空間格局. 生物多樣性, 2011, 19(2): 224- 231.

    [18]常新華, 趙秀海, 曾凡勇, 張春雨. 長白山針闊混交林主要樹種空間分布及其環(huán)境解釋. 北京林業(yè)大學(xué)學(xué)報, 2009, 31(1): 7- 12.

    [19]杜志, 亢新剛, 孟京輝, 孔雷, 郭韋韋, 岳剛. 長白山楊樺次生林主要樹種的空間分布格局及其關(guān)聯(lián)性. 東北林業(yè)大學(xué)學(xué)報, 2013, 22(4): 36- 42.

    [20]焦文婧, 郭東罡, 張婕, 上官鐵梁. 靈空山自然保護(hù)區(qū)油松-遼東櫟林建群種關(guān)聯(lián)性. 生態(tài)學(xué)雜志, 2012, 31(12): 3050- 3057.

    [21]岳永杰, 余新曉, 劉彥, 樊登星, 史宇, 王小平, 陳俊崎. 華北落葉松林不同發(fā)育階段種群分布格局研究. 北京林業(yè)大學(xué)學(xué)報, 2008, 31(Z2): 171- 176.

    [22]田麗, 王孝安. 黃土高原馬欄林區(qū)遼東櫟種群動態(tài)分析. 江西農(nóng)業(yè)學(xué)報, 2011, 23(2): 39- 42.

    [23]曹旭平, 郭其強(qiáng), 張文輝. 黃龍山油松林和油松+遼東櫟混交林物種組成及優(yōu)勢種群動態(tài). 西北植物學(xué)報, 2010, 30(5): 1012- 1019.

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

    [25]霍萌萌, 郭東罡, 張婕, 王治明, 上官鐵梁, 李潤強(qiáng), 劉衛(wèi)華. 靈空山油松-遼東櫟林喬木樹種群落學(xué)特征及空間分布格局研究. 生態(tài)學(xué)報, 2014, 34(20): 5925- 5935.

    [26]王巍, 劉燦然, 馬克平, 于順利. 東靈山兩個落葉闊葉林中遼東櫟種群結(jié)構(gòu)和動態(tài). 植物學(xué)報, 1999, 41(4): 425- 432.

    [27]伊力塔, 韓海榮, 程小琴, 康峰峰, 張志杰. 靈空山林區(qū)遼東櫟(Quercusliaotungensis)種群空間分布格局. 生態(tài)學(xué)報, 2008, 28(7): 3254- 3261.

    [28]孟憲宇. 測樹學(xué). 北京: 中國林業(yè)出版社, 2006: 63.

    [29]牛麗麗, 余新曉, 劉淑燕, 劉彥, 岳永杰, 王小平, 陳俊崎. 北京松山自然保護(hù)區(qū)油松種群生活史及其空間分布格局. 北京林業(yè)大學(xué)學(xué)報, 2008, 30(Z2): 17- 21.

    [30]胡秀娟, 程積民, 楊曉梅, 范文娟, 孟蕾, 韓娟娟. 黃土高原子午嶺林區(qū)遼東櫟種群分布格局及群落特征研究. 西北林學(xué)院學(xué)報, 2010, 25(5): 1- 6.

    [31]杜志, 亢新剛, 包昱君, 楊鑫霞. 長白山云冷杉林不同演替階段的樹種空間分布格局及其關(guān)聯(lián)性. 北京林業(yè)大學(xué)學(xué)報, 2012, 34(2): 14- 19.

    [32]張莉, 張小平, 陸暢, 李曉紅. 安徽瑯琊山青檀種群空間格局. 林業(yè)科學(xué), 2012, 48(2): 9- 15.

    [33]郭垚鑫, 康冰, 李剛, 王得祥, 楊改河, 王大偉. 小隴山紅樺次生林物種組成與立木的點格局分析. 應(yīng)用生態(tài)學(xué)報, 2011, 22(10): 2574- 2580.

    [34]馬莘, 夏穎, 馬慶慶, 宋會興. 四川牡丹群落物種生態(tài)位及空間關(guān)聯(lián)度分析. 重慶師范大學(xué)學(xué)報: 自然科學(xué)版, 2011, 28(4): 57- 61.

    [35]伊力塔, 豪樹奇, 徐涌, 康峰峰. 遼東櫟群落特征研究及其進(jìn)展. 遼寧林業(yè)科技, 2010, (5): 25- 30, 39- 39.

    [36]龔直文, 亢新剛, 顧麗, 楊華. 長白山楊樺次生林生長過程與演替動向分析. 林業(yè)科學(xué)研究, 2009, 22(3): 74- 79.

    [37]宋于洋, 李園園, 張文輝. 梭梭種群不同發(fā)育階段的空間格局與關(guān)聯(lián)性分析. 生態(tài)學(xué)報, 2010, 30(16): 4317- 4327.

    [38]孫書存, 陳靈芝. 遼東櫟植冠的構(gòu)型分析. 植物生態(tài)學(xué)報, 1999, 23(5): 433- 440.

    [39]龔直文, 顧麗, 亢新剛, 楊華, 趙浩彥. 長白山森林次生演替過程中林木空間格局研究. 北京林業(yè)大學(xué)學(xué)報, 2010, 32(2): 92- 99.

    Spatial pattern ofQuercuswutaishanicain natural secondary forest of Huanglong Mountain

    WANG Daoliang1, LI Weizhong1,*, CAO Zhu1, WANG Qing1, WANG Guangru2

    1CollegeofForestry,NorthwestA&FUniversity,Yangling712100,China2ForestBureauofHuanglongshaninYan′anCity,Huanglong715700,China

    Abstract:Spatial pattern analysis is an important method for studying population characteristics, interspecies relationships, as well as the relationships between population and environment. In the present paper, population characteristics of Quercus wutaishanica in natural secondary forest of Huanglong Mountain were investigated. One natural secondary forest plot (with area of 1 hm2) was selected and set up in each Q. wutaishanica-Pinus tabuliformis-Betula platyphylla natural secondary mixed forest and Q. wutaishanica natural secondary pure forest. The data and coordinates of each tree in the two plots were measured and recorded following the standards for plot investigation. Using the method of diameter class (instead of the age), we converted the classification of different growth stages of Q. wutaishanica into size classes for diameter at breast height (DBH) according to the static life table. Thus, Q. wutaishanica was divided into four different growth stages: the immature stage A, 3 cm ≤ DBH < 10 cm; mid-mature stage B, 10 cm ≤ DBH < 20 cm; near-mature stage C, 20 cm ≤ DBH < 30cm; mature and over-mature stage D, DBH ≥ 30 cm. The spatial pattern of Q. wutaishanica was analyzed through O-ring statistical means for each growth stage by using single variable method and the interconnectedness between the growth stages by bivariate-analysis method. The results showed the following: (1) The diameter structure of Q. wutaishanica in the two plots was in line with inverted "J" type, indicating that the natural regeneration in the two plots is better from the perspective of diameter class distribution. (2) There was a close relationship between the spatial pattern and the scale of Q. wutaishanica. It showed random distribution in the two selected plots at a large and medium scale and aggregated distribution at a small scale, especially at the scale of 0—6 m in plot I. (3) The spatial patterns differed at different growth stages, young trees exhibited aggregated distribution at a small scale, especially at the scale of 0—11 m in plot II. (4) The correlation between different growth stages was not identical in the plots. There was no significant correlation between different growth stages within plot I. The immature and mid-mature individuals showed a significant positive correlation at a small scale of plot II, whereas the near-mature, mature, and over-mature individuals showed a significant negative correlation at a small scale. The mid-mature, mature, and over-mature individuals were negatively correlated at the 2—5 m scale. Q. wutaishanica stand in the natural secondary forest on Huanglong Mountain was one of the important constituent species. The study of spatial distribution at different growth stages of individual organisms is critical to capture the growth situation of Q. wutaishanica in Huanglong Mountain. Spatial correlation of different growth stages is vital to understand the development trend of the community as well as to provide an important basis for the development of structured forest management for greater population density.

    Key Words:natural secondary forest; spatial distribution pattern; spatial correlation; O-ring statistical analysis; Quercus wutaishanica

    基金項目:國家自然科學(xué)基金項目(31170587)

    收稿日期:2014- 06- 14; 網(wǎng)絡(luò)出版日期:2015- 08- 26

    *通訊作者

    Corresponding author.E-mail: wzhli6465@163.com

    DOI:10.5846/stxb201406141233

    王道亮,李衛(wèi)忠,曹鑄,王慶,王廣儒.黃龍山天然次生林遼東櫟種群空間格局.生態(tài)學(xué)報,2016,36(9):2677- 2685.

    Wang D L, Li W Z, Cao Z, Wang Q, Wang G R.Spatial pattern ofQuercuswutaishanicain natural secondary forest of Huanglong Mountain.Acta Ecologica Sinica,2016,36(9):2677- 2685.

    一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 波多野结衣av一区二区av| 伊人久久大香线蕉亚洲五| 成年动漫av网址| 女性生殖器流出的白浆| 久久99精品国语久久久| 国产亚洲午夜精品一区二区久久| av电影中文网址| 飞空精品影院首页| 欧美老熟妇乱子伦牲交| 最新的欧美精品一区二区| av福利片在线| av国产精品久久久久影院| 免费黄色在线免费观看| 亚洲精品一区蜜桃| 欧美最新免费一区二区三区| 少妇人妻 视频| 国产一区二区 视频在线| 美女国产高潮福利片在线看| 美女午夜性视频免费| 国产成人精品在线电影| av福利片在线| 日本av手机在线免费观看| 人妻少妇偷人精品九色| avwww免费| 99国产综合亚洲精品| 99香蕉大伊视频| 久久午夜综合久久蜜桃| 一边摸一边做爽爽视频免费| 男女下面进入的视频免费午夜 | 久久久久久人人人人人| 亚洲久久久国产精品| 国产精品 国内视频| 久久精品国产99精品国产亚洲性色 | 免费高清在线观看日韩| 日本wwww免费看| 夜夜躁狠狠躁天天躁| 国产精品免费视频内射| 少妇被粗大的猛进出69影院| a级毛片在线看网站| aaaaa片日本免费| 国产99久久九九免费精品| 国产亚洲av高清不卡| 久久亚洲精品不卡| 久久久久久久久久久久大奶| 伦理电影免费视频| 一级a爱视频在线免费观看| 亚洲av成人不卡在线观看播放网| 久久久国产欧美日韩av| 久久久国产欧美日韩av| 成人三级黄色视频| av视频免费观看在线观看| 欧美老熟妇乱子伦牲交| av有码第一页| 少妇裸体淫交视频免费看高清 | 久久中文字幕一级| 满18在线观看网站| 欧美成狂野欧美在线观看| 亚洲人成电影免费在线| 久久精品国产99精品国产亚洲性色 | 性色av乱码一区二区三区2| 成人手机av| tocl精华| 两性夫妻黄色片| 女人被狂操c到高潮| 国产精品 欧美亚洲| 国产在线精品亚洲第一网站| 国产欧美日韩一区二区精品| 亚洲一区中文字幕在线| 亚洲欧美日韩另类电影网站| 嫩草影视91久久| 黄片播放在线免费| 精品第一国产精品| 我的亚洲天堂| 黑人欧美特级aaaaaa片| 日日夜夜操网爽| 久久久久精品国产欧美久久久| 精品久久久久久电影网| 欧美一级毛片孕妇| 久久人人爽av亚洲精品天堂| 黄片播放在线免费| 国产亚洲精品综合一区在线观看 | 一级毛片高清免费大全| 国产亚洲精品一区二区www| 亚洲精品在线美女| 亚洲国产精品sss在线观看 | 国产精品自产拍在线观看55亚洲| 欧美在线一区亚洲| 国产欧美日韩一区二区精品| 亚洲欧洲精品一区二区精品久久久| 丝袜美足系列| 亚洲 欧美一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 热re99久久国产66热| 欧美国产精品va在线观看不卡| 日本欧美视频一区| 久久精品亚洲精品国产色婷小说| 在线国产一区二区在线| 在线观看www视频免费| 亚洲 国产 在线| 亚洲男人天堂网一区| 一级,二级,三级黄色视频| 在线观看日韩欧美| 欧美中文日本在线观看视频| 91大片在线观看| 两个人看的免费小视频| 日韩国内少妇激情av| 久热爱精品视频在线9| 波多野结衣av一区二区av| 国产一卡二卡三卡精品| 五月开心婷婷网| 中亚洲国语对白在线视频| 91国产中文字幕| 午夜免费成人在线视频| 一级毛片精品| 9热在线视频观看99| 欧美日韩av久久| 老司机靠b影院| 巨乳人妻的诱惑在线观看| 欧美激情 高清一区二区三区| 99国产综合亚洲精品| 男人的好看免费观看在线视频 | 十分钟在线观看高清视频www| 国产精品一区二区免费欧美| 超碰成人久久| 999精品在线视频| 精品一区二区三卡| 久久久国产精品麻豆| 一边摸一边抽搐一进一小说| 两人在一起打扑克的视频| 丰满迷人的少妇在线观看| 搡老熟女国产l中国老女人| 日本黄色视频三级网站网址| 欧美av亚洲av综合av国产av| 日韩三级视频一区二区三区| 人成视频在线观看免费观看| 99热国产这里只有精品6| 日本三级黄在线观看| 国产高清videossex| 91大片在线观看| av福利片在线| 美女午夜性视频免费| 波多野结衣av一区二区av| 亚洲中文字幕日韩| 黄色丝袜av网址大全| 色播在线永久视频| 亚洲欧美激情在线| 99re在线观看精品视频| 久久精品人人爽人人爽视色| 又大又爽又粗| 久久久久国产一级毛片高清牌| 真人做人爱边吃奶动态| 午夜两性在线视频| 搡老岳熟女国产| 美女高潮喷水抽搐中文字幕| 大型av网站在线播放| 老司机午夜十八禁免费视频| 怎么达到女性高潮| 亚洲精品久久午夜乱码| 久久久久久久久免费视频了| 在线观看午夜福利视频| 国产精品久久视频播放| av福利片在线| 亚洲国产精品999在线| 91成人精品电影| 欧美日本亚洲视频在线播放| 国产av在哪里看| 岛国视频午夜一区免费看| 亚洲欧美一区二区三区黑人| 制服诱惑二区| 成熟少妇高潮喷水视频| 黄色毛片三级朝国网站| 亚洲五月色婷婷综合| 午夜影院日韩av| 在线天堂中文资源库| 欧美成人性av电影在线观看| 黄色片一级片一级黄色片| 久久午夜综合久久蜜桃| 激情在线观看视频在线高清| 午夜影院日韩av| 精品久久久久久久毛片微露脸| 69精品国产乱码久久久| 日本撒尿小便嘘嘘汇集6| 黄片播放在线免费| 国产免费男女视频| 亚洲第一av免费看| 欧美人与性动交α欧美精品济南到| 高清在线国产一区| 成年版毛片免费区| av天堂在线播放| 宅男免费午夜| 国产一区二区三区在线臀色熟女 | 高清av免费在线| 午夜91福利影院| 韩国精品一区二区三区| 老汉色av国产亚洲站长工具| 真人一进一出gif抽搐免费| 欧美一级毛片孕妇| 脱女人内裤的视频| 一边摸一边做爽爽视频免费| 99re在线观看精品视频| 国产av精品麻豆| 久久香蕉激情| 一边摸一边抽搐一进一小说| 欧美成狂野欧美在线观看| 男女高潮啪啪啪动态图| 老司机深夜福利视频在线观看| 欧美日韩黄片免| 宅男免费午夜| 久久九九热精品免费| svipshipincom国产片| 国产熟女xx| 嫩草影视91久久| 一二三四社区在线视频社区8| 少妇裸体淫交视频免费看高清 | 久久久国产精品麻豆| 两性夫妻黄色片| www.熟女人妻精品国产| 久久久久国产精品人妻aⅴ院| 三上悠亚av全集在线观看| 亚洲成人国产一区在线观看| 人成视频在线观看免费观看| 少妇的丰满在线观看| 国产亚洲精品综合一区在线观看 | 国产精品香港三级国产av潘金莲| 久久中文字幕一级| 日本一区二区免费在线视频| 黄色毛片三级朝国网站| www.自偷自拍.com| 亚洲av日韩精品久久久久久密| 国产av又大| 欧美中文综合在线视频| www日本在线高清视频| 老熟妇仑乱视频hdxx| 免费av中文字幕在线| 亚洲欧美日韩高清在线视频| 欧美 亚洲 国产 日韩一| 国产真人三级小视频在线观看| 久久久久久久精品吃奶| 亚洲avbb在线观看| a级毛片黄视频| 午夜视频精品福利| 亚洲精品美女久久av网站| 99在线视频只有这里精品首页| 亚洲欧美日韩高清在线视频| 国产99白浆流出| 人人妻人人爽人人添夜夜欢视频| 色精品久久人妻99蜜桃| 色综合欧美亚洲国产小说| 中文字幕精品免费在线观看视频| 久久国产乱子伦精品免费另类| 国内毛片毛片毛片毛片毛片| 十分钟在线观看高清视频www| 亚洲欧美一区二区三区久久| 最好的美女福利视频网| 欧美日韩亚洲综合一区二区三区_| 在线av久久热| 校园春色视频在线观看| 精品久久久久久久久久免费视频 | 啦啦啦在线免费观看视频4| 亚洲午夜精品一区,二区,三区| 成人国产一区最新在线观看| 亚洲久久久国产精品| 亚洲精华国产精华精| 老司机福利观看| 久久久久国产精品人妻aⅴ院| 亚洲国产看品久久| 亚洲一码二码三码区别大吗| 99精品久久久久人妻精品| 久久影院123| 亚洲一区中文字幕在线| 一a级毛片在线观看| 国产精品爽爽va在线观看网站 | 国产99久久九九免费精品| 纯流量卡能插随身wifi吗| 国产亚洲欧美在线一区二区| 午夜老司机福利片| 亚洲午夜理论影院| 午夜免费成人在线视频| 成在线人永久免费视频| 免费不卡黄色视频| 亚洲va日本ⅴa欧美va伊人久久| 免费少妇av软件| 中亚洲国语对白在线视频| 亚洲国产精品sss在线观看 | 老司机在亚洲福利影院| 美女高潮喷水抽搐中文字幕| 男人操女人黄网站| 黄色怎么调成土黄色| 丝袜在线中文字幕| 欧美日韩亚洲综合一区二区三区_| 十分钟在线观看高清视频www| 一区二区三区国产精品乱码| 黄网站色视频无遮挡免费观看| 亚洲成人国产一区在线观看| 香蕉国产在线看| 看免费av毛片| 悠悠久久av| 露出奶头的视频| 亚洲成a人片在线一区二区| 国产成+人综合+亚洲专区| 国产在线观看jvid| 久久亚洲真实| 国产aⅴ精品一区二区三区波| 悠悠久久av| 黑丝袜美女国产一区| 国产av精品麻豆| 久9热在线精品视频| 免费一级毛片在线播放高清视频 | 又紧又爽又黄一区二区| 丰满的人妻完整版| 午夜免费激情av| 脱女人内裤的视频| 乱人伦中国视频| 亚洲av第一区精品v没综合| 熟女少妇亚洲综合色aaa.| 久久国产精品影院| 午夜两性在线视频| 在线天堂中文资源库| 日韩一卡2卡3卡4卡2021年| 国产伦一二天堂av在线观看| 欧美av亚洲av综合av国产av| 这个男人来自地球电影免费观看| 中文字幕精品免费在线观看视频| 叶爱在线成人免费视频播放| 黑人操中国人逼视频| 长腿黑丝高跟| 黄色怎么调成土黄色| 日韩欧美一区二区三区在线观看| 欧美另类亚洲清纯唯美| 在线视频色国产色| 两人在一起打扑克的视频| 亚洲精品国产区一区二| 中文欧美无线码| av免费在线观看网站| 男女之事视频高清在线观看| 高清欧美精品videossex| 亚洲中文av在线| 亚洲情色 制服丝袜| 午夜福利一区二区在线看| 中文字幕色久视频| 午夜影院日韩av| 国产成人av教育| 老司机靠b影院| 高清毛片免费观看视频网站 | 欧美激情极品国产一区二区三区| 高潮久久久久久久久久久不卡| 日韩欧美一区二区三区在线观看| 色尼玛亚洲综合影院| 巨乳人妻的诱惑在线观看| 97超级碰碰碰精品色视频在线观看| 老熟妇仑乱视频hdxx| 青草久久国产| 超碰97精品在线观看| 国产黄色免费在线视频| 九色亚洲精品在线播放| 久久天躁狠狠躁夜夜2o2o| 亚洲一区中文字幕在线| 精品福利永久在线观看| 欧美精品亚洲一区二区| 久久久久久久久久久久大奶| 九色亚洲精品在线播放| 麻豆av在线久日| 欧美日韩av久久| 亚洲av第一区精品v没综合| 欧美日韩亚洲国产一区二区在线观看| 村上凉子中文字幕在线| 久久精品国产99精品国产亚洲性色 | www国产在线视频色| 久久久久国产精品人妻aⅴ院| 午夜福利,免费看| 无人区码免费观看不卡| www日本在线高清视频| 欧美激情 高清一区二区三区| 别揉我奶头~嗯~啊~动态视频| 超碰成人久久| 91精品三级在线观看| 欧美丝袜亚洲另类 | 亚洲九九香蕉| 国产成人啪精品午夜网站| 999久久久国产精品视频| 久久久久久久精品吃奶| 精品一区二区三区av网在线观看| 亚洲精品一区av在线观看| 超色免费av| 老司机深夜福利视频在线观看| www.999成人在线观看| 免费在线观看黄色视频的| 亚洲一区中文字幕在线| 国产精品 国内视频| 曰老女人黄片| 露出奶头的视频| 免费av中文字幕在线| 自线自在国产av| 亚洲av第一区精品v没综合| 亚洲精品久久午夜乱码| 国产av在哪里看| 亚洲精华国产精华精| 可以在线观看毛片的网站| 午夜a级毛片| 一级片免费观看大全| 50天的宝宝边吃奶边哭怎么回事| 男女午夜视频在线观看| 免费看十八禁软件| 欧美乱码精品一区二区三区| 亚洲精品久久成人aⅴ小说| 久久精品亚洲av国产电影网| 久久精品亚洲精品国产色婷小说| 欧美老熟妇乱子伦牲交| 久久久久久亚洲精品国产蜜桃av| 午夜两性在线视频| 免费观看人在逋| 亚洲欧美日韩高清在线视频| 巨乳人妻的诱惑在线观看| 欧美日韩福利视频一区二区| 国产精品99久久99久久久不卡| 亚洲av电影在线进入| 男女高潮啪啪啪动态图| 久久久国产欧美日韩av| 亚洲色图 男人天堂 中文字幕| 成年版毛片免费区| 欧美日本中文国产一区发布| 亚洲欧美精品综合久久99| 亚洲精品中文字幕在线视频| 精品免费久久久久久久清纯| 免费观看精品视频网站| 亚洲精品久久成人aⅴ小说| 国产麻豆69| 亚洲三区欧美一区| 亚洲aⅴ乱码一区二区在线播放 | 国产精品亚洲一级av第二区| 丝袜美腿诱惑在线| 怎么达到女性高潮| 日韩欧美一区视频在线观看| 免费不卡黄色视频| 欧美亚洲日本最大视频资源| 男女下面插进去视频免费观看| 啦啦啦在线免费观看视频4| av视频免费观看在线观看| 亚洲免费av在线视频| 午夜精品国产一区二区电影| 亚洲成a人片在线一区二区| 亚洲国产毛片av蜜桃av| 久久伊人香网站| 1024香蕉在线观看| 看免费av毛片| 午夜免费观看网址| www.精华液| 久久久久亚洲av毛片大全| 国产色视频综合| 曰老女人黄片| 国内毛片毛片毛片毛片毛片| 久久国产精品人妻蜜桃| 在线视频色国产色| 波多野结衣一区麻豆| 亚洲va日本ⅴa欧美va伊人久久| 国产三级黄色录像| 国产av一区在线观看免费| 91成人精品电影| 免费高清视频大片| 午夜福利在线观看吧| 国产激情欧美一区二区| 他把我摸到了高潮在线观看| 欧美中文日本在线观看视频| 天天躁夜夜躁狠狠躁躁| a在线观看视频网站| 亚洲一区中文字幕在线| 俄罗斯特黄特色一大片| 久久青草综合色| 久久精品国产清高在天天线| 老司机福利观看| 男女午夜视频在线观看| 欧美日韩一级在线毛片| 亚洲中文字幕日韩| 久久天躁狠狠躁夜夜2o2o| 变态另类成人亚洲欧美熟女 | 欧美色视频一区免费| 国产精品亚洲一级av第二区| 最近最新中文字幕大全免费视频| 亚洲午夜理论影院| 欧美日韩中文字幕国产精品一区二区三区 | 日韩 欧美 亚洲 中文字幕| 色婷婷av一区二区三区视频| 欧美亚洲日本最大视频资源| 在线观看午夜福利视频| 亚洲国产精品合色在线| 久久人妻福利社区极品人妻图片| 真人做人爱边吃奶动态| 亚洲久久久国产精品| 黄色丝袜av网址大全| 精品一品国产午夜福利视频| 亚洲第一欧美日韩一区二区三区| 亚洲 国产 在线| 一区在线观看完整版| 国产亚洲av高清不卡| 一个人观看的视频www高清免费观看 | 亚洲成人精品中文字幕电影 | 国产1区2区3区精品| 亚洲男人天堂网一区| 制服人妻中文乱码| 国产精品久久久人人做人人爽| 国产激情久久老熟女| 亚洲国产精品sss在线观看 | 免费一级毛片在线播放高清视频 | 女生性感内裤真人,穿戴方法视频| 国产日韩一区二区三区精品不卡| 久久精品国产亚洲av香蕉五月| 久久久国产一区二区| 巨乳人妻的诱惑在线观看| 成人国语在线视频| 神马国产精品三级电影在线观看 | 成人三级做爰电影| 亚洲欧美日韩无卡精品| 人人妻人人澡人人看| 亚洲全国av大片| 欧美黄色淫秽网站| 欧美午夜高清在线| 国产成人免费无遮挡视频| 亚洲av第一区精品v没综合| 中文亚洲av片在线观看爽| 99久久国产精品久久久| 欧美在线黄色| 黄色怎么调成土黄色| 啦啦啦 在线观看视频| 免费在线观看日本一区| 亚洲熟妇中文字幕五十中出 | 80岁老熟妇乱子伦牲交| 亚洲avbb在线观看| 在线观看一区二区三区| 一级片'在线观看视频| 久久国产精品男人的天堂亚洲| 成人国语在线视频| 亚洲成人免费av在线播放| 女人精品久久久久毛片| 日本五十路高清| 亚洲成a人片在线一区二区| 亚洲av成人不卡在线观看播放网| 99久久国产精品久久久| 天堂俺去俺来也www色官网| 又大又爽又粗| 亚洲熟女毛片儿| 国产精品成人在线| 激情视频va一区二区三区| 青草久久国产| 亚洲一区二区三区色噜噜 | 国产精品一区二区免费欧美| 成人免费观看视频高清| av免费在线观看网站| 日本一区二区免费在线视频| 久久亚洲精品不卡| 精品高清国产在线一区| 亚洲五月天丁香| 好男人电影高清在线观看| 国产精品秋霞免费鲁丝片| 91字幕亚洲| 夜夜夜夜夜久久久久| av中文乱码字幕在线| 国产精品亚洲av一区麻豆| 国产精品久久久av美女十八| 国产97色在线日韩免费| 久久久国产欧美日韩av| 女人被狂操c到高潮| 在线视频色国产色| 人妻久久中文字幕网| 男女下面进入的视频免费午夜 | 好男人电影高清在线观看| 高清av免费在线| 国产一区二区在线av高清观看| 色老头精品视频在线观看| 99精品久久久久人妻精品| 日韩 欧美 亚洲 中文字幕| 国产精品电影一区二区三区| 亚洲黑人精品在线| e午夜精品久久久久久久| 天天躁夜夜躁狠狠躁躁| 操出白浆在线播放| 欧美日本亚洲视频在线播放| 中文字幕色久视频| 久久精品aⅴ一区二区三区四区| 日韩三级视频一区二区三区| 久久久久久久久免费视频了| 一区二区三区激情视频| 久久久久久久久久久久大奶| 亚洲三区欧美一区| 99在线视频只有这里精品首页| 日韩免费av在线播放| 精品久久久久久电影网| 乱人伦中国视频| 国产成人欧美在线观看| 国产又色又爽无遮挡免费看| 黑人猛操日本美女一级片| 在线免费观看的www视频| 免费观看精品视频网站| 亚洲av日韩精品久久久久久密| 黄片播放在线免费| 国产真人三级小视频在线观看| svipshipincom国产片| 国产欧美日韩精品亚洲av| 91在线观看av| 亚洲精品一二三| 亚洲色图综合在线观看| 在线免费观看的www视频| 国产亚洲精品综合一区在线观看 | 国产黄色免费在线视频| 麻豆久久精品国产亚洲av | 乱人伦中国视频| 久久久久国产一级毛片高清牌| 在线永久观看黄色视频| 一级毛片女人18水好多| 国产av精品麻豆| 老司机在亚洲福利影院| 香蕉国产在线看| 精品一区二区三区av网在线观看| 长腿黑丝高跟| 香蕉久久夜色|