土壤微生物群落對麻櫟-刺槐混交林凋落物分解的影響

董學德,高 鵬,*,李 騰,張佳辰,董金偉 ,許景偉,囤興建

1 山東農業大學林學院,泰山森林生態系統國家定位觀測研究站, 泰安 271018 2 山東省林業科學研究院, 濟南 250014

森林凋落物分解是一個由生物因素和非生物因素引起的物理、化學性質的變化[1-2]。影響凋落物分解的因素主要包括凋落物質量、氣候因素等非生物因素和生物因素,其中,生物因素是影響凋落物分解的主導因子,非生物因素一部分是通過影響生物因素發揮作用[3]。土壤微生物作為連接植物和土壤的重要生物因子,不僅能通過產生各種酶水解和氧化凋落物,而且能夠進行氨化、硝化、固氮等過程,改變凋落物的結構和化學組成[4]。土壤微生物多樣性是指生命體在遺傳、種類和生態系統層次上的變化,代表著微生物群落的穩定性。在凋落物分解過程中,土壤微生物群落會對外界環境的變化迅速作出響應,微生物多樣性也會發生變化,進而影響凋落物的分解速率[5]。因此,關于土壤微生物群落多樣性對森林凋落物分解過程的影響機制研究受到人們普遍關注[1,6]。

近些年,許多學者針對凋落物質量和環境因素對凋落物分解速率的影響進行了較為全面的研究[7-8]。但是,土壤微生物作為影響凋落物分解的關鍵因子尚缺少系統的研究,混交林土壤微生物多樣性與凋落物分解速率的關系尚不明確。山東泰安市黃前流域是國家級水土流失重點治理區,麻櫟林及其刺槐混交林是該區域重要的生態造林類型。開展麻櫟-刺槐混交林凋落物分解過程中土壤微生物群落多樣性變化特征及其對分解速率的影響的研究,可為該區域生態修復樹種選擇和綜合效果評價提供重要理論依據。因此,以麻櫟-刺槐混交林和麻櫟純林為研究對象,采用野外定點采樣、室內分析和高通量測序相結合的方法,對麻櫟-刺槐混交林和麻櫟純林凋落物分解過程中,土壤細菌、真菌群落多樣性特征及其對凋落物分解速率的影響進行研究,有助于深入了解混交林凋落物分解速率與土壤微生物群落的關系。

1 材料與方法

1.1 研究區概況

研究區位于山東泰安市黃前流域(117°04′—117°22′E,36°17′—36°27′N),屬于黃河下游大汶河流域,是國家級水土流失重點治理區,總面積292.2 km2,海拔310—950 m,平均海拔530 m。屬暖溫帶大陸性半濕潤季風氣候,年平均氣溫為18.5℃,年平均降水量758 mm,75%的降水集中在6—9月,無霜期197 d,≥10℃的活動積溫3819℃,多年平均日照時數2640 h,土壤類型為棕壤(表1)。研究區位于國家級重點生態防護林區,森林覆蓋率為96%,植被類型為暖溫帶落葉闊葉林,主要樹種為麻櫟(QuercusacutissimaCarruth.)、刺槐(RobiniapseudoacaciaL.)、赤松(PinusdensifloraSieb. et Zucc.)、日本落葉松(Larixkaempferi(Lamb.)Carr)等。

表1 研究區標準地基本狀況

1.2 研究方法

1.2.1凋落物采樣

于2019年4月,選擇立地條件相近的麻櫟-刺槐林混交林和麻櫟純林林地,在每個林地隨機設置3個20 m×20 m的標準樣地,并在每個標準樣地隨機設置3個2 m×2 m采樣樣方。采樣時為保證凋落物的新鮮程度一致,只收集樣方內現存在的尚未分解的完整凋落葉和樹枝(樹枝數量極少且為細枝)作為試驗樣品。將試驗樣品帶回實驗室,剔除其中的泥土等雜質,并將樹枝剪成5 cm左右,置于室內風干,然后在85℃下烘干至恒重。將凋落葉和樹枝混合均勻后稱取16 g(精確到0.001 g)烘干的原狀樣品,裝入帶有編號的分解袋(孔徑0.2 mm,規格為16 cm×26 cm)中。取上述麻櫟-刺槐和麻櫟凋落物,測定凋落物初始TOC、TN、TP濃度。2019年5月將麻櫟-刺槐混合凋落物(混合比例為：麻櫟∶刺槐=12∶4)置于麻櫟-刺槐混交林采樣樣方內,麻櫟凋落物置于麻櫟純林采樣樣方內,用竹簽固定,讓其在自然條件下分解。在凋落物分解30、60、90、120、150、180、210、240 d時,在各個采樣樣方內隨機取3個分解袋,帶回實驗室,將凋落物烘干、稱重,測定分解速率；然后,用粉碎機粉碎,過0.149 mm篩子,測定TOC、TN、TP濃度。

