环境科学  2022, Vol. 43 Issue (1): 550-559   PDF    
秦岭中段撂荒地植被恢复过程中土壤微生物代谢特征
薛悦1,2, 康海斌1, 杨航1, 冰德叶1, 晁志1, 张凯1, 王得祥1,2     
1. 西北农林科技大学林学院, 杨凌 712100;
2. 国家林业和草原局黄土高原林木培育重点实验室, 杨凌 712100
摘要: 明确撂荒地植被恢复过程中土壤微生物养分限制特征及其驱动机制对揭示土壤养分循环和维持生态系统稳定具有重要意义.以秦岭中段不同年限撂荒地为研究对象,通过测定土壤理化性质、5种胞外酶活性[β-1,4-葡萄糖苷酶(BG)、纤维二糖水解酶(CBH)、β-1,4-N-乙酰基氨基葡萄糖苷酶(NAG)、亮氨酸氨基肽酶(LAP)和酸性磷酸酶(AP)],运用胞外酶化学计量学模型,探讨撂荒地植被恢复过程中的土壤微生物养分限制状况及其与土壤性质之间的关系.结果表明,耕地撂荒显著增加了BG、CBH、NAG、LAP和AP的活性;随撂荒年限增加,(BG+CBH)/(NAG+LAP)和(BG+CBH)/AP显著降低,而(NAG+LAP)/AP则相反,相关性分析表明大多数土壤性质与胞外酶活性及其计量比显著相关.酶化学计量向量长度随撂荒年限的增加而减小,表明土壤微生物受碳(C)限制减缓,向量角度均大于45°且呈现降低趋势,表明微生物受到磷(P)限制且逐渐减轻.回归分析表明土壤C和P限制与土壤全量养分、速效养分、养分计量比和物理性质显著相关.偏最小二乘路径模型结果显示,养分计量比是微生物C和P限制的主要驱动因素.综上所述,秦岭植被恢复过程中微生物代谢特征反映了微生物参与下的土壤养分循环机制,为揭示该地区植被恢复过程中群落动态及稳定性,维持区域生态环境安全提供了理论依据.
关键词: 撂荒地      酶化学计量      微生物养分限制      植被恢复      秦岭     
Extracellular Enzyme Stoichiometry and Microbial Metabolism Limitation During Vegetation Restoration Process in the Middle of the Qinling Mountains, China
XUE Yue1,2 , KANG Hai-bin1 , YANG Hang1 , BING De-ye1 , CHAO Zhi1 , ZHANG Kai1 , WANG De-xiang1,2     
1. College of Forestry, Northwest A&F University, Yangling 712100, China;
2. Key Laboratory of Forest Cultivation on the Loess Plateau, State Forestry and Grassland Administration, Yangling 712100, China
Abstract: Clarifying the characteristics of soil microbial nutrient limitation and its driving mechanisms during vegetation restoration after farmland abandonment has important implications for revealing soil nutrient cycling and maintaining ecosystem stability. To determine the limitation of soil microbial nutrients and its relationship with soil properties along a chronosequence of abandoned farmland in the middle of the Qinling Mountains, the soil physicochemical properties and five enzyme activities (β-1, 4-glucosidase (BG), cellobiohydrolase (CBH), β-1, 4-N-acetylglucosaminidase (NAG), leucine aminopeptidase (LAP), and acid phosphatase (AP)) were measured, and models of extracellular enzymatic activity were applied. The results showed that the activities of BG, CBH, NAG, LAP, and AP were significantly increased following farmland abandonment. With the increasing years of abandonment, the ratios of (BG+CBH)/(NAG+LAP) and (BG+CBH)/AP significantly decreased, whereas the ratio of (NAG+LAP)/AP increased. Correlation analysis showed that most soil physicochemical properties were significantly correlated with extracellular enzyme activities and extracellular enzymatic stoichiometry. The vector length of extracellular enzymatic stoichiometry decreased with the increase in abandonment years, indicating that the limitation of soil microorganisms on carbon (C) was reduced. Moreover, the vector angles (>45°) showed a decreasing trend, indicating that microbial metabolisms were limited by phosphorus (P) and gradually decreased. Regression analysis showed that the C and P limitations were significantly related to total nutrients, available nutrients, nutrient ratio, and soil physical properties. Partial least squares path modeling (PLS-PM) revealed that the C and P limitations were directly regulated by nutrient ratio. PLS-PM further showed that soil total nutrients indirectly affected soil microbial C and P limitations by affecting nutrient ratio, and nutrient ratio affected the soil metabolism limitation via available nutrients and pH. Our study suggests that the characteristics of microbial metabolism during the vegetation restoration process reflect the mechanism of microorganism-mediated soil nutrient cycling, which provides a theoretical basis for revealing the community dynamics and stability during the vegetation restoration process and maintaining the regional ecological environment security in the Qinling Mountains.
Key words: farmland abandonment      enzymatic stoichiometry      microbial nutrient limitation      vegetation restoration      Qinling Mountains     

