模擬酸雨對我國亞熱帶毛竹林土壤呼吸及微生物多樣性的影響

王 楠, 潘小承, 白尚斌

1 東北林業大學林學院,哈爾濱 150040 2 浙江農林大學浙江省森林生態系統碳循環與固碳減排重點實驗室,杭州 311300

酸雨作為全球性環境問題之一,已受到國際社會普遍關注[1]。我國酸雨區面積仍不斷擴大,危害程度逐漸加重,已超過國土面積的40%,成為繼歐美之后的第三大酸雨區域,其中亞熱帶地區酸雨污染尤為嚴重[2]。浙江位于我國中亞熱帶,全省平均酸雨率超過80%,是酸雨發生典型區域[3]。酸雨給森林生態系統帶來嚴重影響和破壞,不僅對植物葉片造成直接損害,使葉片變黃甚至脫落,影響植物光合作用,而且可導致土壤酸化,造成鉀、鈣、鎂等鹽基離子淋溶,從而使土壤貧瘠,影響植物正常生長發育[4]。同時,土壤酸化使Al3+增加,引起鋁中毒,毒害植物根系[5]。土壤酸化還能促進H+濃度增加,引起土壤微生物群落組成和結構發生變化,不可避免地影響著土壤碳循環過程[6-7]。

土壤呼吸作為陸地生態系統碳循環的一個主要組分,在調控大氣CO2濃度和氣候變化方面起著關鍵作用[8]。土壤碳庫是大氣碳庫的4倍,土壤呼吸釋放碳的速率比人類活動釋放碳的速率要大一個數量級,其微小變化能極大改變大氣CO2濃度[9]。森林是陸地生態系統的主體,其碳庫占全球總碳庫的46%,占土壤呼吸總量的69%,因此,森林土壤呼吸成為陸地生態系統土壤呼吸的重要組分,其動態變化也對全球碳平衡影響深遠[10]。土壤呼吸主要包括植物根系呼吸和土壤微生物分解兩個過程,呼吸速率受植被特征、氣候、水熱條件、環境變化等諸多因素影響。酸雨可通過影響植物根系生長及影響微生物組成和活性變化來改變森林土壤呼吸速率[7]。一些研究認為酸雨促進了土壤呼吸[11],也有研究認為酸雨會抑制或不影響土壤呼吸速率[12]。目前關于酸雨影響森林土壤呼吸還未有明確定論,有待進一步深入研究。

竹林是世界重要的森林類型之一,蘊含著巨大的固碳增匯潛力。我國竹子有39屬,500余種,面積約601萬公頃,是竹子資源最豐富的國家,其中毛竹林(Phyllostachysedulis)面積約占73.8%。毛竹具有較高的社會和經濟效益,是應對全球變化的重要樹種,已得到廣泛推廣種植[13-14]。影響毛竹林土壤呼吸的環境因子諸多,其中酸雨脅迫是一個主要因素,然而關于酸雨對毛竹林土壤呼吸的影響及其微生物學響應機制研究卻相對較少[15]。因此,本文以我國熱帶區域分布較多的毛竹林為研究對象,采用PCR-DGGE技術,探討酸雨對毛竹林土壤細菌和真菌群落結構及多樣性的影響,及其與土壤呼吸的關聯,揭示土壤碳釋放模式的變化,從而為進一步研究亞熱帶毛竹林土壤碳釋放對環境問題的響應提供理論基礎。

1 研究地區與研究方法

1.1 研究地區概況