1.2.2土壤微生物測定樣品采樣

土壤微生物測定樣品的采樣時間為凋落物初次放置及分解120、240 d時,采樣點位于凋落物采樣樣方內。在每個采樣樣方內,采用五點取樣法取0—20 cm的土壤混合均勻后裝于試管,立即放于液氮內保存。每次采樣麻櫟-刺槐混交林和麻櫟純林各得9個(每種林分包括3塊標準樣地,每塊標準樣地包括3個采樣區域)土壤樣本,用于測定土壤細菌和真菌群落結構。

1.2.3樣品測定

凋落物：TOC濃度測定采用K2Cr2O7容量法—外加熱法；TN濃度測定采用半微量開氏法消煮法；TP濃度測定采用鉬銻抗比色法測定[9]。

土壤微生物：采用CTAB方法對樣本的基因組DNA進行提取,之后利用瓊脂糖凝膠電泳檢測DNA的純度和濃度。土壤細菌多樣性的檢測：引物擴增區域為16S V3-V4；引物名稱及序列：341F：CCTACGGGRBGCASCAG；806R：GGACTACNNGGGTATCTAAT[10]。土壤真菌多樣性檢測：引物擴增區域為ITS1；引物名稱及序列：ITS5-1737F：GGAAGTAAAAGTCGTAACAAGG；ITS2-2043R：GCTGCGTTCTTCATCGA TGC[11]。PCR擴增體系：2*taq PCR mix：25.0 μL,Primer F (10 μM)：1.0 μL,Primer FR(10 μM)：1.0 μL,gDNA：2.5 μL,H2O：8.0 μL。PCR反應程序為：95℃預變性5 min,34個循環 (94℃,1 min; 57℃,45 s；72℃,1 min)；72℃,10 min；16℃,5 min[10]。使用New England Biolabs 公司的Phusion? High-Fidelity PCR Master Mix with GC Buffer和高效高保真酶進行PCR。產物純化試劑盒使用的是Thermo Scientific 公司GeneJET 膠回收試劑盒。使用Thermofisher 公司的Ion Plus Fragment Library Kit 48 rxns 建庫試劑盒進行文庫的構建,構建好的文庫經過Qubit 定量和文庫檢測合格后,使用Thermofisher的Ion S5TMXL進行上機測序。

1.3 數據處理與分析

凋落物質量殘留率：凋落物分解一段時間后剩余的質量與凋落物初始質量比值[12],即：

RM=Mt/M0×100%

式中,RM為凋落物質量殘留率(%)；Mt為凋落物分解t時刻的質量(g)；M0為凋落物的初始質量(g)。

凋落物分解的Olson模型：凋落物分解是一個動態過程,凋落物質量殘留量與分解時間存在負指數函數關系[13],即：

Mt/M0=ae-kt

式中,a為修正參數；k為凋落物分解速率；t為凋落物分解時間(a)。凋落物分解50%(t50%)和分解95%(t95%)所需的時間分別為：

t50%=-ln0.5/k

t95%=-ln0.05/k

凋落物元素釋放率：凋落物分解一段時間內釋放的元素量與凋落物初始元素含量比值的百分數,即：

E=[(E0-Et)/E0]×100%

式中,E為凋落物元素釋放率(%)；Et為凋落物分解t時刻凋落物質量與元素濃度的乘積(g)；E0為凋落物的初始質量與元素濃度的乘積(g)。當百分數大于0時,說明凋落物中的元素發生了凈釋放；當百分數小于0時,說明凋落物中的元素發生了凈積累[14]。