随着经济社会的快速发展, 一些立地环境差的山区和林区, 出现了大面积的撂荒地.耕地在撂荒后恢复过程中, 原有的植被类型格局发生了变化, 立地环境也随之得到改善[1, 2], 这种动态变化进程影响着陆地生态系统的养分平衡.植被恢复会引起地上植物的物种组成、生物量和地下土壤养分、微生物群落等的变化[3, 4].土壤养分状况是影响植物群落建成和演替的主要因素之一, 地上植物的生长, 针阔叶树种比例的变化影响了植物凋落物及根系分泌物的输入从而改变了土壤养分状况[5~7].土壤微生物通过矿化和分解土壤有机质提供更多的可利用性养分, 在植物生长和土壤养分循环中发挥了至关重要的作用[8, 9].目前, 植被恢复的研究主要集中在植物多样性、土壤养分和微生物群落等方面[10~14], 而撂荒地植被恢复过程中, 微生物代谢特征及其影响因素等研究较少.因此, 明确土壤微生物代谢对植被恢复的响应特征及其驱动因素对了解植被恢复过程中土壤养分现状和维持区域生态系统土壤养分平衡具有重要意义.

土壤微生物主要通过胞外酶分解土壤有机质获取碳和养分, 这一过程也会为植物提供更多的有效养分[15, 16], 对土壤物质循环和能量流动起着一定的调控作用[17, 18].其中, 主要有与碳(C)循环相关的β-1, 4-葡萄糖苷酶(β-1, 4-glucosidase, BG)、纤维二糖水解酶(cellobiohydrolase, CBH), 与氮(N)循环相关的β-1, 4-N-乙酰基氨基葡萄糖苷酶(β-1, 4-N-acetylglucosa-minidase, NAG)、亮氨酸氨基肽酶(leucine aminopeptidase, LAP)和与磷(P)循环相关的酸性磷酸酶(acid phosphatase, AP), 这些酶主要参与高分子有机聚合物降解末端的催化反应, 能有效反映土壤微生物的能量(C)和养分(N、P)代谢水平[17, 19].

微生物代谢限制相比于土壤养分和微生物群落对环境变化更为灵敏[20].酶化学计量能反映微生物代谢与养分需求间的生态化学平衡, 常用于评价土壤微生物能量和养分的限制特征[17, 18, 21].当土壤微生物受到养分限制时, 为满足自身生长和代谢需求会分泌相关的胞外酶从而引起酶计量比的变化[22].Sinsabaugh等[18]的研究首先发现在全球尺度上土壤C、N和P获取酶的生态计量比对数转换后收敛于1 ∶1 ∶1, 微生物受到了C的限制, 此后大量的学者对不同生态系统的酶计量比进行了研究.Waring等[23]的研究发现热带森林较温带森林土壤酶C ∶P、N ∶P偏低, 受到了P限制.曾泉鑫等[24]的研究发现N添加刺激微生物分泌P获取酶, 加剧了毛竹林微生物C、P限制.此外, Moorhead等[25]的研究提出了酶计量矢量分析模型, 用向量长度和向量角度分别量化微生物受到的C限制和养分(N、P)限制.基于此, 王强等[26]的研究发现过量的生物炭会导致土壤微生物P限制的加剧, 进而使得C限制增强.Deng等[27]的研究发现黄土高原退耕后植被恢复过程中土壤微生物主要受到C和P限制.因此胞外酶计量比和酶计量矢量分析模型有助于进一步理解植被恢复过程中微生物的养分限制特征及其影响因素.受环境因素的影响, 不同生态系统的胞外酶化学计量特征存在较大的差异[23, 28].Xu等[29]的研究发现我国东部南北样带不同森林生态系统土壤酶计量比与气温和降水显著相关外, 还与土壤pH、土壤养分计量比和叶片化学计量比有关.Zhang等[30]对黄土高原撂荒地和人工刺槐林研究发现, 土壤速效养分是影响酶化学计量的主要因素.虽然, 影响特定生态系统酶计量比的土壤因素受到关注, 但各因子调控的相对大小和途径仍不明确.

秦岭地区植被群落演替及稳定对调节陆地生态系统平衡和改善生态环境具有重要意义[31].目前, 该区域植被恢复研究主要集中在地上植物多样性变化和地下土壤微生物群落结构组成等方面[32], 较少关注植被恢复过程中土壤微生物代谢限制特征及其影响因素.本研究以秦岭中段不同撂荒年限恢复的典型植被为对象, 从酶计量比和酶计量矢量分析模型角度探究恢复过程中土壤微生物代谢的主要养分限制, 以及植被恢复过程中, 微生物代谢限制的变化特征及关键驱动因素, 旨在为秦岭地区植被恢复和森林经营提供理论依据.

1 材料与方法 1.1 研究区概况

研究区位于陕西省宁陕县火地塘林区(33°18′~33°28′N, 108°21′~108°39′E), 该林区地处秦岭中段南坡, 属于北亚热带和暖温带过渡地带, 海拔1 450~2 473 m, 年平均气温8~10℃, 年降水量1 000~1 200 mm, 年蒸发量800~950 mm, 年日照时数1 100~1 300 h, 无霜期170 d.土壤以山地棕壤为主, 土层厚度45~105 cm.上世纪林区人口大量流动, 集中出现了不同年限的撂荒地, 撂荒地植被逐渐恢复, 现植被恢复类型有草丛、灌丛和天然林.乔木优势种有锐齿槲栎(Quercus aliena var. acuteserrata)、油松(Pinus tabuliformis)、华山松(P. armandii)、栓皮栎(Q. variabilis)、铁杉(Tsuga chinensis)和红桦(Betula albosinensis)等, 灌木优势种有白檀(Symplocos paniculata)、榛子(Corylus heterophylla)、南方六道木(Abelia dielsii)、胡枝子(Lespedeza bicolor)、忍冬(Lonicera japonica)和卫矛(Euonymus alatus), 草本优势种有大披针叶苔草(Carex lanceolata)、艾篙(Artemisia lavandulaefolia)、铁杆篙(A. gmelinii)、唐松草(Thalictrum aquilegiifolium var. sibiricum)和茜草(Rubia cordifolia)等.

