趙愛花,劉 蕾,付 偉,武 慧,陳保冬

1 中國科學(xué)院生態(tài)環(huán)境研究中心城市與區(qū)域生態(tài)國家重點(diǎn)實(shí)驗(yàn)室, 北京 100085 2 中國科學(xué)院大學(xué), 北京 100049 3 南京信息工程大學(xué)生態(tài)研究院,江蘇省農(nóng)業(yè)氣象重點(diǎn)實(shí)驗(yàn)室, 南京 210044

自工業(yè)革命以來,大量化石燃料的燃燒,以及農(nóng)業(yè)和工業(yè)生產(chǎn)活動(dòng)中氮(N)素消耗大量增加,加速了全球大氣N沉降并導(dǎo)致一系列生態(tài)環(huán)境問題。據(jù)估計(jì),1995年全球平均N沉降速率大約是1860年的3倍,而2050全球平均N沉降速率將增加到1860年的6倍[1]。我國陸地生態(tài)系統(tǒng)N沉降速率已達(dá)21.1 kg hm-2a-1[2],其中,濕沉降速率已經(jīng)由20紀(jì)90 年代的11.11 kg hm-2a-1增加到13.87 kg hm-2a-1[3],增長了近25%,在今后幾十年中我國N沉降的形勢將更加嚴(yán)峻[2]。N沉降增加一方面可以提高受N限制的生態(tài)系統(tǒng)的生產(chǎn)力[4],另一方面過多的N素供應(yīng)也會(huì)給生態(tài)系統(tǒng)帶來諸多負(fù)面影響,如土壤酸化[5]、植物物種多樣性喪失[6]、水體富營養(yǎng)化[7]等。同時(shí),N沉降還可以通過改變土壤環(huán)境(如N的有效性、土壤酸化等)直接影響土壤微生物多樣性[8- 9],或通過地上植被的生理生態(tài)響應(yīng)間接作用于土壤微生物[10]。

叢枝菌根(Arbuscular mycorrhiza,AM)真菌是一類在陸地生態(tài)系統(tǒng)中廣泛存在的土壤真菌,能與絕大多數(shù)陸地植物形成互惠共生體系,具有許多重要的生態(tài)功能,如幫助植物獲取N、P等礦質(zhì)養(yǎng)分和水分,增強(qiáng)植物對(duì)生物和非生物脅迫的適應(yīng)能力[11],改變植物群落的組成、多樣性,還能影響生態(tài)系統(tǒng)的生產(chǎn)力和結(jié)構(gòu)功能穩(wěn)定性[12- 13]。在陸地生態(tài)系統(tǒng)N素循環(huán)過程中,AM真菌除了能夠吸收、轉(zhuǎn)運(yùn)N素以外,在N素生物固定、硝化、反硝化以及N素淋洗過程中均具有潛在重要作用[14- 15]。AM真菌一方面通過改善植物磷(P)營養(yǎng),提高豆科植物的固N(yùn)能力[14],另一方面可能通過影響固N(yùn)微生物功能基因的表達(dá)來調(diào)控N素的生物固定[16]。AM真菌還可能通過與氨氧化微生物競爭底物[17],改變氨氧化微生物的群落組成[18]及植物根系營養(yǎng)生理[19]調(diào)控硝化過程,通過調(diào)控反硝化細(xì)菌的數(shù)量[20]及相關(guān)功能基因的表達(dá)進(jìn)而影響反硝化過程[16]。此外,AM真菌還能夠通過改善土壤結(jié)構(gòu)[21],促進(jìn)宿主植物對(duì)無機(jī)N的吸收等多種途徑減少因淋洗而造成的土壤N素?fù)p失。