土壤微生物Alpha多樣性：微生物菌群Alpha多樣性可以分為菌群豐富度指數(Community richness index)和菌群多樣性指數(Community diversity index)。菌群豐富度指數主要包括Chao1和ACE,用于評估群落中含有的微生物物種總數,菌群多樣性指數主要包括Simpson和Shannon,用于評估微生物的多樣性[15]。

使用Uparse軟件對所有樣品的全部有效數據在97%相似度水平下進行聚類,獲得Operational Taxonomic Units(OTU),統計各個土壤樣品每個OTU中的豐度信息。使用Qiime 1.9.1軟件計算兩種林分土壤細菌、真菌群落的Alpha多樣性指數值。使用R 3.5.1進行數據分析,麻櫟-刺槐混交林和麻櫟純林分解過程中土壤微生物優勢類群間的差異采用單因素方差分析；采用T檢驗分析微生物菌群豐富度指數和菌群多樣性指數值在麻櫟-刺槐混交林和麻櫟純林之間的差異；對麻櫟-刺槐混交林和麻櫟純林凋落物的分解速率與土壤微生物菌群豐富度指數和菌群多樣性指數值之間的相關性進行Spearman相關分析。其余相關圖的制作使用R 3.5.1的“ggplot2”包完成。

2 結果與分析

2.1 麻櫟-刺槐混交林和麻櫟純林凋落物分解動態特征

圖1 麻櫟-刺槐混交林和麻櫟純林凋落物分解過程中的質量損失量動態變化Fig.1 Dynamics changes of mass loss during litter decomposition of QR and QAC圖中數據為平均值±標準差(n=3)

在凋落物分解期間(2019年5月至2020年1月),麻櫟-刺槐混交林和麻櫟純林凋落物質量損失量均表現出先增加后減小的趨勢(圖1)。麻櫟-刺槐混交林凋落物質量損失量在0—90 d逐漸增大,此后逐漸降低,累積質量損失6.32 g；麻櫟純林凋落物質量損失在0—120 d逐漸增大,120 d后逐漸降低,累積質量損失6.21 g。麻櫟-刺槐混交林凋落物質量損失量在0—150 d高于麻櫟純林凋落物質量損失量,此后低于麻櫟純林。

麻櫟-刺槐混交林和麻櫟純林凋落物質量殘留率與分解時間符合Olson負指數模型(圖2)。麻櫟-刺槐混交林凋落物分解速率高于麻櫟純林的,麻櫟-刺槐混交林凋落物分解速率0.70,分解240 d后質量殘留率為60.2%,分解50%所需時間為0.99 a,分解95%所需時間為4.28a；麻櫟純林凋落物分解速率0.62,分解240 d后質量殘留率62.8%,分解50%所需時間為1.12 a,分解95%所需時間為4.83 a。

圖2 麻櫟-刺槐混交林和麻櫟純林凋落物殘留率與時間的Olson負指數回歸關系Fig.2 Olson negative index regression relation between litter residues and time of QR and QAC

在凋落物分解過程中,凋落物的TOC、TN、TP濃度在麻櫟-刺槐混交林和麻櫟純林之間存在差異(圖3)。麻櫟-刺槐混交林凋落物TOC濃度低于麻櫟純林的,麻櫟純林平均TOC濃度是麻櫟-刺槐混交林的1.2倍,并且兩種林分TOC濃度總體呈波動降低的趨勢(圖3)。麻櫟-刺槐混交林凋落物TN濃度高于麻櫟純林的,麻櫟-刺槐混交林凋落物平均TN濃度是麻櫟純林的1.2倍,兩種林分TN濃度總體呈波動升高趨勢(圖3)。麻櫟-刺槐混交林凋落物TP濃度高于麻櫟純林的,麻櫟-刺槐混交林凋落物TP濃度是麻櫟純林的1.1倍,兩種林分TP濃度在0—120 d逐漸升高,此后至240 d逐漸降低(圖3)。

圖3 麻櫟-刺槐混交林和麻櫟純林凋落物分解過程中TOC、TN、TP濃度變化特征Fig.3 Dynamics of TOC, TN and TP concentrations during litter decomposition of QR and QAC圖中數據為平均值±標準差(n=3)