1.2 样地设置

2019年7~9月, 在火地塘林区选择恢复1、15和30 a的典型撂荒植被作为样地(表 1图 1), 选择种植玉米的耕地作为对照样地, 同时选择未砍伐过的天然林代表恢复至顶级群落的样地, 共计5个恢复阶段.每个阶段设置5个样方, 共25块.林地阶段设置样方面积为20 m×20 m, 灌木阶段样方面积为5 m×5 m, 草地阶段样方面积为2 m×2 m.所有样方背景信息(海拔、坡向和坡度等)相似, 其中, 样方间距大于20 m, 不同恢复阶段样地间距大于300 m, 所有样地均在同一流域.

表 1 样地植被特征 Table 1 Vegetation characteristics of sampling sites

图 1 研究区位置示意 Fig. 1 Location of the study area

1.3 样品采集

在样方内用土钻(口径4.5 cm)按“S”形取土样, 其中, 林地阶段选取12个采样点, 灌草阶段选取8个采样点.取样前去除枯落物, 取样深度为0~20 cm, 将同一样方内的土样混合, 过筛(孔径2 mm), 剔除石头、根等杂质装入冰盒带回实验室.将土样分为两部分, 一部分在实验室内自然风干, 用于后续土壤理化性质的测定; 一部分放入4℃冰箱, 用于胞外酶活性的测定.在样方对角线设置3个样点, 用体积100 cm3环刀随机取0~10 cm土层的土样, 并装入铝盒带回实验室立即称重, 随后放入烘箱105℃烘干至恒重, 以测定土壤容重.

1.4 样品分析 1.4.1 土壤理化性质测定

采用重铬酸钾容量法测定土壤有机碳(SOC), 采用凯氏定氮法测定土壤全氮(TN), 土壤全磷(TP)采用H2SO4-HClO4消解, 钼锑抗比色法测定.土壤可溶性有机碳(DOC)用K2SO4浸提, 总有机碳分析仪(TOC-VCPH, 日本)测定, 土壤硝铵态氮(NO3--N和hNH4+-N)用KCl浸提, 连续流动分析仪(Auto Analyzer 3-AA3, 德国)测定, 速效磷(Olsen-P)用NaHCO3浸提, 钼锑抗比色法测定.使用pH计(奥豪斯ST2100, 上海)测定土壤pH值, 水土比为2.5 ∶1.土壤含水量(SWC)采用烘干法测定.土壤容重(SBD)采用环刀法测定.

1.4.2 胞外酶活性测定

采用96微孔板荧光法测定相关酶: β-1, 4-葡萄糖苷酶(BG)、纤维二糖水解酶(CBH)、β-1, 4-N-乙酰基氨基葡萄糖苷酶(NAG)、亮氨酸氨基肽酶(LAP)和酸性磷酸酶(AP).酶底物为4-甲基-7-羟基香豆素(MUB), 称取1 g鲜土加入125 mL蒸馏水在25℃恒温摇床上振荡2 h(240 r ·min-1), 使之混合均匀, 用移液器吸取1 mL悬浊液到离心管内, 分别加入0.25 mL蒸馏水作为空白管、0.25 mL标准物质溶液为淬火管以及0.25 mL对应荧光底物为测定样品管, 摇匀后25℃黑暗条件下培养4 h, 随后每个离心管加入50 μL 0.5 mol ·L-1的NaOH终止反应, 最后分别吸取250 μL样品、空白、淬火、标准和阴性溶液到96孔酶标板中, 用多功能酶标仪(M200PRO, 瑞士)在365 nm激发波长和450 nm检测波长下进行测定.

1.5 数据处理

通过计算酶化学计量的向量长度(vector length, VL)和向量角度(vector angle, VA)分析土壤微生物受到养分限制的情况, 即[25]:

式中, VL越大表明微生物受C限制越大; VA < 45°时表明微生物受N限制; VA>45°时表明微生物受P限制, VA越大微生物P限制越大.

采用R 3.6.2软件对实验数据进行统计分析.运用单因素方差分析法(one-way ANOVA)分析不同撂荒年限土壤理化性质、胞外酶活性及其化学计量比、酶化学计量向量长度及角度的差异, 并采用Tukey's HSD进行事后多重比较.对酶活性、酶化学计量比及其向量特征与土壤性质进行Pearson相关性分析.采用广义线性模型(GLM)确定微生物代谢限制与环境变量的关系.选择重要的环境变量利用R中的“plspm”程序包构建偏最小二乘路径模型(PLS-PM), 进一步探究控制植被恢复过程中微生物代谢限制的可能路径.

2 结果与分析 2.1 不同撂荒年限耕地土壤理化性质特征

不同撂荒年限耕地的土壤理化性质变化特征见表 2, 土壤SOC、TN、C ∶N、C ∶P、N ∶P和NO3-+NH4+-N随撂荒年限的增加呈显著上升的趋势(p < 0.05), 在撂荒30 a时达到最大, 撂荒30 a时土壤C ∶N显著小于天然林, 而TN显著大于天然林.土壤TP、DOC和Olsen-P呈先增大后减小的趋势(p < 0.05), 在撂荒15 a时最大, 天然林土壤Olsen-P和TP显著小于撂荒15 a.天然林土壤pH、SBD和SWC显著小于耕地和1、15和30 a撂荒地(p < 0.05).

表 2 不同撂荒年限耕地土壤理化性质1) Table 2 Soil physicochemical properties in abandoned farmlands in different years