鑒于N沉降形勢的嚴(yán)峻性和AM真菌重要的生態(tài)功能,越來越多的土壤生態(tài)學(xué)家開始關(guān)注N沉降對(duì)AM真菌的影響并開展了試驗(yàn)研究。目前,已有相關(guān)研究多采用N素添加的方式模擬N沉降,受到具體試驗(yàn)條件的影響,研究結(jié)果不甚一致：N沉降對(duì)AM真菌多樣性可能有積極作用[22],沒有顯著影響[23],或者具有消極作用[24]；類似地,對(duì)于AM真菌群落組成而言,有的研究顯示N沉降有顯著影響[25- 28],而另外一些研究則顯示無顯著影響[29- 30]。影響N沉降生態(tài)效應(yīng)的因素有很多,如施N速率[22]、N素形態(tài)[31]、試驗(yàn)方法[29]、試驗(yàn)周期[32]等。此外,還可能受到環(huán)境中N的背景值及生態(tài)系統(tǒng)類型的影響[33]。因此,客觀認(rèn)識(shí)和評(píng)價(jià)N沉降對(duì)AM真菌的影響,還需要開展更為系統(tǒng)的研究工作。

總體上,目前N沉降對(duì)AM真菌多樣性和群落組成影響的研究多集中在草地生態(tài)系統(tǒng),對(duì)森林生態(tài)系統(tǒng)的研究相對(duì)匱乏。在森林生態(tài)系統(tǒng)中開展的有限研究中,施N方式又多為林下施N[24- 25,29],這種方式忽略了冠層發(fā)生的一系列生態(tài)過程,如N的吸收、固定、揮發(fā)、轉(zhuǎn)化等[34-38],而這些過程都會(huì)影響到達(dá)地面的N的形態(tài)和數(shù)量[39],因此林下施N可能無法準(zhǔn)確反映自然條件下N沉降對(duì)AM真菌多樣性和群落組成的影響。Huang等[40]研究發(fā)現(xiàn),林冠施N對(duì)土壤微生物生物量和群落組成的影響小于林下施N,那么對(duì)AM真菌來說是否也是類似的情況,目前尚無相關(guān)報(bào)道?；诖?我們依托雞公山林冠模擬氮沉降野外試驗(yàn)平臺(tái),通過連續(xù)4 a的樣品采集和分析,嘗試回答這個(gè)問題。我們提出的研究假設(shè)是：林下施N模擬N沉降會(huì)高估自然N沉降對(duì)森林土壤AM真菌多樣性和群落組成的影響。

1 材料與方法

1.1 研究區(qū)概況

雞公山林冠模擬N沉降野外試驗(yàn)平臺(tái)于2013年建于河南省信陽市以南 38 km的雞公山國家級(jí)自然保護(hù)區(qū)(31°46′—31°52′N, 114°01′—114°06′E)的原生林中。研究區(qū)由于地處北亞熱帶邊緣秦嶺山系西端的淺山區(qū),受東亞季風(fēng)氣候的影響,具有典型的北亞熱帶向暖溫帶過渡的季風(fēng)氣候和山地氣候特征,四季分明。研究區(qū)內(nèi)土壤以黃棕壤、黃褐土為主,pH 值在 5.0到6.0 之間。植被類型為落葉闊葉混交林,林齡約45 a。優(yōu)勢喬木樹種為麻櫟(QuercusacutissimaCarruth.)、栓皮櫟(QuercusvariabilisBl.)和楓香(LiquidambarformosanaHance),林下樹種和草本豐富。

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

實(shí)驗(yàn)采用完全隨機(jī)區(qū)組設(shè)計(jì),包括4個(gè)區(qū)組(相當(dāng)于4個(gè)試驗(yàn)重復(fù)),每個(gè)區(qū)組隨機(jī)設(shè)置5個(gè)半徑為17 m的樣方,對(duì)應(yīng)5個(gè)試驗(yàn)處理：對(duì)照(CK)、林冠施N 25 kg hm-2a-1(CN25)或50 kg hm-2a-1(CN50)、林下施N 25 kg hm-2a-1(UN25)或50 kg hm-2a-1(UN50)。為防止處理間的干擾,各樣方之間留有至少20 m的緩沖帶,每個(gè)緩沖帶中間都加裝深度為 1 m 的PVC隔離板。

