• <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中文字幕在线观看| 亚洲精品日本国产第一区| 国产爽快片一区二区三区| 97超碰精品成人国产| 男女下面插进去视频免费观看 | 日韩一区二区三区影片| 97超碰精品成人国产| 伊人亚洲综合成人网| 免费观看a级毛片全部| 国产亚洲最大av| av电影中文网址| 久久人人97超碰香蕉20202| 久久久久久久久久成人| 婷婷色综合大香蕉| 国产1区2区3区精品| 激情视频va一区二区三区| 久热这里只有精品99| 满18在线观看网站| 亚洲综合色网址| 最近手机中文字幕大全| 欧美97在线视频| 插逼视频在线观看| 国产一区亚洲一区在线观看| 亚洲国产色片| 欧美变态另类bdsm刘玥| 汤姆久久久久久久影院中文字幕| 国产视频首页在线观看| 欧美3d第一页| 在线天堂中文资源库| 成人手机av| 青春草亚洲视频在线观看| 国产黄色视频一区二区在线观看| kizo精华| av播播在线观看一区| 男女国产视频网站| 久久久精品区二区三区| 看免费成人av毛片| 国产精品秋霞免费鲁丝片| 黄色毛片三级朝国网站| 中文字幕制服av| 不卡视频在线观看欧美| 久久精品久久久久久噜噜老黄| 天美传媒精品一区二区| 中文字幕制服av| 看免费成人av毛片| 99热网站在线观看| 丰满少妇做爰视频| 男女午夜视频在线观看 | 国产精品一区www在线观看| 亚洲第一av免费看| a级毛片在线看网站| 97在线视频观看| 午夜激情av网站| 丁香六月天网| 亚洲av日韩在线播放| 亚洲成国产人片在线观看| 午夜老司机福利剧场| 精品亚洲成a人片在线观看| 久久国产精品大桥未久av| 亚洲国产av影院在线观看| 国产成人免费无遮挡视频| 男人舔女人的私密视频| 国产有黄有色有爽视频| 少妇的丰满在线观看| 高清欧美精品videossex| 一个人免费看片子| 男女免费视频国产| 伦理电影免费视频| 国产精品久久久久久久久免| 久久鲁丝午夜福利片| 最黄视频免费看| 日日爽夜夜爽网站| 国产色爽女视频免费观看| 全区人妻精品视频| 亚洲综合色网址| 黄色视频在线播放观看不卡| 亚洲精华国产精华液的使用体验| 2022亚洲国产成人精品| 观看美女的网站| 亚洲精品中文字幕在线视频| 亚洲经典国产精华液单| 在线观看三级黄色| 黄网站色视频无遮挡免费观看| kizo精华| 成人国产麻豆网| 精品国产国语对白av| 久久久a久久爽久久v久久| 久久女婷五月综合色啪小说| 啦啦啦在线观看免费高清www| 国产精品一区二区在线不卡| 男女边吃奶边做爰视频| 新久久久久国产一级毛片| 亚洲美女搞黄在线观看| av在线老鸭窝| 国产免费又黄又爽又色| 宅男免费午夜| 久久国产精品大桥未久av| 国产男女超爽视频在线观看| 五月开心婷婷网| 91精品三级在线观看| 亚洲成色77777| 成人亚洲欧美一区二区av| av.在线天堂| 最新中文字幕久久久久| 啦啦啦在线观看免费高清www| 91在线精品国自产拍蜜月| 久久人人爽av亚洲精品天堂| 亚洲情色 制服丝袜| 在线观看美女被高潮喷水网站| 国产免费视频播放在线视频| 高清欧美精品videossex| 亚洲精品中文字幕在线视频| 精品国产一区二区三区久久久樱花| 成人国语在线视频| a级毛色黄片| 精品少妇黑人巨大在线播放| 久久久久视频综合| 国产精品一二三区在线看| 亚洲在久久综合| 日韩精品有码人妻一区| 男人爽女人下面视频在线观看| 两个人免费观看高清视频| 欧美日韩成人在线一区二区| 一区二区三区四区激情视频| 亚洲第一av免费看| 最近中文字幕2019免费版| 99国产综合亚洲精品| 亚洲美女视频黄频| 日韩 亚洲 欧美在线| 午夜视频国产福利| 岛国毛片在线播放| 日韩成人伦理影院| a级片在线免费高清观看视频| 亚洲成av片中文字幕在线观看 | 青春草国产在线视频| 国产日韩欧美在线精品| 日本av免费视频播放| 91精品国产国语对白视频| 国产日韩欧美亚洲二区| 男女啪啪激烈高潮av片| 一级片'在线观看视频| 91在线精品国自产拍蜜月| 中文乱码字字幕精品一区二区三区| 亚洲欧美日韩另类电影网站| 热re99久久精品国产66热6| a级毛色黄片| 国产色爽女视频免费观看| 久久精品国产亚洲av涩爱| 欧美精品一区二区大全| 2018国产大陆天天弄谢| 纯流量卡能插随身wifi吗| 国产av码专区亚洲av| 男女边吃奶边做爰视频| 精品亚洲乱码少妇综合久久| 欧美日韩一区二区视频在线观看视频在线| 午夜激情久久久久久久| 亚洲国产精品一区三区| 深夜精品福利| 91精品三级在线观看| 黄色配什么色好看| a级毛片在线看网站| 久久精品久久久久久噜噜老黄| 免费av中文字幕在线| 夜夜骑夜夜射夜夜干| 男人添女人高潮全过程视频| 中文字幕人妻熟女乱码| 黄色视频在线播放观看不卡| 精品人妻熟女毛片av久久网站| 美女国产视频在线观看| 国产69精品久久久久777片| 国产精品一区二区在线观看99| 国产成人精品婷婷| 高清毛片免费看| 国产成人精品无人区| 国产亚洲最大av| 久久久久久人妻| 国产免费福利视频在线观看| 青青草视频在线视频观看| 精品亚洲乱码少妇综合久久| 国产黄频视频在线观看| a级毛片黄视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日本-黄色视频高清免费观看| 91久久精品国产一区二区三区| 香蕉精品网在线| 中文字幕另类日韩欧美亚洲嫩草| 伦理电影免费视频| 精品亚洲成a人片在线观看| 午夜日本视频在线| 99久久精品国产国产毛片| 亚洲欧美成人精品一区二区| 五月玫瑰六月丁香| 黑人巨大精品欧美一区二区蜜桃 | 香蕉精品网在线| 午夜福利在线观看免费完整高清在| 搡老乐熟女国产| 十八禁高潮呻吟视频| 一级,二级,三级黄色视频| 91精品伊人久久大香线蕉| 午夜福利乱码中文字幕| 精品一区二区三区四区五区乱码 | 狠狠婷婷综合久久久久久88av| av一本久久久久| 国产深夜福利视频在线观看| 91午夜精品亚洲一区二区三区| 免费观看无遮挡的男女| 亚洲,欧美精品.| 国产成人91sexporn| 99热6这里只有精品| 久久女婷五月综合色啪小说| 精品第一国产精品| 亚洲精品久久成人aⅴ小说| 欧美日韩精品成人综合77777| 久久久久网色| 捣出白浆h1v1| 亚洲内射少妇av| 中文字幕制服av| 亚洲精品aⅴ在线观看| 国精品久久久久久国模美| 美女中出高潮动态图| 26uuu在线亚洲综合色| 亚洲国产色片| 一级毛片黄色毛片免费观看视频| 一二三四中文在线观看免费高清| 亚洲国产精品999| 18禁裸乳无遮挡动漫免费视频| 青春草国产在线视频| av.在线天堂| 校园人妻丝袜中文字幕| kizo精华| 多毛熟女@视频| 日日啪夜夜爽| 哪个播放器可以免费观看大片| 成年人午夜在线观看视频| 精品酒店卫生间| 激情五月婷婷亚洲| 国产成人精品一,二区| 中文字幕人妻熟女乱码| 亚洲人与动物交配视频| 男女免费视频国产| av电影中文网址| 99re6热这里在线精品视频| 国产精品人妻久久久影院| 丝袜在线中文字幕| 亚洲精品久久久久久婷婷小说| 又黄又粗又硬又大视频| 成人亚洲精品一区在线观看| www.熟女人妻精品国产 | 中文字幕最新亚洲高清| 日韩精品免费视频一区二区三区 | 蜜桃国产av成人99| 丝袜人妻中文字幕| 欧美人与性动交α欧美软件 | av不卡在线播放| 亚洲国产日韩一区二区| 美女xxoo啪啪120秒动态图| 午夜影院在线不卡| 少妇人妻久久综合中文| 国产一级毛片在线| 欧美日韩成人在线一区二区| 国产黄色视频一区二区在线观看| 免费观看在线日韩| 啦啦啦视频在线资源免费观看| 免费看不卡的av| 亚洲一级一片aⅴ在线观看| 香蕉精品网在线| 五月开心婷婷网| 亚洲国产欧美在线一区| 午夜福利视频在线观看免费| 亚洲伊人久久精品综合| 性色av一级| 边亲边吃奶的免费视频| 热99国产精品久久久久久7| 岛国毛片在线播放| 看免费av毛片| 色5月婷婷丁香| 一级片'在线观看视频| 国产乱人偷精品视频| 中文欧美无线码| 欧美日韩av久久| 精品国产一区二区三区久久久樱花| xxxhd国产人妻xxx| 欧美人与善性xxx| 一级毛片电影观看| 日韩精品免费视频一区二区三区 | 欧美日韩亚洲高清精品| 两性夫妻黄色片 | 国产精品久久久久久久电影| 大片免费播放器 马上看| 99re6热这里在线精品视频| 欧美国产精品一级二级三级| 亚洲丝袜综合中文字幕| 午夜激情久久久久久久| 男人添女人高潮全过程视频| 亚洲一级一片aⅴ在线观看| 国产爽快片一区二区三区| 国产黄频视频在线观看| 久久99热6这里只有精品| 一级片免费观看大全| 国内精品宾馆在线| 你懂的网址亚洲精品在线观看| 国产av国产精品国产| a级片在线免费高清观看视频| 五月玫瑰六月丁香| 免费少妇av软件| 亚洲国产av新网站| 国产成人欧美| 日韩中字成人| 久久久久久人人人人人| 日韩中字成人| 久久久久久久国产电影| 免费观看av网站的网址| 日本猛色少妇xxxxx猛交久久| 欧美 亚洲 国产 日韩一| 丁香六月天网| 一级片免费观看大全| 黄色怎么调成土黄色| 亚洲伊人色综图| 人成视频在线观看免费观看| 天美传媒精品一区二区| av黄色大香蕉| 亚洲av国产av综合av卡| 91精品伊人久久大香线蕉| 男人添女人高潮全过程视频| 香蕉国产在线看| 18禁在线无遮挡免费观看视频| 男女无遮挡免费网站观看| 亚洲熟女精品中文字幕| 日韩欧美精品免费久久| 国产淫语在线视频| 美女主播在线视频| 国产精品嫩草影院av在线观看| 亚洲久久久国产精品| 一本色道久久久久久精品综合| 亚洲精品成人av观看孕妇| 黄色一级大片看看| 国产视频首页在线观看| 欧美激情 高清一区二区三区| 国产片内射在线| 国产av一区二区精品久久| 2018国产大陆天天弄谢| 午夜免费观看性视频| 一级片'在线观看视频| 天天影视国产精品| 高清av免费在线| 国产色婷婷99| 国产亚洲午夜精品一区二区久久| 另类亚洲欧美激情| 丰满迷人的少妇在线观看| 18+在线观看网站| √禁漫天堂资源中文www| 美女中出高潮动态图| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产毛片av蜜桃av| 国产欧美另类精品又又久久亚洲欧美| 日韩视频在线欧美| 欧美国产精品一级二级三级| 熟女av电影| 一级a做视频免费观看| 国产老妇伦熟女老妇高清| 男女国产视频网站| 欧美国产精品va在线观看不卡| 老女人水多毛片| 在线观看免费日韩欧美大片| 一区二区三区乱码不卡18| 亚洲精品一二三| 日韩制服骚丝袜av| 日韩一区二区三区影片| 国产av码专区亚洲av| 综合色丁香网| 精品一区二区三卡| 少妇被粗大猛烈的视频| 国产精品一区www在线观看| 久久久久国产网址| 国产精品三级大全| 亚洲精品一区蜜桃| 国产精品国产三级专区第一集| 大香蕉97超碰在线| 国产一区二区激情短视频 | 亚洲精品乱久久久久久| 日本欧美视频一区| 欧美国产精品一级二级三级| 少妇的逼水好多| www.熟女人妻精品国产 | 亚洲高清免费不卡视频| 国产精品久久久久久久电影| 亚洲国产色片| 免费高清在线观看视频在线观看| 欧美丝袜亚洲另类| 日韩精品有码人妻一区| 久久久国产欧美日韩av| 国产激情久久老熟女| 国产av一区二区精品久久| 欧美精品av麻豆av| 在线精品无人区一区二区三| 国产成人a∨麻豆精品| 亚洲成人av在线免费| 夜夜骑夜夜射夜夜干| 男女免费视频国产| 黄色视频在线播放观看不卡| 在线观看免费日韩欧美大片| 在线亚洲精品国产二区图片欧美| 久久久久久伊人网av| 国产精品麻豆人妻色哟哟久久| 亚洲综合色网址| 日韩一区二区视频免费看| 伦理电影免费视频| 国产探花极品一区二区| 18禁在线无遮挡免费观看视频| 国产av码专区亚洲av| 老司机亚洲免费影院| 男的添女的下面高潮视频| 国产精品99久久99久久久不卡 | 纵有疾风起免费观看全集完整版| 亚洲精品第二区| 国产精品无大码| 18禁观看日本| 亚洲一级一片aⅴ在线观看| 国产xxxxx性猛交| √禁漫天堂资源中文www| 亚洲精品一二三| 日本爱情动作片www.在线观看| 五月开心婷婷网| 最后的刺客免费高清国语| 免费女性裸体啪啪无遮挡网站| 国产精品久久久久久精品电影小说| 亚洲精品视频女| 热re99久久精品国产66热6| 少妇被粗大猛烈的视频| 亚洲少妇的诱惑av| 九色成人免费人妻av| 久久久欧美国产精品| 黑人欧美特级aaaaaa片| 90打野战视频偷拍视频| 2021少妇久久久久久久久久久| 久久久久久久久久久久大奶| 欧美3d第一页| 少妇人妻久久综合中文| 巨乳人妻的诱惑在线观看| www.熟女人妻精品国产 | 国产在视频线精品| 欧美精品一区二区大全| 亚洲av综合色区一区| 亚洲欧洲国产日韩| 国产精品久久久av美女十八| 中文字幕最新亚洲高清| 国产成人精品福利久久| 九色成人免费人妻av| 热99久久久久精品小说推荐| 日本爱情动作片www.在线观看| 免费av不卡在线播放| 大香蕉97超碰在线| 日产精品乱码卡一卡2卡三| 女的被弄到高潮叫床怎么办| 各种免费的搞黄视频| 51国产日韩欧美| 亚洲av中文av极速乱| 久热久热在线精品观看| 成人国语在线视频| 国产成人午夜福利电影在线观看| 免费观看无遮挡的男女| 制服丝袜香蕉在线| 黄网站色视频无遮挡免费观看| 精品国产乱码久久久久久小说| 日本91视频免费播放| 男女边吃奶边做爰视频| 成年人免费黄色播放视频| av电影中文网址| 日韩av不卡免费在线播放| 欧美3d第一页| 成人国产av品久久久| 国产日韩欧美视频二区| 久久精品久久久久久噜噜老黄| 一个人免费看片子| 制服人妻中文乱码| 成人亚洲精品一区在线观看| 欧美日韩av久久| 欧美日韩视频精品一区| 丝袜喷水一区| 精品第一国产精品| 色网站视频免费| 熟女人妻精品中文字幕| 久久精品国产自在天天线| 在线观看人妻少妇| 在线观看国产h片| 成人免费观看视频高清| 曰老女人黄片| 91午夜精品亚洲一区二区三区| 亚洲人与动物交配视频| 伊人亚洲综合成人网| 一级,二级,三级黄色视频| 国产免费视频播放在线视频| 精品人妻熟女毛片av久久网站| 成人18禁高潮啪啪吃奶动态图| 久久久久久久亚洲中文字幕| 亚洲综合色惰| 久久精品国产a三级三级三级| 视频区图区小说| 飞空精品影院首页| 午夜av观看不卡| 婷婷成人精品国产| 波多野结衣一区麻豆| 亚洲国产欧美日韩在线播放| 一区在线观看完整版| 婷婷成人精品国产| 日本色播在线视频| 欧美日韩亚洲高清精品| 国产一区二区在线观看日韩| 亚洲国产成人一精品久久久| av免费观看日本| 纯流量卡能插随身wifi吗| 日韩伦理黄色片| 99久久人妻综合| 亚洲av男天堂| 久久久a久久爽久久v久久| 精品久久蜜臀av无| 精品视频人人做人人爽| 有码 亚洲区| 亚洲精华国产精华液的使用体验| 9191精品国产免费久久| freevideosex欧美| 新久久久久国产一级毛片| 搡老乐熟女国产| 亚洲成人av在线免费| 少妇熟女欧美另类| av免费在线看不卡| 国产毛片在线视频| 亚洲欧美日韩卡通动漫| 女人久久www免费人成看片| 男女边吃奶边做爰视频| 日韩免费高清中文字幕av| 性高湖久久久久久久久免费观看| 日本爱情动作片www.在线观看| 人妻一区二区av| xxxhd国产人妻xxx| kizo精华| 国产综合精华液| 亚洲国产日韩一区二区| 91精品国产国语对白视频| 久久女婷五月综合色啪小说| videossex国产| 你懂的网址亚洲精品在线观看| av不卡在线播放| 美女大奶头黄色视频| 欧美成人午夜精品| 国产精品无大码| 精品99又大又爽又粗少妇毛片| 国产成人91sexporn| 亚洲欧美色中文字幕在线| 一级a做视频免费观看| 亚洲精华国产精华液的使用体验| 久久人人爽av亚洲精品天堂| 黄色配什么色好看| 99九九在线精品视频| 免费女性裸体啪啪无遮挡网站| av女优亚洲男人天堂| a级毛片在线看网站| 天堂俺去俺来也www色官网| 制服诱惑二区| 男人舔女人的私密视频| 日日撸夜夜添| 精品视频人人做人人爽| 天堂8中文在线网| 成人亚洲欧美一区二区av| 国内精品宾馆在线| 1024视频免费在线观看| 婷婷成人精品国产| 国产精品秋霞免费鲁丝片| 国产高清三级在线| 亚洲av国产av综合av卡| 中文乱码字字幕精品一区二区三区| 久久免费观看电影| 国产精品久久久久久久久免| www.色视频.com| 日韩一区二区视频免费看| 国产综合精华液| 老司机影院毛片| 欧美日韩av久久| 久久99精品国语久久久| 女人被躁到高潮嗷嗷叫费观| 七月丁香在线播放| 91成人精品电影| 女人精品久久久久毛片| 精品一区二区三卡| 国产精品三级大全| 午夜免费观看性视频| 综合色丁香网| 国产高清不卡午夜福利| 免费日韩欧美在线观看| 丝袜人妻中文字幕| 黄色 视频免费看| 精品久久久久久电影网| 亚洲国产毛片av蜜桃av| 男女高潮啪啪啪动态图| 久久精品国产亚洲av涩爱| 好男人视频免费观看在线| 午夜福利网站1000一区二区三区| 99热这里只有是精品在线观看| 久久精品久久久久久噜噜老黄| 久久久a久久爽久久v久久| 免费av中文字幕在线| 精品人妻熟女毛片av久久网站| 王馨瑶露胸无遮挡在线观看| 午夜福利网站1000一区二区三区| 国产69精品久久久久777片| 精品一区二区三区四区五区乱码 | 婷婷成人精品国产| 亚洲一级一片aⅴ在线观看| 成年人免费黄色播放视频| 秋霞在线观看毛片| 自拍欧美九色日韩亚洲蝌蚪91| 免费av不卡在线播放| 国产一级毛片在线|