• <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| 亚洲在线观看片| 国产一级毛片在线| 一本一本综合久久| 水蜜桃什么品种好| 亚洲av免费高清在线观看| 精品人妻熟女av久视频| 久久6这里有精品| 91在线精品国自产拍蜜月| 久久久精品欧美日韩精品| 免费av毛片视频| 久久久成人免费电影| 亚洲精品久久久久久婷婷小说 | 国内精品美女久久久久久| 人人妻人人澡人人爽人人夜夜 | 久久精品久久精品一区二区三区| 精品久久久久久电影网 | 成人亚洲精品av一区二区| 在线天堂最新版资源| 国产精品野战在线观看| 99热这里只有是精品在线观看| 亚洲五月天丁香| 有码 亚洲区| 欧美又色又爽又黄视频| 日韩国内少妇激情av| 久久精品熟女亚洲av麻豆精品 | 成人综合一区亚洲| 欧美+日韩+精品| 亚洲精品aⅴ在线观看| 国产免费一级a男人的天堂| 精品久久久久久成人av| 久久久精品大字幕| 看十八女毛片水多多多| 夫妻性生交免费视频一级片| 亚洲丝袜综合中文字幕| 国产激情偷乱视频一区二区| 国产欧美日韩精品一区二区| 国产高清三级在线| 小蜜桃在线观看免费完整版高清| 少妇熟女aⅴ在线视频| 亚洲成人av在线免费| 国产三级中文精品| 日韩精品青青久久久久久| eeuss影院久久| 国产精品久久久久久精品电影| 我的女老师完整版在线观看| 又粗又硬又长又爽又黄的视频| 狂野欧美激情性xxxx在线观看| 一级毛片电影观看 | 99热这里只有精品一区| 国产色爽女视频免费观看| 人人妻人人看人人澡| 日韩国内少妇激情av| 成人亚洲欧美一区二区av| 欧美一区二区精品小视频在线| 日本与韩国留学比较| 国产一区二区在线av高清观看| 国产片特级美女逼逼视频| 大又大粗又爽又黄少妇毛片口| 免费观看在线日韩| 亚州av有码| 亚洲欧美清纯卡通| 成人美女网站在线观看视频| 中文天堂在线官网| 精品久久久久久久久av| 51国产日韩欧美| 天堂av国产一区二区熟女人妻| 国产女主播在线喷水免费视频网站 | 精品人妻视频免费看| 青春草亚洲视频在线观看| 最近的中文字幕免费完整| 天堂av国产一区二区熟女人妻| 欧美一区二区精品小视频在线| 国产一区二区在线av高清观看| 人妻少妇偷人精品九色| 日韩欧美在线乱码| 亚洲乱码一区二区免费版| 听说在线观看完整版免费高清| 少妇熟女欧美另类| 国产成人一区二区在线| 婷婷色综合大香蕉| 美女内射精品一级片tv| 老女人水多毛片| 国产熟女欧美一区二区| 亚洲自偷自拍三级| 国产精品美女特级片免费视频播放器| 国产成人午夜福利电影在线观看| 国产三级在线视频| 日日干狠狠操夜夜爽| 亚洲人成网站在线播| 国产白丝娇喘喷水9色精品| 成人国产麻豆网| 性插视频无遮挡在线免费观看| 亚洲美女搞黄在线观看| 最近中文字幕2019免费版| 久久久久久伊人网av| 亚洲综合精品二区| 大又大粗又爽又黄少妇毛片口| 成人毛片a级毛片在线播放| 亚洲内射少妇av| 51国产日韩欧美| 中国美白少妇内射xxxbb| 亚洲精品久久久久久婷婷小说 | 国产淫语在线视频| 久久久午夜欧美精品| 欧美性猛交╳xxx乱大交人| 亚洲一区高清亚洲精品| 国产精品一区二区三区四区免费观看| 免费av毛片视频| 免费观看的影片在线观看| 免费观看精品视频网站| 男人舔奶头视频| 欧美性感艳星| 成年av动漫网址| 人体艺术视频欧美日本| 国产黄a三级三级三级人| 国产中年淑女户外野战色| 观看免费一级毛片| 日韩人妻高清精品专区| 嫩草影院新地址| 欧美另类亚洲清纯唯美| 全区人妻精品视频| 色哟哟·www| 少妇人妻精品综合一区二区| 不卡视频在线观看欧美| 久久久久精品久久久久真实原创| 久久久久九九精品影院| 成人欧美大片| 日本与韩国留学比较| 岛国在线免费视频观看| 最近最新中文字幕大全电影3| 中文亚洲av片在线观看爽| 中文精品一卡2卡3卡4更新| 亚洲欧美中文字幕日韩二区| 欧美日韩国产亚洲二区| 日本wwww免费看| 2021少妇久久久久久久久久久| 亚洲成人中文字幕在线播放| av卡一久久| 中文字幕人妻熟人妻熟丝袜美| av.在线天堂| 精品午夜福利在线看| 亚洲国产精品成人久久小说| 特级一级黄色大片| 久久久久久久久大av| 国产片特级美女逼逼视频| 国产精品人妻久久久久久| 欧美日韩在线观看h| 毛片一级片免费看久久久久| 亚洲在线自拍视频| 蜜臀久久99精品久久宅男| 欧美激情在线99| 国产大屁股一区二区在线视频| 国产免费福利视频在线观看| 亚洲av电影在线观看一区二区三区 | 精品午夜福利在线看| 久久精品夜色国产| 国产精品一区二区三区四区免费观看| 中文精品一卡2卡3卡4更新| 非洲黑人性xxxx精品又粗又长| 亚洲欧美日韩东京热| 99久久精品热视频| 好男人在线观看高清免费视频| 日韩欧美在线乱码| 99久久精品热视频| 在线免费观看不下载黄p国产| 伊人久久精品亚洲午夜| 国产精品一及| 非洲黑人性xxxx精品又粗又长| www.av在线官网国产| 国产私拍福利视频在线观看| 少妇人妻精品综合一区二区| 九草在线视频观看| 国产成人精品婷婷| 亚洲精品456在线播放app| 精品国产一区二区三区久久久樱花 | 在现免费观看毛片| 国产精品不卡视频一区二区| 国产伦精品一区二区三区视频9| 亚洲最大成人手机在线| 熟妇人妻久久中文字幕3abv| 日日撸夜夜添| 精品无人区乱码1区二区| www.av在线官网国产| 国产精品,欧美在线| 51国产日韩欧美| 大香蕉久久网| 国产成人91sexporn| 中文乱码字字幕精品一区二区三区 | 国产午夜精品论理片| 草草在线视频免费看| 成人av在线播放网站| 搞女人的毛片| 国产男人的电影天堂91| 国产精品一二三区在线看| 国产精品.久久久| 国产成人精品婷婷| 欧美极品一区二区三区四区| 我的老师免费观看完整版| 观看美女的网站| 级片在线观看| 村上凉子中文字幕在线| 久久久久九九精品影院| 内射极品少妇av片p| 丰满人妻一区二区三区视频av| 最近的中文字幕免费完整| 女人久久www免费人成看片 | 国产精品.久久久| 久久精品影院6| 嘟嘟电影网在线观看| 日本免费在线观看一区| 青青草视频在线视频观看| 级片在线观看| 欧美日本亚洲视频在线播放| 国产亚洲av嫩草精品影院| 国产伦在线观看视频一区| 国产私拍福利视频在线观看| 久久久国产成人精品二区| 欧美日韩一区二区视频在线观看视频在线 | 不卡视频在线观看欧美| 日本猛色少妇xxxxx猛交久久| 中文字幕制服av| 国产免费一级a男人的天堂| av在线播放精品| 啦啦啦韩国在线观看视频| 精华霜和精华液先用哪个| 好男人在线观看高清免费视频| 亚洲欧美日韩东京热| 天堂√8在线中文| 亚洲av电影不卡..在线观看| 欧美精品一区二区大全| 色播亚洲综合网| 精品久久久久久久人妻蜜臀av| 日本wwww免费看| 国产中年淑女户外野战色| 色网站视频免费| 欧美性猛交╳xxx乱大交人| 久久久a久久爽久久v久久| 日本与韩国留学比较| 99久国产av精品国产电影| 搞女人的毛片| 国产精品日韩av在线免费观看| 九九爱精品视频在线观看| 午夜日本视频在线| 国产亚洲精品久久久com| 国产精品乱码一区二三区的特点| 成年女人永久免费观看视频| videos熟女内射| 中文欧美无线码| 1024手机看黄色片| 欧美最新免费一区二区三区| 免费看美女性在线毛片视频| 国产精品电影一区二区三区| 亚洲欧美成人综合另类久久久 | 亚洲欧美成人综合另类久久久 | 婷婷六月久久综合丁香| 久久精品国产99精品国产亚洲性色| av在线蜜桃| 99国产精品一区二区蜜桃av| 波多野结衣巨乳人妻| 免费av毛片视频| 超碰97精品在线观看| 三级国产精品欧美在线观看| 18禁在线无遮挡免费观看视频| 亚洲国产日韩欧美精品在线观看| 国产av不卡久久| 亚洲图色成人| 哪个播放器可以免费观看大片| 国产精品一区www在线观看| 在线a可以看的网站| 青春草亚洲视频在线观看| 国产成人a区在线观看| 3wmmmm亚洲av在线观看| 麻豆精品久久久久久蜜桃| 精品熟女少妇av免费看| 亚洲高清免费不卡视频| 久久久久国产网址| 18禁在线播放成人免费| 99久久精品一区二区三区| 亚洲伊人久久精品综合 | 91久久精品电影网| 男人的好看免费观看在线视频| 两性午夜刺激爽爽歪歪视频在线观看| 欧美一区二区国产精品久久精品| 国产亚洲av片在线观看秒播厂 | 国产麻豆成人av免费视频| 最近中文字幕2019免费版| 免费观看性生交大片5| 看片在线看免费视频| 美女黄网站色视频| 亚洲图色成人| 久久精品夜夜夜夜夜久久蜜豆| 建设人人有责人人尽责人人享有的 | 黑人高潮一二区| 免费黄网站久久成人精品| av播播在线观看一区| 久久久色成人| 亚洲第一区二区三区不卡| 丝袜美腿在线中文| 日韩国内少妇激情av| 国产精品爽爽va在线观看网站| 国产精品麻豆人妻色哟哟久久 | 国产成人a∨麻豆精品| 人人妻人人看人人澡| 国产男人的电影天堂91| 国产精品麻豆人妻色哟哟久久 | 18禁在线无遮挡免费观看视频| 国产免费男女视频| 一级二级三级毛片免费看| 美女被艹到高潮喷水动态| 在线免费观看的www视频| 老司机福利观看| 国产三级在线视频| 桃色一区二区三区在线观看| 国产三级中文精品| 国产麻豆成人av免费视频| 国产免费福利视频在线观看| 99久久精品国产国产毛片| .国产精品久久| av视频在线观看入口| 国产午夜精品论理片| 久久久精品94久久精品| 国产欧美日韩精品一区二区| 精品久久久久久久久av| av专区在线播放| 亚洲av福利一区| 中文字幕人妻熟人妻熟丝袜美| 免费播放大片免费观看视频在线观看 | 亚洲精华国产精华液的使用体验| 天堂√8在线中文| 十八禁国产超污无遮挡网站| 亚洲av免费在线观看| 精品久久久久久久久久久久久| 欧美bdsm另类| 国产精品久久久久久久久免| av国产免费在线观看| 色视频www国产| 午夜日本视频在线| 人妻系列 视频| 国产乱人视频| 精品无人区乱码1区二区| 天天躁夜夜躁狠狠久久av| 日韩亚洲欧美综合| 国产极品精品免费视频能看的| 一区二区三区乱码不卡18| 一个人看的www免费观看视频| 日日干狠狠操夜夜爽| 搞女人的毛片| 亚洲av熟女| 一个人看视频在线观看www免费| 乱人视频在线观看| 麻豆成人av视频| 18禁在线播放成人免费| 熟妇人妻久久中文字幕3abv| 美女xxoo啪啪120秒动态图| 亚洲成av人片在线播放无| a级毛片免费高清观看在线播放| 免费人成在线观看视频色| 亚洲自偷自拍三级| 一级毛片aaaaaa免费看小| 亚洲国产欧美在线一区| 可以在线观看毛片的网站| 国产男人的电影天堂91| 秋霞伦理黄片| 男女边吃奶边做爰视频| 人人妻人人看人人澡| 久久久久久大精品| 日日摸夜夜添夜夜添av毛片| 一级毛片久久久久久久久女| 两性午夜刺激爽爽歪歪视频在线观看| 欧美激情在线99| 国产精品国产高清国产av| 青青草视频在线视频观看| 麻豆成人av视频| 久久综合国产亚洲精品| 天天一区二区日本电影三级| 免费搜索国产男女视频| 最近中文字幕2019免费版| 插阴视频在线观看视频| 日本免费在线观看一区| 禁无遮挡网站| 91久久精品电影网| 免费在线观看成人毛片| 亚洲国产精品sss在线观看| 国产一级毛片在线| 亚洲欧美成人综合另类久久久 | 精品久久久久久久久av| 看黄色毛片网站| 久久久精品欧美日韩精品| 国产伦理片在线播放av一区| 日日干狠狠操夜夜爽| 大香蕉久久网| 日韩高清综合在线| 国内精品美女久久久久久| 国产片特级美女逼逼视频| 中文天堂在线官网| 三级国产精品片| 黄片无遮挡物在线观看| 午夜精品国产一区二区电影 | 色综合色国产| 性色avwww在线观看| 亚洲中文字幕日韩| 午夜爱爱视频在线播放| av在线观看视频网站免费| 成年免费大片在线观看| av免费观看日本| 国产91av在线免费观看| 亚洲伊人久久精品综合 | 久久精品夜夜夜夜夜久久蜜豆| 亚洲美女搞黄在线观看| 97超视频在线观看视频| 亚洲丝袜综合中文字幕| 亚洲av电影在线观看一区二区三区 | 色综合亚洲欧美另类图片| 久久99热6这里只有精品| 国产免费一级a男人的天堂| 国产成人a∨麻豆精品| 成年版毛片免费区| 一本久久精品| 国产精品综合久久久久久久免费| a级毛色黄片| 成人三级黄色视频| 亚洲图色成人| 亚洲av福利一区| 午夜视频国产福利| 搡老妇女老女人老熟妇| 国产精品1区2区在线观看.| 亚洲国产最新在线播放| av在线蜜桃| 少妇人妻精品综合一区二区| 久久国产乱子免费精品| 在线播放国产精品三级| 亚洲国产欧美在线一区| 99热精品在线国产| 国产亚洲最大av| 欧美日韩国产亚洲二区| 99久久无色码亚洲精品果冻| 日本-黄色视频高清免费观看| 简卡轻食公司| 91aial.com中文字幕在线观看| 午夜福利成人在线免费观看| 国产精品国产三级国产av玫瑰| 99热精品在线国产| 国产精品女同一区二区软件| 国产黄片美女视频| 老司机影院成人| 亚洲精品日韩在线中文字幕| 18禁动态无遮挡网站| 精品国产三级普通话版| 天堂网av新在线| 成人特级av手机在线观看| 日日摸夜夜添夜夜添av毛片| 久久韩国三级中文字幕| 国产极品精品免费视频能看的| 国产美女午夜福利| 日本色播在线视频| 国产 一区 欧美 日韩| 精品久久国产蜜桃| 一个人看视频在线观看www免费| 天天躁日日操中文字幕| 两个人的视频大全免费| 国产精品国产三级国产专区5o | av女优亚洲男人天堂| 欧美日韩综合久久久久久| 亚洲五月天丁香| 国语自产精品视频在线第100页| 一级毛片久久久久久久久女| 97超碰精品成人国产| 欧美一区二区国产精品久久精品| 精品国内亚洲2022精品成人| 国产亚洲91精品色在线| 亚洲在线自拍视频| 国产精品,欧美在线| 久久精品国产鲁丝片午夜精品| 草草在线视频免费看| 国产日韩欧美在线精品| 狂野欧美激情性xxxx在线观看| 看免费成人av毛片| 毛片女人毛片| 人人妻人人澡欧美一区二区| 国产一级毛片在线| 丰满人妻一区二区三区视频av| 看片在线看免费视频| 最近中文字幕高清免费大全6| 最近中文字幕2019免费版| 久99久视频精品免费| 91aial.com中文字幕在线观看| 99久久人妻综合| 亚洲av成人av| 欧美人与善性xxx| 亚洲高清免费不卡视频| 久久久久久九九精品二区国产| 女人被狂操c到高潮| 哪个播放器可以免费观看大片| 日韩欧美国产在线观看| 久久久精品94久久精品| 国产精品.久久久| 1024手机看黄色片| 亚洲精品日韩av片在线观看| 又黄又爽又刺激的免费视频.| 国国产精品蜜臀av免费| www.av在线官网国产| 午夜精品在线福利| 亚洲乱码一区二区免费版| 直男gayav资源| 国产亚洲午夜精品一区二区久久 | 亚洲三级黄色毛片| 日韩成人av中文字幕在线观看| 亚洲经典国产精华液单| 久99久视频精品免费| www日本黄色视频网| 国产精品麻豆人妻色哟哟久久 | 午夜亚洲福利在线播放| 欧美日韩一区二区视频在线观看视频在线 | av播播在线观看一区| 尤物成人国产欧美一区二区三区| 国产视频首页在线观看| av福利片在线观看| 国产中年淑女户外野战色| 热99在线观看视频| 国产激情偷乱视频一区二区| 99久久无色码亚洲精品果冻| 日韩人妻高清精品专区| 国产91av在线免费观看| 麻豆国产97在线/欧美| 两性午夜刺激爽爽歪歪视频在线观看| 大香蕉97超碰在线| 亚洲国产精品sss在线观看| 国产高清三级在线| 我要搜黄色片| 国产午夜福利久久久久久| 韩国高清视频一区二区三区| 久久婷婷人人爽人人干人人爱| 99久久精品一区二区三区| 色播亚洲综合网| 日本免费一区二区三区高清不卡| 中文字幕亚洲精品专区| 高清av免费在线| 精品少妇黑人巨大在线播放 | 国产精品国产三级专区第一集| 亚洲精品乱码久久久v下载方式| 成年女人永久免费观看视频| h日本视频在线播放| 国产美女午夜福利| 久久婷婷人人爽人人干人人爱| 赤兔流量卡办理| 国产精品久久视频播放| 国产三级中文精品| 日韩欧美国产在线观看| 欧美一区二区亚洲| 高清午夜精品一区二区三区| 免费看av在线观看网站| 国产精品人妻久久久影院| 国产探花极品一区二区| 亚洲国产精品成人久久小说| 久久精品夜夜夜夜夜久久蜜豆| 麻豆一二三区av精品| 午夜福利在线观看吧| 日产精品乱码卡一卡2卡三| 国产综合懂色| 搡老妇女老女人老熟妇| 91久久精品电影网| 亚洲不卡免费看| 91久久精品国产一区二区三区| av福利片在线观看| 国产亚洲一区二区精品| 亚洲18禁久久av| 亚洲国产精品合色在线| 最近2019中文字幕mv第一页| 亚洲国产精品专区欧美| av在线亚洲专区| 久久久色成人| 午夜a级毛片| 日韩欧美在线乱码| 免费观看的影片在线观看| av国产免费在线观看| 99久久精品一区二区三区| 一区二区三区免费毛片| 国产亚洲5aaaaa淫片| 十八禁国产超污无遮挡网站| 最近中文字幕高清免费大全6| 欧美精品一区二区大全| 国产在线一区二区三区精 | 久久这里有精品视频免费| 18禁在线无遮挡免费观看视频| 亚洲精品,欧美精品| 久久久精品大字幕| 少妇猛男粗大的猛烈进出视频 | 边亲边吃奶的免费视频| 亚洲不卡免费看| 18禁在线播放成人免费| 人妻夜夜爽99麻豆av| 内地一区二区视频在线| 乱人视频在线观看| 国产精品一区二区性色av| 亚洲美女搞黄在线观看| 99九九线精品视频在线观看视频| 身体一侧抽搐| 高清毛片免费看| 美女xxoo啪啪120秒动态图| 我的老师免费观看完整版| av线在线观看网站| 偷拍熟女少妇极品色| 边亲边吃奶的免费视频| 99久久成人亚洲精品观看| 国产黄色小视频在线观看| 成人毛片a级毛片在线播放| 精品久久久久久久久久久久久| 日本av手机在线免费观看| 日韩一本色道免费dvd| 国产高清三级在线| 一个人观看的视频www高清免费观看| 蜜臀久久99精品久久宅男|