研究區域位于浙江省杭州臨安天目山國家級自然保護區(30°18′30″—30°21′37″N,119°24′11″—119°27′11″ E),保護區海拔300—1500 m,屬北亞熱帶季風氣候,四季分明,溫暖濕潤,年均氣溫14.5℃,年年平均降雨量1400 mm。其森林植被大體分為6種類型：常綠闊葉林,常綠落葉闊葉混交林,落葉闊葉林、針闊混交林,針葉林、毛竹林。毛竹林主要分布在海拔350—950 m。林冠層樹種主要有青岡(Cyclobalanopsisglance)、木荷(Schimasuperba)、苦櫧(Castanopsissclerophylla)、毛竹(Phyllostachyspubscens)等；灌木層有毛花連蕊茶(Camelliafraternal)、山礬(Symplocoscaudate)和山胡椒(Linderglauce)等；草本層有菊科(Compositae)和禾本科(Gramineae)等植物[16]。

1.2 試驗樣地設計

實驗樣地位于海拔485 m的毛竹林內,共設置3個樣地,在每個實驗樣地分別設置2個模擬酸雨小區和1個對照小區,在處理1(T1)小區噴灑pH 4.0的模擬酸雨；在處理2(T2)小區中噴灑pH 2.0的模擬酸雨；在對照(CK)小區噴灑pH值約為5.8的當地湖水。3個樣地呈三角形分布,間隔距離100 m,同一樣地的各處理在同一水平地上,緩沖間隔5 m,每個小區面積為1 m×1 m,3個樣地共計9個小區。參照臨安市以往的常規酸雨監測資料[6],按照H2SO4∶HNO3=8∶1的摩爾比配置模擬酸雨母液,再用母液與去離子水配成相應pH值的模擬酸雨。自2016年9月開始模擬酸雨處理,噴淋頻率為每兩周一次(如遇雨天或土壤濕度過大情況則順延,另擇時間進行噴淋,以防止噴灑的酸雨被地表徑流輸出)。每個小區每次噴灑酸雨量為1 mm,對照小區等量噴灑pH 5.8的湖水。

1.3 土壤呼吸速率及溫度測定

從2016年9月開始,進行了為期2年的土壤呼吸觀測試驗。預先把直徑20 cm、高10 cm的土壤呼吸底座埋入每個小區,將PVC環的一端沿坡面壓入土中約5 cm深,盡量減少布置PVC環對土壤的鎮壓作用,基座埋設過程中盡量不破壞土壤,以減少土壤擾動及根系損傷對測量結果的影響,PVC環在整個測量期間位置不變[17]。試驗期間,定期去除掉底座內的植物,以避免測定的土壤呼吸包含植株呼吸作用。

于每個觀測日8:00—12:00采用Li-8100土壤碳通量觀測系統(Li-COR,Lincoln,Nebraska,USA)測定土壤呼吸,測定頻率為每月2次,每次每個PVC環重復測定3次,取其平均值代表該次測定值,土壤呼吸測定在每次模擬酸雨噴灑前進行。在測定土壤呼吸的同時,采用Li-8100土壤碳通量觀測系統自帶的溫度探頭測5 cm深處的土壤溫度。

1.4 土壤生化指標的測定

于2018年8月(模擬酸雨二年后)進行土壤采集。分別在每個小區的四角與中心布設5個取樣點,采用無菌土鉆取表層大約5—20 cm深度的土壤樣品,更換樣地取樣前對工具采取清洗殺菌措施。每個取樣點中大約取出0.5 kg的土壤樣品,然后將每個土樣分為2份,一份土樣裝入無菌自封樣品袋中并置于4℃電子恒溫箱中,9個小區中共獲取45份土樣,用于測定土壤化學性質。將同一小區剩余5份土樣充分混合,過2 mm篩網,去除雜草與石塊后裝入無菌自封樣品袋中,并儲存與-80℃的干冰保溫箱中,用于提取DNA,共獲取9個土壤微生物樣品。

土壤pH值采用1∶2.5土水比,電位法測定；土壤溶解性有機碳(Soil dissolved organic carbon, DOC)使用高純度水浸提,取上清液過0.45m濾膜后使用島津TOC-VCPH 分析儀測定；堿解氮(Available nitrogen, AN)采用堿解蒸餾法測定；有效磷(Available phosphorous, AP)采用鉬銻抗分光光度法測定；速效鉀(Available potassium, AK)采用火焰原子吸收法測定；土壤微生物生物量碳(Soil microbial biomass of carbon, MBC)采用氯仿熏蒸-K2SO4浸提法測定[18]。

