抚育间伐是常用的森林管理措施^[1]，因伐除林冠相对密集的部分树木，增加了太阳辐射，改变了森林小气候和土壤微生境，必然影响森林生态系统的养分和生物地球化学循环过程，以及该循环过程的核心环节——土壤微生物活动和酶活性。目前，土壤胞外酶研究更多关注于碳、氮和磷循环相关的降解酶，如碳酶[β-葡糖苷酶(BG)、纤维二糖水解酶(CBH)、β-木糖苷酶(BX)]，氮酶[β-1,4-N-乙酰氨基葡萄糖苷酶(NAG)、亮氨酸氨基肽酶(LAP)]和磷酶[酸性或碱性磷酸酶(AcP)]，其活性可作为微生物资源分配的代理指标^[2]。在养分循环期间酶活性的相对丰度变化可反映微生物群落的代谢水平。SINSABAUGH等^[3]最先通过整合分析发现：在全球尺度上碳、氮和磷循环相关酶计量比接近1∶1∶1，表明土壤酶化学计量比呈稳态性。但也有研究发现：土壤酶化学计量比呈非稳态性^[4−6]，说明微生物可能受到能量或关键营养物质(即碳、氮和磷)的限制^[7]。

间伐措施对土壤胞外酶活性和酶化学计量的影响仍不确定。如土壤酶活性在森林间伐后会增加^[8]、减少^[9]或保持不变^[10]。大多数研究主要围绕不同间伐强度对酶活性的影响^[11]。间伐措施的影响效果还会随森林恢复过程而发生改变。如QIU等^[12]对塞罕坝林场内华北落叶松Larix principis-rupprechtii人工林进行间伐恢复9 a后的结果显示：间伐措施显著增加了土壤BG、NAG+LAP和AcP活性。而LULL等^[13]对地中海栎Quercus ilex林间伐后5个月至7 a内，氮和磷循环酶的活性并未发生显著改变。间伐处理和林下移除可在短时间内减少微生物对土壤资源的竞争，进而改变酶的活性^[14]。但随树木生长速度和土壤养分含量的变化，微生物资源利用策略也发生改变，可能造成微生物受到不同养分的限制^[15]。

目前，关于间伐处理对土壤胞外酶活性的研究大多侧重于间伐强度和人工林生态系统的研究，而对天然林生态系统间伐后不同恢复阶段土壤酶活性的研究较少。鉴于此，本研究采用空间代替时间的方法，探讨北亚热带秦岭松栎混交林在抚育间伐后不同恢复时间内林地表层土壤酶活性、酶化学计量比的变化规律，为制定森林可持续经营方案及合理的生态恢复措施提供理论依据。

1.   研究地区与研究方法

1.1   研究区概况

研究区位于陕西省安康市宁东林业局新矿林场(33°20′~33°26′N，108°32′~108°34′E)，地处秦岭山脉，海拔为1 400.0~1 800.0 m。该区属于北亚热带与温带过渡区，年均气温为8.5 ℃，年平均降水量为908.0 mm，土壤为山地棕壤。研究区域为20世纪70年代末采伐后天然更新形成的次生针阔混交林^[16]，采取的是低强度间伐和林冠下补植等保护经营作业法。林内主要以油松Pinus tabuliformis、锐齿槲栎Quercus aliena var. acutiserrata、华山松Pinus armandii为主要建群种，伴生有漆树Toxicodendron vernicifluum、小叶女贞Ligustrum quihoui、青榨槭Acer davidii等树种。林下植被以卫矛Euonymus alatus、木姜子Litsea pungens、披针叶薹草Carex lanceolata、龙牙草Agrimonia pilosa、茜草Rubia cordifolia为主。

