刘宗悦，吕时新，徐钧杰，陈 旭，楼一恺,3，祁祥斌，余树全

（1. 浙江农林大学 林业与生物技术学院，浙江 杭州 311300；2. 永康市自然资源和规划局，浙江 金华321300；3. 国家林业和草原局 华东调查规划设计院，浙江 杭州 311201；4. 浙江天目山国家级自然保护区管理局，浙江 杭州 311311）

森林土壤的物理、化学及生物学性质综合表现为土壤肥力。土壤肥力的高低影响着森林植物的分布、生长以及群落的组成、结构和生产力，植物又通过凋落物、根系分泌物等反作用于土壤，改变着土壤的物理、化学及生物学性质，由此形成一个动态的相互作用机制，提升森林生态系统结构与功能[1−3]

毛竹Phyllostachys edulis属于禾本科Gramineae刚竹属Phyllostachys植物，原产于中国亚热带地区，以地下竹鞭进行无性繁殖，具有扩展蔓延速度快和竞争力强的特点。从20世纪90年代开始，毛竹因其具有良好的竹材、竹芛等多种用途而作为笋材两用的经济树种在中国亚热带地区大面积种植，为山区农民脱贫致富做出贡献[4]。随着中国经济的发展，毛竹林经营成本不断提高，而竹材价格反而降低，因此，很多竹农已经放弃了对毛竹林的经营管理，导致失去管理的毛竹林无序扩张蔓延。近年来，这种趋势不断加剧。在许多地方由于毛竹地下竹鞭扩张渗透到阔叶林等相邻森林中，迅速成林并侵占其他物种的生境，形成毛竹单优势群落，导致植物种类减少，群落层次结构简化，生物多样性显著降低，土壤物理性质、化学性质和土壤微生物发生显著变化，整个区域生态系统退化[5−9]，威胁到生物多样性保护和区域生态安全。因此，控制毛竹林无序蔓延扩张，开展毛竹林皆伐后生态恢复效果研究，已受到全社会广泛关注。2015年，为控制和减少毛竹入侵对天目山国家级自然保护区生态破坏，经审批，对该区域部分入侵毛竹林进行皆伐。本研究在毛竹林皆伐区和未伐区设置样地进行比较研究，探讨毛竹林皆伐后土壤自然恢复效果，为天目山毛竹林入侵区域生态修复提供技术支持，为亚热带毛竹林入侵区生态修复提供借鉴。

1 研究方法与数据来源

1.1 研究区概况

研究区位于浙江省杭州市临安区天目山国家级自然保护区，地处中亚热带北缘，浙皖两省交界处，是南北植物汇流之区，也是中国最主要的野生植物类型自然保护区之一，面积为4 300 hm2，海拔300 ～1 556 m。研究区土壤类型随着海拔变化具有较大差异，600 m以下为红壤，600 ～850 m主要为黄红壤，850 ～1 200 m为黄壤，1 200 m以上为黄棕壤。植被分布也有着明显的垂直带谱特征：在海拔230 ～850 m分布常绿阔叶林，850 ～1 100 m分布常绿落叶阔叶混交林，1 100 ～1 350 m分布落叶阔叶林，1 400 m以上为落叶矮林。

在天目山国家级自然保护区，毛竹林作为一种特殊的森林类型，主要分布在海拔350 ～900 m。由于保护区不允许开展森林经营活动，毛竹林完全处于自然状态，立竹度高，并不断向周边蔓延扩张，面积从55 hm2蔓延扩张到105 hm2，对保护区原有植被蚕食的现象日益严重。为控制和减少毛竹入侵对天目山国家级自然保护区的生态破坏，促进生物多样性保护。2015年经上级批准，对部分入侵毛竹林进行皆伐作业，并于皆伐后连续3 a进行人工除笋，防止毛竹复占。毛竹林皆伐后经过5 a自然恢复，物种组成有了明显的变化，且物种种类明显增加，与未皆伐毛竹林相比植物多样性有显著增加[10−11]。

1.2 试验设计