2.2 不同撂荒年限耕地胞外酶活性及酶计量比特征

图 2可见, 不同撂荒年限耕地对土壤BG、CBH、NAG、LAP和AP活性有显著的影响(p < 0.05).随着撂荒年限的增加, 土壤BG、CBH、NAG、LAP和AP呈上升趋势, 并且均在撂荒30 a时达到最大, 与耕地土壤相比分别增加了63.6%、193.4%、274.7%、196.2%和97.9%.天然林土壤CBH活性和LAP活性显著大于撂荒30 a土壤, 而土壤NAG活性则相反.

不同小写字母表示不同撂荒年限耕地之间差异显著, 下同 图 2 不同撂荒年限耕地胞外酶活性变化 Fig. 2 Changes in extracellular enzymatic activities in abandoned farmlands in different years

图 3可见, 随着撂荒年限的增加, (BG+CBH)/(NAG+LAP)和(BG+CBH)/AP呈显著下降趋势(p < 0.05), 在撂荒30 a时达到最低, 与耕地土壤相比分别降低了55.2%, 24.0%.而(NAG+LAP)/AP呈相反趋势(p < 0.05), 在撂荒30 a时达到最高, 比耕地土壤增加了81.8%.天然林(NAG+LAP)/AP显著大于撂荒30 a, 相比增加了52.5%, 而(BG+CBH)/(NAG+LAP)和(BG+CBH)/AP则相反分别降低了45.2%和18.4%.天然林土壤胞外酶C ∶N ∶P[用ln(BG+CBH) ∶ln(NAG+LAP) ∶ln AP表示]接近于1 ∶1 ∶1.13.此外, 由表 3可见, 土壤胞外酶活性及酶计量比与土壤全量养分、速效养分、养分计量比及物理性质有显著的相关性(表 3), 如(BG+CBH)/(NAG+LAP)与SOC、C ∶N、C ∶P、DOC和NO3-+NH4+-N呈显著负相关, 与TP、Olsen-P、SBD和pH呈显著正相关(p < 0.05).

图 3 不同撂荒年限耕地胞外酶计量比变化 Fig. 3 Changes in extracellular enzymatic stoichiometry in abandoned farmlands in different years

表 3 土壤理化性质与微生物代谢特征的Pearson相关性1) Table 3 Pearsoncorrelations between the soil physicochemical properties and microbial metabolic characteristics

2.3 不同撂荒年限耕地胞外酶计量向量特征

不同撂荒年限耕地向量特征差异显著(p < 0.05), 向量长度及角度变化范围分别为0.66~0.95和52.92°~59.71°, 变化趋势如图 4.随撂荒年限的增加向量长度显著降低, 撂荒30 a时最小(0.81), 表明随着撂荒年限的增加土壤微生物受到的C限制逐渐减缓, 天然林向量长度(0.66)显著小于撂荒30 a.向量角度变化趋势与长度相同, 撂荒30 a最小(56.36°), 天然次生林向量角度为52.92°显著小于撂荒30 a, 这表明虽然土壤微生物受到P限制逐渐减缓, 但仍然受到较强P限制.此外, 由表 3可见, 土壤理化性质(如SOC、C ∶N、C ∶P、DOC、NO3-+NH4+-N、SBD和pH)与向量长度和角度均显著相关(p < 0.05).

图 4 不同撂荒年限耕地胞外酶计量比向量长度和角度变化 Fig. 4 Changes in vector length and angle of extracellular enzymatic stoichiometry in abandoned farmlands in different years

2.4 微生物代谢与土壤理化性质的关系

线性回归分析表明(图 5图 6), 土壤微生物C限制与土壤SOC、C ∶N、C ∶P、DOC和NO3-+NH4+-N呈显著负相关, 与TP、Olsen-P、pH、SMC和SBD呈显著正相关(p < 0.05).P限制与土壤SOC、C ∶N、C ∶P、N ∶P、DOC和NO3-+NH4+-N呈显著负相关, 与TP、pH和SBD呈显著正相关(p < 0.05).偏最小二乘路径模型分析显示(图 7), 土壤物理性质、总养分、速效养分、养分计量比和pH对微生物C、P限制均存在直接和间接的影响, 如土壤总养分直接影响微生物对C和P的获取, 也可以通过影响养分计量比间接对微生物C和P限制产生影响.养分计量比(-0.96)、土壤总养分(-0.51)和速效养分(-0.48) 对微生物C限制表现为负效应, 而土壤物理性质(0.29)和pH(0.13)则表现为正效应, 其中土壤养分计量比影响最大.对于微生物P限制, 土壤养分计量比(-0.71)、土壤总养分(-0.59)、pH(-0.14)和速效养分(-0.09)表现为负效应, 土壤物理性质(0.36)为正效应, 其中受到土壤养分计量比的影响最大.

实线表示向量长度与土壤理化性质之间的模型拟合, 灰色区域表示模型95%置信区间, 下同 图 5 向量长度与土壤理化性质线性回归 Fig. 5 Linear regression between the vector length and soil physicochemical properties

图 6 向量角度与土壤理化性质线性回归 Fig. 6 Linear regression between the vector angle and soil physicochemical properties

1表示土壤物理性质, 2表示土壤总养分, 3表示土壤养分计量比, 4表示土壤速效养分, 5表示土壤pH; 实线表示正向因果关系, 虚线表示负向因果关系; 箭头上的数字表示标准化路径系数; R2表示模型解释的因变量的方差 图 7 偏最小二乘路径建模分析 Fig. 7 Partial least-squares path modeling (PLS-PM)

