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

    河南省自然林草本層植物群落β多樣性及其影響因素

    2024-12-31 00:00:00孟偉杜曉軍焦志華高賢明劉龍昌王宇
    廣西植物 2024年10期
    關(guān)鍵詞:多樣性河南省

    摘 要:" 對反映不同群落間物種組成的時空差異或變化的β多樣性及其組分開展研究,有助于提高對群落結(jié)構(gòu)及其維持機(jī)制的理解和認(rèn)識。森林草本層是森林群落和森林生態(tài)系統(tǒng)的重要組成部分,盡管對林下草本層植物群落β多樣性已有一些相關(guān)研究,但關(guān)于β多樣性的組分(如物種周轉(zhuǎn)和嵌套組分)的相對比例、影響因素及其影響程度等方面仍缺乏統(tǒng)一認(rèn)識。該研究基于河南省168個自然林樣地草本層植物群落調(diào)查數(shù)據(jù)資料,通過Jaccard相異性指數(shù)、方差分解等方法分析河南省自然林草本層植物群落β多樣性格局,并量化空間距離、環(huán)境因素和物種豐富度等影響因素對其β多樣性的相對貢獻(xiàn)。結(jié)果表明:(1)自然林草本層植物群落總β多樣性的Jaccard相異性指數(shù)為0.94,主要來源于物種周轉(zhuǎn)組分(占總β多樣性的96.8%)。(2)自然林草本層植物群落總β多樣性及物種周轉(zhuǎn)組分隨空間距離增加呈顯著增加趨勢,而嵌套組分隨空間距離增加呈顯著減少趨勢。(3)Mantel檢驗(yàn)顯示空間距離、環(huán)境因素和物種豐富度對自然林草本層植物群落總β多樣性均具有顯著影響,其中環(huán)境因素影響最大,空間距離影響次之,物種豐富度影響最?。环讲罘纸饨Y(jié)果進(jìn)一步顯示空間距離、環(huán)境因素和物種豐富度共分別解釋了總β多樣性及其周轉(zhuǎn)組分27.14%、20.35%的變異,其中環(huán)境因素可以單獨(dú)解釋的變異最多(分別占總β多樣性的10.62%和周轉(zhuǎn)組分的9.35%)。該研究有助于增強(qiáng)對自然林草本層植物組成變化規(guī)律及其影響因素的認(rèn)識,進(jìn)而為森林植被保護(hù)管理和生態(tài)修復(fù)提供堅實(shí)的科學(xué)基礎(chǔ)。

    關(guān)鍵詞: β多樣性, 物種周轉(zhuǎn), 物種嵌套, 空間距離, 自然林草本層, 河南省

    中圖分類號:" Q948

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

    文章編號:" 1000-3142(2024)10-1931-12

    收稿日期:" 2024-05-03

    接受日期:" 2024-07-02

    基金項(xiàng)目:" 中國科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)(XDA0505020305); 國家自然科學(xué)基金(41471045,31370586)。

    第一作者: 孟偉(1996—),碩士研究生,主要從事植物多樣性研究,(E-mail)2393092463@qq.com。

    *通信作者:" 杜曉軍,博士,主要從事恢復(fù)生態(tài)學(xué)與生物多樣性研究,(E-mail)xjdu@ibcas.ac.cn;劉龍昌,博士,主要從事植物學(xué)研究,(E-mail)snowliu@126.com。

    β-diversity of the herbaceous layer plant communities in

    natural forests and its influencing factors in Henan Province

    MENG Wei1,2, DU Xiaojun2,3*, JIAO Zhihua4, GAO Xianming2,

    LIU Longchang1*, WANG Yu1

    ( 1. College of Horticulture and Plant Protection, Henan University of Science and Technology, Luoyang 471003, Henan, China; 2. State Key

    Laboratory of Vegetation and Environment Change, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China;

    3. China National Botanical Garden, Beijing 100093, China; 4. School of Chemical amp; Environmental

    Engineering, China University of Mining amp; Technology-Beijing, Beijing 100083, China )

    Abstract:" Conducting research on β-diversity (beta diversity) and its components, which reflect spatial or temporal variations in species composition among communities, helps to improve the understanding and knowledge of community structure and its maintenance mechanisms. The herbaceous layer of forest is an important part of forest community and forest ecosystem. Although some studies have been conducted on β-diversity of the understory herbaceous plant communities," there is still a lack of consensus on the proportion of β-diversity components such as turnover and nestedness components, as well as the influencing factors and the extent of their impact on β-diversity. This study was based on the survey data of herbaceous layer plant communities from 168 natural forest plots in Henan Province. It analyzed the β-diversity patterns of herbaceous layer in natural forests in Henan Province using Jaccard dissimilarity index, variance partitioning and other methods, and quantified the relative contributions of influencing factors such as spatial distance, environmental factor, and species richness to its β-diversity. The results were as follows: (1) The total β-diversity (as measured by the Jaccard dissimilarity index) of the herbaceous layer plant communities in natural forests was 0.94, which was mainly derived from the species turnover component (accounting for 96.8% of the total β-diversity). (2) The total β-diversity and its species turnover component of plant communities in the herbaceous layer of natural forests showed a significant increasing trend with increasing spatial distance, while the nestedness component showed a significant decreasing trend with increasing spatial distance. (3) Mantel test showed that spatial distance, environmental factor, and species richness had significant effects on the total β-diversity of plant communities in the herbaceous layer of natural forests, with environmental factor, spatial distance, and species richness having the greatest, second, and least effect, respectively; variance partitioning results further showed that spatial distance, environmental factor, and species richness accounting for 27.14% and 20.35% of the variation in total β-diversity and its turnover component, respectively, with environmental factor alone accounting for the most variation (10.62% of total β-diversity and 9.35% of turnover component). This study helps to enhance the understanding of the variation rules of plant composition in the herbaceous layers of natural forests and its influencing factors, thereby providing a solid scientific basis for the protection and management of forest vegetation and ecological restoration.

    Key words: β-diversity, species turnover, species nestedness, spatial distance, herbaceous layer of natural forest, Henan Province

    β多樣性反映不同群落間物種組成在時空上的差異或變化(陳圣賓等,2010; Jiang et al., 2021;曲夢君等,2022;Qian amp; Qian, 2023),它聯(lián)系著局域尺度的α多樣性和區(qū)域尺度的γ多樣性(Whittaker, 1972;李大標(biāo)等,2023)。β多樣性可分解成物種周轉(zhuǎn)和物種嵌套兩個組成部分,物種周轉(zhuǎn)描述了物種在空間或環(huán)境梯度中替換的現(xiàn)象,而物種嵌套描述的是在特定環(huán)境中物種的出現(xiàn)和消失(Baselga, 2010)。對β多樣性的周轉(zhuǎn)和嵌套組分及其影響因素研究,有助于對群落結(jié)構(gòu)及其運(yùn)行機(jī)制等的理解和認(rèn)識(斯幸峰等,2017;Soininen et al., 2018;姜小蕾等,2020;Jiang et al., 2022)。

    植物群落β多樣性分布格局受到多種因素的影響(陳圣賓等,2010;翁昌露等,2019;Shi et al., 2021),一般包括生物因素、非生物環(huán)境因素、空間距離等。其中,非生物環(huán)境因素如氣候(Wang et al., 2021;李星等,2022;曲夢君等,2022;楊欣等,2023)、海拔梯度(盧訓(xùn)令等,2010;Han et al., 2023;He et al., 2023)、坡位(譚珊珊等,2013;Hu et al., 2024)、坡度(譚珊珊等,2013;He et al., 2023)、凹凸度(楊潔等,2014;何存存等,2022)和土壤(王丹等,2013;Jiang et al., 2022;Xu et al., 2022)等對植物β多樣性有影響作用;空間距離與擴(kuò)散過程有關(guān),有研究報道植物群落間β多樣性隨著空間距離的增加而增大(Soininen et al., 2007;Chapman amp; McEwan, 2013;李新輝等,2016;Fang et al., 2024);冠層組成(Singh et al., 2023)、林冠大小(沙歡等,2024)等也能影響植物群落的β多樣性;除此之外,人類活動干擾(Schulz et al., 2019;高璐鑫等,2022)也會影響植物群落β多樣性格局。目前,對植物群落β多樣性形成和維持機(jī)制的研究大多關(guān)注生態(tài)位理論和中性理論:生態(tài)位理論認(rèn)為,物種組成同時受環(huán)境因素和物種相互作用的影響,物種組成隨環(huán)境梯度的變化而變化(Chase amp; Leibold, 2009;Fang et al., 2024),一般情況下當(dāng)群落間的環(huán)境差異增加時它們之間的物種組成差異也會增加,從而導(dǎo)致β多樣性的升高(Tuomisto et al., 2003;譚珊珊等,2013;Jiang et al., 2021);中性理論認(rèn)為擴(kuò)散限制是影響群落β多樣性格局的主要因素之一(Hubbell, 2001;Soininen et al., 2007;??瞬?,2009),群落β多樣性一般會因擴(kuò)散限制隨空間距離的增加而增大(Cacciatori et al., 2020;Jiang et al., 2021;楊歡等,2021;李大標(biāo)等,2023)。已有一些研究發(fā)現(xiàn)生境過濾和擴(kuò)散限制共同決定β多樣性等多樣性格局,二者的相對作用在不同群落及不同尺度存在差異(McClain amp; Rex, 2015;米湘成等,2022;Mugnai et al., 2022;Peguero et al., 2022;Zhao et al., 2023),在β多樣性相關(guān)研究中對生態(tài)位理論與中性理論等進(jìn)行整合考慮,從物種-功能-譜系等多維度綜合考慮也越來越受到更多重視。

    關(guān)于森林群落β多樣性的研究以往對高大喬木或木本植物關(guān)注的更多(翁昌露等,2019;Shi et al., 2021;楊欣等,2023),隨著對林下草本層在整個森林生態(tài)系統(tǒng)中的重要作用,如促進(jìn)生物多樣性,加速枯落物質(zhì)的分解以及推動養(yǎng)分循環(huán)等(Gilliam, 2007)的理解認(rèn)識,林下草本層也越發(fā)受到重視和關(guān)注(Nilsson amp; Wardle, 2005)。關(guān)于林下植物β多樣性的研究也陸續(xù)有報道,如塞罕壩主要森林類型林下植被(楊曉勤和田國恒,2012)、美國肯塔基州東南部的老齡林林下草本層(Chapman amp; McEwan, 2013)、北美東部的溫帶森林林下草本層(Buskey et al., 2020)、廣西憑祥紅椎林林下植物(尤業(yè)明等,2016)、黃土高原腹地人工林下草本層(施晶晶等,2019)、黃土高原油松和遼東櫟林下植物(王世雄等,2021)、雅魯藏布江下游高山峽谷區(qū)林下草本植物(杜方雪等,2022)、大興安嶺林區(qū)白樺天然林林下植物(張軼超,2023)等。林下草本層植物群落β多樣性通常受到多種因素影響,如Deng等(2023)認(rèn)為溫度是影響側(cè)柏以及濕地松林下植物群落多樣性重要因子,施晶晶等(2019)研究發(fā)現(xiàn),對樣地間物種組成變化解釋率最高的因子是地理距離,其次是年降水量;張軼超(2023)研究發(fā)現(xiàn),緯度差異影響著大興安嶺地區(qū)白樺天然林林下草本層物種組成。此外,還有研究發(fā)現(xiàn)海拔梯度影響著林下植物群落的組成(王世雄等,2021;杜方雪等,2022;Han et al., 2023),姜有為等(2021)發(fā)現(xiàn)林齡的增長會減弱白梭梭林下植被群落的異質(zhì)性,而毛亦楊等(2021)卻發(fā)現(xiàn)林齡對紅錐林下植被物種多樣性沒有顯著影響。盡管對林下草本層植物群落β多樣性已有一些相關(guān)研究,但對林下草本層植物群落β多樣性組分,如周轉(zhuǎn)組分和嵌套組分的占比、β多樣性的影響因素及其影響程度等問題仍缺乏統(tǒng)一認(rèn)識。

    河南省地處中原,同時作為一個農(nóng)業(yè)大省和人口大省,加強(qiáng)對該區(qū)自然植被的保護(hù)和修復(fù)的科學(xué)基礎(chǔ)研究(包括β多樣性等)尤為重要。本研究將以河南省168個自然林樣地的草本層植物群落為研究對象,對其總β多樣性及其組分的格局特征、影響因素等進(jìn)行分析研究,試圖回答如下科學(xué)問題:(1)自然林草本層的草本植物群落總β多樣性及其周轉(zhuǎn)組分、嵌套組分的格局有何規(guī)律;(2)自然林草本層的草本植物群落β多樣性隨空間距離的增大有何變化趨勢;(3)環(huán)境因素、空間距離、生物因素等對自然林草本層植物群落β多樣性及其組分的相對影響程度有何特點(diǎn)。研究結(jié)果有望為森林植被保護(hù)管理和生態(tài)修復(fù)提供科學(xué)基礎(chǔ)。

    1 材料與方法

    1.1 研究區(qū)概況

    河南省位于中國中東部其地理坐標(biāo)為110°21′—116°39′ E、31°23′—36°22′ N,河南省東西和南北直線距離分別約為580 km和550 km。大部分地區(qū)處于暖溫帶,南部跨亞熱帶,屬于北亞熱帶向暖溫帶過渡的大陸性季風(fēng)氣候。近十年河南省年平均氣溫為15.1~15.9 ℃,年均降水量為512.6~1 129.1 mm,年均日照時數(shù)為1 774.5~2 024.1 h,全年無霜期為207.9~271.7 d(https://www.henan.gov.cn/2023/06-13/2760285.html),具有由平原向丘陵山地氣候過渡的特征, 這種氣候四季分明,雨熱同期,復(fù)雜多樣。河南省的主要植被類型包括森林、灌叢、草地和農(nóng)田等,其中森林是本區(qū)自然植被類型的主體,其自然林資源主要分布在的大別山、伏牛山、太行山,占全省總自然林面積的90.09%(張莉等,2009)。林下草本層植物在森林群落結(jié)構(gòu)穩(wěn)定性的調(diào)節(jié)和生態(tài)功能的發(fā)揮上具有重要作用,是森林生態(tài)系統(tǒng)不可或缺的組成部分。

    1.2 樣地設(shè)置與物種豐富度獲取

    在河南省內(nèi)選取168個自然林環(huán)境下的草本層植物群落樣地進(jìn)行調(diào)查(圖1),記錄所有草本植物的物種名、高度、多度等基本信息,同時對每個樣地胸徑≥5 cm的喬木個體進(jìn)行每木調(diào)查,記錄其密度、平均高度和胸高斷面積和等結(jié)構(gòu)指標(biāo)。每個自然林草本層植物群落樣地面積為50 m × 20 m,隨后在喬木樣地中設(shè)立3個10 m × 10 m的灌木樣方,并在每個灌木樣方中設(shè)立1個2 m × 2 m草本層樣方,共布置了504個2 m × 2 m草本小樣方。用GPS儀測定168個樣地的經(jīng)緯度坐標(biāo),樣地之間南北距離最小和最大分別為0.12 km和485.40 km,樣地之間東西距離最小和最大分別為0.39 km和555.16 km,海拔高度52~1 720 m。

    1.3 環(huán)境因素獲取與分析

    本研究以野外調(diào)查時記錄的每個樣地的海拔高度數(shù)據(jù)為地形因子的代表。同時,從世界氣候數(shù)據(jù)網(wǎng)站(http://www.worldclim.org/)獲取所需要的19個氣候因子數(shù)據(jù)。根據(jù)樣地經(jīng)緯度坐標(biāo),用ArcGIS提取,包括年平均溫度(BIO1)、平均溫度日較差(BIO2)、等溫性(BIO3)、溫度季節(jié)性變化(BIO4)、最暖月最高溫度(BIO5)、最冷月最低溫度(BIO6)、氣溫年較差(BIO7)、最濕季度平均溫度(BIO8)、最干季度平均溫度(BIO9)、最暖季度平均溫度(BIO10)、最冷季度平均溫度(BIO11)、年降水量(BIO12)、最濕潤月降水量(BIO13)、最干月降水量(BIO14)、降水量季節(jié)性變化(BIO15)、最濕季度降雨量(BIO16)、最干季度降雨量(BIO17)、最暖季度降水量(BIO18)以及最冷季度降水量(BIO19)。

    為了避免多重共線性的影響,使用R軟件中內(nèi)置函數(shù)“scale”對海拔及上述19個氣候因子數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化,以消除量綱的差異,再用“Hmisc”包(Harrell, 2023)中“varclus”函數(shù)進(jìn)行聚類分析,計算環(huán)境變量之間的斯皮爾曼(Spearman)相關(guān)系數(shù)來評估共線性,并去除Spearman ρ2gt;0.7的變量,最終保留溫度季節(jié)性變化、最冷月最低溫度、最暖月最高溫度、最暖季度平均溫度、氣溫年較差、平均溫度日較差和最濕潤月降水量7個氣候因子指標(biāo)。林分結(jié)構(gòu),如喬木密度、喬木平均高度以及總胸高斷面積通過影響陽光的透過程度,進(jìn)而影響下層植物的生長。本研究將喬木密度、喬木平均高度、總胸高斷面積、海拔以及氣候因子(包括溫度季節(jié)性變化、最冷月最高溫、最暖月最高溫、最暖季度平均溫度、氣溫年較差、平均溫度月較差和最濕潤月降水量)共11個指標(biāo)作為環(huán)境因子納入分析,用“vegan”包(Oksanen et al., 2019)計算環(huán)境距離,即歐氏距離(Euclidean distance),用“geosphere”包(Hijmans et al., 2022)基于樣地位置經(jīng)緯度坐標(biāo)計算獲得空間距離。

    1.4 數(shù)據(jù)分析

    本研究通過計算不同樣地間β多樣性指數(shù)(Jaccard相異性指數(shù),Jaccard dissimilarity index)(Jaccard, 1912;陳圣賓等,2010;Fang et al., 2024)來表示物種組成的差異,該方法在計算時僅需考慮物種的有無(1代表有,0代表無)(斯幸峰等,2017;Fang et al., 2024)。使用R中“betapart”包(Baselga amp; Orme, 2012;Baselga et al., 2023)來計算總β多樣性(βtotal)及其周轉(zhuǎn)組分(βturnover)和嵌套組分(βnestedness),分析這兩個過程對群落多樣性差異有何影響。本研究主要用線性方程和對數(shù)方程對總β多樣性及其組分與空間距離之間的關(guān)系進(jìn)行回歸模擬分析。

    Jaccard相異性指數(shù)計算公式如下:

    βtotal=b+ca+b+c;

    βturnover=2min(b,c)a+2min(b,c);

    βnestedness=βtotal-βturnover。

    模擬方程如下:

    y=a1x+b1;

    y=a2lnx+b2。

    式中: a為兩個群落共有物種數(shù);b和c分別是兩個群落各自特有的物種數(shù);x為空間距離;y為β多樣性;a1、a2、b1、b2為系數(shù)。

    以空間距離、11個環(huán)境因子的歐氏距離和草本層植物物種豐富度的歐氏距離作為解釋變量,以河南省植物群落β多樣性及其組分作為響應(yīng)變量,利用Mantel分析來檢驗(yàn)空間距離、環(huán)境因子、物種豐富度及其兩兩配對以及三者之間共同作用與植物群落總β多樣性及其組分的相關(guān)性;使用偏Mantel檢驗(yàn),通過控制地理因素、環(huán)境因素或物種豐富度的影響,來檢驗(yàn)空間距離、環(huán)境因素、物種豐富度兩兩之間共同作用與植物群落總β多樣性及其組分的相關(guān)性。對每個檢驗(yàn)采用Pearson方法、設(shè)置999次重復(fù)的顯著性檢驗(yàn)來獲得相關(guān)性r值以及顯著性。

    先用“vegan”包(Oksanen et al., 2019)計算環(huán)境因素、物種豐富度的歐氏距離,再與地理距離基于相異性矩陣的多元回歸(multiple regression on dissimilarity matrices, MRM)分析來檢驗(yàn)其對植物群落組成的影響(Lichstein, 2007)。采用“vegan”包中“varpart”函數(shù)對空間距離、環(huán)境因素和物種豐富度進(jìn)行方差分解分析,將β多樣性分解為由空間距離、環(huán)境因素和物種豐富度3個解釋的部分以及未解釋的部分,分析各因素對植物群落β多樣性及其組分相對貢獻(xiàn)。以上數(shù)據(jù)分析及作圖均用R 4.2.3軟件。

    2 結(jié)果與分析

    2.1 自然林草本層植物群落β多樣性及其組分

    河南省自然林草本層植物群落β多樣性(此處用Jaccard相異性指數(shù))及其組分分析如圖2所示,研究區(qū)域自然林草本層植物群落間總β多樣性指數(shù)在0.94左右,利用“betapart”包將總β多樣性分解成周轉(zhuǎn)組分和嵌套組分后顯示植物群落調(diào)查樣地整體上是其周轉(zhuǎn)組分(0.91)遠(yuǎn)大于嵌套組分(0.03),周轉(zhuǎn)組分和嵌套組分分別占總β多樣性的96.8%和3.2%。

    2.2 自然林草本層植物群落β多樣性隨空間距離變化特征

    對河南省自然林草本層植物群落β多樣性及其組分隨空間距離變化分析結(jié)果如圖3所示:其中對數(shù)回歸模型擬合效果好于線性回歸擬合;自然林草本層植物群落總β多樣性及其周轉(zhuǎn)組分隨著空間距離的增加均呈增加趨勢,即呈現(xiàn)顯著的正相關(guān)關(guān)系(Plt;0.001);嵌套組分則隨著空間距離的增加呈降低趨勢,即呈現(xiàn)顯著的負(fù)相關(guān)關(guān)系(Plt;0.001)。

    2.3 空間距離、環(huán)境因素和物種豐富度對植物群落β多樣性及其組分的影響

    基于Mantel和偏Mantel檢驗(yàn)空間距離、環(huán)境因素和物種豐富度對植物群落β多樣性及其組分的影響,由表1結(jié)果可知:(1)對于研究區(qū)域的自然林草本層植物群落總β多樣性來說,按顯著相關(guān)性強(qiáng)弱排序?yàn)榭臻g距離-環(huán)境因素-物種豐富度三者共同作用、空間距離和環(huán)境因素共同作用、環(huán)境因素、環(huán)境因素和物種豐富度共同作用、空間距離、空間距離和物種豐富度共同作用、物種豐富度(Plt;0.001);(2)對于自然林草本層植物群落多樣性周轉(zhuǎn)組分來說,按顯著相關(guān)性強(qiáng)弱排序?yàn)榭臻g距離和環(huán)境因素共同作用、環(huán)境因素、空間距離-環(huán)境因素-物種豐富度三者共同作用、環(huán)境因素和物種豐富度共同作用、空間距離(Plt;0.001),物種豐富度以及空間距離和物種豐富度共同作用與其沒有顯著的相關(guān)性(Pgt;0.05);(3)對于自然林草本層植物群落多樣性嵌套組分來說,按顯著相關(guān)性強(qiáng)弱排序?yàn)榭臻g距離和物種豐富度共同作用、物種豐富度(Plt;0.001),而空間距離、環(huán)境因素、環(huán)境因素和物種豐富度共同作用、空間距離和環(huán)境因素共同作用、空間距離-環(huán)境因素-物種豐富度三者共同作用均與其沒有顯著的相關(guān)性(Pgt;0.05)。

    基于距離的多元回歸(MRM)來分析空間距離、環(huán)境因素和物種豐富度對β多樣性及其組分的影響,結(jié)果也進(jìn)一步表明空間距離、環(huán)境因素和物種豐富度對研究區(qū)植物群落的總β多樣性存在顯著影響(Plt;0.001,表2),其中空間距離對研究區(qū)植物群落的總β多樣性的單獨(dú)影響最大,而單獨(dú)的環(huán)境因素和物種豐富度的影響均低于空間距離的影響。除物種豐富度外的其他影響因素對于物種周轉(zhuǎn)組分也存在顯著影響。

    將自然林草本層植物群落β多樣性依據(jù)空間距離、環(huán)境因素和物種豐富度進(jìn)行方差分解其結(jié)果(圖4)如下:(1)對于自然林草本層植物群落β多樣性總體部分,空間距離、環(huán)境因素和物種豐富度共同解釋了27.14%,其余72.86%無法被現(xiàn)有的變量解釋;環(huán)境因素的單獨(dú)解釋率(10.62%)高于空間距離的單獨(dú)解釋率(7.88%)和物種豐富度的單獨(dú)解釋率(2.13%)。(2)對于自然林草本層植物群落β多樣性的周轉(zhuǎn)部分,空間距離、環(huán)境因素和物種豐富度共同解釋了20.35%,其余79.65%無法被現(xiàn)有的變量解釋;環(huán)境因素的單獨(dú)解釋率(9.35%)高于空間距離的單獨(dú)解釋率(6.01%)和物種豐富度的單獨(dú)解釋率(0.21%)。(3)對于自然林草本層植物群落β多樣性的嵌套組分,空間距離、環(huán)境因素和物種豐富度共同解釋了12.19%,其余87.81%無法被現(xiàn)有的變量解釋;物種豐富度的單獨(dú)解釋率(6.05%)高于環(huán)境因素的單獨(dú)解釋率(3.43%)和空間距離的單獨(dú)解釋率(1.57%)。

    3 討論與結(jié)論

    3.1 自然林草本層植物群落β多樣性及其組分

    關(guān)于將β多樣性進(jìn)一步分解成物種周轉(zhuǎn)組分和物種嵌套組分,先前一些研究表明周轉(zhuǎn)組分在β多樣性中起主導(dǎo)作用(李星等,2022;曲夢君等,2022;Fang et al., 2024)。本研究使用Jaccard相異性指數(shù)來衡量植物群落的β多樣性,研究結(jié)果也顯示,物種周轉(zhuǎn)組分對河南省自然林草本層植物群落總β多樣性的貢獻(xiàn)更大,占比96.8%(相較于嵌套組分占比3.2%),這表明該地區(qū)自然林草本層植物群落β多樣性主要由空間上或群落間的物種周轉(zhuǎn)組分主導(dǎo),物種組成差異主要來源于樣地間的物種替換,而物種嵌套的貢獻(xiàn)相對較低;河南省自然林草本層植物群落周轉(zhuǎn)組分主要受到空間距離、海拔、溫度季節(jié)性變化等的影響。由于研究發(fā)現(xiàn)植物群落間物種組成差異主要來源于周轉(zhuǎn)組分,因此對本區(qū)保護(hù)工作來說,如要盡可能保護(hù)更多自然林草本層物種,在多區(qū)域建保護(hù)區(qū)或多區(qū)域保護(hù)會更有幫助。

    3.2 自然林草本層植物群落β多樣性隨空間距離變化特征

    空間距離一般認(rèn)為是度量擴(kuò)散限制的相關(guān)指標(biāo)(Blundo et al., 2016)。本研究結(jié)果顯示,河南省自然林草本層植物群落總β多樣性及其物種周轉(zhuǎn)組分隨著空間距離增加而增大,這種變化趨勢與一些關(guān)于植物群落相似性與空間距離之間關(guān)系的研究結(jié)果相似(Qian et al., 2005;Cacciatori et al., 2020;楊歡等,2021;曲夢君等,2022;Fang et"" al., 2024),群落總β多樣性和周轉(zhuǎn)組分隨著空間距離增加而增大,可能是受到生態(tài)位過程和中性過程以及其他過程共同作用的影響。河南省自然林草本層植物群落總β多樣性的物種嵌套組分隨著空間距離的增加呈降低趨勢,這種格局可能和擴(kuò)散能力的差異性、環(huán)境異質(zhì)性、生物之間的相互作用等多種因素有關(guān)。在其他研究中也有類似發(fā)現(xiàn)或不同格局,例如:Fang等(2024)研究發(fā)現(xiàn),中國落葉松林草本層β多樣性的嵌套組分也隨著空間距離的增加呈降低趨勢;曲夢君等(2022)在對阿拉善戈壁植物群落β多樣性研究中發(fā)現(xiàn),嵌套組分隨空間距離的增加而降低(Srensen相異性指數(shù))或嵌套組分與空間距離沒有顯著相關(guān)關(guān)系(Bray-Curtis相異性指數(shù));楊歡等(2021)對庫姆塔格沙漠植物群落β多樣性(Bray-Curtis相似性指數(shù))的研究發(fā)現(xiàn),β多樣性嵌套組分無顯著的地理衰減格局。

    3.3 空間距離、環(huán)境因素和物種豐富度對β多樣性及其組分的影響

    越來越多的研究支持β多樣性受環(huán)境過濾和擴(kuò)散限制共同影響(王世雄等,2013;趙鳴飛等,2017;Peguero et al., 2022;Fang et al., 2024)。本研究結(jié)果表明,空間距離、環(huán)境因素和物種豐富度與河南省自然林草本層植物群落總β多樣性都有顯著相關(guān)性,意味著生境過濾、擴(kuò)散限制、群落結(jié)構(gòu)因素都對河南省自然林草本層植物群落β多樣性格局有影響作用;但三者對總β多樣性及其周轉(zhuǎn)組分的影響存在差異,對總β多樣性和周轉(zhuǎn)組分的單獨(dú)解釋率由大到小依次為環(huán)境因素、空間距離、物種豐富度;空間距離對群落總β多樣性和周轉(zhuǎn)組分的影響均低于環(huán)境因素說明在本研究中相比于空間距離因素來說環(huán)境因素影響更大,與曲夢君等(2022)和王健銘等(2022)的研究結(jié)果相似。本研究也分析了每個因子對植物群落β多樣性及其周轉(zhuǎn)組分的獨(dú)立貢獻(xiàn),發(fā)現(xiàn)環(huán)境因素中海拔因子對于總β多樣性及周轉(zhuǎn)組分影響最高,平均溫度月較差、最暖季度平均溫度和最暖月最高溫度對群落總β多樣性及周轉(zhuǎn)組分的影響也較大,其他研究也分析了海拔以及氣候因子對解釋多樣性格局的重要性(曲夢君等,2022;李星等,2022)。海拔可通過對氣候、土壤等因素的作用對植物群落產(chǎn)生影響,研究區(qū)域樣地之間最大海拔高差在1 600 m以上,這可能是海拔因子對草本層植物群落總β多樣性和周轉(zhuǎn)組分影響較高的原因。本研究結(jié)果顯示,環(huán)境因素和空間距離是河南省自然林草本層植物群落總β多樣性和周轉(zhuǎn)組分的重要影響因素且環(huán)境因素的影響更大,這提醒我們在進(jìn)行本區(qū)自然林草本層植物多樣性保護(hù)和恢復(fù)時首先應(yīng)重視環(huán)境因素,同時也應(yīng)重視空間因素及其與環(huán)境因素的交互作用。

    另外,方差分解結(jié)果顯示,本研究所考慮的三種因素之外仍有較多未解釋的比率,這一定程度也說明群落物種組成是一個多因素多生態(tài)過程綜合影響的格局,本研究所考慮的三種因素并不能完全解釋其變化,這也是后續(xù)研究需要著重注意的一個地方。有研究表明土壤理化性質(zhì)、土壤養(yǎng)分等因子與植物β多樣性相關(guān)(王丹等,2013;Jiang et al., 2022;Xu et al., 2022;李大標(biāo)等,2023),人為干擾也會對植物β多樣性產(chǎn)生一定的影響(Schulz et al., 2019;高璐鑫等,2022),后續(xù)如能補(bǔ)充土壤理化性質(zhì)、土壤養(yǎng)分等數(shù)據(jù)可能有助于提高本研究的β多樣性解釋度,以致提高對β多樣性格局及其生態(tài)機(jī)制的理解認(rèn)識程度。

    致謝 感謝數(shù)十位大學(xué)生、研究生、科研人員和師傅們參加河南森林樣地野外調(diào)查。

    參考文獻(xiàn):

    BASELGA A, 2010. Partitioning the turnover and nestedness components of beta diversity [J]. Global Ecol Biogeogr, 19(1): 134-143.

    BASELGA A, ORME CDL, 2012. Betapart: an R package for the study of beta diversity [J]. Meth Ecol Evol, 3(5): 808-812.

    BASELGA A, ORME D, VILLEGER S, et al., 2023. Betapart: Partitioning beta diversity into turnover and nestedness components [CP]. R package version 1.6. https://CRAN.R-project.org/package=betapart.

    BLUNDO C, GONZALEZ-ESPINOSA M, MALIZIA LR, 2016. Relative contribution of niche and neutral processes on tree species turnover across scales in seasonal forests of NW Argentina [J]. Plant Ecol, 217: 359-368.

    BUSKEY TM, MALONEY ME, CHAPMAN JI, et al., 2020. Herb-layer dynamics in an old-growth forest: vegetation-environment relationships and response to invasion-related perturbations [J]. Forests, 11(10): 1069.

    CACCIATORI C, TORDOIN E, PETRUZZELLIS F, et al., 2020. Drivers of distance-decay in bryophyte assemblages at multiple spatial scales: dispersal limitations or environmental control?" [J]. J Veg Sci, 31(2): 293-306.

    CHAPMAN JI, MCEWAN RW, 2013. Spatiotemporal dynamics of α- and β-diversity across topographic gradients in the herbaceous layer of an old-growth deciduous forest [J]. Oikos, 122(12): 1679-1686.

    CHASE JM, LEIBOLD MA, 2009. Ecological niches: linking classical and contemporary approaches [M]. Chicago: University of Chicago Press.

    CHEN SB, OUYANG ZY, XU WH, et al., 2010. A review of beta diversity studies [J]. Biodivers Sci, 18(4): 323-335. [陳圣賓, 歐陽志云, 徐衛(wèi)華, 等, 2010. Beta多樣性研究進(jìn)展 [J]. 生物多樣性, 18(4): 323-335.]

    DENG N, CAIXIA L, MA F, et al., 2023. Understory vegetation diversity patterns of Platycladus orientalis and Pinus elliottii communities in Central and Southern China [J]. Open Life Sci, 18: 20220791.

    DU FX, LI JJ, WU JH, et al., 2022. Study on the variation characteristics of understory herb diversity along the altitude gradient in the alpine valley area of the lower reaches of the Yarlung Zangbo River [J]. Environ Ecol, 4(12): 37-45. [杜方雪, 李景吉, 吳健輝等, 2022. 雅魯藏布江下游高山峽谷區(qū)林下草本植物多樣性沿海拔梯度變化特征研究 [J]. 環(huán)境生態(tài)學(xué), 4(12): 37-45.]

    FANG WJ, CAI Q, JI CJ, et al., 2024. Life forms affect beta-diversity patterns of larch forests in China [J]. Plant Divers, 46(1): 49-58.

    GAO LX, LAN TY, ZHAO ZX, et al., 2022. Spatial turnover of shrubland communities and underlying factors in northern mid-subtropical China [J]. Chin J Plant Ecol, 46(11): 1411-1421. [高璐鑫, 蘭天元, 趙志霞, 等, 2022. 中國中亞熱帶北部灌叢群落植物空間周轉(zhuǎn)及其驅(qū)動因素 [J]. 植物生態(tài)學(xué)報, 46(11): 1411-1421.]

    GILLIAM FS, 2007. The ecological significance of the herbaceous layer in temperate forest ecosystems [J]. Bioscience, 57(10): 845-858.

    HAN J, YIN H, XUE J, et al., 2023. Vertical distribution differences of the understory herbs and their driving factors on shady and sunny slopes in high altitude mountainous areas [J]. Front For Glob Change, 6: 1138317.

    HARRELL JRFE, 2023. Hmisc: Harrell Miscellaneous [CP]. R package version 5.1-1. https://CRAN.R-project.org/package=Hmisc.

    HE CC, ZHANG CY, ZHAO XH, 2022. Factors influencing β diversity of coniferous and broad-leaved mixed forest in Jiaohe, Jilin Province, China" [J]. Chin J Appl Environ Biol, 28 (2): 325-330." [何存存, 張春雨, 趙秀海, 2022. 吉林蛟河針闊混交林群落β多樣性的影響因素 [J]. 應(yīng)用與環(huán)境生物學(xué)報, 28(2): 325-330.]

    HE X, ARIF M, ZHENG J, et al., 2023. Plant diversity patterns along an elevation gradient: the relative impact of environmental and spatial variation on plant diversity and assembly in arid and semi-arid regions [J]. Front Environ Sci, 11: 1021157.

    HIJMANS RJ, WILLIAMS E, VENNES C, et al., 2022. Geosphere: Spherical Trigonometry [CP]. R package version 1.5-18. https://CRAN.R-project.org/package=geosphere.

    HU G, PANG Q, HU C, et al., 2024. Beta diversity patterns and determinants among vertical layers of tropical seasonal rainforest in karst peak-cluster depressions [J]. Forests, 15: 365.

    HUBBELL SP, 2001. The unified neutral theory of biodiversity and biogeography [M]. New Jersey: Princeton University Press.

    JACCARD P, 1912. The distribution of the flora in the alpine zone [J]. New Phytol, 11(2): 37-50.

    JIANG L, LV G, GONG Y, et al., 2021. Characteristics and driving mechanisms of species beta diversity in desert plant communities [J]. PLoS ONE, 16(1): e0245249.

    JIANG LM, HU D, WANG HF, et al., 2022. Discriminating ecological processes affecting different dimensions of α- and β-diversity in desert plant communities [J]. Ecol Evol, 12(3): e8710.

    JIANG XL, SUN ZY, HAO Q, et al., 2020. Interpretation of environmental factors affecting beta diversity and its components of secondary forest in Lao Mountain [J]. Chin J Ecol, 39(10): 3211-3220" [姜小蕾, 孫振元, 郝青, 等, 2020. 嶗山次生林群落β多樣性格局及其組分的驅(qū)動因素 [J]. 生態(tài)學(xué)雜志, 39(10): 3211-3220.]

    JIANG YW, ZHANG H, TAO HF, et al., 2021. Study on succession pattern of undergrowth plants during the recovery process of vegetation in the Haloxylon persicum plantation [J]. J Xinjiang Agric Univ, 44(3): 213-222." [姜有為, 張恒, 陶洪飛, 等, 2021. 白梭梭人工林林下植被恢復(fù)的演替規(guī)律研究 [J]. 新疆農(nóng)業(yè)大學(xué)學(xué)報, 44(3): 213-222.]

    LI DB, ZHONG YP, GONG XF, et al., 2023. Analysis of species and phylogenetic β diversity drivers in the Masson pine forests in Suichang, Zhejiang Province [J]. Guihaia, 43(7): 1258-1267. [李大標(biāo), 鐘毓萍, 龔笑飛, 等, 2023. 浙江遂昌馬尾松林物種和譜系β多樣性驅(qū)動因子分析 [J]. 廣西植物, 43(7): 1258-1267.]

    LI X, XIN ZM, DONG X, et al., 2022. β diversity and explanation of plant communities in Beishan and Mane Mountain areas of Dunhuang [J]. Arid Zone Res, 39(5): 1464-1472. [李星, 辛智鳴, 董雪, 等, 2022. 敦煌北山及馬鬃山地區(qū)植物群落β多樣性及其解釋 [J]. 干旱區(qū)研究, 39(5): 1464-1472.]

    LI XH, LIU YH, LIU Y, et al., 2016. Impact of geographical distance and environmental differences on the beta diversity of plant communities in the dry-hot valley of the Yuanjiang River [J]. Biodivers Sci, 24(4): 399-406. [李新輝, 劉延虹, 劉曄, 等, 2016. 地理距離及環(huán)境差異對云南元江干熱河谷植物群落beta多樣性的影響 [J]. 生物多樣性, 24(4): 399-406.]

    LICHSTEIN JW, 2007. Multiple regression on distance matrices: a multivariate spatial analysis tool [J]. Plant Ecol, 188: 117-131.

    LU XL, HU N, DING SY, et al., 2010. The pattern of plant species diversity of Funiu Mountain nature reserve [J]. Acta Ecol Sin, 30(21): 5790-5798. [盧訓(xùn)令, 胡楠, 丁圣彥, 等, 2010. 伏牛山自然保護(hù)區(qū)物種多樣性分布格局 [J]. 生態(tài)學(xué)報, 30(21): 5790-5798.]

    MAO YY, CHEN FQ, GUO Y, et al., 2021. Restoration dynamics of understory species diversity in Castanopsis hystrix plantations at different ages [J]. Chin J Appl Environ Biol, 27(4): 930-937." [毛亦楊, 陳富強(qiáng), 郭勇, 等, 2021. 不同林齡紅錐人工林林下植被物種多樣性恢復(fù)動態(tài) [J]. 應(yīng)用與環(huán)境生物學(xué)報, 27(4): 930-937.]

    MCCLAIN CR, REX MA, 2015. Toward a conceptual understanding of β-diversity in the deep-sea benthos [J]. Ann Rev Ecol Evol Syst, 46: 623-642.

    MI XC, WANG XG, SHEN GC, et al., 2022. Chinese forest biodiversity monitoring network (CForBio): Twenty years of exploring community assembly mechanisms and prospects for future research [J]. Biodivers Sci, 30(10): 211-233." [米湘成, 王緒高, 沈國春, 等, 2022. 中國森林生物多樣性監(jiān)測網(wǎng)絡(luò): 二十年群落構(gòu)建機(jī)制探索的回顧與展望 [J]. 生物多樣性, 30(10): 211-233.]

    MUGNAI M, TRINDADE DPF, THIERRY M, et al., 2022. Environment and space drive the community assembly of Atlantic European grasslands: Insights from multiple facets [J]. J Biogeogr, 49(4): 699-711.

    NILSSON MC, WARDLE DA, 2005. Understory vegetation as a forest ecosystem driver: evidence from the northern Swedish boreal forest [J]. Front Ecol Environ, 3(8): 421-428.

    NIU KC, LIU YN, SHEN ZH, et al., 2009. Community assembly: the relative importance of neutral theory and niche theory [J]. Biodivers Sci, 17(6): 579-593." [??瞬?, 劉懌寧, 沈澤昊, 等, 2009. 群落構(gòu)建的中性理論和生態(tài)位理論 [J]. 生物多樣性, 17(6): 579-593.]

    OKSANEN J, BLANCHET FG, KINDT R, et al., 2019. Vegan: community ecology package [CP]. R package version 2.5-4. https://CRAN.R-project.org/package=vegan.

    PEGUERO G, FERRN M, SARDANS J, et al., 2022. Decay of similarity across tropical forest communities: integrating spatial distance with soil nutrients [J]. Ecology, 103(2): e03599.

    QIAN H, QIAN SH, 2023. Geographic patterns of taxonomic and phylogenetic β-diversity of angiosperm genera in regional floras across the world [J]. Plant Divers, 45: 491-500.

    QIAN H, RICKLEFS RE, WHITE PS, 2005. Beta diversity of angiosperms in temperate floras of eastern Asia and eastern North America [J]. Ecol Lett, 8(1): 15-22.

    QU MJ, ABABAIKE N, ZOU XG, et al., 2022. Influence of geographic distance and environmental factors on beta diversity of plants in the Alxa Gobi region in northern China [J]. Biodivers Sci, 30(11): 22029." [曲夢君, 努爾依拉·阿巴拜克, 鄒旭閣, 等, 2022. 地理距離和環(huán)境因子對阿拉善戈壁植物群落β多樣性的影響 [J]. 生物多樣性, 30(11): 22029.]

    SCHULZ K, GUSCHAL M, KOWARIK I, et al., 2019. Grazing reduces plant species diversity of Caatinga dry forests in northeastern Brazil [J]. Appl Veg Sci, 22(2): 348-359.

    SHA H, LIANG YL, PENG N, et al., 2024. Effect of forest gap size on species diversity of vascular plants in shrub and herb layer of Picea crassifolia in Helan Mountains [J]. Chin J Ecol, 43(8): 2382-2394." [沙歡, 梁詠亮, 彭妞, 等, 2024. 林窗大小對青海云杉林灌草層維管植物物種多樣性的影響 [J]. 生態(tài)學(xué)雜志, 43(8): 2382-2394.]

    SHI JJ, ZHAO MF, WANG YH, et al., 2019. Community assembly of herbaceous layer of the planted forests in the central Loess Plateau, China [J]. Chin J Plant Ecol, 43(9): 834-842." [施晶晶, 趙鳴飛, 王宇航, 等, 2019. 黃土高原腹地人工林下草本層群落構(gòu)建機(jī)制 [J]. 植物生態(tài)學(xué)報, 43(9): 834-842.]

    SHI W, WANG YQ, XIANG WS, et al., 2021. Environmental filtering and dispersal limitation jointly shaped the taxonomic and phylogenetic beta diversity of natural forests in southern China [J]. Ecol Evol, 11(13): 8783-8794.

    SI XF, ZHAO YH, CHEN CW, et al., 2017. Beta-diversity partitioning: methods, applications and perspectives [J]. Biodivers Sci, 25(5): 464-480." [斯幸峰, 趙郁豪, 陳傳武, 等, 2017. Beta多樣性分解: 方法、應(yīng)用與展望 [J]. 生物多樣性, 25(5): 464-480.]

    SINGH A, WAGNER B, KASEL S, et al., 2023. Canopy composition and spatial configuration influences beta diversity in temperate regrowth forests of southeastern Australia [J]. Drones, 7(3): 155.

    SOININEN J, HEINO J, WANG J, 2018. A meta-analysis of nestedness and turnover components of beta diversity across organisms and ecosystems [J]. Global Ecol Biogeogr, 27: 96-109.

    SOININEN J, MCDONALD R, HILLEBRAND H, 2007. The distance decay of similarity in ecological communities [J]. Ecography, 30(1): 3-12.

    TAN SS, YE ZL, YUAN LB, at al., 2013. Beta diversity of plant communities in Baishanzu Nature Reserve [J]. Acta Ecol Sin, 33(21): 6944-6956. [譚珊珊, 葉珍林, 袁留斌,等, 2013. 百山祖自然保護(hù)區(qū)植物群落Beta多樣性 [J]. 生態(tài)學(xué)報, 33(21): 6944-6956.]

    TUOMISTO H, RUOKOLAINEN K, YLI-HALLA M, 2003. Dispersal, environment, and floristic variation of western Amazonian forests [J]. Science, 299(5604): 241-244.

    WANG D, WANG XA, GUO H, et al., 2013. Effect of species dispersal and environmental factors on species assemblages in grassland communities [J]. Acta Ecol Sin, 33(14): 4409-4415." [王丹, 王孝安, 郭華, 等, 2013. 環(huán)境和擴(kuò)散對草地群落構(gòu)建的影響 [J]. 生態(tài)學(xué)報, 33(14): 4409-4415.]

    WANG JM, QU MJ, WANG Y, et al., 2022. The drivers of plant taxonomic, functional, and phylogenetic β-diversity in the gobi desert of northern Qinghai-Tibet Plateau [J]. Biodivers Sci, 30(6): 21503." [王健銘, 曲夢君, 王寅, 等, 2022. 青藏高原北部戈壁植物群落物種、功能與系統(tǒng)發(fā)育β多樣性分布格局及其影響因素 [J]. 生物多樣性, 30(6): 21503.]

    WANG JM, WANG Y, LI MX, et al., 2021. Divergent roles of environmental and spatial factors in shaping plant β-diversity of different growth forms in drylands [J]. Global Ecol Conserv, 26: e01487.

    WANG SX, GUO H, WANG XA, et al., 2013. Dispersal limitation versus environment filtering in the assembly of plant communities in the Ziwu Mountains [J]. Chin Agric Sci, 46(22): 4733-4744. [王世雄, 郭華, 王孝安, 等, 2013. 擴(kuò)散限制和環(huán)境篩選對子午嶺森林群落構(gòu)建的相對貢獻(xiàn) [J]. 中國農(nóng)業(yè)科學(xué), 46(22): 4733-4744.]

    WANG SX, XIA TT, WANG XA, 2021. Understory β diversity differences and influential factors between Pinus tabulaeformis plantation and natural Quercus wutaishanica forest on Loess Plateau [J]. Guihaia, 41(3): 362-371." [王世雄, 夏婷婷, 王孝安, 2021. 黃土高原油松和遼東櫟林下植物β多樣性差異及影響因素 [J]. 廣西植物, 41(3): 362-371.]

    WENG CL, ZHANG TT, WU DH, et al., 2019. Drivers and patterns of α- and β-diversity in ten main forest community types in Gutianshan, eastern China [J]. Biodivers Sci, 27(1): 33-41. [翁昌露, 張?zhí)锾铮?巫東豪, 等, 2019. 古田山10種主要森林群落類型的α和β多樣性格局及影響因素 [J]. 生物多樣性, 27(1): 33-41.]

    WHITTAKER RH, 1972. Evolution and measurement of species diversity [J]. Taxon, 21(2/3): 213-251.

    XU W, GONZALEZ-RODRIGUEZ MA, LI Z, et al., 2022. Effects of edaphic factors at different depths on β-diversity patterns for subtropical plant communities based on MS-GDM in southern China [J]. Forests, 13(12): 2184.

    YANG H, WANG Y, WANG JM, et al., 2021. Effects of environmental filtering and dispersal limitation on the β-diversity of plant communities in the south fringe of Kumtag Desert [J]. J Desert Res, 41(3): 147-154." [楊歡, 王寅, 王健銘, 等, 2021. 環(huán)境過濾和擴(kuò)散限制對庫姆塔格沙漠南緣植物群落β-多樣性的影響 [J]. 中國沙漠, 41(3): 147-154.]

    YANG J, LU MM, CAO M, et al., 2014. Phylogenetic and functional alpha and beta diversity in mid-mountain humid evergreen broad-leaved forest [J]. Chin Sci Bull, 59(24): 2349-2358. [楊潔, 盧孟孟, 曹敏, 等, 2014. 中山濕性常綠闊葉林系統(tǒng)發(fā)育和功能性狀的α及β多樣性 [J]. 科學(xué)通報, 59(24): 2349-2358.]

    YANG X, YAO ZL, WANG B, et al., 2023. Driving effects of forest stand structure of a subtropical evergreen broad-leaved forest on species composition variation: From local to regional scales [J]. Biodivers Sci, 31(2): 22139." [楊欣, 姚志良, 王彬, 等, 2023. 亞熱帶常綠闊葉林林分結(jié)構(gòu)對物種組成變異的驅(qū)動作用: 從局域到區(qū)域尺度 [J]. 生物多樣性, 31(2): 22139.]

    YANG XQ, TIAN GH, 2012. Research on diversity of understory vegetation of main forest community in Saihanba [J]. Hebei J For Orchard Res, 27(1): 10-13. [楊曉勤, 田國恒, 2012. 塞罕壩主要森林類型林下植被多樣性研究 [J]. 河北林果研究, 27(1): 10-13.]

    YOU YM, XU JY, CAI DX, et al., 2016. Environmental factors affecting plant species diversity of understory plant communities in a Castanopsis hystrix plantation chronosequence in Pingxiang, Guangxi, China [J]. Acta Ecol Sin, 36(1): 164-172." [尤業(yè)明, 徐佳玉, 蔡道雄, 等, 2016. 廣西憑祥不同年齡紅椎林林下植物物種多樣性及其環(huán)境解釋 [J]. 生態(tài)學(xué)報, 36(1): 164-172.]

    ZHANG L, GUO JR, TANG YL, et al., 2009. Roles of natural forest in eco-construction of Henan Province [J]. J Henan For Sci Technol, 29(3): 62-64. [張莉, 郭際榮, 湯雨露, 等, 2009. 論河南省天然林在生態(tài)體系建設(shè)中的作用 [J]. 河南林業(yè)科技, 29(3): 62-64.]

    ZHANG YC, 2023. Understory plant species diversity and its influencing factors in Betula platyphylla natural forest at different latitudes in the Greater Khingan Mountains [D]. Hohhot: Inner Mongolia Agricultural University: 37-42." [張軼超, 2023. 大興安嶺林區(qū)不同緯度白樺天然林林下植物物種多樣性及其影響因素研究 [D]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學(xué): 37-42.]

    ZHAO MF, WANG GY, XING KX, et al., 2017. Patterns and determinants of species similarity decay of forest communities in the western Qinling Mountains [J]. Biodivers Sci, 25(1): 3-10. [趙鳴飛, 王國義, 邢開雄, 等, 2017. 秦嶺西部森林群落相似性遞減格局及其影響因素 [J]. 生物多樣性, 25(1): 3-10.]

    ZHAO Z, HE L, LI G, et al., 2023. Partitioning beta diversity of dry and hot valley vegetation in the Nujiang River in Southwest China [J]. Front Ecol Evol, 11: 1199874.

    (責(zé)任編輯 周翠鳴)

    猜你喜歡
    多樣性河南省
    河南省樹人教育交流中心
    河南省樹人教育交流中心
    河南省樹人教育交流中心
    河南省樹人教育交流中心
    河南省樹人教育交流中心
    河南省樹人教育交流中心
    由古典戲曲看“代言體”在中國的前世今生
    戲劇之家(2016年22期)2016-11-30 15:13:39
    淺談新時期群文輔導(dǎo)工作的特征
    新時期群文輔導(dǎo)工作的特征
    海洋微生物次生代謝的生物合成機(jī)制
    科技資訊(2016年19期)2016-11-15 10:39:12
    欧美 亚洲 国产 日韩一| 禁无遮挡网站| 美国免费a级毛片| 成年版毛片免费区| 国产高清videossex| 亚洲自偷自拍图片 自拍| 在线十欧美十亚洲十日本专区| 欧美黄色淫秽网站| 色哟哟哟哟哟哟| 国产成人精品久久二区二区免费| 久热这里只有精品99| 1024香蕉在线观看| 99国产精品一区二区三区| 欧美成人免费av一区二区三区| 亚洲精品国产一区二区精华液| 久久久久国产一级毛片高清牌| 久久久国产成人精品二区| 欧美日韩一级在线毛片| 欧美日韩福利视频一区二区| 99久久国产精品久久久| 欧美色欧美亚洲另类二区 | 色婷婷久久久亚洲欧美| 欧美绝顶高潮抽搐喷水| 欧美+亚洲+日韩+国产| 美女高潮喷水抽搐中文字幕| 18美女黄网站色大片免费观看| 性色av乱码一区二区三区2| 女警被强在线播放| 黄色a级毛片大全视频| 欧美激情极品国产一区二区三区| 日本一区二区免费在线视频| 色播亚洲综合网| a级毛片在线看网站| 女警被强在线播放| 色av中文字幕| 国产1区2区3区精品| 国产一级毛片七仙女欲春2 | 欧美成人午夜精品| 老熟妇仑乱视频hdxx| 一区二区三区高清视频在线| 精品人妻1区二区| avwww免费| 亚洲av美国av| 变态另类丝袜制服| 在线免费观看的www视频| 午夜亚洲福利在线播放| 别揉我奶头~嗯~啊~动态视频| 国内精品久久久久久久电影| 国产单亲对白刺激| 精品欧美一区二区三区在线| 国产成人一区二区三区免费视频网站| 日韩有码中文字幕| 成人亚洲精品av一区二区| 韩国av一区二区三区四区| 日本黄色视频三级网站网址| 免费在线观看影片大全网站| 日本欧美视频一区| 亚洲av成人av| 老汉色av国产亚洲站长工具| 国产一区在线观看成人免费| 视频区欧美日本亚洲| 搡老妇女老女人老熟妇| 欧美日韩福利视频一区二区| 亚洲人成电影免费在线| 亚洲av电影不卡..在线观看| 国产亚洲精品综合一区在线观看 | tocl精华| 亚洲人成电影免费在线| 成年人黄色毛片网站| 日本vs欧美在线观看视频| 亚洲男人的天堂狠狠| 免费在线观看日本一区| 无限看片的www在线观看| 9热在线视频观看99| 亚洲国产高清在线一区二区三 | 9热在线视频观看99| 久久婷婷成人综合色麻豆| 亚洲精品久久成人aⅴ小说| 欧美激情高清一区二区三区| 久久天躁狠狠躁夜夜2o2o| 中亚洲国语对白在线视频| 青草久久国产| 嫩草影院精品99| 叶爱在线成人免费视频播放| 亚洲狠狠婷婷综合久久图片| av在线天堂中文字幕| 人成视频在线观看免费观看| 丝袜人妻中文字幕| √禁漫天堂资源中文www| svipshipincom国产片| 日本在线视频免费播放| x7x7x7水蜜桃| 久久人人爽av亚洲精品天堂| 亚洲专区字幕在线| 日本黄色视频三级网站网址| 日韩精品青青久久久久久| 不卡一级毛片| 变态另类成人亚洲欧美熟女 | 久久香蕉精品热| 色哟哟哟哟哟哟| 亚洲av成人av| 香蕉国产在线看| 色av中文字幕| 一区二区三区国产精品乱码| 悠悠久久av| 国产一区在线观看成人免费| 黄片小视频在线播放| 性色av乱码一区二区三区2| bbb黄色大片| 午夜福利视频1000在线观看 | 国产成人系列免费观看| 婷婷精品国产亚洲av在线| 国产精品二区激情视频| 精品久久久久久久毛片微露脸| 两性夫妻黄色片| 久久久久久亚洲精品国产蜜桃av| 69av精品久久久久久| 可以在线观看的亚洲视频| 人妻久久中文字幕网| 久久热在线av| 好看av亚洲va欧美ⅴa在| 12—13女人毛片做爰片一| 久久精品aⅴ一区二区三区四区| 久久久久久国产a免费观看| 90打野战视频偷拍视频| 久久人妻福利社区极品人妻图片| 国产欧美日韩精品亚洲av| 老鸭窝网址在线观看| 国产精品 国内视频| 亚洲人成电影观看| 别揉我奶头~嗯~啊~动态视频| 多毛熟女@视频| 欧美乱色亚洲激情| 91成人精品电影| 欧美激情 高清一区二区三区| 日日干狠狠操夜夜爽| 国产成人一区二区三区免费视频网站| 少妇的丰满在线观看| 黑人欧美特级aaaaaa片| 亚洲欧美激情在线| 美女午夜性视频免费| 久久香蕉国产精品| 欧美激情 高清一区二区三区| 亚洲精品久久国产高清桃花| 可以免费在线观看a视频的电影网站| 亚洲天堂国产精品一区在线| av免费在线观看网站| 久久久久亚洲av毛片大全| 19禁男女啪啪无遮挡网站| 免费在线观看日本一区| 欧美精品啪啪一区二区三区| 国产男靠女视频免费网站| 亚洲av电影在线进入| 性色av乱码一区二区三区2| 国产高清有码在线观看视频 | 亚洲午夜精品一区,二区,三区| 亚洲va日本ⅴa欧美va伊人久久| 天天躁夜夜躁狠狠躁躁| 久久天堂一区二区三区四区| 亚洲狠狠婷婷综合久久图片| av福利片在线| 黄色毛片三级朝国网站| 身体一侧抽搐| 露出奶头的视频| a级毛片在线看网站| x7x7x7水蜜桃| 亚洲精品av麻豆狂野| 国产午夜福利久久久久久| 波多野结衣av一区二区av| 亚洲精品美女久久av网站| 日本 欧美在线| 国产精品久久电影中文字幕| 女性被躁到高潮视频| 日韩 欧美 亚洲 中文字幕| 亚洲久久久国产精品| 日韩高清综合在线| 欧美中文综合在线视频| 国产伦一二天堂av在线观看| 涩涩av久久男人的天堂| 国产午夜福利久久久久久| 欧美色欧美亚洲另类二区 | 免费观看精品视频网站| 一夜夜www| 50天的宝宝边吃奶边哭怎么回事| 九色亚洲精品在线播放| 女人高潮潮喷娇喘18禁视频| 亚洲熟妇熟女久久| 亚洲,欧美精品.| 两个人视频免费观看高清| 亚洲精品久久国产高清桃花| 十八禁人妻一区二区| 欧美日韩亚洲综合一区二区三区_| 91麻豆精品激情在线观看国产| 男女之事视频高清在线观看| 成年人黄色毛片网站| 亚洲人成电影免费在线| 免费在线观看亚洲国产| 99精品久久久久人妻精品| av在线天堂中文字幕| 色播亚洲综合网| 午夜精品在线福利| 欧美一级a爱片免费观看看 | 久久午夜亚洲精品久久| 亚洲av电影不卡..在线观看| 美国免费a级毛片| 日韩成人在线观看一区二区三区| 亚洲无线在线观看| 91国产中文字幕| 两个人看的免费小视频| 亚洲欧美精品综合一区二区三区| ponron亚洲| 麻豆成人av在线观看| 午夜激情av网站| 国产不卡一卡二| 中亚洲国语对白在线视频| 免费在线观看黄色视频的| 色尼玛亚洲综合影院| 亚洲国产欧美网| 777久久人妻少妇嫩草av网站| 国产xxxxx性猛交| 国产精品久久久久久精品电影 | 69av精品久久久久久| 人妻久久中文字幕网| 怎么达到女性高潮| 久久欧美精品欧美久久欧美| 一个人免费在线观看的高清视频| 日本撒尿小便嘘嘘汇集6| 欧美黄色淫秽网站| 欧美最黄视频在线播放免费| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久久久久精品电影 | 日韩 欧美 亚洲 中文字幕| 宅男免费午夜| 高清黄色对白视频在线免费看| 国产亚洲欧美精品永久| 午夜福利影视在线免费观看| 亚洲美女黄片视频| 宅男免费午夜| 亚洲九九香蕉| 18美女黄网站色大片免费观看| 亚洲人成电影免费在线| av天堂久久9| 少妇 在线观看| 国产一区在线观看成人免费| 满18在线观看网站| av片东京热男人的天堂| 50天的宝宝边吃奶边哭怎么回事| 99国产精品一区二区三区| 久久久久久免费高清国产稀缺| 国产高清激情床上av| 色播亚洲综合网| 国产精品亚洲av一区麻豆| 国产激情欧美一区二区| 黄色 视频免费看| 成人亚洲精品一区在线观看| 成年女人毛片免费观看观看9| 国产精品香港三级国产av潘金莲| 男女做爰动态图高潮gif福利片 | 精品熟女少妇八av免费久了| 人人澡人人妻人| 老司机深夜福利视频在线观看| 精品国产超薄肉色丝袜足j| 日韩欧美国产一区二区入口| 亚洲少妇的诱惑av| 国产男靠女视频免费网站| 色综合站精品国产| 亚洲人成伊人成综合网2020| 伦理电影免费视频| 一区二区三区激情视频| 在线观看66精品国产| 午夜精品国产一区二区电影| 亚洲av美国av| 亚洲国产欧美一区二区综合| 级片在线观看| 大型黄色视频在线免费观看| 国产一区二区三区综合在线观看| 香蕉国产在线看| 丰满人妻熟妇乱又伦精品不卡| 人人澡人人妻人| 午夜激情av网站| 在线天堂中文资源库| 国产精品亚洲av一区麻豆| 日韩大尺度精品在线看网址 | 少妇被粗大的猛进出69影院| 久久中文字幕一级| 一区二区三区精品91| 亚洲欧洲精品一区二区精品久久久| 丁香六月欧美| 黄色毛片三级朝国网站| 国产午夜福利久久久久久| 久久国产精品男人的天堂亚洲| 亚洲性夜色夜夜综合| av福利片在线| 日韩高清综合在线| 法律面前人人平等表现在哪些方面| 99国产精品一区二区三区| 一进一出抽搐gif免费好疼| 成在线人永久免费视频| 狂野欧美激情性xxxx| 亚洲欧美精品综合一区二区三区| 精品乱码久久久久久99久播| 亚洲色图 男人天堂 中文字幕| 国产精品电影一区二区三区| 久久国产乱子伦精品免费另类| 免费人成视频x8x8入口观看| 久久午夜综合久久蜜桃| 悠悠久久av| or卡值多少钱| 亚洲国产中文字幕在线视频| 国语自产精品视频在线第100页| 亚洲成av片中文字幕在线观看| 如日韩欧美国产精品一区二区三区| 窝窝影院91人妻| 99久久久亚洲精品蜜臀av| а√天堂www在线а√下载| 黑人操中国人逼视频| 日韩欧美三级三区| 黄色丝袜av网址大全| 99re在线观看精品视频| 午夜福利成人在线免费观看| 男人舔女人下体高潮全视频| 人妻久久中文字幕网| 99国产极品粉嫩在线观看| 国产色视频综合| 精品福利观看| 欧美成人性av电影在线观看| 色婷婷久久久亚洲欧美| 日本免费a在线| 丝袜美足系列| 亚洲无线在线观看| 99久久99久久久精品蜜桃| 中文字幕人妻丝袜一区二区| 高清在线国产一区| 国产高清激情床上av| 成人三级做爰电影| 国产亚洲精品第一综合不卡| 精品一区二区三区四区五区乱码| 一区二区三区高清视频在线| 可以在线观看的亚洲视频| 丝袜在线中文字幕| 制服诱惑二区| 老汉色∧v一级毛片| 视频在线观看一区二区三区| 国产视频一区二区在线看| 99国产精品一区二区蜜桃av| 欧美日韩精品网址| 国产成人精品无人区| 搞女人的毛片| av视频在线观看入口| 久久精品亚洲熟妇少妇任你| 欧美日韩乱码在线| 少妇熟女aⅴ在线视频| 久久中文字幕一级| 香蕉久久夜色| 丝袜美腿诱惑在线| 99久久99久久久精品蜜桃| 久久久水蜜桃国产精品网| 国产人伦9x9x在线观看| 亚洲精品一区av在线观看| 女警被强在线播放| 999久久久国产精品视频| 国产av又大| 精品少妇一区二区三区视频日本电影| cao死你这个sao货| 激情在线观看视频在线高清| 国产成人精品久久二区二区91| 国产精品亚洲美女久久久| 国产av精品麻豆| netflix在线观看网站| 亚洲精品美女久久av网站| 免费在线观看日本一区| 国产午夜精品久久久久久| 亚洲色图综合在线观看| videosex国产| 国产在线精品亚洲第一网站| 91成人精品电影| 亚洲 欧美一区二区三区| 亚洲人成网站在线播放欧美日韩| 天天添夜夜摸| 色尼玛亚洲综合影院| 神马国产精品三级电影在线观看 | svipshipincom国产片| 一级毛片高清免费大全| 国产野战对白在线观看| 视频区欧美日本亚洲| 久久草成人影院| 久久久久久久久免费视频了| 变态另类成人亚洲欧美熟女 | xxx96com| 自线自在国产av| 两个人视频免费观看高清| 国产日韩一区二区三区精品不卡| 亚洲精品一区av在线观看| 欧美成人性av电影在线观看| 亚洲国产中文字幕在线视频| 欧美一区二区精品小视频在线| 人成视频在线观看免费观看| 中文字幕最新亚洲高清| 亚洲精品在线观看二区| 男女午夜视频在线观看| 90打野战视频偷拍视频| 国产aⅴ精品一区二区三区波| 色在线成人网| 久久香蕉激情| 女性被躁到高潮视频| 国产成人精品在线电影| 黑人欧美特级aaaaaa片| 巨乳人妻的诱惑在线观看| 国产精品秋霞免费鲁丝片| 精品一区二区三区视频在线观看免费| 69av精品久久久久久| 精品国产国语对白av| 日本精品一区二区三区蜜桃| 久久草成人影院| 乱人伦中国视频| 真人一进一出gif抽搐免费| 岛国在线观看网站| 美女高潮喷水抽搐中文字幕| 亚洲欧美激情综合另类| 久久欧美精品欧美久久欧美| 精品久久蜜臀av无| 国产精品久久久人人做人人爽| 午夜久久久在线观看| 少妇被粗大的猛进出69影院| 成人18禁高潮啪啪吃奶动态图| 人人妻人人爽人人添夜夜欢视频| 中文字幕高清在线视频| 久久中文字幕人妻熟女| 9热在线视频观看99| 这个男人来自地球电影免费观看| 动漫黄色视频在线观看| 久久久国产成人免费| 国产精品免费视频内射| 大香蕉久久成人网| 女人被躁到高潮嗷嗷叫费观| av片东京热男人的天堂| aaaaa片日本免费| 日本五十路高清| 亚洲成a人片在线一区二区| 国产精品久久久av美女十八| 法律面前人人平等表现在哪些方面| 欧美老熟妇乱子伦牲交| 国产激情欧美一区二区| 欧美日韩乱码在线| 麻豆成人av在线观看| 精品欧美国产一区二区三| www.自偷自拍.com| e午夜精品久久久久久久| 亚洲成人精品中文字幕电影| 欧美国产日韩亚洲一区| 久久精品影院6| netflix在线观看网站| 香蕉丝袜av| 成人国产综合亚洲| 欧美乱妇无乱码| 韩国精品一区二区三区| 国产亚洲欧美98| 久久狼人影院| 精品国产乱子伦一区二区三区| 国产亚洲欧美在线一区二区| 久久精品国产综合久久久| 极品教师在线免费播放| 老熟妇乱子伦视频在线观看| 青草久久国产| 精品电影一区二区在线| 国产成人欧美在线观看| 欧美午夜高清在线| 午夜免费鲁丝| 精品国产乱码久久久久久男人| 色婷婷久久久亚洲欧美| 成人手机av| 97碰自拍视频| 国产成人啪精品午夜网站| 亚洲成人国产一区在线观看| 可以在线观看毛片的网站| 中文字幕另类日韩欧美亚洲嫩草| 久久国产精品影院| 制服诱惑二区| 亚洲自拍偷在线| 黄色丝袜av网址大全| 亚洲精品国产精品久久久不卡| 午夜日韩欧美国产| 性色av乱码一区二区三区2| 免费无遮挡裸体视频| 制服诱惑二区| 亚洲欧美日韩无卡精品| 久久人妻av系列| 亚洲中文av在线| 亚洲自偷自拍图片 自拍| 人成视频在线观看免费观看| 丝袜在线中文字幕| 亚洲第一欧美日韩一区二区三区| 欧美黄色片欧美黄色片| 亚洲中文日韩欧美视频| 两个人视频免费观看高清| 精品国产亚洲在线| 国内精品久久久久久久电影| 国产av一区二区精品久久| 亚洲欧美日韩高清在线视频| 久久中文字幕一级| 亚洲欧美精品综合一区二区三区| 一区在线观看完整版| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品国产色婷婷电影| 久久久久久国产a免费观看| 超碰成人久久| 丝袜人妻中文字幕| 久久久久亚洲av毛片大全| 免费久久久久久久精品成人欧美视频| 88av欧美| 老司机深夜福利视频在线观看| 久久久久九九精品影院| 男女午夜视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲久久久国产精品| 亚洲一区二区三区不卡视频| 一区二区三区精品91| 国产男靠女视频免费网站| 亚洲精品国产色婷婷电影| 亚洲无线在线观看| 一二三四在线观看免费中文在| 久久午夜亚洲精品久久| 97人妻精品一区二区三区麻豆 | 热re99久久国产66热| 99国产精品一区二区蜜桃av| 一区二区三区高清视频在线| 欧美日韩亚洲国产一区二区在线观看| 日韩三级视频一区二区三区| 好看av亚洲va欧美ⅴa在| 亚洲精品久久国产高清桃花| 久久久久亚洲av毛片大全| 高清在线国产一区| 久久久久久久午夜电影| 欧美日韩精品网址| 亚洲一卡2卡3卡4卡5卡精品中文| 精品久久久久久,| 黑人欧美特级aaaaaa片| 少妇裸体淫交视频免费看高清 | 黄片播放在线免费| 国产亚洲精品久久久久5区| 在线观看66精品国产| 久久这里只有精品19| 黄色女人牲交| 欧美日本视频| 99riav亚洲国产免费| 国产三级在线视频| 久久国产精品人妻蜜桃| 国产成人精品在线电影| 亚洲国产精品合色在线| 大型av网站在线播放| 国产精品乱码一区二三区的特点 | 99国产精品99久久久久| 淫秽高清视频在线观看| 欧美绝顶高潮抽搐喷水| 丝袜美腿诱惑在线| 精品高清国产在线一区| 天堂√8在线中文| 一级毛片高清免费大全| 久久中文字幕人妻熟女| 色婷婷久久久亚洲欧美| a级毛片在线看网站| 侵犯人妻中文字幕一二三四区| 91国产中文字幕| 久久久久久人人人人人| aaaaa片日本免费| 久久国产精品男人的天堂亚洲| 91麻豆av在线| 日韩精品免费视频一区二区三区| 国产精品久久电影中文字幕| 国产一卡二卡三卡精品| 黑人欧美特级aaaaaa片| av有码第一页| 欧美乱码精品一区二区三区| 免费观看人在逋| 黑人操中国人逼视频| 日韩欧美免费精品| 久久精品亚洲精品国产色婷小说| 亚洲aⅴ乱码一区二区在线播放 | 精品无人区乱码1区二区| 美女免费视频网站| 黄色丝袜av网址大全| 在线天堂中文资源库| 亚洲专区字幕在线| 91国产中文字幕| 精品国产一区二区三区四区第35| 一级毛片高清免费大全| 十八禁人妻一区二区| 国产三级在线视频| 啦啦啦观看免费观看视频高清 | 成人av一区二区三区在线看| 午夜久久久在线观看| 精品一区二区三区四区五区乱码| 亚洲精品久久成人aⅴ小说| 亚洲午夜理论影院| 国产亚洲av嫩草精品影院| 国产又色又爽无遮挡免费看| 色播在线永久视频| 18禁裸乳无遮挡免费网站照片 | 午夜福利,免费看| 亚洲一区中文字幕在线| 国产精品美女特级片免费视频播放器 | 免费女性裸体啪啪无遮挡网站| 色哟哟哟哟哟哟| 亚洲国产精品久久男人天堂| 久久人妻福利社区极品人妻图片| 国产视频一区二区在线看| 操美女的视频在线观看| 男人舔女人下体高潮全视频| 一级毛片精品| 久久亚洲真实| 亚洲自拍偷在线| 久久国产亚洲av麻豆专区| 麻豆av在线久日| av中文乱码字幕在线| 精品国内亚洲2022精品成人|