在凋落物分解240 d后,麻櫟-刺槐混交林和麻櫟純林凋落物TOC、TN、TP均發生凈釋放,且麻櫟-刺槐混交林釋放率高于麻櫟純林(表2)。麻櫟-刺槐混交林和麻櫟純林凋落物TOC釋放率表現為逐漸增大的趨勢,在分解第240 d時釋放率達到最大值,分別為47.4%和46.3%。TN釋放率在麻櫟-刺槐混交林凋落物中表現為逐漸增加的趨勢,在240 d達到最大值,為26.3%；TN釋放率在麻櫟純林凋落物中表現為先增加后減小的趨勢,在120 d時達到最大值,為22.8%。TP在麻櫟-刺槐混交林和麻櫟純林凋落物中表現出先積累后釋放的趨勢,麻櫟-刺槐混交林發生積累的時間比麻櫟純林的長,并且均在240 d時達到最大值。

表2 麻櫟-刺槐混交林和麻櫟純林凋落物分解過程中TOC、TN、TP釋放率變化特征 (平均值±標準差)

2.2 麻櫟-刺槐混交林和麻櫟純林凋落物分解過程中土壤微生物群落變化特征

土壤微生物群落結構的變化主要是優勢類群的變化引起的[16],因此,本試驗只對土壤細菌和真菌優勢類群進行重點分析。在麻櫟-刺槐混交林和麻櫟純林凋落物的分解過程中,土壤細菌優勢類群為放線菌門(Acidobacteria)、變形菌門(Proteobacteria)、酸桿菌門(Actinobacteria)和疣微菌門(Verrucomicrobia)(圖4)。兩種林分四種細菌優勢類群相對含量之間存在顯著差異(P<0.05),放線菌門的相對含量顯著高于變形菌門、酸桿菌門和疣微菌門。土壤細菌優勢類群的相對含量在麻櫟-刺槐混交林和麻櫟純林之間存在顯著差異(P<0.05),麻櫟-刺槐混交林的放線菌門、變形菌門和疣微菌門的相對含量高于麻櫟純林,酸桿菌門相對含量低于麻櫟純林。土壤真菌優勢類群為擔子菌門(Basidiomycota)、子囊菌門(Ascomycota)和被孢霉門(Moritierellomycota)(圖4)。兩種林分3種真菌優勢類群相對含量之間存在顯著差異,擔子菌門的相對含量顯著高于子囊菌門和被孢霉門(P<0.05)。土壤真菌優勢類群的相對含量在麻櫟-刺槐混交林和麻櫟純林間存在顯著差異(P<0.05)。麻櫟-刺槐混交林中擔子菌門的相對含量高于麻櫟純林的,兩種林分土壤擔子菌門的相對含量為70.4%—85.1%和41.9%—64.2%；麻櫟-刺槐混交林土壤子囊菌門的相對含量顯著低于麻櫟純林的(P<0.05),兩種林分土壤子囊菌門的相對含量為6.6%—8.0%和24.5%—29.8%；麻櫟-刺槐混交林被孢霉門的相對含量顯著高于麻櫟純林,其土壤被孢霉門的相對含量為1.4%—2.8%和1.5%—3.9%。

圖4 麻櫟-刺槐混交林和麻櫟純林凋落物分解過程中土壤細菌和真菌優勢類群相對含量變化Fig.4 Relative content changes of soil bacterial and fungal dominant communities during litter decomposition of QR and QACQR：麻櫟-刺槐混交林Mixed Quercus acutissima Carruth. and Robinia pseudoacacia L. forest; QAC：麻櫟純林Pure Quercus acutissima Carruth forest; QR1、QR2、QR3分別代表麻櫟-刺槐混交林凋落物初次放置、分解120、240 d時土壤微生物取樣,QAC1、QAC2、QAC3分別代表麻櫟純林凋落物初次放置、分解120、240 d時土壤微生物取樣;圖中數據為樣本平均值(n=3)