1.5 土壤細菌和真菌的變性梯度電泳(DGGE)

采用DGGE方法對土壤細菌和真菌進行檢測,細菌采用16S rDNA-PCR擴增所用一對通用引物338F-GC和518R,擴增片段長約250 bp；真菌18S rDNA-PCR擴增所用一對通用引物FR1-GC和FF390,擴增片段長約390 bp。20 μL反應體系如下：10×Premix Tap 10 L,引物(10 μmol/L)各0.2 μL,模板DNA 0.3 μL,牛血清蛋白(BSA,20 mg/mL)0.2 μL。細菌的PCR反應參數：94℃預變性5 min,94℃變性1 min,50℃退火1 min,72℃延伸1 min,30個循環,最后72℃延伸10 min。真菌的PCR反應參數：95℃預變性8 min,94℃變性30 s,50℃退火45 s,72℃延伸2 min,30個循環,最后72℃延伸10 min。取3 μL PCR產物經1.0%瓊脂糖凝膠電泳檢測產物及其濃度[19]。

采用The DCodeTMUniversal Mutation Detection System(Bio-Rad),細菌擴增產物在變性梯度為45%—60%的6%聚丙烯酰胺凝膠上進行電泳。真菌擴增產物在變性梯度為25%—40%的6% 聚丙烯酰胺凝膠上進行電泳。細菌和真菌的電泳過程條件一致：在60℃、80 V條件下電泳13 h,使用SYBR greenⅠ染色0.5 h,染色結果于Gel Doc TM EQ(Bio-Rad)凝膠成像系統成像,使用Quantity One 4.4軟件(Bio-Rad)進行圖像分析。

1.6 數據分析

土壤呼吸速率與土壤5 cm深溫度的關系采用如下指數模型：

R=αebT

式中,R為土壤呼吸速率；T為土壤5 cm深溫度；α是溫度為0℃時土壤呼吸速率；b為溫度反應系數。

Q10值的計算公式[20]：

Q10=e10b

Shannon多樣性指數(H)、Margalef豐富度指數(D)以及Pielou均勻度指數(E)對土壤細菌微生物群落多樣性進行分析[21],其計算公式為：

H= -∑(ni/N) ln(ni/N)

D= (S-1)/lnN

E=H/lnS

式中,ni為每條電泳條帶光密度峰值,N為同一泳道中所有條帶光密度峰值總和。電泳條帶光密度的峰值通過分析軟件Quantity one進行讀取。

采用Origin Pro8.0軟件對數據進行處理和繪圖,采用Pearson法進行相關性分析,采用t檢驗進行不同處理土壤呼吸速率的均值比較,采用單因素方差分析(one-way ANOVA)和最小顯著差異法(LSD)對3種處理間土壤呼吸速率進行方差分析和多重比較,用Curve Estimation對土壤呼吸速率與土壤溫度關系進行回歸分析,R語言vegan包進行冗余分析(Redundancy Analysis, RDA)和作圖。

2 結果與分析

2.1 不同酸雨處理土壤呼吸的季節變化

土壤呼吸速率具有明顯的季節動態規律,最高值出現在7—8月,最低值出現在12月份至次年2月期間。如圖1所示,觀測期可分為5個階段：階段1(2016年9月—2017年2月)模擬酸雨對土壤呼吸速率的影響沒有顯著差異；階段2(2017年3月—2017年6月)模擬酸雨對土壤呼吸速率產生了促進作用,模擬強酸降雨樣地與對照之間具有顯著差異(P<0.05)；階段3 (2017年7月—2017年10月)模擬強酸降雨與對照樣地、模擬弱酸降雨樣地之間具有顯著差異(P<0.05),而對照樣地與模擬弱酸降雨樣地之間的差異不顯著,階段4(2017年11月—2018年2月)模擬酸雨對土壤呼吸速率的影響沒有顯著差異；階段5(2018年3月—2018年8月)模擬酸雨對土壤呼吸速率產生了抑制作用,且模擬強酸降雨樣地與對照之間具有及顯著差異(P<0.05)。