2021年10月，根据研究区内实际间伐处理、林木生长和分布状况，选择立地条件基本一致的林分，设置3种间伐处理，即未间伐(ck)、间伐恢复5 a (5 a，2018年间伐)和间伐恢复13 a (13 a，2010年间伐)。每个间伐处理设置4块面积为20 m×30 m的样地，共计12块样地。为防止样地之间相互干扰，样方之间的间隔不小于100 m。进行间伐处理后林下物种数量增加，更新了枫杨Pterocarya stenoptera、栗Castanea mollissima、桤木Alnus cremastogyne、灯台树Cornus controversa和胡桃楸Juglans mandshurica等树种。其中各样地内物种丰富度和Shannon-Wiener指数参照刘思泽等^[17]的方法计算。样地调查基本概况见表1。

表 1  试验样地基本概况

Table 1  Basic survey of test plots

间伐后恢 复时间/a	海拔/ m	株数密度/ (株·hm−2)	胸径/ cm	郁闭度	物种 丰富度	Shannon-Wiener 指数	凋落物量/ (t·hm−2·a−1)	林内主要树种
ck	1 585.00±61.85	1 420±88	14.60±0.49	0.7	25	2.48	7.01±0.37	油松、锐齿槲栎、华山松、毛樱桃、垂柳、 　木姜子、三桠乌药
5	1 457.32±13.14	1 208±355	13.80±0.84	0.5	32	2.78	5.69±0.26	锐齿槲栎、栗、油松、白桦、垂柳、 　榆树、桤木
13	1 757.57±20.17	1 254±207	13.80±1.19	0.6	29	2.68	6.55±0.29	毛樱桃、油松、锐齿槲栎、漆树、水蜡树、 　木姜子、灯台树
说明：毛樱桃Prunus tomentosa，垂柳Salix babylonica，三桠乌药Lindera obtusiloba，白桦Betula platyphylla，榆树Ulmus pumila，水蜡树Ligustrum obtusifolium。

下载: 导出CSV 

| 显示表格

1.2   采样设计

2023年7月，根据S型取样方法，在ck、5 a、13 a间伐样地内，用直径为3.6 cm的土钻采集0~10 cm的表层土样，为避免样品受到污染，将土壤混合储存于灭菌自封袋中，再用便携冷藏箱带回实验室。在室内充分混匀后过2 mm筛。一份新鲜土样于4 ℃冰箱保存，用于有效氮、土壤酶活性和土壤微生物生物量的测定；另一份土壤样品自然风干，用于其他土壤理化性质的测定。

1.3   测定指标及方法

1.3.1   土壤理化性质测定

土壤含水率采用105 ℃烘干法；土壤pH采用电位法(土水体积质量比为1.0∶2.5)；土壤总氮采用元素分析仪测定；土壤有机碳采用重铬酸钾氧化-外加热法；土壤有效氮指铵态氮和硝态氮的总和，分别采用2 mol·L⁻¹氯化钾浸提-靛酚蓝比色法、氯化钾提取-双波长紫外分光光度法测定；土壤总磷采用硫酸-高氯酸-钼锑抗比色法^[18]。微生物生物量碳、氮采用氯仿熏蒸法，使用岛津总有机碳分析仪测定。

1.3.2   土壤胞外酶活性及酶计量的测定与计算

参照SAIYA-CORK等^[19]的方法，测定与碳、氮、磷循环密切相关的酶活性，各种土壤酶的名称、简称及底物见表2。其中，水解酶(BG、BX、CBH、NAG、LAP、AcP)活性采用微孔板荧光法，用多功能酶标仪在365 nm波长处激发，450 nm波长处荧光测定；氧化酶(POX、PER)活性采用DOPA-紫外分光光度法，用多功能酶标仪在450 nm波长处测定。

表 2  土壤胞外酶的简称及所用底物

Table 2  Soil enzyme along with their enzyme abbreviation and substrate of soil enzyme