3 讨论 3.1 耕地撂荒对胞外酶活性的影响

胞外酶活性在有机质分解和养分循环中起着重要的作用, 可以在一定程度上反映微生物的代谢活动[17, 18].本研究发现, 随着撂荒年限增加土壤C、N和P代谢酶活性增强(图 2).钟泽坤等[22]对黄土高原农田撂荒过程中酶活性变化特征的研究同样发现, 土壤NAG、LAP和AP随着撂荒年限增加而增加.Cui等[33]的研究发现黄土高原农田转变为灌木的植被恢复过程中C、N和P代谢酶活性逐渐增强, 与本研究的结果一致.有研究发现胞外酶活性与土壤有机物的分解有关[15, 16].耕地撂荒后的恢复过程中, 土壤有机质逐渐增多, 使得土壤微生物有了更多的养分来源, 从而刺激微生物分泌大量的胞外酶[34, 35].不同生态系统胞外酶活性的影响因素不同, 本研究表明胞外酶活性与土壤SOC、TP、C ∶P和pH等显著相关(表 3).孙彩丽等[35]的研究发现土壤TN、C ∶N、N ∶P和pH是引起酶活性变化的主要驱动因素.罗攀等[36]的研究发现喀斯特山区土壤酶活性主要受土壤pH、总氮和无机氮的影响.因此, 随着撂荒年限的增加, 土壤物理结构和化学性质逐渐得到改善, 微生物活性增强, 且土壤养分水平总体较低, 使得微生物分泌更多的胞外酶分解有机物.

3.2 耕地撂荒对土壤微生物代谢限制的影响

胞外酶计量比可以为生物的代谢需求和环境养分可利用性提供信息, 进而反映土壤的生物化学过程和养分限制情况[21].结合资源分配学说, 土壤微生物会通过调节自身酶的合成, 提高限制元素的有效性, 从而缓解土壤养分的限制.酶计量矢量化可以更加直观分析微生物所受的能量(C)和养分(N、P)限制[37].本研究中, 胞外酶计量向量分析发现, 微生物受到C和P的共同限制[撂荒后向量角度均大于45°且(NAG+LAP)/AP均低于全球尺度水平0.44[18], 图 4]. Xiao等[38]对农田植被恢复过程中酶计量矢量化研究发现, 微生物受到C和P的共同限制, 与本研究的结果一致.Cui等[39]的研究同样发现黄土高原的不同植被恢复过程中微生物受到C和P的共同限制.本研究还发现, 随着撂荒时间的延长(BG+CBH)/(NAG+LAP)和(BG+CBH)/AP呈下降趋势(图 3), 表明微生物对N、P的需求大于对C的需求.同样, 胞外酶化学计量的向量长度逐渐减小, 核实了土壤微生物受到的C限制减弱, 逐渐转化为养分限制.钟泽坤等[22]的研究发现农田撂荒后向量长度及土壤BG ∶(NAG+LAP)和BG ∶ALP逐渐减小, 微生物C限制逐渐减弱, 与本研究的结果一致.撂荒地经过植被恢复过程, 凋落物及根系分泌物等增加, 使得土壤有机质含量显著增加, 从而缓解了微生物C限制.本研究中, 撂荒后土壤N ∶P随撂荒时间呈增大趋势(表 2), 表明植被恢复中P是生态系统主要的限制因素[38, 40].这可能是由于土壤中的P元素主要来自岩石的风化, 长时间的植被生长已经将本来不富足的P消耗殆尽, 且研究区的土壤呈酸性, P元素易被络合[41], 另外撂荒后豆科植物(如圆锥山蚂蝗)的出现和有机物的输入增加了N素的含量, 因此土壤微生物N限制得到逐步缓解, 表现为受P限制.本研究中, (NAG+LAP)/AP随撂荒时间的增加而增加(图 3).这可能是由于耕地撂荒后主要受到P元素的限制, 但植物多样性和丰富度下降, 植物根系与微生物的养分竞争缓解, 导致微生物P限制逐渐降低, 且N元素也不富裕.综上, 胞外酶计量比及其矢量分析模型揭示了撂荒地植被恢复过程中土壤微生物主要受到C和P限制, 且均随撂荒年限的增加而逐渐减缓.

3.3 土壤微生物C、P限制与土壤理化性质的关系

虽然随撂荒年限延长, 土壤质量得到改善, 有机物含量增加, 进一步影响土壤微生物C和养分获取, 但限制微生物C和P的直接因素仍不清楚.本研究发现土壤C ∶N和C ∶P均与向量长度有显著相关性(图 5), 土壤C ∶N、C ∶P和N ∶P与向量角度显著相关(图 6), 表明养分计量比显著影响微生物C和P获取.偏最小二乘路径模型分析同样显示养分计量比对微生物C和P限制影响最大.植被恢复过程中土壤养分计量比调控着微生物C和P限制[20, 42].在不同的植被类型下, 土壤养分水平不同, 而养分化学计量比是否平衡是影响微生物活性的主要原因.Deng等[27]在黄土高原植被恢复过程中土壤微生物代谢限制的驱动因素研究中同样发现, C ∶P对微生物C和P限制影响最大.Cui等[20]的研究也发现土壤C ∶P和N ∶P与微生物C和P限制显著相关.这可能是因为土壤养分化学计量比对微生物群落的结构和活性有很大影响, 进而影响微生物的养分限制.本研究中, 土壤SOC、DOC和NO3--N+NH4+-N与向量长度及角度呈负相关, 而TP、SBD和pH与向量长度和向量角度呈正相关(图 5图 6), 表明土壤全量养分、速效养分、物理性质和pH也显著影响微生物C和P限制.Xiao等[38]的研究发现黄土高原次生演替过程中微生物C和P限制与土壤全量养分显著相关.Barantal等[43]对亚马孙热带森林有机质分解的能量与养分限制研究中发现, 分解者的能量限制受土壤速效养分的影响.偏最小二乘路径模型进一步说明土壤全量养分直接或通过影响养分计量比间接影响微生物C和P限制, 养分计量比直接或通过影响速效养分、pH, 从而影响C和P限制(图 7).本研究中, 土壤SOC、DOC、C ∶N和C ∶P随撂荒年限延长显著增大(表 2)且与微生物C限制呈显著负相关(图 5), 表明植被恢复为土壤微生物提供了更多的碳源, 从而缓解了微生物的C限制.综上, 通过回归分析和偏最小二乘路径模型分析, 本研究发现植被恢复过程中土壤养分计量比是影响微生物C和P限制主要因素.