2.2 不同酸雨處理土壤呼吸與土壤溫度的關系

根據對土壤溫度檢測數據分析(圖2),土壤呼吸速率與土壤溫度的季節變化趨勢相近,呈顯著正相關關系,方差分析表明不同模擬酸雨處理的土壤溫度無顯著差異,年均土壤溫度為16.6℃。

回歸分析表明,CK、T1、T2處理的土壤呼吸與土壤溫度之間呈極顯著的指數回歸關系(表1),這3個小區的回歸方程的決定系數R2分別為 0.960、0.937、0.845,這表明土壤溫度可分別解釋這 3個小區土壤呼吸96.0%、93.7%、84.5%的季節變異。根據指數回歸方程估算得到的 CK、T1、T2處理的土壤呼吸的溫度敏感系數(Q10)分別為2.14、2.10、2.86。T2與T1和CK間差異顯著,T1與CK間無顯著性差異,這一結果表明,重度酸雨處理增加了溫度敏感性。

圖1 不同酸雨處理下毛竹林土壤呼吸動態變化Fig.1 Dynamic change of soil respiration in moso bamboo under different acid rain treatments during 2016—2018CK：對照處理, control treatment；T1：處理一,treatment 1；T2： 處理二,treatment 2

表1 不同pH值酸雨處理的毛竹林土壤呼吸速率與溫度相關性

不同字母表示同一列內各均值存在顯著差異(P<0.05)；*表示顯著相關(P< 0.05),**表示極顯著相關(P < 0.01); CK：對照處理, control treatment；T1：處理一,treatment 1；T2： 處理二,treatment 2

2.3 不同酸雨處理土壤理化性質的差異

模擬酸雨對土壤pH值、養分及碳、氮含量的影響差異顯著(表2)。T2處理土壤的平均pH值比CK處理低17.31%,比T1處理低14.00%；T2處理土壤有機質含量比CK低47.55%,比T1處理低30.67%；CK處理土壤DOC、AK、MBC的含量分別是T1處理的1.22、1.16、1.01倍,分別是T2處理的1.39、1.36、1. 06。T2處理土壤AN、AP的含量分別是CK處理的1.09、1.67倍,T1的1.05、1.62倍。由此可見,不同強度的酸雨對土壤理化性質和微生物量具有顯著的影響。

圖2 土壤呼吸與土壤溫度的關系Fig.2 Relationship between soil respiration and soil temperature

表2 不同酸雨處理下土壤理化性質和微生物量

不同字母表示同一列內各均值存在顯著差異(P<0.05); CK：對照處理, control treatment；T1：處理一,treatment 1；T2： 處理二,treatment 2

2.4 不同酸雨處理土壤細菌和土壤真菌多樣性

模擬酸雨實驗2年后,對不同處理的土壤微生物群落采用DGGE技術檢測,確定了共30條DGGE帶。結果表明,模擬酸雨處理的所有樣品中均含有8個DGGE譜帶(3、7、11、13、14、23、28和30)。在這些條帶中,3、30條帶被鑒定為細菌,7、11、13、14、23、28條帶被分別鑒定為放線菌、分枝桿菌、硝化螺、酸性微生物和鞘氨醇菌。樣品中含有一些共有菌群,如細菌、放線菌、厚壁菌門、硝化螺旋菌門和α-變形菌。酸雨也造成了某些菌群的減少或消失,如帶4(莖菌屬)、9(孢囊桿菌屬)、12(慢生根瘤菌屬)和25(念珠菌屬),同時,酸性條件也導致某些種類的增加,如帶5(肉食桿菌屬)、6(丙酸桿菌屬)、8(梭形芽孢桿菌屬)。這說明在酸雨作用下這類微生物對土壤有機質降解和釋放土壤二氧化碳中起著更大作用。