酶名称	简称	底物
β-葡糖苷酶β-glucosidase	BG	4-MUB-β-D-glucoside
β-木糖苷酶β-xalosidase	BX	4-MUB-β-D-xylopyranoside
纤维二糖水解酶Cellobiohydrolase	CBH	4-MUB-β-D-cellobioside
β-N-乙酰氨基葡萄糖苷酶β-N-acetylglucosaminidase	NAG	4-MUB-N-acetyl-β-D-glucosaminde
亮氨酸氨基肽酶Leucine aminopeptidase	LAP	L-leucine-7-amido-4 methylcounarin
酸性磷酸酶Acid phosphatase	AcP	4-MUB-phosphatase
酚氧化物酶Phenol oxidase	POX	L-dihydroxyphenylalanine(L-DOPA)
过氧化物酶Peroxidase	PER	L-dihydroxyphenylalanine(L-DOPA) and H2O2
说明：MUB为甲基伞形酮酰Methylumbelliferyl。

下载: 导出CSV 

| 显示表格

通过计算碳、氮和磷酶活性的比值研究土壤胞外酶化学计量^[20]，同时，采用酶计量的载体分析，即用矢量长度(V_L)及矢量角(V_A)分析间伐处理对微生物能量和营养的相对限制状况^[21]，计算公式如下。

式(1)~(5)中：$ {E}_\text{C/N} $、$ {E}_\text{C/P} $、$ {E}_\text{N/P} $分别为土壤碳获取酶/氮获取酶比值、土壤碳获取酶/磷获取酶比值、土壤氮获取酶/磷获取酶比值；$ {H}_{\mathrm{B}\mathrm{G}}\mathrm{、}{H}_{\text{NAG}}\mathrm{、}{H}_{\text{LAP}}\mathrm{、}{H}_{{\mathrm{Ac}}\mathrm{P}} $分别为BG、NAG、LAP、AcP的酶活性；SQRT为平方根函数，Degrees为角度转换函数，ATAN2为反正切函数。V_L越大，表明碳限制越严重。V_A以45°为分界线，＞45°为磷限制，＜45°为氮限制。偏离程度越大，限制程度越强。

1.4   数据分析

使用SPSS 25.0对不同间伐恢复时间下的土壤理化性质、胞外酶活性、酶化学计量比、酶矢量长度和角度的差异进行单因素方差分析(one-way ANOVA)和最小显著性差异法(LSD)(P＜0.05)；利用Sperman检验分析与土壤酶活性和酶矢量变化显著相关的土壤因子，利用Origin 2021绘图。以酶活性及其矢量作为物种因子，土壤理化性质作为环境因子，利用Canoco 5.0进行冗余分析。通过方差膨胀因子(VIF)判断解释变量之间的线性关系，剔除共线性较强(VIF＞5)的变量，对剩余的pH、有效氮、有机碳和全磷共4个变量进行研究。

2.   结果与分析

2.1   间伐恢复对土壤理化性质的影响

从表3可见：间伐恢复对土壤pH、有效氮、全磷、碳氮比、氮磷比、有机碳、微生物量碳、微生物量氮和微生物量碳氮比均有显著影响(P＜0.05)。恢复5 a的土壤pH显著高于ck (P＜0.05)。恢复13 a的土壤全磷、微生物量碳和微生物量氮均显著高于ck (P＜0.05)，分别是ck的1.28、1.19和1.15倍。土壤有效氮、碳氮比和氮磷比均显著低于ck (P＜0.05)。恢复5 a的土壤有机碳显著降低了25.93% (P＜0.05)，但恢复13 a的土壤有机碳质量分数逐渐恢复至未间伐前水平。间伐恢复对土壤含水率和全氮无显著影响。

表 3  不同间伐恢复时间下土壤理化特性状况

Table 3  Soil physical and chemical properties under different thinning treatments