N添加在每年4到10月份進(jìn)行,每月一次,一年7次。添加的N素形態(tài)為NH4NO3溶液。每次施N時(shí),添加的NH4NO3溶液相當(dāng)于3 mm的降水量,一年的添加量相當(dāng)于21 mm的降水量,小于該地區(qū)年均降雨量的2%,因此水分的影響可以忽略不計(jì)[41]。林冠施N和林下施N都是由位于相應(yīng)樣方中心的自動(dòng)噴灑裝置來實(shí)現(xiàn)的。林冠噴灑裝置高35 m(冠層以上約5 m),林下噴灑裝置距離地面的高度為1.5 m。該裝置頂端安裝有一套搖臂噴頭,有供水管道與樣地外的蓄水池連接,利用變頻調(diào)速恒壓噴灌設(shè)備提供壓力,驅(qū)動(dòng)搖臂噴頭360° 旋轉(zhuǎn),以保證噴灑的均勻性和精確性。塔基處供水管道安裝有水量計(jì),可精準(zhǔn)控制每次處理的用水量。關(guān)于裝置的細(xì)節(jié)、工作方式和效率參考Zhang等[41]。

1.3 樣品采集與分析

從2013年到2016年,每年最后一次N添加處理完成之后,進(jìn)行土壤樣品采集。采樣時(shí)間為11月中下旬到12月初,具體時(shí)間視當(dāng)?shù)靥鞖馇闆r而定。土壤樣品的采集和分析,以及生物信息學(xué)分析均參考Zhao等[32]。為了盡可能減少采樣位置的影響,在各樣方核心區(qū)20 m×20 m范圍內(nèi)選取5棵樹作為目標(biāo)樹。在距每棵目標(biāo)樹1—2 m范圍內(nèi),取2鉆土(直徑3 cm,深度10 cm),每個(gè)樣方共10鉆土,制備成一個(gè)混合土樣,用冰袋運(yùn)送到實(shí)驗(yàn)室?；旌贤翗舆^2 mm篩之后分成兩部分,一部分冷凍干燥后儲(chǔ)存于-80℃冰箱,用于土壤微生物總DNA的提取；另一部分自然風(fēng)干,用于測定土壤pH、有效磷(available phosphorus,AP)等基本理化性質(zhì)。

土壤有效磷用0.5 mol/L的NaHCO3提取,鉬藍(lán)比色法測定[42]；土壤pH值采用水土比2.5∶1,PB- 10 pH計(jì)(Sartorius, G?ttingen, Germany)測定；銨態(tài)N、硝態(tài)N首先用2 mol/L 的KCl浸提(水土比為5∶1),然后用連續(xù)流動(dòng)分析儀進(jìn)行測定(SAN++, Skakar, Breda, Holland)；銨態(tài)氮、硝態(tài)氮含量之和作為有效N含量(available nitrogen, AN)；有效N與有效P之比即為N/P比。

采用Power-Soil RDNA提取試劑盒 (MO BIO Laboratories, San Diego, CA, United States),按照說明書提取土壤DNA。DNA樣品用滅菌超純水稀釋5倍,然后進(jìn)行巢式PCR擴(kuò)增。擴(kuò)增的引物分別為AML1/AML2[43]和AMV4.5NF/AMDGR[44],擴(kuò)增體系及程序設(shè)置與Zhao等[32]相同：第一次PCR反應(yīng)總體積為25 μL,其中含有2.5 μL 10 × Ex-Taq緩沖液(Mg2+Plus)、2.0 μL dNTP混合物、0.25 μL Ex-Taq(5u/mL)(Takara,大連,中國)、1.0 μL(10 mg/mL)BSA(TaKaRa,大連,中國),雙端引物各0.5 μL(10 μmol/L),17.25 μL 無菌水和1.0 μL DNA模板。PCR擴(kuò)增條件如下：94℃預(yù)變性3 min；94℃ 45 s,51℃ 40 s,72℃ 1 min,擴(kuò)增35個(gè)循環(huán)；然后72℃延伸10 min,16℃冷卻2分鐘。第一輪的PCR產(chǎn)物進(jìn)行10倍稀釋,作為第二輪的模板。第二輪擴(kuò)增體系與第一輪相同,擴(kuò)增條件如下：94℃預(yù)變性3 min；94℃ 40 s,58℃ 1 min,72℃ 1 min,擴(kuò)增35個(gè)循環(huán)；然后72℃延伸10 min,4℃冷卻2 min。擴(kuò)增產(chǎn)物用1%瓊脂糖凝膠電泳檢測后,送上海瀚宇科技有限公司切膠回收、純化、建庫及測序(Illumina Miseq PE300測序平臺(tái))。