不同模擬酸雨處理的土壤細菌和土壤真菌多樣性具有顯著差異(表3),模擬酸雨處理的土壤細菌群落多樣性表現為降低,說明酸性脅迫抑制了土壤細菌群落的多樣性。T2處理的土壤細菌多樣性與CK處理相比存在顯著差異(P<0.05),說明重度酸雨處理下土壤細菌群落對酸性環境的耐受性最低。而模擬酸雨處理的土壤真菌群落多樣性表現相反,T1處理的土壤真菌多樣性與CK處理和T2處理存在顯著差異(P<0.05),說明中度酸雨處理條件下土壤真菌群落對酸性環境的耐受性最高,弱酸可促進土壤真菌群落的多樣性。

表3 不同酸雨處理下土壤微生物群落的多樣性指數

不同字母表示同一列內各均值存在顯著差異(P<0.05)

2.5 不同酸雨處理土壤性質與土壤微生物的關系

土壤性質與土壤細菌、土壤真菌α多樣性指數的相關性分析表明(表4),土壤細菌的Shannon多樣性指數、Margalef豐富度指數與堿解氮和有效磷呈顯著或極顯著負相關(P<0.01),與pH值、速效鉀、可溶性有機碳和微生物量碳呈顯著或極顯著正相關(P<0.01)；土壤真菌的Shannon多樣性指數與堿解氮顯著正相關(P<0.05),與可溶性有機碳和微生物量碳呈顯著或極顯著負相關(P<0.01),Margalef豐富度指數與pH值、可溶性有機碳和微生物量碳呈顯著負相關(P<0.05),與堿解氮為顯著正相關(P<0.05)。

表4 土壤微生物α多樣性指數與土壤理化性質、微生物量的相關系數

*表示顯著相關(P< 0.05),**表示極顯著相關(P< 0.01)

圖3 土壤微生物群落與環境影響因子的冗余分析Fig.3 RDA of soil microbial community and environmental factors pH：酸堿度；AN：堿解氮,Available N；AP：有效磷,Available P；AK：速效鉀,Available K；DOC：可溶性有機碳,Dissolved organic carbon；MBC：微生物量碳,Microbial biomass carbon

為了解釋酸雨與土壤微生物群落結構及環境因子之間的關系,對不同模擬酸雨樣地土壤微生物群落與環境因子進行冗余分析,可以看出CK與T1、T2之間具有明顯地距離,這證實了之前的分析,即不同模擬酸雨處理樣地的土壤微生物群落在差異很大程度上與環境變量有關。兩軸累計貢獻率為97.74%,T1和T2分布在第二排序軸正方向,CK分布在第二排序軸負方向,表明不同模擬酸雨處理是第二排序軸的主要影響因素。pH、AK、DOC和MBC與模擬酸雨處理的土壤微生物群落結構關系緊密且正相關,AN、AP與土壤微生物群落結構負相關。排列檢驗顯示 pH、DOC、AP、AN和MBC對土壤細菌群落結構的影響最顯著(P<0.05)(圖3)。

2.6 不同酸雨處理土壤呼吸與土壤性質、土壤微生物的關系

土壤呼吸與土壤理化性質、土壤微生物α多樣性指數的相關性分析表明(表4),土壤呼吸與pH值、速效鉀、可溶性有機碳、微生物量碳和細菌Shannon指數呈顯著正相關(P<0.05),與堿解氮和有效磷為顯著負相關(P<0.05)。

3 討論

3.1 不同酸雨處理土壤呼吸的差異性及控制因子