在麻櫟-刺槐混交林凋落物分解過程中,微生物菌群豐富度指數和菌群多樣性指數發生顯著變化(表3)。細菌菌群豐富度指數(Chao1指數和ACE指數)呈現先增大后減小的趨勢,在分解120 d時達到最大值,在分解240 d時顯著降低。細菌菌群多樣性指數(Shannon指數和Simpson指數)均呈先增大后減小的趨勢,在分解240 d時顯著降低。真菌菌群豐富度指數呈先減小后增大的趨勢,在初次放置時達到最大值,此時顯著高于凋落物初分解120、240 d時,分解中期達到最小值。真菌菌群多樣性指數呈逐漸減小的趨勢,但在不同分解時間差異不顯著。

在麻櫟純林凋落物分解過程中,微生物菌群豐富度指數和菌群多樣性指數發生顯著變化(表3)。細菌菌群豐富度指數呈先增加后減小的趨勢,不同分解時間存在顯著差異,在分解120 d時達到最大值,分解240 d時達到最小值。細菌菌群多樣性指數呈先增大后減小的趨勢,在分解240 d時顯著降低。真菌豐富度指數呈逐漸減小的趨勢,凋落物初次放置時顯著高于分解120、240 d,在分解240 d時達到最小值。真菌多樣性指數呈先減小后增加的趨勢,凋落物初次放置時顯著高于分解120、240 d。

麻櫟-刺槐混交林土壤微生物菌群豐富度指數和菌群多樣性指數變化范圍小于麻櫟純林(表3)。麻櫟-刺槐混交林細菌Chao1和ACE指數的最大值是最小值的2.2倍和2.1倍,麻櫟純林細菌Chao1和ACE指數最大值是最小值的3.0倍和2.9.倍；麻櫟-刺槐混交林細菌Shannon和Simpson指數的最大值是最小值的1.1倍和1.0倍,麻櫟純林細菌Shannon和Simpson指數最大值是最小值的1.3倍和1.0倍。麻櫟-刺槐混交林真菌Chao1和ACE指數最大值均是最小值的1.1倍,麻櫟純林真菌Chao1和ACE最大值均是最小值的1.5倍；麻櫟-刺槐混交林真菌Shannon和Simpson指數最大值是最小值的1.1倍和1.0倍,麻櫟純林細菌Shannon和Simpson指數最大值是最小值的1.3倍和1.1倍。

表3 麻櫟-刺槐混交林和麻櫟純林土壤細菌、真菌菌群豐富度指數和多樣性指數(平均值±標準差)

2.3 麻櫟-刺槐混交林和麻櫟純林凋落物分解速率與土壤微生物多樣性的關系

試驗觀測時間為2019年5月—2020年1月。為研究兩種林分凋落物分解速率與土壤微生物多樣性的關系,以2019年5—7月,2019年8—10月,2019年11月到2020年1月的凋落物月平均質量損失量分別代表3次土壤微生物取樣時凋落物的分解速率,然后與土壤微生物菌群豐富度指數、菌群多樣性指數值之間進行Spearman 相關分析。結果表明,兩種林分凋落物分解速率與細菌菌群豐富度指數(Chao1和ACE指數)呈顯著正相關,與細菌菌群多樣性指數(Shannon和Simpson指數)呈顯著正相關(表4)。兩種林分凋落物分解速率與真菌菌群豐富度指數(Chao1和ACE指數)呈顯著正相關,與真菌菌群多樣性指數(Shannon和Simpson指數)呈正相關但未達到顯著水平。

表4 麻櫟-刺槐混交林和麻櫟純林凋落物分解速率與土壤微生物菌群豐富度指數、多樣性指數的相關系數

3 討論

3.1 麻櫟-刺槐混交林和麻櫟純林凋落物分解過程中元素遷移特征

在凋落物分解過程中會伴隨著元素含量的變化,這與凋落物類型、分解環境、分解時間以及元素本身的特點有關[17-18]。碳是組成植物干物質的重要元素,凋落物分解后,碳一部分會被植物根系吸收后重新利用,一部分會固定在土壤團聚體中[19]。本研究發現,在凋落物分解過程中,麻櫟-刺槐混交林和麻櫟純林凋落物的TOC濃度整體呈現出下降的趨勢(圖3),而且TOC殘留率小于100%(表2)。這是因為在凋落物分解初期,由于可溶性物質的淋溶作用,凋落物TOC含量會迅速減少；隨分解的進行,大分子有機物逐漸降解為小分子有機物,并被微生物和土壤動物等吸收轉換,凋落物TOC含量繼續減少,所以TOC在整個分解過程中表現出凈釋放[20]。