最初的數(shù)據(jù)質(zhì)量控制和序列拼接由上海瀚宇生物科技有限公司完成。原始測序數(shù)據(jù)用 Trimmomatic v 0.32[45]進(jìn)行質(zhì)控,主要步驟如下：(1)刪除帶有N堿基的reads；(2)刪除低質(zhì)量的reads Q value<20)；(3)刪除過短的reads(<50 bp)及其配對(duì)reads。用Mothur v.1.32.1[46]進(jìn)行序列拼接,允許引物有1個(gè)堿基的錯(cuò)配。刪除拼接同聚物長度>8,序列總長<200 bp的序列。嵌合體的過濾由Chimera.uchime v.4.2[47]完成。使用Usearch v.9.0.2132_i86linux32[48]在97%的相似性水平劃分OTU (operational taxonomic unit,操作分類單元)。刪除序列數(shù)<5條的OTUs。其余OTUs的代表序列和NCBI GenBank數(shù)據(jù)庫及MaarjAM數(shù)據(jù)庫進(jìn)行比對(duì),滿足下列條件的OTUs用于后續(xù)分析：Identifications >97% similarity,>90% coverage 及>200 BLAST score value,不滿足條件者被舍去。為了進(jìn)一步確認(rèn)分類信息的可靠性,我們把OTUs的代表序列和在GenBank數(shù)據(jù)庫中比對(duì)到的序列放在一起,采用MEGA v5[49],依據(jù)Kimura 2-parameter模型構(gòu)建了Neighbor-Joining 系統(tǒng)發(fā)育樹。除了bootstrap數(shù)目設(shè)為1000以外,其余參數(shù)均采用默認(rèn)值。

1.4 統(tǒng)計(jì)分析

所有的統(tǒng)計(jì)分析均采用R 3.3.2 (R Core Team, 2016) 完成。采用中位數(shù)稀釋法對(duì)每年的樣品測序數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化[50]。標(biāo)準(zhǔn)化后的OTU表格進(jìn)行Hellinger轉(zhuǎn)化,作為分析 AM真菌群落組成的數(shù)據(jù),也用來計(jì)算AM真菌群落的alpha多樣性指數(shù)(豐富度S、香儂指數(shù) H、Pielou′s 均勻度指數(shù)J和辛普森指數(shù))。用“vegan”包[51]“specnumber”函數(shù)計(jì)算豐富度。香儂指數(shù)和辛普森指數(shù)的計(jì)算采用“vegan”包“diversity”函數(shù)完成。Pielou′s 均勻度指數(shù)J=H/log(S)。采用“aov”函數(shù)進(jìn)行雙因素方差分析,評(píng)估施N方式和施N速率對(duì)土壤因子和AM真菌alpha多樣性指數(shù)的影響,分析之前用“stats”包“shapiro.test”函數(shù)進(jìn)行正態(tài)性檢驗(yàn),不符合要求的數(shù)據(jù)進(jìn)行1次或2次乘方、開方、log或 Box-Cox轉(zhuǎn)化,若轉(zhuǎn)化后仍不滿足正態(tài)分布,則采用“friedman.test” 函數(shù)進(jìn)行Friedman檢驗(yàn)。

施N方式和施N速率對(duì)AM真菌群落組成的影響及不同處理間的兩兩比較,均采用 “vegan” 包“adonis2” 函數(shù)進(jìn)行PERMANOVA分析,距離矩陣為“Bray-Curtis”距離,置換次數(shù)設(shè)為9999 次。為了實(shí)現(xiàn)數(shù)據(jù)的可視化,使用“vegan”包“metaMDS”函數(shù)進(jìn)行NMDS分析并作圖。

2 結(jié)果

2.1 施氮對(duì)土壤理化性質(zhì)的影響

