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

    岷江干旱河谷優(yōu)勢灌叢對(duì)土壤微生物群落組成的影響

    2015-03-13 03:33:44王冰冰曲來葉宋成軍宮淵波
    生態(tài)學(xué)報(bào) 2015年8期
    關(guān)鍵詞:飛虹灌叢空地

    王冰冰, 曲來葉, 宋成軍, 宮淵波

    1 中國科學(xué)院生態(tài)環(huán)境研究中心, 北京 100085 2 中國科學(xué)院大學(xué), 北京 100049 3 農(nóng)業(yè)部規(guī)劃設(shè)計(jì)研究院, 北京 100125 4 四川農(nóng)業(yè)大學(xué), 雅安 625014

    岷江干旱河谷優(yōu)勢灌叢對(duì)土壤微生物群落組成的影響

    王冰冰1,2, 曲來葉1,*, 宋成軍3, 宮淵波4

    1 中國科學(xué)院生態(tài)環(huán)境研究中心, 北京 100085 2 中國科學(xué)院大學(xué), 北京 100049 3 農(nóng)業(yè)部規(guī)劃設(shè)計(jì)研究院, 北京 100125 4 四川農(nóng)業(yè)大學(xué), 雅安 625014

    通過調(diào)查岷江干旱河谷兩河口、飛虹、撮箕和牟托4個(gè)樣地優(yōu)勢灌叢及其灌叢間空地的表土土壤物理化學(xué)性質(zhì)和微生物群落組成,探討植物灌叢群落對(duì)土壤微生物群落組成的影響。研究發(fā)現(xiàn)不同灌叢種類對(duì)土壤微生物群落組成以及土壤物理化學(xué)性質(zhì)并沒有顯著影響,而同一樣地灌叢與空地間的差異卻較為顯著。灌叢下比空地土壤中具有更高的有機(jī)質(zhì)、養(yǎng)分含量,更高的土壤含水量和更低的容重,而灌叢下相對(duì)富集的養(yǎng)分資源是造成灌叢與空地間微生物群落組成差異的主要原因。不同樣地影響微生物群落的主要因子存在一定差異,但與氮相關(guān)的因子(總氮、有效氮、碳/氮比)對(duì)土壤微生物群落著非常重要的影響,特別是對(duì)土壤微生物群落總生物量和細(xì)菌類群(革蘭氏陽性菌、革蘭氏陰性菌、細(xì)菌等)。雖然不同灌叢和空地下土壤中細(xì)菌群落都沒有顯著地變化,但真菌和菌根真菌卻明顯的在灌叢下富集。在飛虹和牟托樣地,總磷和碳/磷比與真菌類群,主要指真菌和菌根真菌,表現(xiàn)出顯著正相關(guān)性,這或許反映了真菌類群對(duì)于該區(qū)域磷循環(huán)的重要作用。研究結(jié)果揭示了灌叢植被在干旱河谷地區(qū)地下生態(tài)系統(tǒng)中的重要作用,以及氮、磷這兩種養(yǎng)分元素對(duì)土壤微生物群落的重要影響。同時(shí),未來對(duì)于干旱河谷地區(qū)植物-土壤關(guān)系的研究應(yīng)該關(guān)注真菌和菌根真菌類群的作用。

    岷江; 干旱河谷; 微生物群落組成; 磷脂脂肪酸分析; 優(yōu)勢灌叢; 灌叢間空地

    全球干旱生態(tài)系統(tǒng)占據(jù)陸地表面的41%,承載著38%的人口和豐富的生物類群[1],但惡劣而脆弱的生態(tài)環(huán)境使其更易受到全球變化的影響[1-2]。與濕潤和半濕潤地區(qū)的森林和草地生態(tài)系統(tǒng)相比,人們對(duì)干旱地區(qū)生態(tài)系統(tǒng)的了解相對(duì)缺乏。干旱生態(tài)系統(tǒng)的植被主要由旱生多年生草本和矮灌木構(gòu)成,并表現(xiàn)出明顯的斑塊狀分布[3-4]。這種斑塊狀分布的形成通常是認(rèn)為受水分、養(yǎng)分資源從灌叢間空地遷移至灌叢下造成的[5]。灌叢下富集的資源有利于植物的生長,從而反過來又進(jìn)一步促進(jìn)了灌叢斑塊的形成,即所謂的“灌叢沃島效應(yīng)”[6]。研究發(fā)現(xiàn),灌叢下富饒的養(yǎng)分和溫和的環(huán)境條件也推動(dòng)了微生物生物量、有機(jī)質(zhì)分解、土壤酶活性等“沃島”的形成[7-9]。

    土壤微生物的群落組成、生物量以及其主導(dǎo)的分解過程與土壤中碳的有效性密切相關(guān)[10-11]。有機(jī)質(zhì)的輸入主要來源于植物的凋落物、根系分泌物和根際沉降等,其數(shù)量和質(zhì)量在很大程度上受到植物種類的調(diào)控。這些不同來源的有機(jī)質(zhì),對(duì)微生物的分解過程的影響差異明顯[12-13]。同時(shí)不同植物對(duì)水分、碳氮磷等養(yǎng)分元素的利用存在差異,強(qiáng)烈影響著土壤的物理、化學(xué)和生物學(xué)性質(zhì),如土壤容重、孔隙度、pH值、氮、磷等主要養(yǎng)分元素的相對(duì)比例,以及與特定植物共生的微生物類群,從而會(huì)對(duì)土壤的微生物的群落組成和功能產(chǎn)生重要影響[14]。

    岷江干旱河谷植被稀疏,水土流失嚴(yán)重,生態(tài)環(huán)境脆弱[15]。其植被以旱生灌草植被類型為主,并呈現(xiàn)出明顯的斑塊狀分布。目前對(duì)岷江干旱河谷生態(tài)系統(tǒng)的研究大多數(shù)關(guān)注了地上植物區(qū)系和一些基本的土壤物理化學(xué)性質(zhì)[15-21],尤其土壤水分的變化方面,而對(duì)地下土壤生態(tài)過程,特別是灌叢植被影響土壤微生物群落的研究匱乏。

    鑒于岷江干旱河谷斑塊狀分布的灌叢植被對(duì)于河谷的生態(tài)恢復(fù)和生態(tài)環(huán)境的改善具有重要意義。因此本研究在干旱河谷沿岸,選取基本處于同一海拔的不同水分梯度下的四個(gè)樣地中優(yōu)勢灌叢進(jìn)行了群落調(diào)查,同時(shí)對(duì)灌叢間空地(無植被覆蓋的裸地)土壤微生物群落及土壤物理化學(xué)性質(zhì)進(jìn)行研究,通過比較該地區(qū)不同優(yōu)勢灌叢植被下和灌叢間空地的土壤微生物群落特性的差異,探討灌叢群落對(duì)土壤微生物群落組成的影響,為理解干旱河谷地區(qū)養(yǎng)分循環(huán)和植被恢復(fù)提供參考。本研究探討以下3個(gè)問題: (1) 灌叢與裸地間土壤微生物群落組成的差異;(2)不同灌叢種類間土壤微生物群落組成的差異;(3)影響土壤微生物群落組成的主要土壤理化因子。

    1 實(shí)驗(yàn)方法

    1.1 研究區(qū)域概況

    岷江干旱河谷地區(qū), 行政區(qū)劃上隸屬于四川省阿壩藏族羌族自治州的汶川、理縣、茂縣、黑水和松潘,主要分布于松潘鎮(zhèn)江關(guān)以下, 經(jīng)茂縣鳳儀鎮(zhèn)至汶川縣的岷江正河,位于東經(jīng)103°41′—103°44′,北緯31°21′—31°44′之間。該區(qū)域年平均氣溫10—11 ℃左右,降水量494.8 mm 左右,年蒸發(fā)量1332.4 mm,氣候干燥。土壤類型以山地褐土和山地棕壤土為主,土體十分干燥[22-23]。植被主要由旱中生小葉落葉具刺灌木或肉質(zhì)具刺灌木及耐旱草本(禾草為主)植物組成,通常具有生長稀疏、覆蓋度低等特征,塊狀分布的耐旱灌叢群落被嚴(yán)重沖刷的裸地分隔開,嚴(yán)重退化[24-25]。

    1.2 土樣采集

    本研究選取了岷江上游干旱河谷沿岸4塊樣地,分別位于兩河口、飛虹撮箕、和牟托的多個(gè)山坡上。樣地基本情況如表1所示。

    表1 樣地基本概況Table 1 Basic characteristics of the four sites

    LHK: 兩河口; FH: 飛虹;CJ: 撮箕;MT: 牟托

    土壤樣本的采集進(jìn)行于2012年7月,沿干旱河谷核心區(qū),在相同海拔不同水分條件的四個(gè)區(qū)域分別設(shè)置樣地。每個(gè)樣地中隨機(jī)設(shè)置3個(gè)10 m×10 m的樣方,共12個(gè)樣方。每個(gè)樣方中沿對(duì)角線分別選取5個(gè)裸地(叢間空地),采集0—15 cm的表土層;再選擇樣方中占優(yōu)勢地位的兩種植物,分別選取3個(gè)個(gè)體大小中等、長勢健康的灌叢,采集下方0—15 cm的土壤,土樣立刻放入便攜式冰盒中。土樣運(yùn)回實(shí)驗(yàn)室后,立即過2 mm篩,并挑出較大的植物根系。經(jīng)冷凍干燥和球磨儀研磨后,儲(chǔ)藏在—80 ℃超低溫冰箱中 (Themo Scientific) 用于磷脂脂肪酸(PLFA)分析。另一部分土樣風(fēng)干過 2 mm 篩用于土壤含水量 (SWC)、 pH 值和電導(dǎo)率 (EC) 的測定,過1 mm篩用于有效磷 (AP) 和有效氮 (AN) 的測定,過 100 目篩用于有機(jī)碳 (SOC) 、總氮 (TN) 、總磷 (TP) 的測定。

    1.3 實(shí)驗(yàn)方法

    土壤理化指標(biāo)測定[26]:容重 (BD) 采用環(huán)刀法;SWC采用105 ℃連續(xù)烘干24 h后計(jì)算得出;pH和EC采用酸度計(jì)和電導(dǎo)率儀分別進(jìn)行測定(土∶水=1∶2. 5);TN采用元素分析儀 (vario MAX CN, Elementar) 測定;SOC測定采用低溫外熱重鉻酸鉀氧化比色法;AP采用NaHCO3浸提,鉬銻抗比色法;TP采用HClO4-H2SO4消煮-鉬銻抗比色法,連續(xù)流動(dòng)分析儀 (AutoAnalyzer 3, Bran+Luebbe) 測定; AN采用堿解擴(kuò)散法測定。

    PLFA的提取和純化參考Bligh和Dyer法[27-28],上機(jī)前用1 ml含內(nèi)標(biāo)物19:0的正己烷溶解吹干的脂肪酸甲酯,然后進(jìn)行GC-FID (Agilent 6850N) 檢測,脂肪酸種類通過脂肪酸Sherlock 微生物鑒定系統(tǒng) (MIDI Inc., Newark, DE, USA) 來鑒定。PLFAs總量 (Total PLFA, PLFAT) 及單體脂肪酸的含量基于19:0內(nèi)標(biāo)計(jì)算,分析中僅包括碳鏈數(shù)低于20、摩爾百分比 (mol%) 含量大于1% 且在多數(shù)樣品中都有出現(xiàn)的脂肪酸。細(xì)菌 (bacteria, B) 以14:0、i14:0、i15:0、a15:0、15:1ω6c、i16:0、16:1ω7c、16:1ω9c、i17:0、a17:0、17:1ω8c、16:1 2OH、cy17:0、18:1ω5c、cy17:0和cy19:0 C11—12等15種脂肪酸表示[28];真菌 (F) 用18:2ω6c脂肪酸來估計(jì)[29];真菌:細(xì)菌比通過F/B×100來計(jì)算;AM真菌 (AMF) 以16:1ω5c表示[30];革蘭氏陽性細(xì)菌(G+) 用支鏈、飽和脂肪酸i14:0、i15:0、a15:0、i16:0、i17:0和a17:0表示;革蘭氏陰性細(xì)菌(G-) 用單不飽和脂肪酸和環(huán)式脂肪酸15:1ω6c、16:1ω7c、17:1ω8c、18:1ω5c、16:1 2OH、cy17:0、cy19:0 C11—12表示[28, 31];革蘭氏陽性與陰性細(xì)菌比通過G+/G-計(jì)算;環(huán)丙烷脂肪酸 (cyclo) 以cy17:0和cy19:0 C11—12表示[32],環(huán)丙烷脂肪酸與其前體脂肪酸比值 (cy/pre) 通過cy17:0/16:1ω7c來實(shí)現(xiàn)[32-33]。

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

    采用參數(shù)檢驗(yàn)單因素方差分析 (One-way ANOVA) 和非參數(shù)方差分析(Kruskal-Wallis One-way ANOVA),比較不同植物種類之間和植物與空地之間,土壤物理化學(xué)性質(zhì)和土壤微生物主要功能類型(PLFAT、F、B、G+、G-等)生物量的差異。參數(shù)檢驗(yàn)的多重比較采用Fisher LSD法,Kruskal-Wallis檢驗(yàn)的多重比較方法參見R軟件統(tǒng)計(jì)包pgirmess的kruskalmc函數(shù)[34]。所有數(shù)據(jù)在分析前都進(jìn)行了Shapiro-Wilk正態(tài)分布和levene方差齊質(zhì)性檢驗(yàn)。

    為從整體上探討灌叢與空地以及不同灌叢種類間土壤理化性質(zhì)和微生物群落組成的差異,對(duì)理化數(shù)據(jù)和PLFA數(shù)據(jù)分別做了基于相關(guān)的主成分分析 (PCA),并對(duì)PCA的前3個(gè)主成分軸進(jìn)行了單因素方差分析。PLFA數(shù)據(jù)在分析前進(jìn)行了Chord距離轉(zhuǎn)化[35]。利用逐步回歸后向消去法 (Backward elimination) 和冗余分析 (RDA) 來探討土壤微生物主要類群與土壤理化性質(zhì)的相關(guān)性。PCA、逐步回歸和RDA都使用R軟件包vegan[36]。

    以上所有分析都在R 3.0.0軟件中進(jìn)行,繪圖使用OriginPro8.6和R 3.0.0軟件。

    2 結(jié)果與分析

    2.1 土壤理化性質(zhì)變異

    四個(gè)樣地中,灌叢與空地之間、不同灌叢之間土壤物理化學(xué)性質(zhì)的差異表現(xiàn)不同(表2)。兩河口和牟托的多數(shù)理化指標(biāo)在灌叢與空地之間、不同灌叢之間存在顯著差異,而撮箕僅有含水率和容重在金花噪草和空地間有明顯差異,飛虹所有因子都沒有顯著差異??傮w來說,灌叢下土壤較空地土壤具有較高的SWC、SOC、TN和AN,較低的土壤pH和BD。TP、C/P比、AP在4個(gè)樣地中均沒有顯著差異,其變幅也很??;C/N比只有在兩河口空地顯著的高于灌叢下。EC在飛虹和兩河口都具有較大的變幅,如飛虹空地的EC比川甘亞菊下高了近30%,但由于組內(nèi)方差較大,并沒有表現(xiàn)出統(tǒng)計(jì)學(xué)上的差異顯著性。

    灌叢與空地土壤理化性質(zhì)在PCA軸上得到較好的分離(圖1),其中兩河口小馬鞍羊蹄甲與空地間土壤理化性質(zhì)在PC2軸上差異顯著(P<0.05);撮箕山金花噪草與空地在PC3軸上差異顯著;牟托光果疣、鐵桿蒿與空地在PC1軸上差異極顯著(P<0.01)。飛虹鐵桿蒿、川甘亞菊兩種植物與空地在PC1、PC2軸上盡管都有較好的分離,但均未達(dá)到顯著差異。

    表2 不同優(yōu)勢灌叢、空間間土壤物理化學(xué)性質(zhì)方差分析Table 2 One-way ANOVA of soil physiochemical properties under dominant shrubs and shrub-interspace

    K、T、B、Y、C、J和G分別表示灌叢間空地、鐵桿蒿、白刺花、小馬鞍羊蹄甲、川甘亞菊、金花噪草和光果蕕,后面表示方法與此一樣;表中數(shù)值為均值±標(biāo)準(zhǔn)誤,多重比較結(jié)果以字母a/b (差異顯著,P<0.05)和A/B (差異極顯著,P<0.01) 表示;斜體數(shù)值為采用了非參數(shù)檢驗(yàn);SWC: 土壤含水量 soil water content;BD: 容重 bulk density;EC: 電導(dǎo)率 electronic conductivity;SOC: 土壤有機(jī)碳 soil organic carbon;TN: 總氮 total nitrogen;TP: 總磷 total phosphorus;C/N: 碳氮比 carbon to nitrogen ratio;C/P: 碳磷比 carbon to phosphorus ratio;AP: 有效磷 available phosphorus;AN: 有效氮 available nitrogen

    圖1 4個(gè)樣地土壤物理化學(xué)性質(zhì)主成分分析(PCA)Fig.1 Principal component analysis of Soil physiochemical properties at the four sites橫豎誤差棒分別為PC1軸和PC2/PC3軸的標(biāo)準(zhǔn)誤;K、T、B、Y、C、J和G分別表示灌叢間空地、鐵桿蒿、白刺花、小馬鞍羊蹄甲、川甘亞菊、金花噪草和光果蕕

    2.2 土壤微生物群落組成變異

    對(duì)不同類群微生物磷脂脂肪酸含量(mol%)的方差分析 (表3) 可知,撮箕各微生物類群在不同灌叢、灌叢與空地間都沒有顯著性差異。兩河口、飛虹和牟托3個(gè)樣地中,灌叢下土壤微生物PLFAT、F、AMF和F/B 都顯著地高于空地。G+僅在兩河口表現(xiàn)出差異顯著性,G-及G+/G-都僅在牟托樣地有顯著地變化,而B在4個(gè)樣地均無明顯變化。兩河口樣地,空地cyclo和cy/pre比都極顯著的高于灌叢下,其他樣地則無顯著變化。不同灌叢種類間,微生物類群基本無顯著差異。

    不同樣地土壤微生物磷脂脂肪酸(PLFAs)主成分分析結(jié)果如圖2所示,PCA前三軸可較好反映PLFAs的變異。兩河口樣地,PCA前三軸解釋了67%的變異,其中PC1軸可解釋36.37%的PLFA變異,并與i15:1G、16:1ω5c、a17:0、16:1 2OH、18:2ω6c、18:1ω9c、18:1ω5c等脂肪酸表現(xiàn)出很強(qiáng)的相關(guān)性(圖2)。灌叢與空地在PC1軸上有著極顯著的差異,3種灌叢明顯的聚為一類。灌叢中18:2ω6c、16:1ω5c、18:1ω5c、i15:1G等脂肪酸比例較高,而空地中則富集a17:0、16:1 2OH、cy19:0、18:1ω9c、 i17:0等脂肪酸。

    飛虹、撮箕和牟托樣地,PCA前三軸共分別解釋了69%、67%和62%的變異。PC1軸上這3個(gè)樣地的空地與灌叢、灌叢與灌叢之間都沒有顯著差異。相比之下,PC2軸較PC1軸更能反映微生物組成的變異(P<0.05)。飛虹樣地,PC2軸可以解釋15.06%的變異,鐵桿蒿與空地在PC2軸上表現(xiàn)出顯著地差異性,鐵桿蒿下18:2ω6c、16:1ω5c、i16:1G、16:0等脂肪酸比例較高,空地中則相對(duì)富集更多的18:1ω9c、16:1 2OH、i17:1G(圖2)。撮箕樣地,PC2軸可解釋17.00%的變異,金花噪草與空地在PC2軸上差異顯著,金花噪草下i14:0、18:0、18:1ω9c、16:1 2OH、cy19:0相對(duì)富集,空地中則具有更多的15:0、16:1ω9c、18:1ω5c等脂肪酸(圖2)。牟托樣地, PC2軸可以解釋22.71%的變異,鐵桿蒿與空地在PC2軸上差異顯著,鐵桿蒿下i17:1G、16:1ω9c、i15:1G、18:1ω5c等脂肪酸比例相對(duì)較高,而空地中i17:0、18:0、16:1 2OH、cy19:0等脂肪酸比例較高(圖2)。

    表3 不同優(yōu)勢灌叢、空地間微生物類群方差分析

    Table 3 One-way ANOVA of PLFAs concentration of different soil microbial groups under dominant shrubs and in the shrub-interspaces

    樣地Sites因子Factors總磷脂脂肪酸濃度PLFAT/(nmol/g)革蘭氏陽性菌G+/%革蘭氏陰性菌G-/%革蘭氏陽性:陰性細(xì)菌比G+/G-細(xì)菌B/%真菌F/%真菌:細(xì)菌比F/B/%菌根真菌AMF/%環(huán)丙烷脂肪酸cyclo/%環(huán)丙烷脂肪酸與其前體脂肪酸比cy/preLHKK52.61±4.37aA25.53±0.67b24.08±0.561.07±0.0450.70±0.672.59±0.17aA5.11±0.33aA2.72±0.11aA6.90±0.20bB0.41±0.01bBT82.59±5.944bB23.31±0.69a23.40±0.511.00±0.0347.90±0.595.01±0.63bB10.50±1.33bB4.04±0.39bB5.83±0.21aB0.31±0.02aAY89.60±5.21bB23.42±0.36a23.85±0.550.98±0.0348.43±1.114.01±0.40bB8.28±0.81bAB4.01±0.49bB5.83±0.25aB0.32±0.02aAB75.84±5.38bB23.41±0.67a23.72±0.660.99±0.0248.30±1.234.67±0.43bB9.65±0.82bB4.15±0.35bB5.43±0.15aA0.28±0.02aAFHK64.38±5.64a24.34±0.6026.97±0.790.92±0.0452.54±0.722.99±0.23a5.70±0.45a3.26±0.30a6.57±0.330.30±0.02T83.62±7.72b25.07±0.9125.96±1.210.97±0.0451.18±1.694.45±0.14b7.73±0.50b4.27±0.28b5.98±0.270.26±0.01C80.20±5.42ab25.14±0.3924.96±1.021.04±0.0451.89±1.323.32±0.36a6.47±0.62ab3.63±0.18ab6.07±0.200.28±0.01CJK75.88±7.5324.27±0.6026.86±0.780.91±0.0252.30±1.274.60±0.379.07±0.994.25±0.167.12±0.290.22±0.00T68.62±4.4923.94±0.7226.87±1.250.90±0.0351.92±1.924.87±0.469.29±0.654.71±0.277.20±0.480.21±0.01J82.74±5.9925.39±1.0029.13±1.410.87±0.0155.59±2.425.85±0.4510.56±0.774.91±0.377.09±0.370.21±0.02MTK53.93±5.09a26.29±0.6223.87±0.26b1.10±0.02ab51.22±0.802.87±0.14a5.63±0.30a3.45±0.11a6.75±0.210.29±0.01T88.07±4.36b25.23±0.6224.57±0.55b1.03±0.03a50.90±1.003.71±0.32b7.23±0.50b4.39±0.32b6.18±0.300.27±0.02G72.89±10.55b26.98±0.9022.05±1.11a1.25±0.11b50.12±0.923.70±0.36b7.39±0.73b4.14±0.33ab6.51±0.230.30±0.02

    PLFAT 總磷脂脂肪酸濃度 total phospholipid fatty acid; G+: 革蘭氏陽性菌 Gran-positive bacteria; G-: 革蘭氏陰性菌 Gran-negative bacteria; B: 細(xì)菌 general bacteria; F: 真菌 fungi; F/B: 真菌細(xì)菌比 fungi to bacteria ratio; AMF: 菌根真菌 arbuscular mycorrhizal fungi; cyclo: 環(huán)丙烷脂肪酸 cyclopropyl; cy/pre: 環(huán)丙烷脂肪酸與其前體脂肪酸比cyclo to precursor ratio

    圖2 4個(gè)樣地土壤微生物磷脂脂肪酸(PLFAs)主成分分析Fig.2 Principal component analysis of soil microbial PLFAs at the four sites橫豎誤差棒分別為PC1軸和PC2軸的標(biāo)準(zhǔn)誤

    2.3 土壤因子與微生物群落組成的相關(guān)性

    RDA分析表明,獲得的所有土壤理化指標(biāo)共同可解釋兩河口樣地微生物類群變異的38.32%,飛虹樣地的55.99%,撮箕樣地的50.59%,牟托樣地的61.93%。兩河口樣地,C/N和AN是影響微生物類群變異的兩個(gè)主要因子,可以解釋27.44%的變異,占全部因子解釋的71.61% (圖3)。對(duì)于飛虹,TP和TN是主導(dǎo)的因子,解釋了39.66%的微生物類群變異量,占全部因子解釋的70.83% (圖3)。在撮箕山,C/N、AN和SWC 3個(gè)因子共同可解釋35.98%的變異,占全部可解釋變異的71.12% (圖3)。在牟托,AN、C/P、SWC和pH值可解釋48.94%的土壤微生物類群變異,為全部因子解釋量的79.02% (圖3)。

    4個(gè)樣地中,PLFAT與RDA1軸表現(xiàn)出明顯的正相關(guān)性(圖3,表4),其與細(xì)菌、真菌類群在RDA排序圖上表現(xiàn)出明顯的不同。F、AMF和F/B在RDA圖上的分布較為相似,而細(xì)菌類群,包括B、G+、G-、cyclo以及G+/G-和cy/pre也較為明顯的聚集在一起。

    3 討論

    3.1 灌叢與空地間的土壤理化性質(zhì)和微生物群落組成差異

    由于根系和凋落物的輸入,灌叢下土壤往往較空地中具有更高的C、N、P等養(yǎng)分含量[37-38]。灌叢可影響土壤的形成過程,緩解降雨和冰雪消融對(duì)地表的沖刷,減少土壤養(yǎng)分的淋失,從而有效的改善土壤的物理和化學(xué)性質(zhì)[39]。本研究中,牟托和兩河口樣地,灌叢下土壤較空地土壤具有更高的有機(jī)碳和養(yǎng)分含量,更為適宜的土壤環(huán)境(更高的土壤含水量和更低的容重),更為豐富的PLFAT、F、AMF含量及F/B(表2,表3)。這些說明灌叢所創(chuàng)造的良好微環(huán)境總體上促進(jìn)了微生物生物量的增加及真菌類群的繁衍。此外,兩河口灌叢下土壤中的cyclo和cy/pre都極顯著低于空地中的。飛虹樣地的盡管未達(dá)到顯著差異,但灌叢下這兩者也都低于空地中。研究表明,cyclo和cy/pre可以指示環(huán)境的脅迫,高溫、饑餓狀態(tài)下,這兩者的值都會(huì)增高,即環(huán)丙烷脂肪酸含量增加而其前體脂肪酸含量下降[32, 40]。兩河口和飛虹樣地是本研究中最為干旱的兩個(gè)樣地,其植被覆蓋度也相對(duì)較低,裸地直接遭受陽光的曝曬,加上有機(jī)質(zhì)輸入?yún)T乏,造成空地相對(duì)不適合微生物的生存。這說明,灌叢的存在有利于降低環(huán)境對(duì)微生物生存的脅迫,而這也從側(cè)面反映了灌叢植被對(duì)于干旱河谷樣地地下生態(tài)系統(tǒng)的重要意義。

    圖3 不同樣地土壤微生物群落主要類群與土壤理化因子的冗余分析(RDA)Fig.3 Redundancy Analysis of the main types of microbial community and soil physiochemical factors at different sitesSWC: 土壤含水量 soil water content; TN: 總氮 total nitrogen; TP: 總磷 total phosphorus; C/N: 碳氮比 carbon to nitrogen ratio; C/P: 碳磷比 carbon to phosphorus ratio; AN: 有效氮 available nitrogen

    表4 通過逐步回歸選擇出來的主要環(huán)境因子與土壤微生物主要類群的Pearson相關(guān)分析Table 4 Pearson correlation between the main types of soil microbial and the environmental variables selected by stepwise model selection

    *P<0.05, **P<0.01

    撮箕和飛虹兩地,灌叢與其灌叢間空地上的土壤理化性質(zhì)和微生物群落組成差異并不總是顯著的(表2,表3)。撮箕灌叢與空地土壤性質(zhì)的差異僅僅表現(xiàn)在pH和BD兩個(gè)理化指標(biāo)上,而飛虹的也主要是體現(xiàn)在PLFAT和真菌類指標(biāo)上。對(duì)理化性質(zhì)的PCA排序顯示,僅僅撮箕樣地金花噪草與空地在PC3軸表現(xiàn)出顯著差異。而對(duì)PLFA數(shù)據(jù)的PCA分析卻顯示,兩個(gè)樣地微生物群落組成在PC2軸上都有顯著差異。這或許能夠說明土壤微生物是一個(gè)相對(duì)更為敏感的指標(biāo),更能反映植物對(duì)土壤性質(zhì)的影響。Zornoza 等也發(fā)現(xiàn),微生物群落組成可以比生物量等指標(biāo)更為快速的反映土地利用的變化,從而可以作為評(píng)價(jià)土壤肥力、擾動(dòng)和管理措施效應(yīng)的有效指標(biāo)[41]。

    傳統(tǒng)上認(rèn)為“灌叢沃島”與空地是兩個(gè)生物地球化學(xué)循環(huán)完全不同的區(qū)域[42-43],然而現(xiàn)在越來越多的研究表明,這兩者之間存在著千絲萬縷的聯(lián)系[44]。這種聯(lián)系的發(fā)生,很大程度上要?dú)w功于龐大的地下“真菌網(wǎng)絡(luò)”。通過“真菌網(wǎng)絡(luò)”,一個(gè)個(gè)獨(dú)立的、斑塊狀分布的“灌叢沃島”,以及灌叢之間分隔的空地都有機(jī)的聯(lián)系起來。本研究中撮箕和飛虹兩地灌叢與空地間或許通過“真菌網(wǎng)絡(luò)”有機(jī)的聯(lián)系起來,從而造成兩者之間差異不甚明顯。

    而另一方面,這種較小的差異,或許指示著灌叢植被對(duì)其周圍空地土壤具有一定的改善作用。灌從的存在有利凋落物的聚集,從而加速養(yǎng)分循環(huán),促進(jìn)植物根系的生長,導(dǎo)致更多的根際沉積,灌叢下土壤微生物代謝活性增強(qiáng),而根際的活躍又會(huì)進(jìn)一步增加灌叢下土壤中的資源量。在自然恢復(fù)的過程中,利用“土壤-植被-土壤微生物”系統(tǒng)之間的相互促進(jìn)作用改善局部環(huán)境,通過植被的不斷恢復(fù),可以提高相鄰空地的土壤質(zhì)量,有利于植被更新和生長,也有可能形成良性循環(huán),加快區(qū)域內(nèi)植被恢復(fù)進(jìn)程。李新榮報(bào)道了沙坡頭地區(qū)人工固沙植被近50年的演變與區(qū)域生境的恢復(fù),證明了沙漠綴塊狀分布的灌木植被在植被恢復(fù)過程中的重要作用[45]。岷江干旱河谷地區(qū)相比于沙坡頭地區(qū)具有更高的降雨量,在沒有人為干擾(破壞)的情況下,干旱河谷地區(qū)目前的斑塊狀分布的灌叢群落,有可能成為該區(qū)域植被恢復(fù)的基礎(chǔ)。

    3.2 不同灌叢間土壤理化性質(zhì)和微生物群落組成差異

    已有研究表明土壤理化性質(zhì)和微生物群落組成在不同植物下存在一定差異[46-49]。但本研究發(fā)現(xiàn)灌叢間(主要是各個(gè)樣地中優(yōu)勢種鐵桿蒿和其他種灌叢間)的土壤理化性質(zhì)和微生物群落組成都沒有表現(xiàn)出顯著的差異 (表2,表3)。所調(diào)查的4個(gè)樣地中,僅兩河口的小馬羊蹄甲灌叢下土壤pH、SOC和C/P比顯著高于白刺花灌叢下的;真菌生物量僅在鐵桿蒿和川甘亞菊間表現(xiàn)出顯著差異。造成這種結(jié)果的原因,或許存在一定的偶然因素,但不同種類植物在這種干旱條件下的趨同適應(yīng)可能是一個(gè)更為根本的原因[50]。趨同適應(yīng)的結(jié)果使得不同灌叢植物對(duì)土壤物理化學(xué)、生物學(xué)上的影響差異相對(duì)縮小。

    3.3 土壤理化因子對(duì)微生物群落組成的影響

    植物對(duì)微生物群落組成的影響常??赏ㄟ^土壤的物理化學(xué)性質(zhì)來調(diào)節(jié)。土壤pH、SOC被認(rèn)為是影響微生物群落結(jié)構(gòu)的主要因子[10, 49, 51-53],但是養(yǎng)分元素的缺失、有機(jī)質(zhì)的質(zhì)量 (C/N、C/P) 以及SWC都會(huì)限制微生物的生長[10, 49, 53-54]。

    從RDA排序結(jié)果來看,N是影響岷江干旱河谷地區(qū)土壤微生物群落組成的一個(gè)非常重要的因子(圖3,表4)。4個(gè)樣地中,與N相關(guān)的因子 (C/N、TN、AN) 都對(duì)微生物群落產(chǎn)生重要影響,特別是PLFAT和細(xì)菌群落。兩河口樣地的AN和C/N比、飛虹的TN、撮箕的C/N比以及牟托的AN都主要與PLFAT、細(xì)菌群落(B、G+、G-和cyclo)、G+/G-和cy/pre有顯著地相關(guān)性。而另一方面,PLFAT由細(xì)菌貢獻(xiàn)了很大一部分,從而N或許是影響該地區(qū)細(xì)菌群落變異的一個(gè)重要因素。黃元元等對(duì)該地區(qū)灌叢“肥島”效應(yīng)的研究也發(fā)現(xiàn),AN、TN和C/N對(duì)細(xì)菌群落有著顯著地影響[55]。

    對(duì)于飛虹和牟托樣地的微生物群落組成來說,P是另一個(gè)重要的影響因素。飛虹樣地,TP與AMF和F表現(xiàn)出顯著地正相關(guān)性;牟托樣地,隨著C/P比的增加,F(xiàn)、AMF比例減少,F(xiàn)/B比值降低。從而說明,土壤中P的增加或者有機(jī)質(zhì)中P含量的提高有利于F和AMF的增加,反過來則反映了F和AMF對(duì)該區(qū)域P循環(huán)的重要價(jià)值。Burke等、黃元元在研究中也發(fā)現(xiàn)土壤P對(duì)真菌群落的分布的構(gòu)建作用[56]。干旱河谷地區(qū),雖然土壤中磷的固定較多,但土壤有效磷尤其偏低,土壤潛在性肥力向有效性肥力轉(zhuǎn)化的能力較差,致使土壤有效性肥力水平較低[23],從而植物增加對(duì)AMF的投資,有利于自身對(duì)P等限制性養(yǎng)分元素的吸收利用[57]。

    對(duì)于撮箕和牟托樣地,SWC對(duì)土壤微生物類群也有著重要影響,其單獨(dú)分別可以解釋兩樣地17.85%和23.35%的微生物類群變異。不同的是,撮箕樣地,SWC與PLFAT呈顯著的負(fù)相關(guān)關(guān)系,與細(xì)菌類群表現(xiàn)出顯著的正相關(guān)關(guān)系,而牟托樣地則正好相反。此外,兩河口與撮箕的C/N比同微生物類群的關(guān)系也表現(xiàn)出類似的變化。事實(shí)上微生物群落組成受到很多因素的影響,某些本研究未涉及的因素如不同樣地溫度、團(tuán)聚體、土壤有機(jī)酸等的差異可能比SWC、C/N更加強(qiáng)烈的影響著微生物群落組成,從而造成這種差異[54, 58-59]。

    從RDA排序中(圖3,表4)可見,這些具有顯著影響的因子,對(duì)F、F/B和AMF類群的變異解釋的相對(duì)較少。真菌類群,特別是AMF,與植物之間有著很密切的聯(lián)系。本研究中,兩河口、飛虹和牟托3個(gè)樣地中,灌叢、空地間真菌類群有著顯著地差異,因此,其可能更多的受到來自植物的直接調(diào)控。

    pH值通常被認(rèn)為是影響微生物群落的一個(gè)重要因素,但本研究中pH的作用并未突出,4個(gè)樣地中,pH對(duì)微生物類群變異的作用較小,這大概與本研究中pH的較小變異有關(guān)。Pietri和Brookes也發(fā)現(xiàn),干旱土壤中pH對(duì)微生物的作用相對(duì)于濕潤的森林土壤來說并不明顯[60]。事實(shí)上,本研究中,pH也只在濕度較高的牟托樣地表現(xiàn)出一定的影響。

    總體來說,灌叢下相對(duì)更高的養(yǎng)分含量 (AN、TN、TP)、更優(yōu)的物理環(huán)境 (SWC、BD),是造成灌叢下微生物,特別是真菌、AMF類群富集的重要因素。

    4 結(jié)論

    綜上所述,灌叢與空地間土壤理化性質(zhì)和微生物群落組成存在顯著差異,而不同灌叢種類間則無明顯差異。灌叢下土壤相比于空地,具有更高的養(yǎng)分含量、更為適宜的土壤物理?xiàng)l件,以及更為豐富的真菌和菌根真菌群落。N是影響干旱河谷地區(qū)土壤微生物群落的主要因素,其次是P和土壤含水量。此外,P與真菌、AMF類群表現(xiàn)的很強(qiáng)相關(guān)性,表明真菌和AMF對(duì)于該地區(qū)土壤中P的有效性和P循環(huán)的重要作用。

    [1] Maestre F, Quero J, Gotelli N, Escudero A, Ochoa V, Delgado-Baquerizo M, García-Gómez M, Bowker M, Soliveres S, Escolar C, García-Palacios P, Berdugo M, Valencia E, Gozalo B, Gallardo A, Aguilera L, Arredondo T, Blones J, Boeken B, Bran D, Concei??o A, Cabrera O, Chaieb M, Derak M, Eldridge D, Espinosa C, Florentino A, Gaitán J, Gatica M, Ghiloufi W, Gómez-González S, Gutiérrez J, Hernández R, Huang X, Huber-Sannwald E, Jankju M, Miriti M, Monerris J, Mau R, Morici E, Naseri K, Ospina A, Polo V, Prina A, Pucheta E, Ramírez-Collantes D, Rom?o R, Tighe M, Torres-Díaz C, Val J, Veiga J, Wang D, Zaady E. Plant species richness and ecosystem multifunctionality in global drylands. Science, 2012, 335(6065): 214-218.

    [2] Seager R, Ting M F, Held I, Kushnir Y, Lu J, Vecchi G, Huang H P, Harnik N, Leetmaa A, Lau N C, Li C H, Velez J, Naik N. Model projections of an imminent transition to a more arid climate in southwestern North America. Science, 2007, 316(5828): 1181-1184.

    [3] Noy-Meir I. Desert ecosystems: environment and producers. Annual Review of Ecology and Systematics, 1973, 4(1): 25-51.

    [4] Kéfi S, Rietkerk M, Alados C L, Pueyo Y, Papanastasis V P, ElAich A, De Ruiter P C. Spatial vegetation patterns and imminent desertification in Mediterranean arid ecosystems. Nature, 2007, 449(7159): 213-217.

    [5] Schlesinger W H, Raikes J A, Hartley A E, Cross A F. On the spatial pattern of soil nutrients in desert ecosystems. Ecology, 1996, 77(2): 364-374.

    [6] Bolling J D, Walker L R. Fertile island development around perennial shrubs across a Mojave Desert chronosequence. Western North American Naturalist, 2002, 62(1): 88-100.

    [7] Ben-David E A, Zaady E, Sher Y, Nejidat A. Assessment of the spatial distribution of soil microbial communities in patchy arid and semi-arid landscapes of the Negev Desert using combined PLFA and DGGE analyses. FEMS Microbiology Ecology, 2011, 76(3): 492-503.

    [8] Goberna M, Pascual J A, García C, Sánchez J. Do plant clumps constitute microbial hotspots in semiarid Mediterranean patchy landscapes? Soil Biology and Biochemistry, 2007, 39(5): 1047-1054.

    [9] Thompson T L, Zaady E, Huancheng P, Wilson T B, Martens D A. Soil C and N pools in patchy shrublands of the Negev and Chihuahuan Deserts. Soil Biology and Biochemistry, 2006, 38(7): 1943-1955.

    [10] Meril? P, Malmivaara-L?ms? M, Spetz P, Stark S, Vierikko K, Derome J, Fritze H. Soil organic matter quality as a link between microbial community structure and vegetation composition along a successional gradient in a boreal forest. Applied Soil Ecology, 2010, 46(2): 259-267.

    [11] Weintraub S P, Wieder W R, Cleveland C C, Townsend A R. Organic matter inputs shift soil enzyme activity and allocation patterns in a wet tropical forest. Biogeochemistry, 2013, 114(1/3): 313-326.

    [12] Kuzyakov Y, Blagodatskaya E, Blagodatsky S. Comments on the paper by Kemmitt et al. (2008) ‘Mineralization of native soil organic matter is not regulated by the size, activity or composition of the soil microbial biomass-A new perspective′ [Soil Biology & Biochemistry 40, 61-73]: The biology of the Regulatory Gate. Soil Biology and Biochemistry, 2009, 41(2): 435-439.

    [13] Kaiser C, Koranda M, Kitzler B, Fuchslueger L, Schnecker J, Schweiger P, Rasche F, Zechmeister-Boltenstern S, Sessitsch A, Richter A. Belowground carbon allocation by trees drives seasonal patterns of extracellular enzyme activities by altering microbial community composition in a beech forest soil. New Phytologist, 2010, 187(3): 843-858.

    [14] Ayres E, Steltzer H, Berg S, Wallenstein M D, Simmons B L, Wall D H. Tree species traits influence soil physical, chemical, and biological properties in high elevation forests. PLoS One, 2009, 4(6): e5964.

    [15] Ma K M, Fu B J, Liu S L, Guan W B, Liu G H, Lu Y H, Anand M. Multiple-scale soil moisture distribution and its implications for ecosystem restoration in an Arid River valley, China. Land Degradation and Development, 2004, 15(1): 75-85.

    [16] 劉慶, 包維楷, 喬永康, 錢能斌. 岷江上游茂縣半干旱河谷灌叢優(yōu)勢種間關(guān)系的研究. 應(yīng)用與環(huán)境生物學(xué)報(bào), 1996, 2(1): 36-42.

    [17] 楊欽周. 岷江上游干旱河谷灌叢研究. 山地學(xué)報(bào), 2007, 25(1): 1-32.

    [18] 劉文彬. 岷江上游半干旱河谷灌叢植物區(qū)系. 山地研究, 1992, 10(2): 83-88.

    [19] 冶民生, 關(guān)文彬, 吳斌, 馬克明, 劉國華, 汪西林. 岷江干旱河谷主要灌木種群生態(tài)位研究. 北京林業(yè)大學(xué)學(xué)報(bào), 2006, 28(1): 7-13.

    [20] 劉國華, 馬克明, 傅伯杰, 關(guān)文彬, 康永祥, 周建云, 劉世梁. 岷江干旱河谷主要灌叢類型地上生物量研究. 生態(tài)學(xué)報(bào), 2003, 23(9): 1757-1764.

    [21] 陳廣生, 曾德慧, 陳伏生, 范志平, 耿海東. 干旱和半干旱地區(qū)灌木下土壤 “肥島” 研究進(jìn)展. 應(yīng)用生態(tài)學(xué)報(bào), 2003, 14(12): 2295-2300.

    [22] 冶民生, 關(guān)文彬, 譚輝, 馬克明, 劉國華, 汪西林. 岷江干旱河谷灌叢α多樣性分析. 生態(tài)學(xué)報(bào), 2004, 24(6): 1123-1130.

    [23] 王春明, 包維楷, 陳建中, 孫輝, 謝嘉穗. 岷江上游干旱河谷區(qū)褐土不同亞類剖面及養(yǎng)分特征. 應(yīng)用與環(huán)境生物學(xué)報(bào), 2003, 9(3): 230-234.

    [24] 陳泓. 岷江上游干旱河谷灌叢群落物種多樣性與生產(chǎn)力的海拔梯度格局研究 [D]. 雅安: 四川農(nóng)業(yè)大學(xué), 2007.

    [25] 黎燕瓊, 劉興良, 鄭紹偉, 陳泓, 岳永杰, 慕長龍, 劉軍. 岷江上游干旱河谷四種灌木的抗旱生理動(dòng)態(tài)變化. 生態(tài)學(xué)報(bào), 2007, 27(3): 870-878.

    [26] 鮑士旦. 土壤農(nóng)化分析 (第三版). 北京: 中國農(nóng)業(yè)出版社, 2000.

    [27] Bligh E G, Dyer W J. A rapid method of total lipid extraction and purification. Canadian Journal of Biochemistry and Physiology, 1959, 37(8): 911-917.

    [28] Waldrop M P, Balser T C, Firestone M K. Linking microbial community composition to function in a tropical soil. Soil Biology and Biochemistry, 2000, 32(13): 1837-1846.

    [29] Frosteg?rd A, Tunlid A, B??th E. Phospholipid Fatty Acid composition, biomass, and activity of microbial communities from two soil types experimentally exposed to different heavy metals. Applied and Environmental Microbiology, 1993, 59(11): 3605-3617.

    [30] Kourtev P S, Ehrenfeld J G, Haggblom M. Exotic plant species alter the microbial community structure and function in the soil. Ecology, 2002, 83(11): 3152-3166.

    [31] Frosteg?rd A, B??th E. The use of phospholipid fatty acid analysis to estimate bacterial and fungal biomass in soil. Biology and Fertility of Soils, 1996, 22(1/2): 59-65.

    [32] Bossio D A, Scow K M. Impacts of carbon and flooding on soil microbial communities: Phospholipid fatty acid profiles and substrate utilization patterns. Microbial Ecology, 1998, 35(3): 265-278.

    [33] Chen D M, Zhou L X, Wu J P, Hsu J N, Lin Y B, Fu S L. Tree girdling affects the soil microbial community by modifying resource availability in two subtropical plantations. Applied Soil Ecology, 2012, 53: 108-115.

    [34] Giraudoux P. pgirmess: Data analysis in ecology. R package version 1. 5. 7. 2013.

    [35] Legendre P, Gallagher E D. Ecologically meaningful transformations for ordination of species data. Oecologia, 2001, 129(2): 271-280.

    [36] Oksanen J, Blanchet F G, Kindt R, Legendre P, Minchin P R, O′Hara R B, Simpson G L, Solymos P, Stevens M H H, Wagner H. vegan: Community Ecology Package. R package version 2. 0-7. 2013.

    [37] Gallardo A, Schlesinger W H. Carbon and Nitrogen Limitations of Soil Microbial Biomass in Desert Ecosystems. Biogeochemistry, 1992, 18(1): 1-17.

    [38] Schlesinger W H, Pilmanis A M. Plant-soil interactions in deserts. Biogeochemistry, 1998, 42(1/2): 169-187.

    [39] Aanderud Z T, Shuldman M I, Drenovsky R E, Richards J H. Shrub-interspace dynamics alter relationships between microbial community composition and belowground ecosystem characteristics. Soil Biology and Biochemistry, 2008, 40(9): 2206-2216.

    [40] Bossio D A, Scow K M, Gunapala N, Graham K J. Determinants of soil microbial communities: Effects of agricultural management, season, and soil type on phospholipid fatty acid profiles. Microbial Ecology, 1998, 36(1): 1-12.

    [41] Zornoza R, Guerrero C, Mataix-Solera J, Scow K M, Arcenegui V, Mataix-Beneyto J. Changes in soil microbial community structure following the abandonment of agricultural terraces in mountainous areas of Eastern Spain. Applied Soil Ecology, 2009, 42(3): 315-323.

    [42] Schlesinger W H, Reynolds J, Cunningham G L, Huenneke L, Jarrell W, Virginia R, Whitford W. Biological feedbacks in global desertification. Science, 1990, 247(4946): 1043-1048.

    [43] Jackson R B, Caldwell M M. The scale of nutrient heterogeneity around individual plants and its quantification with geostatistics. Ecology, 1993, 74(2): 612-614.

    [44] Collins S L, Sinsabaugh R L, Crenshaw C, Green L, Porras-Alfaro A, Stursova M, Zeglin L H. Pulse dynamics and microbial processes in aridland ecosystems. Journal of Ecology, 2008, 96(3): 413-420.

    [45] 李新榮. 干旱沙區(qū)土壤空間異質(zhì)性變化對(duì)植被恢復(fù)的影響. 中國科學(xué)D輯: 地球科學(xué), 2005, 35(4): 361-370.

    [47] Osono T, Takeda H. Microfungi associated with Abies needles and Betula leaf litter in a subalpine coniferous forest. Canadian Journal of Microbiology, 2007, 53(1): 1-7.

    [48] Priha O, Grayston S J, Hiukka R, Pennanen T, Smolander A. Microbial community structure and characteristics of the organic matter in soils under Pinus sylvestris, Picea abies and Betula pendula at two forest sites. Biology and Fertility of Soils, 2001, 33(1): 17-24.

    [49] Wang M, Qu L Y, Ma K M, Yuan X. Soil microbial properties under different vegetation types on Mountain Han. Science China Life Sciences, 2013, 56(6): 561-570.

    [50] 馬淼, 李博, 陳家寬. 植物對(duì)荒漠生境的趨同適應(yīng). 生態(tài)學(xué)報(bào), 2006, 26(11): 3861-3869.

    [51] Johannes R, Philip C B, Erland B. The microbial PLFA composition as affected by pH in an arable soil. Soil Biology and Biochemistry, 2010, 42(3): 516-520.

    [52] B??th E, Anderson T H. Comparison of soil fungal/bacterial ratios in a pH gradient using physiological and PLFA-based techniques. Soil Biology and Biochemistry, 2003, 35(7): 955-963.

    [53] Drenovsky R E, Vo D, Graham K J, Scow K M. Soil water content and organic carbon availability are major determinants of soil microbial community composition. Microbial Ecology, 2004, 48(3): 424-430.

    [54] Demoling F, Figueroa D, Baath E. Comparison of factors limiting bacterial growth in different soils. Soil Biology and Biochemistry, 2007, 39(10): 2485-2495.

    [55] 黃元元. 岷江干旱河谷灌木恢復(fù)對(duì)土壤肥力質(zhì)量的影響 [D]. 北京: 中國科學(xué)院大學(xué), 2011.

    [56] Burke D J, Lopez-Gutierrez J C, Smemo K A, Chan C R. Vegetation and Soil Environment Influence the Spatial Distribution of Root-Associated Fungi in a Mature Beech-Maple Forest. Applied and Environmental Microbiology, 2009, 75(24): 7639-7648.

    [57] Bolan N S. A critical-review on the role of mycorrhizal fungi in the uptake of phosphorus by plants. Plant and Soil, 1991, 134(2): 189-207.

    [58] Ushio M, Balser T C, Kitayama K. Effects of condensed tannins in conifer leaves on the composition and activity of the soil microbial community in a tropical montane forest. Plant and Soil, 2013, 365(1/2): 157-170.

    [59] Briar S S, Fonte S J, Park I, Six J, Scow K, Ferris H. The distribution of nematodes and soil microbial communities across soil aggregate fractions and farm management systems. Soil Biology and Biochemistry, 2011, 43(5): 905-914.

    [60] Pietri J C A, Brookes P C. Substrate inputs and pH as factors controlling microbial biomass, activity and community structure in an arable soil. Soil Biology and Biochemistry, 2009, 41(7): 1396-1405.

    Effects of dominant shrubs on soil microbial composition in the semi-arid Minjiang River valley

    WANG Bingbing1,2, QU Laiye1,*, SONG Chengjun3, GONG Yuanbo4

    1ResearchCenterforEco-EnvironmentalSciences,ChineseAcademyofSciences,Beijing100085,China2UniversityofChineseAcademyofSciences,Beijing100049,China3ChineseAcademyofAgriculturalEngineering,Beijing100125,China4SichuanAgriculturalUniversity,Ya′an625014,China

    Although soil microorganisms play a critical role in soil organic matter decomposition, nutrient cycling and ecosystem productivity, we still know very little about how microbial communities are affected by environmental factors and how the structure and function of soil microbial communities influence key soil processes, especially for the arid and semi-arid river valley district. Studies focusing on soil microbial communities under this region are still under-represented when compared to moist forest and grassland ecosystems. Thus, the objectives of this work were to quantify the effects of shrubs on soil physicochemical properties and microbial composition (Phospholipid fatty acid, PLFA), and to find out soil factors that significantly affect soil microbial composition. To achieve this, we investigated the variations of soil properties between dominant shrubs and shrub-interspaces at a set of sites along the mountainous semi-arid valley of the upper Minjiang River, Sichuan Province. The dominant shrubs had significant effects on soil physicochemical properties and microbial composition, while the effects of shrub species seem to be trivial. The soil under shrub canopies had more organic carbon and nutrient resources, higher water content and lower bulk density than those of the shrub-interspaces. The bacterial communities′ abundance did not vary significantly between the two types of soil, while fungi and arbuscular mycorrhizal fungi (AMF) tended to dominate under shrubs. Furthermore, the correlation based Principal Component Analysis (PCA) demonstrated a clear dissimilarity of microbial community composition and physicochemical properties between shrub and shrub-interspace soils. Besides, the PCA revealed that the microbial community composition was a more sensitive and effective indicator than soil physiochemical properties in this study. Backward selection procedure and the combined redundancy analysis (RDA) were used to select the most important soil physiochemical variables affecting soil microbial communities. Though the main factors varied among the four sites, the nitrogen (N) associated factors (Total N, Available N and carbon-to-nitrogen ratio) always played a pivotal role in structuring the microbial groups, especially total microbial biomass (total PLFA) and the bacterial groups (Gran-positive (G+), Gran-negative (G-), general bacteria (B), and cyclo bacteria (cyclo)). The correlation between fungi groups (fungi, AMF) and N associated factors was relatively weak. Total phosphorus and carbon-to-phosphorus ratio strongly affected the fungi and AMF groups in Feihong and Mutuo. This may indicate the importance of phosphorus in shaping fungi and AMF communities, as well as the indispensable role that fungi and AMF communities may play in influencing phosphorus availability and cycling. In addition, soil water content exerted effects on the microbial composition in the two relative humid sites, Cuoji and Mutuo. Moreover, pH did not significantly affect microbial communities, though previous studies have pointed out its important role in affecting microbial composition. This may be due to the small variation of pH in this semi-arid river valley. Collectively, the enrichment of soil nutrients and better physical condition under shrub canopies were the main reasons for the difference of microbial community composition between shrubs and shrub-interspaces. Our work demonstrated the important role of shrubs in the semi-arid valley of upper Minjiang River, and provided insights into the nutrient cycling of semi-arid valley. It also revealed the importance of nitrogen and phosphorus for bacterial and fungal communities, respectively. Future studies on the plant-soil relationship in dryland ecosystems should pay more attention to fungi and AMF communities.

    Minjiang River; semi-arid valley; microbial community composition; PLFA; dominant shrubs; shrub-interspaces

    國家自然科學(xué)基金 (31170581); 國家自然科學(xué)基金青年基金(41101270)

    2013-06-10;

    日期:2014-05-16

    10.5846/stxb201306101570

    *通訊作者Corresponding author.E-mail: lyqu@rcees.ac.cn

    王冰冰, 曲來葉, 宋成軍, 宮淵波.岷江干旱河谷優(yōu)勢灌叢對(duì)土壤微生物群落組成的影響.生態(tài)學(xué)報(bào),2015,35(8):2481-2493.

    Wang B B, Qu L Y, Song C J, Gong Y B.Effects of dominant shrubs on soil microbial composition in the semi-arid Minjiang River valley.Acta Ecologica Sinica,2015,35(8):2481-2493.

    猜你喜歡
    飛虹灌叢空地
    “飛虹姑娘”帶百萬網(wǎng)友“云”旅游
    灌叢化對(duì)高寒草甸土壤水力性質(zhì)的影響
    吉蘭泰鹽湖綠洲檉柳灌叢生長與沙堆形態(tài)特征的關(guān)系
    看塔去
    第六章 神秘的空地
    第六章 神秘的空地
    空地一體化技術(shù)在IFTD試飛中的應(yīng)用
    荒漠化草原錦雞兒屬灌叢擴(kuò)增對(duì)牧草產(chǎn)量和植物多樣性的影響
    假如你有一塊空地,千萬別做庭院
    海峽姐妹(2017年11期)2018-01-30 08:57:42
    “石上飛虹”2017慶國慶·迎盛會(huì)雨花蛋白石展
    寶藏(2017年11期)2017-03-24 11:33:01
    美女免费视频网站| АⅤ资源中文在线天堂| 十八禁人妻一区二区| 国产视频内射| 黑人操中国人逼视频| 日韩中文字幕欧美一区二区| 一进一出好大好爽视频| 久久精品91无色码中文字幕| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩国产亚洲二区| 精品久久久久久,| 九九热线精品视视频播放| 日韩高清综合在线| 最近最新中文字幕大全电影3| 99国产精品一区二区三区| 欧美又色又爽又黄视频| 九九热线精品视视频播放| 国产真实乱freesex| 国产 一区 欧美 日韩| 成年版毛片免费区| 久久久成人免费电影| 叶爱在线成人免费视频播放| 免费看日本二区| 中文资源天堂在线| 无遮挡黄片免费观看| 90打野战视频偷拍视频| 国产日本99.免费观看| 美女被艹到高潮喷水动态| av天堂中文字幕网| 日韩欧美一区二区三区在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲av电影不卡..在线观看| 亚洲色图av天堂| 18禁黄网站禁片免费观看直播| 黄色视频,在线免费观看| 非洲黑人性xxxx精品又粗又长| 欧美一级a爱片免费观看看| 久久香蕉精品热| 久久久成人免费电影| 老鸭窝网址在线观看| 婷婷丁香在线五月| 国产精品一区二区精品视频观看| 国产精品一及| 又大又爽又粗| 欧美一级a爱片免费观看看| 国产极品精品免费视频能看的| 国产主播在线观看一区二区| 日韩高清综合在线| 国产精品电影一区二区三区| 神马国产精品三级电影在线观看| 精品国产乱码久久久久久男人| 免费在线观看亚洲国产| 99国产精品99久久久久| 亚洲精品在线观看二区| 美女黄网站色视频| 麻豆国产av国片精品| 久久人妻av系列| 麻豆成人av在线观看| av欧美777| 欧美zozozo另类| 特大巨黑吊av在线直播| 99国产精品99久久久久| 小说图片视频综合网站| 亚洲第一欧美日韩一区二区三区| 亚洲中文av在线| 亚洲熟妇熟女久久| 欧美日韩福利视频一区二区| 91九色精品人成在线观看| 99国产极品粉嫩在线观看| 亚洲av电影在线进入| 免费在线观看成人毛片| 国内精品一区二区在线观看| 国产精品久久久人人做人人爽| 亚洲国产精品成人综合色| 午夜福利成人在线免费观看| 午夜两性在线视频| 综合色av麻豆| 亚洲欧美精品综合久久99| 欧美乱色亚洲激情| 国产高清激情床上av| 男女视频在线观看网站免费| 国产精品一区二区免费欧美| 老熟妇仑乱视频hdxx| 亚洲精品粉嫩美女一区| 亚洲成a人片在线一区二区| 亚洲av成人一区二区三| 老熟妇仑乱视频hdxx| 成年女人永久免费观看视频| 亚洲成a人片在线一区二区| 国产亚洲精品一区二区www| 亚洲精品在线观看二区| 免费搜索国产男女视频| 中文字幕最新亚洲高清| 国产精品综合久久久久久久免费| 九九热线精品视视频播放| 99久久精品国产亚洲精品| 欧美日韩国产亚洲二区| 草草在线视频免费看| 麻豆av在线久日| 级片在线观看| 午夜视频精品福利| 午夜亚洲福利在线播放| 亚洲第一电影网av| 在线a可以看的网站| 国产私拍福利视频在线观看| 精品无人区乱码1区二区| 在线观看舔阴道视频| 久久久国产欧美日韩av| 欧美xxxx黑人xx丫x性爽| 国产在线精品亚洲第一网站| 亚洲精品在线观看二区| 嫩草影院入口| 亚洲 欧美一区二区三区| 精品电影一区二区在线| 18禁国产床啪视频网站| 精品国产三级普通话版| 久久伊人香网站| 精品午夜福利视频在线观看一区| 久久国产精品人妻蜜桃| 欧美乱码精品一区二区三区| 1000部很黄的大片| 美女大奶头视频| 深夜精品福利| 免费大片18禁| 亚洲国产色片| 午夜精品一区二区三区免费看| 少妇熟女aⅴ在线视频| 老熟妇乱子伦视频在线观看| 中文资源天堂在线| 午夜福利在线在线| 久久精品91无色码中文字幕| 少妇的逼水好多| 午夜影院日韩av| 国产高清有码在线观看视频| 夜夜躁狠狠躁天天躁| 午夜福利高清视频| 欧美激情在线99| av黄色大香蕉| 国产精品久久久久久亚洲av鲁大| 最近最新中文字幕大全免费视频| 精品一区二区三区av网在线观看| 国产野战对白在线观看| www.www免费av| 国内精品久久久久精免费| 欧美+亚洲+日韩+国产| 最近视频中文字幕2019在线8| 久久精品国产99精品国产亚洲性色| 国产精品一区二区三区四区免费观看 | 少妇的丰满在线观看| 国产精品一区二区三区四区免费观看 | 老汉色av国产亚洲站长工具| 欧美中文日本在线观看视频| 夜夜看夜夜爽夜夜摸| 久久精品综合一区二区三区| 午夜成年电影在线免费观看| 每晚都被弄得嗷嗷叫到高潮| 91在线精品国自产拍蜜月 | 在线观看免费视频日本深夜| or卡值多少钱| 日韩欧美在线二视频| 麻豆久久精品国产亚洲av| 国产精品爽爽va在线观看网站| 久久草成人影院| 三级男女做爰猛烈吃奶摸视频| 五月玫瑰六月丁香| 国产v大片淫在线免费观看| 三级男女做爰猛烈吃奶摸视频| 床上黄色一级片| 狂野欧美激情性xxxx| 亚洲精品一卡2卡三卡4卡5卡| 九色成人免费人妻av| 最新在线观看一区二区三区| 亚洲 欧美 日韩 在线 免费| 免费观看人在逋| 国内久久婷婷六月综合欲色啪| 三级国产精品欧美在线观看 | 麻豆成人av在线观看| 久久精品综合一区二区三区| 一本精品99久久精品77| 国产野战对白在线观看| 亚洲色图 男人天堂 中文字幕| 99国产精品99久久久久| 村上凉子中文字幕在线| 国产精华一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 1024香蕉在线观看| 欧美中文日本在线观看视频| 亚洲成人免费电影在线观看| a级毛片a级免费在线| 精品日产1卡2卡| 日本成人三级电影网站| 国产精品av久久久久免费| 亚洲avbb在线观看| 日本免费一区二区三区高清不卡| www.www免费av| 99久久99久久久精品蜜桃| 久久婷婷人人爽人人干人人爱| 日韩有码中文字幕| 在线十欧美十亚洲十日本专区| 日本精品一区二区三区蜜桃| 老司机午夜福利在线观看视频| 亚洲中文字幕日韩| 久久久久久大精品| 观看免费一级毛片| 亚洲国产精品sss在线观看| 给我免费播放毛片高清在线观看| 老熟妇乱子伦视频在线观看| 国产精品女同一区二区软件 | 亚洲av电影在线进入| 日本三级黄在线观看| ponron亚洲| 国产成人精品无人区| 成人一区二区视频在线观看| 国产精品久久久人人做人人爽| 精品无人区乱码1区二区| 在线观看66精品国产| 国产精品九九99| 精品日产1卡2卡| 国产精品综合久久久久久久免费| 高清毛片免费观看视频网站| 欧美绝顶高潮抽搐喷水| 床上黄色一级片| 欧美中文综合在线视频| 国产男靠女视频免费网站| 欧美xxxx黑人xx丫x性爽| 久久精品国产综合久久久| 国产蜜桃级精品一区二区三区| 99久久国产精品久久久| 动漫黄色视频在线观看| 免费搜索国产男女视频| 欧美绝顶高潮抽搐喷水| 90打野战视频偷拍视频| 脱女人内裤的视频| 久久天堂一区二区三区四区| 欧美一级毛片孕妇| 久久欧美精品欧美久久欧美| 桃红色精品国产亚洲av| 婷婷精品国产亚洲av| 成年女人永久免费观看视频| 观看免费一级毛片| 国产真人三级小视频在线观看| 操出白浆在线播放| 日本黄色视频三级网站网址| 亚洲黑人精品在线| 欧美日韩国产亚洲二区| 99在线人妻在线中文字幕| 午夜成年电影在线免费观看| 久久久久亚洲av毛片大全| 亚洲五月婷婷丁香| 日韩精品中文字幕看吧| 99国产综合亚洲精品| 久久久久久久久免费视频了| 婷婷六月久久综合丁香| 人妻久久中文字幕网| www.自偷自拍.com| 久99久视频精品免费| 午夜精品一区二区三区免费看| 午夜a级毛片| 国产成人一区二区三区免费视频网站| 在线a可以看的网站| 免费观看人在逋| 美女被艹到高潮喷水动态| 久久亚洲精品不卡| 久久精品国产99精品国产亚洲性色| 国产高清视频在线播放一区| 男人的好看免费观看在线视频| 岛国在线观看网站| 亚洲精品美女久久av网站| 国产激情久久老熟女| 日韩精品青青久久久久久| 18禁裸乳无遮挡免费网站照片| 中文资源天堂在线| 琪琪午夜伦伦电影理论片6080| 国产av一区在线观看免费| 国产三级在线视频| 精品一区二区三区视频在线 | 国产精品av视频在线免费观看| 国产精品美女特级片免费视频播放器 | 一个人免费在线观看的高清视频| 亚洲人成网站高清观看| 国产精品av视频在线免费观看| 亚洲一区高清亚洲精品| 久久欧美精品欧美久久欧美| 久99久视频精品免费| 夜夜看夜夜爽夜夜摸| 免费人成视频x8x8入口观看| www日本在线高清视频| 国产精品九九99| 国产高清三级在线| 1024香蕉在线观看| 国产精品永久免费网站| 国产又色又爽无遮挡免费看| 女人被狂操c到高潮| 国产精品综合久久久久久久免费| 久9热在线精品视频| 午夜两性在线视频| 99riav亚洲国产免费| 亚洲va日本ⅴa欧美va伊人久久| 女生性感内裤真人,穿戴方法视频| 黄色丝袜av网址大全| 国产一区在线观看成人免费| 女人被狂操c到高潮| 伊人久久大香线蕉亚洲五| 两性夫妻黄色片| 国产精品99久久99久久久不卡| 性色av乱码一区二区三区2| 真人做人爱边吃奶动态| 亚洲精品美女久久av网站| 制服人妻中文乱码| 国产精品,欧美在线| 亚洲男人的天堂狠狠| 午夜影院日韩av| 99久久国产精品久久久| 男人舔奶头视频| 后天国语完整版免费观看| 亚洲人成网站在线播放欧美日韩| 性色avwww在线观看| 免费无遮挡裸体视频| 一夜夜www| 在线国产一区二区在线| 欧美性猛交黑人性爽| 午夜精品一区二区三区免费看| 国产91精品成人一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 女人高潮潮喷娇喘18禁视频| 在线免费观看不下载黄p国产 | 国产97色在线日韩免费| 欧美又色又爽又黄视频| aaaaa片日本免费| 国产高潮美女av| 亚洲色图 男人天堂 中文字幕| 亚洲成人免费电影在线观看| 国产亚洲av嫩草精品影院| 国产视频一区二区在线看| 91久久精品国产一区二区成人 | 国产精品久久久久久久电影 | 一卡2卡三卡四卡精品乱码亚洲| 亚洲欧美日韩无卡精品| 99久久综合精品五月天人人| 日韩 欧美 亚洲 中文字幕| 久久精品国产99精品国产亚洲性色| av中文乱码字幕在线| 亚洲成人久久爱视频| 99久久久亚洲精品蜜臀av| 国产亚洲精品久久久com| 成年版毛片免费区| 国产精品久久视频播放| 国产在线精品亚洲第一网站| 亚洲熟妇中文字幕五十中出| 天堂网av新在线| 国产精品一及| 91老司机精品| 亚洲av成人av| 啦啦啦韩国在线观看视频| 日韩欧美三级三区| 国产在线精品亚洲第一网站| АⅤ资源中文在线天堂| 一区二区三区国产精品乱码| 日韩欧美精品v在线| 精品免费久久久久久久清纯| 超碰成人久久| 亚洲国产精品成人综合色| 曰老女人黄片| 一边摸一边抽搐一进一小说| 国产亚洲精品久久久久久毛片| 男女视频在线观看网站免费| 小说图片视频综合网站| 亚洲成av人片免费观看| 麻豆久久精品国产亚洲av| 亚洲欧美日韩高清专用| 亚洲成人久久性| 一本一本综合久久| 男人舔女人的私密视频| 精品无人区乱码1区二区| 欧美又色又爽又黄视频| 欧美中文日本在线观看视频| 国产不卡一卡二| 国产精品久久电影中文字幕| 99久久精品热视频| 国产在线精品亚洲第一网站| 国产午夜精品久久久久久| 脱女人内裤的视频| 午夜福利免费观看在线| 搡老熟女国产l中国老女人| 国产精品女同一区二区软件 | 最近最新免费中文字幕在线| 欧美日韩综合久久久久久 | 亚洲va日本ⅴa欧美va伊人久久| 香蕉丝袜av| 免费在线观看视频国产中文字幕亚洲| 午夜亚洲福利在线播放| 国语自产精品视频在线第100页| 亚洲国产精品成人综合色| 国产精品99久久99久久久不卡| 色播亚洲综合网| 变态另类成人亚洲欧美熟女| 国模一区二区三区四区视频 | 久9热在线精品视频| 麻豆久久精品国产亚洲av| 成人av在线播放网站| 国产爱豆传媒在线观看| xxxwww97欧美| 18禁黄网站禁片免费观看直播| 成人亚洲精品av一区二区| 国产又色又爽无遮挡免费看| 午夜久久久久精精品| 18禁国产床啪视频网站| 窝窝影院91人妻| 亚洲国产欧美人成| 欧美日本视频| а√天堂www在线а√下载| 老熟妇乱子伦视频在线观看| 桃红色精品国产亚洲av| 国内少妇人妻偷人精品xxx网站 | 国产伦精品一区二区三区视频9 | 欧美乱码精品一区二区三区| 岛国在线免费视频观看| 黄色成人免费大全| 99在线视频只有这里精品首页| 男人舔女人的私密视频| 国产精品久久久人人做人人爽| 老司机午夜十八禁免费视频| 老汉色∧v一级毛片| 给我免费播放毛片高清在线观看| 国产高清视频在线播放一区| 成年免费大片在线观看| 国产不卡一卡二| 伦理电影免费视频| 黄片小视频在线播放| 高潮久久久久久久久久久不卡| 床上黄色一级片| 一级毛片高清免费大全| 免费观看人在逋| 757午夜福利合集在线观看| 一个人观看的视频www高清免费观看 | 2021天堂中文幕一二区在线观| 欧美丝袜亚洲另类 | 久久伊人香网站| 久久天堂一区二区三区四区| 美女被艹到高潮喷水动态| 听说在线观看完整版免费高清| 好看av亚洲va欧美ⅴa在| 国产成人精品久久二区二区91| 麻豆成人av在线观看| tocl精华| 午夜免费成人在线视频| 久久草成人影院| 黄色女人牲交| 丁香六月欧美| 我要搜黄色片| 女同久久另类99精品国产91| 国产精品美女特级片免费视频播放器 | 亚洲精华国产精华精| 午夜两性在线视频| 麻豆国产av国片精品| 嫩草影视91久久| 成人鲁丝片一二三区免费| 亚洲精品一区av在线观看| 国产精品美女特级片免费视频播放器 | 久久久久九九精品影院| 欧美一区二区精品小视频在线| 国产成人av教育| 精华霜和精华液先用哪个| 99热这里只有精品一区 | 欧美极品一区二区三区四区| 亚洲av免费在线观看| 一级毛片精品| 欧美成人性av电影在线观看| 九九在线视频观看精品| 中文字幕精品亚洲无线码一区| 啦啦啦免费观看视频1| 国内久久婷婷六月综合欲色啪| 美女高潮喷水抽搐中文字幕| 亚洲国产日韩欧美精品在线观看 | 欧美中文综合在线视频| 一二三四社区在线视频社区8| 91在线观看av| 99久久无色码亚洲精品果冻| 12—13女人毛片做爰片一| 久久久久性生活片| 国产免费av片在线观看野外av| 看免费av毛片| 天天添夜夜摸| 在线观看66精品国产| 成人性生交大片免费视频hd| 五月伊人婷婷丁香| 日本一本二区三区精品| 久久久精品欧美日韩精品| 亚洲色图av天堂| 精品一区二区三区视频在线 | а√天堂www在线а√下载| 午夜两性在线视频| 欧美色欧美亚洲另类二区| 草草在线视频免费看| 可以在线观看毛片的网站| 国内毛片毛片毛片毛片毛片| 国产精品亚洲一级av第二区| 99久久国产精品久久久| 热99在线观看视频| 国产1区2区3区精品| 岛国在线免费视频观看| 免费在线观看影片大全网站| 国产v大片淫在线免费观看| 午夜免费观看网址| 亚洲成av人片在线播放无| a级毛片在线看网站| 国内少妇人妻偷人精品xxx网站 | 哪里可以看免费的av片| 亚洲最大成人中文| 麻豆国产av国片精品| 日韩成人在线观看一区二区三区| av黄色大香蕉| 国产精品精品国产色婷婷| 麻豆成人av在线观看| 国产成人一区二区三区免费视频网站| 一级毛片高清免费大全| 欧美zozozo另类| 悠悠久久av| 91麻豆av在线| 哪里可以看免费的av片| 日韩精品中文字幕看吧| 亚洲人成电影免费在线| 美女午夜性视频免费| 中文字幕久久专区| 搡老妇女老女人老熟妇| 成年女人看的毛片在线观看| 亚洲成a人片在线一区二区| 久久天堂一区二区三区四区| 欧美日本视频| 日本一二三区视频观看| 在线视频色国产色| 99久久成人亚洲精品观看| 90打野战视频偷拍视频| 日本成人三级电影网站| 桃色一区二区三区在线观看| 国产精品久久久久久精品电影| 午夜日韩欧美国产| 99riav亚洲国产免费| 人人妻人人澡欧美一区二区| 在线播放国产精品三级| 亚洲无线在线观看| 亚洲人成网站在线播放欧美日韩| 美女高潮的动态| 亚洲精品色激情综合| 超碰成人久久| 一级作爱视频免费观看| 国产主播在线观看一区二区| 人人妻人人看人人澡| 又大又爽又粗| 无遮挡黄片免费观看| 国产成人av激情在线播放| 亚洲七黄色美女视频| 欧美日韩黄片免| 精品日产1卡2卡| 亚洲精品国产精品久久久不卡| 欧美日韩精品网址| 国产麻豆成人av免费视频| 丰满的人妻完整版| 国产亚洲精品综合一区在线观看| 欧美xxxx黑人xx丫x性爽| 亚洲精品一区av在线观看| 日本一本二区三区精品| 91在线精品国自产拍蜜月 | 午夜福利免费观看在线| 99视频精品全部免费 在线 | 亚洲av美国av| 少妇熟女aⅴ在线视频| 欧美日韩一级在线毛片| 国产蜜桃级精品一区二区三区| 久久精品国产综合久久久| 禁无遮挡网站| 丁香欧美五月| 成人一区二区视频在线观看| 亚洲五月天丁香| 人人妻人人澡欧美一区二区| 亚洲欧美日韩高清在线视频| 啪啪无遮挡十八禁网站| 亚洲真实伦在线观看| 制服人妻中文乱码| 亚洲精品中文字幕一二三四区| 日韩免费av在线播放| 99久久成人亚洲精品观看| 免费大片18禁| 久久久久久人人人人人| 欧美中文综合在线视频| 听说在线观看完整版免费高清| 亚洲精品在线美女| 精品久久久久久成人av| 舔av片在线| 亚洲一区二区三区色噜噜| 老司机福利观看| 天堂√8在线中文| 好看av亚洲va欧美ⅴa在| 母亲3免费完整高清在线观看| 久久精品91蜜桃| 国产精品香港三级国产av潘金莲| 午夜精品久久久久久毛片777| 国产男靠女视频免费网站| 首页视频小说图片口味搜索| 国产免费av片在线观看野外av| 脱女人内裤的视频| 九九热线精品视视频播放| 国产伦一二天堂av在线观看| 99国产精品一区二区蜜桃av| 亚洲中文字幕日韩| 精品无人区乱码1区二区| 国产一区在线观看成人免费| 男女那种视频在线观看| 露出奶头的视频| 手机成人av网站| 99精品久久久久人妻精品| 色噜噜av男人的天堂激情| 成人欧美大片| 久久亚洲真实| 波多野结衣高清无吗|