氮是植物氨基酸、蛋白質、核酸的重要組成元素,凋落物中的氮元素通過土壤微生物的氨化和硝化作用轉化成可被植物吸收的NH4+和NO3-,引起凋落物TN含量發生變化[21]。在本研究中,麻櫟-刺槐混交林和麻櫟純林凋落物TN濃度隨分解的進行表現出增加的趨勢(圖3),但是TN含量表現為凈釋放(表2)。凋落物TN的釋放或積累與凋落物中TN濃度有關,在TN濃度滿足微生物生長時表現為凈釋放,反之,表現為凈積累[22]。本研究中,麻櫟-刺槐混交林凋落物TN平均濃度為16.6 g/kg,麻櫟純林凋落物TN平均濃度為13.6 g/kg,兩種林分的凋落物TN濃度均高于全球平均水平10.9 g/kg[23],認為凋落物TN濃度處于較高水平,可以滿足微生物正常生長、發育、繁殖的需求,微生物物會不斷分解凋落物中的含氮有機物,所以TN表現為凈釋放[20]。

磷是組成植物遺傳物質的重要元素,磷的釋放和積累對維持植物正常的生長代謝具有重要意義[24]。本研究發現,凋落物分解過程中,麻櫟-刺槐混交林和麻櫟純林凋落物TP濃度表現出先升高后降低的趨勢(圖3),TP在前期表現出凈積累的過程中,隨后表現出凈釋放的過程(表2)。這是因為在凋落物分解前期TP濃度較低,還不能被微生物利用,所以會發生暫時的性積累。隨著分解的進行,凋落物內部結構發生變化,TP濃度逐漸升高,此時凋落物內的TP在微生物作用下被礦化,再表現出TP凈釋放[20,25]。

研究發現,麻櫟-刺槐混交林凋落物TOC、TN、TP釋放率高于麻櫟純林(表2)。這是因為混交林與純林相比,凋落物的分解環境往往會發生改變,使混交林土壤微生物群落具有較高的代謝強度和豐富度導致的[26]。李英花等[27]對落葉松人工混交林凋落物分解過程進一步進行了研究,發現落葉松和紅松混交林對凋落物養分釋放表現出促進作用,這與本研究研究結果一致,她認為混交林對凋落物養分釋放的影響與樹種組成和凋落物混合比例有關。因此,在今后的研究中應加大不同混交樹種及混交比例對凋落物養分釋放影響的研究。

3.2 麻櫟-刺槐混交林和麻櫟純林土壤微生物多樣性變化分析

森林生態系統會因林分結構、凋落物性質、季節更替產生差異,這些差異可能會對土壤細菌和真菌群落多樣性產生影響[28]。表征土壤微生物多樣性常用的指標有菌群豐富度指數(Chao1和ACE)和菌群多樣性指數(Shannon和Simpson)[29]。在麻櫟-刺槐混交林和麻櫟純林凋落物分解過程中,土壤微生物菌群豐富度指數和菌群多樣性指數發生顯著變化。土壤溫度的差異可能是導致土壤微生物群落結構變化最重要的原因[30]。本試驗土壤微生物三次采樣時間分別為2019年5月,2019年9月和2020年1月,土壤溫度在2019年9月最高,2020年1月最低。一般來說,每種微生物都有其生存的最適合溫度,土壤溫度的變化會影響微生物的生長、繁殖過程,進而影響菌群豐富度和多樣性[31]。有研究表明,溫度在3—25℃范圍內,溫度升高后會使部分微生物群落代謝低溫時不能代謝的物質,增強微生物活性、菌群豐富度和多樣性[32]。另外,低溫可能會降低微生物膜脂質流動性,導致細胞內液冰凍、細胞破裂和死亡,從而降低微生物多樣性[33]。研究發現,麻櫟-刺槐混交林土壤微生物菌群豐富度指數和菌群多樣性指數變化范圍小于麻櫟純林,這可能是混交林能為土壤微生物生長提供更加穩定的微生態環境引起的。王凌云等[34]對楓香-木荷混交林和馬尾松純林土壤微生態環境差異進行了研究,發現香楓-木荷混交林對改善土壤微生物生存環境的改善效果更明顯。黃雅麗等[35]研究認為混交林能從宏觀上改變植被類型和林內植被群落, 也從微觀上改變土壤微生物群落, 使混交林較純林更加穩定,進而為微生物生長提供一個穩定的微生態環境。