间伐后恢复时间/a	pH	含水率/%	有效氮/(mg·kg−1)	全氮/(g·kg−1)	全磷/(g·kg−1)	碳氮比
ck	5.48±0.10 b	37.28±4.01 a	21.34±1.96 a	4.58±0.86 a	0.60±0.08 b	10.02±1.16 a
5	5.98±0.13 a	35.10±6.81 a	17.19±0.48 ab	3.28±0.68 a	0.52±0.10 b	9.34±1.41 ab
13	5.76±0.17 ab	40.37±1.67 a	16.56±0.58 b	3.93±0.44 a	0.77±0.07 a	8.55±1.32 b
间伐后恢复时间/a	氮磷比	有机碳/(g·kg−1)	微生物量碳/(g·kg−1)	微生物量氮/(g·kg−1)	微生物量碳氮比
ck	7.49±0.71 a	35.94±3.84 a	1.14±0.04 b	0.20±0.01 b	5.97±0.37 ab
5	6.45±0.95 ab	26.62±2.79 b	1.14±0.09 b	0.22±0.01 ab	5.09±0.13 b
13	5.04±0.34 b	33.33±2.27 ab	1.36±0.02 a	0.23±0.01 a	6.11±0.33 a
说明：数据均为平均值±标准误。不同小写字母表示不同处理间差异显著 (P＜0.05)。

下载: 导出CSV 

| 显示表格

2.2   间伐恢复对土壤酶活性及胞外酶计量比的影响

从图1可见：间伐恢复对不同土壤酶活性的影响并不一致。恢复13 a时土壤BX、AcP和NAG+LAP活性显著下降(P＜0.05)，较ck分别降低了25.39%、22.92%和46.25%，同时土壤BG活性还显著提高(P＜0.05)，是ck的1.34倍(P＜0.05)。土壤氧化酶(POX、PER)和CBH活性变化趋势与前4种酶不同，在恢复5 a时活性最低，在恢复13 a时活性最高。

图 1  间伐恢复对土壤酶活性的影响

Figure 1  Effect of thinning treatment on soil enzyme activity

下载: 全尺寸图片 幻灯片

通过矢量分析发现：V_A＞45°，且E_N/P＜1、E_C/N＞1 (图2A)，表明研究区土壤微生物生长代谢主要受碳和磷共同限制。森林土壤E_C/P和E_N/P显著偏离1，且随间伐后时间的持续而逐渐恢复或显著增大(P＜0.05，图2B)。V_A和V_L在3个间伐恢复间均有明显差异(图2C~D)。与ck相比，间伐恢复5 a的V_A显著降低了4.42%，13 a的V_L是ck的1.13倍(P＜0.05)。表明间伐措施在恢复初期能够缓解土壤微生物受磷限制的状况，而后随恢复时间的持续，微生物受碳限制程度显著增加(P＜0.05)。

图 2  间伐恢复对土壤酶化学计量及酶矢量的影响

Figure 2  Effects of thinning treatment on soil enzyme stoichiometry and enzyme vector

下载: 全尺寸图片 幻灯片

2.3   土壤酶整体变化和土壤理化性质的相关性分析及冗余分析

相关性分析(表4)表明：水解酶活性与有效氮、有机碳和微生物量碳氮比均呈正相关关系。其中土壤碳获取酶(BG、CBH)与土壤全磷、有机碳、微生物量碳呈显著(P＜0.05)或极显著(P＜0.01)正相关，BX活性与土壤有效氮、微生物量碳氮比呈显著正相关(P＜0.05)。土壤氮获取酶(NAG+LAP)和磷获取酶(AcP)均与土壤有效氮呈极显著正相关(P＜0.01)。酚氧化物酶(PER)除与pH呈显著负相关外(P＜0.05)，还与有机碳、微生物量碳氮比呈极显著正相关(P＜0.01)。V_A仅与pH呈极显著负相关(P＜0.01)。V_L与全磷和微生物量碳呈显著正相关外(P＜0.05)，还与氮磷比呈极显著负相关(P＜0.01)。

表 4  土壤酶变化与土壤理化性质的相关性分析

