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

    長(zhǎng)時(shí)間尺度草地物種多樣性-群落穩(wěn)定性關(guān)系及機(jī)制

    2025-05-30 00:00:00祁旭陽(yáng)王欣雨尹子鳴馬曉瑞楊夢(mèng)茹王子臨王亦波田慧慧拓行行張?zhí)碛?/span>李偉
    草地學(xué)報(bào) 2025年4期
    關(guān)鍵詞:物種效應(yīng)

    關(guān)鍵詞:時(shí)間尺度;物種多樣性;群落穩(wěn)定性;異步性;種群穩(wěn)定性;統(tǒng)計(jì)平均效應(yīng)中圖分類(lèi)號(hào): .2 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1007-0435(2025)04-1239-10

    Abstract:Revealing the paterns,relationships,and potential mechanisms ofchanges in grassland species diversity and community stability on long time scales is one of the core objectives of ecological research.Based on a 40-year long-term monitoring experiment in a typical steppe in Yunwu Mountain,Ningxia,this study used Patrick's species richnessindex,Shannon-Wiener diversity index,and other indicators to measure species diversity.By applying mathematical decomposition methods,community stability was broken down into population stability,asynchrony,compensatory efects,and statistical averaging effects.Using generalinear models and random forest models,the relationships and underlying mechanisms between species diversity and community stability were explored.The results showed that from 1982 to 2O21,both grassand species diversity and community stability exhibited a“unimodal\" patern ofchange. Species diversity was predominantly positively correlated with community stability.The diferent mechanisms were allpositively corelated with community stability,with their explanatory power for community stability variation ranked as follows:population stability ( 9 4 % ) gt; asynchrony ( 8 6 % ) gt; statistical averaging effect ( 6 4 % ) compensatory effect ( 4 6 % ) . In summary,allfour mechanisms contributed to the maintenance of community stability,but the statistical averaging effect was more critical than the compensatory effect.

    Key words: Time scale; Species diversity; Community stability; Asynchrony; Population stability; Statisticalaveraging effect

    人類(lèi)活動(dòng)和氣候變化引起的全球物種多樣性的持續(xù)減少顯著改變了草地生態(tài)系統(tǒng)的結(jié)構(gòu)和功能,進(jìn)而降低了草地生態(tài)系統(tǒng)為人類(lèi)提供服務(wù)的能力[1-3,而生態(tài)系統(tǒng)維持其功能和服務(wù)的能力與生態(tài)系統(tǒng)的穩(wěn)定性密切相關(guān)[4-5]。經(jīng)典生態(tài)演替理論認(rèn)為,演替時(shí)間對(duì)物種更迭和種群波動(dòng)具有重要影響[6-7],但長(zhǎng)時(shí)間尺度上物種多樣性-群落穩(wěn)定性的關(guān)系目前仍存在較大爭(zhēng)議。Elton81958年提出高物種多樣性和更復(fù)雜的群落結(jié)構(gòu)會(huì)顯著增強(qiáng)生態(tài)系統(tǒng)穩(wěn)定性。然而,May[9]利用數(shù)學(xué)模型通過(guò)抽象的簡(jiǎn)單系統(tǒng)來(lái)代替自然生態(tài)系統(tǒng),并將物種的相互作用強(qiáng)度隨機(jī)分布,發(fā)現(xiàn)生態(tài)系統(tǒng)穩(wěn)定性隨著物種多樣性增加而下降。Tilman等[10]在CedarCreek進(jìn)行的人工控制物種豐富度實(shí)驗(yàn)表明物種多樣性與群落穩(wěn)定性呈正相關(guān)關(guān)系。然而,現(xiàn)有研究仍有諸多不足之處。第一,無(wú)論是數(shù)學(xué)模型模擬還是人工控制物種多樣性實(shí)驗(yàn)均沒(méi)有以自然群落為研究對(duì)象,未充分考慮自然生態(tài)系統(tǒng)本身的復(fù)雜性[9-11]。第二,這些實(shí)驗(yàn)持續(xù)時(shí)間均較短,大多是 1 ~ 10年,所得出的結(jié)論很難推廣至更長(zhǎng)的時(shí)間尺度[10.12-13]。第三,近期的研究主要集中于外部環(huán)境改變對(duì)群落穩(wěn)定性的影響,例如放牧、氮添加、降水改變等[14-16],并得出異步性是維持群落穩(wěn)定的關(guān)鍵機(jī)制,但對(duì)于驅(qū)動(dòng)異步性的統(tǒng)計(jì)平均效應(yīng)或補(bǔ)償效應(yīng)動(dòng)態(tài)變化及其相對(duì)重要性仍不明確,維持群落穩(wěn)定性的關(guān)鍵機(jī)制仍有很大爭(zhēng)論[17-19]。因此,研究長(zhǎng)時(shí)間尺度下自然群落的穩(wěn)定性變化及維持機(jī)制具有重要意義。

    群落穩(wěn)定性可以分解為種群穩(wěn)定性和不同物種之間的同步性乘積,該研究結(jié)果有效的促進(jìn)了自然群落和理想群落研究結(jié)果的整合[20]。物種多樣性-種群穩(wěn)定性的關(guān)系可能是正相關(guān)或負(fù)相關(guān)[21],但由于異步性具有關(guān)鍵的正向作用,物種多樣性-群落穩(wěn)定性關(guān)系總體上為正相關(guān)[22]。然而,這種正向作用是生物學(xué)機(jī)制還是統(tǒng)計(jì)學(xué)機(jī)制目前仍沒(méi)有一致結(jié)論,可應(yīng)用數(shù)學(xué)分解方法將異步性進(jìn)一步分解為補(bǔ)償效應(yīng)和統(tǒng)計(jì)平均效應(yīng)來(lái)明確其內(nèi)在機(jī)制[23]。補(bǔ)償效應(yīng)是一種生物學(xué)機(jī)制,物種對(duì)環(huán)境變化的響應(yīng)不同導(dǎo)致不同物種之間的產(chǎn)生性狀差異,一些物種生產(chǎn)力下降會(huì)被其他生產(chǎn)力上升的物種所補(bǔ)償,從而維持整個(gè)群落的穩(wěn)定[24]。統(tǒng)計(jì)平均效應(yīng)描述了一個(gè)群落內(nèi)獨(dú)立物種或種群越多則整個(gè)群落越穩(wěn)定的現(xiàn)象,是一種統(tǒng)計(jì)學(xué)機(jī)制[23]。Cardinale等[25]分析了44個(gè)控制物種豐富度實(shí)驗(yàn)結(jié)果發(fā)現(xiàn)補(bǔ)償效應(yīng)的作用被嚴(yán)重低估,在生產(chǎn)力低于平均水平物種在群落占主導(dǎo)地位時(shí),補(bǔ)償效應(yīng)對(duì)生物量的貢獻(xiàn)是選擇效應(yīng)的1.9倍。內(nèi)蒙古草原24年的物種多樣性實(shí)驗(yàn)表明,種群或功能群水平的相互補(bǔ)償對(duì)維持群落穩(wěn)定起到了關(guān)鍵作用[24]。但相關(guān)理論和實(shí)驗(yàn)研究卻發(fā)現(xiàn),統(tǒng)計(jì)平均效應(yīng)可以對(duì)群落穩(wěn)定產(chǎn)生積極或消極影響,具有不容忽視的作用[26]。因此,辨析長(zhǎng)時(shí)間尺度下自然群落中補(bǔ)償效應(yīng)和統(tǒng)計(jì)平均效應(yīng)的相對(duì)重要性對(duì)深入理解草地物種多樣性-群落穩(wěn)定性關(guān)系至關(guān)重要。

    位于黃土高原西部的寧夏云霧山草原國(guó)家級(jí)自然保護(hù)區(qū)自1982年以來(lái)嚴(yán)格實(shí)施封山禁牧和補(bǔ)播優(yōu)良牧草等措施[27-28],有效促進(jìn)了草地的正向演替[29-30],是自前黃土高原保存最完整,最具代表性的長(zhǎng)芒草(Stipabungeana)天然草原[31]。本研究以云霧山典型草原生態(tài)系統(tǒng)為研究對(duì)象,基于40年長(zhǎng)期定位監(jiān)測(cè)實(shí)驗(yàn),深入分析長(zhǎng)時(shí)間尺度上物種多樣性和群落穩(wěn)定性的變化模式、相關(guān)關(guān)系及其維持機(jī)制,可為草地生態(tài)系統(tǒng)可持續(xù)健康管理提供理論依據(jù),也可為區(qū)域生物多樣性保護(hù)和草牧業(yè)發(fā)展提供重要參考。

    1 材料與方法

    1.1 研究區(qū)域概況

    研究區(qū)位于寧夏固原市云霧山國(guó)家級(jí)自然保護(hù)區(qū)內(nèi) ,依托西北農(nóng)林科技大學(xué)寧夏云霧山草原生態(tài)系統(tǒng)國(guó)家長(zhǎng)期科研基地。該區(qū)域海拔高度為 處于溫帶半干旱氣候區(qū),主要受溫帶大陸性季風(fēng)氣候的影響;年平均氣溫 月氣溫最高 ;1月氣溫最低 ;年大于零度積溫 ;年日照時(shí)長(zhǎng) ;年降水量 左右, 60 % ~ 7 5 % 的降水集中分布于夏季7一9月;山地灰褐土和黑瀘土為主要土壤類(lèi)型, 。植被類(lèi)型為暖溫帶典型草原。目前保護(hù)區(qū)內(nèi)探明野生物種297種,優(yōu)勢(shì)種和建群種主要有大針茅(Stipagrandis)、白蓮蒿(Arte-misiasacrorum)、長(zhǎng)芒草(S.bungeana)甘菊(Chrysanthemumlauandulifolium)甘青針茅(S.przewal-skyi)和百里香(Thymus mongolicus)等[33]。

    1. 2 試驗(yàn)設(shè)計(jì)

    自1982年保護(hù)區(qū)成立之初,實(shí)施封山禁牧措施。在保護(hù)區(qū)核心區(qū)內(nèi)隨機(jī)確定15個(gè) 的小區(qū)并利用鐵樁固定四角,用于長(zhǎng)期觀測(cè),小區(qū)之間距離 以上。每個(gè)小區(qū)分成兩部分:一部分為 ,用于植被調(diào)查;另一部分為 18m,用于植物個(gè)體采樣[34]

    1.3 研究方法

    1982一2021年每年8月中旬在每個(gè)定位小區(qū)的 部分,隨機(jī)選取1個(gè) 的樣方,調(diào)查每個(gè)樣方內(nèi)植物群落的物種豐富度、多度、高度、蓋度。為降低人工采樣誤差,樣方位置逐年更換;物種豐富度和物種多度以群落中每年出現(xiàn)的物種數(shù)及物種株叢數(shù)來(lái)衡量;每個(gè)物種小于5株全部測(cè)量,大于等于5株則隨機(jī)選取5株來(lái)測(cè)量高度,取其均值記錄為物種高度[34]。對(duì)樣方內(nèi)每個(gè)物種分別齊地劉割,裝人帶編號(hào)紙袋內(nèi)封口并帶回實(shí)驗(yàn)室。放入烘箱之內(nèi), 烘48小時(shí)至恒重,所得干重即為物種地上生物量。

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

    1.4.1物種多樣性指標(biāo)本研究采用Patrick物種豐富度指數(shù)、Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù)、Simpson指數(shù)作為物種多樣性的衡量指標(biāo)[35],計(jì)算公式為:

    Patrick物種豐富度指數(shù): (1-1)

    Shannon-Wiener多樣性指數(shù)

    Pielou均勻度指數(shù):

    Simpson指數(shù)

    式中, 為物種 i 的個(gè)體數(shù), N 為群落中所有物種總個(gè)體數(shù), S 為總物種數(shù)。

    1.4.2群落結(jié)構(gòu)根據(jù)植物的生長(zhǎng)型和生活型,將群落物種組成劃分為一年生或二年生(AB)多年生叢生禾草(PB)多年生根莖禾草(PR)多年生雜草(PF)和灌木或半灌木(SS)5個(gè)功能群,以此探究群落結(jié)構(gòu)的變化模式。

    1.4.3 群落穩(wěn)定性和效應(yīng)值根據(jù)Thibautamp;Connolly20]理論,群落穩(wěn)定性可分解為種群穩(wěn)定性和不同種群之間的異步性乘積,即

    群落穩(wěn)定性計(jì)算公式為[23]:

    式中: 表示在時(shí)間 時(shí)物種 …·S的生物量; , 分別表示物種i隨時(shí)間計(jì)算的地上生物量平均值,標(biāo)準(zhǔn)差,方差; 表示物種i,j隨時(shí)間變化的協(xié)方差動(dòng)態(tài);則 ;當(dāng)種群時(shí)間序列不相關(guān)時(shí), 為了在時(shí)間尺度上描述群落穩(wěn)定性變化,以?xún)赡隇闈L動(dòng)窗口(即1982—1983,1983—1984,1984—1985··2020—2021)分別計(jì)算群落穩(wěn)定性指數(shù)及相關(guān)指標(biāo)。

    種群穩(wěn)定性計(jì)算公式為:

    群落穩(wěn)定性與種群穩(wěn)定性的比率是Loreau-deMazancourt異步指數(shù) ,其計(jì)算公式為:

    假設(shè)一個(gè)理想群落物種水平方差和均值與原始群落相同,但物種協(xié)方差等于零,因此公式 ( 2 - 1)變?yōu)? ,在不相關(guān)或獨(dú)立種群的情況下, 是群落穩(wěn)定性的值。SAE是統(tǒng)計(jì)平均效應(yīng)中 Φ 的分量(\"ip”代表“獨(dú)立種群”),其計(jì)算公式為:

    補(bǔ)償效應(yīng)是 Φ 的其余部分,其計(jì)算公式為:

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

    首先建立長(zhǎng)時(shí)間序列變化模式,分別擬合一般線(xiàn)性模型和二次方模型,利用赤池信息準(zhǔn)則(Akaike'sinformationcriterion,AIC)選擇最優(yōu)模型,AIC越小說(shuō)明擬合程度越高;利用MOS(Mitch-ell-OldsandShaw)確定變化模式是否為典型“單峰狀”,綜合兩項(xiàng)分析選擇時(shí)間尺度上物種多樣性、群落穩(wěn)定性及其效應(yīng)值的最優(yōu)變化模式。其次,利用一元線(xiàn)性回歸法擬合物種多樣性和群落穩(wěn)定性及相關(guān)效應(yīng)值關(guān)系,通過(guò)隨機(jī)森林模型和方差分解方法進(jìn)一步量化不同穩(wěn)定性機(jī)制對(duì)群落穩(wěn)定性的解釋率。2007年和2014年因定位監(jiān)測(cè)小區(qū)圍欄內(nèi)闖入野生動(dòng)物,部分啃食了草地地上部分,因此我們剔除了這兩年的數(shù)據(jù)。所有統(tǒng)計(jì)分析采用R-4.3.3完成,統(tǒng)計(jì)檢驗(yàn)的概率顯著性水平為 。

    2 結(jié)果與分析

    2.1時(shí)間尺度上物種多樣性和群落結(jié)構(gòu)變化模式在長(zhǎng)時(shí)間尺度上,物種豐富度和Shannon

    Wiener多樣性指數(shù)隨時(shí)間變化模式為“單峰狀”;Pielou均勻度指數(shù)變化模式為“倒單峰狀”,Simpson指數(shù)變化模式為“波動(dòng)上升狀”(圖1,表1)。群落結(jié)構(gòu)早期主要以多年生叢生禾草和多年生雜草為主,占據(jù)整個(gè)群落 73 % 以上的生物量,但2013—2021出現(xiàn)多年生叢生禾草、多年生雜草、灌木或半灌木三個(gè)功能群生物量接近的變化(圖2A,B)。

    物種豐富度具體變化模式為“急劇上升 維持穩(wěn)定 緩慢下降”,變化范圍為 個(gè),平均為18.67個(gè),從起初封育時(shí)的平均10個(gè)物種急劇上升至1995年25個(gè),維持高物種豐富度的狀態(tài)到2012年,隨后急劇下降,到2021年平均物種數(shù)為6個(gè);Pielou均勻度指數(shù)變化范圍為 ,平均為0.86,Simpson指數(shù)平均為0.1。

    圖1物種多樣性變化模式Fig.1Change patterns of species diversity

    注:圖1A,1B,1C,1D分別為物種多樣性指數(shù)物種豐富度、香農(nóng)-威納多樣性指數(shù)、Pielou均勻度指數(shù)、Simpson指數(shù)變化模式 Note:Figdspiellateaepfsessto index,and Simpson index

    2.2時(shí)間尺度上群落穩(wěn)定性及效應(yīng)值變化模式

    40年時(shí)間尺度上,群落穩(wěn)定性、種群穩(wěn)定性、統(tǒng)計(jì)平均效應(yīng)隨時(shí)間變化模式為“單峰狀”,異步性為“波動(dòng)狀”,補(bǔ)償效應(yīng)為“波動(dòng)上升狀”(圖3,表1)。

    2.3時(shí)間尺度上物種多樣性與群落穩(wěn)定性及相關(guān)機(jī)制關(guān)系

    物種豐富度和Shannon-Wiener多樣性指數(shù)與群落穩(wěn)定性、種群穩(wěn)定性、異步性、統(tǒng)計(jì)平均效應(yīng)均有顯著的正相關(guān)關(guān)系 ,與補(bǔ)償效應(yīng)無(wú)顯

    圖2群落功能群變化模式Fig.2Change patterns of community functional group

    注:圖2A、2B分別為植物功能群隨年份變化模式和功能群地上生物量平均值。圖2A中AB,PB,PF,PR,SS分別表示一年生或二年生、多年生叢生禾草、多年生雜草、多年生根莖禾草、灌木或半灌木5個(gè)功能群

    Note:FigAandBoteaalaatiopaeoflantfuctioalgrousndteaveaeabovegodassoffuctioaousnFigure2A,ABesptielyeprtiefuclosasalserealcPerennialForbs,Perennial Rhizomatous Grassesand Shrub or Semi-shrubs

    表11982一2021年年份與物種多樣性、群落穩(wěn)定性之間回歸分析Table1 Regression analysis between 1982—2O2l years and species diversity,community stability
    注:AIC表示赤池信息準(zhǔn)則;MOS表示Mitchell-OldsandShawNote:AIC indicates Akaike information criterion;MOS indicates Mitchell-Olds and Shaw

    著相關(guān)關(guān)系;Pielou均勻度指數(shù)與異步性、補(bǔ)償效應(yīng)顯著正相關(guān)( ? lt; 0 . 0 5 ) ,與種群穩(wěn)定性顯著負(fù)相關(guān)( ,與群落穩(wěn)定性、統(tǒng)計(jì)平均效應(yīng)無(wú)顯著相關(guān)關(guān)系;Simpson指數(shù)與群落穩(wěn)定性、異步性、統(tǒng)計(jì)平均效應(yīng)、補(bǔ)償效應(yīng)顯著負(fù)相關(guān)( ? lt; 0 . 0 5 ) ,與種群穩(wěn)定性無(wú)顯著相關(guān)關(guān)系(圖4)??傮w而言,群落中物種多樣性越高,群落物種分布越均勻,整個(gè)群落穩(wěn)定性越高。

    2.4時(shí)間尺度上群落穩(wěn)定性維持機(jī)制

    回歸分析表明,種群穩(wěn)定性、異步性、統(tǒng)計(jì)平均效應(yīng)與群落穩(wěn)定性呈極顯著正相關(guān)關(guān)系 ( P lt; 0.001),群落穩(wěn)定性與補(bǔ)償效應(yīng)顯著正相關(guān)( ( P lt; 0.05)(圖5)。隨機(jī)森林模型結(jié)果表明,種群穩(wěn)定性、統(tǒng)計(jì)平均效應(yīng)和異步性是影響群落穩(wěn)定性的主要機(jī)制(圖6A)。方差分解結(jié)果表明,種群穩(wěn)定性、異步性、統(tǒng)計(jì)平均效應(yīng)、補(bǔ)償效應(yīng)總解釋 9 5 % 的群落穩(wěn)定性變異,解釋度排序?yàn)榉N群穩(wěn)定性 ( 9 4 % ) gt; 異步性 ( 8 6 % ) gt; 統(tǒng)計(jì)平均效應(yīng) ( 6 4 % ) gt; 補(bǔ)償效應(yīng)( 46 % ) ,其中種群穩(wěn)定性、異步性和統(tǒng)計(jì)平均效應(yīng)共同解釋 42 % 的群落穩(wěn)定性變異,異步性和補(bǔ)償效應(yīng)共同解釋 31 % 的群落穩(wěn)定性變異(圖6B)。

    圖3群落穩(wěn)定性及效應(yīng)值變化模式Fig.3 Change patterns of community stability and effect value

    注:圖3A為群落穩(wěn)定性隨時(shí)間變化模式,圖3B,3C,3D,3E分別群落穩(wěn)定性效應(yīng)值種群穩(wěn)定性、異步性、補(bǔ)償效應(yīng)、統(tǒng)計(jì)平均效應(yīng)變化模式 Note:Fig3AilstratesthtoralarationpsofompesatyabiltigBCDnd3Espiveepictal variationpatemityabilituulaiabilityodmpestortatisticalet

    3討論

    云霧山典型草原自1982年實(shí)施封山禁牧措施后,物種多樣性和群落穩(wěn)定性表現(xiàn)為“單峰狀\"變化模式[36-37],產(chǎn)生這種模式的原因可能是:在群落演替前期,封山禁牧措施降低了人為干擾,動(dòng)物對(duì)草地的影響程度急劇減少,自然資源(光照、水分、養(yǎng)分等)對(duì)植物生長(zhǎng)的限制被逐步解除[38],同時(shí)植物養(yǎng)分吸收策略也發(fā)生顯著改變[18],主要以生長(zhǎng)迅速且擴(kuò)散能力強(qiáng)的r對(duì)策植物為主,因此有效促進(jìn)了群落內(nèi)物種的更新和物種多樣性的增加[39],種群穩(wěn)定性和統(tǒng)計(jì)平均效應(yīng)也相應(yīng)增加,且統(tǒng)計(jì)平均效應(yīng)和補(bǔ)償效應(yīng)均通過(guò)增加物種間異步性來(lái)提高群落穩(wěn)定性(圖3,圖5);在群落演替后期,草地發(fā)生灌叢化現(xiàn)象,從早期以多年生叢生禾草和多年生雜草為主向以多年生叢生禾草、多年生雜草、灌木或半灌木(矮腳錦雞兒、白蓮蒿等)為主轉(zhuǎn)變[35,40]。物種組成由r對(duì)策植物逐漸向少數(shù)光競(jìng)爭(zhēng)能力強(qiáng)的k對(duì)策物種轉(zhuǎn)變[41],同時(shí)由于地表凋落物的積累,抑制了下層植物的生長(zhǎng)和發(fā)育[28],種群穩(wěn)定性和種間異步性也逐漸降低,從而降低了群落穩(wěn)定性。總之,禁牧措施對(duì)早期演替物種更新以及功能維持起到積極作用[42],但長(zhǎng)時(shí)間禁牧并不利于物種多樣性和群落穩(wěn)定性的維持[43]

    圖4物種多樣性與群落穩(wěn)定性相關(guān)性矩陣

    自20世紀(jì)70年代以來(lái),物種多樣性-群落穩(wěn)定性的正相關(guān)關(guān)系被多數(shù)生態(tài)學(xué)家所接受[44-45],即一個(gè)群落內(nèi)物種豐富度越大,種群分布越均勻,則整個(gè)群落的穩(wěn)定性越高[18]。MacArthur等[46-48]提出的多樣性與穩(wěn)定性關(guān)系假說(shuō)認(rèn)為,一個(gè)群落內(nèi)物種組成越豐富且種群大小越穩(wěn)定,那么整個(gè)群落的穩(wěn)定性就越高。本研究中物種豐富度和Shannon-wiener多樣性指數(shù)與群落穩(wěn)定性呈極顯著正相關(guān)關(guān)系為該假說(shuō)提供了有力支撐(圖4)。Houlahan等[21]分析了北美洲91個(gè)群落數(shù)據(jù)也發(fā)現(xiàn)物種豐富度與種群穩(wěn)定性、群落穩(wěn)定性呈正相關(guān)。Downing等49]研究了水生食物網(wǎng)也發(fā)現(xiàn)物種豐富群落的種群穩(wěn)定性、群落穩(wěn)定性也較高。一項(xiàng)10年尺度的控制草地物種多樣性實(shí)驗(yàn)結(jié)果支持競(jìng)爭(zhēng)理論的預(yù)測(cè),即更大的物種多樣性導(dǎo)致更高的群落穩(wěn)定性和較低的種群穩(wěn)定性[10]。研究結(jié)果產(chǎn)生差異的主要原因可能是研究對(duì)象不同,本研究對(duì)象為野外封育狀態(tài)下的自然植物群落,保持了原有群落的物種豐富度和結(jié)構(gòu),而Tilman等[10對(duì)直接控制物種豐富度的群落進(jìn)行研究。此外,本研究結(jié)果中,Pielou均勻度指數(shù)與補(bǔ)償效應(yīng)顯著正相關(guān)、Simpson指數(shù)與補(bǔ)償效應(yīng)顯著負(fù)相關(guān)以及補(bǔ)償效應(yīng)與群落穩(wěn)定性顯著正相關(guān)也為全面理解物種共存和生態(tài)位理論提供了佐證[50-51]

    圖5群落穩(wěn)定性和效應(yīng)值回歸分析Fig.5Regression analysis of community stability and effect value

    Note:FigABdspectielntteiaspulaiblitcatistalea pensatoryeffect,and communitystability

    4,圖5)。自然群落中不同物種之間具有生態(tài)位差異,種群分布越均勻,補(bǔ)償效應(yīng)發(fā)揮的作用就越大,物種對(duì)資源的分化利用就越強(qiáng)[52],從而促進(jìn)了物種共存和群落穩(wěn)定。

    圖6"群落穩(wěn)定性和效應(yīng)值的隨機(jī)森林模型(A)和方案分解{B}結(jié)果

    物種多樣性主要通過(guò)增強(qiáng)種間異步性來(lái)維持群落穩(wěn)定[53-54],種間異步性可進(jìn)一步分解為補(bǔ)償效應(yīng)和統(tǒng)計(jì)平均效應(yīng)[23]。本研究結(jié)果表明,統(tǒng)計(jì)平均效應(yīng)和補(bǔ)償效應(yīng)均對(duì)群落穩(wěn)定性起關(guān)鍵作用,并且統(tǒng)計(jì)平均效應(yīng)比補(bǔ)償效應(yīng)更重要,這是解釋物種多樣性-群落穩(wěn)定性正相關(guān)關(guān)系的關(guān)鍵機(jī)制(圖6)。本研究也發(fā)現(xiàn)物種豐富度主要通過(guò)增加群落生物量和降低種群隨機(jī)性來(lái)穩(wěn)定群落,種群隨機(jī)性比環(huán)境隨機(jī)性更關(guān)鍵,統(tǒng)計(jì)平均效應(yīng)發(fā)揮了關(guān)鍵作用。Gonzalez等[56]同樣認(rèn)為物種或種群之間的隨機(jī)組合在物種多樣性-群落穩(wěn)定性正相關(guān)關(guān)系中發(fā)揮重要作用[55],且統(tǒng)計(jì)平均效應(yīng)可能更重要。Hallett等[15]利用9個(gè)草地物種多樣性實(shí)驗(yàn)研究了降水變異對(duì)群落穩(wěn)定性的影響,認(rèn)為群落穩(wěn)定性機(jī)制的重要性會(huì)隨環(huán)境梯度變化而變化,補(bǔ)償效應(yīng)對(duì)群落穩(wěn)定性的貢獻(xiàn)隨降水變異性增加而增大。隨著研究空間尺度的增加,環(huán)境隨機(jī)性也隨之增強(qiáng),這可能導(dǎo)致種群隨機(jī)性對(duì)群落穩(wěn)定性的影響相對(duì)降低[57]。相應(yīng)地,補(bǔ)償效應(yīng)在維持群落穩(wěn)定性中的作用可能會(huì)變得更加顯著。因此,未來(lái)的研究應(yīng)更多關(guān)注不同環(huán)境梯度下各種群落穩(wěn)定性機(jī)制的相對(duì)重要性,這對(duì)于預(yù)測(cè)在全球氣候變化背景下生物多樣性的減少及其對(duì)生態(tài)系統(tǒng)穩(wěn)定性的潛在影響至關(guān)重要。

    4結(jié)論

    云霧山40年的草地生態(tài)系統(tǒng)演替表明,物種多樣性與群落穩(wěn)定性以正相關(guān)關(guān)系為主,群落物種越豐富且分布越均勻,則群落穩(wěn)定性越高。在較小的空間尺度上,種群穩(wěn)定性和異步性是維持群落穩(wěn)定的主要機(jī)制,且統(tǒng)計(jì)平均效應(yīng)比補(bǔ)償效應(yīng)更重要。本研究為草地生態(tài)系統(tǒng)生物多樣性保護(hù)和可持續(xù)管理提供了重要理論依據(jù)。

    參考文獻(xiàn)

    [1]ISBELLFI,POLLEYHW,WILSEYBJ.Biodiversity,Pro ductivity and the Temporal Stability of Productivity:Patterns and Processes[J].EcologyLetters,2009,12(5):443-451

    [2] MCCANNKS.Thediversity-stability debate[J].Nature, 2000,405(6783):228-233

    [3] 李奇,朱建華,肖文發(fā).生物多樣性與生態(tài)系統(tǒng)服務(wù)一關(guān) 系、權(quán)衡與管理[J].生態(tài)學(xué)報(bào),2019,39(8):2655-2666

    [4] WHITEHJ,CAPLATP,EMMERSONMC,etal.Predictingfuture stability of ecosystem functioning under climate change[J].Agriculture,Ecosystemsamp;Environment,2021, 320:107600

    [5] 李周園,葉小洲,王少鵬,生態(tài)系統(tǒng)穩(wěn)定性及其與生物多樣性 的關(guān)系[J].植物生態(tài)學(xué)報(bào),2021,45(10):1127-1139

    [6] WALKERLR,WARDLED A.Plant succession asaninte gratorof contrasting ecological time scales[J]. Trends in Ecologyamp;Evolution,2014,29(9):504-510

    [7] CONNELLJH,SLATYERRO.Mechanismsof succession innatural communitiesand their rolein community stability and organization[J].TheAmericanNaturalist,1977,111(982): 1119-1144

    [8] ELTON C S. The ecology of invasions by animals and plants [M].London:Methuen,1958:143-153

    [9] MAYR M.Will a large complex system be stable?[J]. Nature,1972,238(5364):413-414

    [10]TILMAN D,REICHP B,KNOPS J M H. Biodiversity and ecosystem stability in a decade-long grassland experiment[J]. Nature,2006,441(7093):629-632

    [11]RIVA F,GRACO-ROZA C,DASKALOVA G N,et al. Toward a cohesive understanding of ecological complexity[J]. Science Advances,2023,9(25):eabq4207

    [12]LI W J,ZHOU X,XIANG Z Q,et al. Biomass temporal stabilityincreases at two spatial scales during secondary succession [J].JournalofEcology,2023,111(7):1575-1586

    [13]沈艷,馬紅彬,趙菲,等.荒漠草原土壤養(yǎng)分和植物群落穩(wěn)定 性對(duì)不同管理方式的響應(yīng)[J].草地學(xué)報(bào),2015,23(2): 264-270

    [14]劉晶,楊雪,張博,等.長(zhǎng)期施肥對(duì)黃土高原典型草原群落穩(wěn) 定性的影響及機(jī)制研究[J].西北植物學(xué)報(bào),2021,41(2): 310-316

    [15]HALLETTLM,HSUJS,CLELANDEE,etal.Biotic mechanisms of community stability shift along a precipitation gradient[J]. Ecology,2014,95(6):1693-1700

    [16]LIU J,YANG X,GHANIZADEH H,et al. Long-term enclosure can benefit grassland community stability on the Loess Plateau of China[J].Sustainability,2021,13(1):213

    [17]TILMAN D.The ecological consequences of changes in biodiversity:a search for general principles[J].Ecology,1999,80 (5):1455

    [18] TILMAN D,ISBELL F,COWLES JM. Biodiversity and ecosystem functioning[J]. Annual Review of Ecology,Evolution, and Systematics,2014,45:471-493

    [19] TREDENNICK A T,DE MAZANCOURT C,LOREAU M, et al.Environmental responses,not species interactions,deter mine synchrony of dominant species in semiarid grasslands[J]. Ecology,2017,98(4):971-981

    [20]THIBAUT L M,CONNOLLY S R.Understanding diversitystability relationships:towards a unified model of portfolio effects[J].Ecology Letters,2013,16(2):140-150

    [21]HOULAHAN JE,CURRIE D J,COTTENIE K,et al. Negative relationships between species richness and temporal variability arecommon but weak in natural systems[J]. Ecology, 2018,99(11):2592-2604

    [22] CRAVEN D,EISENHAUER N,PEARSE W D,et al. Multiplefacetsofbiodiversity drive the diversity-stabilityrelationship[J].Nature Ecologyamp; Evolution,2018,2(10):1579-1587

    [23]ZHAOL,WANGSP,SHENRH,etal.Biodiversitystabilizes plant communities through statistical-averaging effects rather than compensatory dynamics[J].Nature Communications,2022,13(1):7804

    [24]BAI Y F,HAN XG,WU JG,et al. Ecosystem stability and compensatory effects in the Inner Mongolia grassland[J]. Nature,2004,431(7005):181-184

    [25] CARDINALE B J,WRIGHT J P,CADOTTE M W,et al. Impacts of plant diversity on biomass production increase through time because of species complementarity[J].Proceedings of the National Academy of Sciences of the United States ofAmerica,2007,104(46):18123-18128

    [26]DOAKDF,BIGGERD,HARDINGEK,etal.The statistical inevitability of stability-diversity relationships in community ecology[J]. The American Naturalist,1998,151(3):264-276

    [27]鄒厚遠(yuǎn),關(guān)秀琦,張信,等.云霧山草原自然保護(hù)區(qū)的管理途 徑探討[J].草業(yè)科學(xué),1997,14(1):3-4

    [28]趙凌平,白欣,譚世圖,等.不同年限封育對(duì)黃土高原典型草 原地上植被的影響[J].草業(yè)科學(xué),2018,35(1):27-35

    [29]程積民,井趙斌,金晶煒,等.黃土高原半干旱區(qū)退化草地恢 復(fù)與利用過(guò)程研究[J].中國(guó)科學(xué):生命科學(xué),2014,44(3): 267-279

    [30]史曉曉,程積民,于飛,等.云霧山天然草地30年恢復(fù)演替過(guò) 程中優(yōu)勢(shì)草種生態(tài)位動(dòng)態(tài)[J].草地學(xué)報(bào),2014,22(4): 677-684

    [31]張譯文,常小峰,張玉薇.黃土高原云霧山長(zhǎng)期封育草地灌叢 化空間分布格局[J].草地學(xué)報(bào),2023,31(12):3802-3808

    [32]楊雪.時(shí)間尺度上黃土高原典型草原物種多樣性和地上生產(chǎn) 力的關(guān)系研究[D].楊凌:西北農(nóng)林科技大學(xué),2022:8-11

    [33]井趙斌.云霧山退化草地生態(tài)系統(tǒng)恢復(fù)與調(diào)控機(jī)理[D].楊 凌:西北農(nóng)林科技大學(xué),2014:14

    [34]楊雪,晏昕輝,李美慧,等.時(shí)間尺度上草地物種多樣性和地 上生產(chǎn)力的關(guān)系[J].草地學(xué)報(bào),2022,30(2):259-268

    [35]李美慧,李玉華,晏昕輝,等.半灌木擴(kuò)張驅(qū)動(dòng)的草地植物多 樣性與地上生產(chǎn)力特征及其關(guān)系研究[J].草業(yè)學(xué)報(bào),2023,32 (5):27-39

    [36] BONET A,PAUSAS JG. Species richness and cover along a 60-yearchronosequenceinold-fields of southeastern Spain[J]. Plant Ecology,2004,174(2):257-270

    [37]LI W J,LIJH,LIU S S,et al. Magnitude of species diversity effect on aboveground plant biomass increases through successional time of abandoned farmlands on the eastern Tibetan Plateau of China[J]. Land Degradation amp; Development,2017,28 (1):370-378

    [38]LOREAU M. Ecosystem development explained by competition within and between material cycles[J]. Proceedings of the Royal Society of London Series B: Biological Sciences,1998, 265(1390):33-38

    [39]趙凌平,譚世圖,白欣,等.封育年限對(duì)云霧山典型草原植物 繁殖與植被更新的影響[J].草業(yè)學(xué)報(bào),2017,26(10):1-9

    [40]拓行行,黃建弟,白玥,等.半灌木擴(kuò)張草原生長(zhǎng)季內(nèi)土壤氮 礦化特征及其影響因素[J].草地學(xué)報(bào),2023,31(9):2730- 2739

    [41]PICKETT S T A. Population patterns through twenty years of oldfield succession[J].Vegetatio,1982,49(1):45-59

    [42]單貴蓮,徐柱,寧發(fā),等.圍封年限對(duì)典型草原群落結(jié)構(gòu)及物 種多樣性的影響[J].草業(yè)學(xué)報(bào),2008,17(6):1-8

    [43]王黎黎,張克斌,程中秋,等.圍欄封育對(duì)半干旱區(qū)植物群落 穩(wěn)定性的影響[J].甘肅農(nóng)業(yè)大學(xué)學(xué)報(bào),2011,46(5):86-90

    [44]SHE Y D,LI X L,LI C Y,et al. Relationship between species diversityand communitystabilityindegraded alpinemeadows during bare patch succession[J]. Plants,2023,12(20):3582

    [45]VAN PLASF.Biodiversity and ecosystem functioning in naturallyassembled communities[J].Biological Reviews of the Cambridge Philosophical Society,2019,94(4):1220-1245

    [46]MACARTHUR R.Fluctuations of animal populations and a measureof community stability[J].Ecology,1955,36(3): 533-536

    [47]黃建輝,韓興國(guó).生物多樣性和生態(tài)系統(tǒng)穩(wěn)定性[J].生物多 樣性,1995,3(1):31-37

    [48]李永政.多樣性—穩(wěn)定性關(guān)系與環(huán)境保護(hù)[D].鄭州:河南大 學(xué),2012:3-18

    [49]DOWNING AL,BROWNBL,LEIBOLDMA.Multiple diversity-stability mechanisms enhance population and communitystabilityin aquatic food webs[J].Ecology,20l4,95(1): 173-184

    [50]李偉.施肥導(dǎo)致高寒草甸物種多樣性喪失機(jī)制研究[D].蘭 州:蘭州大學(xué),2011:1-3

    [51]儲(chǔ)誠(chéng)進(jìn),王酉石,劉宇,等.物種共存理論研究進(jìn)展[J].生物 多樣性,2017,25(4):345-354

    [52]WANGCY,WEI M,WUBD,etal.Alpine grassland degradationreduced plant species diversity and stability of plant communities in theNorthern TibetPlateau[J].Acta Oecologica, 2019,98:25-29

    [53]SCHNABELF,LIUXJ,KUNZM,etal.Species richness stabilizes productivity via asynchrony and drought-tolerance diversityina large-scale tree biodiversityexperiment[J].Science Advances,2021,7(51):eabk1643

    [54]LOREAUM,DEMAZANCOURTC.Species synchrony and itsdrivers:neutral and nonneutral community dynamicsin fluctuatingenvironments[J].The American Naturalist,2Oo8,172 (2):E48-E66

    [55]DEMAZANCOURTC,ISBELLF,LAROCQUEA,etal. Predicting ecosystem stability from community composition and biodiversity[J].EcologyLetters,2013,16(5):617-625

    [56]GONZALEZA,GERMAINRM,SRIVASTAVADS,etal. Scaling-up biodiversity-ecosystem functioning research [J]. EcologyLetters,2020,23(4):757-776

    [57]CHALCRAFTDR. Changes in ecological stability across realisticbiodiversity gradientsdepend on spatial scale[J].Global EcologyandBiogeography,2013,22(1):19-28

    (責(zé)任編輯彭露茜)

    猜你喜歡
    物種效應(yīng)
    物種大偵探
    物種大偵探
    吃光入侵物種真的是解決之道嗎?
    鈾對(duì)大型溞的急性毒性效應(yīng)
    懶馬效應(yīng)
    場(chǎng)景效應(yīng)
    回首2018,這些新物種值得關(guān)注
    電咖再造新物種
    應(yīng)變效應(yīng)及其應(yīng)用
    瘋狂的外來(lái)入侵物種
    嫩草影院精品99| 夜夜夜夜夜久久久久| 99精品在免费线老司机午夜| 久久久久性生活片| 成年人黄色毛片网站| 成人特级av手机在线观看| 欧美成人一区二区免费高清观看| 久久伊人香网站| 国产熟女xx| 国产免费一级a男人的天堂| 色噜噜av男人的天堂激情| 美女黄网站色视频| 午夜日韩欧美国产| 毛片女人毛片| 国产亚洲精品久久久com| 91久久精品国产一区二区成人| 国产精品av视频在线免费观看| 少妇丰满av| 亚洲国产精品999在线| 欧美一区二区精品小视频在线| 国产乱人视频| 久久欧美精品欧美久久欧美| 免费电影在线观看免费观看| 日本 欧美在线| 国产v大片淫在线免费观看| 99久久99久久久精品蜜桃| 欧美+亚洲+日韩+国产| 日本精品一区二区三区蜜桃| 日韩欧美精品免费久久 | 99热这里只有是精品50| 日韩中文字幕欧美一区二区| 日韩精品青青久久久久久| 欧美黄色片欧美黄色片| 真人一进一出gif抽搐免费| 国产在视频线在精品| 老熟妇乱子伦视频在线观看| 国产私拍福利视频在线观看| 成人无遮挡网站| 亚洲av不卡在线观看| 成人永久免费在线观看视频| 美女大奶头视频| 99国产综合亚洲精品| 无人区码免费观看不卡| 国产伦在线观看视频一区| 男人舔奶头视频| 床上黄色一级片| 午夜精品在线福利| 国产欧美日韩精品一区二区| 神马国产精品三级电影在线观看| 国产一区二区激情短视频| 久久久色成人| a级一级毛片免费在线观看| 免费在线观看影片大全网站| 久久久久久久午夜电影| 日本黄色片子视频| 日韩精品中文字幕看吧| 国产av在哪里看| 国产精品久久久久久亚洲av鲁大| 亚洲精品影视一区二区三区av| 亚洲自拍偷在线| 深夜a级毛片| 国产伦精品一区二区三区视频9| 少妇的逼好多水| 性色av乱码一区二区三区2| 久久久久国内视频| 窝窝影院91人妻| 国产免费一级a男人的天堂| 一二三四社区在线视频社区8| 禁无遮挡网站| 中文字幕久久专区| 成人毛片a级毛片在线播放| 国产精品三级大全| 美女高潮的动态| 好看av亚洲va欧美ⅴa在| 国产av在哪里看| 日韩欧美免费精品| 欧洲精品卡2卡3卡4卡5卡区| 在线观看免费视频日本深夜| 哪里可以看免费的av片| 三级毛片av免费| 国产激情偷乱视频一区二区| 精华霜和精华液先用哪个| 亚洲av第一区精品v没综合| 国产精品亚洲av一区麻豆| 熟妇人妻久久中文字幕3abv| 九九久久精品国产亚洲av麻豆| 性插视频无遮挡在线免费观看| 亚洲美女搞黄在线观看 | 精品一区二区三区av网在线观看| 国产一区二区三区视频了| 日本精品一区二区三区蜜桃| 一进一出抽搐动态| 免费看光身美女| 亚洲片人在线观看| 男人和女人高潮做爰伦理| 免费观看精品视频网站| 深夜a级毛片| 免费人成视频x8x8入口观看| av在线老鸭窝| 神马国产精品三级电影在线观看| 91久久精品国产一区二区成人| 尤物成人国产欧美一区二区三区| 成人av一区二区三区在线看| 日韩欧美免费精品| 校园春色视频在线观看| 精品一区二区三区人妻视频| 国产精品女同一区二区软件 | 欧美乱色亚洲激情| 深夜a级毛片| 老鸭窝网址在线观看| 一本精品99久久精品77| 亚洲精品亚洲一区二区| 搡老妇女老女人老熟妇| 中文字幕人妻熟人妻熟丝袜美| 在线a可以看的网站| 精品乱码久久久久久99久播| 国产人妻一区二区三区在| 免费在线观看日本一区| 亚洲真实伦在线观看| 免费看美女性在线毛片视频| 欧美国产日韩亚洲一区| 又粗又爽又猛毛片免费看| 欧美性猛交╳xxx乱大交人| 麻豆av噜噜一区二区三区| 日本黄色片子视频| 亚洲色图av天堂| 国产欧美日韩一区二区三| 麻豆国产av国片精品| 丁香六月欧美| 中文字幕人成人乱码亚洲影| 麻豆国产av国片精品| 国产精品人妻久久久久久| 亚洲欧美日韩无卡精品| 性色avwww在线观看| 国内精品久久久久精免费| 大型黄色视频在线免费观看| 99国产精品一区二区蜜桃av| 国产在线精品亚洲第一网站| 日韩欧美精品v在线| 亚洲欧美日韩卡通动漫| 夜夜躁狠狠躁天天躁| 日本免费一区二区三区高清不卡| 啪啪无遮挡十八禁网站| 啦啦啦韩国在线观看视频| 成人国产一区最新在线观看| 亚洲精华国产精华精| 怎么达到女性高潮| 久久国产乱子免费精品| 欧美成狂野欧美在线观看| 国产伦精品一区二区三区四那| 久久久久亚洲av毛片大全| 搡老妇女老女人老熟妇| 99久久精品国产亚洲精品| 亚洲久久久久久中文字幕| 哪里可以看免费的av片| 国产黄色小视频在线观看| 简卡轻食公司| 真实男女啪啪啪动态图| 亚洲专区国产一区二区| 婷婷精品国产亚洲av在线| avwww免费| 性色av乱码一区二区三区2| 美女高潮的动态| 观看免费一级毛片| 亚洲欧美日韩东京热| 国产毛片a区久久久久| 欧美+亚洲+日韩+国产| 色哟哟·www| 麻豆成人午夜福利视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 十八禁国产超污无遮挡网站| 此物有八面人人有两片| 啪啪无遮挡十八禁网站| 日韩欧美在线乱码| 丰满的人妻完整版| 久久久国产成人精品二区| 国产黄片美女视频| 亚洲无线在线观看| 国产aⅴ精品一区二区三区波| 1000部很黄的大片| 免费大片18禁| 九色国产91popny在线| 91久久精品电影网| 欧美成人a在线观看| 十八禁国产超污无遮挡网站| 亚洲国产精品999在线| 亚洲五月婷婷丁香| 日本熟妇午夜| 变态另类丝袜制服| 天天躁日日操中文字幕| 十八禁人妻一区二区| av在线观看视频网站免费| 久久99热这里只有精品18| 亚洲最大成人av| 麻豆成人午夜福利视频| 免费观看精品视频网站| 国产一级毛片七仙女欲春2| 深夜精品福利| 日韩欧美精品免费久久 | 少妇的逼好多水| 如何舔出高潮| 国产免费一级a男人的天堂| 一个人免费在线观看电影| 国产久久久一区二区三区| 免费在线观看影片大全网站| 国产高清有码在线观看视频| 怎么达到女性高潮| 国产欧美日韩一区二区三| 国产日本99.免费观看| 国产精品98久久久久久宅男小说| 有码 亚洲区| 免费在线观看成人毛片| 91麻豆av在线| 欧美成人性av电影在线观看| 欧美xxxx性猛交bbbb| 精品午夜福利视频在线观看一区| 亚洲国产日韩欧美精品在线观看| 亚洲av成人不卡在线观看播放网| 久久精品综合一区二区三区| 精品熟女少妇八av免费久了| 看黄色毛片网站| 国产三级中文精品| 丰满乱子伦码专区| 亚洲人成电影免费在线| 在线观看美女被高潮喷水网站 | 老司机午夜十八禁免费视频| 国产欧美日韩精品一区二区| 国产三级中文精品| 成人午夜高清在线视频| 中文字幕精品亚洲无线码一区| 国内精品久久久久精免费| 免费在线观看影片大全网站| a级毛片免费高清观看在线播放| 成人av在线播放网站| 日本一本二区三区精品| 神马国产精品三级电影在线观看| 91久久精品国产一区二区成人| 人妻久久中文字幕网| 国内精品一区二区在线观看| 人妻制服诱惑在线中文字幕| av女优亚洲男人天堂| 午夜激情欧美在线| 国产高清三级在线| 一级黄色大片毛片| 午夜福利欧美成人| 久久久久国内视频| 国产精品永久免费网站| 精品国产亚洲在线| 国产亚洲欧美在线一区二区| 欧美3d第一页| 黄色日韩在线| 国产精品免费一区二区三区在线| 日韩 亚洲 欧美在线| 看黄色毛片网站| 久久精品人妻少妇| www.色视频.com| 成人无遮挡网站| 波多野结衣巨乳人妻| 日韩欧美精品v在线| 嫩草影院入口| 亚洲欧美清纯卡通| ponron亚洲| 999久久久精品免费观看国产| 亚洲av电影不卡..在线观看| 成年女人毛片免费观看观看9| 在线免费观看不下载黄p国产 | 亚洲第一区二区三区不卡| 麻豆久久精品国产亚洲av| 12—13女人毛片做爰片一| 国产免费一级a男人的天堂| 一区二区三区高清视频在线| 久久精品人妻少妇| 久久伊人香网站| 一本一本综合久久| 国产精华一区二区三区| 嫩草影院入口| 午夜影院日韩av| 直男gayav资源| 热99re8久久精品国产| АⅤ资源中文在线天堂| 我的女老师完整版在线观看| av天堂在线播放| 91在线观看av| 亚洲精品成人久久久久久| 亚洲专区国产一区二区| 午夜免费成人在线视频| 国产国拍精品亚洲av在线观看| 51午夜福利影视在线观看| 日韩亚洲欧美综合| 在线播放国产精品三级| 久久久久久久久久成人| 国产乱人视频| 精华霜和精华液先用哪个| 欧美区成人在线视频| 国产午夜福利久久久久久| 老熟妇仑乱视频hdxx| 精品人妻偷拍中文字幕| 99久久九九国产精品国产免费| 亚洲第一欧美日韩一区二区三区| 男人的好看免费观看在线视频| 最近视频中文字幕2019在线8| 日韩欧美三级三区| 搡老岳熟女国产| 波多野结衣高清无吗| 亚洲,欧美,日韩| 欧美最新免费一区二区三区 | 亚洲中文字幕日韩| 亚洲av二区三区四区| 一个人看的www免费观看视频| 婷婷丁香在线五月| 国产精品av视频在线免费观看| 中亚洲国语对白在线视频| 我要搜黄色片| 听说在线观看完整版免费高清| 九九在线视频观看精品| 成年版毛片免费区| 99久久久亚洲精品蜜臀av| 亚洲成人免费电影在线观看| 国产一区二区在线观看日韩| 免费看美女性在线毛片视频| 成人鲁丝片一二三区免费| 亚洲第一欧美日韩一区二区三区| 亚洲成人久久性| 国产亚洲av嫩草精品影院| 99国产精品一区二区三区| 午夜福利免费观看在线| 免费人成在线观看视频色| 又粗又爽又猛毛片免费看| a级毛片a级免费在线| 十八禁人妻一区二区| 老司机午夜十八禁免费视频| 国产又黄又爽又无遮挡在线| 又黄又爽又免费观看的视频| 国产一级毛片七仙女欲春2| 亚洲中文字幕日韩| 日韩精品青青久久久久久| 国产精品免费一区二区三区在线| 淫秽高清视频在线观看| 亚洲最大成人手机在线| 天堂影院成人在线观看| 欧美日韩瑟瑟在线播放| 无遮挡黄片免费观看| 偷拍熟女少妇极品色| 亚洲一区高清亚洲精品| 欧美最黄视频在线播放免费| 日韩国内少妇激情av| 国内少妇人妻偷人精品xxx网站| 国产高清视频在线播放一区| 婷婷色综合大香蕉| 免费人成视频x8x8入口观看| 无人区码免费观看不卡| 久久久久精品国产欧美久久久| 三级男女做爰猛烈吃奶摸视频| 国产精品一及| 俺也久久电影网| 国产精品野战在线观看| 精品人妻偷拍中文字幕| 欧美一区二区国产精品久久精品| 在线免费观看不下载黄p国产 | 91麻豆av在线| 婷婷亚洲欧美| 久久欧美精品欧美久久欧美| 精品久久国产蜜桃| 1024手机看黄色片| 成人欧美大片| 国产精品影院久久| 国产精品一区二区三区四区久久| 欧美成人性av电影在线观看| 欧美另类亚洲清纯唯美| 日韩精品中文字幕看吧| 国产精品一区二区三区四区久久| 美女 人体艺术 gogo| 99国产精品一区二区三区| 五月玫瑰六月丁香| 精品日产1卡2卡| 日本免费一区二区三区高清不卡| 一个人观看的视频www高清免费观看| 亚洲欧美日韩高清在线视频| 2021天堂中文幕一二区在线观| 嫩草影院入口| 欧洲精品卡2卡3卡4卡5卡区| 久久国产精品影院| 日韩精品青青久久久久久| 亚洲av一区综合| 国产成人啪精品午夜网站| 天天一区二区日本电影三级| 99国产极品粉嫩在线观看| 亚洲狠狠婷婷综合久久图片| 成人鲁丝片一二三区免费| netflix在线观看网站| 最好的美女福利视频网| 啪啪无遮挡十八禁网站| 欧美一级a爱片免费观看看| 午夜福利18| 最近最新免费中文字幕在线| 免费一级毛片在线播放高清视频| 国产一区二区在线观看日韩| 国产精品亚洲av一区麻豆| 亚洲av不卡在线观看| 日本黄大片高清| 国产一区二区三区在线臀色熟女| 精品熟女少妇八av免费久了| 亚洲av中文字字幕乱码综合| 日韩中字成人| 18禁裸乳无遮挡免费网站照片| 亚洲不卡免费看| 99久久久亚洲精品蜜臀av| 精品久久久久久久久av| 亚洲五月天丁香| 久久久精品大字幕| 国产午夜精品久久久久久一区二区三区 | 网址你懂的国产日韩在线| 国产一区二区在线观看日韩| 午夜亚洲福利在线播放| 精品久久久久久成人av| 成人毛片a级毛片在线播放| 18禁在线播放成人免费| 国产成人a区在线观看| 日本与韩国留学比较| 亚洲精华国产精华精| 99精品在免费线老司机午夜| 又紧又爽又黄一区二区| 两个人的视频大全免费| 亚洲,欧美精品.| 简卡轻食公司| 99精品在免费线老司机午夜| 国产一级毛片七仙女欲春2| 乱人视频在线观看| 亚洲av不卡在线观看| 国内少妇人妻偷人精品xxx网站| 18禁在线播放成人免费| 精品午夜福利视频在线观看一区| 长腿黑丝高跟| 久久人人精品亚洲av| 毛片一级片免费看久久久久 | 中国美女看黄片| 在线观看舔阴道视频| 国产黄a三级三级三级人| 中文在线观看免费www的网站| 精品久久久久久久久亚洲 | 日本三级黄在线观看| 97人妻精品一区二区三区麻豆| 午夜精品久久久久久毛片777| 禁无遮挡网站| 伊人久久精品亚洲午夜| 亚洲av免费在线观看| 欧美另类亚洲清纯唯美| 色在线成人网| 别揉我奶头~嗯~啊~动态视频| 在线a可以看的网站| 国产一区二区在线观看日韩| 国产精品一及| 狂野欧美白嫩少妇大欣赏| 国产精品1区2区在线观看.| 国产精品三级大全| 麻豆av噜噜一区二区三区| av天堂中文字幕网| 久久人人爽人人爽人人片va | 在现免费观看毛片| 赤兔流量卡办理| 欧美+亚洲+日韩+国产| 久久天躁狠狠躁夜夜2o2o| 国产高清三级在线| 日本与韩国留学比较| 精品人妻熟女av久视频| 黄色一级大片看看| 露出奶头的视频| 精品国内亚洲2022精品成人| 一区福利在线观看| x7x7x7水蜜桃| 一级黄色大片毛片| 亚洲av一区综合| 国产精品野战在线观看| 波多野结衣高清作品| 亚洲第一电影网av| 五月伊人婷婷丁香| 一级黄片播放器| 精品一区二区三区视频在线观看免费| 97超视频在线观看视频| 欧美色欧美亚洲另类二区| 国产精品人妻久久久久久| 国产高清激情床上av| 女人被狂操c到高潮| 国产男靠女视频免费网站| 日本免费a在线| 久久国产精品人妻蜜桃| 久久精品人妻少妇| 久久午夜福利片| 国产精品不卡视频一区二区 | 中文字幕精品亚洲无线码一区| 午夜福利高清视频| 可以在线观看的亚洲视频| 亚洲专区中文字幕在线| 能在线免费观看的黄片| 国内精品美女久久久久久| 在现免费观看毛片| 不卡一级毛片| 欧美乱妇无乱码| 亚洲一区二区三区色噜噜| 亚洲最大成人手机在线| 亚洲人成网站在线播放欧美日韩| 露出奶头的视频| 国产 一区 欧美 日韩| 国产午夜精品久久久久久一区二区三区 | 亚洲欧美清纯卡通| av福利片在线观看| 最近最新免费中文字幕在线| 国语自产精品视频在线第100页| 99国产精品一区二区蜜桃av| 亚洲,欧美,日韩| 欧美性猛交╳xxx乱大交人| a级一级毛片免费在线观看| 亚洲av电影在线进入| 欧美高清性xxxxhd video| 在线观看av片永久免费下载| 少妇裸体淫交视频免费看高清| 免费黄网站久久成人精品 | 最近在线观看免费完整版| 欧美一区二区精品小视频在线| 色精品久久人妻99蜜桃| 亚洲av二区三区四区| 国产黄色小视频在线观看| 欧美日韩国产亚洲二区| 少妇裸体淫交视频免费看高清| 日韩国内少妇激情av| 一区二区三区四区激情视频 | 少妇丰满av| 特大巨黑吊av在线直播| 一二三四社区在线视频社区8| 伊人久久精品亚洲午夜| 精品一区二区三区视频在线| 国产不卡一卡二| 国产 一区 欧美 日韩| 偷拍熟女少妇极品色| 亚洲中文字幕日韩| 久久久国产成人精品二区| 国产91精品成人一区二区三区| 小蜜桃在线观看免费完整版高清| 欧美xxxx黑人xx丫x性爽| 国产免费一级a男人的天堂| 简卡轻食公司| 精品人妻1区二区| 国产视频内射| 精品不卡国产一区二区三区| 欧美成人性av电影在线观看| 亚洲av中文字字幕乱码综合| 丰满人妻一区二区三区视频av| 国产aⅴ精品一区二区三区波| 深夜精品福利| 一本精品99久久精品77| 亚洲精品日韩av片在线观看| 成人性生交大片免费视频hd| 精品一区二区三区人妻视频| 美女cb高潮喷水在线观看| 日本成人三级电影网站| 亚洲aⅴ乱码一区二区在线播放| 熟妇人妻久久中文字幕3abv| 丝袜美腿在线中文| 欧美区成人在线视频| 最近最新免费中文字幕在线| 欧美+日韩+精品| 美女xxoo啪啪120秒动态图 | 亚洲国产欧洲综合997久久,| 久久久久精品国产欧美久久久| 我的女老师完整版在线观看| 怎么达到女性高潮| 床上黄色一级片| 日日摸夜夜添夜夜添小说| 亚洲av成人av| 午夜久久久久精精品| 欧美bdsm另类| 久久99热6这里只有精品| 最好的美女福利视频网| bbb黄色大片| 日本熟妇午夜| 琪琪午夜伦伦电影理论片6080| 搞女人的毛片| 18美女黄网站色大片免费观看| 精品福利观看| 一级毛片久久久久久久久女| 久久久久久久久久成人| 51午夜福利影视在线观看| .国产精品久久| 男插女下体视频免费在线播放| 国产极品精品免费视频能看的| 热99re8久久精品国产| 亚洲欧美日韩无卡精品| 直男gayav资源| 变态另类丝袜制服| 观看免费一级毛片| 亚洲激情在线av| 国产激情偷乱视频一区二区| 国产精品女同一区二区软件 | 国产成人a区在线观看| 男女做爰动态图高潮gif福利片| 一个人免费在线观看的高清视频| 99久久99久久久精品蜜桃| 日本免费a在线| 国产白丝娇喘喷水9色精品| 精品久久久久久久末码| 国产免费一级a男人的天堂| 午夜福利在线观看吧| 午夜精品久久久久久毛片777| 欧美又色又爽又黄视频| 人妻久久中文字幕网| 一进一出好大好爽视频| 啦啦啦观看免费观看视频高清| 亚洲无线在线观看| 伦理电影大哥的女人| 在线观看午夜福利视频| 一进一出抽搐gif免费好疼| 欧美乱妇无乱码| 国内精品美女久久久久久| 亚洲第一电影网av| 日本a在线网址|