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

    羌塘高寒草地物種多樣性與生態(tài)系統(tǒng)多功能關(guān)系格局

    2016-08-09 01:22:02熊定鵬趙廣帥武建雙石培禮張憲洲
    生態(tài)學報 2016年11期
    關(guān)鍵詞:物種多樣性

    熊定鵬,趙廣帥,武建雙,石培禮,*,張憲洲

    1 中國科學院地理科學與資源研究所生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測與模擬重點實驗室, 北京 100101 2 中國科學院大學, 北京 100049

    ?

    羌塘高寒草地物種多樣性與生態(tài)系統(tǒng)多功能關(guān)系格局

    熊定鵬1,2,趙廣帥1,2,武建雙1,石培禮1,*,張憲洲1

    1 中國科學院地理科學與資源研究所生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測與模擬重點實驗室, 北京100101 2 中國科學院大學, 北京100049

    摘要:傳統(tǒng)的生物多樣性-生態(tài)系統(tǒng)功能研究大多側(cè)重于單一生態(tài)系統(tǒng)功能與物種多樣性的關(guān)系,忽略了生態(tài)系統(tǒng)的重要價值在于其能夠同時提供多種功能或服務(wù),即生態(tài)系統(tǒng)的多功能性。基于藏北羌塘高寒草地樣帶調(diào)查數(shù)據(jù),選取植被地上生物量、地下生物量、土壤全氮、硝態(tài)氮及銨態(tài)氮含量、土壤全磷含量、土壤有機碳儲量等7個與植物生長、養(yǎng)分循環(huán)、土壤有機碳蓄積相關(guān)的參數(shù)來表征生態(tài)系統(tǒng)多功能性。采用上述參數(shù)轉(zhuǎn)換為Z分數(shù)后的平均值計算多功能性指數(shù)(M)。分析了不同生物多樣性指數(shù)與生態(tài)系統(tǒng)多功能指數(shù)的關(guān)系以及年降水量和年均溫度對物種多樣性和生態(tài)系統(tǒng)多功能性指數(shù)的影響。結(jié)果表明,物種豐富度指數(shù)與生態(tài)系統(tǒng)多功能性之間呈極顯著的正相關(guān)關(guān)系,Shannon-wiener和Simpson物種多樣性指數(shù)也與多功能性指數(shù)間呈顯著的正相關(guān),但多功能性指數(shù)與Pielou均勻度指數(shù)沒有表現(xiàn)出明顯的相關(guān)關(guān)系。物種豐富度與表征植物生長、養(yǎng)分循環(huán)以及土壤有機碳蓄積的生態(tài)系統(tǒng)功能指數(shù)間也均呈極顯著的正相關(guān)關(guān)系。降水格局顯著影響羌塘高原物種豐富度和生態(tài)系統(tǒng)多功能指數(shù),二者均隨年降雨量的增加而顯著增加,但物種多樣性指數(shù)并未與年降水量呈現(xiàn)顯著相關(guān)關(guān)系。研究強調(diào)了群落物種豐富度即群落物種數(shù)量對維持生態(tài)系統(tǒng)多功能性的重要意義,這意味著由于人類活動導(dǎo)致的物種喪失可能會給藏北高寒草地生態(tài)系統(tǒng)多功能和生態(tài)服務(wù)帶來更為嚴重的后果。就退化草地恢復(fù)或草地可持續(xù)管理而言,在藏北羌塘地區(qū),本地植物種的物種豐富度恢復(fù)和維持應(yīng)作為重要目標之一。

    關(guān)鍵詞:羌塘高原;高寒草地;物種多樣性;物種豐富度;生態(tài)系統(tǒng)多功能性

    在全球范圍內(nèi),由于人類活動導(dǎo)致的物種滅絕和生態(tài)系統(tǒng)結(jié)構(gòu)簡化的趨勢不斷加劇,生態(tài)系統(tǒng)中物種喪失可能帶來的潛在生態(tài)后果備受關(guān)注。近20年以來,生物多樣性與生態(tài)系統(tǒng)過程(功能)的關(guān)系成為生態(tài)學領(lǐng)域內(nèi)的核心科學問題。大量相關(guān)實驗研究在全球范圍內(nèi)廣泛開展,多數(shù)實驗結(jié)果表明生物多樣性與單一的生態(tài)系統(tǒng)功能(例如:生產(chǎn)力,氮素累積)之間存在因果關(guān)系[1- 3],多樣性導(dǎo)致更高的群落生產(chǎn)力、更高的系統(tǒng)穩(wěn)定性和更高的抗入侵能力[1,4-5]。對生物多樣性-生態(tài)系統(tǒng)功能關(guān)系實驗結(jié)果的整合分析也表明,物種多樣性與生態(tài)系統(tǒng)功能之間呈正相關(guān)關(guān)系,隨物種多樣性的增加這種關(guān)系通常會趨于飽和,但相關(guān)的機理性解釋仍存在較大爭議[6-8]。然而,這些生物多樣性-生態(tài)系統(tǒng)功能關(guān)系研究僅考慮了單一生態(tài)系統(tǒng)功能與物種多樣性的關(guān)系,卻忽略了生態(tài)系統(tǒng)最為重要價值在于其能夠同時維系或提供諸如凈初級生產(chǎn)力、碳氮循環(huán)、養(yǎng)分維持等多種不同的功能或服務(wù)[9-10],即生態(tài)系統(tǒng)的多功能性(multifunctionality)[11]。

    近年來,一些學者認為只有將各類的生態(tài)系統(tǒng)系統(tǒng)功能作為整體考慮,才能更利于全面地理解生態(tài)系統(tǒng)功能與生物多樣性的關(guān)系。此外,從生態(tài)系統(tǒng)多功能性角度出發(fā),也為生態(tài)系統(tǒng)管理提供了一個全新的、綜合性的視角。因為生態(tài)系統(tǒng)功能間存在著權(quán)衡(trade-off)關(guān)系,例如,植物生產(chǎn)力與抗脅迫能力間通常表現(xiàn)為負相關(guān),即不太可能同時最大化此兩種生態(tài)系統(tǒng)功能[12],擁有較高的某種生態(tài)功能的物種,其性狀通常限制了其他功能水平。因此僅考慮某一生態(tài)系統(tǒng)功能的管理策略可能會削弱生態(tài)系統(tǒng)維持或提供其他功能的能力。在此背景下,生物多樣性-生態(tài)系統(tǒng)多功能性關(guān)系研究逐漸成為一個新的熱點議題[9,12]。例如,Hector[9]基于歐洲草地BIODEPTH實驗數(shù)據(jù),選擇了7種生態(tài)系統(tǒng)過程參數(shù)包括地上、地下生產(chǎn)力、地上植被碳庫、土壤無機氮含量、光合有效輻射傳輸比率以及木質(zhì)素和纖維素降解速率,研究了物種多樣性與生態(tài)系統(tǒng)多功能性的關(guān)系。Zavaleta[12]等根據(jù)美國Cedar Creek長期實驗研究數(shù)據(jù),綜合考慮了地上、地下生產(chǎn)力、土壤無機氮含量、群落抗入侵能力、昆蟲物種多樣性、土壤碳庫和植物氮素含量8種生態(tài)系統(tǒng)功能,分析了不同時間尺度上物種多樣性與生態(tài)系統(tǒng)多功能性的關(guān)系。Maestre[10]等基于全球干旱地區(qū)的實驗結(jié)果研究生態(tài)系統(tǒng)多功能性與物種豐富度的關(guān)系,其選取了與生態(tài)系統(tǒng)碳、氮、磷循環(huán)相關(guān)的14種生態(tài)系統(tǒng)功能來表征生態(tài)系統(tǒng)多功能。以上生物多樣性-生態(tài)系統(tǒng)多功能性研究都認為維持或提供生態(tài)系統(tǒng)多功能性需要更多的物種數(shù)或功能群組合[9-10,12-13],可能原因是不同物種(或功能型)對不同的生態(tài)系統(tǒng)功能(過程)都具有重要影響。然而,目前國內(nèi)對生態(tài)系統(tǒng)的多功能研究還很少,只有對農(nóng)林生態(tài)系統(tǒng)多功能的概念框架進行了介紹[14-16],而對多功能性與生物多樣性的關(guān)系研究尤其缺乏。毫無疑問,這會局限了人們對生態(tài)系統(tǒng)多功能性的認識,不利于生態(tài)系統(tǒng)管理和生態(tài)服務(wù)的發(fā)揮。

    藏北高原地處青藏高原西北腹地,是西藏重要畜牧業(yè)生產(chǎn)基地,也是生態(tài)環(huán)境極為脆弱的典型區(qū)域[17]。草地是藏北地區(qū)最主要的生態(tài)系統(tǒng)類型,各類天然草地面積達4.8×105km2, 占藏北高原土地面積的81%[18]。高寒草地生態(tài)系統(tǒng)為藏北地區(qū)牧民提供了凈初級物質(zhì)生產(chǎn)、碳蓄積與碳匯、調(diào)節(jié)氣候、涵養(yǎng)水源、水土保持、防風固沙和維持生物多樣性等生態(tài)服務(wù)功能,是區(qū)域牧業(yè)經(jīng)濟發(fā)展的基礎(chǔ),也是青藏高原生態(tài)安全屏障的重要組成部分[19-21]。該區(qū)域以往生物多樣性-生態(tài)系統(tǒng)功能關(guān)系研究多集中于樣點尺度,且僅針對單一的生態(tài)系統(tǒng)功能,基于樣帶尺度的研究也僅限于分析地上生產(chǎn)力與生物多樣性的關(guān)系格局[22]。尚缺乏藏北天然草地生態(tài)系統(tǒng)多功能性和生物多樣性關(guān)系的定量研究。據(jù)此,本研究采用樣帶分析方法,基于樣地調(diào)查數(shù)據(jù),綜合7種與植被生產(chǎn)力、養(yǎng)分循環(huán)、土壤有機碳蓄積相關(guān)的生態(tài)系統(tǒng)參數(shù),分析了藏北高寒草地物種多樣性和功能多樣性與生態(tài)系統(tǒng)多功能性之間的關(guān)系。研究結(jié)果將加深人們對高寒草地生態(tài)系統(tǒng)多功能性與生物多樣性關(guān)系的認識,為草地生物多樣性保護以及可持續(xù)利用管理提供有益科學參考。

    1材料與方法

    1.1研究區(qū)概況

    藏北高原,又稱“羌塘”,意為北方高平地,位于西藏自治區(qū)北部。整個羌塘地區(qū),東迄91°E左右的內(nèi)外流水系之分水嶺,西至國界并與克什米爾地區(qū)相鄰接,南達岡底斯-念青唐古拉山脈,北到昆侖山脈,面積約為5.95×105km2[18]。研究區(qū)屬亞寒帶氣候,氣候寒冷干燥,海拔在4600—5100m之間,大部分地區(qū)年平均氣溫低于0℃, 最暖月均溫不及14℃, 最冷月均溫-10℃以下。受西南季風運動路徑及地形因素影響,降水自東向西呈遞減趨勢,其中65%—85%的年降水量集中在植物生長旺盛的6—8月。研究區(qū)內(nèi)部主要包括高寒草甸,高寒草原,高寒荒漠草原等地帶性高寒草地生態(tài)系統(tǒng)類型[23]。土壤類型主要包括栗鈣土、寒凍鈣土、雛漠土及寒漠土等。

    1.2野外取樣和室內(nèi)分析

    2011年7月下旬—8月初(植物生長旺盛期),在藏北羌塘高原東西樣帶上選取了20個圍欄樣地作為調(diào)查樣點(表1),研究樣地均為西藏自治區(qū)退牧還草工程劃定的禁牧草地,于2004或2005年起實施全年禁牧管理。在每個樣地選取一塊地勢相對平坦,植物生長均勻的區(qū)域(100m×100m)進行取樣調(diào)查,在該區(qū)域?qū)蔷€上設(shè)置一條 100m 樣線,沿樣線設(shè)置5個0.5m×0.5m的樣方,相當于5次重復(fù)。在各樣方中,分別記錄樣方內(nèi)植物種類、分物種蓋度、平均高度和群落總蓋度,采用齊地刈割法收獲分種地上生物量;同時在每個樣點隨機拋0.1m2樣圈30次,記錄所出現(xiàn)的物種和次數(shù),用以計算樣地的物種豐富度和各物種的相對頻度。每個樣方,地上部分刈割后,用根鉆(直徑=7cm)取4鉆(0—20cm深度)帶根系的土樣,混合后裝入布袋中,樣品帶回住處后立即沖洗漂凈并將根系分離;隨后用直徑4cm的土鉆獲取5鉆(0—20cm深度)土壤樣品,土壤樣品混合后一部分裝入塑封袋中放入冰箱冷凍,一部分裝入布袋待測。此外,在樣線旁挖掘一個容重坑,用環(huán)刀法獲取土壤容重樣品,用以測定土壤容重及礫石比。

    1.3數(shù)據(jù)處理1.3.1物種重要值

    重要值(IV)是群落物種的重要定量指標, 本研究采用草地調(diào)查樣方內(nèi)物種相對蓋度(Cr),相對高度(Hr)和樣圈法統(tǒng)計的物種相對頻度(Fr)計算重要值[23],公式如下:

    IV=(Cr+Hr+Fr)/3 (1)

    1.3.2物種多樣性

    本研究采用物種豐富度指數(shù)(Species richness,S),Shannon-wiener多樣性指數(shù)(H′),Simpson優(yōu)勢度指數(shù)(P),Pielou均勻度指數(shù)(E)來表征物種多樣性[26]。

    (2)

    (3)

    (4)

    式中,S為群落物種數(shù)目,Pi為第i種植物的相對重要值。

    1.3.3生態(tài)系統(tǒng)多功能性指數(shù)

    選擇與植物生長(地上生物量、地下生物量)、養(yǎng)分循環(huán)(土壤全氮含量、土壤氨態(tài)氮含量、土壤硝態(tài)氮含量和土壤全磷含量)和土壤有機碳蓄積(土壤有機碳密度)相關(guān)的7種重要生態(tài)系統(tǒng)參數(shù)來表征生態(tài)系統(tǒng)多功能性(M)。為量化生態(tài)系統(tǒng)多功能性,首先計算20個調(diào)查樣地7種生態(tài)系統(tǒng)參數(shù)的Z分數(shù),單樣本Kolmogorov-Smirnov檢驗結(jié)果表明各參數(shù)均服從正態(tài)分布(P> 0.05)(表2),因此直接使用原始數(shù)據(jù)進行計算。Z分數(shù)計算公式如下:

    Zij=(xij-μj)/σj

    (5)

    式中,Zij為樣地i第j種生態(tài)系統(tǒng)參數(shù)的Z分數(shù),i的范圍在1—7之間,j的范圍在1—20之間;xij為樣地i第j種生態(tài)系統(tǒng)參數(shù)的數(shù)值;μj為第j種生態(tài)系統(tǒng)參數(shù)在20個樣地間的平均值;σj為第j種生態(tài)系統(tǒng)參數(shù)在20個樣地間的標準差。

    對每一樣地7種不同生態(tài)系統(tǒng)參數(shù)的Z分數(shù)求平均值,即為各樣地的多功能多樣性指數(shù)(M)。計算公式如下:

    (6)

    式中,Mi為樣地i的生態(tài)系統(tǒng)多功能性指數(shù)。該方法的優(yōu)點是,其計算得到多功能性指數(shù)具有較好的統(tǒng)計學特性,且不會改變原始數(shù)據(jù)的變異程度,計算簡便,能夠直觀的表征植物群落同時提供多種生態(tài)系統(tǒng)功能或服務(wù)的能力[10]。

    表2 羌塘樣帶各樣地調(diào)查指標描述性統(tǒng)計

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

    通過一元線性回歸方法分析物種多樣性和功能多樣性指數(shù)與生態(tài)系統(tǒng)多功能性指數(shù)間的關(guān)系,此外,將7個變量劃分為與植物生長、養(yǎng)分循環(huán)、土壤有機碳蓄積相關(guān)的3種重要生態(tài)系統(tǒng)功能,并分別與不同生物多樣性指數(shù)進行線性回歸。同時分析了年降水量、年均溫度與物種多樣性和生態(tài)系統(tǒng)功能多樣性指數(shù)間的線性回歸關(guān)系。以上分析采用最小二乘法擬合回歸方程。所有數(shù)據(jù)分析在SPSS 20.0中進行。

    2結(jié)果與分析

    線性回歸分析結(jié)果表明(表3,圖1),藏北高寒草地生態(tài)系統(tǒng)多功能性指數(shù)隨物種豐富度的增加而呈極為顯著的增加趨勢,物種豐富度解釋了約55%的草地生態(tài)系統(tǒng)多功能性指數(shù)的變化。Shannon-wiener多樣性指數(shù)和Simpson優(yōu)勢度指數(shù)也與藏北高寒草地生態(tài)系統(tǒng)多功能性指數(shù)間呈顯著的正相關(guān)關(guān)系(P<0.05),Pielou均勻度指數(shù)與多功能性指數(shù)間擬合結(jié)果不顯著(表3,圖1)。

    表3 生物多樣性與生態(tài)系統(tǒng)多功能性線性回歸的F值和顯著性

    *P<0.05, **P<0.01, ***P<0.001

    圖1 高寒草地物種多樣性與生態(tài)系統(tǒng)多功能性的關(guān)系Fig.1 Relationship between species diversity and ecosystem multifunctionality in alpine grasslands

    物種豐富度和表征植物生長、養(yǎng)分循環(huán)以及土壤有機碳蓄積3種重要的生態(tài)系統(tǒng)功能指數(shù)間均呈極顯著的正相關(guān)關(guān)系,且線性擬合的P值均小于0.01,在統(tǒng)計上達到了99%的置信度水平(表3,圖2)。與植物生長和土壤有機碳蓄積相關(guān)的生態(tài)系統(tǒng)功能指數(shù)呈現(xiàn)隨Shannon-wiener指數(shù)和Simpson指數(shù)增加而顯著增加的趨勢,但與養(yǎng)分循環(huán)相關(guān)的生態(tài)系統(tǒng)功能指數(shù)與兩種多樣性指數(shù)間無明顯相關(guān)關(guān)系(表3,圖2)。Pielou均勻度指數(shù)與3種生態(tài)系統(tǒng)功能指數(shù)間均未表現(xiàn)出明顯的相關(guān)關(guān)系。兩種功能多樣性指數(shù)與3種生態(tài)系統(tǒng)功能指數(shù)間的線性擬合結(jié)果均不顯著(表3)。

    圖2 高寒草地物種多樣性與不同生態(tài)系統(tǒng)功能的關(guān)系Fig.2 Relationship between species diversity and different ecosystem functions in the alpine grasslands

    羌塘高原草地群落的物種豐富度隨年降雨量的增加顯著增加(R2=0.27,P<0.05),但年降雨量與Shannon-wienner指數(shù)、Simpson指數(shù)和Pielou指數(shù)間沒有呈現(xiàn)出明顯的相關(guān)關(guān)系。年降水量顯著的影響生態(tài)系統(tǒng)多功能性及其相關(guān)組分,其與多功能性綜合指數(shù)和表征植物生長、土壤養(yǎng)分循環(huán)相關(guān)的功能性指數(shù)間均呈極其顯著的正相關(guān)關(guān)系(P<0.001),與土壤有機碳蓄積指數(shù)間呈顯著的正相關(guān)關(guān)系(圖3)。而物種多樣性指數(shù)和多功能性指數(shù)與溫度之間均沒有相關(guān)關(guān)系??梢?藏北羌塘高原降水梯度上多功能性及其組分指數(shù)隨年降水量增加而顯著增加的格局主要是降水因子決定的。

    圖3 高寒草地物種多樣性與年降水量的關(guān)系Fig.3 Relationship between ecosystem multifunctionality and mean annual precipitation in the alpine grasslands

    3討論

    在藏北羌塘高寒草地,不同物種多樣性指數(shù)與草地生態(tài)系統(tǒng)多功能性指數(shù)間關(guān)系格局存在較大差異。本研究發(fā)現(xiàn)物種豐富度指數(shù)與多功能性指數(shù)間線性回歸關(guān)系最顯著,而Pielou均勻度指數(shù)與生態(tài)系統(tǒng)多能性指數(shù)關(guān)系不顯著。物種豐富度指數(shù)代表了植物群落內(nèi)物種數(shù)目的多少,均勻度指數(shù)則體現(xiàn)了各物種個體數(shù)目分配的均勻程度。因此,這一結(jié)果說明,藏北羌塘高寒草地生態(tài)系統(tǒng)多功能性主要受群落物種數(shù)量而非物種分布的均勻程度影響。在羌塘高原,草地生態(tài)系統(tǒng)同時維持和提供多種生態(tài)系統(tǒng)功能和服務(wù)需要較高的植物種數(shù)量來支撐。本研究結(jié)果與諸多基于區(qū)域或全球尺度的研究相一致,例如,Hector[9]等基于歐洲BIODEPTH草地生物多樣性-生態(tài)系統(tǒng)功能實驗數(shù)據(jù)的研究結(jié)果,Zavaleta[12]等基于全球干旱地區(qū)實驗數(shù)據(jù)的研究以及Lefcheck[27]等在全球尺度上跨越不同營養(yǎng)級和生境的整合分析研究。這些研究都表明物種豐富度的增加對生態(tài)系統(tǒng)多功能性具有正效應(yīng)。Shannon-wienner和Simpson多樣性指數(shù)指數(shù)通過對群落的豐富度和均勻度進行綜合評價來表征群落的復(fù)雜程度,本研究發(fā)現(xiàn)兩種指數(shù)也與生態(tài)系統(tǒng)多功能性指數(shù)間成顯著的正相關(guān)關(guān)系(P<0.05)。這種正相關(guān)關(guān)系表明,生態(tài)系統(tǒng)的多功能性主要受物種豐富度的影響。

    本研究中,植物初級生產(chǎn)力、土壤碳蓄積,以及養(yǎng)分循環(huán)這3類重要的生態(tài)系統(tǒng)功能指數(shù)均隨物種多樣性指數(shù)特別是物種豐富度的增加而顯著增加。研究結(jié)果與大多數(shù)的物種多樣性-生產(chǎn)力試驗結(jié)果相一致,即物種豐富的群落具有更高的群落生產(chǎn)力[1,3]。對土壤有機碳固持而言,物種豐富度的增加通常通過提高地上、地下生物量生產(chǎn)力[28]和根系凋落物碳庫來增加土壤碳累積[29],另一方面,物種多樣性能夠通過改變凋落物的構(gòu)成和多樣性以及土壤的生物/非生物屬性來影響土壤有機質(zhì)分解[29],從而提高土壤有機碳蓄積。有研究表明,植物多樣性越高,會增加根際微生物碳源的輸入,從而導(dǎo)致微生物活性增強和凋落物分解,有利于碳儲量增加[30]。物種豐富度與土壤養(yǎng)分循環(huán)指數(shù)間的正相關(guān)關(guān)系可能是由于較高的物種豐富度增加了土壤潛在氮礦化效率[31],同時較高的土壤有機碳蓄積有利于土壤中氨離子的交換,因此較高的土壤儲量通常與較高土壤N礦化效率相聯(lián)系[31-32]。這表明本研究中3種重要的生態(tài)系統(tǒng)功能/過程間存在一定的協(xié)同效應(yīng),維持較高的物種多樣性水平不但能夠提高藏北羌塘地區(qū)高寒草地生態(tài)系統(tǒng)多功能性,而且可以同時最大化這些重要的生態(tài)系統(tǒng)功能和服務(wù)。Wu[22]等人關(guān)于藏北羌塘高原草地地上生物量-物種豐富度關(guān)系的研究揭示了與以往諸多研究相一致的經(jīng)典的單峰曲線格局。這一結(jié)果意味著:僅針對維持地上生產(chǎn)力水平而言,草地群落中某些植物種可能在功能上表現(xiàn)冗余,也就是對群落生產(chǎn)力的貢獻不大。然而,從生態(tài)系統(tǒng)多功能性角度出發(fā),本研究發(fā)現(xiàn)羌塘高寒草地生態(tài)系統(tǒng)多功能性隨物種豐富度的增加呈顯著線性增加趨勢。這說明在羌塘高原,高寒草地生態(tài)系統(tǒng)維持多種功能或服務(wù)需要更多物種的參與,因為不同物種可能會在對維系不同種類的生態(tài)系統(tǒng)功能都產(chǎn)生貢獻[9,13,33]。因次,孤立的考慮某一生態(tài)系統(tǒng)功能的研究可能會低估生物多樣性的貢獻[9,27]。由于不同生態(tài)系統(tǒng)功能之間可能存在的權(quán)衡關(guān)系,使用生態(tài)系統(tǒng)多功能性這一綜合性概念有助于從整體上理解和認識生物多樣性與生態(tài)系統(tǒng)功能間的關(guān)系格局[34]。

    氣候因子對于草地植物群落結(jié)構(gòu)和物種豐富度的空間分布格局的形成具有重要意義[35],藏北羌塘草地樣帶橫跨十個經(jīng)度,年降水量從東至西由550mm遞減至106mm,而在緯度上跨度很小,溫度的影響很小。因此水分梯度變化可能是本研究中影響物種多樣性變化的重要因素。研究結(jié)果表明,物種豐富度與年降水量間呈顯著的正相關(guān)關(guān)系,這與楊元合等[36]對整個青藏高原高寒草地樣帶和胡云峰等[37]對內(nèi)蒙古東北-西南草地樣帶多樣性-降水量關(guān)系的研究結(jié)果相符。其他與均勻度緊密相關(guān)的多樣性指數(shù)均與年降雨量間無明顯的相關(guān)關(guān)系,表明均勻度更大程度上受到群落結(jié)構(gòu)自身的影響,而對降水梯度變化不太敏感[37]。本研究中生態(tài)系統(tǒng)多功能性指數(shù)以及植被生產(chǎn)力、土壤碳儲量和養(yǎng)分循環(huán)相關(guān)的參數(shù)均與年降水量間呈極顯著的正相關(guān)關(guān)系,表明水分條件是決定羌塘高寒草原植物生長、養(yǎng)分循環(huán)和碳蓄積的主要限制因子。

    藏北羌塘地區(qū)氣候條件惡劣,環(huán)境脆弱,其草地對氣候變化和人類活動的干擾響應(yīng)極為敏感,生態(tài)系統(tǒng)脆弱。因此,藏北高原高寒草地由于過度放牧利用和氣候變化的影響,威脅到群落的多樣性,從而對生態(tài)系統(tǒng)多功能性造成顯著地影響[38]。本研究結(jié)果強調(diào)了物種豐富度,即群落的物種數(shù)目對于維持藏北高寒草地生態(tài)系統(tǒng)多功能性的重要意義。較高的物種多樣性水平有利于提高植物生產(chǎn)力、土壤有機碳固定以及土壤養(yǎng)分循環(huán)。相應(yīng)的,物種喪失可能會給此區(qū)域帶來更加嚴重的生態(tài)后果,威脅當?shù)啬撩竦馁囈陨娴奈镔|(zhì)經(jīng)濟基礎(chǔ),以及野生動物的生存。因此,在退化草地恢復(fù)過程中,生物多樣性水平特別是本地植物種豐富度的恢復(fù)應(yīng)當作為一個重要目標。就草地可持續(xù)管理而言,應(yīng)多從生態(tài)系統(tǒng)多功能性角度出發(fā),而不是單純考慮某一生態(tài)系統(tǒng)功能或服務(wù)的實現(xiàn)。

    4結(jié)論

    本研究表明,物種豐富度水平對于維持羌塘高寒草地生態(tài)系統(tǒng)多功能性具有重要意義,物種豐富度可以同時最大化藏北高寒草地的主要生態(tài)系統(tǒng)功能和服務(wù)。孤立的考慮某一生態(tài)系統(tǒng)功能的研究可能會低估生物多樣性的貢獻。由于人類活動導(dǎo)致的物種喪失可能會給藏北羌塘草原帶來更加嚴重的生態(tài)后果,生物多樣性水平(特別是本地物種豐富度)的恢復(fù)和維持應(yīng)當作為該地區(qū)退化草地恢復(fù)的重要目標之一。

    參考文獻(References):

    [1]Hector A, Schmid B, Beierkuhnlein C, CaldeiraM C, Diemer M, Dimitrakopoulos P G, Finn J A, Freitas H, Giller P S, Good J, Harris R, H?gberg P, Huss-Danell K, Joshi J, Jumpponen A, K?rner C, LeadleyP W, Loreau M, Minns A, Mulder C P H, O′Donovan G, Otway S J, Pereira J S, Prinz A, Read D J, Scherer-Lorenzen M, Schulze E D, Siamantziouras A S D, Spehn E M, Terry A C, Troumbis A Y, Woodward F I, Yachi S, Lawton J H. Plant diversity and productivity experiments in European grasslands. Science, 1999, 286(5442): 1123- 1127.

    [2]Hooper D U, Chapin F S, Ewel J J, Hector A, InchaustiP, Lavorel S, Lawton J H, Lodge D M, Loreau M, Naeem S, Schmid B, Set?l? H, SymstadA J, Vandermeer J, Wardle D A. Effects of biodiversity on ecosystem functioning: A consensus of current knowledge. Ecological Monographs, 2005, 75(1): 3- 35.

    [3]Flombaum P, Sala O E. Higher effect of plant species diversity on productivity in natural than artificial ecosystems. Proceedings of the National Academy of Sciences of the United States of America, 2008, 105(16): 6087- 6090.

    [4]張全國, 張大勇. 生物多樣性與生態(tài)系統(tǒng)功能: 最新的進展與動向. 生物多樣性, 2003, 11(5): 351- 363.

    [5]Kennedy T, Naeem S, Howe K M, Knops J M H, Tilman D, Reich P. Biodiversity as a barrier to ecologicalinvasion. Nature, 2002, 417(6889): 636- 638.

    [6]Srivastava D S, Velland M. Biodiversity-ecosystem function research: is it relevant toconservation? Annual Review of Ecology, Evolution, and Systematics, 2005, 36(1): 267- 294.

    [7]Balvanera P, Pfisterer A B, Buchmann N, He JS, Nakashizuka T, Raffaelli D, Schmid B. Quantifying the evidence for biodiversity effects on ecosystem functioning and services. Ecology Letters, 2006, 9(10): 1146- 1156.

    [8]Cardinale B J, Srivastava D S, Duffy J E,Wright J P, Downing A L, Sankaran M, Jouseau C. Effects of biodiversity on the functioning of trophic groups and ecosystems. Nature, 2006, 443(7114): 989- 992.

    [9]Hector A, Bagchi R. Biodiversity and ecosystem multifunctionality. Nature, 2007, 448(7150): 188- 190.

    [10]Maestre F T, Quero J L, Gotelli N J, Escudero A, Ochoa V, Delgado-Baquerizo M, Garcia-Gomez M, Bowker M A, Soliveres S, Escolar C, Garcia-Palacios P, Berdugo M, Valencia E, Gozalo B, Gallardo A, Aguilera L, Arredondo T, Blones J, Boeken B, Bran D, Conceicao A A, Cabrera O, Chaieb M, Derak M, Eldridge D J, Espinosa C I, Florentino A, Gaitan J, Gatica M G, Ghiloufi W, Gomez-Gonzalez S, Gutierrez J R, Hernandez R M, Huang X, Huber-Sannwald E, Jankju M, Miriti M, Monerris J, Mau R L, Morici E, Naseri K, Ospina A, Polo V, Prina A, Pucheta E, Ramirez-Collantes D A, Romao R, Tighe M, Torres-Diaz C, Val J, Veiga J P, Wang D, Zaady E. Plant species richness and ecosystem multifunctionality in global drylands. Science, 2012, 335(6065): 214- 218.

    [11]Byrnes J E K, Gamfeldt L, Isbell F, Lefcheck J S, Griffin J N, Hector A, Cardinale B J, Hooper D U, Dee L E, Duffy J E, Freckleton R. Investigating the relationship between biodiversity and ecosystem multifunctionality: Challenges and solutions. Methods in Ecology and Evolution, 2014, 5(2): 111- 124.

    [12]Zavaleta E S, Pasari J R, Hulvey K B, Tilman G D. Sustaining multiple ecosystem functions in grassland communities requires higher biodiversity. Proceedings of the National Academy of Sciences of the United States of America, 2010, 107(4): 1443- 1446.

    [13]Gamfeldt L, Hillebrand H, Jonsson PR. Multiple functions increase the importance of biodiversity for overall ecosystem functioning. Ecology, 2008, 89(5): 1223- 1231.

    [14]江娟. 北京市延慶縣森林資源及其多功能效益評價[D]. 北京: 中國林業(yè)科學研究院, 2010.

    [15]彭建, 劉志聰, 劉焱序. 農(nóng)業(yè)多功能性評價研究進展. 中國農(nóng)業(yè)資源與區(qū)劃, 2014, 35(6): 1- 8.

    [16]吳晶晶. 林農(nóng)復(fù)合經(jīng)營系統(tǒng)多功能性分析及評價研究[D]. 南京: 南京林業(yè)大學, 2013.

    [17]王景升, 張憲洲, 趙玉萍, 沈振西, 石培禮, 余成群. 藏北羌塘高原氣候變化的時空格局. 資源科學, 2008, 30(12): 1852- 1859.

    [18]李明森. 藏北高原草地資源合理利用. 自然資源學報, 2000, 15(4): 335- 339.

    [19]劉興元, 龍瑞軍, 尚占環(huán). 草地生態(tài)系統(tǒng)服務(wù)功能及其價值評估方法研究. 草業(yè)學報, 2011, 20(1): 167- 174.

    [20]宋春橋, 游松財, 柯靈紅, 劉高煥, 鐘新科. 藏北高原地表覆蓋時空動態(tài)及其對氣候變化的響應(yīng). 應(yīng)用生態(tài)學報, 2011, 22(8): 2091- 2097.

    [21]孫鴻烈, 鄭度, 姚檀棟, 張鐿鋰. 青藏高原國家生態(tài)安全屏障保護與建設(shè). 地理學報, 2012, 67(1): 3- 12.

    [22]Wu J S, Shen Z X, Zhang X Z. Precipitation and species composition primarily determinethe diversity-productivity relationship of alpine grasslandson the Northern Tibetan Plateau. Alpine Botany, 2014, 124(1): 13- 25.

    [23]武建雙, 李曉佳, 沈振西, 張憲洲, 石培禮, 余成群, 王景升, 周宇庭. 藏北高寒草地樣帶物種多樣性沿降水梯度的分布格局. 草業(yè)學報, 2012, 21(3): 17- 25.

    [24]中國科學院南京土壤研究所. 土壤理化分析. 上海: 上海科技出版社, 1978.

    [25]陳杰. 貴州省植煙土壤礫石含量分布. 安徽農(nóng)業(yè)科學,2007, 35(32): 10334- 10335.

    [26]Tilman D, Elhaddi A. Drought and biodiversity in grasslands. Oecologia, 1992, 89(2): 257- 264.

    [27]Lefcheck J S, Byrnes J E K, Isbell F, Gamfeldt L, Griffin J N, Eisenhauer N, Hensel M J S, Hector A, CardinaleB J, Duffy J E. Biodiversity enhances ecosystem multifunctionality across trophic levels and habitats. Nature Communication, 2015, 6: 6936.

    [28]FornaraDA, Tilman D. Plant functional composition influences rates of soilcarbon and nitrogen accumulation. Journal of Ecology, 2008, 96(2): 314- 322.

    [29]Cong W F, Van Ruijven J, Van Der Werf W, De Deyn G B, Mommer L, Berendse F, Hoffland E. Plant species richness leaves a legacy of enhanced root litter-induceddecomposition in soil. Soil Biology and Biochemistry, 2015, 80: 341- 348.

    [30]Lange M, Eisenhauer N, Sierra CA, Bessler H, Engels C, Griffiths RI, Mellado-VazquezPG, Malik AA, Roy J, Scheu S, Steinbeiss S,Thomson B C,Trumbore S E, Gleixner G. Plant diversity increases soil microbial activity and soil carbon storage. Nature Communications, 2014, 6: 6707.

    [31]Mueller K, Hobbie S E, Tilman D, Reich P B. Effects of plant diversity, N fertilization, and elevatedcarbon dioxide on grassland soil N cycling in along-term experiment. Global Change Biology, 2013, 19(4): 1249- 1261.

    [32]Cong W F, Van Ruijven J, Mommer L, De Deyn G B, Berendse F, Hoffland E, Lavorel S. Plant species richness promotes soil carbon andnitrogen stocks in grasslands without legumes. Journal of Ecology, 2014, 102(5): 1163- 1170.

    [33]Isbell F, Calcagno V, Hector A, Connolly J, Harpole W S, Reich P B, Scherer-Lorenzen M, Schmid B, Tilman D, Van RuijvenJ, Weigelt A, Wilsey B J, Zavaleta E S, Loreau M.High plant diversity is needed to maintain ecosystem services. Nature, 2011, 477(7363): 199- 202.

    [34]Byrnes J, Lefcheck J S, Gamfeldt L, Griffin JN, Isbell F, Hector A. Multifunctionality does not imply that all functions are positively correlated. Proceedings of the National Academy of Sciences of the United States of America, 2014, 111(51): E5490.

    [35]Cleland EE, CollinsSL, Dickson TL, FarrerEC, GrossKL, Gherardi L A, Hallett L M, Hobbs R J, Hsu J S, Turnbull L, Suding K N. Sensitivity of grassland plant community compositionto spatial vs. temporal variation in precipitation. Ecology, 2013, 94(8): 1687- 1696.

    [36]楊元合, 饒勝, 胡會峰, 陳安平, 吉成均, 朱彪, 左聞韻, 李軒然, 沈?;? 王志恒, 唐艷鴻, 方精云. 青藏高原高寒草地植物物種豐富度及其與環(huán)境因子和生物量的關(guān)系. 生物多樣性, 2004, 12(1): 200- 205.

    [37]胡云鋒, 艷燕, 阿拉騰圖雅, 畢力格吉夫. 內(nèi)蒙古東北-西南草地樣帶植物多樣性變化. 資源科學, 2012, 34(6): 1024- 1031.

    [38]Gao Q Z, Li Y, Wan Y F, Jiangcun W Z, Qin X B, Wang B S. Significant achievements in protection and restoration of alpine grassland ecosystem in northern Tibet, China. Restoration Ecology, 2009, 17(3): 320- 323.

    基金項目:國家自然科學基金項目(41271067); 中國科學院戰(zhàn)略性先導(dǎo)專項課題(XDA05060700)

    收稿日期:2015- 06- 30;

    修訂日期:2015- 11- 13

    *通訊作者

    Corresponding author.E-mail: shipl@igsnrr.ac.cn

    DOI:10.5846/stxb201506301361

    The relationship between species diversity and ecosystem multifunctionality in alpine grasslands on the Tibetan Changtang Plateau

    XIONG Dingpeng1,2, ZHAO Guangshuai1,2, WU Jianshuang1, SHI Peili1,*,ZHANG Xianzhou1

    1KeyLaboratoryofEcosystemNetworkObservationandModeling,InstituteofGeographicSciencesandNaturalResourcesResearch,Chinese,Beijing100101,China2UniversityofChineseAcademyofSciences,Beijing100049,China

    Abstract:The classic species diversity-ecosystem function relationship studies attached more importance on single function and its correlation with diversity despite the fact that most ecosystems are managed or valued for their ability to maintain multiple goods and services simultaneously, so called multifunctionality. Therefore, multifunctionality is an important perspective for biodiversity-function research in point of view of comprehensive function analysis. We conducted an ecosystem multifunctionality vs. biodiversity study on the Tibetan Changtang Plateau, in terms of seven indicators of ecosystem functioning related to plant growth (aboveground biomass, belowground biomass), nutrient cycling (soil total nitrogen content, soil and content, soil phosphorus content), and soil carbon storage (soil organic carbon density). Z scores (standardized deviates) of the seven indicators were assessed to represent multifunctionality. The multifunctionality index (M) was calculated as the average Z-score for all indicators measured in each site. Linear regression was used to analyze the relationship between multifunctionality and species diversity indices (Species richness, Shannon-wiener index, Simpson dominance index, Pielou evenness index). We also tested whether ecosystem multifunctionality and species diversity were associated with climatic factors (mean annual temperature and mean annual precipitation) by using linear regression. It was indicated that ecosystem multifunctionality indices were correlated most significantly with species richness, notable positive with Shannon-wiener and Simpson index, but not with Pielou evenness index. Separate analyses of functions related to the plant growth, soil carbon storage and nutrient cycling also yield significantly positive relationships with species richness. Precipitation plays a positive role in affecting species diversity and ecosystem multifunctionality indices, as the species diversity and multiple functions increase linearly with mean annual precipitation. The present study highlights the importance of species richness in maintaining multiple ecosystem functions, goods and services in Changtang Plateau. Because each species contribute to different functions, studies with focus on single function (i.e. net production) might underestimate levels of species richness required to maintain multiple ecosystem goods and services. The alpine grassland ecosystem is vulnerable to climate change and human induced biodiversity loss. It is therefore implied that species loss may cause severe consequence for ecosystem functioning, which would threaten ecosystem health and human well-being in this region. We suggest that the species richness of native species must be maintained and sustained in order to meet the goal of degraded grassland restoration and sustainable management on the Changtang Plateau.

    Key Words:Tibetan Changtang Plateau; alpine grasslands; species diversity; species richness; ecosystem multifunctionality

    熊定鵬,趙廣帥,武建雙,石培禮,張憲洲.羌塘高寒草地物種多樣性與生態(tài)系統(tǒng)多功能關(guān)系格局.生態(tài)學報,2016,36(11):3362- 3371.

    Xiong D P, Zhao G S, Wu J S, Shi P L,Zhang X Z.The relationship between species diversity and ecosystem multifunctionality in alpine grasslands on the Tibetan Changtang Plateau.Acta Ecologica Sinica,2016,36(11):3362- 3371.

    猜你喜歡
    物種多樣性
    小麥—苜蓿鄰作帶步甲科物種多樣性及其時空動態(tài)
    關(guān)于我國海洋生物物種多樣性的探討
    山東棉田雜草組成及其群落特征
    增益生防在大田作物害蟲生物防治中的應(yīng)用
    江西大崗山木荷天然林群落物種多樣性分析
    綠色科技(2016年21期)2016-12-27 10:53:26
    納板河保護區(qū)不同林分群落喬木物種多樣性研究
    用固定樣地法監(jiān)測銀瓶山闊葉林物種多樣性
    輝騰錫勒風電場區(qū)域內(nèi)物種多樣性研究
    桂南地區(qū)植物群落多樣性分析
    綠色科技(2015年6期)2015-08-05 18:12:16
    大亮子河國家森林公園林喬灌草不同層次植物物種組成分析
    久久久成人免费电影| 欧美丝袜亚洲另类 | aaaaa片日本免费| 国产黄色小视频在线观看| 国产黄色小视频在线观看| 一边摸一边抽搐一进一小说| 最后的刺客免费高清国语| 又爽又黄无遮挡网站| 99久久成人亚洲精品观看| 久久人妻av系列| 成人美女网站在线观看视频| 亚洲午夜理论影院| 亚洲真实伦在线观看| 亚洲自拍偷在线| 他把我摸到了高潮在线观看| a级毛片免费高清观看在线播放| 午夜亚洲福利在线播放| 99在线视频只有这里精品首页| x7x7x7水蜜桃| 乱码一卡2卡4卡精品| 国内精品久久久久精免费| 91狼人影院| 国内精品一区二区在线观看| 男人舔女人下体高潮全视频| 九九热线精品视视频播放| 69人妻影院| 亚洲精品影视一区二区三区av| 国产精品久久久久久久久免| 欧美潮喷喷水| .国产精品久久| 99久久久亚洲精品蜜臀av| 日本色播在线视频| 亚洲熟妇熟女久久| 亚洲精品影视一区二区三区av| 国产成人av教育| 亚洲avbb在线观看| 久久久久免费精品人妻一区二区| 欧美三级亚洲精品| www日本黄色视频网| 亚洲av美国av| 久久6这里有精品| 99riav亚洲国产免费| 亚洲欧美日韩高清专用| 99久久精品一区二区三区| 日本熟妇午夜| 久久精品国产亚洲av香蕉五月| 美女高潮喷水抽搐中文字幕| 欧美zozozo另类| 亚洲国产色片| 久久精品国产鲁丝片午夜精品 | 成年女人毛片免费观看观看9| 性色avwww在线观看| 久久精品影院6| 精品久久久久久久久亚洲 | 深夜精品福利| 精品人妻视频免费看| 免费在线观看影片大全网站| 人妻制服诱惑在线中文字幕| 99热精品在线国产| 国产久久久一区二区三区| 亚洲图色成人| 亚洲成av人片在线播放无| 亚洲av五月六月丁香网| 嫩草影视91久久| 国产不卡一卡二| 国产国拍精品亚洲av在线观看| 日韩欧美免费精品| 亚洲国产日韩欧美精品在线观看| 美女大奶头视频| 欧美精品国产亚洲| 色播亚洲综合网| 久久午夜亚洲精品久久| 老司机午夜福利在线观看视频| 乱码一卡2卡4卡精品| 韩国av一区二区三区四区| 亚洲av二区三区四区| 久久久午夜欧美精品| 亚洲av.av天堂| 欧美一区二区国产精品久久精品| 日本一二三区视频观看| 精品久久久久久成人av| av在线亚洲专区| 亚洲人成网站在线播放欧美日韩| 少妇被粗大猛烈的视频| 午夜福利在线在线| 亚洲国产欧美人成| 99riav亚洲国产免费| 亚洲欧美激情综合另类| 精品一区二区三区av网在线观看| 少妇高潮的动态图| 久久欧美精品欧美久久欧美| 91在线精品国自产拍蜜月| 久久久国产成人免费| 日本黄色片子视频| 天堂影院成人在线观看| 内射极品少妇av片p| 欧美三级亚洲精品| 色综合亚洲欧美另类图片| 国产精品自产拍在线观看55亚洲| 日韩高清综合在线| 日本在线视频免费播放| 在线a可以看的网站| 神马国产精品三级电影在线观看| 美女黄网站色视频| 国模一区二区三区四区视频| 99热这里只有精品一区| 色5月婷婷丁香| 亚洲国产高清在线一区二区三| 久久人人精品亚洲av| 波多野结衣高清作品| 成人午夜高清在线视频| 国产69精品久久久久777片| 国产色婷婷99| 动漫黄色视频在线观看| 久久中文看片网| 亚洲av一区综合| 给我免费播放毛片高清在线观看| 久久精品国产鲁丝片午夜精品 | 波多野结衣高清作品| 一夜夜www| 久久久色成人| 亚洲第一电影网av| 国产精品一区二区三区四区免费观看 | 乱码一卡2卡4卡精品| 成人高潮视频无遮挡免费网站| 国产免费一级a男人的天堂| 亚洲美女视频黄频| АⅤ资源中文在线天堂| 国模一区二区三区四区视频| 真实男女啪啪啪动态图| 看十八女毛片水多多多| 精华霜和精华液先用哪个| 国产精品日韩av在线免费观看| 免费观看的影片在线观看| 女的被弄到高潮叫床怎么办 | 国产精品免费一区二区三区在线| 欧美成人性av电影在线观看| 免费看av在线观看网站| 韩国av一区二区三区四区| 免费在线观看影片大全网站| 亚洲av免费在线观看| 嫩草影院新地址| 国产精品一区www在线观看 | 成人毛片a级毛片在线播放| 男女啪啪激烈高潮av片| 国产精品人妻久久久影院| 色播亚洲综合网| 99热这里只有是精品在线观看| 欧美日韩亚洲国产一区二区在线观看| 少妇的逼水好多| 国产老妇女一区| eeuss影院久久| 听说在线观看完整版免费高清| 春色校园在线视频观看| 久久久久久久久久久丰满 | 日本与韩国留学比较| 亚洲自偷自拍三级| 精品无人区乱码1区二区| 成人一区二区视频在线观看| 12—13女人毛片做爰片一| 久久久久久久久大av| 黄片wwwwww| 变态另类成人亚洲欧美熟女| 国产精品不卡视频一区二区| 18+在线观看网站| 熟妇人妻久久中文字幕3abv| 在线播放无遮挡| 啦啦啦观看免费观看视频高清| 日韩欧美在线乱码| 成人欧美大片| 国产精品野战在线观看| 99久久精品一区二区三区| 中文字幕av在线有码专区| 人人妻人人看人人澡| 亚洲,欧美,日韩| 国产在线男女| 国产精品99久久久久久久久| 波多野结衣高清无吗| 搡老岳熟女国产| 国产精品国产三级国产av玫瑰| 欧美人与善性xxx| 免费看a级黄色片| 国产午夜精品久久久久久一区二区三区 | 欧美xxxx性猛交bbbb| 啦啦啦韩国在线观看视频| 成人特级黄色片久久久久久久| 草草在线视频免费看| 久久久久久久精品吃奶| 尤物成人国产欧美一区二区三区| 亚洲色图av天堂| 日本与韩国留学比较| 色哟哟·www| 成人特级av手机在线观看| 1024手机看黄色片| 特大巨黑吊av在线直播| 亚洲,欧美,日韩| 国产精品无大码| 听说在线观看完整版免费高清| 免费在线观看成人毛片| 中文字幕高清在线视频| 欧美+亚洲+日韩+国产| 变态另类丝袜制服| 成人av在线播放网站| 日日夜夜操网爽| 久久精品91蜜桃| 色吧在线观看| 一本久久中文字幕| 老师上课跳d突然被开到最大视频| 精品福利观看| 久久这里只有精品中国| 最近最新免费中文字幕在线| 亚洲久久久久久中文字幕| 欧美绝顶高潮抽搐喷水| 成人三级黄色视频| 国产高清有码在线观看视频| 校园春色视频在线观看| av天堂中文字幕网| 亚洲自偷自拍三级| 18禁黄网站禁片午夜丰满| 男女那种视频在线观看| 精品一区二区三区人妻视频| 一边摸一边抽搐一进一小说| 黄色欧美视频在线观看| 成人av一区二区三区在线看| 亚洲内射少妇av| 在现免费观看毛片| 99热精品在线国产| 国产精品亚洲美女久久久| 亚洲精品亚洲一区二区| 国产精品无大码| 18禁在线播放成人免费| 性欧美人与动物交配| 中国美白少妇内射xxxbb| 久久欧美精品欧美久久欧美| 九九热线精品视视频播放| 午夜激情福利司机影院| 不卡视频在线观看欧美| 男女下面进入的视频免费午夜| 99国产精品一区二区蜜桃av| 亚洲一区高清亚洲精品| 欧美日韩国产亚洲二区| 色噜噜av男人的天堂激情| 亚洲黑人精品在线| 欧美成人免费av一区二区三区| 国内精品一区二区在线观看| 久久中文看片网| 久久国内精品自在自线图片| 美女被艹到高潮喷水动态| 22中文网久久字幕| 日本爱情动作片www.在线观看 | 国产成人a区在线观看| 欧美日韩综合久久久久久 | 国产精品免费一区二区三区在线| 黄色视频,在线免费观看| 亚洲av成人精品一区久久| 麻豆久久精品国产亚洲av| 久久精品综合一区二区三区| 老司机福利观看| 网址你懂的国产日韩在线| 亚洲综合色惰| 国产人妻一区二区三区在| 成人综合一区亚洲| 国产高清视频在线播放一区| 欧美成人一区二区免费高清观看| 久久久精品大字幕| 亚洲精品在线观看二区| 一夜夜www| 精品日产1卡2卡| 日日摸夜夜添夜夜添av毛片 | 88av欧美| 亚洲成人久久爱视频| 国内精品久久久久精免费| 色噜噜av男人的天堂激情| 禁无遮挡网站| 麻豆一二三区av精品| 亚洲av中文av极速乱 | 人人妻人人看人人澡| 免费看a级黄色片| 国产高清有码在线观看视频| 中出人妻视频一区二区| 日韩欧美国产一区二区入口| 在线观看午夜福利视频| 校园春色视频在线观看| 女同久久另类99精品国产91| 神马国产精品三级电影在线观看| 欧美bdsm另类| 国产精品国产三级国产av玫瑰| 美女xxoo啪啪120秒动态图| 嫁个100分男人电影在线观看| 成年女人毛片免费观看观看9| 小说图片视频综合网站| 久久久久久国产a免费观看| 免费看av在线观看网站| 99视频精品全部免费 在线| 九色成人免费人妻av| 热99re8久久精品国产| 91在线观看av| 97人妻精品一区二区三区麻豆| 欧美成人性av电影在线观看| 精品久久久久久久久久免费视频| 亚洲 国产 在线| 我要搜黄色片| 日本三级黄在线观看| 国产免费一级a男人的天堂| 噜噜噜噜噜久久久久久91| 亚洲一区高清亚洲精品| 日韩中文字幕欧美一区二区| 国产一区二区亚洲精品在线观看| 午夜福利成人在线免费观看| 精品一区二区三区人妻视频| 在线免费观看的www视频| 亚洲av电影不卡..在线观看| 精品一区二区三区视频在线观看免费| a级毛片免费高清观看在线播放| 欧美3d第一页| 成人特级黄色片久久久久久久| 麻豆国产97在线/欧美| 美女 人体艺术 gogo| 91狼人影院| www日本黄色视频网| www.www免费av| 欧美最黄视频在线播放免费| 无遮挡黄片免费观看| 丰满乱子伦码专区| 内地一区二区视频在线| 国产精品一区二区性色av| 国模一区二区三区四区视频| 国产亚洲精品av在线| 男人舔女人下体高潮全视频| 久久国内精品自在自线图片| 久久久久久久久久成人| 亚洲av.av天堂| 国产成人aa在线观看| 日韩精品中文字幕看吧| 深夜精品福利| 亚洲欧美精品综合久久99| 免费av不卡在线播放| 国产精品爽爽va在线观看网站| 老司机福利观看| 国产午夜福利久久久久久| 免费看美女性在线毛片视频| 国产综合懂色| 丰满人妻一区二区三区视频av| 99精品久久久久人妻精品| 亚洲最大成人中文| 国产精品福利在线免费观看| 97热精品久久久久久| 中文字幕高清在线视频| 国产不卡一卡二| 麻豆成人午夜福利视频| 国产免费男女视频| 高清日韩中文字幕在线| 热99re8久久精品国产| 一个人看的www免费观看视频| 美女高潮喷水抽搐中文字幕| 蜜桃亚洲精品一区二区三区| 永久网站在线| 亚洲七黄色美女视频| 乱码一卡2卡4卡精品| 久久久午夜欧美精品| 成人无遮挡网站| av专区在线播放| av天堂在线播放| 两个人视频免费观看高清| 乱码一卡2卡4卡精品| 久久精品人妻少妇| 日韩av在线大香蕉| 免费在线观看日本一区| 日日撸夜夜添| x7x7x7水蜜桃| 免费观看在线日韩| 在线观看美女被高潮喷水网站| 国产 一区 欧美 日韩| 两个人的视频大全免费| 欧美在线一区亚洲| 小蜜桃在线观看免费完整版高清| 国产精品免费一区二区三区在线| 啦啦啦啦在线视频资源| 999久久久精品免费观看国产| 免费大片18禁| 最近最新中文字幕大全电影3| av中文乱码字幕在线| www.www免费av| 午夜日韩欧美国产| 日日撸夜夜添| 干丝袜人妻中文字幕| 老司机深夜福利视频在线观看| 国产乱人视频| 91在线精品国自产拍蜜月| 亚洲在线观看片| 欧美日本亚洲视频在线播放| 男人和女人高潮做爰伦理| 久久久久久久亚洲中文字幕| 欧美性感艳星| 一级av片app| 国产乱人视频| 1024手机看黄色片| 国产乱人视频| 久9热在线精品视频| 韩国av一区二区三区四区| 国产精品自产拍在线观看55亚洲| 可以在线观看毛片的网站| 免费高清视频大片| 久久精品国产亚洲av天美| 精品人妻1区二区| 久久99热6这里只有精品| 两个人的视频大全免费| 俺也久久电影网| 久久久久久久久久成人| 国内毛片毛片毛片毛片毛片| 国产激情偷乱视频一区二区| 免费电影在线观看免费观看| 欧洲精品卡2卡3卡4卡5卡区| 免费看a级黄色片| 一本久久中文字幕| 三级国产精品欧美在线观看| 91久久精品国产一区二区成人| 波多野结衣巨乳人妻| 他把我摸到了高潮在线观看| 日韩 亚洲 欧美在线| 97碰自拍视频| 久久久色成人| 少妇高潮的动态图| 中亚洲国语对白在线视频| 伊人久久精品亚洲午夜| 亚洲狠狠婷婷综合久久图片| 国产 一区 欧美 日韩| 久久亚洲精品不卡| 最新在线观看一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 国模一区二区三区四区视频| 亚洲av二区三区四区| 日本一二三区视频观看| 中国美白少妇内射xxxbb| 久久久久久伊人网av| 啦啦啦观看免费观看视频高清| 男女做爰动态图高潮gif福利片| 99热这里只有是精品在线观看| 老熟妇仑乱视频hdxx| 俺也久久电影网| 波多野结衣巨乳人妻| 精品午夜福利在线看| 成人av在线播放网站| 99久久精品一区二区三区| 日日啪夜夜撸| 日韩欧美免费精品| 91麻豆av在线| 亚洲最大成人中文| 美女大奶头视频| 一夜夜www| 中文字幕人妻熟人妻熟丝袜美| 我的老师免费观看完整版| 91久久精品国产一区二区三区| 国产v大片淫在线免费观看| 国产激情偷乱视频一区二区| 亚洲精品在线观看二区| 日本a在线网址| 亚洲久久久久久中文字幕| 欧美三级亚洲精品| 色综合亚洲欧美另类图片| 久久午夜亚洲精品久久| 一区二区三区激情视频| 18禁黄网站禁片免费观看直播| 99热6这里只有精品| 精品不卡国产一区二区三区| 日本爱情动作片www.在线观看 | 精品欧美国产一区二区三| 国产精品国产高清国产av| 亚洲狠狠婷婷综合久久图片| 韩国av在线不卡| ponron亚洲| 成年人黄色毛片网站| 女同久久另类99精品国产91| 久久久久久久精品吃奶| 少妇的逼水好多| 哪里可以看免费的av片| 啦啦啦韩国在线观看视频| 精华霜和精华液先用哪个| 精品一区二区三区视频在线| 中国美白少妇内射xxxbb| 18禁黄网站禁片免费观看直播| 久久99热6这里只有精品| 亚洲美女视频黄频| 国产av麻豆久久久久久久| 一夜夜www| 亚洲中文字幕日韩| 日本撒尿小便嘘嘘汇集6| 国产主播在线观看一区二区| 国产欧美日韩精品一区二区| 在线观看美女被高潮喷水网站| 在线a可以看的网站| 国产精品美女特级片免费视频播放器| 少妇人妻一区二区三区视频| 尤物成人国产欧美一区二区三区| 亚洲av日韩精品久久久久久密| 免费观看精品视频网站| 99热只有精品国产| 日韩欧美国产在线观看| 搡女人真爽免费视频火全软件 | 久久亚洲真实| 亚洲无线观看免费| 亚洲精品粉嫩美女一区| 最近最新中文字幕大全电影3| 丰满的人妻完整版| 国产69精品久久久久777片| 日韩精品青青久久久久久| av女优亚洲男人天堂| av福利片在线观看| 色播亚洲综合网| 欧美+日韩+精品| 久久久精品欧美日韩精品| 免费在线观看日本一区| 国产成年人精品一区二区| 国产在视频线在精品| 国产精品久久久久久精品电影| 校园春色视频在线观看| 午夜a级毛片| 蜜桃亚洲精品一区二区三区| 自拍偷自拍亚洲精品老妇| 无人区码免费观看不卡| 国内精品美女久久久久久| 亚洲精华国产精华精| 亚洲av二区三区四区| 国产精品久久久久久av不卡| 又黄又爽又免费观看的视频| 美女xxoo啪啪120秒动态图| 国产亚洲精品综合一区在线观看| 88av欧美| 男女那种视频在线观看| 精品乱码久久久久久99久播| 成人二区视频| 亚洲真实伦在线观看| 麻豆精品久久久久久蜜桃| 22中文网久久字幕| 午夜爱爱视频在线播放| 12—13女人毛片做爰片一| 精品欧美国产一区二区三| 亚洲中文字幕一区二区三区有码在线看| 亚洲成a人片在线一区二区| 久久精品国产亚洲av香蕉五月| 99久国产av精品| 久久精品国产自在天天线| 久久99热6这里只有精品| 国产一区二区在线观看日韩| 亚洲三级黄色毛片| 欧美性猛交黑人性爽| 日韩精品中文字幕看吧| 欧美成人一区二区免费高清观看| 国产亚洲精品久久久久久毛片| 亚洲aⅴ乱码一区二区在线播放| 国内少妇人妻偷人精品xxx网站| 国产成人aa在线观看| 男人舔奶头视频| 好男人在线观看高清免费视频| 精品久久久久久久久亚洲 | 午夜福利高清视频| 毛片女人毛片| 黄色日韩在线| 亚洲va在线va天堂va国产| 国产视频一区二区在线看| 一边摸一边抽搐一进一小说| 2021天堂中文幕一二区在线观| 亚洲成a人片在线一区二区| 国产一区二区三区av在线 | 丰满的人妻完整版| 又黄又爽又刺激的免费视频.| 亚洲熟妇中文字幕五十中出| 日韩一区二区视频免费看| 欧美精品国产亚洲| 国内精品一区二区在线观看| 中国美女看黄片| 九九在线视频观看精品| 国产亚洲av嫩草精品影院| 91午夜精品亚洲一区二区三区 | 婷婷精品国产亚洲av| 日本免费a在线| 日本黄大片高清| 国产麻豆成人av免费视频| 波多野结衣高清无吗| 亚洲avbb在线观看| 精品乱码久久久久久99久播| 大又大粗又爽又黄少妇毛片口| 国产不卡一卡二| 午夜日韩欧美国产| 日本免费a在线| av在线亚洲专区| 精品一区二区三区视频在线| 又粗又爽又猛毛片免费看| 国产美女午夜福利| 我的女老师完整版在线观看| 亚洲三级黄色毛片| 国产成年人精品一区二区| 中文字幕av在线有码专区| 国产精品免费一区二区三区在线| 欧美精品啪啪一区二区三区| 亚洲黑人精品在线| 亚洲精品一卡2卡三卡4卡5卡| 麻豆av噜噜一区二区三区| 美女大奶头视频| 亚洲av电影不卡..在线观看| 国产一区二区三区在线臀色熟女| 日韩国内少妇激情av| 我要看日韩黄色一级片| 亚洲电影在线观看av| 可以在线观看的亚洲视频| 亚洲中文字幕一区二区三区有码在线看| 美女高潮的动态| 午夜a级毛片| 欧美高清性xxxxhd video| 97超视频在线观看视频| 久久久久久久久久久丰满 | 亚洲精品色激情综合| 男人狂女人下面高潮的视频| 亚洲精品粉嫩美女一区| 久久久久免费精品人妻一区二区|