Table 4  Correlation analysis between soil enzyme changes and soil physical and chemical properties

指标	pH	IN	TP	SOC	MBC	MBC/MBN	N/P
POX	−0.54	−0.29	−0.04	−0.07	0.26	0.30	−0.04
PER	−0.65*	0.19	0.32	0.45*	0.22	0.52**	0.21
BG	0.28	0.35	0.73**	0.55**	0.63**	0.38	−0.25
BX	−0.53	0.54**	−0.01	0.27	0.10	0.45*	0.56
CBH	−0.01	0.24	0.46*	0.43*	0.53**	0.65**	0.17
AcP	−0.72*	0.57**	−0.38	0.06	−0.13	0.22	0.85**
NAG+LAP	0.17	0.66**	−0.08	0.14	−0.01	0.00	0.60
VA	−0.95**	0.01	−0.30	−0.06	−0.04	0.35	0.43
VL	0.45	−0.28	0.70**	0.31	0.48*	0.15	−0.63*
说明：IN为土壤有效氮，TP为土壤全磷，SOC为土壤有机碳，MBC为微生物量碳，MBN为微生物量氮，N/P为氮磷比。POX为酚氧化物酶，PER为过氧化物酶，BG为β-葡糖苷酶，BX为β-木糖苷酶，CBH为纤维二糖水解酶，AcP为酸性磷酸酶，NAG+LAP为氮获取酶(β-N-乙酰氨基葡萄糖苷酶和亮氨酸氨基肽酶总和)，VA为酶矢量角度，VL为酶矢量长度。*表示显著相关(P＜0.05)，**表示极显著相关(P＜0.01)。

下载: 导出CSV 

| 显示表格

冗余分析(图3)表明：剔除存在共线性关系的变量后，pH、有效氮、有机碳和全磷共解释了酶活性和酶矢量变异的73.71%。其中pH和有机碳是对土壤酶整体变化解释度最高的因子，分别解释了变量的48.80%和13.10%，且pH与酶指标变化显著相关(P＜0.05)。

图 3  土壤酶活性与理化性质关系的冗余分析

Figure 3  Redundancy analysis of soil enzyme activity and physical and chemical properties

下载: 全尺寸图片 幻灯片

3.   讨论

3.1   间伐恢复年限对土壤理化性质及微生物生物量的影响

间伐改变了秦岭松栎混交林表层土壤pH和养分质量分数，但在不同恢复阶段规律不一致。在本研究中，间伐导致pH提高，尤其是间伐恢复5 a后，这与许多学者的研究结果一致。如对云杉Picea crassifolia^[22]林和火炬松Pinus taeda^[23]林研究表明：间伐减少了针叶凋落物作为有机酸主要输入组分的产生，从而显著提高土壤表层pH。本研究中针叶树种的胸高断面积占比在间伐后有所降低，这在一定程度上能缓解土壤酸化。同时，间伐后土壤含水率、全氮、全磷和有机碳质量分数均呈先减少后逐步恢复的趋势。这可能是因为间伐短期内树冠层郁闭度减小，导致土壤蒸发增强的同时，也促进林下植被的快速生长，加快了土壤水分的消耗^[24]。凋落物作为土壤最主要的有机碳源，通过微生物转化为腐殖质的同时也改变了土壤pH，影响凋落物的分解，改变土壤养分水平^[25]。相较于ck，间伐恢复5、13 a后，凋落物量分别恢复至81.16%和93.41%，间伐恢复13 a的土壤全氮、全磷和有机碳质量分数有所提高，表明随时间的持续，林分结构及相关生态过程在一定程度上得到恢复。此外，本研究中微生物量碳、氮和土壤有效氮在间伐恢复13 a后的变化趋势不一致，可能因为间伐后林地内出现了栗、桤木和水蜡树等阳性植物，以及毛樱桃、白桦和漆树等阔叶树种，林地内相对多度增加，根系密度和根系分泌物增多，有利于土壤微生物生物量的积累^[26]。而林下喜光物种的快速生长^[27]，对土壤游离态氮的需求增大，导致土壤有效氮质量分数有所降低。这与周璇等^[28]对8年生柳杉Cryptomeria japonica人工林进行间伐后的研究结果一致。