在2013年,施N速率對(duì)土壤pH值有顯著影響(P=0.024)。兩個(gè)施N速率與對(duì)照相比,pH值均顯著降低,但是兩個(gè)施N速率之間沒有顯著差異。施N方式對(duì)土壤pH沒有顯著影響(P=0.826),而且施N方式和施N速率之間也沒有顯著交互作用(P=0.929)(表1)。對(duì)于其他土壤性質(zhì)而言,施N方式、施N速率的影響及二者交互作用均不顯著(表1)。在2014年,施N方式和施N速率對(duì)所有土壤性質(zhì)參數(shù)均無顯著影響(表1)。在2015年,施N方式顯著影響了土壤pH值(P=0.021)(表1),林下施N土壤pH值(4.21)顯著低于林冠施N處理(4.35)。對(duì)于其他土壤性質(zhì)而言,施N方式、施N速率的影響以及二者交互作用均不顯著(表1)。2016年和2013年結(jié)果基本相同。

表1 施N方式和施N速率對(duì)土壤性質(zhì)的影響

2.2 高通量測序基本信息及AM真菌OTUs分類情況

對(duì)連續(xù)4 a共80個(gè)土壤DNA樣品進(jìn)行高通量測序,初步得到761939條序列,劃分為661個(gè)OTUs。經(jīng)過嚴(yán)格篩選以后,最終確認(rèn)AM真菌的OTUs為181個(gè),共604890條序列,隸屬于1門,1綱,4目,8科,9屬,50VT(virtual taxa)。其中,絕大多數(shù)的OTUs和序列屬于Glomus屬,其次為Acaulospora屬(圖1)。OTUs的代表序列已上傳到Genbank數(shù)據(jù)庫(MN559107-MN559287)。序列數(shù)過少的3個(gè)樣品(2014年區(qū)組3的CN50處理 33條；2015年區(qū)組1的UN50處理 194條； 2016 年區(qū)組4的UN25處理 60條),在后續(xù)的分析中沒有被包括進(jìn)來,其余樣品的序列數(shù)均多于1957條。

圖1 所有土壤樣品AM真菌序列和操作分類單元的分布情況 Fig.1 The proportional distributions of AM fungal sequences and derived AM fungal OTUs detected in all soil samples

2.3 施氮對(duì)AM真菌alpha多樣性指數(shù)和群落組成的影響

從2013年到2016年,施N方式、施N速率及二者之間的交互作用對(duì)AM真菌alpha多樣性指數(shù)都沒有顯著影響(表2)。隨著處理時(shí)間的延長,不同處理下AM真菌的alpha多樣性指數(shù)有趨同的趨勢(表2)。

表2 不同年份下施N方式和施N速率對(duì)AM真菌alpha多樣性指數(shù)的影響

在2013年,施N方式對(duì)AM真菌群落組成有輕微影響(P=0.052),而施N速率有極顯著的影響(P=0.004),且二者之間有顯著交互作用(P=0.045)(表3)：施N速率為25 kg hm-2a-1時(shí),林冠施N(CN25)與對(duì)照差異不顯著(P=0.115),林下施N(UN25)與對(duì)照差異極顯著(P=0.005),且林冠施N(CN25)和林下施N(UN25)相比,群落組成差異極顯著(P=0.002)(表4,圖2)；當(dāng)施N速率為50 kg hm-2a-1時(shí),林冠施N(CN50)、林下施N(UN50)及對(duì)照三者之間均無顯著差異(表4,圖2)。從2014年到2016年,施N方式和施N速率對(duì)AM真菌群落組成都沒有顯著影響,雙因素交互作用也不顯著(表3,圖2)。

表3 不同年份施N方式和施N速率對(duì)AM真菌群落組成的影響

圖2 不同處理下AM真菌群落組成Fig.2 The variation of AM fungal community composition among different treatments不同顏色的橢圓代表不同處理下95%的置信區(qū)間;CK: 對(duì)照；CN25: 林冠施N 25 kg hm-2 a-1；CN50: 林冠施N 50 kg hm-2 a-1；UN25: 林下施N 25 kg hm-2 a-1；UN50: 林下施N 50 kg hm-2 a-1

2.4 AM真菌alpha多樣性指數(shù)和群落組成的變化與土壤因子變化之間的關(guān)系