由于各地區的氣候、土壤養分及植被類型等條件的不同,模擬酸雨對土壤呼吸的影響結果存在較大差異。Reth等[22]研究發現模擬酸雨抑制了土壤呼吸,這種抑制作用可能是氫離子高載荷對土壤微生物產生毒性所致。Fangueiro等[23]報道,低模擬酸雨強度反而促進了土壤呼吸。本研究結果表明,實驗初期模擬酸雨對土壤呼吸無明顯影響,這主要是由于土壤擁有巨大的緩沖體系,對酸性物質的輸入具有很強的緩沖作用,毛竹林根系短期內對酸雨的耐受性較強。隨著時間推進,在模擬酸雨作用下,土壤呼吸呈現出先增大后減小的趨勢,酸雨的促進土壤呼吸作用可能是根細胞中維持蛋白質周轉和溶質梯度所消耗的能量占根系維持呼吸的絕大部分,如果根系組織中的氮含量保持不變,那么根系呼吸也不會發生變化,而酸雨中的NO3-的施肥作用恰恰增加了毛竹林根系細胞的氮含量,使得地下的土壤呼吸速率加強[12]。模擬酸雨下森林的土壤酸化是個逐漸累積的過程,隨著酸雨強度的增大和時間的推移,模擬酸雨對土壤呼吸產生抑制作用,這與其脅迫下土壤酸化從而導致土壤微生物活性及土壤養分流失有關。相關性分析表明(表4),土壤呼吸與土壤細菌α多樣性指數呈顯著正相關(P<0.05),土壤微生物活性下降,從而延緩了土壤有機物的礦化和分解速率,抑制了土壤CO2的釋放,并呈現逐漸增強的趨勢[24]。

表4 土壤呼吸與土壤微生物α多樣性指數、土壤理化性質的相關系數

*表示顯著相關(P< 0.05),**表示極顯著相關(P< 0.01)

酸雨對土壤微生物量碳的影響與酸雨脅迫下土壤持續酸化有關,通常pH值較低的土壤中微生物量碳含量相應較低[25]。本研究在模擬酸雨的持續作用下,土壤微生物量碳的含量及土壤呼吸速率呈顯著下降特征,這是由于酸雨導致土壤 pH改變,引起土壤微生物活性、功能及土壤微生物量碳的變化。該結果與王國兵等[26]研究結論一致,土壤呼吸產生的CO2部分是由于微生物的自養呼吸和分解有機質和植物凋落物的異養呼吸而產生的,微生物生物量碳的減少間接地導致了土壤呼吸速率的減少[27]。另外,酸雨使土壤pH值降低,有機質礦化速率下降,土壤腐殖質層DOC的釋放量減少[28]。本研究中毛竹林土壤溶解性有機碳呈隨著酸雨強度增加而下降的趨勢,可能由于DOC來源于凋落物的分解以及微生物代謝,微生物活減少,產生的水溶性有機物質也減少[29-30]。因此,酸雨對土壤微生物活性的影響被認為是影響DOC釋放的關鍵因素。

土壤溫度是影響土壤呼吸季節性動態變化的最主要的因素,土壤呼吸為酶促反應,而土壤溫度是影響酶促反應的重要因素,溫度升高往往會導致酶促反應加快,從而促進土壤呼吸。本研究表明,酸雨各處理的土壤呼吸季節變化具有明顯規律性,在不受土壤水分脅迫的條件下酸雨各處理的土壤呼吸與土壤溫度之間遵循指數回歸關系,這與Nikolova等[31]以及Karhu等[32]的研究結果一致。梁國華等[33]發現酸雨降低了土壤呼吸的溫度敏感性,并且酸性越強抑制作用越顯著。Chen等[34]研究表明酸雨對土壤呼吸的溫度敏感性沒有影響。本研究發現模擬酸雨改變了土壤呼吸對溫度的響應,表現為Q10值在酸處理下有上升的趨勢,說明模擬酸雨增加了土壤呼吸的溫度敏感性,這可能是因為DOC隨酸處理的增加而減少,不同生態系統土壤呼吸季節變異的Q10與DOC呈負相關[35]。毛竹林生態系統中毛竹林凋落物分解因酸雨而促進,凋落物有機碳的質量逐漸下降,土壤呼吸的Q10值越來越高[36]。Q10值是在整個測定期間進行計算的,分高溫階段或低溫不同階段分析可能有更為具體直觀的結果,今后可進一步研究酸雨脅迫下土壤呼吸在不同溫度時段的溫度敏感性,以及進一步測定土壤呼吸的日變化也是今后研究值得考慮的方面。