3.2   间伐恢复年限对土壤酶活性的影响

在本研究中，间伐恢复年限导致土壤BX、AcP和NAG+LAP活性显著降低，但对其他土壤酶活性影响趋势不同，如POX、PER、BG和CBH通常在间伐恢复5 a时活性最低，在13 a时恢复到间伐前水平或高于未间伐处理(如BG)。这与其他研究结果相似，但并不完全一致^[29−30]。这种结果可能是由于不同的林分环境以及微生物利用资源多寡的差异，导致土壤酶活性对同一干扰方式的不同改变^[31]。随着间伐恢复时间的持续，易分解有机物质减少而难降解的碳相对较多^[32]，POX、PER和BG、CBH作为土壤中主要的木质素降解酶和纤维素降解酶，其活性得到显著提高，以增强微生物利用顽固性有机碳的能力。这与MEISAM等^[33]的研究结果一致。而以分解几丁质和蛋白质、半纤维素等易分解物质为主的NAG+LAP、BX活性的显著降低也映证了SINSABAUGH等^[34]的资源分配理论。

土壤胞外酶与土壤养分输入和微生物量等密切相关^[35]。通过相关分析发现：BG和CBH活性与微生物量碳、全磷显著正相关，表明土壤微生物数量的变化与碳循环土壤酶活性的变化关系最为密切，而全磷则是磷素限制环境中影响微生物生长的主要因素^{[7, 16]}。有效氮质量分数的减少虽然在一定程度上能促使氮获取酶的产生，但同样也会降低土壤微生物的活性和限制酶促反应底物供应，从而减少部分酶的释放^[36]，这与孙鹏跃等^[37]的研究结果一致。冗余分析发现：土壤pH也是影响土壤酶活性的主要因素，并与部分酶变化表现出负相关关系，这与多数研究结果是一致的^[3]。有研究表明：大多数土壤酶在特定的pH范围(最适值在4.0~5.5)内表现出最大的活性和稳定性，当pH超过这个范围时，酶活性会降低^[38]。

3.3   间伐恢复年限对酶化学计量比和微生物养分限制的影响

本研究中所有处理的土壤酶矢量角度均＞45°，符合亚热带地区森林土壤微生物更受磷素限制的理论^[39]。同时参与土壤碳、氮和磷循环相关酶计量比偏离了表层土壤中接近1∶1∶1的平均水平^[3]，也在一定程度上反映了秦岭区域松栎混交林间伐恢复过程中微生物受碳和磷的共同限制，这与薛悦等^[40]对安康市火池塘林区撂荒地恢复过程的研究结果相一致。与未间伐样地相比，间伐后恢复5 a时显著降低的酶矢量角度表征了微生物受到的磷限制减弱，随时间进程减弱效应逐渐消失，林内物种丰富度的提高和凋落物量的增加，促使土壤微生物分泌更多碳获取酶(如BG)来降解有机质，释放磷以供给微生物活动，以缓解磷限制，这些过程都会导致微生物碳限制的进一步增加。相关性分析结果中，酶矢量长度与微生物量碳呈显著正相关，证实了微生物需要更多的碳源来满足代谢活动所耗的能量，这与CUI等^[41]的研究结果相似。

4.   结论

间伐改变了松栎混交林区域内的年凋落物总量及针叶与阔叶的凋落量比例，同时改变了林内物种丰富度和林分郁闭度，从而影响了土壤基本理化性质。抚育间伐在一定程度上能够缓解土壤微生物受磷限制的状况，但随恢复时间持续，林内凋落物量逐渐增加使土壤微生物受碳限制更为严重。