在目前的時(shí)間尺度下,AM真菌alpha多樣性指數(shù)與土壤因子之間沒有顯著的相關(guān)性(表5)。在2013年,AM真菌群落組成的變化與土壤pH呈顯著正相關(guān)關(guān)系(表6),與其他土壤因子之間無顯著相關(guān)關(guān)系,而從2014年到2016年,所有土壤因子與AM真菌群落組成之間的相關(guān)性均未達(dá)到顯著水平(表6)。

表5 所有樣品AM真菌alpha多樣性指數(shù)與土壤因子之間的Pearson 相關(guān)性分析

3 討論

3.1 氮沉降對(duì)土壤理化性質(zhì)的影響

之前的meta分析表明,不論是在區(qū)域尺度上[52],還是在全球范圍內(nèi)[53],溫帶森林生態(tài)系統(tǒng)土壤pH值對(duì)N沉降的響應(yīng)都十分敏感。對(duì)于我國森林生態(tài)系統(tǒng)來說,N沉降速率<30 kg hm-2a-1也會(huì)顯著降低土壤pH值[52]。我們的研究也支持此觀點(diǎn),因?yàn)樵?013年和2016年,低N處理(25 kg hm-2a-1)與對(duì)照處理相比,pH值有顯著下降。同時(shí),高N(50 kg hm-2a-1)處理與對(duì)照相比,pH值也有顯著下降,但與低N處理(25 kg hm-2a-1)相比,并無顯著差異,這可能是由于本研究區(qū)域土壤pH背景值較低造成的。因?yàn)門ian 等[53]meta分析表明,土壤pH值的背景值越低,對(duì)N沉降的響應(yīng)越不敏感。而本研究區(qū)域的初始pH值為弱酸性(5.0—6.0),所以盡管施N速率增加一倍,土壤pH值并沒有很大的改變。

在2015年,林下施N土壤pH值與林冠施N處理相比顯著降低,這與Huang等[40]研究結(jié)果一致。這一現(xiàn)象說明在特定條件下,林下施N會(huì)高估自然N沉降對(duì)土壤環(huán)境的影響。

3.2 氮沉降對(duì)AM真菌alpha多樣性的影響

在目前的時(shí)間尺度下(4 a),施N方式、施N速率對(duì)AM真菌物種alpha多樣性指數(shù)都沒有顯著影響,且二者也沒有顯著交互作用,說明在目前的N素添加水平和時(shí)間尺度下,林下施N并沒有高估自然N沉降對(duì)AM真菌alpha多樣性的影響。盡管施N以后土壤pH值有顯著下降,但是pH值的變化并沒有對(duì)AM真菌alpha多樣性造成影響。土壤pH值與AM真菌的alpha指數(shù)沒有明顯關(guān)系(表5),這與之前該試驗(yàn)平臺(tái)的研究結(jié)果相符[32]。此外,之前的研究還發(fā)現(xiàn),土壤N/P比與AM真菌的豐富度和香儂指數(shù)關(guān)系密切[32],而在目前的時(shí)間尺度下,施N方式和施N速率都沒有明顯改變土壤N/P比(表1)。我們推測這可能是AM真菌alpha多樣性對(duì)N沉降沒有響應(yīng)的原因之一。隨著處理時(shí)間的延長,幾種試驗(yàn)處理下AM真菌alpha多樣性有趨同的趨勢(表2),這可能是森林生態(tài)系統(tǒng)土壤AM真菌對(duì)持續(xù)N素輸入產(chǎn)生了適應(yīng),也可能和N素添加濃度有關(guān)。例如,劉永俊等[54]研究發(fā)現(xiàn),低量施肥(30 g/m2(NH4)2HPO4)對(duì)垂穗披堿草(Elymusnutans)根系中AM 真菌物種豐富度沒有顯著影響,而當(dāng)施肥量達(dá)到為120 g/m2時(shí),施肥表現(xiàn)出極其顯著的負(fù)面作用。對(duì)于土壤真菌來說,N素添加帶來的影響還受到試驗(yàn)持續(xù)時(shí)間的影響[33],如Van Diepen等[25]研究發(fā)現(xiàn),在美國北部森林中,經(jīng)過連續(xù)12年的N素添加,區(qū)域A的AM真菌香儂指數(shù)有輕微的增加,而區(qū)域C有輕微的下降。本研究時(shí)間跨度為4 a,那么在更大的時(shí)間尺度上,N沉降是否會(huì)對(duì) AM真菌alpha多樣性產(chǎn)生顯著影響還值得進(jìn)一步研究。

