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

    鄂西南五種野生莢蒾屬植物種群結(jié)構(gòu)與動態(tài)特征研究

    2023-11-01 04:51:20趙奐敦李春輝艾訓(xùn)儒薛衛(wèi)星徐來仙向欽周云
    廣西植物 2023年9期

    趙奐敦 李春輝 艾訓(xùn)儒 薛衛(wèi)星 徐來仙 向欽 周云

    摘?要:為探明鄂西南地區(qū)野生莢蒾屬(Viburnum)植物的種群數(shù)量特征及區(qū)域分布狀況,并揭示宜昌莢蒾(V. erosum)、樺葉莢蒾(V. betulifolium)、合軸莢蒾(V. sympodiale)、茶莢蒾(V. setigerum)、莢蒾(V. dilatatum)5種莢蒾屬植物種群的生存現(xiàn)狀及發(fā)展趨勢,在鄂西南金子山國有林場、木林子國家級自然保護(hù)區(qū)及七姊妹山國家級自然保護(hù)區(qū),共布設(shè)27 hm2固定動態(tài)監(jiān)測樣地對莢蒾屬植物進(jìn)行調(diào)查,通過對5種莢蒾屬植物種群進(jìn)行年齡結(jié)構(gòu)、動態(tài)量化分析和靜態(tài)生命表以及相關(guān)曲線等研究,探討莢蒾屬植物種群結(jié)構(gòu)與動態(tài)特征及未來發(fā)展?jié)摿?。結(jié)果表明:(1)鄂西南莢蒾屬植物在分布區(qū)域及種群數(shù)量大小上均存在明顯差異。(2)種群結(jié)構(gòu)與動態(tài)分析顯示,5種莢蒾屬植物種群年齡結(jié)構(gòu)呈現(xiàn)為金字塔型,幼齡階段種群個(gè)體數(shù)量較多,具有較強(qiáng)的增長潛力,但對外界干擾均具有較強(qiáng)的敏感性。(3)靜態(tài)生命表顯示,5種莢蒾屬植物存活量均隨著齡級的增加而單調(diào)遞減;消失率與死亡率曲線變化趨勢相似,但不同植物消失率與死亡率曲線波動具有差異性;種群存活曲線均趨近于Deevey-Ⅱ型。(4)4個(gè)生存函數(shù)曲線表明,5種莢蒾屬植物種群均表現(xiàn)出前期大量減少、中期趨于穩(wěn)定、后期衰退的特點(diǎn)。綜上認(rèn)為,鄂西南莢蒾屬植物資源豐富,對其屬植物的合理開發(fā)具有一定的科學(xué)性和可行性,幼齡階段時(shí)可適當(dāng)進(jìn)行人工疏伐,提高保存率;中齡穩(wěn)定階段時(shí)可篩選觀賞性狀優(yōu)勢種進(jìn)行適當(dāng)開發(fā)利用;部分物種則需要加強(qiáng)保護(hù),以維持該地區(qū)生物多樣性。

    關(guān)鍵詞: 莢蒾屬, 種群結(jié)構(gòu), 動態(tài)量化分析, 靜態(tài)生命表, 生存函數(shù)

    中圖分類號:Q948

    文獻(xiàn)標(biāo)識碼:A

    文章編號:1000-3142(2023)09-1747-12

    收稿日期:2022-11-14

    基金項(xiàng)目:中國生態(tài)環(huán)境部生物多樣性調(diào)查評估項(xiàng)目(2019HJ2096001006); 湖北省林業(yè)科學(xué)研究院委托項(xiàng)目(H20064)。

    第一作者: 趙奐敦(1996-),碩士研究生,主要從事森林生態(tài)系統(tǒng)結(jié)構(gòu)與功能研究,(E-mail)2429989761@qq.com。

    *通信作者:李春輝,高級農(nóng)藝師,主要從事植物生態(tài)及生物資源開發(fā)利用研究,(E-mail)1106125155@qq.com。

    Population structure and dynamic characteristics of five

    wild Viburnum plants in Southwest Hubei, China

    ZHAO Huandun1, LI Chunhui2*, AI Xunru1, XUE Weixing1,

    XU Laixian1, XIANG Qin1, ZHOU Yun1

    ( 1. School of Forestry and Horticulture, Hubei Minzu University, Enshi 445000, Hubei, China; 2. Agriculture and Rural

    Affairs Bureau of Enshi Tujia and Miao Autonomous Prefecture, Enshi 445000, Hubei, China )

    Abstract:To explore the population quantitative characteristics and regional distribution of wild Viburnum plants in Southwest Hubei, and reveal the survival status and development trend of the following five Viburnum species: V. erosum, V. betulifolium, V. sympodiale, V. setigerum and V. dilatatum. A total of 27 hm2 dynamic monitoring stands were set up in Jinzi Mountains National Forest Farm, Mulinzi National Nature Reserve and Qizimei Mountains National Nature Reserve in Southwest Hubei. The age structure, dynamic quantitative analysis and static life table of five Viburnum species were used to explore the population structure, dynamic characteristics and future development potential of Viburnum plants. The results were as follows: (1) The distribution area and population size of Viburnum plants in Southwest Hubei were significantly different. (2) The age structure of the five Viburnum species displayed as the pyramidal type summarized from population structure and dynamic analysis. The number of individuals of each population was large at the young stage (89.48% of V. erosum in Jinzi Mountains, 76.38% of V. betulifolium, 67.24% of V. erosum in Mulinzi, 62.11% of V. dilatatum, 56.53% of V. sympodiale, 58.44% of V. setigerum), and only a small proportion of individuals (0.18% of V. erosum in Jinzi Mountains, 0.85% of V. betulifolium, 3.29% of V. erosum in Mulinzi, 3.97% of V. dilatatum, 8.37% of V. sympodiale, 2.79% of V. setigerum) at the old stage. The young stage plants of five Viburnum species all had strong growth potentials but were sensitive to external disturbances. (3) The survival quantity of the five Viburnum species was decreasing monotonically with increasing age class presented by static life tables. Except for V. setigerum, the life expectancy reached peak at the young age, all the other Viburnum species populations reached the peak at the middle age. The change trend of mortality rate and vanish rate were similar, but the fluctuation of vanish rate and mortality rate curves were different among different species. The survival curves were tended to be Deevey-Ⅱ type. (4) The four survival function curves showed that all five Viburnum species populations exhibited a substantial decrease at the early stage, stabilization at the middle stage, and decline at the late stage. To sum up, plant resources of Viburnum species are abundant in Southwest Hubei. Therefore, it is scientific and feasible to exploit Viburnum plants rationally. At the young stage, artificial thinning can be carried out appropriately to improve the preserving rate. During the stable middle stage, the dominant species which have good ornamental characters can be selected for appropriate development and utilization. However, some species need to be protected to maintain biodiversity in this area.

    Key words: Viburnum, population structure, dynamic quantitative analysis, static life table, survival function

    種群是物種和植物群落的基本單位,對種群的重要屬性即種群結(jié)構(gòu)與動態(tài)進(jìn)行研究也一直是種群生態(tài)學(xué)研究的主要內(nèi)容(Condit et al., 2000;Manabe et al., 2000;牛翠娟等,2015;劉浩棟等,2020)。種群結(jié)構(gòu)與動態(tài)既能反映種群在特定階段的數(shù)量分布和配置情況,又能揭示種群與環(huán)境之間的關(guān)系、種群動態(tài)變化及預(yù)測未來發(fā)展趨勢,對植物資源的保護(hù)和利用具有重要價(jià)值(金慧等,2017;何斌等,2020)。種群統(tǒng)計(jì)是研究種群結(jié)構(gòu)與動態(tài)的基本方法,主要包括種群年齡結(jié)構(gòu)、動態(tài)量化分析、靜態(tài)生命表、存活曲線、生存函數(shù)分析等核心內(nèi)容(王卓等,2009)。其中,種群年齡結(jié)構(gòu)及動態(tài)量化分析不僅可以揭示種群在各齡級階段的現(xiàn)存數(shù)量、種群的更新繁殖能力及種群未來發(fā)展的穩(wěn)定性,還能在一定程度上反映種群在所處環(huán)境中的適應(yīng)能力及種群在該群落中的地位和作用(劉普幸,2011;拓鋒等,2021)。種群動態(tài)量化分析是對種群結(jié)構(gòu)動態(tài)分析中等級歸類劃分粗放等缺點(diǎn)進(jìn)行彌補(bǔ)(陳曉德,1998)。靜態(tài)生命表可以直觀地展現(xiàn)種群的存活率、死亡率及生命期望等重要種群特征,是研究種群結(jié)構(gòu)與動態(tài)變化的主要手段,其與生存函數(shù)結(jié)合對種群進(jìn)行分析,對揭示種群的生命過程、生存現(xiàn)狀及種群在一定條件下生存、繁衍及變化的可能性具有重要意義(羅西等,2021;拓鋒等,2021)。因此,綜合利用各種群結(jié)構(gòu)與動態(tài)的研究方法,更有助于人們深入了解種群的現(xiàn)狀、動態(tài)特征及演替規(guī)律,預(yù)測種群未來發(fā)展趨勢,對種群的保護(hù)與利用具有重要指導(dǎo)意義(解婷婷等,2014)。

    莢蒾屬(Viburnum)植物隸屬于五?;疲ˋdoxaceae),全球約200種,是亞熱帶森林生態(tài)系統(tǒng)重要的組成部分,也是極具開發(fā)價(jià)值的觀賞灌木屬,在國際上被譽(yù)為“萬能綠化灌木”。湖北省分布有莢蒾屬植物43種,僅鄂西南恩施州就有28種(含22種6變種)(艾訓(xùn)儒等,2018)。因?yàn)樵搮^(qū)域最具代表性的植被類型為亞熱帶常綠落葉闊葉混交林,所以目前有關(guān)群落中植物種群結(jié)構(gòu)與動態(tài)的研究主要集中在優(yōu)勢種(姚蘭等,2017;朱強(qiáng)等,2019;王進(jìn)等,2020)、瀕危植物(楊心兵等,2000;林勇等,2017;羅西等,2021)等相關(guān)研究,而極少涉及在群落中較為劣勢卻具有較高觀賞價(jià)值的莢蒾屬植物。據(jù)報(bào)道,國內(nèi)約90%的莢蒾屬植物資源仍處于野生狀態(tài)(呂文君等,2019),并且對該屬的研究主要集中在常見種的種子休眠與萌發(fā)特性(肖月娥等,2007;楊軼華等,2015)、扦插繁殖(黃增艷和宋蕾,2007)、化學(xué)成分及元素含量測定(白長財(cái)?shù)龋?011;孫曉蕾等,2020)等方面,因此查清莢蒾屬植物資源的野生分布現(xiàn)狀,探究優(yōu)勢莢蒾屬植物的生存現(xiàn)狀與更新潛力,這對莢蒾屬植物原生種種質(zhì)資源的調(diào)查、保護(hù)及開發(fā)栽培利用等具有重要意義。

    本研究以鄂西南金子山國有林場、木林子國家級自然保護(hù)區(qū)及七姊妹山國家級自然保護(hù)區(qū)內(nèi)的莢蒾屬植物種群為研究對象,基于莢蒾屬植物資源野外調(diào)查數(shù)據(jù),篩選出5個(gè)種群個(gè)體數(shù)量較為優(yōu)勢的莢蒾屬植物為代表進(jìn)行種群結(jié)構(gòu)與動態(tài)特征分析。通過對5種優(yōu)勢莢蒾屬植物的種群個(gè)體數(shù)量特征、分布現(xiàn)狀、種群結(jié)構(gòu)與動態(tài)特征、野外生存現(xiàn)狀、自然更新與發(fā)展趨勢等進(jìn)行分析討論,以期為國內(nèi)莢蒾屬植物資源的開發(fā)利用等提供一定理論基礎(chǔ),同時(shí)為鄂西南地區(qū)生物多樣性的保護(hù)與維持提供一定科學(xué)依據(jù)。

    1?研究區(qū)概況

    鄂西南地處云貴高原與東部低山丘陵的過渡區(qū)域,與重慶、湖南相鄰,主要組成部分是恩施土家族苗族自治州,地勢起伏大,地形復(fù)雜多樣,主要為喀斯特地貌,境內(nèi)小氣候豐富,主要為亞熱帶大陸性季風(fēng)濕潤氣候,具有豐富的野生動植物資源。研究樣地位于鄂西南金子山國有林場、木林子國家級自然保護(hù)區(qū)和七姊妹山國家級自然保護(hù)區(qū)的常綠落葉闊葉混交林內(nèi),總面積為27 hm2,土壤類型主要為棕壤或黃棕壤,年均氣溫為12.7 ~ 15.5 ℃,年均降水量為1 475.0 ~ 1 667.5 mm,年日照時(shí)數(shù)為1 212~1 300 h,無霜期約為4個(gè)月(王進(jìn)等,2020);樣地中木植被有明顯垂直分層,群落物種組成多樣(姚蘭等,2017;朱強(qiáng)等,2019)(表1,圖1)。

    2?研究方法

    2.1 樣地建設(shè)及野外調(diào)查

    2014年4月至2016年9月,通過實(shí)地調(diào)查,在金子山國有林場、木林子國家級自然保護(hù)區(qū)、七姊妹山國家級自然保護(hù)區(qū)中選擇植被保存較為完整的天然林,利用全站儀按照CTFS(Center for Tropical Forest Science)樣地建設(shè)標(biāo)準(zhǔn)(Condit, 1998)布設(shè)固定監(jiān)測樣地,面積共27 hm2,包含675個(gè)20 m × 20 m的固定樣地,其中木林子大樣地為15 hm2,包含375個(gè)20 m × 20 m的固定樣地;金子山及七姊妹山大樣地均為6 hm2,分別包含150個(gè)20 m × 20 m的固定樣地(王進(jìn)等,2020)。野外調(diào)查以20 m × 20 m的固定樣地為基本單元,用紅色油漆對樣地中所有胸徑(DBH)≥ 1 cm的木本植物個(gè)體在距離地面1.3 m處的樹干上進(jìn)行標(biāo)記,并掛上特制的具有唯一編號的鋁牌;進(jìn)行物種學(xué)名(科、屬、種)鑒定,測定并記錄符合要求的木本植物個(gè)體的胸徑(DBH)、冠幅、樹高(H)及生長狀況等信息;根據(jù)設(shè)定的坐標(biāo)原點(diǎn),記錄樣地中所有木本植物個(gè)體的相對坐標(biāo)值(X,Y)。基于野外調(diào)查數(shù)據(jù),對所有樣地中的莢蒾屬植物進(jìn)行統(tǒng)計(jì),并分析各莢蒾屬植物種群個(gè)體數(shù)量及在各樣地中的分布情況。結(jié)合樣地及各莢蒾屬植物種群個(gè)體數(shù)量等實(shí)際情況,挑選5個(gè)種群個(gè)體數(shù)量較為優(yōu)勢的莢蒾屬植物進(jìn)行種群結(jié)構(gòu)與動態(tài)特征分析。

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

    采用徑級代替齡級的方法,結(jié)合樣地內(nèi)5種莢蒾屬植物的徑級特征,以5 mm為步長,按照莢蒾屬植物胸徑將大小級劃分為9個(gè)徑級:Ⅰ級,1 cm≤ DBH<1.5 cm;Ⅱ級,1.5 cm≤DBH<2 cm,Ⅲ級,2 cm≤DBH<2.5 cm;Ⅳ級,2.5 cm≤DBH<3 cm;Ⅴ級,3 cm≤DBH<3.5 cm;Ⅵ級,3.5 cm≤DBH<4 cm;Ⅶ級,4 cm≤DBH<4.5 cm;Ⅷ級,4.5 cm≤DBH<5 cm;Ⅸ級,DBH≥5 cm。根據(jù)徑級將所有個(gè)體劃分為幼齡(Ⅰ級,Ⅱ級)、低齡(Ⅲ級,Ⅳ級)、中齡(Ⅴ級,Ⅵ級)和高齡(Ⅶ級以上)4個(gè)階段。采用陳曉德(1998)的種群結(jié)構(gòu)動態(tài)量化分析方法,對5種莢蒾屬植物種群結(jié)構(gòu)進(jìn)行定量分析。其中,Vpi和V′pi分別為不考慮和考慮外部干擾因素下整個(gè)種群結(jié)構(gòu)的數(shù)量變化動態(tài)指數(shù),P極大(Pmax)為種群隨機(jī)干擾風(fēng)險(xiǎn)極大值。Vn、Vpi、V′pi取正、0、負(fù)值時(shí)分別反映種群個(gè)體數(shù)量的增長、穩(wěn)定、衰退的結(jié)構(gòu)動態(tài)關(guān)系。

    2.3 種群靜態(tài)生命表及存活曲線

    靜態(tài)生命表及存活曲線在研究種群的結(jié)構(gòu)與動態(tài)特征中具有重要作用,能直觀地展現(xiàn)種群在不同徑級的實(shí)際生存?zhèn)€體數(shù)、死亡數(shù)量及存活趨勢,由于所研究的莢蒾屬植物均為野生種群,為避免生命表中出現(xiàn)死亡率為負(fù)的情況,因此采用勻滑技術(shù)對數(shù)據(jù)進(jìn)行處理(江洪,1992;吳承禎等,2000)。將5種莢蒾屬植物種群全部個(gè)體按照大小級編制靜態(tài)生命表,并以大小級為橫坐標(biāo),繪制5種莢蒾屬植物的死亡率曲線、消失率曲線及存活曲線;采用Hett和Louckst(1976)提出的指數(shù)函數(shù)模型Nx=N0e-bx(代表Deevey-Ⅱ型曲線)與冪函數(shù)模型Nx=N0x-b (代表Deevey-Ⅲ型曲線) 對種群存活曲線進(jìn)行檢驗(yàn)。引入生存函數(shù)S(i)、i表示大小級、累積死亡率函數(shù)F(i)、死亡密度函數(shù)f(ti)和危險(xiǎn)率函數(shù)λ(ti)ti表示種群生存時(shí)間4個(gè)生存分析函數(shù)對5種莢蒾屬植物種群結(jié)構(gòu)進(jìn)行分析,hi為區(qū)間長度(李金昕等,2016)。

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

    采用Excel對文中的數(shù)據(jù)進(jìn)行處理及分析,使用Origin 2021軟件繪制死亡率與消失率曲線、存活曲線及生存函數(shù)曲線,并用SPSS軟件對各種群大小級進(jìn)行Kolomogorov-Smirnov檢驗(yàn)(簡稱KS檢驗(yàn))。

    3?結(jié)果與分析

    3.1 莢蒾屬植物種群數(shù)量特征及分布現(xiàn)狀

    據(jù)調(diào)查所得(表2),鄂西南3個(gè)區(qū)域莢蒾屬植物共13種,分別為巴東莢蒾(Viburnum henryi)、茶莢蒾(V. setigerum)、樺葉莢蒾(V. betulifolium)、莢蒾(V. dilatatum)、宜昌莢蒾(V. erosum)、水紅木(V. cylindricum)、合軸莢蒾(V. sympodiale)、球核莢蒾(V. propinquum)、短序莢蒾(V. brachybotryum)、粉團(tuán)(V. plicatum)、蝴蝶戲珠花(V. plicatum)、少花莢蒾(V. oliganthum)和陜西莢蒾(V. schensianum),各莢蒾屬植物的分布狀況及種群個(gè)體數(shù)量存在顯著差異。其中,合軸莢蒾分布在金子山、木林子樣地,球核莢蒾分布在金子山、七姊妹山樣地,少花莢蒾和陜西莢蒾僅分布在金子山樣地,短序莢蒾和粉團(tuán)僅分布在木林子樣地,蝴蝶戲珠花僅分布在七姊妹山樣地,其余莢蒾屬植物則在3個(gè)調(diào)查樣地均有分布。樣地中莢蒾屬植物共有8 198株,個(gè)體總數(shù)為金子山 > 木林子 > 七姊妹山,其中種群個(gè)體數(shù)量較為優(yōu)勢的有樺葉莢蒾、茶莢蒾、莢蒾、宜昌莢蒾、合軸莢蒾等,其余莢蒾屬植物種群個(gè)體數(shù)量均比較少,因此本研究以這5種優(yōu)勢種為主要研究對象。

    3.2 5種莢蒾屬植物種群結(jié)構(gòu)特征及動態(tài)變化分析

    對3個(gè)區(qū)域中各莢蒾屬植物種群大小級個(gè)體數(shù)量進(jìn)行分析,結(jié)果見圖2。由圖2可知,5種莢蒾屬植物種群個(gè)體數(shù)量均隨大小級的增加而總體呈現(xiàn)下降趨勢。5種莢蒾屬植物種群在幼齡階段的個(gè)體數(shù)量占比均較大,但相較于其他區(qū)域中莢蒾屬植物種群各大小級個(gè)體數(shù)量分布情況,金子山宜昌莢蒾、樺葉莢蒾種群在幼齡階段個(gè)體數(shù)量變化更明顯,說明在該樣地中宜昌莢蒾和樺葉莢蒾種群在幼齡階段對環(huán)境的適應(yīng)能力較弱。除了木林子樣地中合軸莢蒾種群個(gè)體數(shù)量在高齡階段占比大于中齡階段且種群大小級結(jié)構(gòu)呈現(xiàn)為不規(guī)則的金字塔型以外,其他莢蒾屬植物種群在高齡階段的個(gè)體數(shù)量占比均最小且種群大小級結(jié)構(gòu)均呈現(xiàn)為高度規(guī)則的金字塔型。由此可知,5種莢蒾屬植物種群均具有大量的幼齡樹,種群的增長狀態(tài)穩(wěn)定,表明鄂西南地區(qū)5種莢蒾屬植物種群具有長期穩(wěn)定的自我發(fā)展能力及維持持續(xù)繁衍更新的能力。對各莢蒾屬植物種群的大小級進(jìn)行精確性KS檢驗(yàn)得到,宜昌莢蒾(金子山)、樺葉莢蒾、宜昌莢蒾(木林子)、莢蒾、合軸莢蒾及茶莢蒾的KS值分別為0.160、0.429、0.466、0.537、0.405和0.748,不同種群呈現(xiàn)為正態(tài)分布,顯著性值均大于0.05。

    對3個(gè)區(qū)域中各莢蒾屬植物進(jìn)行動態(tài)量化分析,結(jié)果見表3。由表3可知,5種莢蒾屬植物種群均在Ⅸ級出現(xiàn)負(fù)值的情況,呈現(xiàn)衰退狀態(tài),但5種莢蒾屬植物種群數(shù)量動態(tài)變化指數(shù)Vpi(不考慮外部干擾)和V′pi(考慮外部干擾)均顯示Vpi>V′pi>0,說明5種莢蒾屬植物種群總體屬于增長型,這與種群大小級結(jié)構(gòu)分析結(jié)果一致。木林子的莢蒾、合軸莢蒾種群及七姊妹山的茶莢蒾種群V′pi指數(shù)均不大于1,說明在考慮外界干擾情況下,這3個(gè)莢蒾屬植物種群的動態(tài)變化趨勢處于增長型與穩(wěn)定型之間,各種群動態(tài)變化指數(shù)值由大到小為茶莢蒾、莢蒾、合軸莢蒾種群;金子山宜昌莢蒾、樺葉莢蒾種群及木林子宜昌莢蒾種群V′pi指數(shù)均大于1。種群結(jié)構(gòu)對隨機(jī)干擾的敏感性指數(shù)(Pmax)均大于0,說明5種莢蒾屬植物種群穩(wěn)定性受到外界環(huán)境干擾的敏感性強(qiáng),并且敏感性最高的是金子山宜昌莢蒾,最低的是木林子合軸莢蒾。

    3.3 5種莢蒾植物種群靜態(tài)生命表及存活曲線

    對3個(gè)區(qū)域中5種莢蒾屬植物種群的數(shù)據(jù)進(jìn)一步整理并計(jì)算得到種群靜態(tài)生命表,結(jié)果見表4。由表4可知,隨著大小級的增加,5種莢蒾屬植物種群的標(biāo)準(zhǔn)化存活個(gè)數(shù)(lx)逐漸減小,符合物種的生物學(xué)特性。種群個(gè)體的期望壽命(ex)除七姊妹山茶莢蒾種群外均在Ⅴ級處出現(xiàn)最高峰值, 其中合軸莢蒾種群的期望壽命值最大,且5種莢蒾屬植物種群的期望壽命值均在Ⅷ級時(shí)偏低,說明在Ⅷ級時(shí)各種群已處于生理衰退期。5種莢蒾屬植物種群消失率(Kx)與死亡率(qx)曲線變化趨勢一致(圖3),但不同莢蒾屬植物種群的死亡率與存活率曲線波動具有較大差異。金子山宜昌莢蒾、樺葉莢蒾種群的死亡率與存活率分別在Ⅱ級、Ⅰ級時(shí)段達(dá)到第一次峰值,其中,宜昌莢蒾種群死亡率與消失率在Ⅰ~Ⅴ級呈現(xiàn)先上升后下降的趨勢,而樺葉莢蒾種群死亡率與消失率在Ⅰ~Ⅴ級呈現(xiàn)先降后升的趨勢,但2種莢蒾屬植物種群的期望壽命值總體均呈現(xiàn)上升的趨勢。木林子3種莢蒾屬植物的死亡率與存活率則在Ⅰ級和Ⅱ級達(dá)到第一次峰值,3種莢蒾屬植物種群死亡率與消失率在Ⅰ~Ⅴ級呈現(xiàn)先升后降的趨勢,期望壽命值則先減后增;在Ⅴ級后,3種莢蒾屬植物種群死亡率與消失率變化較為緩慢,但期望壽命值下降明顯。七姊妹山樣地中茶莢蒾種群死亡率與消失率在Ⅰ~Ⅶ級的變化較為平緩,期望壽命值在Ⅱ級時(shí)達(dá)到峰值,之后整體呈下降趨勢。綜上認(rèn)為,5種莢蒾屬植物在幼齡、低齡階段種群易受到競爭、環(huán)境篩選等作用的影響,但在中齡階段后,除茶莢蒾以外,其他莢蒾屬植物受到的環(huán)境篩選和競爭壓力相對減弱,環(huán)境適應(yīng)能力增強(qiáng),并且合軸莢蒾、莢蒾在樣地中有更好的生存能力。

    不同樣地的莢蒾屬植物種群的存活曲線趨勢變化不一致,結(jié)果見圖4。由圖4可知,其中木林子3種莢蒾屬植物種群及金子山宜昌莢蒾種群的變化趨勢整體呈現(xiàn)為“凸-凹-凸”型,但七姊妹山茶莢蒾種群則在整體呈現(xiàn)出“凹-凸-凹”的變化趨勢,而金子山樺葉莢蒾種群則呈現(xiàn)為“凹-凸-凹-凸-凹”的變化趨勢。總體而言,5種莢蒾屬植物種群的存活曲線介于“Deevey-Ⅱ”和“Deevey-Ⅲ”型之間。結(jié)合Helt和Louckst提出的指數(shù)函數(shù)和對數(shù)函數(shù),對5種莢蒾屬植物種群存活曲線進(jìn)行檢驗(yàn)(表5),對比擬合的兩種模型的R2,各莢蒾屬植物種群的指數(shù)模型R2值均大于冪函數(shù)模型,且P值均小于0.001,結(jié)合靜態(tài)生命表和存活曲線分析認(rèn)為,可以明顯得出5種莢蒾屬植物種群存活曲線與“Deevey-Ⅱ”型相近,種群存活曲線的變化趨勢和模型檢驗(yàn)的結(jié)果保持一致,說明各莢蒾屬植物種群結(jié)構(gòu)總體處于穩(wěn)定狀態(tài),但不同種群存在一定差異性。

    3.4 5種莢蒾屬植物種群生存分析

    基于種群靜態(tài)生命表,對5種莢蒾屬植物種群的生存分析結(jié)果(圖5)顯示,各莢蒾屬植物種群累積死亡率隨大小級單調(diào)遞增,生存率隨大小級單調(diào)遞減,并且均表現(xiàn)為幼齡、低齡階段變化幅度大,在中齡及高齡階段變化幅度小,說明天然條件下5種莢蒾屬植物幼齡階段死亡率都較高,這與靜態(tài)生命表的分析結(jié)果一致。木林子3種莢蒾屬植物種群及七姊妹山樣地茶莢蒾種群的生存率與累積死亡率均在Ⅱ級時(shí)達(dá)到平衡,而金子山宜昌莢蒾、樺葉莢蒾種群的累積死亡率則一直高于生存率,這說明金子山宜昌莢蒾、樺葉莢蒾種群在幼齡階段之前已經(jīng)進(jìn)入衰退,木林子及七姊妹山莢蒾屬植物種群進(jìn)入衰退期的時(shí)間則在幼齡階段時(shí)期。5種莢蒾屬植物的死亡密度曲線均隨著大小級的增加而呈現(xiàn)波動下降的趨勢且在幼齡及低齡階段的值相對較高,說明5種莢蒾屬植物種群個(gè)體數(shù)量在幼齡及低齡階段大量銳減。5種莢蒾屬植物種群的危險(xiǎn)率在Ⅰ~Ⅴ級隨大小的增長而單調(diào)遞增,進(jìn)一步說明5種莢蒾屬植物種群個(gè)體數(shù)量在幼齡及低齡階段隨著齡級的增加而減少,之后危險(xiǎn)率變化趨勢則隨著大小級的增加而逐漸趨于平緩。4個(gè)生存函數(shù)曲線表明,5種莢蒾屬植物種群均表現(xiàn)為前期個(gè)體數(shù)大量減少,中期趨于穩(wěn)定,后期衰退的狀態(tài)。

    4?討論

    4.1 樣地中莢蒾屬植物種群數(shù)量特征及分布現(xiàn)狀

    種群數(shù)量特征可以表征種群個(gè)體數(shù)量及分配狀況,對于揭示群落的形成和維持機(jī)制具有重要意義(費(fèi)世民等,2008)。對鄂西南3個(gè)區(qū)域莢蒾屬植物個(gè)體數(shù)量對比分析得到,單位面積內(nèi)莢蒾屬植物個(gè)體數(shù)量排序依次為金子山 > 七姊妹山 > 木林子。3個(gè)區(qū)域中莢蒾屬植物均為灌木或小喬木,在群落中基本處于中下層,各區(qū)域地理位置及氣候相近,受水熱條件差異的影響較小,但群落演替階段存在差異性,群落演替程度為木林子 > 七姊妹山 > 金子山(艾訓(xùn)儒等,2018;王進(jìn)等,2020;李瑋宜等,2021),這說明各區(qū)域莢蒾屬植物個(gè)體數(shù)量受到群落演替階段的影響。從鄂西南各莢蒾屬植物種群個(gè)體數(shù)量來看,單位面積內(nèi)種群個(gè)體數(shù)量大小為樺葉莢蒾 > 茶莢蒾 > 莢蒾 > 宜昌莢蒾 > 合軸莢蒾等。由此可見,不同莢蒾屬植物在同一地區(qū)的種群大小存在明顯差異,這可能與不同種群對環(huán)境適應(yīng)能力、繁殖及生存能力的差異性有關(guān)。對5種莢蒾屬植物優(yōu)勢種群的分布情況進(jìn)行分析,樺葉莢蒾主要分布于金子山,茶莢蒾主要分布于七姊妹山,莢蒾與合軸莢蒾主要分布在木林子,宜昌莢蒾則在金子山與木林子均有分布。一般而言,同屬物種具有相似的進(jìn)化特征和生態(tài)適應(yīng)性,導(dǎo)致其對環(huán)境資源的利用也可能具有相似性,因此它們可能會面臨更強(qiáng)烈的種間競爭,進(jìn)而影響共存(Ackerly & Donoghue,1998)。綜上所述,鄂西南野生莢蒾屬植物資源豐富,但各莢蒾屬植物種群個(gè)體數(shù)量及分布情況均存在明顯的差異性。因此,對于種群較為優(yōu)勢的物種,可考慮進(jìn)行開發(fā)利用,并加強(qiáng)對其自身生物學(xué)性狀、競爭能力及環(huán)境影響因素等的進(jìn)一步研究;而較為劣勢或稀少的種群則需加強(qiáng)其野生種質(zhì)資源的保護(hù)。

    4.2 5種莢蒾屬植物種群結(jié)構(gòu)特征

    種群的徑級結(jié)構(gòu)可以替代種群年齡結(jié)構(gòu)來反映種群發(fā)展對生境的適應(yīng)與種群的生存狀態(tài)(陳育等,2014)。本研究對3個(gè)區(qū)域中5種莢蒾屬植物的種群結(jié)構(gòu)與動態(tài)特征分析顯示,各莢蒾屬植物種群的個(gè)體數(shù)量均隨著大小級的增大而逐漸減小,總體均呈現(xiàn)為金字塔結(jié)構(gòu),這與李功麟等(2021)對固沙灌木種群年齡結(jié)構(gòu)的研究結(jié)果相似,種群幼齡階段數(shù)量多,可以適應(yīng)當(dāng)?shù)厣常ㄟ^大量繁殖幼苗的策略來維持種群穩(wěn)定性,但高齡階段種群個(gè)體數(shù)量相對較少,反映各莢蒾屬植物種群在高齡階段的植株生存較為困難的現(xiàn)象。對各種群的大小級進(jìn)行KS檢驗(yàn),得出各種群均呈現(xiàn)正態(tài)分布。種群中個(gè)體對環(huán)境的適應(yīng)能力可以通過種群數(shù)量動態(tài)來反映(Bierzychudek, 2014)。本研究對各莢蒾屬植物種群動態(tài)量化分析顯示,5種莢蒾屬植物都對外部環(huán)境具有較高的敏感性,表現(xiàn)為種群整體均呈現(xiàn)為穩(wěn)定增長、后期衰退的趨勢。結(jié)合樣地中實(shí)際調(diào)查情況,莢蒾屬植物在森林群落中主要位于中下層的位置,在群落中屬于相對劣勢的物種,光照等環(huán)境因素也在一定程度上影響其種群增長的穩(wěn)定性,同層種間資源競爭壓力大,并且隨著幼齡個(gè)體的逐步增長,需要從環(huán)境中獲得更多的資源,因此種內(nèi)競爭也逐漸加劇。在不考慮和考慮外部干擾時(shí),各莢蒾屬植物種群增長潛力均為宜昌莢蒾(金子山) > 樺葉莢蒾 > 宜昌莢蒾(木林子)> 茶莢蒾 > 莢蒾 > 合軸莢蒾,但在考慮外部干擾時(shí),宜昌莢蒾及樺葉莢蒾種群具有更大的增長潛力;其他莢蒾屬植物的動態(tài)變化為處于增長型和穩(wěn)定型之間的過渡階段,這進(jìn)一步反映各莢蒾屬植物種群增長及動態(tài)變化趨勢的差異性可能是受到其生境及植物自身性狀的影響。

    4.3 5種莢蒾屬植物種群動態(tài)變化趨勢

    靜態(tài)生命表、存活曲線、生存函數(shù)分析等不僅能反映種群的生存狀況及環(huán)境對其的影響程度,還可以預(yù)測種群未來的發(fā)展趨勢(陳曉德,1998)。本研究對3個(gè)區(qū)域內(nèi)莢蒾屬植物種群數(shù)量進(jìn)行統(tǒng)計(jì),并編制種群靜態(tài)生命表以及相關(guān)曲線等來分析種群的生存狀態(tài)及發(fā)展趨勢。5種莢蒾屬植物種群存活曲線均趨近于Deevey-Ⅱ曲線,說明各莢蒾屬植物種群具有相對穩(wěn)定的發(fā)展趨勢,存活量隨大小級的增加而減少。七姊妹山茶莢蒾種群的期望壽命值在Ⅱ級時(shí)出現(xiàn)最高峰,其他莢蒾屬植物種群的期望壽命值均在Ⅴ級處出現(xiàn)最高峰且合軸莢蒾種群期望壽命值最大,說明七姊妹山茶莢蒾種群在幼齡階段時(shí)能夠更好地適應(yīng)生存環(huán)境;而其他莢蒾屬植物則在中齡階段時(shí)對環(huán)境的適應(yīng)能力較好,其中合軸莢蒾的適應(yīng)能力最強(qiáng),但隨著大小級的增加,各莢蒾屬植物種群在高齡階段均處于生理衰退期。這可能與七姊妹山樣地處于干擾后的中期階段且樣地中主要物種的大徑級個(gè)體在空間分布上表現(xiàn)出特定的生境偏好有關(guān)(姚良錦等,2017)。

    5種莢蒾屬植物種群的死亡率與存活率曲線研究結(jié)果顯示,各莢蒾屬植物種群的死亡率與存活率曲線變化趨勢較為接近,但不同莢蒾屬植物種群的死亡率與存活率曲線波動差異較為明顯,這可能與環(huán)境的篩選作用及各莢蒾屬植物的資源競爭力有關(guān)。各莢蒾屬植物種群的死亡率與存活率均在幼齡或低齡階段達(dá)到第一次峰值,由于各莢蒾屬植物種群在幼齡和低齡個(gè)體密度大,隨著生長過程中對資源和空間的需求增加,種內(nèi)及種間競爭能力增強(qiáng),因此由密度制約效應(yīng)引發(fā)的種群自疏和他疏作用導(dǎo)致種群有較高的死亡率(韓文軒和方精云,2008;何斌等,2020)。

    生存函數(shù)能夠更直觀、具體地說明種群的結(jié)構(gòu)和動態(tài)變化(李文英等,2018)。本研究生存函數(shù)分析結(jié)果與靜態(tài)生命表分析結(jié)果一致,5種莢蒾屬植物種群均表現(xiàn)為前期幼齡及低齡階段個(gè)體數(shù)大量減少,中期趨于穩(wěn)定,后期衰退的狀態(tài),這說明在天然更新下,雖然各莢蒾屬植物種群在幼齡階段死亡率較高,種群個(gè)體數(shù)量急劇下降,但大量的幼齡及低齡樹依舊對種群中齡階段具有較強(qiáng)的貢獻(xiàn),表明通過環(huán)境篩選及密度制約而存活的樹種對環(huán)境具有較強(qiáng)的適應(yīng)及資源競爭能力,有助于維持種群的平穩(wěn)發(fā)展,而在高齡階段,各莢蒾屬植物則進(jìn)入生理衰退時(shí)期(李文英等,2018)。

    目前,國內(nèi)關(guān)于野生莢蒾屬植物野外資源的研究較少,對莢蒾屬植物的認(rèn)知嚴(yán)重不足,雖然本研究掌握了鄂西南部分地區(qū)莢蒾屬植物的基礎(chǔ)數(shù)據(jù),并揭示了部分較為優(yōu)勢的莢蒾屬植物種群結(jié)構(gòu)與動態(tài)特征,但各莢蒾屬植物的功能性狀、繁殖更新能力、種內(nèi)及種間關(guān)系等還有待進(jìn)一步研究。因此,建議查明該地區(qū)其他莢蒾屬植物種群較為弱勢的原因及其種群結(jié)構(gòu)與動態(tài)變化,并加強(qiáng)對該地區(qū)各莢蒾屬植物的監(jiān)測與系統(tǒng)性研究,明確各莢蒾屬植物自身生物學(xué)特性及空間分布特征,為維護(hù)鄂西南地區(qū)的生物多樣性提供幫助,也為國內(nèi)其他地區(qū)的莢蒾屬植物資源的研究及開發(fā)利用提供科學(xué)依據(jù)。

    5?結(jié)論

    鄂西南地區(qū)是中國莢蒾屬植物分布的主要區(qū)域,莢蒾屬植物受到物種自身性狀、種間關(guān)系及群落演替階段的影響,在分布情況及種群個(gè)體數(shù)量方面都具有明顯差異。研究區(qū)內(nèi)個(gè)體數(shù)量較多的5種莢蒾屬植物為樺葉莢蒾、宜昌莢蒾、莢蒾、茶莢蒾和合軸莢蒾;5種莢蒾屬植物種群大小級結(jié)構(gòu)呈現(xiàn)為金字塔型,種群增長穩(wěn)定,幼齡個(gè)體數(shù)量多,自我更新能力較強(qiáng),具有較大的開發(fā)潛力,但種群又易受到外界環(huán)境干擾、物種的更新能力及資源競爭力等因素的影響,尤其是各種群在高齡階段時(shí)均受到較強(qiáng)的影響。因此,考慮該屬植物對環(huán)境資源需求的前提下,在對鄂西南莢蒾屬植物資源保護(hù)利用過程中,幼齡階段時(shí)可適當(dāng)進(jìn)行人工疏伐,提高保存率;中齡穩(wěn)定階段時(shí)可篩選性狀優(yōu)勢種進(jìn)行適當(dāng)開發(fā)利用。

    參考文獻(xiàn):

    AI XR, HUANG S, YAO L, et al., 2018. Flora of Enshi [M]. Beijing: Science Press, 4: 103-116.[艾訓(xùn)儒, 黃升, 姚蘭, 等, 2018. 《恩施植物志》 [M]. 北京: 科學(xué)出版社, 4: 103-116.]

    ACKERLY DD, DONOGHUE MJ, et al., 1998. Leaf size, sapling allometry, and Corners Rules: phylogeny and correlated evolution in maples [J]. Am Nat, 152(6): 767.

    BAI CC, HAN L, ZHAO YS, 2011. Determination of 12 inorganic elements in Viburnum betulifolium of fruit [J]. Asia-Pacific Trad Med, 7(10): 18-19.[白長財(cái), 韓璐, 趙云生, 2011. 樺葉莢蒾果實(shí)中12種無機(jī)元素含量測定 [J]. 亞太傳統(tǒng)醫(yī)藥, 7(10): 18-19.]

    BIERZYCHUDEK P, 2014. Plant biodiversity and population dynamics [M]//MONSON RK. Ecology and the Environment. New York: Springer, 8: 29-65.

    CONDIT R, ASHTON PS, BAKER P, et al., 2000. Spatial patterns in the distribution of tropical tree species [J]. Science, 288(5470): 1414-1418.

    CONDIT R, 1998.Tropical forest census plots: methods and results from Barro Colorado Island, Panama and a comparison with other plots[M]. Berlin: Springer: 23-55.

    CHEN XD, 1998. A study on the method of quantitative analysis for population and community structural dynamics [J]. Acta Ecol Sin, 18(2): 104-107.[陳曉德, 1998. 植物種群與群落結(jié)構(gòu)動態(tài)量化分析方法研究 [J]. 生態(tài)學(xué)報(bào), 18(2): 104-107.]

    CHEN Y, YANG J, ZHANG PJ, et al., 2014. Population structure and spatial point pattern of Helianthemum soongoricum in West Ordos, Inner Mongolia, China [J]. J Desert Res, 34(1): 75-82.[陳育, 楊劼, 張璞進(jìn), 等,2014. 西鄂爾多斯半日花(Helianthemum soongoricum)種群結(jié)構(gòu)和點(diǎn)格局分析 [J]. 中國沙漠, 34(1): 75-82.]

    FEI SM, HE YP, CHEN XM, et al., 2008. Quantitative features of populations of Pinus tabulaeformis and P. armandii regenerated following water damage at Qinling Mountain, China [J]. Chin J Plant Ecol, 32(1): 95-105.[費(fèi)世民, 何亞平, 陳秀明, 等, 2008. 秦嶺水災(zāi)跡地油松和華山松更新種群數(shù)量特征 [J]. 植物生態(tài)學(xué)報(bào), 32(1): 95-105.]

    HE B, LI Q, CHEN QL, et al., 2020. Quantitative characteristics and population dynamics of Pinus massoniana plantation in northwest Guizhou Province[J]. J Centr S Univ For Technol, 40(11): 129-137.[何斌, 李青, 陳群利, 等, 2020. 貴州省西北部馬尾松人工林種群數(shù)量特征與動態(tài) [J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 40(11): 129-137.]

    HETT JM, LOUCKS OL, 1976. Age structure models of balsam fir and eastern hemlock [J]. J Ecol, 3(64): 1029-1044.

    HAN WX, FANG JY, 2008. Self-thinning law in plant population : -3/2 vs -4/3? [J]. J Peiking Univ(Nat Sci Ed), 44(4): 661-668.[韓文軒, 方精云, 2008. 植物種群的自然稀疏規(guī)律—— -3/2還是-4/3? [J]. 北京大學(xué)學(xué)報(bào)(自然科學(xué)版), 44(4): 661-668.]

    JIN H, ZHAO Y, YIN H, et al., 2017. Quautitative characteristics and dynamic analysis of the endangered species Rhododendron chrysamium population in Changbai Mountain [J]. Chin J Ecol, 36(10): 3123-3130.[金慧, 趙瑩, 尹航, 等, 2017. 長白山瀕危植物牛皮杜鵑(Rhododendron chrysanthum)種群數(shù)量特征與動態(tài)分析 [J]. 生態(tài)學(xué)雜志, 36(11): 3123-3130.]

    JIANG H,1992. Picea asperata population ecology [M]. Beijing: China Foresty Press: 11-12.[江洪, 1992. 云杉種群生態(tài)學(xué) [M]. 北京: 中國林業(yè)出版社: 11-12.]

    LIU HD, CHEN Q, XU ZY, et al., 2020. Natural population structure and spatial distribution pattern of rare and endangered species of Dacrydium pectinatum [J]. Acta Ecol Sin, 40(9): 2985-2995.[劉浩棟, 陳巧, 徐志揚(yáng), 等, 2020. 珍稀瀕危樹種陸均松天然種群結(jié)構(gòu)與空間分布格局 [J]. 生態(tài)學(xué)報(bào), 40(9): 2985-2995.]

    LIU PX, 2011. Study on population structure and dynamics of Populus euphratica in the middle and lower reaches of the Shule River Basin Oasis, Hexi Corridor [J]. J Nat Resour, 26(3): 429-439.[劉普幸, 2011. 疏勒河中下游綠洲胡楊種群結(jié)構(gòu)與動態(tài)研究 [J]. 自然資源學(xué)報(bào), 26(3): 429-439.]

    LUO X, GUO QJ, YAO L, et al., 2021. Characteristics of natural population structure of endangered plant Liriodendron chinense [J]. J Centr S Univ For Technol, 41(7): 115-123.[羅西, 郭秋菊, 姚蘭, 等, 2021. 瀕危植物鵝掌楸的天然種群結(jié)構(gòu)特征 [J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 41(7): 115-123.]

    LIN Y, AI XR, YAO L, et al., 2017. Population structure and dynamics of Metasequoia glyptostroboides parent trees [J]. Chin J Ecol, 36(6): 1531-1538.[林勇, 艾訓(xùn)儒, 姚蘭, 等, 2017. 水杉原生母樹種群結(jié)構(gòu)與動態(tài) [J]. 生態(tài)學(xué)雜志, 36(6): 1531-1538.]

    L WJ, LIU HT, XIA BS, et al., 2019. Viburnum resources and their application to landscaping [J]. World For Res, 32(3): 36-41.[呂文君, 劉宏濤, 夏伯順, 等, 2019. 莢蒾屬植物資源及其園林應(yīng)用 [J]. 世界林業(yè)研究, 32(3): 36-41.]

    LI JX, WU DJ, ZHANG SP, et al., 2016. Life table and dynamic analysis of Fagus hayatae population in Micangshan Nature Reserve, Sichuan Province, China [J]. Bull Bot Res, 36(1): 68-74.[李金昕, 吳定軍, 章世鵬, 等, 2016. 四川米倉山自然保護(hù)區(qū)臺灣水青岡種群生命表及動態(tài)分析 [J]. 植物研究, 36(1): 68-74.]

    LI GL, ZHANG DH, ZHANG ZS, et al., 2021. Population dynamics of main sand-fixing shrubs in Gurbantunggut Desert [J]. J Desert Res, 41(2): 129-137.[李功麟, 張定海, 張志山, 等, 2021. 古爾班通古特沙漠沙丘主要灌木的種群數(shù)量動態(tài) [J]. 中國沙漠, 41(2): 129-137.]

    LI WY, LI X, GAN XH, 2018. Population structure and dynamics of endangered plant Tetracentron sinense [J]. Subtrop Plant Sci, 47(3): 222-228.[李文英, 李欣, 甘小洪, 2018. 瀕危植物水青樹的種群結(jié)構(gòu)與數(shù)量動態(tài) [J]. 亞熱帶植物科學(xué), 47(3): 222-228.]

    LI WY, AI XR, YAO L, et al., 2021. Biodiversity evaluation of three national nature reserves in Southwest Hubei [J]. Hubei Fore Sci Technol, 50(2):18-22. [李瑋宜, 艾訓(xùn)儒, 姚蘭, 等, 2021. 鄂西南3個(gè)國家級自然保護(hù)區(qū)生物多樣性評價(jià) [J].湖北林業(yè)科技, 50(2):18-22.]

    MANABE T, NISHIMURA N, MIURA M, et al., 2000. Population structure and spatial patterns for trees in a temperate old-growth evergreen broad-leaved forest in Japan [J]. Plant Ecol, 151(2): 181-197.

    NIU CJ, LOU AR, SUN RY,et al., 2015. Foundations in ecology [M]. Beijing: Higher Education Press: 64-65.[牛翠娟, 婁安如, 孫儒泳, 等, 2015. 基礎(chǔ)生態(tài)學(xué)(第3版) [M]. 北京: 高等教育出版社: 64-65.]

    SUN XL, ZHANG SY, CHEN WH, et al., 2020. Chemical constituents research progress of the genus Viburnum [J]. Chin J Mod Appl Pharm, 37(15): 1907-1913.[孫曉蕾, 張世堯, 陳文華, 等, 2020. 莢蒾屬植物化學(xué)成分的研究進(jìn)展 [J]. 中國現(xiàn)代應(yīng)用藥學(xué), 37(15): 1907-1913.]

    HUANG ZY, SONG L, 2007. The softwood-cutting propagation of Viburnum × carlcehalum [J]. Hubei Agris Sci, 46(6): 950-952.[黃增艷, 宋蕾, 2007. 紅蕾雪球莢蒾嫩枝扦插的繁殖技術(shù) [J]. 湖北農(nóng)業(yè)科學(xué), 46(6): 950-952.]

    TA F, LIU XD, HUANG DL, et al., 2021. Quantitative dynamics of Picea crassifolia population in Dayekou basin of Qilian Mountains [J]. Acta Ecol Sin, 41(17): 6871-6882.[拓鋒, 劉賢德, 黃冬柳, 等, 2021. 祁連山大野口流域青海云杉種群數(shù)量動態(tài) [J]. 生態(tài)學(xué)報(bào), 41(17): 6871-6882.]

    WANG Z, HUANG RF, WANG LH, et al., 2009. Analysis on life table of natural Sabina vulgaris population in MU Us Sandy Land [J]. J Desert Res, 29(1): 118-124.[王卓, 黃榮鳳, 王林和, 等, 2009. 毛烏素沙地天然臭柏種群生命表分析 [J]. 中國沙漠, 29(1): 118-124.]

    WANG J, YAO L, AI XR, et al., 2020. Structure and dynamic characteristics of Betula luminifera populations in different regions of Southwest Hubei Province, China [J]. Chin J Appl Ecol, 31(2): 357-365.[王進(jìn), 姚蘭, 艾訓(xùn)儒, 等, 2020. 鄂西南不同區(qū)域亮葉樺種群結(jié)構(gòu)與動態(tài)特征 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 31(2): 357-365.]

    WU CZ, HONG W, XIE JS, et al., 2000. Life table analysis of Tsuga longibracteata population [J]. Chin J Appl Ecol, 11(3): 333-336.[吳承禎, 洪偉, 謝金壽, 等, 2000. 珍稀瀕危植物長苞鐵杉種群生命表分析 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 11(3): 333-336.]

    XIE TT, SU PX, ZHOU ZJ, et al., 2014. Structure and dynamic characteristics of Calligonum mongolicum population in the desertoasis ecotone [J]. Acta Ecol Sin, 34(15): 4272-4279.[解婷婷, 蘇培璽, 周紫鵑, 等, 2014. 荒漠綠洲過渡帶沙拐棗種群結(jié)構(gòu)及動態(tài)特征 [J]. 生態(tài)學(xué)報(bào), 34(15): 4272-4279.]

    XIAO YE, ZHOU XY, ZHANG XQ, et al., 2007. Advance in research on dormancy and germination of Viburnum seed [J]. Seed, 26(6): 56-59.[肖月娥, 周翔宇, 張憲權(quán), 等, 2007. 莢蒾屬(Viburnum)種子休眠與萌發(fā)特性研究進(jìn)展 [J]. 種子, 26(6): 56-59.]

    YAO L, AI XR, YI YM, et al., 2017. Structure and dynamics of dominant populations in the mixed forest of subtropical evergreen and deciduous broad-leaved tree species in the Southwest of Hubei Province [J]. Sci Silv Sin, 53(2): 10-18.[姚蘭, 艾訓(xùn)儒, 易詠梅, 等, 2017. 鄂西南亞熱帶常綠落葉闊葉混交林優(yōu)勢種群的結(jié)構(gòu)與動態(tài)特征 [J]. 林業(yè)科學(xué), 53(2): 10-18.]

    YANG XB, LIU SX, YANG FS, 2000. A study on the structure of Davidia involucrata var. vilmoriniana population in Houhe, Hubei Province [J]. J Biol, 17(1): 16-18. [楊心兵, 劉勝祥, 楊福生, 2000. 湖北省后河自然保護(hù)區(qū)光葉珙桐種群結(jié)構(gòu)的研究 [J]. 生物學(xué)雜志, 17(1): 16-18.]

    YANG YH, GONG W, LI H, et al., 2015, Study on seed germination characteristics of 4 species of Viburnum [J]. Territ Nat Res Stud, (1): 88-90.[楊軼華, 宮偉, 李虹, 等, 2015. 4種莢蒾種子萌發(fā)特性研究 [J]. 國土與自然資源研究, (1): 88-90.]

    YAO LJ, YAO L, YI YM, et al., 2017. Species composition and communty structure of a 6-ha subtropical evergreen and deciduous broad-leaved mixed forest dynamics plot in the Qizimei Mountains, Hubei Province [J]. Biodivers Sci, 25(3): 275-284.[姚良錦, 姚蘭, 易詠梅, 等, 2017. 湖北七姊妹山亞熱帶常綠落葉闊葉混交林的物種組成和群落結(jié)構(gòu) [J]. 生物多樣性, 25(3): 275-284.]

    ZHU Q, AI XR, YAO L, et al., 2019. Structure and dynamics of Carpinus fargesiana population in southwest Hubei Province [J]. J Centr S Univ For Technol, 39(8): 93-100.[朱強(qiáng), 艾訓(xùn)儒, 姚蘭, 等, 2019. 鄂西南川陜鵝耳櫪種群結(jié)構(gòu)與動態(tài) [J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 39(8): 93-100.]

    (責(zé)任編輯?李?莉?王登惠)

    国产视频内射| 午夜日韩欧美国产| 国产亚洲精品一区二区www| 国产亚洲精品久久久久久毛片| 日韩成人在线观看一区二区三区| 国产精品九九99| 精品久久久久久久末码| 99热这里只有精品一区 | 精品第一国产精品| 亚洲中文字幕日韩| 深夜精品福利| 91国产中文字幕| www.精华液| 日本熟妇午夜| 亚洲国产看品久久| cao死你这个sao货| 欧美中文综合在线视频| 在线观看66精品国产| 男人舔女人下体高潮全视频| 性色av乱码一区二区三区2| 亚洲人成电影免费在线| 一区二区三区国产精品乱码| 国产精品综合久久久久久久免费| 一二三四社区在线视频社区8| 一级片免费观看大全| 国产v大片淫在线免费观看| 好看av亚洲va欧美ⅴa在| 又粗又爽又猛毛片免费看| 国产欧美日韩精品亚洲av| 欧美乱妇无乱码| 最近最新中文字幕大全电影3| 婷婷精品国产亚洲av| 18禁黄网站禁片免费观看直播| 最近视频中文字幕2019在线8| 麻豆成人午夜福利视频| 中文字幕av在线有码专区| 悠悠久久av| 五月玫瑰六月丁香| 精品熟女少妇八av免费久了| 久久久久久久久免费视频了| 亚洲 欧美 日韩 在线 免费| 丰满人妻一区二区三区视频av | 日韩欧美一区二区三区在线观看| 99久久无色码亚洲精品果冻| www.精华液| 精品福利观看| 宅男免费午夜| 在线观看日韩欧美| 久久精品国产综合久久久| 91字幕亚洲| 亚洲无线在线观看| 午夜老司机福利片| 99国产精品一区二区三区| 一级a爱片免费观看的视频| 国内毛片毛片毛片毛片毛片| 成年免费大片在线观看| 欧美激情久久久久久爽电影| 日韩精品免费视频一区二区三区| 色播亚洲综合网| 一卡2卡三卡四卡精品乱码亚洲| 国产免费男女视频| 日韩免费av在线播放| 精品久久久久久,| 村上凉子中文字幕在线| 巨乳人妻的诱惑在线观看| 精华霜和精华液先用哪个| 中文字幕人妻丝袜一区二区| 一二三四社区在线视频社区8| 久久香蕉国产精品| 亚洲精品粉嫩美女一区| 欧美黑人精品巨大| 中文字幕熟女人妻在线| 中出人妻视频一区二区| 久热爱精品视频在线9| 久久天堂一区二区三区四区| 男插女下体视频免费在线播放| 国产片内射在线| 在线观看免费视频日本深夜| 美女高潮喷水抽搐中文字幕| 一级作爱视频免费观看| 男男h啪啪无遮挡| 男人舔女人的私密视频| 亚洲av第一区精品v没综合| 欧美另类亚洲清纯唯美| 欧美成人午夜精品| 国产黄a三级三级三级人| 一级作爱视频免费观看| 熟妇人妻久久中文字幕3abv| 亚洲片人在线观看| 毛片女人毛片| 黑人操中国人逼视频| 女人高潮潮喷娇喘18禁视频| 免费观看人在逋| 成熟少妇高潮喷水视频| 亚洲人成网站在线播放欧美日韩| 亚洲精品美女久久久久99蜜臀| 在线观看日韩欧美| 欧美性长视频在线观看| 亚洲精品久久成人aⅴ小说| 免费看十八禁软件| 人妻久久中文字幕网| 波多野结衣巨乳人妻| 欧美高清成人免费视频www| 久久久久国产精品人妻aⅴ院| 亚洲熟女毛片儿| 亚洲精品一区av在线观看| 国产精品av视频在线免费观看| 女生性感内裤真人,穿戴方法视频| 一区福利在线观看| 国产野战对白在线观看| 亚洲精品中文字幕在线视频| 美女扒开内裤让男人捅视频| 午夜激情福利司机影院| 99在线视频只有这里精品首页| 99riav亚洲国产免费| 亚洲av第一区精品v没综合| videosex国产| 亚洲av美国av| 日韩欧美一区二区三区在线观看| 五月伊人婷婷丁香| 国产爱豆传媒在线观看 | 俄罗斯特黄特色一大片| 长腿黑丝高跟| 在线免费观看的www视频| 日本 av在线| 国产又色又爽无遮挡免费看| 欧美日本亚洲视频在线播放| 欧美+亚洲+日韩+国产| 国产一区在线观看成人免费| avwww免费| 精品久久久久久,| 2021天堂中文幕一二区在线观| 亚洲一区二区三区色噜噜| av福利片在线观看| 国内精品一区二区在线观看| 成人永久免费在线观看视频| 免费一级毛片在线播放高清视频| 精品国产乱子伦一区二区三区| 级片在线观看| 午夜成年电影在线免费观看| 桃色一区二区三区在线观看| 很黄的视频免费| bbb黄色大片| 午夜精品久久久久久毛片777| 禁无遮挡网站| 毛片女人毛片| 欧美成狂野欧美在线观看| 嫁个100分男人电影在线观看| 国产精品久久视频播放| a级毛片在线看网站| 一级片免费观看大全| 99久久无色码亚洲精品果冻| 一级片免费观看大全| a级毛片在线看网站| 亚洲国产欧美一区二区综合| 欧美3d第一页| 国产1区2区3区精品| 国产亚洲精品久久久久久毛片| 国产一区在线观看成人免费| 精品免费久久久久久久清纯| 国产97色在线日韩免费| 亚洲av成人一区二区三| 天天躁夜夜躁狠狠躁躁| 久久精品影院6| 一本一本综合久久| 国产精品一区二区免费欧美| 中文亚洲av片在线观看爽| 久久精品国产亚洲av香蕉五月| 免费人成视频x8x8入口观看| 精品欧美一区二区三区在线| 老司机深夜福利视频在线观看| 午夜福利18| av在线天堂中文字幕| 可以免费在线观看a视频的电影网站| 丁香六月欧美| 欧美一区二区精品小视频在线| 国产av不卡久久| 亚洲人成电影免费在线| 脱女人内裤的视频| 身体一侧抽搐| 18禁国产床啪视频网站| 在线观看免费日韩欧美大片| 50天的宝宝边吃奶边哭怎么回事| 黄色 视频免费看| 久久久久久大精品| 精品无人区乱码1区二区| 久久久久免费精品人妻一区二区| 美女免费视频网站| 国产成人精品无人区| 在线播放国产精品三级| 免费一级毛片在线播放高清视频| 在线观看免费视频日本深夜| 日本精品一区二区三区蜜桃| 最近最新免费中文字幕在线| 天堂影院成人在线观看| 久久香蕉精品热| 日韩欧美三级三区| 精品熟女少妇八av免费久了| netflix在线观看网站| 18禁观看日本| 亚洲av电影在线进入| 国产亚洲欧美98| 亚洲自偷自拍图片 自拍| 久久久久免费精品人妻一区二区| 久久人人精品亚洲av| 大型av网站在线播放| 伦理电影免费视频| 国产精品一区二区精品视频观看| 午夜免费成人在线视频| 天堂√8在线中文| 99在线视频只有这里精品首页| 两个人看的免费小视频| 日本黄大片高清| 国产av不卡久久| 可以在线观看的亚洲视频| 久久久久久久久中文| 精品电影一区二区在线| 亚洲中文av在线| 午夜精品在线福利| 久久午夜综合久久蜜桃| 免费无遮挡裸体视频| 亚洲成人中文字幕在线播放| 男女之事视频高清在线观看| 99热这里只有精品一区 | 人妻久久中文字幕网| 白带黄色成豆腐渣| 在线免费观看的www视频| 这个男人来自地球电影免费观看| 国产人伦9x9x在线观看| 亚洲人成网站在线播放欧美日韩| 一a级毛片在线观看| 18禁美女被吸乳视频| 欧美极品一区二区三区四区| 日韩精品青青久久久久久| 老鸭窝网址在线观看| 午夜日韩欧美国产| 深夜精品福利| 色av中文字幕| 亚洲欧美日韩无卡精品| 国内精品久久久久久久电影| 亚洲 欧美 日韩 在线 免费| 国产精品久久久久久人妻精品电影| 99精品在免费线老司机午夜| 婷婷六月久久综合丁香| 中出人妻视频一区二区| 国产亚洲精品第一综合不卡| 欧美日韩一级在线毛片| 亚洲中文字幕日韩| 久久天躁狠狠躁夜夜2o2o| 国产午夜精品论理片| 美女黄网站色视频| 麻豆国产av国片精品| 午夜久久久久精精品| 亚洲国产精品成人综合色| 1024视频免费在线观看| 国内久久婷婷六月综合欲色啪| 国产亚洲精品久久久久久毛片| 亚洲国产欧美人成| 国产成人一区二区三区免费视频网站| 无限看片的www在线观看| 丁香欧美五月| 亚洲国产精品成人综合色| 男女那种视频在线观看| 国产一区二区三区视频了| 亚洲人成77777在线视频| 18禁黄网站禁片午夜丰满| 观看免费一级毛片| 亚洲精品在线美女| 床上黄色一级片| 99精品久久久久人妻精品| 制服人妻中文乱码| 国产午夜福利久久久久久| 午夜福利免费观看在线| 免费观看精品视频网站| 88av欧美| 校园春色视频在线观看| 亚洲精品一区av在线观看| 麻豆成人午夜福利视频| 精品不卡国产一区二区三区| 人成视频在线观看免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 麻豆av在线久日| 好看av亚洲va欧美ⅴa在| 久久精品影院6| 精品久久久久久久久久久久久| 亚洲乱码一区二区免费版| 别揉我奶头~嗯~啊~动态视频| 国产精品98久久久久久宅男小说| 香蕉丝袜av| 一区二区三区国产精品乱码| 亚洲精品美女久久av网站| 亚洲国产精品久久男人天堂| 午夜福利欧美成人| 成人亚洲精品av一区二区| 2021天堂中文幕一二区在线观| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产高清在线一区二区三| 欧美黑人巨大hd| 国产伦一二天堂av在线观看| 淫秽高清视频在线观看| 久久久久久久久中文| 亚洲熟妇中文字幕五十中出| 日本一本二区三区精品| 午夜免费观看网址| 中文字幕最新亚洲高清| 香蕉av资源在线| 校园春色视频在线观看| avwww免费| www国产在线视频色| 伊人久久大香线蕉亚洲五| 好男人在线观看高清免费视频| 国产精品亚洲一级av第二区| 久久精品成人免费网站| 变态另类成人亚洲欧美熟女| a级毛片在线看网站| 久久久久国产精品人妻aⅴ院| 亚洲狠狠婷婷综合久久图片| 人人妻人人看人人澡| 亚洲精品美女久久久久99蜜臀| 亚洲天堂国产精品一区在线| 欧美午夜高清在线| 中文在线观看免费www的网站 | 脱女人内裤的视频| 变态另类丝袜制服| 九色国产91popny在线| 国产精品久久久人人做人人爽| 午夜福利成人在线免费观看| 色噜噜av男人的天堂激情| 久久久国产成人精品二区| 午夜福利在线在线| 桃色一区二区三区在线观看| 欧美日本亚洲视频在线播放| 欧美丝袜亚洲另类 | 99国产综合亚洲精品| 男男h啪啪无遮挡| 亚洲精品粉嫩美女一区| 国产1区2区3区精品| 亚洲精品美女久久av网站| 亚洲人成网站高清观看| 88av欧美| 亚洲 欧美一区二区三区| 天堂动漫精品| 可以在线观看毛片的网站| 岛国视频午夜一区免费看| 一本精品99久久精品77| 88av欧美| 啦啦啦韩国在线观看视频| 十八禁网站免费在线| 国产在线精品亚洲第一网站| 免费在线观看完整版高清| www.熟女人妻精品国产| 他把我摸到了高潮在线观看| 久久久久精品国产欧美久久久| 国产精品香港三级国产av潘金莲| 欧洲精品卡2卡3卡4卡5卡区| 婷婷六月久久综合丁香| 91大片在线观看| 亚洲成人国产一区在线观看| 国产aⅴ精品一区二区三区波| 91麻豆精品激情在线观看国产| 国产成人啪精品午夜网站| 亚洲片人在线观看| 级片在线观看| 五月玫瑰六月丁香| 黄色毛片三级朝国网站| 国产午夜精品久久久久久| 中文字幕久久专区| 亚洲av日韩精品久久久久久密| 久久 成人 亚洲| 麻豆国产97在线/欧美 | 日本精品一区二区三区蜜桃| 脱女人内裤的视频| 露出奶头的视频| 中文字幕av在线有码专区| 亚洲无线在线观看| 国产伦一二天堂av在线观看| 岛国在线免费视频观看| 国产乱人伦免费视频| 中亚洲国语对白在线视频| 国语自产精品视频在线第100页| 夜夜爽天天搞| 日本成人三级电影网站| 一本一本综合久久| 18禁观看日本| 最新在线观看一区二区三区| 国产一区在线观看成人免费| 高清在线国产一区| 99热只有精品国产| 免费观看精品视频网站| 91成年电影在线观看| 在线观看舔阴道视频| 黄色视频,在线免费观看| 久久精品亚洲精品国产色婷小说| cao死你这个sao货| 熟女电影av网| 亚洲欧美日韩高清专用| 又粗又爽又猛毛片免费看| 老司机午夜十八禁免费视频| 精华霜和精华液先用哪个| 欧美极品一区二区三区四区| 日日爽夜夜爽网站| 舔av片在线| 一本大道久久a久久精品| 国产亚洲精品久久久久5区| 少妇熟女aⅴ在线视频| 亚洲黑人精品在线| 中文亚洲av片在线观看爽| 成年人黄色毛片网站| 每晚都被弄得嗷嗷叫到高潮| 国产成人av教育| 亚洲av成人不卡在线观看播放网| 免费在线观看亚洲国产| 国产一区在线观看成人免费| 中文字幕最新亚洲高清| 国产黄a三级三级三级人| 亚洲中文日韩欧美视频| 搡老妇女老女人老熟妇| 精品福利观看| 欧美成狂野欧美在线观看| 久久久久久九九精品二区国产 | 久久精品夜夜夜夜夜久久蜜豆 | 午夜福利欧美成人| 天堂影院成人在线观看| 国产黄a三级三级三级人| 在线国产一区二区在线| 91av网站免费观看| 免费无遮挡裸体视频| 国产成人影院久久av| 久久婷婷成人综合色麻豆| 女警被强在线播放| 国产成人系列免费观看| 婷婷丁香在线五月| 国产欧美日韩一区二区三| 夜夜躁狠狠躁天天躁| 国产黄a三级三级三级人| 国产亚洲精品综合一区在线观看 | 一区二区三区激情视频| 国产精品美女特级片免费视频播放器 | 在线观看免费午夜福利视频| 神马国产精品三级电影在线观看 | 午夜免费观看网址| 亚洲精品av麻豆狂野| 国产精品精品国产色婷婷| 国内精品一区二区在线观看| 男人的好看免费观看在线视频 | 色av中文字幕| 丰满人妻一区二区三区视频av | 中文字幕精品亚洲无线码一区| 久久人人精品亚洲av| 亚洲avbb在线观看| 天堂动漫精品| 久久久久免费精品人妻一区二区| 一个人观看的视频www高清免费观看 | 婷婷亚洲欧美| 毛片女人毛片| 日本五十路高清| 99久久综合精品五月天人人| 美女黄网站色视频| 在线国产一区二区在线| 亚洲黑人精品在线| 免费在线观看影片大全网站| 国产精品亚洲av一区麻豆| av在线播放免费不卡| 一本精品99久久精品77| 九九热线精品视视频播放| 中文字幕高清在线视频| 久久久水蜜桃国产精品网| 欧美日韩一级在线毛片| 欧美成人免费av一区二区三区| 国产一级毛片七仙女欲春2| 国产成人精品无人区| 午夜日韩欧美国产| 九色成人免费人妻av| 日本一本二区三区精品| 久久中文看片网| 伦理电影免费视频| 99热只有精品国产| 成人一区二区视频在线观看| 一二三四社区在线视频社区8| 久久 成人 亚洲| 舔av片在线| 超碰成人久久| 欧美zozozo另类| 亚洲精品一卡2卡三卡4卡5卡| 18禁裸乳无遮挡免费网站照片| 成人特级黄色片久久久久久久| 给我免费播放毛片高清在线观看| 悠悠久久av| 国产麻豆成人av免费视频| 国产三级中文精品| 日本撒尿小便嘘嘘汇集6| 搡老岳熟女国产| √禁漫天堂资源中文www| 国产精品免费一区二区三区在线| 淫妇啪啪啪对白视频| 久久香蕉国产精品| 欧美日本视频| 搡老妇女老女人老熟妇| 成年版毛片免费区| 一边摸一边抽搐一进一小说| 色哟哟哟哟哟哟| 国产麻豆成人av免费视频| 人人妻,人人澡人人爽秒播| 夜夜夜夜夜久久久久| 97人妻精品一区二区三区麻豆| 欧美日韩亚洲国产一区二区在线观看| 色综合亚洲欧美另类图片| 黄色a级毛片大全视频| 欧美性猛交黑人性爽| 香蕉久久夜色| av视频在线观看入口| 国产亚洲精品久久久久久毛片| 后天国语完整版免费观看| 男人舔奶头视频| 国产v大片淫在线免费观看| 日韩欧美一区二区三区在线观看| 久久久精品大字幕| 国产亚洲欧美在线一区二区| 99国产精品99久久久久| 欧美性猛交黑人性爽| 欧美乱码精品一区二区三区| 国产精品av久久久久免费| 国产免费男女视频| cao死你这个sao货| 麻豆成人av在线观看| 成人手机av| 老熟妇仑乱视频hdxx| 两人在一起打扑克的视频| 国产野战对白在线观看| 婷婷丁香在线五月| 久久婷婷人人爽人人干人人爱| 久久国产精品人妻蜜桃| 99精品欧美一区二区三区四区| 亚洲精品av麻豆狂野| 国产一区二区三区在线臀色熟女| 国产午夜精品论理片| 亚洲一码二码三码区别大吗| 亚洲中文字幕日韩| 女人被狂操c到高潮| 88av欧美| 国内精品一区二区在线观看| 久久久久国内视频| 精品一区二区三区视频在线观看免费| 亚洲,欧美精品.| 欧美成狂野欧美在线观看| 小说图片视频综合网站| a在线观看视频网站| 在线a可以看的网站| 色精品久久人妻99蜜桃| 亚洲专区国产一区二区| 国产成人影院久久av| 中文亚洲av片在线观看爽| 亚洲精品久久国产高清桃花| 亚洲国产高清在线一区二区三| 欧美乱妇无乱码| 日本免费一区二区三区高清不卡| 精品无人区乱码1区二区| 真人做人爱边吃奶动态| 国产欧美日韩一区二区三| 国产人伦9x9x在线观看| 88av欧美| 真人做人爱边吃奶动态| 俄罗斯特黄特色一大片| 欧美绝顶高潮抽搐喷水| 国产精品久久久久久人妻精品电影| 一级黄色大片毛片| 亚洲色图 男人天堂 中文字幕| 悠悠久久av| 非洲黑人性xxxx精品又粗又长| 国产在线观看jvid| 国产激情久久老熟女| 亚洲熟妇熟女久久| 白带黄色成豆腐渣| 免费在线观看成人毛片| 亚洲精品色激情综合| 日韩精品中文字幕看吧| 国产精品美女特级片免费视频播放器 | 日韩免费av在线播放| www.熟女人妻精品国产| 国内毛片毛片毛片毛片毛片| 老汉色∧v一级毛片| 法律面前人人平等表现在哪些方面| 精品久久久久久久毛片微露脸| 又紧又爽又黄一区二区| 国产成人欧美在线观看| 女生性感内裤真人,穿戴方法视频| 在线观看66精品国产| 欧美3d第一页| 国产午夜精品久久久久久| 又紧又爽又黄一区二区| 免费看日本二区| 亚洲国产欧美人成| 日本成人三级电影网站| 狂野欧美激情性xxxx| 99在线人妻在线中文字幕| 亚洲专区国产一区二区| 久久精品国产99精品国产亚洲性色| 欧美乱码精品一区二区三区| 午夜亚洲福利在线播放| 一级黄色大片毛片| 亚洲人成电影免费在线| 欧美日本亚洲视频在线播放| 午夜福利成人在线免费观看| 亚洲,欧美精品.| 欧美高清成人免费视频www| 国产成人精品无人区| 精品国产超薄肉色丝袜足j| 非洲黑人性xxxx精品又粗又长| 国产真人三级小视频在线观看| 91大片在线观看| 一个人免费在线观看的高清视频| 久久久久国产精品人妻aⅴ院| 黄片大片在线免费观看| 国产爱豆传媒在线观看 | 成人三级黄色视频|