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

    不同水鹽梯度下功能多樣性和功能冗余對荒漠植物群落穩(wěn)定性的影響

    2018-01-04 05:34:31王恒方呂光輝周耀治
    生態(tài)學(xué)報 2017年23期
    關(guān)鍵詞:水鹽樣方梯度

    王恒方,呂光輝,周耀治,曹 靖

    1 新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院,烏魯木齊 830046 2 綠洲生態(tài)教育部重點(diǎn)實(shí)驗室,烏魯木齊 830046 3 新疆大學(xué)干旱生態(tài)環(huán)境研究所,烏魯木齊 830046

    不同水鹽梯度下功能多樣性和功能冗余對荒漠植物群落穩(wěn)定性的影響

    王恒方1,2,呂光輝2,3,*,周耀治1,2,曹 靖1,2

    1 新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院,烏魯木齊 830046 2 綠洲生態(tài)教育部重點(diǎn)實(shí)驗室,烏魯木齊 830046 3 新疆大學(xué)干旱生態(tài)環(huán)境研究所,烏魯木齊 830046

    功能多樣性和功能冗余是物種多樣性的兩個組成部分,也是影響群落穩(wěn)定性的兩個重要因素?;诓煌}梯度下植物功能多樣性、功能冗余、物種多樣性和群落穩(wěn)定性及其相關(guān)關(guān)系的計算結(jié)果,分析功能多樣性和功能冗余對群落穩(wěn)定性的影響,結(jié)果表明: (1) 功能多樣性、物種多樣性和群落穩(wěn)定性均表現(xiàn)為高水高鹽和中水中鹽群落顯著高于低水低鹽群落(P<0.05);(2) 高水高鹽群落,功能多樣性與物種多樣性的相關(guān)系數(shù)小于功能冗余與物種多樣性的相關(guān)系數(shù),且功能多樣性與穩(wěn)定性的相關(guān)系數(shù)也小于功能冗余與穩(wěn)定性的相關(guān)系數(shù),而中水中鹽和低水低鹽群落的相關(guān)系數(shù)則呈現(xiàn)相反的規(guī)律;(3) 中水中鹽和低水低鹽群落的功能多樣性的標(biāo)準(zhǔn)化偏回歸系數(shù)均大于功能冗余的標(biāo)準(zhǔn)化偏回歸系數(shù);(4) 典范對應(yīng)分析中,土壤含水量可以解釋總特征根的22.7%,而土壤含鹽量僅可以解釋總特征根的1.3%;(5) 高水高鹽群落的穩(wěn)定比最接近20/80,穩(wěn)定性最高;低水低鹽群落遠(yuǎn)離20/80,穩(wěn)定性最低。改進(jìn)后的Godron穩(wěn)定性測定方法與物種種群密度變異系數(shù)方法得出的結(jié)果相同。綜上可知,功能多樣性和功能冗余兩者中與物種多樣性關(guān)系更為密切者對群落穩(wěn)定性的影響也越大,且兩者均可提高群落穩(wěn)定性,也就證明冗余假說在溫帶干旱荒漠區(qū)域的隱域性植物群落中是成立的;群落穩(wěn)定性、功能多樣性、功能冗余及物種多樣性主要是受土壤含水量的影響,土壤含鹽量對其影響較小。

    物種多樣性;Rao′s指數(shù);荒漠植物群落;冗余假說

    隨著全球生物多樣性的迅速喪失,物種多樣性和生態(tài)系統(tǒng)穩(wěn)定性的關(guān)系以及他們的維持機(jī)制已經(jīng)成為人類面臨的重要科學(xué)問題[1- 4]。影響生態(tài)系統(tǒng)穩(wěn)定性的因子眾多,生態(tài)學(xué)家關(guān)注最多的是物種多樣性對其的影響[4]。人們提出很多關(guān)于物種多樣性和生態(tài)系統(tǒng)功能的假說,包括多樣性-穩(wěn)定性假說[2],鉚釘假說[5],物種冗余假說[6],不確定假說[7]和關(guān)鍵種假說[8]等。研究表明物種多樣性和穩(wěn)定性間可能存在正相關(guān)[1- 3,9]、負(fù)相關(guān)[10]和非線性的關(guān)系[11],并沒有一致的結(jié)論。Bello等認(rèn)為物種多樣性包括功能多樣性和功能冗余兩個部分[12],功能多樣性和功能冗余中任何一個因子的變化都可能導(dǎo)致生態(tài)系統(tǒng)穩(wěn)定性的變化。

    功能多樣性是聯(lián)系物種多樣性和生態(tài)系統(tǒng)功能的關(guān)鍵性因素[13],是表示種間功能屬性相對差異性和生態(tài)互補(bǔ)性的指標(biāo),通常認(rèn)為在決定生態(tài)系統(tǒng)功能上比物種多樣性有更直接的作用。植物的功能多樣性在闡明環(huán)境變化、群落組成及與生態(tài)系統(tǒng)的關(guān)系中起到了重要的作用,功能多樣性包括功能豐富度、功能均勻度和功能分離度。常用的方法是以Rao′s功能多樣性指數(shù)代表功能多樣性,Rao′s二次方程整合了物種豐富度和物種之間的功能特征差異??梢院芎玫慕忉屛锓N間趨異性和物種在群落中的比例的表征[14]。通常認(rèn)為當(dāng)功能多樣性隨物種多樣性增加表現(xiàn)為飽和增加時,就預(yù)示著存在功能冗余。Walker[6]在冗余假說中提出在一個包含許多物種的功能群中,物種通常表現(xiàn)出功能冗余。也就是說,一些物種呈現(xiàn)出類似的或相同的功能,并且對環(huán)境的變化有異步響應(yīng),或顯示出時間生態(tài)位分化的現(xiàn)象[6,15],這些物種在群落受到干擾后有助于保持群落的穩(wěn)定性[16- 18]。Bello等[12]提出了一種計算功能冗余的方法,并認(rèn)為功能冗余是物種多樣性除去功能多樣性,即物種多樣性與功能多樣性的差異。

    本文主要研究了不同水鹽梯度下荒漠區(qū)植被群落功能多樣性和功能冗余的關(guān)系。在干旱區(qū)的空間資源再分配中的研究中,水鹽是影響多樣性和生態(tài)系統(tǒng)功能的重要的環(huán)境因素,分析了研究區(qū)植被群落的多樣性和群落穩(wěn)定性的變化,以及他們在不同水鹽梯度下的相互關(guān)系。試圖回答以下問題:(1)冗余假說(冗余可以提高穩(wěn)定性)在該研究區(qū)是否成立。(2)功能多樣性和功能冗余哪個對穩(wěn)定性的影響更大。(3)水分和鹽分哪個因素對群落的功能多樣性、功能冗余及群落穩(wěn)定性影響更大。該研究除了可以加深人們對多樣性和穩(wěn)定性關(guān)系的了解外,也可以為荒漠區(qū)植被的保護(hù)和恢復(fù)提供科學(xué)依據(jù)。

    1 研究方法

    1.1 研究區(qū)概況

    研究區(qū)位于艾比湖國家級濕地自然保護(hù)區(qū)(79°53′—85°02′ E,43°38′—45°52′ N)。該區(qū)域?qū)俦睖貛Т箨懶愿珊禋夂?年均降水量為107 mm,年均蒸發(fā)量1600 mm以上,日照時數(shù)約2800 h,地處溫帶又位于阿拉山口大風(fēng)的主風(fēng)道上,年平均大風(fēng)(大于17 m/s)日數(shù)多達(dá)164 d。土壤類型豐富多樣,典型的土壤為灰漠土、灰棕漠土和風(fēng)沙土[19],多樣化的土壤類型孕育了豐富的植物群落,荒漠旱生植物物種尤其豐富,主要植物種類有胡楊(Populuseuphratica)、梭梭(Sacsaoul)、蘆葦(Reed)、檉柳(Tamarixchinensis)、鹽穗木(Halostachyscaspica)、鹽爪爪(Kalidiumfoliatum)、琵琶柴(Reaumuriasongonica)和駱駝刺(Alhagisparsifolia)等。

    1.2 實(shí)驗設(shè)計

    1.2.1 樣方設(shè)置

    在艾比湖國家級濕地自然保護(hù)區(qū)的荒漠區(qū)內(nèi),按土壤水分和鹽分可能沿著離河距離的變化梯度,垂直于艾比湖補(bǔ)給河流之一的阿其克蘇河,沿河岸在保護(hù)區(qū)東大橋管護(hù)站以北布設(shè)樣地,樣地東西相距2 km,南北相距3 km,共6 km2。將樣地平均劃分為30個400 m×500 m小樣地,在每個小樣地的中心點(diǎn)設(shè)置一個30 m×30 m的樣方,沿從東到西方向共5條樣帶,每條樣帶6個樣方,總計30個30 m×30 m的樣方。

    圖1 研究區(qū)及樣方分布示意圖Fig.1 The study area and location of the plots

    1.2.2 取樣測定

    在每個30 m×30 m的樣方中調(diào)查喬木,沿對角線取2個5 m×5 m的灌木樣方,在每個喬木樣方的對角線上各取2個1 m×1 m的草本樣方。調(diào)查時主要記錄喬木植物的種名、多度、胸徑、株高、冠幅,灌木和草本植物調(diào)查種名、叢數(shù)、叢幅、高度和蓋度。同時記錄各樣方的經(jīng)緯度。

    在每個樣方中,選取每個物種的5片不同健康葉片,利用游標(biāo)卡尺測定葉片厚度,同時用便攜式葉綠素檢測儀(SPAD- 502)測定植物葉片的相對葉綠素含量。用葉面積儀測定葉面積,梭梭、沙拐棗等形狀近似圓柱體的葉片,利用游標(biāo)卡尺測定葉片直徑和長度,按圓柱體表面積計算公式求出葉片總面積。然后帶回實(shí)驗室烘干,用萬分之一天平(AL204,Mettler Toledo Group,上海)稱重,計算平均葉面積(m2)、比葉重(kg/m2)和比葉面積(m2/kg)。

    每個樣方中,采用“S”采樣法,選取3個土壤樣點(diǎn)(包含一份喬木冠下土、一份灌木冠下土以及一份裸地土)采集0—15cm土樣,除去凋落物,分別裝入已稱重的鋁盒和密封袋。立即稱量鋁盒和土壤的質(zhì)量,帶回實(shí)驗室測定土壤含水量和土壤含鹽量。

    表1 不同水鹽梯度樣地的植物種名錄

    1.3 數(shù)據(jù)計算與統(tǒng)計分析

    1.3.1利用SPSS 21.0,采用系統(tǒng)聚類中的wards方法對樣方進(jìn)行系統(tǒng)聚類,利用多重比較的Duncan檢驗方法,對土壤的含水量及含鹽量等理化性質(zhì)進(jìn)行差異性檢驗(表2)。

    表2 不同水鹽梯度土壤理化性質(zhì)特征

    注:同一列相同字母表示不同水鹽梯度下土壤理化性質(zhì)無顯著差異(P>0.05)

    1.3.2 功能多樣性計算

    通過群落調(diào)查,以19個物種(表1)的4個定性性狀(科、生活周期、葉表面類型、生長型)和6個定量性狀(冠幅、株高、葉厚度、比葉重、比葉面積和相對葉綠素含量)共10個功能特征類型(表3),用于計算功能多樣性。在計算功能多樣性時,對物種×性狀狀態(tài)數(shù)據(jù)矩陣中每個物種定性性狀進(jìn)行賦值分類(表3)。

    采用Rao′s指數(shù)對30個樣方的功能多樣性進(jìn)行計算,由于物種功能特征值是數(shù)量、名稱的混合型數(shù)據(jù),所以計算中用到的距離公式都采用Gower′s距離公式[20-21],使用R(3.2.5)軟件進(jìn)行計算。Rao′s計算公式如下[22]:

    (1)

    S為樣方內(nèi)的物種數(shù);pi和pj分別為第i和第j個物種占總物種個體數(shù)的比例;dij為Gower′s距離,dij在0(兩物種具有完全相同的特征)至1(兩物種具有完全不同的特征)之間,表示物種i,j在一組性狀空間中的相異性。Rao系數(shù)描述了種間功能性狀的相對差異,比其他功能多樣性指數(shù)更適合表征群落的種間趨異性,也能更好的將物種喪失與生態(tài)系統(tǒng)功能變化相聯(lián)系。

    1.3.3 功能冗余FR計算

    采用Bello等[12]的公式計算:

    FR=D-Q

    (2)

    式中,D為Shannon-Wiener多樣性指數(shù),代表物種多樣性;Q為表示功能多樣性的Rao系數(shù),代表功能多樣性。D的計算公式如下[23]:

    (3)

    式中,S為樣方內(nèi)的物種數(shù);pi為第i個物種占總物種個體數(shù)的比例。

    表3 植物功能特征類型及賦值

    1.3.4 群落穩(wěn)定性的兩種計算方法

    (1)群落穩(wěn)定性以物種種群密度變異系數(shù)(coefficient of variation,CV)的倒數(shù)ICV表示[24- 26]:

    (4)

    式中,μ為樣方中各物種的平均密度,σ為各物種密度的標(biāo)準(zhǔn)差。ICV值越大,群落穩(wěn)定性越高,因為相對于平均密度,各物種密度的變異性較小。

    (2)以改進(jìn)后的Godron穩(wěn)定性測定方法[27]來測定群落的穩(wěn)定性。具體方法為:首先把樣地中不同種植物的蓋度由大到小的順序排列,換算成相對蓋度,按相對蓋度由大到小的順序逐步累加起來;然后將樣地內(nèi)植物種類的總和取倒數(shù),按植物種類的順序逐步累積;再將植物種類倒數(shù)百分率同累積相對蓋度一一對應(yīng),畫出散點(diǎn)圖并用平滑曲線連接起來。然后再與直線y=100-x相交,交點(diǎn)即為所求群落的穩(wěn)定性比值。一般地,穩(wěn)定性比值越接近20/80,群落就越穩(wěn)定。

    1.3.5 統(tǒng)計分析

    統(tǒng)計分析采用SPSS 21.0軟件進(jìn)行。采用線性回歸模型擬合功能多樣性-物種多樣性、功能冗余-物種多樣性、功能多樣性-功能冗余、群落穩(wěn)定性-功能冗余、群落穩(wěn)定性-功能多樣性以及群落穩(wěn)定性-物種多樣性間的關(guān)系,以F檢驗顯著(P<0.05)且決定系數(shù)R2來進(jìn)行比較分析。

    采用二元逐步回歸分析,以群落穩(wěn)定性為因變量,以功能多樣性和功能冗余為自變量,確定功能多樣性和功能冗余對穩(wěn)定性的相對影響,影響程度以標(biāo)準(zhǔn)化的偏回歸系數(shù)確定。

    為探討功能多樣性、功能冗余、物種多樣性及群落穩(wěn)定性與土壤含水量及含鹽量之間的關(guān)系,首先利用DCA(除趨勢對應(yīng)分析),查看各排序軸的長度(Axis lengths),其中最大值為1.325,小于3,應(yīng)優(yōu)先選用線性RDA模型,但根據(jù)模型效果最優(yōu)的原則,最終選取CCA排序分析[28-29]。以上計算及分析采用R(3.2.5)軟件進(jìn)行分析。

    在改進(jìn)后的Godron穩(wěn)定性測定方法中繪制平滑曲線時,由于可能存在多種關(guān)系如線性、飽和增加及S型曲線等,我們采用線性、對數(shù)、指數(shù)、冪函數(shù)四種模型擬合了群落中植物種類倒數(shù)百分率與累積相對蓋度的關(guān)系,以F檢驗顯著(P<0.05)且決定系數(shù)R2值最大的模型為最優(yōu)模型。

    2 結(jié)果與分析

    2.1 不同水鹽梯度下,群落的功能多樣性、功能冗余、物種多樣性和群落穩(wěn)定性

    不同水鹽梯度下,功能多樣性、物種多樣性及群落穩(wěn)定性均表現(xiàn)為高水高鹽群落>中水中鹽群落>低水低鹽群落。除功能冗余不同水鹽梯度下無顯著性差異外(P>0.05),功能多樣性、物種多樣性和群落穩(wěn)定性均表現(xiàn)為高水高鹽和中水中鹽群落顯著高于低水低鹽群落(P<0.05)(表4)。

    表4不同水鹽梯度的功能多樣性、功能冗余、物種多樣性及群落穩(wěn)定性變化

    Table4Changesoffunctionaldiversity,functionalredundancy,speciesdiversityandcommunitystabilityunderdifferentwater-salinityhabitats

    變量Variables功能多樣性Functionaldiversity功能冗余Functionalredundancy物種多樣性Speciesdiversity群落穩(wěn)定性Communitystability高水高鹽0.7201±0.0999b0.3267±0.1298a1.0467±0.2125b0.5048±0.0930b中水中鹽0.5763±0.1741b0.3300±0.1533a0.9063±0.2745b0.4351±0.02885b低水低鹽0.3245±0.1974a0.2478±0.0907a0.5723±0.2752a0.3245±0.02541a

    同一列相同字母表示不同水鹽梯度下功能多樣性、功能冗余、物種多樣性及群落穩(wěn)定性無顯著差異(P>0.05)

    2.2 水鹽梯度上功能多樣性、功能冗余、物種多樣性與群落穩(wěn)定性之間的回歸分析

    線性回歸分析顯示(圖2),在高水高鹽梯度上,功能多樣性與物種多樣性的相關(guān)性0.8131(由相關(guān)系數(shù)R2決定)小于功能冗余與物種多樣性的相關(guān)系數(shù)0.8924,功能多樣性與群落穩(wěn)定性的相關(guān)性0.770也小于功能冗余與群落穩(wěn)定性的相關(guān)性0.8453;而中水中鹽和低水低鹽梯度上,功能多樣性與物種多樣性的相關(guān)性大于功能冗余與物種多樣性的相關(guān)性,功能多樣性與穩(wěn)定性的相關(guān)性也大于功能冗余與穩(wěn)定性的關(guān)系。由圖2可知,隨著水鹽梯度的降低,功能多樣性與穩(wěn)定性的相關(guān)系數(shù)增大,功能冗余與高水高鹽穩(wěn)定性的相關(guān)性大于與中水中鹽及低水低鹽的相關(guān)性。不同水鹽梯度下,物種多樣性和群落穩(wěn)定性均呈正相關(guān)關(guān)系。

    圖2 不同水鹽梯度下物種多樣性、功能多樣性、功能冗余與群落穩(wěn)定性之間的線性關(guān)系Fig.2 Linear relations between species diversity, functional diversity, functional redundancy, and community stability under different water-salinity habitats

    逐步回歸分析結(jié)果顯示(表5),中水中鹽和低水低鹽梯度,功能多樣性和功能冗余的標(biāo)準(zhǔn)化偏回歸系數(shù)均大于0,說明增加功能多樣性和功能冗余均能提高群落穩(wěn)定性,決定系數(shù)R2的值大小表明能共同解釋群落穩(wěn)定性77.2%—80.3%的變化。高水高鹽群落,由于功能多樣性的貢獻(xiàn)性相對于功能冗余過小,在逐步回歸過程中剔除。中水中鹽的功能多樣性和功能冗余的標(biāo)準(zhǔn)化偏回歸系數(shù)分別為:0.373和0.235;低水低鹽的功能多樣性和功能冗余的標(biāo)準(zhǔn)化偏回歸系數(shù)分別為:0.357和0.243;可見中水中鹽和低水低鹽的功能多樣性的標(biāo)準(zhǔn)化偏回歸系數(shù)均大于功能冗余的標(biāo)準(zhǔn)化偏回歸系數(shù),因此功能多樣性對群落穩(wěn)定性的影響大于功能冗余的影響。由于各逐步回歸的自變量的容忍度值較大,功能多樣性和功能冗余之間不存在共線性(表5)。

    表5不同水鹽梯度上功能多樣性(x1)和功能冗余(x2)對群落穩(wěn)定性(y)影響的逐步回歸分析

    Table5Stepwiseregressiveanalysisoftheeffectsoffunctionaldiversity(x1)andfuctionalredundancy(x2)oncommunitystability(y)underdifferentwater-salinityhabitats

    水鹽梯度Waterandsaltgradients回歸方程RegressionequationsR2F-testP容忍度Tolerance高水高鹽y=0.292+0.652x20.80124.370<0.050.501中水中鹽y=0.143+0.373x1+0.235x20.803147.014<0.0010.804低水低鹽y=0.148+0.357x1+0.243x20.772332.293<0.0010.363

    2.3 典范對應(yīng)分析(CCA)

    圖3 各樣方的功能多樣性、功能冗余及穩(wěn)定性與水鹽的典范對應(yīng)分析(CCA)排序圖Fig.3 Canonical correspondence analysis (CCA) ordination diagram of 30 plots and species diversity, functional diversity, functional redundancy and community stability□代表高水高鹽群落;○代表中水中鹽群落;△代表低水低鹽群落;W代表水 分;S代表鹽分;D代表物種多樣性;Q代表功能多樣性;FR代表功能冗余;ICV代表群落穩(wěn)定性

    在CCA排序中(圖3),主要分析整體上物種多樣性、功能多樣性、功能冗余及其與群落穩(wěn)定性在空間上的含義,其中各排序軸包含物種多樣性、功能多樣性、功能冗余及穩(wěn)定性和樣方水鹽矩陣中所涵蓋的信息量的大小,可以通過其特征值體現(xiàn)。在排序過程中可以看出兩條排序軸的信息量,CCA排序的第一軸特征根值占總特征根值的22.7%,即第一軸能夠解釋樣方與水分的相關(guān)關(guān)系,它與土壤含水量有極顯著相關(guān)關(guān)系,說明第一軸主要反映土壤含水量對穩(wěn)定性、功能多樣性及功能冗余的影響,第二軸與鹽分存在相關(guān)關(guān)系,但僅能解釋1.3%,說明土壤含水量對功能多樣性、功能冗余和穩(wěn)定性的影響大于土壤含鹽量對其影響。由排序圖也可以看出整體來講物種多樣性、功能多樣性與群落穩(wěn)定性的相關(guān)性較強(qiáng),功能冗余與群落穩(wěn)定性的相關(guān)性較弱。

    2.4 改進(jìn)后的Godron穩(wěn)定性測定方法

    從穩(wěn)定性分析結(jié)果可見(表6),高水高鹽的穩(wěn)定比最接近20/80,其穩(wěn)定性最高,低水低鹽的穩(wěn)定比遠(yuǎn)離20/80,其穩(wěn)定性最低,中水中鹽居于兩者之間。

    3 討論

    3.1 功能多樣性和功能冗余對群落穩(wěn)定性的影響

    大量研究表明,物種多樣性對群落穩(wěn)定性具有正效應(yīng)[24,29],而物種多樣性對群落穩(wěn)定性的正效應(yīng)是通過功能多樣性和功能冗余的增加實(shí)現(xiàn)的,且功能多樣性是起互補(bǔ)效應(yīng)的作用,功能冗余是起保險效應(yīng)的作用[30]。Bello等研究結(jié)果表明,物種多樣性與功能多樣性和功能冗余正相關(guān)[12]。本研究為此提供了新的證據(jù),結(jié)果表明,不同水鹽梯度下,物種多樣性與群落穩(wěn)定性正相關(guān),功能多樣性和功能冗余對群落穩(wěn)定性的影響大小不同。高水高鹽條件下,通過線性回歸和逐步回歸(圖2、表5)的結(jié)果表明功能多樣性對群落穩(wěn)定性的影響很小,而功能冗余對群落穩(wěn)定性的影響占主導(dǎo)地位,這與Pillar得出的結(jié)論群落穩(wěn)定性僅由功能冗余所決定相一致[31],可能是在現(xiàn)有土壤水鹽條件下,物種相對豐富,干旱區(qū)物種多樣性已達(dá)到飽和狀態(tài),從而功能冗余對群落穩(wěn)定性的影響相對來講更大。Petchey等認(rèn)為[14],一個群落的功能冗余程度由內(nèi)在的冗余和外在的冗余共同決定,內(nèi)在冗余指種間的功能相似性,而外在冗余產(chǎn)生的原因是物種組成的非隨機(jī)變化引起的植物功能性狀的非隨機(jī)變化。如駱駝等牲畜選擇性的采食適口性物種以及放牧對競爭優(yōu)勢物種的抑制效應(yīng)則會增加物種的相對多度或多樣性,引起性狀組成的非隨機(jī)變化,產(chǎn)生外在冗余。功能冗余是群落本身的內(nèi)在屬性,因此,只有在群落物種多樣性趨于飽和后,功能冗余才會成為影響群落穩(wěn)定性的主要因素。而中水中鹽和低水低鹽條件下,功能多樣性與物種多樣性的相關(guān)性強(qiáng)于功能冗余與物種多樣性的相關(guān)性,說明物種間通過更強(qiáng)的功能互補(bǔ)性維持群落穩(wěn)定性;功能多樣性與群落穩(wěn)定性的相關(guān)性也強(qiáng)于功能冗余與群落穩(wěn)定性的相關(guān)性更進(jìn)一步說明了互補(bǔ)效應(yīng)的作用大于保險效應(yīng)的作用,這與姚天華等的研究結(jié)果一致[32],即群落物種功能多樣性和功能冗余對群落穩(wěn)定性的影響中,功能多樣性對群落穩(wěn)定性的影響較大。得出種間功能性狀差異的互補(bǔ)效應(yīng)才是維持群落穩(wěn)定性的主要動力。綜合以上結(jié)果及證據(jù),我們認(rèn)為群落物種功能多樣性和功能冗余對擾動發(fā)生響應(yīng)變化的程度確實(shí)與其對群落穩(wěn)定性影響的相對強(qiáng)弱正相關(guān),即功能多樣性和功能冗余兩者中與物種多樣性關(guān)系更為密切者對群落穩(wěn)定性的影響也將越大。

    表6 不同水鹽梯度下Godron曲線和交點(diǎn)坐標(biāo)

    在本研究區(qū),高水高鹽群落在逐步回歸結(jié)果中,功能多樣性對群落穩(wěn)定性的貢獻(xiàn)很小,而且在線性回歸模型中,功能多樣性對群落穩(wěn)定性的相關(guān)性小于功能冗余對群落穩(wěn)定性的相關(guān)性,我們推測可能是因為干旱荒漠區(qū)物種多樣性較低,高水高鹽群落的物種多樣性已經(jīng)趨于飽和,功能冗余開始發(fā)揮主要作用。但目前研究區(qū)物種飽和的精確分析還需進(jìn)一步的試驗研究。

    3.2 CCA分析對穩(wěn)定性及環(huán)境的解釋

    CCA排序中同時具有樣方的功能多樣性、功能冗余、物種多樣性及群落穩(wěn)定性矩陣,也有各樣方的水鹽矩陣,一方面可以解釋整體群落中功能多樣性、功能冗余、物種多樣性和穩(wěn)定性的關(guān)系,同時也可以較好的表達(dá)出各指標(biāo)與水鹽梯度的關(guān)系??傮w而言,樣方間水鹽越接近,樣方間的功能多樣性、功能冗余、物種多樣性及穩(wěn)定性越接近,這是群落的功能多樣性、功能冗余、物種多樣性及穩(wěn)定性與水鹽之間形成的生態(tài)關(guān)系。本CCA排序中,因只有兩個環(huán)境因子即水分和鹽分含量,第一軸與水分含量有關(guān),從左到右表達(dá)了水分遞減的一個梯度;第二軸主要與鹽分含量有關(guān),從下到上表達(dá)了鹽分遞增的梯度。綜合來說,可能是土壤含鹽量的范圍不足以對功能多樣性、功能冗余、物種多樣性及群落穩(wěn)定性產(chǎn)生限制性的影響(表2),土壤含水量是影響功能多樣性、功能冗余、物種多樣性及穩(wěn)定性的主要的環(huán)境因子。這可能與植物對土壤水分的生態(tài)適應(yīng)性相關(guān)。水分脅迫可能是導(dǎo)致不同水鹽梯度下功能多樣性和功能冗余的主要原因。這與張雪妮等[33]的研究結(jié)果一致,即艾比湖荒漠植物多樣性在土壤含鹽量為0.03%—1.24%的范圍內(nèi)植物的生長不受影響,干旱脅迫可能是導(dǎo)致兩種生境下植物多樣性指數(shù)差異顯著的主要原因。在各個樣方的排序圖上沿著CCA第一排序軸出現(xiàn)了較規(guī)律的分布,且出現(xiàn)了交錯的現(xiàn)象,顯示出群落的功能多樣性、功能冗余及穩(wěn)定性與其所處環(huán)境之間相互關(guān)系的復(fù)雜性。除了涉及水分和鹽分差距,可能生物間的相互作用等也會影響這些群落特征。

    CCA排序中,物種多樣性與群落穩(wěn)定性的相關(guān)性最大,功能冗余與群落穩(wěn)定性的相關(guān)性小于功能多樣性與群落穩(wěn)定性的相關(guān)性。功能冗余受水鹽條件的影響較小,可能是物種長期與環(huán)境及種群之間的相互作用,形成了極為復(fù)雜的關(guān)系,不易受水鹽條件的影響而改變。而功能多樣性對水鹽變化較敏感。

    3.3 改進(jìn)后的Godron穩(wěn)定性結(jié)果分析

    干旱區(qū)荒漠植物群落穩(wěn)定性主要是受種間競爭、環(huán)境壓力和干擾活動等的影響[27],本研究表明,利用改進(jìn)的Godron穩(wěn)定性測定方法測得的三類群落的穩(wěn)定性與通過物種種群密度變異系數(shù)得到的結(jié)果一致,可能和種間聯(lián)接及物種多樣性有關(guān),穩(wěn)定性由高到低依次為:高水高鹽群落>中水中鹽群落>低水低鹽群落。可能是高水高鹽梯度土壤含鹽量在植物生長所需的限度內(nèi),水分相對充足,高水高鹽的物種多樣性最高,最適合物種的生存,其穩(wěn)定性最高;中水中鹽和低水低鹽梯度,土壤含水量較低,限制了某些物種的生長,穩(wěn)定性較低。

    4 結(jié)論

    綜上所述,從不同水鹽梯度荒漠植物群落穩(wěn)定性對功能多樣性和功能冗余響應(yīng)的研究得出以下結(jié)論:

    (1)高水高鹽梯度下,功能多樣性對群落穩(wěn)定性的影響小于功能冗余對群落穩(wěn)定性的影響;中水中鹽和低水低鹽梯度下,功能多樣性對群落穩(wěn)定性的影響大于功能冗余對群落穩(wěn)定性的影響。

    (2)功能多樣性和功能冗余均可以提高群落穩(wěn)定性。且兩者中與物種多樣性關(guān)系更為密切者對群落穩(wěn)定性的影響也越大,即冗余假說(冗余可以提高群落穩(wěn)定性)在該干旱荒漠區(qū)域的隱域性植物群落中是成立的。

    (3)群落的穩(wěn)定性、功能多樣性、功能冗余及物種多樣性受水鹽的協(xié)同作用,但土壤含水量大于土壤含鹽量的影響。

    [1] Elton C S. The Ecology of Invasions by Animals and Plants. Chicago: The University of Chicago Press, 2000.

    [2] MacArthur R. Fluctuations of animal populations and a measure of community stability. Ecology. 1955, 36(3): 533- 536.

    [3] Mougi A, Kondoh M. Diversity of interaction types and ecological community stability. Science, 2012, 337(6092): 349- 351.

    [4] 張景慧, 黃永梅. 生物多樣性與穩(wěn)定性機(jī)制研究進(jìn)展. 生態(tài)學(xué)報, 2016, 36(13): 3859- 3870.

    [5] Ehrlich P R, Ehrlich A H. Extinction. New York: Ballantine. 1981: 384.

    [6] Walker B H. Biodiversity and ecological redundancy. Conservation biology, 1992, 6(1): 18- 23.

    [7] Lawton J H. What do species do in ecosystems? Oikos, 1994, 71(3): 367- 374.

    [8] Sala O E, Lauenroth W K, McNaughton S J, Rusch G, Zhang X S, Mooney H A, Cushman J H, Medina E, Schulz E D. Biodiversity and ecosystem functioning in grasslands// Mooney H A, Cushman J H, Medina E, Sala O E, Schulze E D, eds. Functional Roles of Biodiversity: A Global Perspective. Chichester: Wiley, 1996: 129- 149.

    [9] Tilman D, Reich P B, Knops J M H. Biodiversity and ecosystem stability in a decade-long grassland experiment. Nature, 2006, 441(7093): 629- 632.

    [10] Mori A S, Furukawa T, Sasaki T. Response diversity determines the resilience of ecosystems to environmental change. Biological Reviews, 2013, 88(2): 349- 364.

    [11] Schulze E D, Mooney H A. Biodiversity and Ecosystem Function. Berlin Heidelberg: Springer, 1994.

    [12] De Bello F, Lep? J, Lavorel S, Moretti M. Importance of species abundance for assessment of trait composition: an example based on pollinator communities. Community Ecology, 2007, 8(2): 163- 170.

    [13] 江小雷, 張衛(wèi)國. 功能多樣性及其研究方法. 生態(tài)學(xué)報, 2010, 30(10): 2766- 2773.

    [14] Petchey O L, Evans K L, Fishburn I S, Gaston K J. Low functional diversity and no redundancy in British avian assemblages. Journal of Animal Ecology, 2007, 76(5): 977- 985.

    [15] Elmqvist T, Folke C, Nystr?m M, Peterson G, Bengtsson J, Walker B, Norberg J. Response diversity, ecosystem change, and resilience. Frontiers in Ecology and the Environment, 2003, 1(9): 488- 494.

    [16] Naeem S. Species redundancy and ecosystem reliability. Conservation Biology, 1998, 12(1): 39- 45.

    [17] Thibaut L M, Connolly S R, Sweatman H P A. Diversity and stability of herbivorous fishes on coral reefs. Ecology, 2012, 93(4): 891- 901.

    [18] Gonzalez A, Loreau M. The causes and consequences of compensatory dynamics in ecological communities. Annual Review of Ecology, Evolution, and Systematics, 2009, 40: 393- 414.

    [19] 中國科學(xué)院新疆綜合考察隊, 中國科學(xué)院植物研究所. 新疆植被及其利用. 北京: 科學(xué)出版社, 1978.

    [20] Gower J C. A general coefficient of similarity and some of its properties. Biometrics, 1971, 27(4): 857- 871.

    [21] Podani J. Extending Gower′s general coefficient of similarity to ordinal characters. Taxon, 1999, 48(2): 331- 340.

    [22] Ricotta C. A note on functional diversity measures. Basic and Applied Ecology, 2005, 6(5): 479- 486.

    [23] Zhang J T. Quantitative Ecology. Beijing: Science Press, 2004.

    [24] Yang Z L, Van Ruijven J, Du G Z. The effects of long-term fertilization on the temporal stability of alpine meadow communities. Plant and Soil, 2011, 345(1- 2): 315- 324.

    [25] Sasaki T, Lauenroth W K. Dominant species, rather than diversity, regulates temporal stability of plant communities. Oeclogia, 2011, 166(3): 761- 768.

    [26] 王海東, 張璐璐, 朱志紅. 刈割、施肥對高寒草甸物種多樣性與生態(tài)系統(tǒng)功能關(guān)系的影響及群落穩(wěn)定性機(jī)制. 植物生態(tài)學(xué)報, 2013, 37(4): 279- 295.

    [27] 呂光輝, 杜昕, 楊建軍, 馬媛, 孟吉翔. 阜康綠洲——荒漠交錯帶荒漠植被群落穩(wěn)定性. 干旱區(qū)地理, 2007, 30(5): 660- 665.

    [28] 國際生物多樣性計劃中國委員會. 中國生物多樣性保護(hù)與研究進(jìn)展. 北京: 氣象出版社, 2014.

    [29] 張欽弟, 段曉梅, 白玉芳, 王祎玲, 張金屯. 山西太岳山脫皮榆群落的功能多樣性. 植物學(xué)報, 2016, 51(2): 218- 225.

    [30] Tilman D. The ecological consequences of changes in biodiversity: a search for general principles. Ecology, 1999, 80(5): 1455- 1474.

    [31] Pillar V D, Blanco C C, Müller S C, Sosinski E E, Joner F, Duarte L D S. Functional redundancy and stability in plant communities. Journal of Vegetation Science, 2013, 24(5): 963- 974.

    [32] 姚天華, 朱志紅, 李英年, 潘石玉, 孔彬彬, 衛(wèi)欣華, 杜家麗. 功能多樣性和功能冗余對高寒草甸群落穩(wěn)定性的影響. 生態(tài)學(xué)報, 2016, 36(6): 1547- 1558.

    [33] 張雪妮, 呂光輝, 楊曉東, 貢璐, 秦璐, 何學(xué)敏, 劉昊奇. 基于鹽分梯度的荒漠植物多樣性與群落、種間聯(lián)接響應(yīng). 生態(tài)學(xué)報, 2013, 33(18): 5714- 5722.

    Effectsoffunctionaldiversityandfunctionalredundancyonthestabilityofdesertplantcommunitiesunderdifferentwaterandsaltgradients

    WANG Hengfang1,2, Lü Guanghui2,3,*, ZHOU Yaozhi1,2, CAO Jing1,2

    1CollegeofResourcesandEnvironmentScience,XinjiangUniversity,Urumqi830046,China2KeyLaboratoryofOasisEcology,MinistryofEducation,Urumqi830046,China3InstituteofAridEcologyandEnvironment,XinjiangUniversity,Urumqi830046,China

    Functional diversity and functional redundancy are two components of species diversity, and are also two important factors that affect the stability of a community. The effects of functional diversity and functional redundancy on the community stability were analyzed based on the estimation of functional diversity, functional redundancy, species diversity, community stability, and their correlations of plant populations under different water and salt gradients The results revealed that: (1) The functional diversity, species diversity, and community stability of plant populations in high and moderate water and salt gradients were significantly higher than those of in low water and salt gradients (P<0.05). (2) For populations in high water and salinity gradient, the correlation coefficient between functional diversity and species diversity was lower than that of between functional redundancy and species diversity, moreover, the correlation coefficient between functional diversity and stability was lower than that of between functional redundancy and stability. However, the correlation coefficients of populations in moderate and low water and salinity gradients were contrary to those of in high water and salinity gradient. (3) The standardized partial regression coefficients of functional diversity were greater than that of functional redundancy of communities in moderate and low water and salt gradients. (4) In the canonical correspondence analysis, the water and salt content of soil could explain 22.7% and 1.3% of the total characteristic root, respectively. (5) The stability was highest in plant communities in high water and salt gradient, with the stability rate near 20/80, whereas populations in low water and salt gradient had the lowest stability, with the stability rate far from 20/80. These results were consistent with those of estimated by the improved Godron Stability Test Method and Population Density Variance Coefficient Method. In summary, both of the functional diversity and functional redundancy could enhance the stability of plant community, and their effects were dependent on their relationships with species diversity which provided evidence for the redundancy hypothesis of plant communities in the temperate arid desert region. The stability, functional diversity, functional redundancy, and species diversity of plant communities were much more affected by soil moisture, but less by soil salinity.

    species diversity; Rao′s index; desert plant community; redundancy hypothesi

    國家自然科學(xué)基金項目(31560131,41571034);國家重點(diǎn)基礎(chǔ)研究發(fā)展計劃資助項目(2015CB251601)

    2016- 10- 19; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 08- 14

    *通訊作者Corresponding author.E-mail: ler@xju.edu.cn

    10.5846/stxb201610192139

    王恒方,呂光輝,周耀治,曹靖.不同水鹽梯度下功能多樣性和功能冗余對荒漠植物群落穩(wěn)定性的影響.生態(tài)學(xué)報,2017,37(23):7928- 7937.

    Wang H F, Lü G H, Zhou Y Z, Cao J.Effects of functional diversity and functional redundancy on the stability of desert plant communities under different water and salt gradients.Acta Ecologica Sinica,2017,37(23):7928- 7937.

    猜你喜歡
    水鹽樣方梯度
    藥王谷多種因素對漏蘆、野韭生長的影響
    鹽地堿蓬生長生理指標(biāo)和冠層光譜對水鹽脅迫的響應(yīng)
    昆明金殿國家森林公園林業(yè)有害生物調(diào)查研究
    一個改進(jìn)的WYL型三項共軛梯度法
    典型野生刺梨(Rosa roxburghii Tratt.)灌草叢植物多樣性研究
    不同供給源下粗粒土水鹽遷移特性分析
    水鹽體系汽液平衡的熱力學(xué)模型研究
    一種自適應(yīng)Dai-Liao共軛梯度法
    云南九龍山自然保護(hù)區(qū)水青岡群落物種多樣性及群落特征
    艾丁湖流域不同水鹽條件與天然植被關(guān)系
    国产片特级美女逼逼视频| 久久久精品免费免费高清| 伦理电影大哥的女人| 老鸭窝网址在线观看| 国产一区二区在线观看av| 在线观看美女被高潮喷水网站| 少妇人妻精品综合一区二区| 久久这里有精品视频免费| 狠狠婷婷综合久久久久久88av| 视频区图区小说| 国产人伦9x9x在线观看 | 最新的欧美精品一区二区| 亚洲精品乱久久久久久| 亚洲av免费高清在线观看| 新久久久久国产一级毛片| 精品酒店卫生间| 如日韩欧美国产精品一区二区三区| 多毛熟女@视频| 不卡av一区二区三区| 精品亚洲成a人片在线观看| av有码第一页| 日韩制服骚丝袜av| 春色校园在线视频观看| 天天影视国产精品| 成人毛片a级毛片在线播放| 亚洲一级一片aⅴ在线观看| 日本av手机在线免费观看| 黄网站色视频无遮挡免费观看| 日韩精品免费视频一区二区三区| 熟女电影av网| 亚洲国产欧美在线一区| 亚洲精品国产av成人精品| 自线自在国产av| www.av在线官网国产| 国产熟女午夜一区二区三区| 毛片一级片免费看久久久久| 亚洲精品一区蜜桃| 免费女性裸体啪啪无遮挡网站| 精品久久久久久电影网| 人妻人人澡人人爽人人| 国产精品久久久av美女十八| 成人毛片a级毛片在线播放| 亚洲av福利一区| 夜夜骑夜夜射夜夜干| 女人高潮潮喷娇喘18禁视频| 国产有黄有色有爽视频| 一本色道久久久久久精品综合| 久久精品国产综合久久久| 亚洲美女搞黄在线观看| 国产成人精品在线电影| 成年女人毛片免费观看观看9 | av在线老鸭窝| 香蕉精品网在线| 欧美激情 高清一区二区三区| videosex国产| 两性夫妻黄色片| 国产男女内射视频| 成人国产麻豆网| 中文字幕最新亚洲高清| 老司机亚洲免费影院| 99热全是精品| 丝袜脚勾引网站| 日韩中文字幕欧美一区二区 | 精品久久久久久电影网| 男女边吃奶边做爰视频| 久久人人爽av亚洲精品天堂| 中文字幕av电影在线播放| 免费日韩欧美在线观看| 欧美在线一区亚洲| 亚洲男人的天堂狠狠| a级毛片黄视频| av欧美777| 成人免费观看视频高清| 亚洲人成电影观看| 国产精品一区二区在线不卡| 超碰成人久久| 亚洲成人精品中文字幕电影 | 亚洲精品av麻豆狂野| 国产高清视频在线播放一区| av在线天堂中文字幕 | 欧美日韩亚洲高清精品| 黄色视频,在线免费观看| 亚洲一区高清亚洲精品| 亚洲精品一卡2卡三卡4卡5卡| 欧美黑人精品巨大| 久久午夜亚洲精品久久| netflix在线观看网站| 日韩欧美国产一区二区入口| 国产精品av久久久久免费| 亚洲免费av在线视频| 亚洲va日本ⅴa欧美va伊人久久| 久久国产精品男人的天堂亚洲| 精品国产乱码久久久久久男人| 中文字幕色久视频| 一二三四在线观看免费中文在| av视频免费观看在线观看| 十八禁人妻一区二区| 国产日韩一区二区三区精品不卡| 国产高清视频在线播放一区| 最近最新中文字幕大全免费视频| 国内久久婷婷六月综合欲色啪| 在线观看免费视频日本深夜| 国产免费现黄频在线看| 亚洲第一欧美日韩一区二区三区| 国产精品亚洲一级av第二区| 亚洲av日韩精品久久久久久密| 亚洲成a人片在线一区二区| 欧美人与性动交α欧美软件| 亚洲精品一二三| 久久性视频一级片| 9191精品国产免费久久| 国产免费av片在线观看野外av| 嫁个100分男人电影在线观看| 日韩大尺度精品在线看网址 | 国产精品爽爽va在线观看网站 | 亚洲七黄色美女视频| 欧美日韩亚洲国产一区二区在线观看| 一区二区三区激情视频| 啦啦啦 在线观看视频| 老司机午夜十八禁免费视频| 久久精品亚洲精品国产色婷小说| 欧美午夜高清在线| 99久久国产精品久久久| 9热在线视频观看99| 精品国产美女av久久久久小说| 久久性视频一级片| 日韩免费高清中文字幕av| avwww免费| 久久久久九九精品影院| 韩国av一区二区三区四区| 精品国产一区二区三区四区第35| 嫁个100分男人电影在线观看| 亚洲色图综合在线观看| 在线观看一区二区三区| 久99久视频精品免费| 9191精品国产免费久久| 18禁国产床啪视频网站| 日韩一卡2卡3卡4卡2021年| 中文字幕精品免费在线观看视频| 欧美激情 高清一区二区三区| 99热国产这里只有精品6| 午夜福利影视在线免费观看| 日本精品一区二区三区蜜桃| 久久精品国产99精品国产亚洲性色 | 欧美乱妇无乱码| 色综合婷婷激情| 婷婷精品国产亚洲av在线| 男女之事视频高清在线观看| xxx96com| 啦啦啦免费观看视频1| 久久久久国内视频| 亚洲av成人一区二区三| 欧美老熟妇乱子伦牲交| 女人精品久久久久毛片| 身体一侧抽搐| 69av精品久久久久久| 成人18禁高潮啪啪吃奶动态图| 欧美日韩亚洲高清精品| 91av网站免费观看| 亚洲 欧美一区二区三区| 亚洲精品在线观看二区| 两性夫妻黄色片| 韩国av一区二区三区四区| 欧美黑人精品巨大| aaaaa片日本免费| 又黄又爽又免费观看的视频| 视频区欧美日本亚洲| 成人亚洲精品av一区二区 | 国产亚洲欧美98| 欧美日韩精品网址| 一二三四社区在线视频社区8| 极品人妻少妇av视频| 在线免费观看的www视频| 女警被强在线播放| 亚洲精品成人av观看孕妇| 国产一区二区在线av高清观看| 中文字幕色久视频| 三上悠亚av全集在线观看| 欧美激情高清一区二区三区| 我的亚洲天堂| 在线观看一区二区三区激情| 操出白浆在线播放| 久热爱精品视频在线9| 天堂√8在线中文| 亚洲一区高清亚洲精品| 国产99白浆流出| 黄色视频,在线免费观看| 日本三级黄在线观看| 色哟哟哟哟哟哟| 精品国产乱码久久久久久男人| 成人影院久久| 99久久人妻综合| 麻豆av在线久日| 美女福利国产在线| 亚洲男人天堂网一区| av片东京热男人的天堂| 黄色 视频免费看| 黄频高清免费视频| 日日夜夜操网爽| 日韩 欧美 亚洲 中文字幕| 欧美乱妇无乱码| 国产免费男女视频| 国产精品九九99| 亚洲精品国产一区二区精华液| 18禁国产床啪视频网站| 久久精品91无色码中文字幕| 又黄又粗又硬又大视频| 欧美老熟妇乱子伦牲交| 美女扒开内裤让男人捅视频| 亚洲精品在线观看二区| 香蕉国产在线看| 老司机深夜福利视频在线观看| 久久中文看片网| 色婷婷av一区二区三区视频| 不卡一级毛片| 欧美在线一区亚洲| 日韩免费高清中文字幕av| 88av欧美| 在线播放国产精品三级| 欧美在线黄色| 很黄的视频免费| 亚洲五月色婷婷综合| 久久久久久大精品| 黄色怎么调成土黄色| 国产真人三级小视频在线观看| tocl精华| 侵犯人妻中文字幕一二三四区| 久久久国产精品麻豆| 午夜久久久在线观看| 自线自在国产av| 久久精品91无色码中文字幕| 女性被躁到高潮视频| 91精品三级在线观看| 国产亚洲欧美精品永久| 国产亚洲精品第一综合不卡| 一进一出抽搐动态| 丝袜在线中文字幕| 这个男人来自地球电影免费观看| 亚洲欧美一区二区三区久久| 深夜精品福利| 成人影院久久| 精品人妻在线不人妻| 国产麻豆69| 一边摸一边做爽爽视频免费| 国产精品久久久久久人妻精品电影| 热re99久久国产66热| 日韩欧美一区二区三区在线观看| 亚洲第一欧美日韩一区二区三区| 国产欧美日韩一区二区三区在线| 青草久久国产| 亚洲专区国产一区二区| 亚洲欧美一区二区三区黑人| netflix在线观看网站| 亚洲男人天堂网一区| 91麻豆av在线| 新久久久久国产一级毛片| 午夜精品在线福利| 麻豆成人av在线观看| 99国产综合亚洲精品| 黄色a级毛片大全视频| av在线播放免费不卡| av福利片在线| 老汉色∧v一级毛片| 啪啪无遮挡十八禁网站| 超碰成人久久| 久久人妻福利社区极品人妻图片| 日韩大尺度精品在线看网址 | a级毛片在线看网站| 女人爽到高潮嗷嗷叫在线视频| 88av欧美| 多毛熟女@视频| 午夜免费成人在线视频| 99久久国产精品久久久| 亚洲国产欧美日韩在线播放| 久久久久国内视频| 黄色成人免费大全| av网站免费在线观看视频| 91在线观看av| 国产精品日韩av在线免费观看 | 又紧又爽又黄一区二区| 男女下面插进去视频免费观看| 亚洲av成人一区二区三| 五月开心婷婷网| 在线播放国产精品三级| 国产主播在线观看一区二区| www国产在线视频色| 天天躁夜夜躁狠狠躁躁| 操出白浆在线播放| 悠悠久久av| 麻豆av在线久日| 亚洲国产欧美网| 人人妻人人爽人人添夜夜欢视频| av天堂久久9| 啦啦啦免费观看视频1| x7x7x7水蜜桃| 97碰自拍视频| 日本wwww免费看| 妹子高潮喷水视频| 国产黄色免费在线视频| 又紧又爽又黄一区二区| 18美女黄网站色大片免费观看| 欧美性长视频在线观看| 99国产综合亚洲精品| 黄色毛片三级朝国网站| 久久精品影院6| 国产高清国产精品国产三级| 久久精品国产亚洲av香蕉五月| 首页视频小说图片口味搜索| 亚洲成人免费av在线播放| 中文字幕人妻熟女乱码| 麻豆久久精品国产亚洲av | 午夜精品国产一区二区电影| 中文字幕人妻熟女乱码| ponron亚洲| 两性夫妻黄色片| 一级a爱视频在线免费观看| 欧美成人性av电影在线观看| 免费日韩欧美在线观看| 高清黄色对白视频在线免费看| 香蕉久久夜色| 99国产综合亚洲精品| 99国产精品99久久久久| 91成年电影在线观看| 欧美最黄视频在线播放免费 | 欧美中文综合在线视频| 久久精品aⅴ一区二区三区四区| 成年女人毛片免费观看观看9| 亚洲少妇的诱惑av| 最新在线观看一区二区三区| 国产男靠女视频免费网站| 中文字幕人妻熟女乱码| 大型黄色视频在线免费观看| 美女 人体艺术 gogo| 老司机靠b影院| 午夜免费激情av| 久久国产精品影院| 超色免费av| 亚洲欧美精品综合一区二区三区| 99国产精品一区二区蜜桃av| 国产av精品麻豆| 真人做人爱边吃奶动态| 国产精品爽爽va在线观看网站 | 很黄的视频免费| 伊人久久大香线蕉亚洲五| 中文字幕人妻熟女乱码| 自线自在国产av| 久99久视频精品免费| 黄色 视频免费看| 欧美成狂野欧美在线观看| 亚洲熟女毛片儿| 丝袜人妻中文字幕| 亚洲熟女毛片儿| 一个人观看的视频www高清免费观看 | 黄色毛片三级朝国网站| 亚洲av日韩精品久久久久久密| 窝窝影院91人妻| 91精品三级在线观看| 精品久久久久久久久久免费视频 | 老司机福利观看| 精品久久久久久,| 视频区图区小说| 91字幕亚洲| 日韩有码中文字幕| 伊人久久大香线蕉亚洲五| 黄色片一级片一级黄色片| 亚洲国产欧美日韩在线播放| 久99久视频精品免费| 亚洲午夜精品一区,二区,三区| 国产欧美日韩精品亚洲av| 精品一区二区三区av网在线观看| 一区二区日韩欧美中文字幕| 性少妇av在线| 国产精品影院久久| 久久精品国产99精品国产亚洲性色 | 超碰97精品在线观看| 两性夫妻黄色片| 超碰97精品在线观看| 波多野结衣一区麻豆| 亚洲精品国产区一区二| 别揉我奶头~嗯~啊~动态视频| 成人av一区二区三区在线看| 亚洲人成伊人成综合网2020| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久成人av| 午夜免费观看网址| 精品国产乱子伦一区二区三区| 岛国视频午夜一区免费看| 淫妇啪啪啪对白视频| 国产精品二区激情视频| 亚洲色图综合在线观看| 天天添夜夜摸| 超碰成人久久| 国产精品久久视频播放| 国产99白浆流出| 午夜精品久久久久久毛片777| 黄色 视频免费看| 日韩有码中文字幕| 国产三级在线视频| 欧美另类亚洲清纯唯美| 另类亚洲欧美激情| 黄色怎么调成土黄色| 在线天堂中文资源库| 亚洲色图综合在线观看| 亚洲性夜色夜夜综合| 女人爽到高潮嗷嗷叫在线视频| 午夜亚洲福利在线播放| 9热在线视频观看99| 亚洲一区二区三区欧美精品| 日韩中文字幕欧美一区二区| 亚洲人成电影免费在线| 久久人人精品亚洲av| 国产成人欧美在线观看| 两个人免费观看高清视频| 高清在线国产一区| 亚洲男人天堂网一区| 欧美最黄视频在线播放免费 | 精品国产乱子伦一区二区三区| 亚洲精品国产精品久久久不卡| 色综合婷婷激情| 日韩欧美免费精品| 亚洲午夜精品一区,二区,三区| 久久精品91无色码中文字幕| 老汉色∧v一级毛片| 视频区图区小说| 国产三级在线视频| √禁漫天堂资源中文www| av视频免费观看在线观看| 亚洲情色 制服丝袜| 亚洲av成人av| 国产熟女xx| 51午夜福利影视在线观看| 又黄又爽又免费观看的视频| 午夜免费激情av| 一进一出抽搐gif免费好疼 | 丰满的人妻完整版| 亚洲欧美日韩高清在线视频| 长腿黑丝高跟| av电影中文网址| 人成视频在线观看免费观看| 欧美成人午夜精品| 日本一区二区免费在线视频| 无遮挡黄片免费观看| 熟女少妇亚洲综合色aaa.| 日韩欧美一区二区三区在线观看| av国产精品久久久久影院| 一个人免费在线观看的高清视频| 欧美精品啪啪一区二区三区| 久久精品91无色码中文字幕| 国产高清视频在线播放一区| 一级a爱片免费观看的视频| 国产97色在线日韩免费| 午夜福利一区二区在线看| 我的亚洲天堂| 精品久久久久久久毛片微露脸| 99国产精品一区二区三区| 亚洲五月天丁香| 丰满饥渴人妻一区二区三| 操出白浆在线播放| 日本欧美视频一区| 精品第一国产精品| 午夜福利在线免费观看网站| 老司机在亚洲福利影院| 国产熟女午夜一区二区三区| 亚洲欧美日韩无卡精品| av天堂久久9| 美女高潮到喷水免费观看| 国产免费男女视频| 超色免费av| 伊人久久大香线蕉亚洲五| 精品欧美一区二区三区在线| 岛国视频午夜一区免费看| 一级毛片高清免费大全| 男女下面插进去视频免费观看| 夜夜看夜夜爽夜夜摸 | 麻豆久久精品国产亚洲av | 宅男免费午夜| 午夜日韩欧美国产| 99国产精品一区二区蜜桃av| 亚洲精品国产区一区二| 久久久国产成人精品二区 | 亚洲av成人一区二区三| 电影成人av| 不卡一级毛片| 在线av久久热| 亚洲,欧美精品.| 看片在线看免费视频| 久久精品国产亚洲av高清一级| 亚洲中文字幕日韩| 久久这里只有精品19| 午夜a级毛片| 日韩大尺度精品在线看网址 | 中出人妻视频一区二区| 国产成年人精品一区二区 | a级毛片在线看网站| 女人被躁到高潮嗷嗷叫费观| 亚洲欧美精品综合久久99| 91麻豆精品激情在线观看国产 | 人成视频在线观看免费观看| 每晚都被弄得嗷嗷叫到高潮| 欧美一级毛片孕妇| 88av欧美| 欧美人与性动交α欧美软件| 国产精品成人在线| 女人高潮潮喷娇喘18禁视频| 国产精品永久免费网站| 成人亚洲精品一区在线观看| 精品福利观看| 欧美日韩一级在线毛片| 国产精品香港三级国产av潘金莲| 脱女人内裤的视频| 性色av乱码一区二区三区2| 中文亚洲av片在线观看爽| 亚洲人成伊人成综合网2020| 另类亚洲欧美激情| 9191精品国产免费久久| 久久久久国内视频| 国产男靠女视频免费网站| www.熟女人妻精品国产| 久久久久久久午夜电影 | 真人一进一出gif抽搐免费| 亚洲少妇的诱惑av| 久久影院123| 国产精品av久久久久免费| www.自偷自拍.com| 亚洲avbb在线观看| 日韩三级视频一区二区三区| 久久久久国产精品人妻aⅴ院| 国产精品成人在线| 亚洲国产精品999在线| svipshipincom国产片| 国产午夜精品久久久久久| 色老头精品视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利欧美成人| 真人做人爱边吃奶动态| 久久精品91无色码中文字幕| 99国产极品粉嫩在线观看| 精品久久久久久成人av| 久久久国产精品麻豆| 一级a爱视频在线免费观看| 精品国产一区二区久久| 久久精品国产亚洲av香蕉五月| 成人精品一区二区免费| 日韩人妻精品一区2区三区| 纯流量卡能插随身wifi吗| 校园春色视频在线观看| 91精品三级在线观看| 在线十欧美十亚洲十日本专区| 麻豆国产av国片精品| 欧美在线黄色| 国产1区2区3区精品| a级毛片在线看网站| 99久久人妻综合| av天堂久久9| 欧美激情高清一区二区三区| 999精品在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 一二三四社区在线视频社区8| 欧美日韩瑟瑟在线播放| 国产精品美女特级片免费视频播放器 | 亚洲,欧美精品.| 国产免费男女视频| 欧美人与性动交α欧美精品济南到| av在线播放免费不卡| 中文字幕高清在线视频| a在线观看视频网站| 色在线成人网| 精品国内亚洲2022精品成人| 一级片免费观看大全| 一进一出抽搐gif免费好疼 | 亚洲中文av在线| 久久国产精品影院| 久久久久久亚洲精品国产蜜桃av| 岛国视频午夜一区免费看| 欧美精品啪啪一区二区三区| 一进一出抽搐动态| 亚洲成人国产一区在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲色图av天堂| 精品久久久久久电影网| 午夜福利,免费看| 热99国产精品久久久久久7| 99香蕉大伊视频| 国产精品久久久人人做人人爽| 精品久久久精品久久久| 波多野结衣一区麻豆| 国产主播在线观看一区二区| 欧美乱码精品一区二区三区| 后天国语完整版免费观看| 国产主播在线观看一区二区| 咕卡用的链子| 99在线视频只有这里精品首页| 久热这里只有精品99| 香蕉久久夜色| 别揉我奶头~嗯~啊~动态视频| 午夜福利一区二区在线看| 亚洲专区字幕在线| 久99久视频精品免费| 1024香蕉在线观看| 免费日韩欧美在线观看| 国产成人一区二区三区免费视频网站| 精品乱码久久久久久99久播| 欧美性长视频在线观看| 三级毛片av免费| 99精品欧美一区二区三区四区| 色婷婷久久久亚洲欧美| bbb黄色大片| 欧美黄色片欧美黄色片| 手机成人av网站| 桃色一区二区三区在线观看| 老汉色av国产亚洲站长工具| 美女国产高潮福利片在线看| 精品国产一区二区久久| 91精品三级在线观看| 变态另类成人亚洲欧美熟女 |