3.3 土壤微生物群落對麻櫟-刺槐混交林和麻櫟純林凋落物分解速率的影響

森林生態系統中,凋落物分解雖然受多種因子影響,但多是土壤微生物作用下對諸多生態因子綜合效應的反應[3]。在土壤微生物群落中,細菌數量最大,種類最多,許多群落具有分解凋落物木質素、半纖維素和纖維素的能力,細菌菌群豐富度和菌群多樣性會對凋落物的分解進程產生影響[36]。研究發現,凋落物分解速率與土壤細菌菌群豐富度指數和多樣性指數呈顯著正相關(表4),表明高的細菌豐富度和多樣性對凋落物分解具有促進作用。許多學者對凋落物分解速率與細菌群落多樣性的關系進行了研究,得出了類似的結論。李俊等[37]對方枝柏和康定柳凋落物分解進行了研究,發現凋落物分解速率與細菌菌群豐富度呈正相關。路穎等[16〗在研究泰山4種優勢造林樹種凋落物分解對細菌群落結構的影響時,發現凋落物分解速率與細菌菌群多樣性指數具有很大的相關性。土壤真菌是土壤微生物的另一大類,本研究中麻櫟-刺槐混交林和麻櫟純林凋落物的分解速率與真菌菌群豐富度指數呈顯著正相關,表明真菌菌群豐富度影響凋落物的分解速率。真菌的營養方式為分解吸收式,在凋落物分解過程中能夠產生大量促進凋落物分解的酶,進而有效促進凋落物的分解[38]。

混交林作為一種常見的營林類型,能夠改善森林群落結構,增強森林生態服務功能,會對凋落物有機碳、全氮、全磷以及土壤溫度、土壤水分等方面產生一些獨特的變化,這會對土壤微生物多樣性產生影響,進而會影響凋落物的分解速率[39-40]。大多數研究認為,混交林與純林相比,往往能增加環境異質性和為微生物提供多樣性的食物,從而增加土壤微生物多樣性,提高凋落物分解速率[41]。但Hooper等[42]研究發現,一些凋落物混合會與微生物發生拮抗作用,降低微生物多樣性,對分解起抑制作用。而Nielsen等[43]研究發現,混交林對土壤微生物多樣性并無直接影響,對分解速率的影響不明顯。研究發現,麻櫟-刺槐混交林凋落物分解速率高于麻櫟純林,這可能與凋落物質量差異引起的土壤微生物菌群豐富度和多樣性的差異有關。一般認為,當初始C/N高于30時,隨C/N的增高,氮元素會成為限制性因素,凋落物的微生物分解過程因受抑制而減慢[44]。本研究中,凋落物初始C/N在麻櫟-刺槐混交林中為31.3,麻櫟純林中為39.5,麻櫟-刺槐混交林較低的C/N會為微生物生長和繁殖創造有利條件,增加土壤微生物多樣性[45]。與麻櫟純林相比,麻櫟-刺槐混交林可溶性氮源種類及含量更加豐富,能為土壤微生物提供較為豐富和多樣的食物,增加土壤微生物多樣性,提高凋落物分解速率[46]。

4 結論

兩種林分土壤微生物優勢類群相同,細菌優勢類群為放線菌門、變形菌門、酸桿菌門和疣微菌門,真菌優勢類群為擔子菌門、子囊菌門和被孢霉門。凋落物分解過程中,麻櫟-刺槐混交林土壤微生物菌群豐富度指數和菌群多樣性指數變化范圍小于麻櫟純林。

麻櫟-刺槐混交林凋落物分解速率高于麻櫟純林,土壤微生物多樣性對凋落物分解速率具有重要影響。凋落物分解速率與土壤細菌菌群豐富度指數和菌群多樣性指數呈顯著正相關,與土壤真菌菌群豐富度指數呈顯著正相關。