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

    黃土區(qū)封育和放牧草地物種多樣性和地上生物量對(duì)氮素添加的響應(yīng)

    2017-09-03 08:59:10趙潔李偉井光花魏琳程積民
    草業(yè)學(xué)報(bào) 2017年8期
    關(guān)鍵詞:物種

    趙潔,李偉,井光花,魏琳,程積民,*

    (1.西北農(nóng)林科技大學(xué)動(dòng)物科技學(xué)院, 陜西 楊凌 712100; 2.中國(guó)科學(xué)院水利部水土保持研究所, 陜西 楊凌 712100)

    黃土區(qū)封育和放牧草地物種多樣性和地上生物量對(duì)氮素添加的響應(yīng)

    趙潔1,李偉2,井光花2,魏琳2,程積民1,2*

    (1.西北農(nóng)林科技大學(xué)動(dòng)物科技學(xué)院, 陜西 楊凌 712100; 2.中國(guó)科學(xué)院水利部水土保持研究所, 陜西 楊凌 712100)

    氮素添加對(duì)草地生態(tài)系統(tǒng)地上生物量具有普遍促進(jìn)作用,但對(duì)解釋引起物種多樣性喪失的機(jī)制仍存在較大爭(zhēng)議。本研究以黃土高原半干旱區(qū)長(zhǎng)期封育和連續(xù)放牧草地為對(duì)象,通過(guò)設(shè)置不同氮素添加水平來(lái)探索兩種典型的草地管理方式下群落地上生物量、物種多樣性、群落組分結(jié)構(gòu)對(duì)氮素添加的響應(yīng)及其異同以及物種多樣性與地上生物量的關(guān)系,并研究這些響應(yīng)的潛在生態(tài)學(xué)機(jī)制。結(jié)果顯示,1)氮素添加顯著提高長(zhǎng)期封育草地地上生物量的同時(shí)降低了物種多樣性;放牧草地群落地上生物量、物種多樣性對(duì)氮素添加的響應(yīng)與封育草地相同,但其響應(yīng)程度均小于封育施肥草地;2)草地群落物種多樣性隨地表凋落物生物量增加而降低,隨植被透光率增加而升高; 3)封育施肥草地Shannon-Wiener多樣性指數(shù)與地上生物量存在顯著線性負(fù)相關(guān)關(guān)系,而放牧施肥草地的物種多樣性與地上生物量不存在顯著的相關(guān)關(guān)系。以上結(jié)果表明放牧能夠減緩施肥對(duì)物種多樣性的負(fù)效應(yīng),在輕微降低物種多樣性的前提下顯著提高群落地上生物量,是一種合理的草地管理方式。

    放牧;封育;氮素添加;物種多樣性;地上生物量;群落組分

    多數(shù)研究表明氮素添加在提高生物量的同時(shí)會(huì)造成草地生態(tài)系統(tǒng)物種多樣性下降[1-2]。氮素添加對(duì)草地地上生物量的影響在年限和氮素水平上可能存在閾值,草地地上生物量在到達(dá)此閾值時(shí)達(dá)到峰值,之后氮素對(duì)地上生物量促進(jìn)效應(yīng)將減小甚至轉(zhuǎn)為負(fù)向影響[1,3-4]。Stevens等[5]根據(jù)逐步回歸預(yù)測(cè)氮沉降速率為0.25 g N/(m2·yr)時(shí),平均每4 m2將有一個(gè)物種消失。但對(duì)引起物種多樣性下降的原因尚存在爭(zhēng)議,目前人們提出的解釋機(jī)制主要有以下幾種:光競(jìng)爭(zhēng)、凋落物積累、集合水平自疏假說(shuō)、銨中毒與土壤酸化[6-10]。其中光競(jìng)爭(zhēng)的理論認(rèn)為,氮素的增加解除了生態(tài)系統(tǒng)對(duì)氮素的限制,嗜氮植物快速生長(zhǎng),在養(yǎng)分吸收、高度、蓋度上占得優(yōu)勢(shì),蔭蔽一些非嗜氮、生長(zhǎng)較慢或矮小的物種,降低物種多樣性。即種間競(jìng)爭(zhēng)由地下對(duì)養(yǎng)分的競(jìng)爭(zhēng)轉(zhuǎn)為地上對(duì)光的競(jìng)爭(zhēng),植物對(duì)養(yǎng)分和光競(jìng)爭(zhēng)力的差異決定了最終群落物種組成格局[11]。另一方面,氮素增加后伴隨著植物生產(chǎn)力提高,凋落物積累量也增加,影響了一些種子萌發(fā)和幼苗生長(zhǎng),也是物種多樣性降低的原因之一[12]。植物群落與環(huán)境之間的相互作用極其復(fù)雜多面,上述研究及假說(shuō)雖然可以從一定程度上解釋物種多樣性對(duì)施肥的響應(yīng),但是到目前還沒(méi)有較為統(tǒng)一的認(rèn)識(shí),氮沉降導(dǎo)致物種多樣性喪失的機(jī)制還有很大的研究空間。

    盡管?chē)?guó)內(nèi)外學(xué)者對(duì)氮沉降已經(jīng)進(jìn)行了許多研究[13-17],但對(duì)氮素和水分限制的黃土高原半干旱草地生態(tài)系統(tǒng)的研究還相對(duì)較少,且其對(duì)氮素添加的響應(yīng)可能更加敏感[18]。同時(shí)近年來(lái)也有少數(shù)研究表明放牧可能會(huì)減緩施肥對(duì)物種多樣性的改變[3,19]?;谏鲜鲈蚣扒叭搜芯炕A(chǔ)選取黃土高原典型草原長(zhǎng)期封育和連續(xù)放牧草地為研究對(duì)象,通過(guò)氮素添加實(shí)驗(yàn)研究?jī)煞N管理方式下草地地上生物量和物種多樣性對(duì)氮素添加的響應(yīng),為黃土高原已恢復(fù)天然草地的后續(xù)合理利用以及我國(guó)生態(tài)草牧業(yè)的可持續(xù)發(fā)展提供理論依據(jù)。本研究的主要目標(biāo)是:1)比較封育和放牧草地地上生物量和物種多樣性隨著不同氮素添加水平變化的異同;2)研究放牧和施肥對(duì)物種多樣性的影響及原因;3)探索黃土高原典型草原上物種多樣性與地上生物量的關(guān)系。

    1 材料與方法

    1.1 研究地概況

    研究地位于寧夏固原市東北部的云霧山國(guó)家草地自然保護(hù)區(qū)實(shí)驗(yàn)區(qū)(106°21′-106°27′ E, 36°10′-36°17′ N),海拔1800~2100 m。該區(qū)域氣候?qū)儆谥袦貛О敫珊祬^(qū),年平均氣溫5 ℃,最冷月1月(平均氣溫-14 ℃),最熱月7月(22~25 ℃),年降水量400~450 mm,降雨季節(jié)分配不均,全年65%~75%的降水集中在7-9月,蒸發(fā)量1330~1640 mm。年日照時(shí)數(shù)為2500 h,無(wú)霜期112~140 d。

    保護(hù)區(qū)植被主要是溫帶典型草原,其建群種和優(yōu)勢(shì)種以白蓮蒿(Artemisiasacrorum)、大針茅(Stipagrandis)為主,伴生種類(lèi)型以豬毛蒿(A.scoparia)、百里香(Thymusmongolicus)、星毛委陵菜(Potentillaacaulis)為主。

    1.2 實(shí)驗(yàn)設(shè)計(jì)

    實(shí)驗(yàn)選取長(zhǎng)期封育和連續(xù)放牧兩種管理方式下的草地,采用隨機(jī)區(qū)組設(shè)計(jì)[3],自2013年起進(jìn)行長(zhǎng)期氮素添加實(shí)驗(yàn),于每年生長(zhǎng)季初期即4月底施肥,為提高肥料利用率于小雨天施肥,按實(shí)驗(yàn)設(shè)計(jì)的施用量用手撒施。其中封育草地封育30年;放牧草地連續(xù)30年進(jìn)行放牧實(shí)驗(yàn),放牧強(qiáng)度為4只羊/hm2,放牧期為6-8月,每月放牧1次,每次10 d,放牧?xí)r間為早上6點(diǎn)將家畜趕入放牧區(qū)自由采食,下午6點(diǎn)趕回畜圈休息。以尿素作為肥料,6個(gè)氮素添加水平依次為:CK,0 g/(m2·yr);N5,5 g/(m2·yr);N10,10 g/(m2·yr);N20,20 g/(m2·yr);N40,40 g/(m2·yr);N80,80 g/(m2·yr),等價(jià)于氮素分別為0、2.34、4.67、9.34、18.68、37.35 g/(m2·yr),其中N5為低氮水平,N10和N20為中氮水平,N40和 N80為高氮水平。每個(gè)水平6個(gè)重復(fù),每塊樣地上設(shè)置36個(gè)4 m×4 m小區(qū),各區(qū)間有2 m緩沖帶。

    1.3 觀測(cè)項(xiàng)目和測(cè)定方法

    于2015年在植物生長(zhǎng)旺期即8月上旬在每個(gè)小區(qū)內(nèi)隨機(jī)作0.5 m×0.5 m的樣方進(jìn)行植被群落調(diào)查,包括物種組成、數(shù)目、各物種植株平均高度、多度及蓋度。剪取植物地上部分并按物種分類(lèi)裝入信封,同時(shí)收集樣方內(nèi)凋落物帶回實(shí)驗(yàn)室80 ℃烘干至恒重,分別測(cè)定植物地上生物量與凋落物生物量。

    在天氣晴朗的日子用數(shù)字照度儀(LX-1010B,中國(guó))在上午11點(diǎn)至下午13點(diǎn)間測(cè)定群落光輻射。在每個(gè)小區(qū)內(nèi)隨機(jī)選取3處分別測(cè)定冠層光輻射和地面光輻射。

    功能群將具有相似生活型、特性或?qū)Νh(huán)境因子響應(yīng)類(lèi)似的植物進(jìn)行分組,基于功能群的研究方法是評(píng)估物種和群落對(duì)環(huán)境因子改變所做出響應(yīng)的有效手段[20]。參考王謝等[21]的方法按可食性進(jìn)行功能群劃分并計(jì)算各功能群地上生物量占總地上生物量的百分比以表征不同物種在群落中的功能。將植物群落以可食性劃分為4個(gè)組:可食禾草類(lèi)(grass species group,GG),包含莎草科干生苔草(Carexaridula)、非禾本科可食草類(lèi)(forbs species group,F(xiàn)G)、可食性低劣或有毒害草類(lèi)(noxious species group,NG)、豆科類(lèi)(leguminous species group,LG)。

    本研究計(jì)算物種重要值、α多樣性和植被透光率。α多樣性包括物種豐富度指數(shù)、Shannon-Wiener多樣性指數(shù)和Pielou均勻度指數(shù)。計(jì)算公式如下:

    (1) 物種重要值

    物種重要值=(相對(duì)生物量+相對(duì)蓋度+相對(duì)高度)×100/3

    (2) α多樣性指數(shù)

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

    Shannon-Wiener多樣性指數(shù):H′=-∑PilnPi

    Pielou均勻度指數(shù):J=H′/lnS′

    式中:S為樣方的平均物種數(shù);Pi為相對(duì)多度,即物種i的地上生物量占樣方內(nèi)所有物種地上生物量百分比;S′為各樣方中的物種數(shù)。

    (3) 植被透光率

    植被透光率=(地表光合有效輻射/冠層光合有效輻射)×100%

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

    利用R軟件進(jìn)行數(shù)據(jù)整理、分析及作圖。采用雙因素方差(Two-way ANOVA)分析草地管理方式和氮素添加對(duì)草地地上生物量和物種多樣性的影響,用Tukey法進(jìn)行多重比較,進(jìn)一步分析封育草地和放牧草地各施肥水平間及同一施肥水平下封育草地和放牧草地地上生物量與物種多樣性差異。采用Pearson法檢驗(yàn)物種多樣性與地上生物量間相關(guān)性,并做一般線性回歸。用Shapiro-Wilk test和Bartlett test檢驗(yàn)數(shù)據(jù)正態(tài)性和方差齊性。

    2 結(jié)果與分析

    2.1 封育草地和放牧草地地上生物量對(duì)氮素添加的響應(yīng)

    雙因素方差分析結(jié)果顯示(表1),草地管理方式和氮素添加對(duì)草地地上生物量皆有顯著影響,但其交互效應(yīng)并不顯著。與封育草地相比,放牧使地上生物量下降27.7%(236.72 g/m2到171.23 g/m2);施肥后封育草地地上生物量總體上呈低氮急速上升—下降—中氮緩慢上升—高氮下降的雙峰趨勢(shì),各氮素水平下地上生物量皆比未施肥(CK)時(shí)高,且在N5水平時(shí)達(dá)到最大,多重比較結(jié)果得出其與CK、N10、N80水平時(shí)地上生物量差異顯著(P<0.01);在放牧施肥草地,地上生物量的變化趨勢(shì)與封育施肥草地相似,與未施肥(CK)相比,氮素濃度為N5、 N10、 N80水平時(shí)地上生物量顯著增高(P<0.05),且在N10水平下地上生物量達(dá)到最大(圖1);兩樣地地上生物量在N5水平下差異顯著(P<0.05)。

    表1 草地管理方式和氮素添加對(duì)草地地上生物量和物種多樣性影響的雙因素方差分析

    N: 氮素添加Nitrogen addition; M: 草地管理方式Grassland management; 誤差方差自由度為60 There are 60 degrees of freedom for error.

    圖1 封育草地(F)和放牧草地(G)地上生物量對(duì)氮素添加的響應(yīng)(平均值±標(biāo)準(zhǔn)誤差)Fig.1 Responses of community aboveground biomass to nitrogen addition (mean±SE) CK、N5、N10、N20、N40、N80指氮素添加量分別為0、5、10、20、40、80 g CO(NH2)2/(m2·yr)。不同小寫(xiě)字母表示封育施肥草地和放牧施肥草地各施肥水平間差異顯著(P<0.05)。CK, N5, N10, N20, N40, N80 represent nitrogen addition levels of 0, 5, 10, 20, 40, 80 g CO(NH2)2/(m2·yr). Different small letters indicate significant differences among the nitrogen addition levels within the two sites (P<0.05).下同The same below.

    2.2 封育草地和放牧草地群落組分對(duì)氮素添加的響應(yīng)

    依照物種重要值>10%為優(yōu)勢(shì)物種,>5%為亞優(yōu)勢(shì)種,<5%為伴生種的標(biāo)準(zhǔn)劃分群落物種[4]。由圖2可見(jiàn),封育施肥草地的優(yōu)勢(shì)物種甘青針茅、大針茅、白蓮蒿等較高大的嗜氮植物重要值隨氮素水平增加而增加,個(gè)體高度較低的蓬子菜則呈降低趨勢(shì);個(gè)體高度較為中等的亞優(yōu)勢(shì)物種賴(lài)草、披堿草、火絨草、翼莖風(fēng)毛菊等重要值在低、中氮水平增加,在高氮水平下又降低,群落下層的甘菊、多毛并頭黃芩、百里香等低矮的物種重要值則表現(xiàn)為連續(xù)下降;重要值<5%的一些物種,如本氏針茅、冰草、蚓果芥等幾乎可看做偶見(jiàn)種,隨著氮素增加自群落中消失。

    在放牧施肥草地中情況大有不同,優(yōu)勢(shì)物種中不僅出現(xiàn)了光稃茅香,更值得注意的是出現(xiàn)了在封育施肥草地中作為伴生種的百里香和干生苔草;甘菊作為亞優(yōu)勢(shì)種出現(xiàn),與封育施肥草地不同的是其重要值隨氮素增加而連續(xù)增加,在封育施肥草地中為偶見(jiàn)種的甘肅蒿也在各氮素水平出現(xiàn);多種封育施肥草地中不存在的物種,如西山委陵菜、西藏點(diǎn)地梅、星毛委陵菜、青海苜蓿等個(gè)體高度很低的物種成為放牧施肥草地的伴生種。

    圖2 封育草地(F)和放牧草地(G)群落物種組成及重要值對(duì)氮素添加的響應(yīng)Fig.2 Responses of species composition and important value to nitrogen addition 按可食性分組的功能群,其中GG為可食禾草類(lèi)(包括莎草科)、FG為非禾本科可食草類(lèi)、NG為可食性低劣或有毒害草類(lèi)、LG為豆科類(lèi)。Functional group divided by edibility, GG (grass and sedge species group), FG (forbs species group), NG (noxious species group), LG (leguminous species group). 下同 The same below.1: 甘青針茅 Stipa przewalskyi;2: 大針茅 Stipa grandis; 3: 賴(lài)草 Leymus secalinus; 4: 干生苔草 Carex aridula; 5: 光稃茅香 Anthoxanthum glabrum; 6:披堿草 Elymus excelsus; 7: 本氏針茅 Stipa bungeana; 8: 早熟禾 Poa sphondylodes; 9: 冰草 Agropyron cristatum; 10: 甘肅蒿 Artemisia gansuensis; 11: 甘菊 Chrysanthemum lavandulifolium; 12: 火絨草 Leontopodium leontopodioides; 13: 阿爾泰狗娃花 Heteropappus altaicus; 14: 狼毒 Stellera chamaejasme; 15: 黃毛棘豆 Oxytropis ochrantha; 16: 青海苜蓿 Medicago archiducis-nicolaii; 17: 小果黃耆 Astragalus zacharensis; 18: 白蓮蒿 Artemisia sacrorum; 19: 百里香 Thymus mongolicus; 20: 蓬子菜 Viola dissecta; 21: 翼莖風(fēng)毛菊 Saussurea japonica var. pteroclada; 22: 絲毛飛廉 Carduus crispus; 23: 長(zhǎng)柱沙參 Adenophora stenanthina; 24: 細(xì)葉沙參 Adenophora capillaria subsp. paniculata; 25: 多毛并頭黃芩 Scutellaria scordifolia var. villosissima; 26: 二裂委陵菜 Potentilla bifurca; 27: 白花枝子花 Dracocephalum heterophyllum; 28: 蚓果芥 Neotorularia humilis; 29: 急折百蕊草 Thesium refractum; 30: 細(xì)葉遠(yuǎn)志 Polygala sibirica; 31: 野韭 Allium ramosum; 32: 北柴胡 Bupleurum chinense; 33: 西山委陵菜 Potentilla sischanensis; 34: 西藏點(diǎn)地梅 Androsace mariae; 35: 裂葉堇菜 Viola dissecta; 36: 星毛委陵菜 Potentilla acaulis; 37: 早開(kāi)堇菜 Viola prionantha; 38: 北方還陽(yáng)參 Crepis crocea.

    以可食性分組的功能群在兩樣地中皆以可食禾草(包括莎草科)為主,其地上生物量占地上總生物量比例平均值分別為59.36%和50.74%。在封育施肥草地中,可食禾草地上生物量所占比重隨氮素水平的變化增加不明顯,甚至在N20水平下降22%(圖3),結(jié)合圖2可知處于優(yōu)勢(shì)地位的甘青針茅在N20水平下大幅的下降造成了這種結(jié)果;非禾本科可食草類(lèi)和可食性低劣或有毒害草類(lèi)地上生物量皆呈先增加再降低的趨勢(shì),氮素水平間差異不顯著;封育施肥草地中未出現(xiàn)豆科類(lèi)草。放牧施肥草地中,可食禾草地上生物量所占比重隨著氮素水平增加連續(xù)增加,在N80時(shí)與CK形成顯著差異(P<0.01),增加1.39倍,除了在封育施肥草地中也存在的大針茅、甘青針茅等禾本科草類(lèi),干生苔草、本氏針茅、光稃茅香和早熟禾地上生物量亦明顯增加(圖2);非禾本科可食草類(lèi)的比重則呈降低趨勢(shì),CK、N5、N20水平下的地上生物量占比皆與N80差異顯著(P<0.01),分別下降63.13%、65.58%、54.38%;可食性低劣或有毒害草類(lèi)地上生物量占比隨氮素增加而降低,但各處理間差異不顯著;與重要值分析相對(duì)應(yīng),放牧施肥草地出現(xiàn)了可食性?xún)?yōu)良的豆科牧草(圖3)。

    2.3 封育草地和放牧草地物種多樣性對(duì)氮素添加的響應(yīng)

    雙因素方差分析結(jié)果顯示(表1),草地管理方式對(duì)3種多樣性指數(shù)都具有顯著影響,氮素添加對(duì)物種豐富度具有顯著影響,但對(duì)Shannon-Wiener多樣性指數(shù)和Pielou均勻度指數(shù)無(wú)顯著影響,且交互作用對(duì)3種指數(shù)亦無(wú)顯著影響。由圖4可見(jiàn),3種多樣性指數(shù)在封育草地中隨施肥水平的遞增皆呈連續(xù)下降趨勢(shì),但各處理間差異不顯著,與CK相比,N80水平下的物種豐富度指數(shù)、Shannon-Wiener多樣性指數(shù)和Pielou均勻度指數(shù)分別下降31.3%、32.64%、22.58%;在放牧施肥草地中,3種指數(shù)皆為先增加后降低,其中物種豐富度指數(shù)和Shannon-Wiener多樣性指數(shù)在N10水平時(shí)達(dá)到最大,均勻度指數(shù)在N20水平時(shí)升至峰值。兩樣地的3種多樣性指數(shù)分別在N10、N80、N80水平下差異顯著(P<0.05)。

    圖3 封育草地(F)和放牧草地(G)各功能群地上生物量對(duì)氮素添加的響應(yīng)Fig.3 Responses of aboveground biomass of functional group to nitrogen addition

    圖4 封育草地(F)和放牧草地(G)α多樣性對(duì)氮素添加的響應(yīng)Fig.4 Responses of alpha diversity to nitrogen addition R、S、J分別代表物種豐富度、Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù)。* 表示封育草地與放牧草地在同一處理下差異顯著(P<0.05)。 R, S, J represent species richness, Shannon-Wiener index and Pielou index. * indicates significant difference among fenced and grazed plots under the same treatment (P<0.05).

    2.4 封育施肥草地和放牧施肥草地物種多樣性與地上生物量的關(guān)系

    封育施肥草地的Shannon-Wiener多樣性指數(shù)和地上生物量呈顯著負(fù)相關(guān)(r=-0.47,P<0.01),在放牧施肥草地兩變量并無(wú)顯著相關(guān)性(圖5)。

    相關(guān)分析結(jié)果顯示,Shannon-Wiener多樣性指數(shù)與凋落物生物量呈顯著負(fù)相關(guān)(r=-0.51,P<0.001),隨凋落物生物量線性降低(圖6);Shannon-Wiener多樣性指數(shù)與植被透光率呈顯著正相關(guān)(r=0.47,P<0.001),隨植被透光率線性增加(圖6)。

    圖5 氮素添加后封育草地(F)和放牧草地(G)物種多樣性指數(shù)與地上生物量之間的關(guān)系Fig.5 Relationship between biodiversity index and aboveground biomass after nitrogen addiction 以Shannon-Wiener多樣性指數(shù)為解釋變量,地上生物量為響應(yīng)變量建立回歸模型。Regression model was established by using Shannon-Wiener index as explanatory variable and aboveground biomass as response variable. 顯著性水平為Significant differences are reported as: NS, P>0.05;*,P<0.05;**,P<0.01;***,P<0.001;下同 The same below.

    圖6 氮素添加后物種多樣性指數(shù)與凋落物生物量和植被透光率之間的關(guān)系Fig.6 Relationship between biodiversity index and aboveground biomass, vegetation light penetration 以Shannon-Wiener多樣性指數(shù)為響應(yīng)變量,凋落物生物量、植被透光率為解釋變量建立回歸模型。Regression models were established by using Shannon-Wiener index as response variable and litter aboveground biomass, vegetation light penetration as explanatory variables.

    3 討論

    3.1 封育草地和放牧草地地上生物量對(duì)氮素添加的響應(yīng)

    添加氮素后,封育和放牧草地地上生物量與對(duì)照相比皆顯著增加,變化趨勢(shì)大體上都為先增加后降低,但二者在數(shù)量上有一定的差距(圖1),這與Bai等[3]的研究結(jié)果一致,他認(rèn)為兩樣地間生物量的差異源自于不同的植物-土壤反饋系統(tǒng)影響了植物對(duì)水的獲取能力。在濕潤(rùn)的高寒草甸和干旱的荒漠草原進(jìn)行氮素添加實(shí)驗(yàn),結(jié)果顯示兩群落的地上生物量分別在氮素濃度為23.34和0.35 g/m2時(shí)達(dá)到峰值[4,22],表明對(duì)位處干旱半干旱區(qū)的草地生態(tài)系統(tǒng)來(lái)說(shuō),除了氮素,水分也是植物生物量增加的限制因素之一[18]。陳靜等[23]在科爾沁沙質(zhì)草地實(shí)施的增雨和增氮實(shí)驗(yàn)發(fā)現(xiàn),單純地添加水分或氮素對(duì)植被總生物量及各構(gòu)件平均生物量無(wú)顯著影響,但在夏季增雨×氮素添加處理下生物量則顯著增加。

    在本實(shí)驗(yàn)中,封育施肥和放牧施肥草地地上生物量達(dá)到最高值時(shí)的氮素濃度分別為2.34和4.67 g/m2,放牧施肥草地地上生物量到達(dá)峰值時(shí)的氮素濃度比封育施肥草地高,表明封育草地地上生物量對(duì)氮素添加的響應(yīng)更加敏感。封育草地中,氮素添加前處于優(yōu)勢(shì)地位的甘青針茅、大針茅、白蓮蒿等個(gè)體高度較高且根系發(fā)達(dá)的物種在添加氮素后依靠其生理特性可獲得更多的營(yíng)養(yǎng)和水分來(lái)生長(zhǎng)繁殖,其凋亡后成為營(yíng)養(yǎng)豐富的凋落物,一方面為植物、土壤微生物提供營(yíng)養(yǎng),另一方面覆蓋于地面減少水分散失,起到一定的保水作用,形成良性的水—植物—土壤反饋系統(tǒng),使地上生物量在氮素添加后快速增加。而放牧施肥草地則由于牲畜的啃食、踩踏凋落物積累量較少,失去了保水和提供營(yíng)養(yǎng)的凋落物層,地上生物量的增長(zhǎng)較封育施肥草地緩慢,體現(xiàn)在與封育施肥草地地上生物量在數(shù)量上有一定差距且對(duì)氮素添加的響應(yīng)較慢。

    3.2 封育草地和放牧草地群落物種組成和物種多樣性對(duì)氮素添加的響應(yīng)

    伴隨氮素梯度增加帶來(lái)的地上生物量增加同時(shí)也導(dǎo)致了凋落物生物量的積累,而這可能是造成物種多樣性降低的原因之一。Jensen等[24]移除了樣方內(nèi)的凋落物,結(jié)果種子出苗率顯著增加,提高了物種豐富度并改變了物種組成,表明凋落物在植物群落構(gòu)建中起著關(guān)鍵性作用。Clark等[7]的凋落物移除實(shí)驗(yàn)顯示,群落中到達(dá)土壤表面的透射光增加了54%,幼苗數(shù)量增加,最終導(dǎo)致物種豐富度增加。根據(jù)以上及前人實(shí)驗(yàn),目前關(guān)于凋落物影響植物群落物種組成和物種多樣性的機(jī)制有以下幾點(diǎn):1)阻礙種子出苗或幼苗建植(包括改變種子的物理環(huán)境導(dǎo)致萌發(fā)的關(guān)鍵因素缺失;分泌化學(xué)物質(zhì)改變化學(xué)環(huán)境抑制種子萌發(fā);形成機(jī)械屏障阻礙種子擴(kuò)散);2)削弱到達(dá)地面的光照強(qiáng)度,影響幼苗或個(gè)體高度較低的植物獲取光能;3)占據(jù)裸露地面抑制了克隆植物無(wú)性繁殖[25]。

    本實(shí)驗(yàn)封育施肥草地中,隨著氮素增加后凋落物的積累,從功能群層面看,禾本科草類(lèi)受凋落物影響較小,豆科草類(lèi)受影響最大。從物種組成層面看,一些個(gè)體高度較低的物種如蓬子菜、多毛并頭黃芩、百里香等重要值隨氮素增加而降低甚至從樣方中消失(圖2),導(dǎo)致了物種多樣性連續(xù)降低(圖4);放牧施肥草地中則由于牲畜的啃食與踩踏行為減小了凋落物積累量,解除了凋落物對(duì)一些較低矮物種的限制,出現(xiàn)了早熟禾、西藏點(diǎn)地梅、青海苜蓿等物種,但氮素輸入導(dǎo)致了禾本科植物地上生物量連續(xù)增高(圖3),凋落物也隨之增高,從而造成放牧施肥草地物種多樣性先升高后下降的趨勢(shì),但其物種多樣性總體水平依然高于封育施肥草地。

    氮素的增加解除了生態(tài)系統(tǒng)對(duì)氮素的限制,封育施肥草地中對(duì)氮素敏感、根系密度高、處于優(yōu)勢(shì)地位的嗜氮物種,如禾本科的賴(lài)草、大針茅,非禾本科可食草類(lèi)的白蓮蒿快速生長(zhǎng),進(jìn)一步在養(yǎng)分吸收、高度、蓋度上占得優(yōu)勢(shì),降低植被透光率,蔭蔽一些對(duì)氮素不敏感生長(zhǎng)較慢或矮小的物種,從而降低了物種多樣性;而放牧施肥草地較高的植被透光率利于較低矮植物的拓殖和存活,一些封育施肥草地中不存在的豆科牧草如青海苜蓿和小果黃耆也相繼出現(xiàn),故隨著植被透光率增加,更多的光能到達(dá)地表,物種多樣性也隨之增加(圖5)。表明種間競(jìng)爭(zhēng)由植物對(duì)地下養(yǎng)分的競(jìng)爭(zhēng)轉(zhuǎn)為地上對(duì)光資源的競(jìng)爭(zhēng),植物對(duì)養(yǎng)分和光競(jìng)爭(zhēng)力的差異造成的生態(tài)位分化決定了最終群落物種組成格局[11,26]。

    故結(jié)合光競(jìng)爭(zhēng)和凋落物積累對(duì)群落物種組成和物種多樣性的影響可見(jiàn),氮素添加后封育施肥草地可能由于凋落物積累和植被對(duì)光資源競(jìng)爭(zhēng)強(qiáng)度的增強(qiáng)導(dǎo)致了物種多樣性下降,而放牧施肥草地凋落物積累量較少,也沒(méi)有改變光競(jìng)爭(zhēng)強(qiáng)度,其物種多樣性下降速度低于封育施肥草地物種,故推論:放牧可減緩施肥對(duì)物種多樣性的降低作用。更有研究指出適度的干擾可以提高物種豐富度,改善種群的生態(tài)位配置及最大限度地利用環(huán)境資源[27]。

    3.3 封育和放牧草地地上生物量和物種多樣性的關(guān)系

    張杰琦等[22]和于麗等[28]在封育草地添加氮素的實(shí)驗(yàn)結(jié)果顯示,物種豐富度和Shannon-Wiener多樣性指數(shù)、Simpson優(yōu)勢(shì)度指數(shù)、Pielou均勻度指數(shù)皆與地上生物量呈顯著的線性負(fù)相關(guān)關(guān)系。與其相同,本實(shí)驗(yàn)封育施肥草地中,Shannon-Wiener多樣性指數(shù)與地上生物量存在顯著負(fù)相關(guān)關(guān)系;在放牧施肥草地中,Shannon-Wiener多樣性指數(shù)與地上生物量無(wú)顯著的相關(guān)關(guān)系(圖6),但Shannon-Wiener多樣性指數(shù)隨地上上生物量增加總體上呈下降趨勢(shì),且比封育施肥草地緩慢。這一結(jié)果與Bai等[3]的結(jié)論一致。

    自20世紀(jì)60年代起,物種多樣性與生產(chǎn)力的關(guān)系就是植物生態(tài)學(xué)家們極具興趣且充滿(mǎn)爭(zhēng)議性的話(huà)題。一些學(xué)者認(rèn)為物種多樣性與生產(chǎn)力之間的關(guān)系隨著研究尺度的變化有著較為普遍的模式,并利用Meta分析總結(jié)出5種被廣泛接受的物種豐富度—生產(chǎn)力關(guān)系(SRPR):正相關(guān),負(fù)相關(guān),單峰,U型和不相關(guān)[29]。研究者們認(rèn)為模式的多變性取決于研究尺度,即SPRP是“尺度依賴(lài)性”的[30]。Meta分析的結(jié)果顯示,在較大的尺度上這種關(guān)系更傾向于表現(xiàn)為正相關(guān),而在較小的尺度上單峰與U型模式則更加普遍。但也有學(xué)者對(duì)這種綜合性的結(jié)果持質(zhì)疑態(tài)度,Adler等[31]通過(guò)在五大洲的標(biāo)準(zhǔn)化取樣方法分別在局域、區(qū)域和全球尺度探索物種豐富度和初級(jí)生產(chǎn)力之間的關(guān)系,結(jié)果并未找到明確的、具有普遍性的關(guān)系模式,認(rèn)為Meta分析忽略了不同研究間方法論的差異,不能夠作為可靠的結(jié)論,生產(chǎn)力并不能作為一個(gè)很好的指標(biāo)來(lái)預(yù)測(cè)物種多樣性的變化,即在局域、區(qū)域和全球尺度不存在一致的生產(chǎn)力與物種多樣性關(guān)系。

    本實(shí)驗(yàn)亦未得出物種多樣性與地上生物量間特定的模型。因?yàn)樵趶?fù)雜多變的生態(tài)環(huán)境中,影響物種多樣性的因素是多元化的,在后續(xù)的研究中應(yīng)綜合分析降雨量,土壤營(yíng)養(yǎng)元素、pH、含水量,微生物活動(dòng),溫度等因子。

    4 結(jié)論

    通過(guò)對(duì)黃土高原半干旱區(qū)長(zhǎng)期封育和連續(xù)放牧草地進(jìn)行氮素添加實(shí)驗(yàn),研究結(jié)果表明:1)施肥顯著提高了封育草地地上生物量,但當(dāng)?shù)厣仙锪窟_(dá)到一定閾值時(shí)也導(dǎo)致了物種多樣性降低,改變?nèi)郝湮锓N組成,這些改變可能會(huì)反過(guò)來(lái)影響地上生物量;2)放牧可在顯著提高地上生物量的同時(shí)減小物種多樣性降低程度。作為生態(tài)系統(tǒng)穩(wěn)定性的核心指標(biāo),較高的地上生物量與高物種多樣性、高穩(wěn)定性是相伴而生的[32]。所以適當(dāng)?shù)姆拍潦呛侠淼奶烊徊莸毓芾矸绞?,這也符合中度干擾能增加物種多樣性的干擾理論。3)未發(fā)現(xiàn)物種多樣性與地上生物量間特定的關(guān)系模型,在以后的研究中應(yīng)加入更加多元化的因子探究各指標(biāo)對(duì)物種多樣性差異的解釋率。另外,相對(duì)于其他類(lèi)型草地(如高寒草甸、荒漠草原等),黃土高原典型草原對(duì)氮素添加的響應(yīng)是否更加緩慢或不敏感,仍需進(jìn)行長(zhǎng)期研究,可從物種個(gè)體功能性狀響應(yīng)等方面來(lái)進(jìn)一步解釋其內(nèi)在機(jī)制。

    References:

    [1] Li L J, Zeng D H, Yu Z Y,etal. Effects of nitrogen addition on grassland species diversity and productivity in Keerqin sandy land. Chinese Journal of Applied Ecology, 2009, 20(8): 1838-1844. 李祿軍, 曾德慧, 于占源, 等. 氮素添加對(duì)科爾沁沙質(zhì)草地物種多樣性和生產(chǎn)力的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2009, 20(8): 1838-1844.

    [2] Stevens C J, Duprè C, Dorland E,etal. Nitrogen deposition threatens species richness of grasslands across Europe. Environmental Pollution, 2010, 158(9): 2940-2945.

    [3] Bai Y, Wu J, Clark C M,etal. Tradeoffs 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.

    [4] He Y H, Liu X P, Xie Z K. Effect of nitrogen addition on species diversity and plant productivity of herbaceous plants in desert grassland of the Loess Plateau. Journal of Desert Research, 2015, 35(1): 66-71. 何玉惠, 劉新平, 謝忠奎. 氮素添加對(duì)黃土高原荒漠草原草本植物物種多樣性和生產(chǎn)力的影響. 中國(guó)沙漠, 2015, 35(1): 66-71.

    [5] Stevens C J, Dise N B, Mountford J O,etal. Impact of nitrogen deposition on the species richness of grasslands. Science, 2004, 303: 1876-1879.

    [6] Hautier Y, Niklaus P A, Hector A. Competition for light causes plant biodiversity loss after eutrophication. Science, 2009, 324: 636-638.

    [7] Clark C M, Tilman D. Recovery of plant diversity following N cessation: effects of recruitment, litter, and elevated N cycling. Ecology, 2010, 91(12): 3620-3630.

    [8] Stevens M H H, Carson W P. Plant density determines species richness along an experimental fertility gradient. Ecology, 1999, 80(2): 455-465.

    [9] Stevens C J, Thompson K, Grime J P,etal. Contribution of acidification and eutrophication to declines in species richness of calcifuge grasslands along a gradient of atmospheric nitrogen deposition. Functional Ecology, 2010, 24(2): 478-484.

    [10] Van Den Berg L J L, Dorland E, Vergeer P,etal. Decline of acid-sensitive plant species in heathland can be attributed to ammonium toxicity in combination with low pH. New Phytologist, 2005, 166(2): 551-564.

    [11] Newman E. Competition and diversity in herbaceous vegetation. Nature, 1973, 244: 310.

    [12] Tilman D. Species richness of experimental productivity gradients: how important is colonization limitation. Ecology, 1993, 74(8): 2179-2191.

    [13] Bowman W D, Cleveland C C, Halada,etal. Negative impact of nitrogen deposition on soil buffering capacity. Nature Geoscience, 2008, 1(11): 767-770.

    [14] Clark C M, Tilman D. Loss of plant species after chronic low-level nitrogen deposition to prairie grasslands. Nature, 2008, 451: 712-715.

    [15] Tilman D. Secondary succession and the pattern of plant dominance along experimental nitrogen gradients. Ecological Monographs, 1987, 57(3): 189-214.

    [16] Li D J, Mo J M, Fang Y T,etal. Impact of nitrogen deposition on forest plants. Acta Ecologica Sinica, 2003, 23(9): 1891-1900. 李德軍, 莫江明, 方運(yùn)霆, 等. 氮沉降對(duì)森林植物的影響. 生態(tài)學(xué)報(bào), 2003, 23(9): 1891-1900.

    [17] Zheng L X, Liu X J, Zhang F S. Atomosphere deposition of organic nitrogen: A review. Acta Ecologica Sinica, 2007, 27(9): 3828-3834. 鄭利霞, 劉學(xué)軍, 張福鎖. 大氣有機(jī)氮沉降研究進(jìn)展. 生態(tài)學(xué)報(bào), 2007, 27(9): 3828-3834.

    [18] Zhou X B, Zhang Y M. Review on the ecological of N deposition in arid and semi-arid areas. Acta Ecologica Sinica, 2009, 29(7): 3835-3845. 周曉兵, 張?jiān)? 干旱半干旱區(qū)氮沉降生態(tài)效應(yīng)研究進(jìn)展. 生態(tài)學(xué)報(bào), 2009, 29(7): 3835-3845.

    [19] Yang Z L, Su F L, Miao Y,etal. Effects of fertilization and grazing on species richness in an alpine meadow of Qinghai-Xizang Plateau. Chinese Journal of Plant Ecology, 2014, 38(10): 1074-1081. 楊中領(lǐng), 蘇芳龍, 苗原, 等. 施肥和放牧對(duì)青藏高原高寒草甸物種豐富度的影響. 植物生態(tài)學(xué)報(bào), 2014, 38(10): 1074-1081.

    [20] Li W, Wen S, Hu W,etal. Root-shoot competition interactions cause diversity loss after fertilization: a field experiment in an alpine meadow on the Tibetan Plateau. Journal of Plant Ecology, 2011, 4(3): 138-146.

    [21] Wang X, Xiang C H, Li X W,etal. Effects of a winter wildfire on plant community structure and forage quality in subalpine grassland of western Sichuan, China. Chinese Journal of Plant Ecology, 2013, 37(10): 922-932. 王謝, 向成華, 李賢偉, 等. 冬季火對(duì)川西亞高山草地植物群落結(jié)構(gòu)和牧草質(zhì)量的影響. 植物生態(tài)學(xué)報(bào), 2013, 37(10): 922-932.

    [22] Zhang J Q, Li Q, Ren Z W,etal. Effects of nitrogen addition on species richness and relationship between species richness and aboveground productivity of alpine meadow of the Qinghai-Tibetan Plateau, China. Chinese Journal of Plant Ecology, 2010, 34(10): 1125-1131. 張杰琦, 李奇, 任正煒, 等. 氮素添加對(duì)青藏高原高寒草甸植物群落物種豐富度及其與地上生產(chǎn)力關(guān)系的影響. 植物生態(tài)學(xué)報(bào), 2010, 34(10): 1125-1131.

    [23] Chen J, Li Y L, Cui D,etal. Response of aboveground biomass allocation in four dominant species to water and nitrogen addition in the Horqin sandy land. Journal of Desert Research, 2014, 34(3): 696-703. 陳靜, 李玉霖, 崔奪, 等. 氮素及水分添加對(duì)科爾沁沙地4種優(yōu)勢(shì)植物地上生物量分配的影響. 中國(guó)沙漠, 2014, 34(3): 696-703.

    [24] Jensen K, Meyer C. Effects of light competition and litter on the performance ofViolapalustrisand on species composition and diversity of an abandoned fen meadow. Plant Ecology, 2001, 155(2): 169-181.

    [25] Bosy J L, Reader R J. Mechanisms underlying the suppression of forb seedling emergence by grass (Poapratensis) litter. Functional Ecology, 1995, 9(4): 635-639.

    [26] Li L J, Yu Z Y, Zeng D H,etal. Effects of fertilizations on species composition and diversity of grassland in Keerqin Sandy Lands. Acta Prataculturae Sinica, 2010, 19(2): 109-115. 李祿軍, 于占源, 曾德慧, 等. 施肥對(duì)科爾沁沙質(zhì)草地群落物種組成和多樣性的影響. 草業(yè)學(xué)報(bào), 2010, 19(2): 109-115.

    [27] Jing G H, Cheng J M, Su J S,etal. Response of dominant niche breadths and niche overlaps to various disturbance factors in typical steppe fenced grassland of China’s Loess Plateau region. Acta Prataculturae Sinica, 2015, 24(9): 43-52. 井光花, 程積民, 蘇紀(jì)帥, 等. 黃土區(qū)長(zhǎng)期封育草地優(yōu)勢(shì)物種生態(tài)位寬度與生態(tài)位重疊對(duì)不同干擾的響應(yīng)特征. 草業(yè)學(xué)報(bào), 2015, 24(9): 43-52.

    [28] Yu L, Zhao J N, Wang H,etal. Effects of nutrient addition on plant diversity and productivity in aStipabaicalensisgrassland in Inner Mongolia, China. Acta Ecologica Sinica, 2015, 35(24): 8165-8173. 于麗, 趙建寧, 王慧, 等. 養(yǎng)分添加對(duì)內(nèi)蒙古貝加爾針茅草原植物多樣性與生產(chǎn)力的影響. 生態(tài)學(xué)報(bào), 2015, 35(24): 8165-8173.

    [29] Whittaker R J. Meta-analyses and mega-mistakes: calling time on meta-analysis of the species richness-productivity relationship. Ecology, 2010, 91(9): 2522-2533.

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

    [31] Adler P B, Seabloom E W, Borer E T,etal. Productivity is a poor predictor of plant species richness. Science, 2011, 333: 1750-1753.

    [32] Bai Y F, Li L H, Huang J H,etal. The influence of plant diversity and functional composition on ecosystem stability of fourStipacommunities in the Inner Mongolia Plateau. Acta Botanica Sinica, 2001, 43(3): 280-287.

    Responses of species diversity and aboveground biomass to nitrogen addition in fenced and grazed grassland on the Loess Plateau

    ZHAO Jie1, LI Wei2, JING Guang-Hua2, WEI Lin2, CHENG Ji-Min1,2*

    1.CollegeofAnimalScienceandTechnology,NorthwestA&FUniversity,Yangling712100,China; 2.InstituteofSoilandWaterConservation,ChineseAcademyofSciences&MinistryofWaterResource,Yangling712100,China

    Nitrogen addition can increase the aboveground biomass of the plant community in grasslands. However, it can also lead to a loss of species diversity, but the mechanisms involved in this process are unclear. The objective of this study was to explore the effects of nitrogen addition on aboveground biomass, species diversity, community composition, and the relationship between species richness and aboveground biomass in long-term fenced grassland and grazed grassland on the semi-arid Loess Plateau. We focused on the differences and similarities in the grassland responses under these two typical management systems, and tried to identify the the ecological mechanisms underlying these responses. We added nitrogen at different rates [0, 5, 10, 20, 40, 80 g CO(NH2)2/(m2·yr)] to long-term fenced grassland and grazed grassland. After 3 years of treatment, we measured aboveground biomass, calculated species diversity indexes (species richness, Shannon-Wiener index, and Pielou index) and analyzed the relationship between species diversity and productivity using regression analyses. The results showed that, in the long-term fenced grassland, nitrogen addition significantly increased the aboveground biomass and decreased species diversity, while the grazed grassland showed similar but weaker responses to nitrogen addition. The species diversity in fenced and fertilized grassland decreased as the litter biomass increased, and also increased with greater light penetration. In comparison, species diversity in grazed and fertilized grassland decreased slowly and did not change significantly with greater light penetration. The Shannon-Wiener index was significantly negatively correlated with above-ground productivity in fenced and fertilized grassland, but no such relationship existed in grazed grassland. These results indicated that grazing can relieve the negative effect of fertilization on species diversity. Therefore, grazing is a reasonable grassland management method to retain species diversity while increasing the community aboveground biomass by fertilization.

    grazing; fencing; nitrogen addition; species diversity; aboveground biomass; community composition

    10.11686/cyxb2017064

    http://cyxb.lzu.edu.cn

    2017-02-24;改回日期:2017-04-17

    國(guó)家重點(diǎn)專(zhuān)項(xiàng)(2016YFC0500700),國(guó)家自然科學(xué)基金(41601586)和陜西省自然科學(xué)基金(2016JQ3008)資助。

    趙潔(1992-),女,陜西西安人,碩士。E-mail: 1440039188@qq.com*通信作者Corresponding author. E-mail: gyzcjm@ms.iswc.ac.cn

    趙潔, 李偉, 井光花, 魏琳, 程積民. 黃土區(qū)封育和放牧草地物種多樣性和地上生物量對(duì)氮素添加的響應(yīng). 草業(yè)學(xué)報(bào), 2017, 26(8): 54-64.

    ZHAO Jie, LI Wei, JING Guang-Hua, WEI Lin, CHENG Ji-Min. Responses of species diversity and aboveground biomass to nitrogen addition in fenced and grazed grassland on the Loess Plateau. Acta Prataculturae Sinica, 2017, 26(8): 54-64.

    猜你喜歡
    物種
    物種大偵探
    物種大偵探
    物種大偵探
    吃光入侵物種真的是解決之道嗎?
    生日禮物種草合集
    物種大滅絕
    麗水發(fā)現(xiàn)新物種
    誰(shuí)在“摧毀”澳大利亞——可怕的物種入侵
    一億年后,地球上可能出現(xiàn)哪些新物種
    回首2018,這些新物種值得關(guān)注
    亚洲精华国产精华精| 免费在线观看影片大全网站| 欧美黄色淫秽网站| 岛国在线免费视频观看| 久久亚洲精品不卡| 精品一区二区三区视频在线观看免费| 国产黄片美女视频| 蜜桃久久精品国产亚洲av| 床上黄色一级片| 999久久久国产精品视频| xxx96com| 国产麻豆成人av免费视频| 日韩av在线大香蕉| 欧美一区二区国产精品久久精品| 他把我摸到了高潮在线观看| 国产精品香港三级国产av潘金莲| 日本熟妇午夜| 欧美一级a爱片免费观看看| 国产成人av教育| 精品久久久久久久末码| 19禁男女啪啪无遮挡网站| 亚洲国产精品999在线| 国产午夜福利久久久久久| 亚洲一区二区三区不卡视频| 国产av一区在线观看免费| 亚洲成人久久爱视频| 午夜亚洲福利在线播放| 精品久久久久久久末码| 90打野战视频偷拍视频| 九九久久精品国产亚洲av麻豆 | 少妇熟女aⅴ在线视频| 欧美大码av| 成年免费大片在线观看| 老司机福利观看| 国产97色在线日韩免费| 日韩欧美免费精品| 一进一出好大好爽视频| 久久久久亚洲av毛片大全| 九九在线视频观看精品| 国产亚洲精品av在线| 中文在线观看免费www的网站| 一卡2卡三卡四卡精品乱码亚洲| 在线看三级毛片| av女优亚洲男人天堂 | 中文字幕久久专区| 欧美国产日韩亚洲一区| 国产日本99.免费观看| 国产亚洲精品久久久com| 男插女下体视频免费在线播放| 久久久久久久精品吃奶| 亚洲九九香蕉| 床上黄色一级片| 免费电影在线观看免费观看| 国产真人三级小视频在线观看| 悠悠久久av| 国产亚洲精品久久久久久毛片| 激情在线观看视频在线高清| 一级毛片精品| 中文在线观看免费www的网站| 亚洲av电影不卡..在线观看| 在线看三级毛片| 免费在线观看亚洲国产| 国产97色在线日韩免费| 亚洲欧美激情综合另类| 在线免费观看不下载黄p国产 | a级毛片a级免费在线| 男女下面进入的视频免费午夜| 亚洲自拍偷在线| 搡老岳熟女国产| 91麻豆精品激情在线观看国产| 18禁国产床啪视频网站| 丁香欧美五月| 久久久成人免费电影| 国产精品一及| 亚洲真实伦在线观看| 精品国产亚洲在线| 精品福利观看| 色老头精品视频在线观看| 日本熟妇午夜| 日本熟妇午夜| 亚洲最大成人中文| 国产精品国产高清国产av| 亚洲精品一卡2卡三卡4卡5卡| 免费一级毛片在线播放高清视频| 久久精品人妻少妇| 久久香蕉精品热| 国产伦一二天堂av在线观看| 少妇的逼水好多| 国产精品一区二区免费欧美| 亚洲在线观看片| 亚洲avbb在线观看| 国产成人啪精品午夜网站| 桃色一区二区三区在线观看| 成年女人毛片免费观看观看9| 在线免费观看不下载黄p国产 | 99久国产av精品| 成人鲁丝片一二三区免费| 男人的好看免费观看在线视频| 中文字幕熟女人妻在线| 午夜免费观看网址| 国产男靠女视频免费网站| 欧美三级亚洲精品| 精品福利观看| 岛国在线观看网站| 亚洲午夜精品一区,二区,三区| 黑人巨大精品欧美一区二区mp4| 黄色成人免费大全| 亚洲精品一卡2卡三卡4卡5卡| 国产乱人伦免费视频| 网址你懂的国产日韩在线| 超碰成人久久| 亚洲成av人片在线播放无| 99久久精品一区二区三区| 啦啦啦韩国在线观看视频| 99热这里只有是精品50| 精品久久久久久久久久久久久| 亚洲国产中文字幕在线视频| 精品日产1卡2卡| 免费看日本二区| 免费看日本二区| 成熟少妇高潮喷水视频| 欧美另类亚洲清纯唯美| 狂野欧美激情性xxxx| 美女扒开内裤让男人捅视频| 亚洲欧美日韩高清在线视频| 久久国产精品人妻蜜桃| 国产aⅴ精品一区二区三区波| 99热精品在线国产| 天堂网av新在线| 国产精品野战在线观看| 亚洲专区字幕在线| 午夜日韩欧美国产| 日本在线视频免费播放| 国产亚洲欧美98| bbb黄色大片| 性欧美人与动物交配| aaaaa片日本免费| 狂野欧美激情性xxxx| 午夜成年电影在线免费观看| 国产欧美日韩精品一区二区| 丰满人妻熟妇乱又伦精品不卡| 一级毛片高清免费大全| 在线国产一区二区在线| 欧美高清成人免费视频www| 成人永久免费在线观看视频| 一区二区三区国产精品乱码| 久久久久久久久免费视频了| 久久国产精品人妻蜜桃| 久久久久久久久中文| 熟妇人妻久久中文字幕3abv| 国产精品国产高清国产av| 亚洲人成网站在线播放欧美日韩| 男人舔奶头视频| 好男人电影高清在线观看| 欧美又色又爽又黄视频| 亚洲成人精品中文字幕电影| 在线观看一区二区三区| 日本 av在线| 久久久成人免费电影| 久久精品国产亚洲av香蕉五月| 中文字幕人成人乱码亚洲影| 婷婷六月久久综合丁香| 在线免费观看的www视频| 精品久久久久久久毛片微露脸| 欧美在线一区亚洲| tocl精华| 99久久精品国产亚洲精品| 精品日产1卡2卡| 亚洲av成人精品一区久久| 在线国产一区二区在线| 亚洲最大成人中文| 国产亚洲精品久久久久久毛片| 男人舔女人的私密视频| 男人舔女人下体高潮全视频| 国产精品久久久久久精品电影| www.999成人在线观看| 长腿黑丝高跟| 黄色视频,在线免费观看| 色综合站精品国产| 极品教师在线免费播放| 欧美在线一区亚洲| 中文在线观看免费www的网站| 亚洲人成伊人成综合网2020| 精品午夜福利视频在线观看一区| 国产aⅴ精品一区二区三区波| 小蜜桃在线观看免费完整版高清| 亚洲av五月六月丁香网| 日本免费一区二区三区高清不卡| 999久久久精品免费观看国产| 大型黄色视频在线免费观看| 国产三级在线视频| 热99在线观看视频| www.999成人在线观看| 午夜福利免费观看在线| 51午夜福利影视在线观看| 中亚洲国语对白在线视频| 久久九九热精品免费| 熟女人妻精品中文字幕| 久久久精品大字幕| 九九在线视频观看精品| 伊人久久大香线蕉亚洲五| 亚洲成人精品中文字幕电影| 黄色女人牲交| 女生性感内裤真人,穿戴方法视频| 成在线人永久免费视频| 亚洲七黄色美女视频| 嫁个100分男人电影在线观看| 三级毛片av免费| av黄色大香蕉| 好看av亚洲va欧美ⅴa在| 在线a可以看的网站| 久久久久亚洲av毛片大全| 色吧在线观看| 久久久国产精品麻豆| 国产成人影院久久av| 首页视频小说图片口味搜索| 国产高清videossex| 欧美绝顶高潮抽搐喷水| 精品一区二区三区视频在线观看免费| 亚洲成人精品中文字幕电影| 亚洲国产欧洲综合997久久,| 99国产精品一区二区蜜桃av| 精品欧美国产一区二区三| 欧美丝袜亚洲另类 | 国产高清有码在线观看视频| 一本一本综合久久| 婷婷丁香在线五月| 香蕉久久夜色| 大型黄色视频在线免费观看| 19禁男女啪啪无遮挡网站| 免费人成视频x8x8入口观看| 午夜精品在线福利| 午夜久久久久精精品| 一区二区三区国产精品乱码| 精品不卡国产一区二区三区| 男人舔女人的私密视频| 国产真人三级小视频在线观看| 国产野战对白在线观看| 丰满的人妻完整版| 久久香蕉国产精品| 日本 欧美在线| 在线免费观看不下载黄p国产 | 午夜亚洲福利在线播放| 久久九九热精品免费| 国产精品 欧美亚洲| 国产成+人综合+亚洲专区| 亚洲自拍偷在线| www.999成人在线观看| 岛国在线免费视频观看| 亚洲欧美激情综合另类| 色尼玛亚洲综合影院| 一区二区三区国产精品乱码| 亚洲av五月六月丁香网| 久久国产乱子伦精品免费另类| 男女床上黄色一级片免费看| 高清毛片免费观看视频网站| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成av人片在线播放无| 日本熟妇午夜| 亚洲国产高清在线一区二区三| 国产三级中文精品| 国产欧美日韩一区二区三| 90打野战视频偷拍视频| 亚洲五月天丁香| 狂野欧美激情性xxxx| 成人精品一区二区免费| 国产亚洲av高清不卡| 国产真人三级小视频在线观看| 久久久久久久午夜电影| 黄片大片在线免费观看| 亚洲欧洲精品一区二区精品久久久| 欧美最黄视频在线播放免费| 欧美不卡视频在线免费观看| 成人国产一区最新在线观看| 深夜精品福利| 欧美xxxx黑人xx丫x性爽| 成年女人毛片免费观看观看9| 亚洲专区字幕在线| 少妇人妻一区二区三区视频| 午夜福利免费观看在线| 欧美xxxx黑人xx丫x性爽| 夜夜夜夜夜久久久久| 国产成人av激情在线播放| 久久精品亚洲精品国产色婷小说| 日本黄色视频三级网站网址| av欧美777| 欧美日韩乱码在线| 操出白浆在线播放| 男女那种视频在线观看| 亚洲色图av天堂| 欧美又色又爽又黄视频| 久久这里只有精品19| 男人和女人高潮做爰伦理| 在线播放国产精品三级| av福利片在线观看| 此物有八面人人有两片| 久久精品国产清高在天天线| 叶爱在线成人免费视频播放| 亚洲人与动物交配视频| 99精品在免费线老司机午夜| 国内精品美女久久久久久| 丰满人妻一区二区三区视频av | 国产伦精品一区二区三区四那| 人人妻人人看人人澡| 亚洲人成网站在线播放欧美日韩| 亚洲成人中文字幕在线播放| 极品教师在线免费播放| 在线免费观看不下载黄p国产 | 精品日产1卡2卡| 久久精品亚洲精品国产色婷小说| 男人舔奶头视频| 噜噜噜噜噜久久久久久91| 亚洲成人久久爱视频| xxx96com| 在线a可以看的网站| 国产又色又爽无遮挡免费看| 狠狠狠狠99中文字幕| 在线免费观看不下载黄p国产 | www国产在线视频色| 久久久久久久久中文| 国产精品综合久久久久久久免费| 黄片小视频在线播放| 午夜福利在线观看吧| www.熟女人妻精品国产| 黑人巨大精品欧美一区二区mp4| 亚洲美女视频黄频| 男女视频在线观看网站免费| 亚洲国产精品久久男人天堂| 白带黄色成豆腐渣| 国产又色又爽无遮挡免费看| 最近最新中文字幕大全免费视频| 女人被狂操c到高潮| 亚洲午夜精品一区,二区,三区| 国产三级在线视频| 国产精品自产拍在线观看55亚洲| 国产精品美女特级片免费视频播放器 | 欧美xxxx黑人xx丫x性爽| www国产在线视频色| 欧美不卡视频在线免费观看| 国产精品久久久人人做人人爽| 欧美成人一区二区免费高清观看 | 好男人电影高清在线观看| 精品国产美女av久久久久小说| 日日干狠狠操夜夜爽| 一进一出抽搐动态| 精品国产乱子伦一区二区三区| 免费av不卡在线播放| 亚洲欧美日韩无卡精品| 一个人免费在线观看电影 | 精品福利观看| 最新中文字幕久久久久 | e午夜精品久久久久久久| 亚洲午夜理论影院| 亚洲乱码一区二区免费版| 久久中文字幕一级| 精品久久蜜臀av无| 免费无遮挡裸体视频| 亚洲av美国av| 国产成人影院久久av| 久久人妻av系列| 女同久久另类99精品国产91| 午夜免费观看网址| 日韩高清综合在线| 午夜精品一区二区三区免费看| 中文亚洲av片在线观看爽| 午夜免费观看网址| 欧美成狂野欧美在线观看| 99热这里只有是精品50| 久久久久久九九精品二区国产| 成在线人永久免费视频| 一卡2卡三卡四卡精品乱码亚洲| 国产免费男女视频| 啦啦啦免费观看视频1| a级毛片在线看网站| 91麻豆av在线| 757午夜福利合集在线观看| 成人国产一区最新在线观看| 久久精品综合一区二区三区| 禁无遮挡网站| 久久久久久久久免费视频了| 亚洲电影在线观看av| 99在线人妻在线中文字幕| x7x7x7水蜜桃| 少妇熟女aⅴ在线视频| 久久久成人免费电影| 男插女下体视频免费在线播放| 国内精品久久久久精免费| 亚洲成人久久性| 久久久精品欧美日韩精品| 欧洲精品卡2卡3卡4卡5卡区| 久久精品亚洲精品国产色婷小说| 亚洲成人中文字幕在线播放| 69av精品久久久久久| 黄色丝袜av网址大全| 精品国产亚洲在线| 亚洲中文av在线| 欧美丝袜亚洲另类 | 日本三级黄在线观看| 午夜两性在线视频| 变态另类丝袜制服| 国产真实乱freesex| 亚洲午夜理论影院| 国产男靠女视频免费网站| 制服丝袜大香蕉在线| 国产成人欧美在线观看| 在线国产一区二区在线| 美女大奶头视频| 又大又爽又粗| 午夜视频精品福利| 久久久久久久午夜电影| 成年女人毛片免费观看观看9| 亚洲精品一卡2卡三卡4卡5卡| 老汉色av国产亚洲站长工具| 国产97色在线日韩免费| 熟女电影av网| 天堂√8在线中文| 国内精品久久久久久久电影| 不卡av一区二区三区| 精品国产美女av久久久久小说| 99在线视频只有这里精品首页| 久久精品国产综合久久久| 国产一区二区在线观看日韩 | 亚洲 国产 在线| 在线观看美女被高潮喷水网站 | 搞女人的毛片| 亚洲av免费在线观看| 桃色一区二区三区在线观看| 亚洲色图av天堂| 九色成人免费人妻av| 色综合欧美亚洲国产小说| av天堂中文字幕网| 麻豆国产97在线/欧美| 欧美成人免费av一区二区三区| 久久精品91蜜桃| 久久中文字幕一级| 亚洲欧洲精品一区二区精品久久久| 日韩欧美 国产精品| 日本a在线网址| 91老司机精品| 亚洲国产看品久久| 嫁个100分男人电影在线观看| 国模一区二区三区四区视频 | 精品福利观看| 国产精品日韩av在线免费观看| 亚洲欧美精品综合久久99| 成人特级黄色片久久久久久久| 国产精品爽爽va在线观看网站| 午夜福利成人在线免费观看| 夜夜看夜夜爽夜夜摸| 国产精品精品国产色婷婷| 久久午夜综合久久蜜桃| 人人妻,人人澡人人爽秒播| 久久性视频一级片| 国产亚洲精品av在线| 国产精品一及| 色av中文字幕| 精品久久久久久,| 国产精品美女特级片免费视频播放器 | 99视频精品全部免费 在线 | 91老司机精品| 日韩成人在线观看一区二区三区| 国产成人精品久久二区二区91| 狂野欧美白嫩少妇大欣赏| 熟妇人妻久久中文字幕3abv| 亚洲精品在线观看二区| 国产精品精品国产色婷婷| 精品久久久久久久久久免费视频| 丁香欧美五月| 给我免费播放毛片高清在线观看| 一区二区三区国产精品乱码| 亚洲片人在线观看| 一级作爱视频免费观看| 日韩精品中文字幕看吧| 女生性感内裤真人,穿戴方法视频| 日本在线视频免费播放| 亚洲精品在线观看二区| 老鸭窝网址在线观看| 亚洲av成人一区二区三| 女生性感内裤真人,穿戴方法视频| 两个人看的免费小视频| 日本 av在线| 老司机午夜十八禁免费视频| 国产黄片美女视频| 久久香蕉国产精品| 丰满人妻熟妇乱又伦精品不卡| 91字幕亚洲| 成人无遮挡网站| 亚洲国产精品合色在线| 国产黄a三级三级三级人| 99久久精品国产亚洲精品| 欧美日韩精品网址| 在线观看美女被高潮喷水网站 | 色综合婷婷激情| 久久中文看片网| 色综合欧美亚洲国产小说| av中文乱码字幕在线| 高清毛片免费观看视频网站| 天天一区二区日本电影三级| 国产精品九九99| 国产高清videossex| 亚洲欧美日韩无卡精品| 18禁黄网站禁片午夜丰满| 亚洲一区二区三区色噜噜| 亚洲av熟女| 99国产精品99久久久久| 深夜精品福利| 午夜成年电影在线免费观看| 国产三级在线视频| 久9热在线精品视频| 黄色日韩在线| 观看美女的网站| 天堂av国产一区二区熟女人妻| 精品国产美女av久久久久小说| 老司机午夜福利在线观看视频| 午夜福利高清视频| 久久精品91无色码中文字幕| 特级一级黄色大片| 综合色av麻豆| 99热这里只有是精品50| 亚洲国产看品久久| 免费在线观看亚洲国产| 精品国产亚洲在线| 亚洲国产精品成人综合色| 中文资源天堂在线| 亚洲国产欧美人成| 观看免费一级毛片| 久久亚洲精品不卡| 日本撒尿小便嘘嘘汇集6| 人妻夜夜爽99麻豆av| 俺也久久电影网| 国产精品国产高清国产av| 免费观看人在逋| 在线国产一区二区在线| АⅤ资源中文在线天堂| 黄色丝袜av网址大全| 巨乳人妻的诱惑在线观看| x7x7x7水蜜桃| 久久久久久九九精品二区国产| 久久九九热精品免费| 麻豆成人av在线观看| 国产成人精品无人区| 亚洲午夜理论影院| 日韩三级视频一区二区三区| 黑人欧美特级aaaaaa片| 精品国产美女av久久久久小说| 亚洲av熟女| 国产1区2区3区精品| 九色成人免费人妻av| 中文字幕久久专区| 在线免费观看的www视频| 级片在线观看| 国产伦精品一区二区三区视频9 | 亚洲av日韩精品久久久久久密| 日韩大尺度精品在线看网址| 美女黄网站色视频| 国产成人啪精品午夜网站| 啦啦啦免费观看视频1| 亚洲欧美日韩东京热| 精品乱码久久久久久99久播| av黄色大香蕉| 久久精品夜夜夜夜夜久久蜜豆| 亚洲欧美精品综合久久99| 欧美在线黄色| 特大巨黑吊av在线直播| 少妇人妻一区二区三区视频| 国产69精品久久久久777片 | 88av欧美| 欧美国产日韩亚洲一区| 搞女人的毛片| 亚洲熟妇熟女久久| 欧美黑人巨大hd| 日本黄大片高清| 可以在线观看的亚洲视频| 亚洲欧美精品综合一区二区三区| 亚洲av美国av| 久久久久久久久中文| 亚洲av成人av| 亚洲国产日韩欧美精品在线观看 | 亚洲人成网站高清观看| 久久精品影院6| 国产欧美日韩精品一区二区| 亚洲欧美日韩卡通动漫| 午夜两性在线视频| 亚洲人与动物交配视频| 丰满的人妻完整版| 真人做人爱边吃奶动态| 欧美xxxx黑人xx丫x性爽| 欧美日韩黄片免| 精品日产1卡2卡| 中文字幕人妻丝袜一区二区| 超碰成人久久| 成人av在线播放网站| 久久精品aⅴ一区二区三区四区| 欧美日本视频| 久久九九热精品免费| 亚洲自拍偷在线| 久久久久久久久免费视频了| 国产高清有码在线观看视频| 两个人的视频大全免费| 精品日产1卡2卡| 亚洲精品粉嫩美女一区| 岛国在线观看网站| 99久久国产精品久久久| 18禁国产床啪视频网站| 少妇的逼水好多| 一本久久中文字幕| 亚洲成人中文字幕在线播放| 日韩大尺度精品在线看网址| 99久久无色码亚洲精品果冻| 亚洲精品在线观看二区| 国产在线精品亚洲第一网站| 最近最新中文字幕大全免费视频| 国产乱人视频| 最新美女视频免费是黄的| 亚洲av第一区精品v没综合|