3.3 氮沉降對(duì)AM真菌群落組成的影響

在試驗(yàn)早期(2013年),施N方式對(duì)AM真菌的群落組成有輕微影響,而且受到施N速率的制約：當(dāng)施N速率較低(25 kg hm-2a-1)時(shí),林下施N對(duì)AM真菌群落組成的影響大于林冠施N,支持了我們的假設(shè)；而當(dāng)施N速率較高(50 kg hm-2a-1)時(shí),兩種施N方式都沒有顯著改變AM真菌群落組成,與我們的假設(shè)不一致。這些結(jié)果說明在特定試驗(yàn)條件下,林冠施N比林下施N能更好地反應(yīng)自然N 沉降對(duì)AM真菌群落組成的影響,與Huang等[40]的研究結(jié)論相似。冠層對(duì)N的截留會(huì)影響到達(dá)地面的N的數(shù)量[39],截留的百分比和生態(tài)系統(tǒng)類型、N沉降量及研究區(qū)域土壤特征有關(guān),少則10%以下,多則70%—80%[41]。在本試驗(yàn)樣地內(nèi),林冠噴灑的N被冠層截留大概占到一半左右[55],在短時(shí)間內(nèi),林冠噴灑的N只有一半左右直接到達(dá)林下。所以林冠施N 25 kg hm-2a-1(CN25),大致相當(dāng)于林下施N 12.5 kg hm-2a-1,林冠施N速率為50 kg hm-2a-1(CN50),則大概相當(dāng)于林下施N 25 kg hm-2a-1(UN25),如此施N方式的問題就基本簡化為林下不同施N速率的問題。在Jiang等[56]的研究中,隨著施N速率的增加,N添加對(duì)AM真菌群落組成的影響也越大。在本研究中,CN25、CN50、UN25與對(duì)照相比,P值依次為0.115,0.080,0.005,基本上與Jiang等[56]的研究結(jié)果相一致,但是UN50與對(duì)照相比差異卻不顯著(P=0.480),這有可能是因?yàn)樵囼?yàn)處理重復(fù)次數(shù)較少,組內(nèi)的差異可能掩蓋了組間差異(圖2),所以在空間異質(zhì)性較大的環(huán)境條件下進(jìn)行控制試驗(yàn)時(shí),需要盡可能地增加重復(fù)次數(shù)。

從2014年到2016年,不同試驗(yàn)處理之間AM真菌群落組成差異不顯著,可能有兩方面的原因：(1)AM真菌對(duì)N素添加產(chǎn)生了一定程度的適應(yīng)；(2)試驗(yàn)處理時(shí)間還不夠長。Van Diepen等[25]研究還發(fā)現(xiàn),經(jīng)過連續(xù)12年的慢性N沉降,楓木 (Acerspp.)根中AM真菌的群落組成發(fā)生了顯著的改變,因此在更長的時(shí)間尺度上,施N是否會(huì)對(duì)土壤AM真菌群落組成產(chǎn)生影響還需要后續(xù)研究。

4 結(jié)論

本試驗(yàn)條件下,施氮對(duì)落葉闊葉混交林土壤AM真菌alpha多樣性沒有顯著影響,但在試驗(yàn)早期施氮對(duì)AM真菌的群落組成產(chǎn)生一定的影響,林下施氮高估了自然氮沉降對(duì)AM真菌群落組成的影響。隨著處理時(shí)間的延長,不同試驗(yàn)處理下AM真菌群落有趨同趨勢,這可能是因?yàn)锳M真菌群落對(duì)N沉降產(chǎn)生了適應(yīng)。本研究評(píng)估了不同施N方式對(duì)森林生態(tài)系統(tǒng)土壤AM真菌的影響,考慮到實(shí)驗(yàn)條件的局限性,在未來的研究中需要考慮更多的N素添加梯度和更長的時(shí)間跨度,才能夠更全面的認(rèn)識(shí)氮沉降的生態(tài)效應(yīng)。