3.2 不同模擬酸雨條件下土壤微生物群落多樣性的比較及控制因子

土壤呼吸與土壤有機物的轉化與循環密切相關,而土壤有機物的轉化與循環有賴于微生物活化才能完成。酸化的土壤由于H+的毒害作用,從而改變了土壤分解者微生物的種類,結構以及生物活性[37]。不同pH的模擬酸雨均能抑制細菌、放線菌、氮素生理類群的生長,其數量隨著酸雨pH的降低而不斷減少,真菌、纖維素分解菌的數量隨酸雨pH的降低呈現先升高后下降的趨勢[38]。本研究結果表明,模擬酸雨對毛竹林土壤細菌和真菌群落多樣性、豐富度具有顯著影響,通過比較指數值的大小可以發現,酸雨能夠在一定程度上降低土壤細菌群落多樣性,而弱酸能夠顯著提高土壤細菌群落多樣性,同時對酸雨具有抗性或耐性的微生物種群會表現出更高的活性。

有關不同地區土壤微生物多樣性的研究表明,土壤pH值是影響土壤微生物多樣性的重要因子[39]。Pearson相關性分析表明,土壤細菌的多樣性指數和豐富度指數與pH值呈顯著或極顯著正相關,這與孫巖和全炳武[40]的研究結果一致,這可能是酸雨脅迫會抑制大多數細菌的生長和繁殖,使土壤細菌的多樣性和豐富度下降。同時,我們利用冗余分析可以揭示不同酸雨處理條件下土壤性質對微生物結構的影響,發現其結果與土壤細菌相近,這也說明土壤細菌群落是毛竹林土壤環境中的主導微生物。本研究中,土壤微生物與pH值、有效鉀、可溶性有機碳和微生物量碳呈顯著正相關,同時表現出與土壤呼吸呈正相關關系。這可能是與南方酸性土壤有關,酸雨淋溶增強,土壤的分解與礦化過程較快,土壤元素的遷移較強,尤其是鉀的流失量較大,土壤含鉀量低,抑制土壤微生物的生長。土壤微生物與土壤堿解氮和有效磷含量呈顯著負相關關系,可能原因pH值降低,土壤對氮、磷酸根的吸附量降低。酸雨對土壤生態系統的影響是個長期累積的過程,開展酸雨脅迫下毛竹林土壤呼吸的長期定位觀測具有較高研究價值,盡管本次研究經歷兩年的實驗,但酸雨處理區域面積偏小,今后可增大實驗區面積,以更好地深入研究酸雨脅迫引起的土壤酸化對土壤生態系統過程影響,這些工作開展將有利于我國毛竹林的科學經營及亞熱帶常綠闊葉林的生態保護。

4 結論

(1)模擬酸雨對土壤呼吸速率的影響表現為先促進后抑制趨勢,土壤呼吸速率呈季節性規律變化,不同模擬酸雨處理改變了土壤呼吸對溫度的敏感性。

(2)模擬酸雨改變土壤微生物群落結構及多樣性,重度酸雨處理抑制了土壤細菌群落生長,而中度酸雨對土壤真菌多樣性和豐富度具有促進作用。

(3)pH值、速效鉀、可溶有機碳等土壤理化性質與土壤微生物群落結構、土壤呼吸密切相關。