于2017年在毛竹林皆伐区域设置了清除(CR)和保留(UR)采伐剩余物样地各3个，在未伐区毛竹林内设置对照样地(ck) 3个。样地面积20 m×20 m，在每个样地四角和中间分别设置1个5 m×5 m样方。

在2020年11月样地复查的同时采集土壤样品，每个样地内采用5点取样法，均匀采取0 ～10、10 ～20 cm土层中3个点位土壤的混合样品，并过2 mm筛，去除样品中的植物碎屑及石块，经过研磨过筛处理后测定土壤化学性质；使用200 cm3环刀于样方中心采集0 ～10、10 ～20 cm土层土壤用于测定土壤物理性质；新鲜土壤样品保存在4 ℃环境中用于测定土壤胞外酶活性。

1.3 测定方法

土壤理化性质参考《土壤农业化学分析方法》进行测定[12]。土壤容重、孔隙度、最大持水量、最小持水量、毛管持水量及含水量采用环刀法测定。土壤pH采用pH计测定(水土质量比为2.5∶1.0)；土壤有机质采用高温外热重铬酸钾氧化-容量法测定；全氮用凯氏定氮法测定；全磷采用氢氧化钠熔融法-钼锑抗比色测定；全钾采用氢氧化钠熔融法测定；碱解氮采用碱解扩散法测定；有效磷采用盐酸-氟化铵浸提，钼锑抗比色法测定；速效钾采用火焰光度计测定。土壤胞外酶活性参考GERMAN等[13]使用鲜土进行测定，脲酶采用比色法测定，水解酶活性使用甲基伞形酮(4-MUB)作为底物标示，氧化酶活性使用L-二羟苯丙氨酸(L-DOPA)作为底物标示。利用微孔板荧光法、使用多功能酶标仪测定单位时间内特定底物被各种酶催化释放出的荧光产物的量来计算酶活性。

1.4 数据分析

采用 Excel 2019 进行数据整理；用 SPSS 26.0 进行统计分析。对毛竹林样地和砍伐后自然恢复迹地土壤物理、化学性质及胞外酶活性差异做单因素方差分析(one-way ANOVA)，显著性水平设置为0.05，并采用Tukey 方法进行显著性检验。统计结果用Origin 2018制图。

1.5 综合恢复效果评价

1.5.1 构建评价指标 综合考虑土壤物理性质、化学性质及胞外酶活性，本研究选取了34 个候选指标，分别是：0 ～10、10 ～20 cm土层土壤容重、孔隙度、毛管持水量、田间持水量、饱和持水量、pH、有机质质量分数、全氮质量分数、全磷质量分数、全钾质量分数、碱解氮质量分数、有效磷质量分数、速效钾质量分数、过氧化物酶活性、脲酶活性、酸性磷酸酶和β-葡萄糖苷酶活性。在运用主成分分析时，首先，根据特征值不小于 1 和累积方差贡献率不小于 80%，提取主成分。然后，在每个主成分中，选出不小于最大载荷值 90%的指标，对其进行 Pearson 相关分析，将相关系数之和最大的指标作为评价指标；若相关系数之和最大的指标不止1个，则选取载荷值较大的指标；若不小于最大载荷值 90%的指标只有1个，则该指标直接作为评价指标。

1.5.2 指标权重 根据主成分分析得到所选取指标的公因子方差，计算各指标的权重(Wi)：其中：n是综合评价的指标个数，Ci是第i个指标的公因子方差。

1.5.3 综合评价模型 基于模糊数学模型，计算土壤性质的综合得分(S) ：其中：Fi是指标i的离差标准化后的值。S越大，表示毛竹林皆伐后自然恢复的土壤综合恢复效果越好，反之则越差。

2 结果与分析

2.1 土壤物理性质

由表1可见：在0 ～10 cm表层土壤中，ck处理的土壤容重分别比CR、UR处理高31%和 14%；UR处理土壤容重略高于CR处理, 但无显著差异；UR、CR处理的土壤总孔隙度、毛管持水量、田间持水量和饱和持水量均高于ck处理，但无显著差异。在10 ～20 cm土壤中，ck处理土壤容重分别比CR、UR处理高39%和37%；ck、CR和UR处理的总孔隙度无显著差异；UR处理的毛管持水量、田间持水量和饱和持水量均高于CR和ck处理(P＜0.05)。

