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

    荒漠草原物種多樣性、生產(chǎn)力與氣候因子和土壤養(yǎng)分之間關(guān)系的研究

    2015-06-01 12:15:17孫小麗康薩如拉張慶常昌明馬文靜牛建明
    草業(yè)學(xué)報(bào) 2015年12期
    關(guān)鍵詞:荒漠生產(chǎn)力生物量

    孫小麗,康薩如拉,張慶,2,常昌明,馬文靜,牛建明,2*

    (1.內(nèi)蒙古大學(xué)生命科學(xué)學(xué)院,內(nèi)蒙古 呼和浩特 010021;2.中美生態(tài)、能源及可持續(xù)性科學(xué)研究中心,內(nèi)蒙古 呼和浩特 010021)

    ?

    荒漠草原物種多樣性、生產(chǎn)力與氣候因子和土壤養(yǎng)分之間關(guān)系的研究

    孫小麗1,康薩如拉1,張慶1,2,常昌明1,馬文靜1,牛建明1,2*

    (1.內(nèi)蒙古大學(xué)生命科學(xué)學(xué)院,內(nèi)蒙古 呼和浩特 010021;2.中美生態(tài)、能源及可持續(xù)性科學(xué)研究中心,內(nèi)蒙古 呼和浩特 010021)

    物種多樣性對(duì)生態(tài)系統(tǒng)功能的影響是生態(tài)學(xué)領(lǐng)域的核心問(wèn)題之一。本研究利用內(nèi)蒙古荒漠草原80個(gè)野外植被、土壤調(diào)查樣地資料,結(jié)合氣候數(shù)據(jù),探討物種多樣性、生產(chǎn)力與氣候因子和土壤養(yǎng)分的關(guān)系,為進(jìn)一步認(rèn)識(shí)多樣性-生產(chǎn)力關(guān)系提供參考。結(jié)果表明,1)物種多樣性、生物量與年平均降水量均呈線性正相關(guān),與年平均溫度均呈線性負(fù)相關(guān),與干燥指數(shù)均呈線性正相關(guān)。2)物種豐富度與土壤全磷含量呈線性負(fù)相關(guān)(R2=0.072,P<0.05),與速效氮含量呈線性正相關(guān)(R2=0.183,P<0.01),Shannon-Wiener指數(shù)與速效氮也呈線性正相關(guān)(R2=0.061,P<0.05);生物量與土壤全氮、速效氮和有機(jī)碳含量都呈線性正相關(guān),相關(guān)系數(shù)分別為0.150,0.177,0.089(P<0.01)。3)群落生物量與物種多樣性呈線性正相關(guān)關(guān)系。4)干燥指數(shù)是影響多樣性和生物量的主要環(huán)境因子。

    荒漠草原;生產(chǎn)力;物種多樣性;土壤養(yǎng)分;水熱因子

    物種多樣性對(duì)生態(tài)系統(tǒng)功能的影響是生態(tài)學(xué)領(lǐng)域的核心科學(xué)問(wèn)題之一[1-2]。生產(chǎn)力水平的高低是生態(tài)系統(tǒng)功能的重要表現(xiàn)形式,所以探討物種多樣性與生產(chǎn)力之間的關(guān)系對(duì)于研究物種多樣性對(duì)生態(tài)系統(tǒng)功能的影響具有重要意義[3-4]。

    物種多樣性不僅能夠度量群落的組成結(jié)構(gòu)和功能的復(fù)雜性,而且也能指示環(huán)境[5]。近年來(lái),圍繞生物多樣性、生產(chǎn)力與環(huán)境因子的關(guān)系,國(guó)內(nèi)外開展了一系列討論[6-9],其中較多的是關(guān)于物種多樣性、生產(chǎn)力與緯度和海拔導(dǎo)致的光照和水熱因子變化的關(guān)系[6-7]。尤其是在干旱、半干旱地區(qū),多樣性和生產(chǎn)力極大地受到降水條件的影響[10]。植物的生長(zhǎng)發(fā)育與土壤因子的關(guān)系也十分密切,土壤理化性質(zhì)的不同以及土壤母質(zhì)條件等,都可能影響群落結(jié)構(gòu)、物種多樣性和生產(chǎn)力水平[11-12]。關(guān)于土壤養(yǎng)分與物種多樣性的關(guān)系已有一些報(bào)道,然而由于土壤養(yǎng)分對(duì)多樣性的影響復(fù)雜多樣,研究結(jié)果缺乏一致的規(guī)律性。目前關(guān)于多樣性、生產(chǎn)力與土壤養(yǎng)分關(guān)系的研究多涉及森林[13-14]、草原[15-16]、高寒草甸[17-18]等,但在本研究區(qū),荒漠草原的理論探討還不夠全面。

    內(nèi)蒙古荒漠草原是草原向荒漠過(guò)渡的旱生性最強(qiáng)的草原生態(tài)系統(tǒng)[19]。前人在該地區(qū)積累了大量的研究成果及文獻(xiàn)資料,如張慶等[20]對(duì)群落分布格局與環(huán)境因子關(guān)系的探討,楊雪棟等[21]對(duì)小針茅(Stipaklemenzii)草原土壤養(yǎng)分化學(xué)計(jì)量學(xué)特征的研究,馬文靜等[22]闡述了短花針茅(Stipabreviflora)草原物種多樣性、功能多樣性與生產(chǎn)力之間的關(guān)系。但是該地區(qū)尚缺乏聯(lián)系氣候因子和土壤養(yǎng)分等環(huán)境因子的多樣性、生產(chǎn)力研究。本研究通過(guò)對(duì)內(nèi)蒙古荒漠草原進(jìn)行植被與土壤調(diào)查,結(jié)合氣候資料,在群落水平上探討內(nèi)蒙古荒漠草原物種多樣性、生產(chǎn)力及其與氣候因子和土壤養(yǎng)分的關(guān)系,為進(jìn)一步認(rèn)識(shí)物種多樣性-生產(chǎn)力關(guān)系提供參考。

    1 材料與方法

    1.1研究區(qū)概況

    研究區(qū)位于內(nèi)蒙古中溫型荒漠草原區(qū),地處107°18′-115°35′ E,41°11′-45°5′ N,海拔905~1725 m。該地區(qū)年平均降水量為150~250 mm,多集中在夏季;年平均溫度2~5℃,≥10℃年積溫2200~2500℃。地帶性土壤為棕鈣土[19,23],主要植被類型為小針茅荒漠草原和短花針茅荒漠草原。

    圖1 樣點(diǎn)分布圖

    1.2數(shù)據(jù)獲取

    植被調(diào)查選擇在群落生物量高峰期進(jìn)行,于2012年7月下旬至8月中旬完成。共設(shè)置了80個(gè)具有代表性的樣地(圖1)。植被調(diào)查采用的是樣方法,每個(gè)樣地在10 m×10 m的范圍內(nèi),隨機(jī)設(shè)置3個(gè)1 m×1 m測(cè)產(chǎn)樣方,記錄樣方中出現(xiàn)的每一物種名、株叢數(shù),并按物種將植株地上部分齊地面刈割,稱取鮮重,然后分別用報(bào)紙包好保存,帶回實(shí)驗(yàn)室置于65℃烘箱烘干至恒重,稱取干重。另做7個(gè)頻度樣方,記錄樣方中出現(xiàn)的物種名。同時(shí),用土鉆鉆取0~30 cm混合土樣,3次重復(fù)。帶回實(shí)驗(yàn)室風(fēng)干后過(guò)1 mm和0.15 mm篩后用于測(cè)定土壤基本指標(biāo)。測(cè)試項(xiàng)目為土壤的全氮(total nitrogen,TN)(硒粉、硫酸銅、硫酸鉀消化-蒸餾法[24])、全磷(total phosphorus,TP)(堿熔-鉬銻抗比色法[24])、速效氮(available nitrogen,AN)(凱氏定氮法[24])、速效磷(available phosphorus,AP)(碳酸氫鈉浸提-鉬銻抗比色法[24])和有機(jī)碳(organic carbon,OC)(重鉻酸鉀氧化外加熱法[24])。群落生產(chǎn)力(productivity)用3個(gè)測(cè)產(chǎn)樣方地上生物量干重的平均值代替。群落的物種豐富度由10個(gè)樣方中出現(xiàn)的物種數(shù)決定。

    年平均降水量(mean annual precipitation, MAP)和年平均溫度(mean annual temperature, MAT)數(shù)據(jù)采用WORLDCLIM網(wǎng)站(http://www.worldclim.org/)下載的當(dāng)前(1950-2000年)監(jiān)測(cè)數(shù)據(jù),數(shù)據(jù)空間分辨率為arc 30 s。干燥指數(shù)(aridity index, AI)數(shù)據(jù)源于CGIAR-CSI GeoPortal網(wǎng)站(http://www.csi.cgiar.org)的數(shù)據(jù)產(chǎn)品,這一產(chǎn)品是根據(jù)WORLDCLIM網(wǎng)站的數(shù)據(jù)建模生成的[25]:

    AI=MAP/MAE

    式中,AI為干燥指數(shù),MAP為年平均降水量,MAE(mean annual potential evapo-transpiration)為年平均潛在蒸散量。AI值越大代表氣候條件越濕潤(rùn),AI值越小代表氣候條件越干旱。

    1.3植物群落物種多樣性的指標(biāo)

    物種豐富度指數(shù):R=S

    式中,R為群落物種豐富度指數(shù),S為群落中的總物種數(shù)。

    式中,H′為Shannon-Wiener指數(shù),Pi為種i的相對(duì)重要值=(相對(duì)生物量+相對(duì)多度+相對(duì)頻度)/3,S為種i所在樣地所有物種數(shù)。

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

    采用線性回歸方法分析物種多樣性、生物量與氣候因子和土壤養(yǎng)分之間的關(guān)系以及物種多樣性與生物量之間的關(guān)系。通過(guò)多元逐步回歸分析確定氣候因子及土壤因子對(duì)群落豐富度、多樣性和生物量變異的貢獻(xiàn)。將引入變量的P值定為0.05,剔除變量的P值定為0.10,判定系數(shù)(R2)揭示了回歸模型中自變量能夠解釋的因變量變異比例。

    使用SPSS 21.0軟件進(jìn)行統(tǒng)計(jì)分析,采用SigmaPlot 12.5繪散點(diǎn)圖制。

    2 結(jié)果與分析

    2.1物種多樣性、生物量與氣候因子的關(guān)系

    根據(jù)研究區(qū)80個(gè)樣地的物種豐富度、Shannon-Wiener指數(shù)和生物量與氣候因子建立的回歸模型(圖2)可以看出,物種豐富度、Shannon-Wiener指數(shù)和生物量均隨著年平均降水量的增加而顯著上升(圖2A,B,C),其中物種豐富度和生物量與年平均降水量呈極顯著的線性正相關(guān);Shannon-Wiener指數(shù)和生物量與年平均溫度均呈顯著的線性負(fù)相關(guān)(圖2E,F(xiàn)),物種豐富度與年平均溫度相關(guān)性不顯著(圖2D);物種豐富度、Shannon-Wiener指數(shù)和生物量均隨著干燥指數(shù)的增大呈顯著的上升趨勢(shì)(圖2G,H,I)。

    圖2 物種豐富度、多樣性和生物量與水熱因子和干燥指數(shù)的關(guān)系

    2.2物種多樣性、生物量與土壤養(yǎng)分的關(guān)系

    由表1可知,物種豐富度與土壤全磷含量呈線性負(fù)相關(guān)(P<0.05),物種豐富度、Shannon-Wiener指數(shù)與土壤速效氮含量具有顯著的線性正相關(guān)關(guān)系(P<0.01、P<0.05),生物量與土壤全氮、速效氮以及有機(jī)碳含量均呈極顯著的正相關(guān)(P<0.01)。其他土壤養(yǎng)分與物種豐富度、Shannon-Wiener指數(shù)和生物量關(guān)系不顯著。

    表1 物種豐富度、多樣性和生物量與土壤養(yǎng)分含量的關(guān)系Table 1 The regression analyses of species richness, species diversity, and biomass with soil nutrient content

    NS 代表關(guān)系不顯著;PL代表正線性關(guān)系;NL 代表負(fù)線性關(guān)系。

    NS indicates non-significant relationship; PL indicates positive linear relationship; NL indicates negative linear relationship.

    2.3物種多樣性與生物量的關(guān)系

    群落生物量與物種豐富度和Shannon-Wiener指數(shù)呈顯著的線性正相關(guān),相關(guān)性系數(shù)分別為0.076(P<0.05)和0.074(P<0.05),即物種多樣性對(duì)生產(chǎn)力存在正效應(yīng)(圖3)。

    圖3 物種豐富度、多樣性與生物量的關(guān)系

    2.4影響物種多樣性、生物量的主導(dǎo)環(huán)境因子

    通過(guò)多元逐步回歸分析環(huán)境因子對(duì)物種豐富度的影響發(fā)現(xiàn),該模型引入了AI、TP和AN三個(gè)變量,剔除了其他變量,且這3個(gè)因子與物種豐富度均存在顯著相關(guān)關(guān)系,解釋了物種豐富度的30.9%,其中AI作為第一解釋變量進(jìn)入方程,解釋了總變異的21.8%,成為影響群落物種豐富度最主要的環(huán)境因子(表2)。另外,MAT作為最重要的因子被引入解釋Shannon-Wiener指數(shù)的方程中,解釋率為9.9%,TP和AN同樣作為第二和第三因子被引入,說(shuō)明MAT是對(duì)Shannon-Wiener指數(shù)影響最大的環(huán)境因子。環(huán)境因子對(duì)生物量影響的多元逐步回歸分析,AI作為唯一因子被引入方程,解釋了群落生產(chǎn)力的27.2%。

    2.5干燥指數(shù)對(duì)多樣性-生產(chǎn)力關(guān)系的影響

    為進(jìn)一步探討干燥指數(shù)對(duì)物種多樣性-生產(chǎn)力關(guān)系的影響,按聯(lián)合國(guó)環(huán)境規(guī)劃署基于干燥指數(shù)的全球氣候分類標(biāo)準(zhǔn)[26]將研究區(qū)劃分為干旱、半干旱兩種干濕氣候類型。研究區(qū)80個(gè)樣地,其中25個(gè)分布在干旱區(qū),55個(gè)分布在半干旱區(qū)。利用線性回歸方法分析兩種不同氣候類型(干旱區(qū)、半干旱區(qū))下多樣性-生產(chǎn)力的關(guān)系,依據(jù)R2分析干燥指數(shù)對(duì)兩者關(guān)系的影響(圖4)。發(fā)現(xiàn)無(wú)論物種豐富度還是Shannon-Wiener指數(shù),當(dāng)干燥指數(shù)高時(shí)(半干旱區(qū))對(duì)生產(chǎn)力的貢獻(xiàn),均比干燥指數(shù)低時(shí)(干旱區(qū))大。

    3 結(jié)論與討論

    3.1水熱因子對(duì)物種多樣性、生產(chǎn)力的重要性

    表2 多元逐步回歸分析結(jié)果Table 2 Results of stepwise regression analysis

    AI:干燥指數(shù) Aridity Index; TP:全磷 Total phosphorus;AN:速效氮 Available nitrogen;MAT:年平均溫度 Mean annual temperature.

    內(nèi)蒙古草原處于干旱半干旱區(qū),其生物多樣性和生產(chǎn)力深受降水條件的限制,大量研究表明,隨著降水量的增加,生物多樣性和生產(chǎn)力呈顯著增加的趨勢(shì)[27-28]。有學(xué)者對(duì)美國(guó)[29]、委內(nèi)瑞拉[30]、阿根廷[31]干旱區(qū)草原生態(tài)系統(tǒng)的研究均證明物種豐富度與降水量呈顯著的正相關(guān)關(guān)系。本研究中物種多樣性和生產(chǎn)力與年平均降水量呈顯著正相關(guān),與年平均溫度呈負(fù)相關(guān),與干燥指數(shù)呈顯著正相關(guān)。物種豐富度、Shannon-Wiener指數(shù)和生產(chǎn)力與水熱因子相關(guān)性的比較結(jié)果顯示,水熱因子與群落生產(chǎn)力的相關(guān)性最高,可以解釋其變異的13.0%~28.1%,降水與物種豐富度的相關(guān)性高于Shannon-Wiener指數(shù),可以解釋其變異的22.1%,3個(gè)氣候要素中, 作為水熱綜合因子的干燥指數(shù)同物種豐富度和初級(jí)生產(chǎn)力的相關(guān)性最高。白永飛等[8]對(duì)錫林河流域植物群落的研究也表明物種豐富度、多樣性和生產(chǎn)力與年降水量正相關(guān),與年平均氣溫負(fù)相關(guān),而水熱因子可以解釋生產(chǎn)力的變異高達(dá)71.31%~77.58%,可以解釋物種豐富度變異的40.2%~61.39%;Kang等[28]對(duì)整個(gè)內(nèi)蒙古草原的研究發(fā)現(xiàn)地上生物量與年平均降水量的相關(guān)系數(shù)高達(dá)65%;然而馬文靜等[22]對(duì)短花針茅荒漠草原的研究顯示年平均降水量對(duì)物種豐富度和生產(chǎn)力的解釋率也不高,分別為12.6%和10.5%。造成荒漠草原這種解釋率偏低的原因可能與其獨(dú)特的生境特征有關(guān),該區(qū)年平均降水量普遍偏低,導(dǎo)致研究區(qū)內(nèi)降水梯度不明顯;此外土壤、放牧等因素也可能起重要作用。在本研究區(qū)內(nèi),存在著隨年平均溫度升高,年平均降水量顯著降低的特點(diǎn),二者呈極顯著的負(fù)相關(guān),因此,物種多樣性和生產(chǎn)力與年平均溫度呈負(fù)相關(guān),極有可能是由于水分的限制,導(dǎo)致植物無(wú)法利用這些能量進(jìn)行光合作用,從而呈隨著能量增加群落物種多樣性和生產(chǎn)力降低的關(guān)系。

    3.2土壤養(yǎng)分對(duì)物種多樣性、生產(chǎn)力的重要性

    在土壤-植被體系中,土壤和植被是兩個(gè)相互依存的因子,植被影響土壤,土壤制約植被[32]。探討土壤養(yǎng)分與物種多樣性的關(guān)系,對(duì)了解生態(tài)系統(tǒng)的功能與過(guò)程具有重要意義。目前二者之間關(guān)系的研究結(jié)果缺乏一致的規(guī)律性。多數(shù)氮素添加實(shí)驗(yàn)研究發(fā)現(xiàn),增加氮素會(huì)導(dǎo)致物種多樣性下降[33-34],Perroni-Ventura等[12]對(duì)熱帶半干旱生態(tài)系統(tǒng)的研究顯示物種豐富度隨著土壤氮含量的增加而升高。Olsen等[35]則認(rèn)為在貧瘠的土壤里增加氮含量有利于提高物種豐富度,在養(yǎng)分豐富的土壤增加氮含量會(huì)使物種豐富度降低。究其原因,對(duì)退化的、資源貧瘠的群落添加氮素可以減輕氮素的限制,從而提高物種豐富度;相反,對(duì)相對(duì)成熟穩(wěn)定的群落添加氮素可能會(huì)破壞原有生態(tài)系統(tǒng)的共存機(jī)制,導(dǎo)致物種豐富度下降[36]。就磷元素而言,張林靜等[37]認(rèn)為土壤全磷含量與物種多樣性沒(méi)有顯著相關(guān)性,Jouany等[38]對(duì)歐洲溫性草原和Perroni等[39]對(duì)墨西哥的草原-荒漠疏林生態(tài)系統(tǒng)研究均表明物種豐富度隨著土壤中全磷含量的升高而降低。本研究發(fā)現(xiàn)在荒漠草原這一生態(tài)系統(tǒng)類型中,物種多樣性沿著速效氮含量增加的梯度呈升高的趨勢(shì),沿著土壤全磷含量增加的梯度呈降低的趨勢(shì)。地上生物量與有機(jī)碳、全氮和速效氮含量呈顯著的線性正相關(guān),這與左小安等[5]對(duì)科爾沁沙地的研究是一致的。土壤中的養(yǎng)分,在空間和時(shí)間上的分布是具有異質(zhì)性的,作為對(duì)土壤養(yǎng)分異質(zhì)性的反應(yīng),植物在養(yǎng)分豐富的局部環(huán)境中能選擇性地改變其根系的生長(zhǎng),從而增加養(yǎng)分的吸收[40]。土壤中養(yǎng)分含量的高低直接影響群落的生產(chǎn)力,土壤養(yǎng)分越豐富,群落生產(chǎn)力越高[41]。由表1可知,速效氮與物種多樣性和生產(chǎn)力的相關(guān)系數(shù)均最高,這可能是因?yàn)椴莸厣鷳B(tài)系統(tǒng)普遍缺氮,物種多樣性和生產(chǎn)力對(duì)氮含量的增加反應(yīng)比較敏感,這在一定程度上支持了Vitousek等[42]的觀點(diǎn),即氮素是草地生態(tài)系統(tǒng)凈初級(jí)生產(chǎn)力的限制因子,同時(shí)也是衡量草地生產(chǎn)力的重要指標(biāo)。Tilman等[43]認(rèn)為高的植物多樣性能夠更完全地利用土壤中的限制性養(yǎng)分,減少生態(tài)系統(tǒng)中N素的淋溶損失。同時(shí),多樣性高的生態(tài)系統(tǒng)能夠保持更高比例的限制性養(yǎng)分,使土壤養(yǎng)分的總儲(chǔ)量增加,從而使生產(chǎn)力增加。有學(xué)者指出氮素對(duì)草地生產(chǎn)力的有效性受到水分的限制,尤其在干旱區(qū),水分的缺乏會(huì)導(dǎo)致氮素的有效性偏低[44]。因此,降水的缺乏極有可能是導(dǎo)致荒漠草原氮素對(duì)生產(chǎn)力貢獻(xiàn)(R2=0.177)低于典型草原和草甸草原(R2=0.65)[8]的主要原因。

    3.3水熱因子是影響物種多樣性和生產(chǎn)力的主要環(huán)境因子

    為進(jìn)一步確定影響荒漠草原物種多樣性和生產(chǎn)力的主導(dǎo)環(huán)境因子,建立了多元逐步回歸模型,發(fā)現(xiàn)干燥指數(shù)作為解釋物種豐富度的第一因子和解釋群落生產(chǎn)力的唯一最重要的因子被引入,而且都呈顯著的相關(guān)關(guān)系,年平均溫度作為解釋Shannon-Wiener指數(shù)的第一因子被引入模型。干燥指數(shù)是水分與溫度綜合作用的體現(xiàn),這充分說(shuō)明了水熱因子在內(nèi)蒙古荒漠草原區(qū)對(duì)物種多樣性和生產(chǎn)力作用的主導(dǎo)地位。白永飛等[8]發(fā)現(xiàn)作為水熱綜合因子的干燥度與錫林河流域植物群落多樣性和生產(chǎn)力的相關(guān)系數(shù)是最高的;鄭曉翾等[45]對(duì)呼倫貝爾草甸草原、典型草原和荒漠草原3種草原類型的研究也證明水熱因子是影響物種豐富度的主要環(huán)境因子。但是同時(shí)也發(fā)現(xiàn),這種線性關(guān)系的相關(guān)系數(shù)并不高,這可能與研究尺度和放牧利用等因素有關(guān)。在大的空間尺度下,氣候因子的梯度變化更為突出,稀釋了局地尺度人類活動(dòng)的作用,因此氣候?qū)ξ锓N多樣性分布格局解釋量相對(duì)較高,在小的空間尺度下,氣候因子的梯度變化不明顯,氣候解釋量反而下降[46-47]。此外,放牧是荒漠草原區(qū)的主要利用方式,牲畜對(duì)牧草的采食和土壤的踐踏以及牲畜體內(nèi)營(yíng)養(yǎng)物質(zhì)的歸還和轉(zhuǎn)化,對(duì)草地的微環(huán)境和土壤條件會(huì)造成顯著影響[48],進(jìn)而影響植物群落物種組成[49]、豐富度[50]、多樣性[51]等。本研究區(qū)內(nèi)群落結(jié)構(gòu)較為簡(jiǎn)單,更易受到人類活動(dòng)的干擾。目前,人類活動(dòng)很可能是影響內(nèi)蒙古荒漠草原多樣性和生產(chǎn)力格局的一個(gè)重要因素。

    3.4物種多樣性與生產(chǎn)力之間的關(guān)系受到環(huán)境因子的制約

    許多研究證實(shí)物種多樣性與生產(chǎn)力主要存在4種關(guān)系:正相關(guān)[43];負(fù)相關(guān)[52];單峰曲線[53];不相關(guān)[54]。其中呈單峰曲線和正相關(guān)關(guān)系的結(jié)論居多。就荒漠草原而言,物種多樣性與生產(chǎn)力之間呈線性正相關(guān),這一結(jié)論與馬文靜等[22]對(duì)內(nèi)蒙古短花針茅荒漠草原的研究相一致。在不同尺度下,物種多樣性和生產(chǎn)力的關(guān)系存在差異。在單一的群落類型內(nèi),多樣性與生產(chǎn)力之間的關(guān)系常常并不顯著[55],而當(dāng)研究區(qū)包括幾個(gè)群落類型的時(shí)候,二者之間的關(guān)系常常表現(xiàn)出單峰曲線的形式[56]。Ni等[57]在不同尺度上對(duì)內(nèi)蒙古不同類型草原物種多樣性與生產(chǎn)力關(guān)系的研究結(jié)果也發(fā)現(xiàn)在局地尺度上,物種豐富度與生物量呈線性正相關(guān),隨著尺度的增大,物種豐富度與生產(chǎn)力的關(guān)系可能會(huì)出現(xiàn)單峰曲線關(guān)系。

    干燥指數(shù)、年平均溫度、全磷和速效氮進(jìn)入了解釋物種豐富度、Shannon-wiener指數(shù)的方程中,表明水熱因子和土壤養(yǎng)分對(duì)物種多樣性有重要影響,進(jìn)而推測(cè),荒漠草原多樣性-生產(chǎn)力關(guān)系同時(shí)受到上述主要?dú)夂蛞蜃优c土壤養(yǎng)分的調(diào)控。探討干燥指數(shù)對(duì)物種多樣性-生產(chǎn)力關(guān)系的影響,發(fā)現(xiàn)無(wú)論物種豐富度還是Shannon-Wiener指數(shù),當(dāng)干燥指數(shù)高時(shí)(半干旱區(qū))對(duì)生產(chǎn)力的貢獻(xiàn),均比干燥指數(shù)低時(shí)(干旱區(qū))大(圖4)。He等[58]研究了在土壤養(yǎng)分N、P、K含量存在差異的情況下, 物種豐富度和生產(chǎn)力之間的關(guān)系,結(jié)果表明,在土壤N、P、K含量高的群落中,物種豐富度和生產(chǎn)力呈顯著的正相關(guān)關(guān)系,而在土壤養(yǎng)分N、P、K含量低的處理中, 兩者之間沒(méi)有明顯的關(guān)系。出現(xiàn)上述現(xiàn)象的原因可能是:植物在資源利用上具有互補(bǔ)性,不同的植物具有不同的資源利用形式,資源供給水平高可以使得物種間對(duì)資源的利用差異充分體現(xiàn)出來(lái),多樣性的作用表現(xiàn)更加充分,這樣就加強(qiáng)了多樣性與生產(chǎn)力之間正相關(guān)關(guān)系的表現(xiàn)。

    [1] Loreau M, Naeem S, Inchausti P,etal.Biodiversity and ecosystem functioning:Current knowledge and future challenges.Science, 2001, 294:804-808.

    [2] Hooper D U, Chapin F S III, Ewel J J,etal.Effects of biodiversity on ecosystem functioning:A consensus of current knowledge.Ecological Monographs, 2005, 75(1):3-35.

    [3] Naeem S.Ecosystem consequences of biodiversity loss:The evolution of a paradigm.Ecology, 2002, 83(6):1537-1552.

    [4] Tilman D, Reich P B, Knops J,etal.Diversity and productivity in a long-term grassland experiment.Science, 2001, 294:843-845.

    [5] Zuo X A, Zhao X Y, Zhao H L,etal.Changes of species diversity and productivity in relation to soil properties in sandy grassland in Horqin sand land.Environmental Science, 2007, 28(5):945-951.

    [6] Wilson J B, Lee W G, Mark A F.Species diversity in relation to ultramafic substrate and to altitude in southwestern New Zealand.Vegetatio, 1990, 86:15-20.

    [7] Zhang Q, Hou X Y, Li F Y H,etal.Alpha, beta and gamma diversity differ in response to precipitation in the Inner Mongolia Grassland.PLoS ONE, 9(3):e93518.doi:10.1371/journal.pone.0093518.

    [8] Bai Y F, Li L H, Wang Q B,etal.Changes in plant species diversity and productivity along gradients of precipitation and elevation in the Xilin River Basin, Inner Mongolia.Acta Phytoecologica Sinica, 2000, 24(6):667-673.

    [9] Zhang Z Z, Wu G L, Wang D,etal.Plant community structure and soil moisture in the semi-arid natural grassland of the Loess Plateau.Acta Prataculturae Sinica, 2014, 23(6):313-319.

    [10] Abramsky Z, Rosenzweig M L.Tilman’s predicted productivity-diversity relationship shown by desert rodents.Nature, 1984, 309:150-151.

    [11] Long W X, Yang X B, Li D H.Patterns of species diversity and soil nutrients along a chronosequence of vegetation recovery in Hainan Island, South China.Ecological Research, 2012, 27(3):561-568.

    [12] Perroni-Ventura Y, Montana C, Garcia-Oliva F.Relationship between soil nutrient availability and plant species richness in a tropical semi-arid environment.Journal of Vegetation Science, 2006, 17(6):719-728.

    [13] An S Q, Wang Z F, Zhu X L,etal.Effect of soil factors on species diversity in secondary forest communities.Journal of Wuhan Botanical Research, 1997, 15(2):143-150.

    [14] Firn J, Erskine P D, Lamb D.Woody species diversity influences productivity and soil nutrient availability in tropical plantations.Oecologia, 2007, 154(3):521-533.

    [15] Kahmen A, Perner J, Audorff V,etal.Effects of plant diversity, community composition and environmental parameters on productivity in montane european grasslands.Oecologia, 2005, 142(4):606-615.

    [16] Ceulemans T, Stevens C J, Duchateau L,etal.Soil phosphorus constrains biodiversity across european grasslands.Global Change Biology, 2014, 20(12):3814-3822.

    [17] Liu M X, Wang G.Responses of plant community diversity and soil factors to slope in alpine meadow.Chinese Journal of Ecology, 2013, 32(2):259-265.

    [18] Wang C T, Long R J, Cao G M,etal.The relationship between soil nutrients and diversity-productivity of different type grasslands in alpine meadow.Chinese Journal of Soil Science, 2008, 39(1):1-8.

    [19] Inner Mongolia-Ningxia jiont inspection group of Chinese Academy of Sciences.Vegetation of Inner Mongolia[M].Beijing:Science Press, 1985.

    [20] Zhang Q, Niu J M, Buyantuyev A,etal.Ecological analysis and classification ofStipabrevifloracommunities in the Inner Mongolia region:the role of environmental factors.Acta Prataculturae Sinica, 2012, 21(1):83-92.

    [21] Yang X D, Sun W G, Baoyin T.Characteristics of soil stoichiometry ofStipaklemeziicommunity in desert steppe of Inner Mongolia.Chinese Journal of Grassland, 2012, 34(5):30-34.

    [22] Ma W J, Zhang Q, Niu J M,etal.Relationship of ecosystem primary productivity to species diversity and functional group diversity:evidence fromStipabrevifloragrassland in Nei Mongol.Chinese Journal of Plant Ecology, 2013, 37(7):620-630.

    [23] Wei Z J, Han G D, Zhao G,etal.Research on Desert Grassland Ecosystem in China[M].Beijing:Science Press, 2013.

    [24] Zhang G L, Gong Z T.Soil Survey Laboratory Methods[M].Beijing:Science Press, 2012.

    [25] Trabucco A, Zomer R J, Bossio D A,etal.Climate change mitigation through afforestation/reforestation:A global analysis of hydrologic impacts with four case studies.Agriculture, Ecosystems and Environment, 2008, 126:81-97.

    [26] Middleton N J, Thomas D.World Atlas of Desertification[M].London:Arnold, United Nations Environment Programme, 1997.

    [27] Han B, Fan J W, Zhong H P.Grassland and biomass of communities along gradients of the Inner Mongolia grassland transect.Journal of Plant Ecology, 2006, 30(4):553-562.

    [28] Kang M, Dai C, Ji W,etal.Biomass and its allocation in relation to temperature, precipitation, and soil nutrients in Inner Mongolia Grasslands, China.PLoS ONE, 8(7):e69561.Doi:10.1371/journal.pone.0069561.

    [29] Daubenmire R.Steppe vegetation of Washington.Washington Agricultural Experiment Station Technical Bulletin, 1970, 62:1-131.

    [30] Sarmiento G.Patterns of specific and phenological diversity in the grass community of the venezuelan tropical savannas.Journal of Biogeography, 1983, 10:373-391.

    [31] Jobbagy E G, Paruelo J M, Leon R J C.Vegetation heterogeneity and diversity in flat and mountain landscapes ofPatagonia(Argentina).Journal of Vegetation Science, 1996, 7(4):599-608.

    [32] Tilman D.Causes, consequences and ethics of biodiversity.Nature, 2000, 405:208-211.

    [33] Zhang J Q, Li Q, Ren Z W,etal.Effect of nitrogen addition on species richness and relationship between species and aboveground productivity of alpine meadow of Qinghai-Tibetan Plateau, China.Chinese Journal of Plant Ecology, 2010, 34(10):1125-1131.

    [34] Cornwell W K, Grubb P J.Regional and local patterns in plant species richness with respect to resource availability.Oikos, 2003, 100(3):417-428.

    [35] Olsen S L, Sandvik S M, Totland O.Influence of two N-fixing legumes on plant community properties and soil nutrient levels in an Alpine.Arctic Antarctic and Alpine Research, 2013, 45(3):363-371.

    [36] Bai Y F, Wu J G, Clark C M,etal.Trade offs and thresholds in the effects of nitrogen addition on biodiversity and ecosystem functioning:Evidence from Inner Mongolia Grasslands.Global Change Biology, 2010, 16(1):358-372.

    [37] Zhang L J, Yue M, Gu F X,etal.Coupling relationship between plant communities’ species diversity and soil factors in ecotone between desert and oasis in Fukang, Xinjiang.Chinese Journal of Applied Ecology, 2002, 13(6):658-662.

    [38] Jouany C, Cruz P, Daufresne T,etal.Biological phosphorus cycling in grasslands:Interactions with nitrogen.In:Bunemann E K, Oberson A, Frossard E.Phosphorus in Action:Biological Processes in Soil Phosphorus Cycling[M].Berlin:Springer, 2011, 275-294.

    [39] Perroni Y, Garcia-Oliva F, Souza V.Plant species identity and soil P forms in an oligotrophic grassland-desert scrub system.Journal of Arid Environments, 2014, 108:29-37.

    [40] Grime J P.Biodiversity and ecosystem function:The debate deepens.Science, 1997, 277:1260-1261.

    [41] Tilman D.Secondary succession and the pattern of plant dominance along experiment nitrogen gradients.Ecological Monographs, 1987, 57:189-214.

    [42] Vitousek P M, Aber J D, Howarth R W,etal.Human alteration of the global nitrogen cycle:Sources and consequences.Ecological Applications, 1997, 7(3):737-750.

    [43] Tilman D, Knops J, Wedin D,etal.The influence of functional diversity and composition on ecosystem processes.Science, 1997, 277:1300-1302.

    [44] Lee M, Manning P, Rist J,etal.A global comparison of grassland biomass responses to CO2and nitrogen enrichment.Philosophical Transactions of the Royal Society B-Biological Sciences, 2010, 365:2047-2056.

    [45] Zhen X X, Jin T T, Mu L F,etal.The relationship between plant species richness in Hulunbeier grassland and biomass and environment factors.Chinese Journal of Grassland, 2008, 30(6):74-81.

    [46] Field R, Hawkins B A, Cornell H V,etal.Spatial species-richness gradients across scales:A meta-analysis.Journal of Biogeography, 2009, 36(1):132-147.

    [47] Whittaker R J, Nogues-Bravo D, Araujo M B.Geographical gradients of species richness:A test of the water-energy conjecture of hawkinsetal.(2003) using european data for five taxa.Global Ecology and Biogeography, 2007, 16(1):76-89.

    [48] Tasker E M, Bradstock R A.Influence of cattle grazing practices on forest understory structure in north-eastern New South Wales.Austral Ecology, 2006, 31(4):490-502.

    [49] Milchunas D G, Lauenroth W K.Quantitative effects of grazing on vegetation and soils over a global range of environments.Ecological Monographs, 1993, 63(4):327-366.

    [50] Pykala J.Effects of restoration with cattle grazing on plant species composition and richness of semi-natural grasslands.Biodiversity and Conservation, 2003, 12(11):2211-2226.

    [51] Proulx M, Mazumder A.Reversal of grazing impact on plant species richness in nutrient-poor vs.nutrient-rich ecosystems.Ecology, 1998, 79(8):2581-2592.

    [52] Thompson K, Askew A P, Grime J P,etal.Biodiversity, ecosystem function and plant traits in mature and immature plant communities.Functional Ecology, 2005, 19(2):355-358.

    [53] Kassen R, Buckling A, Bell G,etal.Diversity peaks at intermediate productivity in a laboratory microcosm.Nature, 2000, 406:508-512.

    [54] Grace J B, Michael Anderson T, Smith M D,etal.Does species diversity limit productivity in natural grassland communities.Ecology Letters, 2007, 10(8):680-689.

    [55] Moore D R J, Keddy P A.The relationship between species richness and standing crop in wetlands the importance of scale.Vegetatio, 1988, 79(1-2):99-106.

    [56] Waide R B, Willig M R, Steiner C F,etal.The relationship between productivity and species richness.Annual Review of Ecology and Systematics, 1999, 30:257-300.

    [57] Ni J, Wang G H, Bai Y F,etal.Scale-dependent relationships between plant diversity and above-ground biomass in temperate grasslands, south-eastern Mongolia.Journal of arid environments, 2007, 68(1):132-142.

    [58] He J S, Bazzaz F A, Schmid B.Interactive effects of diversity, nutrients and elevated CO2on experimental plant communities.Oikos, 2002, 97, 337-348.

    參考文獻(xiàn):

    [5] 左小安, 趙學(xué)勇, 趙哈林, 等.科爾沁沙質(zhì)草地群落物種多樣性、生產(chǎn)力與土壤特性的關(guān)系.環(huán)境科學(xué), 2007, 28(5):945-951.

    [8] 白永飛, 李凌浩, 王其兵, 等.錫林河流域草原群落植物多樣性和初級(jí)生產(chǎn)力沿水熱梯度變化的樣帶研究.植物生態(tài)學(xué)報(bào), 2000, 24(6):667-673.

    [9] 張志南, 武高林, 王冬, 等.黃土高原半干旱區(qū)天然草地群落結(jié)構(gòu)與土壤水分關(guān)系.草業(yè)學(xué)報(bào), 2014, 23(6):313-319.

    [13] 安樹青, 王崢?lè)? 朱學(xué)雷, 等.土壤因子對(duì)次生森林群落物種多樣性的影響.武漢植物學(xué)研究, 1997, 15(2):143-150.

    [17] 劉旻霞, 王剛.高寒草甸植物群落多樣性及土壤因子對(duì)坡向的響應(yīng).生態(tài)學(xué)雜志, 2013, 32(2):259-265.

    [18] 王長(zhǎng)庭, 龍瑞軍, 曹廣民, 等.高寒草甸不同類型草地土壤養(yǎng)分與物種多樣性——生產(chǎn)力關(guān)系.土壤通報(bào), 2008, 39(1):1-8.

    [19] 中國(guó)科學(xué)院內(nèi)蒙寧夏綜合考察隊(duì).內(nèi)蒙古植被[M].北京:科學(xué)出版社, 1985.

    [20] 張慶, 牛建明, Buyantuyev, 等.內(nèi)蒙古短花針茅群落數(shù)量分類及環(huán)境解釋.草業(yè)學(xué)報(bào), 2012, 21(1):83-92.

    [21] 楊雪棟, 孫衛(wèi)國(guó), 寶音陶格濤.內(nèi)蒙古荒漠草原小針茅群落土壤養(yǎng)分的化學(xué)計(jì)量學(xué)特征.中國(guó)草地學(xué)報(bào), 2012, 34(5):30-34.

    [22] 馬文靜, 張慶, 牛建明, 等.物種多樣性和功能群多樣性與生態(tài)系統(tǒng)生產(chǎn)力的關(guān)系——以內(nèi)蒙古短花針茅草原為例.植物生態(tài)學(xué)報(bào), 2013, 37(7):620-630.

    [23] 衛(wèi)智軍, 韓國(guó)棟, 趙鋼, 等.中國(guó)荒漠草原生態(tài)系統(tǒng)研究[M].北京:科學(xué)出版社, 2013.

    [24] 張甘霖, 龔子同.土壤調(diào)查實(shí)驗(yàn)室分析方法[M].北京:科學(xué)出版社, 2012.

    [27] 韓彬, 樊江文, 鐘華平.內(nèi)蒙古草地樣帶植物群落生物量的梯度研究.植物生態(tài)學(xué)報(bào), 2006, 30(4):553-562.

    [33] 張杰琦, 李奇, 任正煒, 等.氮素添加對(duì)青藏高原高寒草甸植物群落物種豐富度及其與地上生產(chǎn)力關(guān)系的影響.植物生態(tài)學(xué)報(bào), 2010, 34(10):1125-1131.

    [37] 張林靜, 岳明, 顧峰雪, 等.新疆阜康綠洲荒漠過(guò)渡帶植物群落物種多樣性與土壤環(huán)境因子的耦合關(guān)系.應(yīng)用生態(tài)學(xué)報(bào), 2002, 13(6):658-662.

    [45] 鄭曉翾, 靳甜甜, 木麗芬, 等.呼倫貝爾草原物種多樣性與生物量、環(huán)境因子的關(guān)系.中國(guó)草地學(xué)報(bào), 2008, 30(6):74-81.

    Relationship between species diversity, productivity, climatic factors and soil nutrients in the desert steppe

    SUN Xiao-Li1, KANG Sarula1, ZHANG Qing1,2, CHANG Chang-Ming1, MA Wen-Jing1,NIU Jian-Ming1,2*

    1.CollegeofLifeSciences,InnerMongoliaUniversity,Hohhot010021,China; 2.Sino-USCenterforConservation,EnergyandSustainabilityScienceinInnerMongolia,Hohhot010021,China

    The effect of species diversity on ecosystem function is one of the core issues in ecology.Our objective was to explore the relationship between species diversity, productivity, climatic factors, and soil nutrients in the desert steppe, to enhance understanding of the diversity-productivity relationship.Species richness and Shannon-Wiener index were used to measure species diversity.Climatic, vegetation and soil data were collected for 80 sampling sites across the desert steppe of Inner Mongolia.Regression analysis and stepwise regression analysis were used for determining the relationship between environmental factors and species diversity and productivity.The effects of environmental factors on species diversity and productivity were evaluated, based on the coefficient of determination from regression analysis.It was found that species diversity and biomass were both positively and linearly correlated with mean annual precipitation and aridity index.However, Shannon-Wiener index and biomass were both negatively and linearly correlated with mean annual temperature.Species richness was negatively and linearly correlated with soil total phosphorus content (R2=0.072,P<0.05), but positively and linearly correlated with soil available nitrogen content (R2=0.183,P<0.01); Shannon-Wiener index was also positively and linearly correlated with soil available nitrogen content (R2=0.061,P<0.05), and biomass was positively and linearly correlated with total nitrogen, soil available nitrogen, and soil organic carbon content.The coefficients of determination were 0.150, 0.177, 0.089 (P<0.01), respectively.Biomass increased significantly with the increase of species richness and Shannon-Wiener index.These results indicate that aridity index is the most influential factor affecting species diversity and biomass of desert steppe.

    desert steppe; productivity; species diversity; soil nutrient; precipitation and temperature

    10.11686/cyxb2015014

    http://cyxb.lzu.edu.cn

    2015-01-08;改回日期:2015-04-08

    國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973)項(xiàng)目(2012CB722201),國(guó)家科技支撐計(jì)劃課題(2011BAC07B01)和國(guó)家自然科學(xué)基金(31060320)資助。

    孫小麗(1988-),女,內(nèi)蒙古赤峰人,在讀碩士。 E-mail:sxl890325@163.com*通信作者Corresponding author.E-mail:jmniu2005@163.com

    孫小麗,康薩如拉,張慶,常昌明,馬文靜,牛建明.荒漠草原物種多樣性、生產(chǎn)力與氣候因子和土壤養(yǎng)分之間關(guān)系的研究.草業(yè)學(xué)報(bào), 2015, 24(12):10-19.

    SUN Xiao-Li, KANG Sarula, ZHANG Qing, CHANG Chang-Ming, MA Wen-Jing, NIU Jian-Ming.Relationship between species diversity, productivity, climatic factors and soil nutrients in the desert steppe.Acta Prataculturae Sinica, 2015, 24(12):10-19.

    猜你喜歡
    荒漠生產(chǎn)力生物量
    向荒漠宣戰(zhàn)
    中國(guó)石化(2021年8期)2021-11-05 07:00:14
    輪牧能有效促進(jìn)高寒草地生物量和穩(wěn)定性
    荒漠生態(tài)系統(tǒng)的演變
    把“風(fēng)景”轉(zhuǎn)化成“生產(chǎn)力
    紅土地(2016年3期)2017-01-15 13:45:17
    荒漠之路上的三首情歌
    文史春秋(2016年3期)2016-12-01 05:42:20
    生物量高的富鋅酵母的開發(fā)應(yīng)用
    人口紅利與提高生產(chǎn)力
    海峽姐妹(2015年9期)2015-02-27 15:13:07
    基于SPOT-5遙感影像估算玉米成熟期地上生物量及其碳氮累積量
    展望2014:推動(dòng)科技創(chuàng)新轉(zhuǎn)化為第一生產(chǎn)力
    杭州科技(2014年3期)2014-02-27 15:26:46
    色94色欧美一区二区| 国产精品成人在线| 国产精品秋霞免费鲁丝片| 国产日韩一区二区三区精品不卡 | 精品一区在线观看国产| 国产淫片久久久久久久久| 丝袜在线中文字幕| kizo精华| 最新的欧美精品一区二区| 免费不卡的大黄色大毛片视频在线观看| 成人18禁高潮啪啪吃奶动态图 | 熟女人妻精品中文字幕| 一级爰片在线观看| 精品酒店卫生间| 亚洲真实伦在线观看| 2018国产大陆天天弄谢| 插阴视频在线观看视频| 伊人久久国产一区二区| 三级经典国产精品| 精华霜和精华液先用哪个| 99九九线精品视频在线观看视频| 下体分泌物呈黄色| 黑人巨大精品欧美一区二区蜜桃 | 亚洲三级黄色毛片| 边亲边吃奶的免费视频| 日韩中文字幕视频在线看片| av不卡在线播放| 亚洲av在线观看美女高潮| 日本黄色日本黄色录像| 69精品国产乱码久久久| 一二三四中文在线观看免费高清| 亚洲精品国产av蜜桃| 欧美97在线视频| .国产精品久久| 成人黄色视频免费在线看| 91久久精品国产一区二区三区| 美女cb高潮喷水在线观看| 成年av动漫网址| 亚洲,欧美,日韩| 久久免费观看电影| 亚洲av成人精品一二三区| 一级片'在线观看视频| av一本久久久久| 九色成人免费人妻av| 国国产精品蜜臀av免费| 日韩成人av中文字幕在线观看| 国产成人免费无遮挡视频| 热re99久久国产66热| 91精品伊人久久大香线蕉| 男人爽女人下面视频在线观看| 国产精品久久久久成人av| 能在线免费看毛片的网站| 春色校园在线视频观看| 亚洲va在线va天堂va国产| 亚洲四区av| 亚洲欧美日韩另类电影网站| 日本av免费视频播放| av播播在线观看一区| 国产在线男女| 十八禁高潮呻吟视频 | 女人久久www免费人成看片| 边亲边吃奶的免费视频| 日韩中字成人| 有码 亚洲区| 亚洲av电影在线观看一区二区三区| 偷拍熟女少妇极品色| 91精品伊人久久大香线蕉| 一区二区三区乱码不卡18| 麻豆成人午夜福利视频| 亚洲精品第二区| 18+在线观看网站| 日韩av不卡免费在线播放| 免费久久久久久久精品成人欧美视频 | 欧美日韩av久久| 久久狼人影院| 久久久亚洲精品成人影院| 亚洲欧洲日产国产| 老司机影院成人| 国产中年淑女户外野战色| 一级毛片aaaaaa免费看小| 两个人的视频大全免费| 黑人猛操日本美女一级片| 少妇人妻精品综合一区二区| 日本黄色日本黄色录像| 丰满人妻一区二区三区视频av| 中文字幕精品免费在线观看视频 | 三上悠亚av全集在线观看 | 欧美 亚洲 国产 日韩一| 老女人水多毛片| 亚洲电影在线观看av| 亚洲一级一片aⅴ在线观看| 春色校园在线视频观看| 久久精品熟女亚洲av麻豆精品| 男的添女的下面高潮视频| 色哟哟·www| 黄色视频在线播放观看不卡| 伦精品一区二区三区| 国产 一区精品| 综合色丁香网| 国产乱人偷精品视频| 性高湖久久久久久久久免费观看| 人人妻人人澡人人看| 亚洲婷婷狠狠爱综合网| 日韩一区二区视频免费看| 好男人视频免费观看在线| 在现免费观看毛片| 少妇被粗大的猛进出69影院 | 观看免费一级毛片| 日本黄大片高清| 777米奇影视久久| 国产亚洲精品久久久com| 日本欧美国产在线视频| 蜜臀久久99精品久久宅男| av福利片在线| 少妇的逼好多水| 校园人妻丝袜中文字幕| 人妻少妇偷人精品九色| 97在线人人人人妻| 人妻制服诱惑在线中文字幕| 欧美日韩视频精品一区| 丰满迷人的少妇在线观看| 色婷婷久久久亚洲欧美| 国产片特级美女逼逼视频| 久久久亚洲精品成人影院| 最黄视频免费看| 亚洲伊人久久精品综合| 亚洲内射少妇av| 嫩草影院新地址| 婷婷色综合www| 麻豆精品久久久久久蜜桃| 丝瓜视频免费看黄片| 99久久精品一区二区三区| 成人免费观看视频高清| 22中文网久久字幕| 如日韩欧美国产精品一区二区三区 | 久久精品国产亚洲av天美| 亚洲第一av免费看| 男女免费视频国产| 国产黄色免费在线视频| 高清av免费在线| 免费观看的影片在线观看| 91成人精品电影| 女性生殖器流出的白浆| 观看av在线不卡| 我要看日韩黄色一级片| 老司机亚洲免费影院| 久久久a久久爽久久v久久| 国产精品99久久久久久久久| 久久热精品热| 欧美精品国产亚洲| 高清欧美精品videossex| 十八禁高潮呻吟视频 | 一级毛片电影观看| 热99国产精品久久久久久7| 国产亚洲av片在线观看秒播厂| 久久鲁丝午夜福利片| 国产精品成人在线| 精品久久国产蜜桃| 亚洲四区av| 久久久久久久国产电影| 蜜桃在线观看..| 亚洲人成网站在线观看播放| 香蕉精品网在线| 亚洲国产精品国产精品| 人体艺术视频欧美日本| 少妇人妻久久综合中文| 大片免费播放器 马上看| 日韩 亚洲 欧美在线| 国产av一区二区精品久久| 边亲边吃奶的免费视频| 国产成人精品婷婷| 免费观看的影片在线观看| 99热这里只有精品一区| 色94色欧美一区二区| 日产精品乱码卡一卡2卡三| 国产精品三级大全| 午夜福利影视在线免费观看| 亚洲精品成人av观看孕妇| 国产伦精品一区二区三区四那| 天堂中文最新版在线下载| 日本wwww免费看| 毛片一级片免费看久久久久| 自线自在国产av| 日韩熟女老妇一区二区性免费视频| 日日啪夜夜撸| 精华霜和精华液先用哪个| 一级二级三级毛片免费看| 一本久久精品| 麻豆乱淫一区二区| 精品一区二区三卡| 欧美激情极品国产一区二区三区 | 国产免费视频播放在线视频| 久久精品夜色国产| 又爽又黄a免费视频| 免费大片黄手机在线观看| 交换朋友夫妻互换小说| 又大又黄又爽视频免费| 久久6这里有精品| 在线观看www视频免费| 亚洲av在线观看美女高潮| 一区二区三区乱码不卡18| 成年人午夜在线观看视频| 免费黄频网站在线观看国产| 一区二区三区免费毛片| 观看免费一级毛片| 国产熟女欧美一区二区| 少妇被粗大猛烈的视频| 成人18禁高潮啪啪吃奶动态图 | 日本wwww免费看| 久久久久久久精品精品| 国产乱来视频区| kizo精华| 亚洲精品,欧美精品| 少妇的逼好多水| 91aial.com中文字幕在线观看| 9色porny在线观看| 亚洲国产成人一精品久久久| 国产色爽女视频免费观看| 国产永久视频网站| 欧美另类一区| 国内精品宾馆在线| tube8黄色片| 高清黄色对白视频在线免费看 | 久久久久人妻精品一区果冻| 国产精品成人在线| 国产精品一区www在线观看| 久久人人爽人人片av| 男女无遮挡免费网站观看| 精品国产露脸久久av麻豆| 极品人妻少妇av视频| 免费观看无遮挡的男女| 欧美 日韩 精品 国产| 99国产精品免费福利视频| 欧美亚洲 丝袜 人妻 在线| 男女国产视频网站| 国产精品一区www在线观看| 少妇精品久久久久久久| 亚洲欧美一区二区三区黑人 | 好男人视频免费观看在线| 午夜精品国产一区二区电影| 黄色配什么色好看| 国产精品秋霞免费鲁丝片| 久久6这里有精品| 纯流量卡能插随身wifi吗| 日本黄大片高清| 啦啦啦中文免费视频观看日本| 精品久久久噜噜| 五月伊人婷婷丁香| 国产精品一区www在线观看| 色94色欧美一区二区| 欧美国产精品一级二级三级 | 亚洲欧洲日产国产| 精品亚洲成国产av| 边亲边吃奶的免费视频| 啦啦啦中文免费视频观看日本| 精品一品国产午夜福利视频| 亚洲av.av天堂| 天堂俺去俺来也www色官网| 精品熟女少妇av免费看| 王馨瑶露胸无遮挡在线观看| 黄色一级大片看看| 在线观看国产h片| 亚洲第一av免费看| 日韩大片免费观看网站| 黑人巨大精品欧美一区二区蜜桃 | 在线看a的网站| 卡戴珊不雅视频在线播放| av视频免费观看在线观看| 99久久精品热视频| 国产一级毛片在线| 国产精品国产三级专区第一集| 免费观看的影片在线观看| 九草在线视频观看| 亚洲情色 制服丝袜| 精品视频人人做人人爽| 97超视频在线观看视频| 黑人巨大精品欧美一区二区蜜桃 | 交换朋友夫妻互换小说| 日本欧美国产在线视频| 国产有黄有色有爽视频| 91精品伊人久久大香线蕉| 日日摸夜夜添夜夜添av毛片| 国国产精品蜜臀av免费| 中文字幕久久专区| 亚州av有码| 99热这里只有精品一区| 五月伊人婷婷丁香| 亚洲精品456在线播放app| 精品人妻偷拍中文字幕| 日韩亚洲欧美综合| 你懂的网址亚洲精品在线观看| 精品少妇久久久久久888优播| 午夜视频国产福利| 一本色道久久久久久精品综合| 丝袜在线中文字幕| 自拍偷自拍亚洲精品老妇| 亚洲欧美中文字幕日韩二区| 国产精品一区二区性色av| 国产免费一区二区三区四区乱码| 人妻制服诱惑在线中文字幕| 日本爱情动作片www.在线观看| 国产精品国产av在线观看| 欧美精品国产亚洲| 日本欧美视频一区| 五月玫瑰六月丁香| 午夜影院在线不卡| 国产精品国产三级国产av玫瑰| 亚洲人成网站在线播| 亚洲真实伦在线观看| 国产 精品1| 夜夜看夜夜爽夜夜摸| 黄色毛片三级朝国网站 | 最后的刺客免费高清国语| 不卡视频在线观看欧美| 久久久久久久精品精品| 日日摸夜夜添夜夜添av毛片| 中国国产av一级| 国产欧美日韩精品一区二区| 精品久久国产蜜桃| 成人午夜精彩视频在线观看| 有码 亚洲区| 日本91视频免费播放| 又大又黄又爽视频免费| 自线自在国产av| 五月伊人婷婷丁香| 纯流量卡能插随身wifi吗| 亚洲精品乱码久久久久久按摩| 久久久久久久精品精品| 最近的中文字幕免费完整| 亚洲国产最新在线播放| 一级爰片在线观看| 亚洲欧美一区二区三区黑人 | 观看av在线不卡| 国产亚洲最大av| 一级黄片播放器| 亚洲久久久国产精品| 在线观看免费日韩欧美大片 | 国产亚洲欧美精品永久| 黄色毛片三级朝国网站 | 久久国产乱子免费精品| 男男h啪啪无遮挡| 久久婷婷青草| freevideosex欧美| 精品久久国产蜜桃| 亚州av有码| 久久精品国产亚洲网站| 人妻一区二区av| 七月丁香在线播放| 狂野欧美白嫩少妇大欣赏| 看免费成人av毛片| 日韩三级伦理在线观看| 超碰97精品在线观看| 国产免费一级a男人的天堂| 国产成人免费无遮挡视频| 亚洲精品一二三| 最近最新中文字幕免费大全7| 久久人人爽人人片av| 欧美日韩一区二区视频在线观看视频在线| 欧美性感艳星| 中文字幕亚洲精品专区| 国产黄色免费在线视频| 国产在视频线精品| 91精品一卡2卡3卡4卡| 丰满人妻一区二区三区视频av| 国产精品国产三级国产专区5o| 久久久国产精品麻豆| 中文精品一卡2卡3卡4更新| 亚洲精品中文字幕在线视频 | 亚洲精品国产av成人精品| 日韩制服骚丝袜av| 一本大道久久a久久精品| 国产成人一区二区在线| 国产成人免费观看mmmm| 免费在线观看成人毛片| 最近中文字幕2019免费版| 国产国拍精品亚洲av在线观看| 欧美3d第一页| 熟女电影av网| 久久ye,这里只有精品| 欧美成人午夜免费资源| 欧美精品人与动牲交sv欧美| 嫩草影院新地址| 成人综合一区亚洲| 日本色播在线视频| 亚洲欧美清纯卡通| 午夜免费男女啪啪视频观看| 十分钟在线观看高清视频www | 十八禁网站网址无遮挡 | 91成人精品电影| 七月丁香在线播放| 国产色婷婷99| 免费av中文字幕在线| 午夜视频国产福利| 人人澡人人妻人| 国产探花极品一区二区| 丰满迷人的少妇在线观看| 亚洲性久久影院| 黄色怎么调成土黄色| 精品人妻熟女毛片av久久网站| 视频中文字幕在线观看| 欧美激情极品国产一区二区三区 | 成人18禁高潮啪啪吃奶动态图 | 久久久久久久久久人人人人人人| 我要看黄色一级片免费的| 免费观看在线日韩| 丝袜脚勾引网站| 久久久亚洲精品成人影院| 少妇熟女欧美另类| 日本av免费视频播放| 国产精品蜜桃在线观看| 精品国产一区二区久久| 国产亚洲av片在线观看秒播厂| 日韩三级伦理在线观看| 国产av一区二区精品久久| 久久影院123| 久久av网站| 色视频在线一区二区三区| 九草在线视频观看| 人妻人人澡人人爽人人| 亚洲av欧美aⅴ国产| 午夜激情福利司机影院| 日本午夜av视频| 久久精品国产亚洲av天美| 十分钟在线观看高清视频www | 精品久久久久久电影网| av在线老鸭窝| 丰满乱子伦码专区| 欧美激情国产日韩精品一区| 亚洲欧美日韩卡通动漫| 欧美日韩综合久久久久久| 婷婷色综合www| 午夜福利网站1000一区二区三区| 国产精品国产三级国产专区5o| 亚洲熟女精品中文字幕| 高清在线视频一区二区三区| 中文精品一卡2卡3卡4更新| 久久久精品免费免费高清| 丰满迷人的少妇在线观看| 国产成人freesex在线| 黄色毛片三级朝国网站 | 亚洲精品第二区| 日本av免费视频播放| 久久久久久人妻| 18禁裸乳无遮挡动漫免费视频| 亚洲丝袜综合中文字幕| 成人漫画全彩无遮挡| www.av在线官网国产| 嫩草影院入口| 国产成人精品无人区| 最近的中文字幕免费完整| 深夜a级毛片| 九九爱精品视频在线观看| 一本色道久久久久久精品综合| 精品一品国产午夜福利视频| 综合色丁香网| 大陆偷拍与自拍| 国产片特级美女逼逼视频| 汤姆久久久久久久影院中文字幕| 亚洲人成网站在线播| 18禁裸乳无遮挡动漫免费视频| 午夜免费观看性视频| 欧美性感艳星| 夫妻性生交免费视频一级片| 色婷婷久久久亚洲欧美| 国产极品天堂在线| 中文字幕亚洲精品专区| 97超视频在线观看视频| 亚洲精华国产精华液的使用体验| 黄色毛片三级朝国网站 | 亚洲国产精品成人久久小说| 欧美亚洲 丝袜 人妻 在线| 丝袜在线中文字幕| 国产精品成人在线| 十八禁高潮呻吟视频 | 欧美激情国产日韩精品一区| 成人无遮挡网站| 成人毛片60女人毛片免费| 久久人人爽人人片av| 26uuu在线亚洲综合色| 乱码一卡2卡4卡精品| 国产亚洲午夜精品一区二区久久| 十八禁网站网址无遮挡 | 亚洲精品自拍成人| 极品教师在线视频| 国产精品秋霞免费鲁丝片| 夜夜爽夜夜爽视频| 久久久久人妻精品一区果冻| 黄色日韩在线| 水蜜桃什么品种好| 亚洲熟女精品中文字幕| 少妇精品久久久久久久| 国产日韩一区二区三区精品不卡 | 高清黄色对白视频在线免费看 | 国产成人91sexporn| 久久精品熟女亚洲av麻豆精品| 天堂8中文在线网| 黄色日韩在线| 最近的中文字幕免费完整| 3wmmmm亚洲av在线观看| 欧美一级a爱片免费观看看| 黑人高潮一二区| 丝袜喷水一区| 国产精品成人在线| 最近最新中文字幕免费大全7| 女的被弄到高潮叫床怎么办| 2022亚洲国产成人精品| 黄色视频在线播放观看不卡| 国产高清不卡午夜福利| 美女xxoo啪啪120秒动态图| 成人综合一区亚洲| 国产成人免费观看mmmm| 亚洲图色成人| 亚洲av成人精品一二三区| 91精品国产九色| 欧美日韩在线观看h| 国产69精品久久久久777片| 日日啪夜夜撸| 黄色视频在线播放观看不卡| 中文天堂在线官网| 麻豆成人午夜福利视频| 亚洲精品aⅴ在线观看| 高清欧美精品videossex| 九九在线视频观看精品| 午夜老司机福利剧场| 久久久久久久久久人人人人人人| 久久午夜福利片| 亚洲经典国产精华液单| 国产av国产精品国产| 一个人看视频在线观看www免费| 免费大片18禁| 人人妻人人爽人人添夜夜欢视频 | 一区二区三区免费毛片| 熟女电影av网| 我要看黄色一级片免费的| 人妻夜夜爽99麻豆av| 在线 av 中文字幕| 高清黄色对白视频在线免费看 | 欧美精品高潮呻吟av久久| 大片电影免费在线观看免费| 亚洲精品国产av蜜桃| 日韩成人伦理影院| 国产一区二区在线观看日韩| 涩涩av久久男人的天堂| 久久久久视频综合| 国产日韩欧美在线精品| 777米奇影视久久| 国产在线男女| 久久综合国产亚洲精品| 免费大片黄手机在线观看| 大片免费播放器 马上看| 亚洲av二区三区四区| 欧美xxⅹ黑人| 人妻人人澡人人爽人人| 熟女电影av网| 午夜福利影视在线免费观看| 日本-黄色视频高清免费观看| 老司机影院成人| 亚洲人与动物交配视频| 老司机影院毛片| 日日爽夜夜爽网站| 欧美激情极品国产一区二区三区 | 内射极品少妇av片p| 嫩草影院新地址| 看十八女毛片水多多多| 精品卡一卡二卡四卡免费| 中文字幕人妻丝袜制服| 曰老女人黄片| 国产黄色免费在线视频| 国产精品伦人一区二区| 少妇人妻 视频| 一个人免费看片子| 免费久久久久久久精品成人欧美视频 | 成人午夜精彩视频在线观看| 丰满迷人的少妇在线观看| 国产精品99久久久久久久久| 黑人高潮一二区| 国产在线一区二区三区精| 亚洲精品乱久久久久久| 人妻一区二区av| 在线 av 中文字幕| 少妇人妻 视频| 亚洲欧美成人精品一区二区| 夜夜骑夜夜射夜夜干| 少妇人妻一区二区三区视频| 亚洲丝袜综合中文字幕| 国产免费一区二区三区四区乱码| 自拍偷自拍亚洲精品老妇| 国产免费又黄又爽又色| 91久久精品国产一区二区成人| 国产永久视频网站| videos熟女内射| 国产综合精华液| 91久久精品电影网| 亚洲精品乱码久久久v下载方式| 狂野欧美激情性xxxx在线观看| 你懂的网址亚洲精品在线观看| 如何舔出高潮| 激情五月婷婷亚洲| 免费观看在线日韩| 精品99又大又爽又粗少妇毛片| 成人特级av手机在线观看| 一级毛片黄色毛片免费观看视频| 如何舔出高潮| 国产视频首页在线观看| 你懂的网址亚洲精品在线观看| 国产精品蜜桃在线观看| 91久久精品国产一区二区三区| 中文字幕人妻丝袜制服| 狠狠精品人妻久久久久久综合| 激情五月婷婷亚洲| 午夜福利在线观看免费完整高清在| 国产欧美日韩一区二区三区在线 | 51国产日韩欧美| 久久久久久久久大av| 久久久久久久久久成人|