• <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.

    国产男人的电影天堂91| 久久久久精品国产欧美久久久| 免费av观看视频| 国产午夜福利久久久久久| 亚洲成人中文字幕在线播放| 日韩强制内射视频| 能在线免费观看的黄片| 十八禁国产超污无遮挡网站| av福利片在线观看| www日本黄色视频网| 久久精品综合一区二区三区| 最近2019中文字幕mv第一页| 成人欧美大片| 国产女主播在线喷水免费视频网站 | 波野结衣二区三区在线| 精华霜和精华液先用哪个| 91久久精品国产一区二区三区| 在线免费十八禁| 久久久a久久爽久久v久久| 日本熟妇午夜| 五月伊人婷婷丁香| 69av精品久久久久久| 国产69精品久久久久777片| 99久久精品一区二区三区| 不卡一级毛片| 亚洲精品粉嫩美女一区| 国产人妻一区二区三区在| 又黄又爽又刺激的免费视频.| .国产精品久久| 色噜噜av男人的天堂激情| 午夜亚洲福利在线播放| 欧美极品一区二区三区四区| 亚洲欧美日韩东京热| 2021天堂中文幕一二区在线观| 国产一区二区在线av高清观看| 欧美色视频一区免费| av在线老鸭窝| 国产女主播在线喷水免费视频网站 | 啦啦啦啦在线视频资源| 日本熟妇午夜| 国产蜜桃级精品一区二区三区| 亚洲内射少妇av| 欧美性猛交黑人性爽| av在线播放精品| 国产麻豆成人av免费视频| 国产激情偷乱视频一区二区| 亚洲精品在线观看二区| 激情 狠狠 欧美| 国产麻豆成人av免费视频| 99久久精品一区二区三区| 22中文网久久字幕| 欧美性猛交╳xxx乱大交人| 热99re8久久精品国产| 日韩一区二区视频免费看| 两个人的视频大全免费| 长腿黑丝高跟| 少妇熟女欧美另类| 伊人久久精品亚洲午夜| 搡女人真爽免费视频火全软件 | 18+在线观看网站| 在线观看av片永久免费下载| 美女xxoo啪啪120秒动态图| 免费在线观看成人毛片| av在线亚洲专区| 最好的美女福利视频网| 欧美激情在线99| 特大巨黑吊av在线直播| 黄色配什么色好看| 香蕉av资源在线| 成人无遮挡网站| 亚洲国产精品国产精品| 日本-黄色视频高清免费观看| 香蕉av资源在线| 男人的好看免费观看在线视频| 国产大屁股一区二区在线视频| 女同久久另类99精品国产91| 我的女老师完整版在线观看| 成人av在线播放网站| 99精品在免费线老司机午夜| 午夜福利18| 在线观看美女被高潮喷水网站| 亚洲精品乱码久久久v下载方式| 青春草视频在线免费观看| 免费观看在线日韩| 国产一区二区在线av高清观看| 最近的中文字幕免费完整| 精品久久久久久久人妻蜜臀av| 精品一区二区三区视频在线观看免费| 亚洲自偷自拍三级| 少妇熟女欧美另类| 国产黄片美女视频| 久久精品影院6| 亚洲欧美中文字幕日韩二区| 免费人成在线观看视频色| 国产精品综合久久久久久久免费| 日韩强制内射视频| 亚洲国产精品久久男人天堂| 99久久九九国产精品国产免费| 男女做爰动态图高潮gif福利片| 亚洲精品影视一区二区三区av| 亚洲精品乱码久久久v下载方式| 日本成人三级电影网站| 国产成年人精品一区二区| 亚洲av免费高清在线观看| 啦啦啦韩国在线观看视频| 日本免费一区二区三区高清不卡| 午夜亚洲福利在线播放| 99热这里只有是精品50| av天堂在线播放| 一级毛片久久久久久久久女| 乱码一卡2卡4卡精品| 国产午夜福利久久久久久| 少妇的逼好多水| 欧美xxxx性猛交bbbb| 日韩在线高清观看一区二区三区| 热99在线观看视频| 亚洲国产精品成人久久小说 | 亚洲精品亚洲一区二区| 我要看日韩黄色一级片| av黄色大香蕉| 人妻夜夜爽99麻豆av| 综合色av麻豆| 亚洲国产精品久久男人天堂| 大型黄色视频在线免费观看| 两个人视频免费观看高清| 欧洲精品卡2卡3卡4卡5卡区| 国产高清视频在线播放一区| 天美传媒精品一区二区| 免费黄网站久久成人精品| 国产精品久久久久久久电影| 久久久精品欧美日韩精品| 成人二区视频| 久久精品91蜜桃| 久久久久久久午夜电影| 午夜精品国产一区二区电影 | 国产爱豆传媒在线观看| 综合色丁香网| av在线老鸭窝| 97碰自拍视频| 夜夜爽天天搞| 看黄色毛片网站| 国产精品无大码| 亚洲婷婷狠狠爱综合网| 99久久精品一区二区三区| 国产午夜精品久久久久久一区二区三区 | 最新在线观看一区二区三区| 成人亚洲欧美一区二区av| 精品午夜福利在线看| 久99久视频精品免费| 91在线观看av| 久久久久久久久久黄片| 国产精品一区www在线观看| 亚洲欧美成人精品一区二区| 精品一区二区免费观看| 精品一区二区三区视频在线| 日本熟妇午夜| 国产国拍精品亚洲av在线观看| 国产精品一二三区在线看| 国产精品乱码一区二三区的特点| 哪里可以看免费的av片| 极品教师在线视频| 国产精品一区二区性色av| 国产精品爽爽va在线观看网站| 亚洲精品456在线播放app| 九九在线视频观看精品| 亚洲性夜色夜夜综合| 日本三级黄在线观看| 国产成人精品久久久久久| a级毛片免费高清观看在线播放| 欧美3d第一页| 国产午夜福利久久久久久| 黄色视频,在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲第一区二区三区不卡| 我要看日韩黄色一级片| 亚洲在线自拍视频| 美女免费视频网站| 国产精华一区二区三区| 伦精品一区二区三区| 久久精品国产亚洲av香蕉五月| 在现免费观看毛片| 日日摸夜夜添夜夜添av毛片| 97超视频在线观看视频| 免费看日本二区| 在线观看免费视频日本深夜| 俄罗斯特黄特色一大片| 国内精品久久久久精免费| 99国产极品粉嫩在线观看| 亚洲高清免费不卡视频| 国产在线男女| 日日撸夜夜添| 晚上一个人看的免费电影| 夜夜看夜夜爽夜夜摸| 乱人视频在线观看| 国内久久婷婷六月综合欲色啪| 身体一侧抽搐| 精品免费久久久久久久清纯| 国产 一区 欧美 日韩| 男人舔奶头视频| 欧美zozozo另类| 国产精品一区二区免费欧美| av国产免费在线观看| 日韩亚洲欧美综合| 伦理电影大哥的女人| 天堂√8在线中文| 老司机福利观看| 国产午夜福利久久久久久| 欧美日韩乱码在线| 99精品在免费线老司机午夜| 久久久久久久久大av| 午夜激情欧美在线| 国产午夜福利久久久久久| 欧美日韩乱码在线| 小蜜桃在线观看免费完整版高清| 99久久中文字幕三级久久日本| 99热这里只有是精品在线观看| 国产精品亚洲一级av第二区| 丝袜美腿在线中文| 日韩在线高清观看一区二区三区| 国产亚洲91精品色在线| 亚洲国产精品久久男人天堂| 亚洲精品日韩av片在线观看| 中文字幕精品亚洲无线码一区| h日本视频在线播放| 色哟哟·www| 校园春色视频在线观看| 国产三级在线视频| 亚洲国产色片| 久久精品91蜜桃| 日本免费a在线| 亚洲av成人精品一区久久| 一区二区三区高清视频在线| 欧美在线一区亚洲| 日韩成人伦理影院| 国国产精品蜜臀av免费| 国产大屁股一区二区在线视频| 亚洲美女黄片视频| 久久久色成人| 成人亚洲欧美一区二区av| 成人漫画全彩无遮挡| 国产精品一区二区免费欧美| 国产蜜桃级精品一区二区三区| 国产日本99.免费观看| 在线观看免费视频日本深夜| 免费搜索国产男女视频| 97超碰精品成人国产| 久久久久久伊人网av| 麻豆精品久久久久久蜜桃| 欧美最新免费一区二区三区| 男女之事视频高清在线观看| 久久亚洲国产成人精品v| 久久中文看片网| 极品教师在线视频| 亚洲乱码一区二区免费版| 免费观看的影片在线观看| 国产在线男女| 少妇的逼好多水| 亚洲图色成人| 久久精品久久久久久噜噜老黄 | 午夜福利视频1000在线观看| 亚洲专区国产一区二区| 色播亚洲综合网| 赤兔流量卡办理| 国产成人91sexporn| 欧美色欧美亚洲另类二区| 亚洲国产高清在线一区二区三| 一区二区三区四区激情视频 | 神马国产精品三级电影在线观看| 最后的刺客免费高清国语| 国产成人一区二区在线| 国产一区二区亚洲精品在线观看| 精品久久久久久久久久免费视频| 搡女人真爽免费视频火全软件 | 精品人妻偷拍中文字幕| 亚洲欧美日韩高清专用| 成人精品一区二区免费| 草草在线视频免费看| 在线看三级毛片| 午夜福利在线观看免费完整高清在 | 性色avwww在线观看| 免费看光身美女| 九色成人免费人妻av| 久久精品国产亚洲av天美| 国产精品无大码| 狂野欧美白嫩少妇大欣赏| 亚洲第一区二区三区不卡| 毛片一级片免费看久久久久| 老熟妇乱子伦视频在线观看| 天美传媒精品一区二区| 亚洲,欧美,日韩| 97在线视频观看| 老女人水多毛片| 国产不卡一卡二| 韩国av在线不卡| 国产女主播在线喷水免费视频网站 | 国产精品99久久久久久久久| 色综合站精品国产| 亚洲av中文字字幕乱码综合| 天堂动漫精品| 国产色爽女视频免费观看| 一区福利在线观看| 久久久国产成人精品二区| 亚洲精品日韩在线中文字幕 | 欧美潮喷喷水| 国产精品一区二区三区四区久久| 中文在线观看免费www的网站| www日本黄色视频网| 国产亚洲欧美98| 欧美一区二区精品小视频在线| 亚洲欧美清纯卡通| 精品不卡国产一区二区三区| 六月丁香七月| 白带黄色成豆腐渣| av卡一久久| 亚洲成av人片在线播放无| 免费高清视频大片| 天堂网av新在线| 99久久中文字幕三级久久日本| 欧美高清性xxxxhd video| 村上凉子中文字幕在线| 桃色一区二区三区在线观看| 最好的美女福利视频网| 亚洲一区高清亚洲精品| 亚洲国产精品sss在线观看| 日本撒尿小便嘘嘘汇集6| 天天躁夜夜躁狠狠久久av| 亚洲四区av| 亚洲成人av在线免费| 久久韩国三级中文字幕| 少妇高潮的动态图| 99国产精品一区二区蜜桃av| 国产伦在线观看视频一区| 国产精品爽爽va在线观看网站| 成人欧美大片| 97在线视频观看| 一级毛片aaaaaa免费看小| 99热只有精品国产| 国产高清三级在线| 国产欧美日韩精品亚洲av| 国产又黄又爽又无遮挡在线| 亚洲av免费在线观看| 在线播放国产精品三级| 日韩欧美一区二区三区在线观看| 色尼玛亚洲综合影院| 一区福利在线观看| 久久久久久久久久久丰满| 看非洲黑人一级黄片| 免费在线观看成人毛片| 久久人人精品亚洲av| 午夜激情欧美在线| 免费看a级黄色片| 欧美色欧美亚洲另类二区| 婷婷精品国产亚洲av在线| 精品人妻熟女av久视频| 欧美又色又爽又黄视频| 最新在线观看一区二区三区| av天堂在线播放| 99热全是精品| 国产真实伦视频高清在线观看| 欧美丝袜亚洲另类| 国产三级在线视频| 91精品国产九色| 村上凉子中文字幕在线| 久久久国产成人精品二区| 村上凉子中文字幕在线| 精品日产1卡2卡| 亚洲精品日韩av片在线观看| 精品99又大又爽又粗少妇毛片| 3wmmmm亚洲av在线观看| 国内久久婷婷六月综合欲色啪| 色噜噜av男人的天堂激情| 午夜精品国产一区二区电影 | АⅤ资源中文在线天堂| 精品久久久久久久人妻蜜臀av| 国产精品日韩av在线免费观看| 日韩一区二区视频免费看| 女生性感内裤真人,穿戴方法视频| 国产精品野战在线观看| 乱人视频在线观看| 亚洲最大成人中文| 国产精品久久视频播放| 成年女人看的毛片在线观看| 国产人妻一区二区三区在| 日韩人妻高清精品专区| 成年免费大片在线观看| 黄片wwwwww| 成年av动漫网址| 男女做爰动态图高潮gif福利片| 国内久久婷婷六月综合欲色啪| 菩萨蛮人人尽说江南好唐韦庄 | 日本一二三区视频观看| 国产高潮美女av| av女优亚洲男人天堂| 国产乱人视频| 欧美色欧美亚洲另类二区| 2021天堂中文幕一二区在线观| 国产69精品久久久久777片| 深夜a级毛片| 欧美成人免费av一区二区三区| 精品人妻视频免费看| 老司机福利观看| 无遮挡黄片免费观看| 麻豆国产av国片精品| 婷婷色综合大香蕉| 99热全是精品| 神马国产精品三级电影在线观看| 免费看日本二区| 久久精品久久久久久噜噜老黄 | 中文字幕精品亚洲无线码一区| 久久久午夜欧美精品| 国产69精品久久久久777片| 亚洲成人中文字幕在线播放| 亚洲国产精品国产精品| eeuss影院久久| 国产麻豆成人av免费视频| 中文亚洲av片在线观看爽| 综合色丁香网| 欧美激情国产日韩精品一区| 欧美性猛交╳xxx乱大交人| 久久精品国产亚洲av涩爱 | 色尼玛亚洲综合影院| 色吧在线观看| 99精品在免费线老司机午夜| 欧美性猛交╳xxx乱大交人| 伊人久久精品亚洲午夜| 一边摸一边抽搐一进一小说| 精品一区二区三区人妻视频| aaaaa片日本免费| 亚洲乱码一区二区免费版| 老司机福利观看| 免费av毛片视频| 男女那种视频在线观看| 欧美性猛交黑人性爽| 成人永久免费在线观看视频| 欧美zozozo另类| 国产精品亚洲一级av第二区| 能在线免费观看的黄片| 亚洲欧美日韩无卡精品| 精品久久久久久久久久久久久| 国产精品,欧美在线| 国产精品一区二区免费欧美| 少妇人妻精品综合一区二区 | 日日啪夜夜撸| 国产精品伦人一区二区| 一进一出好大好爽视频| 大又大粗又爽又黄少妇毛片口| a级毛片a级免费在线| 日日摸夜夜添夜夜添av毛片| 男女视频在线观看网站免费| 精品一区二区三区人妻视频| 国产视频内射| 99久久精品热视频| 人妻丰满熟妇av一区二区三区| 日韩强制内射视频| 欧美不卡视频在线免费观看| 一本一本综合久久| 自拍偷自拍亚洲精品老妇| 校园春色视频在线观看| 精品午夜福利在线看| 午夜日韩欧美国产| 亚洲欧美日韩卡通动漫| 永久网站在线| 久久99热6这里只有精品| 国内揄拍国产精品人妻在线| 啦啦啦观看免费观看视频高清| 日韩强制内射视频| 精品福利观看| 久久久色成人| 国产三级在线视频| 精品久久久久久久久久久久久| 变态另类丝袜制服| 91狼人影院| 18禁在线无遮挡免费观看视频 | 一级av片app| 久久久久久久久大av| 成人欧美大片| 成人午夜高清在线视频| 成人特级黄色片久久久久久久| 国产淫片久久久久久久久| 日本与韩国留学比较| 欧美另类亚洲清纯唯美| 国产成年人精品一区二区| 日韩在线高清观看一区二区三区| av卡一久久| 12—13女人毛片做爰片一| 欧美3d第一页| 精品乱码久久久久久99久播| 日韩成人伦理影院| 日韩欧美精品免费久久| 神马国产精品三级电影在线观看| 亚洲精品乱码久久久v下载方式| 免费高清视频大片| 伊人久久精品亚洲午夜| 给我免费播放毛片高清在线观看| av在线观看视频网站免费| 久久精品国产亚洲av天美| 亚洲内射少妇av| 两性午夜刺激爽爽歪歪视频在线观看| 一边摸一边抽搐一进一小说| 国产精品久久久久久亚洲av鲁大| 中文字幕免费在线视频6| 亚洲熟妇中文字幕五十中出| 搞女人的毛片| 久久久久九九精品影院| 国产男靠女视频免费网站| 在线播放国产精品三级| 九九热线精品视视频播放| 日韩,欧美,国产一区二区三区 | 性插视频无遮挡在线免费观看| 寂寞人妻少妇视频99o| 麻豆乱淫一区二区| 草草在线视频免费看| 99久久精品一区二区三区| 女人十人毛片免费观看3o分钟| av.在线天堂| 久久人人爽人人片av| 久久婷婷人人爽人人干人人爱| 国产精品野战在线观看| 国产成人a∨麻豆精品| 久久久久久九九精品二区国产| 老司机福利观看| 久久精品影院6| 少妇熟女aⅴ在线视频| 男插女下体视频免费在线播放| 久久欧美精品欧美久久欧美| 亚洲第一电影网av| 色综合站精品国产| 在线观看66精品国产| 非洲黑人性xxxx精品又粗又长| 久久精品国产亚洲av香蕉五月| 国产激情偷乱视频一区二区| av卡一久久| 欧美+亚洲+日韩+国产| 亚洲国产精品久久男人天堂| 91久久精品国产一区二区三区| 黄色日韩在线| 国产麻豆成人av免费视频| 久久这里只有精品中国| 人人妻人人看人人澡| 国产精品久久视频播放| 婷婷精品国产亚洲av在线| 一本精品99久久精品77| 超碰av人人做人人爽久久| 久久精品91蜜桃| 免费看日本二区| 久久精品影院6| 黄色一级大片看看| 国产淫片久久久久久久久| 色噜噜av男人的天堂激情| 国产伦精品一区二区三区四那| 国产成人a∨麻豆精品| 日韩精品青青久久久久久| 久久久久精品国产欧美久久久| 最近视频中文字幕2019在线8| 欧美最新免费一区二区三区| 日本撒尿小便嘘嘘汇集6| 美女黄网站色视频| 一个人看视频在线观看www免费| 国产亚洲精品久久久com| 美女大奶头视频| 亚洲在线观看片| 女人十人毛片免费观看3o分钟| 国产爱豆传媒在线观看| 亚洲第一电影网av| 偷拍熟女少妇极品色| 不卡一级毛片| 精品99又大又爽又粗少妇毛片| 成人特级av手机在线观看| 女同久久另类99精品国产91| 成人鲁丝片一二三区免费| 日韩欧美三级三区| 久久鲁丝午夜福利片| 国产黄色视频一区二区在线观看 | 亚洲美女黄片视频| 黄色欧美视频在线观看| 久久韩国三级中文字幕| 日产精品乱码卡一卡2卡三| 一本精品99久久精品77| 国产黄色小视频在线观看| ponron亚洲| 尤物成人国产欧美一区二区三区| 超碰av人人做人人爽久久| 两个人的视频大全免费| 久久精品国产亚洲av天美| 18禁在线播放成人免费| 国产免费一级a男人的天堂| 欧美极品一区二区三区四区| a级毛片a级免费在线| 免费观看在线日韩| 精品99又大又爽又粗少妇毛片| 综合色av麻豆| 日本-黄色视频高清免费观看| 午夜精品国产一区二区电影 | videossex国产| 国产一区二区三区av在线 | 久久欧美精品欧美久久欧美| 国产高清三级在线| 亚洲国产精品成人综合色| 日产精品乱码卡一卡2卡三| 欧美一级a爱片免费观看看| 欧美激情在线99| 欧美又色又爽又黄视频| 亚洲在线自拍视频| 夜夜爽天天搞| 午夜福利18| 免费一级毛片在线播放高清视频| 性色avwww在线观看| eeuss影院久久| 国产精品不卡视频一区二区| 1024手机看黄色片| 菩萨蛮人人尽说江南好唐韦庄 | 一卡2卡三卡四卡精品乱码亚洲| 亚洲最大成人av| av在线蜜桃| 日本三级黄在线观看|