4 结论

(1) 秦岭中段撂荒地植被恢复提高了土壤SOC、TN和NO3--N+NH4+-N的含量, 改变了土壤物理性质, 改善了土壤质量.

(2) 秦岭中段耕地撂荒后BG、CBH、NAG、LAP和AP的活性显著增加.

(3) 撂荒地植被恢复过程中土壤微生物主要受到了C和P限制, 且随撂荒年限的增加微生物C、P限制逐渐减缓.

(4) 土壤养分计量比取决于全量养分, 土壤养分计量比直接介导了微生物C和P限制.

(5) 本研究通过明确胞外酶化学计量比与土壤养分计量比的耦合关系, 揭示了撂荒地植被恢复过程中微生物养分限制特征, 为进一步明确秦岭中段撂荒地植被恢复过程中微生物代谢特征和土壤养分循环提供了理论依据.

参考文献
[1] Knelman J E, Graham E B, Trahan N A, et al. Fire severity shapes plant colonization effects on bacterial community structure, microbial biomass, and soil enzyme activity in secondary succession of a burned forest[J]. Soil Biology and Biochemistry, 2015, 90: 161-168. DOI:10.1016/j.soilbio.2015.08.004
[2] 高涵, 肖礼, 牛丹, 等. 宁南山区退耕还林还草对土壤氮素组成及其转化酶活的影响[J]. 环境科学, 2019, 40(8): 3825-3832.
Gao H, Xiao L, Niu D, et al. Effects of converting farmland into forest and grassland on soil nitrogen component and conversion enzyme activity in the mountainous area of southern Ningxia[J]. Environmental Science, 2019, 40(8): 3825-3832.
[3] van Hall R, Cammeraat L H, Keesstra S D, et al. Impact of secondary vegetation succession on soil quality in a humid Mediterranean landscape[J]. GETENA, 2017, 149: 836-843.
[4] Cline L C, Zak D R. Soil microbial communities are shaped by plant-driven changes in resource availability during secondary succession[J]. Ecology, 2015, 96(12): 3374-3385. DOI:10.1890/15-0184.1
[5] Zhang W, Zhao J, Pan F J, et al. Changes in nitrogen and phosphorus limitation during secondary succession in a karst region in southwest China[J]. Plant and Soil, 2015, 391(1-2): 77-91. DOI:10.1007/s11104-015-2406-8
[6] Wu Y, Chen W J, Li Q, et al. Ecoenzymatic stoichiometry and nutrient limitation under a natural secondary succession of vegetation on the Loess Plateau, China[J]. Land Degradation & Development, 2020, 32(1): 399-409.
[7] 陈金磊, 张仕吉, 李雷达, 等. 亚热带不同植被恢复阶段林地凋落物层现存量和养分特征[J]. 生态学报, 2020, 40(12): 4073-4086.
Chen J L, Zhang S J, Li L D, et al. Stock and nutrient characteristics of litter layer at different vegetation restoration stages in subtropical region, China[J]. Acta Ecologica Sinica, 2020, 40(12): 4073-4086.
[8] Bitas V, Kim H S, Bennett J W, et al. Sniffing on microbes: diverse roles of microbial volatile organic compounds in plant health[J]. Current Review, 2013, 26(8): 835-843.
[9] Lozano Y M, Hortal S, Armas C, et al. Interactions among soil, plants, and microorganisms drive secondary succession in a dry environment[J]. Soil Biology and Biochemistry, 2014, 78: 298-306. DOI:10.1016/j.soilbio.2014.08.007
[10] Sun C L, Chai Z Z, Liu G B, et al. Changes in species diversity patterns and spatial heterogeneity during the secondary succession of grassland vegetation on the Loess Plateau, China[J]. Frontiers in Plant Science, 2017, 8. DOI:10.3389/fpls.2017.01465
[11] 海旭莹, 董凌勃, 汪晓珍, 等. 黄土高原退耕还草地C、N、P生态化学计量特征对植物多样性的影响[J]. 生态学报, 2020, 40(23): 8570-8581.
Hai X Y, Dong L B, Wang X Z, et al. Effects of carbon, nitrogen, and phosphorus ecological stoichiometry characteristics on plant diversity since returning farmland to grassland on the Loess Plateau[J]. Acta Ecologica Sinica, 2020, 40(23): 8570-8581.
[12] 孙沛沛, 钱朝菊, 尹晓月, 等. 半干旱区沙地蒿类植被建成对土壤细菌的影响[J]. 生态学报, 2020, 40(16): 5783-5792.
Sun P P, Qian C J, Yin X Y, et al. Effects of artemisia vegetation built on soil bacteria in semi-arid sandy land[J]. Acta Ecologica Sinica, 2020, 40(16): 5783-5792.
[13] 张芳, 杜虎, 曾馥平, 等. 西南喀斯特峰丛洼地木本植物群落结构与多样性变化[J]. 生态学报, 2020, 40(12): 4094-4104.
Zhang F, Du H, Zeng F P, et al. Changes of woody community structure and diversity in karst peak-cluster depressions in southwest China[J]. Acta Ecologica Sinica, 2020, 40(12): 4094-4104.
[14] 陈孟立, 曾全超, 黄懿梅, 等. 黄土丘陵区退耕还林还草对土壤细菌群落结构的影响[J]. 环境科学, 2018, 39(4): 1824-1832.
Chen M L, Zeng Q C, Huang Y M, et al. Effects of the farmland-to-forest/grassland conversion program on the soil bacterial community in the Loess Hilly Region[J]. Environmental Science, 2018, 39(4): 1824-1832.
[15] Cui Y X, Fang L C, Guo X B, et al. Ecoenzymatic stoichiometry and microbial nutrient limitation in rhizosphere soil in the arid area of the northern Loess Plateau, China[J]. Soil Biology and Biochemistry, 2018, 116: 11-21. DOI:10.1016/j.soilbio.2017.09.025
[16] Manzon S. Flexible carbon-use efficiency across litter types and during decomposition partly compensates nutrient imbalances-results from analytical stoichiometric models[J]. Frontiers in Plant Science, 2017, 8. DOI:10.3389/fmicb.2017.00661
[17] Sinsabaugh R L, Hill B H, Shah J J F. Ecoenzymatic stoichiometry of microbial organic nutrient acquisition in soil and sediment[J]. Nature, 2009, 462(7274): 795-798. DOI:10.1038/nature08632
[18] Sinsabaugh R L, Lauber C L, Weintraub M N, et al. Stoichiometry of soil enzyme activity at global scale[J]. Ecology Letters, 2008, 11(11): 1252-1264. DOI:10.1111/j.1461-0248.2008.01245.x
[19] Merino C, Godoy R, Matus F. Soil microorganisms and enzyme activity at different levels of organic matter stability[J]. Journal of soil science and plant nutrition, 2016, 16: 14-30.
[20] Cui Y X, Bing H J, Fang L C, et al. Extracellular enzyme stoichiometry reveals the carbon and phosphorus limitations of microbial metabolisms in the rhizosphere and bulk soils in alpine ecosystems[J]. Plant and Soil, 2021, 458(1-2): 7-20. DOI:10.1007/s11104-019-04159-x
[21] 左宜平, 张馨月, 曾辉, 等. 大兴安岭森林土壤胞外酶活力的时空动态及其对潜在碳矿化的影响[J]. 北京大学学报(自然科学版), 2018, 54(6): 1311-1324.
Zuo Y P, Zhang X Y, Zeng H, et al. Spatiotemporal dynamics of soil extracellular enzyme activity and its influence on potential mineralization rate of soil organic carbon in forests of Daxing'an mountain range[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2018, 54(6): 1311-1324.
[22] 钟泽坤, 杨改河, 任成杰, 等. 黄土丘陵区撂荒农田土壤酶活性及酶化学计量变化特征[J]. 环境科学, 2021, 42(1): 411-421.
Zhong Z K, Yang G H, Ren C J, et al. Effects of farmland abandonment on soil enzymatic activity and enzymatic stoichiometry in the Loess Hilly Region, China[J]. Environmental Science, 2021, 42(1): 411-421.
[23] Waring B G, Weintraub S R, Sinsabaugh R L. Ecoenzymatic stoichiometry of microbial nutrient acquisition in tropical soils[J]. Biogeochemistry, 2014, 117(1): 101-113. DOI:10.1007/s10533-013-9849-x
[24] 曾泉鑫, 张秋芳, 林开淼, 等. 酶化学计量揭示5年氮添加加剧毛竹林土壤微生物碳磷限制[J]. 应用生态学报, 2021, 32(2): 521-528.
Zeng Q X, Zhang Q F, Lin K M, et al. Enzyme stoichiometry evidence revealed that five years nitrogen addition exacerbated the carbon and phosphorus limitation of soil microorganisms in a Phyllostachys pubescens forest[J]. Chinese Journal of Applied Ecology, 2021, 32(2): 521-528.
[25] Moorhead D L, Sinsabaugh R L, Hill B H, et al. Vector analysis of ecoenzyme activities reveal constraints on coupled C, N and P dynamics[J]. Soil Biology and Biochemistry, 2016, 93: 1-7. DOI:10.1016/j.soilbio.2015.10.019
[26] 王强, 耿增超, 许晨阳, 等. 施用生物炭对塿土土壤微生物代谢养分限制和碳利用效率的影响[J]. 环境科学, 2020, 41(5): 2425-2433.
Wang Q, Geng Z C, Xu C Y, et al. Effects of biochar application on soil microbial nutrient limitations and carbon use efficiency in Lou Soil[J]. Environmental Science, 2020, 41(5): 2425-2433.
[27] Deng L, Peng C H, Huang C B, et al. Drivers of soil microbial metabolic limitation changes along a vegetation restoration gradient on the Loess Plateau, China[J]. Geoderma, 2019, 353: 188-200. DOI:10.1016/j.geoderma.2019.06.037
[28] Peng X Q, Wang W. Stoichiometry of soil extracellular enzyme activity along a climatic transect in temperate grasslands of northern China[J]. Soil Biology and Biochemistry, 2016, 98: 74-84. DOI:10.1016/j.soilbio.2016.04.008
[29] Xu Z W, Yu G R, Zhang X Y, et al. Soil enzyme activity and stoichiometry in forest ecosystems along the North-South Transect in eastern China (NSTEC)[J]. Soil Biology and Biochemistry, 2017, 104: 152-163. DOI:10.1016/j.soilbio.2016.10.020
[30] Zhang W, Xu Y D, Gao D X, et al. Ecoenzymatic stoichiometry and nutrient dynamics along a revegetation chronosequence in the soils of abandoned land and Robinia pseudoacacia plantation on the Loess Plateau, China[J]. Soil Biology and Biochemistry, 2019, 134: 1-14. DOI:10.1016/j.soilbio.2019.03.017
[31] 马姜明, 李昆. 森林生态系统稳定性研究的现状与趋势[J]. 世界林业研究, 2004, 17(1): 15-19.
Ma J M, Li K. Current situation of research and prospects on forest ecosystem stability[J]. World Forestry Research, 2004, 17(1): 15-19.
[32] 张洪武, 于飞, 王得祥, 等. 秦岭南坡松栎林群落演替过程中物种组成及多样性动态分析[J]. 西北植物学报, 2014, 34(1): 169-176.
Zhang H W, Yu F, Wang D X, et al. Changes of species composition and diversity in the process of community succession of pine oak forests on the south-facing slopes in Qinling mountains[J]. Acta Botanica Boreali-Occidentalia Sinica, 2014, 34(1): 169-176.
[33] Cui Y X, Fang L C, Guo X B, et al. Responses of soil bacterial communities, enzyme activities, and nutrients to agricultural-to-natural ecosystem conversion in the Loess Plateau, China[J]. Journal of Soils and Sediments, 2019, 19(3): 1427-1440. DOI:10.1007/s11368-018-2110-4
[34] 郭天斗, 于露, 孙忠超, 等. 宁夏东部荒漠草原向灌丛地人为转变过程中土壤胞外酶活性响应[J]. 应用生态学报, 2020, 31(8): 2541-2548.
Guo T D, Yu L, Sun Z C, et al. Responses of soil extracellular enzyme activities to the anthropogenic transition from desert grassland to shrubland in eastern Ningxia, China[J]. Chinese Journal of Applied Ecology, 2020, 31(8): 2541-2548.
[35] 孙彩丽, 王艺伟, 王从军, 等. 喀斯特山区土地利用方式转变对土壤酶活性及其化学计量特征的影响[J]. 生态学报, 2021, 41(10): 4140-4149.
Sun C L, Wang Y W, Wang C J, et al. Effects of land use conversion on soil extracellular enzyme activity and its stoichiometric characteristics in karst mountainous areas[J]. Acta Ecologica Sinica, 2021, 41(10): 4140-4149.
[36] 罗攀, 陈浩, 肖孔操, 等. 地形、树种和土壤属性对喀斯特山区土壤胞外酶活性的影响[J]. 环境科学, 2017, 38(6): 2577-2585.
Luo P, Chen H, Xiao K C, et al. Effects of topography, tree species and soil properties on soil enzyme activity in karst regions[J]. Environmental Science, 2017, 38(6): 2577-2585.
[37] Sinsabaugh R L, Moorhead D L. Resource allocation to extracellular enzyme production: a model for nitrogen and phosphorus control of litter decomposition[J]. Soil Biology and Biochemistry, 1994, 26(10): 1305-1311. DOI:10.1016/0038-0717(94)90211-9
[38] Xiao L, Li P, Shi P, et al. Soil nutrient stoichiometries and enzymatic activities along an elevational gradient in the dry-hot valley region of southwestern China[J]. Archives of Agronomy and Soil Science, 2019, 65(3): 322-333. DOI:10.1080/03650340.2018.1502882
[39] Cui Y X, Fang L C, Guo X B, et al. Natural grassland as the optimal pattern of vegetation restoration in arid and semi-arid regions: evidence from nutrient limitation of soil microbes[J]. Science of the Total Environment, 2019, 648: 388-397. DOI:10.1016/j.scitotenv.2018.08.173
[40] Jiao F, Wen Z M, An S S, et al. Successional changes in soil stoichiometry after land abandonment in Loess Plateau, China[J]. Ecological Engineering, 2013, 58: 249-254. DOI:10.1016/j.ecoleng.2013.06.036
[41] 殷爽, 王传宽, 金鹰, 等. 东北地区大秃顶子山土壤-微生物-胞外酶C: N: P化学计量特征沿海拔梯度的变化[J]. 植物生态学报, 2019, 43(11): 999-1009.
Yin S, Wang C K, Jin Y, et al. Changes in soil-microbe-exoenzyme C: N: P stoichiometry along an altitudinal gradient in Mt. Datudingzi, Northeast China[J]. Chinese Journal of Plant Ecology, 2019, 43(1): 9999-1009.
[42] Sinsabaugh R L, Shah J J F. Ecoenzymatic stoichiometry and ecological theory[J]. Annual Review of Ecology, Evolution, and Systematics, 2012, 43: 313-343. DOI:10.1146/annurev-ecolsys-071112-124414
[43] Barantal S, Schimann H, Fromin N, et al. Nutrient and carbon limitation on decomposition in an amazonian moist forest[J]. Ecosystems, 2012, 15(7): 1039-1052. DOI:10.1007/s10021-012-9564-9