表1 不同处理的土壤物理性质Table 1 Soil physical properties of different treatments

2.2 土壤化学性质

由表2可见：在0 ～10 cm土壤中，UR处理的pH显著低于 CR和ck处理(P＜0.05)，UR处理的有机碳、全氮、全磷和碱解氮质量分数均显著高于CR和ck处理，UR处理的有机碳质量分数比CR和ck处理高99%和117%，全氮高93%和123%，全磷高57%和87%，碱解氮高95%和107%。ck处理的速效钾质量分数显著低于UR和CR(P＜0.05)。在10 ～20 cm土壤中，UR处理的pH略低于CR和ck处理，UR处理的有机碳、全氮、全磷和碱解氮质量分数均显著高于CR和ck处理(P＜0.05)；UR处理的全钾质量分数显著低于CR和ck处理，ck处理的速效钾质量分数显著低于UR和CR处理 (P＜0.05)。

表2 不同处理的土壤化学性质Table 2 Soil chemical properties of different treatments

2.3 土壤胞外酶活性

由表3可知：整体来看，除土壤酸性磷酸酶外，毛竹林皆伐后UR、CR处理的土壤脲酶、β-葡萄糖苷酶和过氧化物酶活性高于ck处理；UR处理土壤的4种胞外酶活性均高于CR处理，幅度为46% ～98%。在0 ～10 cm土壤中，UR处理的脲酶、酸性磷酸酶和β-葡萄糖苷酶活性均显著高于CR和ck处理(P＜0.05)，3种处理的过氧化物酶活性无显著差异；ck处理土壤的酸性磷酸酶活性比CR处理高47%。在10 ～20 cm土壤中，UR处理的脲酶、酸性磷酸酶和β-葡萄糖苷酶活性均显著高于CR和ck处理(P＜0.05)，3种处理的过氧化物酶活性无显著差异。

表3 不同处理的土壤胞外酶活性Table 3 Soil extracellar enzyme activity of different treatments

2.4 综合恢复效应评价

通过表层土壤34个指标的主成分分析，提取7个主成分，累积贡献率为83.82%；第1主成分特征值为14.100，贡献率41.5%，主要影响因子为0 ～10 cm 全氮(xTNⅠ)、10 ～20 cm全氮(xTNⅡ)、0 ～10 cm 有机质(xSOCⅠ)、10 ～20 cm有机质(xSOCⅡ)、0 ～10 cm碱解氮(xANⅠ)、10 ～20 cm碱解氮(xANⅡ)等6个指标；第2主成分特征值为4.646，贡献率为13.7%，主要影响因子为0 ～10 cm孔隙度(xPⅠ)、10 ～20 cm孔隙度(xPⅡ)和10 ～20 cm饱和持水量(xShcⅡ)等3个指标；第3主成分特征值为3.100，贡献率为9.1%，主要影响因子为0 ～10 cm速效钾(xAKⅠ)和10 ～20 cm速效钾(xAKⅡ)；第4主成分特征值为2.500，贡献率为7.3%，主要影响因子为0 ～10 cm有效磷(xAPⅠ)、10 ～20 cm有效磷(xAPⅡ)和10 ～20 cm β-葡萄糖苷酶活性(xBDⅡ)；第5主成分特征值为1.610，贡献率为4.735%，主要影响因子为0 ～10 cm孔隙度(xPⅠ)和10 ～20 cm过氧化物酶活性(xPERⅡ)；第6主成分特征值为1.300，贡献率为3.8%，主要影响因子为10 ～20 cm β-葡萄糖苷酶活性(xBDⅡ)；第7主成分特征值为1.300，贡献率为3.7%，主要影响因子为10 ～20 cm酸性磷酸酶活性(xACPⅡ)。各主成分影响因子经相关分析后最终选取7个评价指标，分别为xTNⅠ、xPⅡ、xAKⅠ、xAPⅠ、xPERⅡ、xBGⅡ、xACPⅡ。

根据选定指标的离差标准化值和指标权重(表4)构建土壤性质综合评价模型：

表4 综合评价指标及其权重Table 4 Comprehensive evaluation indices and weights

从0 ～10、10 ～20 cm等2层土壤综合评分结果来看(图1)，毛竹林皆伐后恢复迹地得分均高于毛竹林得分，且UR处理得分高于CR处理。CR、UR、ck处理综合得分分别为0.41、0.61、0.37，UR处理综合得分分别比CR和ck处理高49%和 65%。UR处理的0 ～10 cm全氮得分为0.11，高于CR处理(0.05)和ck处理(0.04)；CR处理和UR处理的0 ～10 cm速效钾得分均为0.10，高于ck处理(0.04)；UR处理的10 ～20 cm β-葡萄糖苷酶得分为0.08，高于CR和ck处理(0.04)；UR处理的10 ～20 cm酸性磷酸酶得分为0.10，高于CR (0.05)和ck处理(0.06)。

图1 各样区的土壤性质综合得分Figure 1 Synthesis scores of soil properties in the sample areas

3 结论与讨论

3.1 毛竹林皆伐后植物多样性恢复对土壤的影响

自然条件下植物和土壤有着最直接的关联，具有密切的相互反馈作用，植物群落的变化会影响其生长环境中的土壤，最终改变土壤结构和土壤养分的有效性[14]。毛竹的入侵性极强，可通过其独特的竞争策略影响入侵林分其他物种的生长，改变入侵地的群落结构[11]，导致植物物种多样性明显下降直至形成毛竹纯林，相应的土壤也会发生变化[15]。本研究结果表明：毛竹林皆伐后经5 a自然恢复，林地土壤与ck相比，土壤容重极显著降低，孔隙度和持水能力增加，土壤有机碳、全氮、全磷、碱解氮和速效钾质量分数提高，土壤胞外酶活性增强。主要原因是毛竹林皆伐后大量物种进入，植物多样性明显提高。这一结果与毛竹入侵相邻森林引起的植物多样性下降[16]、土壤退化是相反的，因此，通过毛竹林皆伐控制其快速扩张，可作为目前毛竹扩张蔓延区域的生态修复措施。

3.2 毛竹林皆伐后采伐剩余物管理方式对土壤的影响

毛竹林采伐后会产生大量的采伐剩余物，处理方式不同会对迹地更新方式和生态修复产生影响。保留采伐剩余物覆盖在地表，腐烂分解后可增加土壤有机质和养分，减少土壤侵蚀，但会对植物种子萌发、幼苗生长和植被恢复产生一定的影响[17−18]。本研究表明：在0 ～10和10 ～20 cm土层土壤中，采伐剩余物保留样区的土壤有机碳、全氮、全磷、碱解氮及速效钾质量分数均显著高于采伐剩余物清理样区。这与池鑫晨等[15]、黄启堂等[19]、ZHAO等[20]、吴家森等[21]的研究结果一致，且采伐剩余物保留样区的脲酶、酸性磷酸酶、β-葡萄糖苷酶活性显著高于采伐剩余物清理样区。有研究表明：土壤酶活性与土壤有机质含量密切相关[22]，土壤有机质含量高可以维持较高的土壤胞外酶活性[23]，保留采伐剩余物后提高了土壤中的有机质。因此，建议在利用毛竹林皆伐手段进行区域生态修复时，应尽量保留采伐剩余物，但应让其分散均匀，不宜局部过厚堆积，减少对植被恢复影响。

3.3 结论

毛竹林皆伐后自然恢复迹地土壤恢复良好，保留采伐剩余物的管理方式相较于清除采伐剩余物的管理方式可以更有效地提高土壤有机碳、全氮、全磷、碱解氮和速效钾质量分数，以及土壤胞外酶活性。在未来人工毛竹林弃管后实施生态修复时，保留采伐剩余物可以节省更多人力及资金成本，同时也可以达成更好的土壤恢复效果。