林下植被去除对杉木人工林土壤酶活性及其化学计量比的影响

刘 仁，袁小兰，刘 俏，张春阳，郭春兰，陈伏生，王辉民，李建军＊

(1. 江西农业大学林学院，鄱阳湖流域森林生态系统保护与修复国家林业和草原局重点实验室，江西 南昌 330045；2. 江西农业大学南昌商学院，江西 南昌 332020；3. 中国科学院地理科学与资源研究所生态系统网络观测与模拟重点实验室，北京 100101)

林下植被可以增加人工林生态系统生物多样性和养分周转率，但其在维持森林生态系统结构和功能方面的作用通常被忽视[1]，特别是在传统森林管理方式中，常采用林下植被去除的方式来减少树木和下层群落之间的空间资源竞争[2]。然而，林下植被可通过改变土壤生物化学过程，如微生物活性和组成、碳和其他营养物质的周转，进而影响森林生态系统长期的生产力和稳定性[3]。以往研究表明，林下植被完整时，土壤养分、微生物量、酶活性和土壤呼吸均高于林下植被移除[4-5]。因此，在人工林的经营管理中，林下植被的保留比去除可能更有利于森林土壤地力的提高，但目前林下植被管理对土壤生化过程的影响发生机制仍不清楚，需要开展更多的研究工作。

土壤酶是由植物根系或微生物产生，对土壤碳（C）、氮（N）和磷（P）等养分变化具有一定的灵敏性，其常被用作衡量土壤质量状况的关键指标[6]。土壤酶活性可以反映底物利用率、微生物和植物对养分的需求，同时也在一定程度上影响植物和微生物在土壤环境变化时维持其养分平衡的策略[7]。Xu 等[8]通过研究我国东部南北样带沿线的9 个森林生态系统中的土壤酶活性及其化学计量比，结果发现，土壤酶化学计量比能很好的反映微生物的营养需求和环境中营养物质的供应状况。因此，开展土壤酶化学计量的研究为理解土壤中养分的周转、循环和平衡过程提供了新途径[9]。目前，有关林下植被去除的相关研究主要集中于土壤酶活性，且研究结果存在较大差异。有研究表明，林下植被去除显著降低了土壤C 含量和微生物生物量，同时也促使土壤水解酶活性降低[10]；而Liu 等[11]发现，去除林下植被对桉树人工林生态系统土壤酶活性无显著影响。因此，林下植被去除对人工林生态系统土壤酶活性的影响研究还需进一步探索，而且，从土壤酶化学计量的角度探究林下植被去除对土壤养分特征的相关研究相对较少，相关过程的研究可能对人工林养分维持及调控具有重要的指导意义。

杉木（Cunninghamia lanceolata(Lamb.) Hook.）广泛分布于我国南方，是主要经济树种之一，种植面积约占全国人工林面积的30%[12]。杉木能够高效吸收土壤养分进而决定了其速生特性[13]，而林下植被去除可能会进一步改变杉木人工林中土壤水分和养分有效性，最终影响土壤酶活性[10]，以往研究更多关注土壤酶活性变化，而从土壤酶化学计量比的角度探讨林下植被去除对土壤养分特征的影响研究较少。因此，本研究以养分较贫瘠的杉木人工林为研究对象，通过测定土壤C、N、P 含量及5 种C、N、P 水解酶活性，计算养分-酶化学计量比，并进一步分析土壤理化性质与酶活性及其化学计量比的相关性，进而揭示亚热带杉木人工林林下植被去除和保留的土壤养分和酶活性的响应特征，为林下植被管理对土壤质量影响的评估提供参考。

1 试验地概况

试验样地位于江西省吉安市泰和县中国科学院千烟洲石溪林场（26°44'52″ N，115°04'13″ E），该区地处亚热带季风气候区，雨量充沛，四季分明，海拔102 m，年平均气温17.9℃，年平均无霜期280 d，年平均降水量1 489 mm。本研究试验样地杉木活立木平均树高13.6 m，平均胸径11.4 cm，林分密度1 278 株·hm－2，郁闭度76%。林下灌木层平均高1.3 m，盖度11%，林龄15 a，坡向为东南方向，坡度为15°；主要物种有毛冬青（Ilex pubescensHook. et Arn.）、杨桐（Adinandra millettii(Hook. et Arn.) Benth. et Hook. f. ex Hance）、紫珠（Phyllostachys nigra(Lodd.) Munro） 和黄栀子（Gardenia jasminoidesEllis）等；林下草本层茂密，平均高度0.8 m，盖度80%，优势种为铁芒萁（Dicranopteris pedata(Houttuyn) Nakaike）和狗脊（Woodwardia japonica(L. F.) Sm.），伴生有淡竹叶（Lophatherum gracileBrongn.）、菝葜（Smilax chinaL.）和土茯苓（Smilax glabraRoxb.）等。

2 研究方法

2.1 试验设计

2013 年1 月，在江西省吉安市泰和县石溪林场杉木人工林内设置8 个15 m×15 m 样方，将样方地面上的所有凋落物移除（单独研究林下植被的影响），采用配对试验设计方法进行2 种处理：（1）林下植被保留（UP）：保留样方中的林下植被；（2）林下植被去除（UR）：在保留林下植被处理的样方附近划分与样方相同面积相近条件的杉木林，去除其中的林下植被。每个月对UP 和UR 样方内的凋落物进行人工徒手清除，以避免凋落物对处理结果的影响。

2.2 样品采集

2018 年7 月，在每个样方内采用 “S” 型布点取样，共计9 个布点，首先用三角铲收集50 cm×50 cm方格内矿质层以上的腐殖质层，然后再利用直径5 cm 的土钻在该区域采集0～5、5～10 cm 土层样品，将相同土层的样品混合后放入无菌自封袋内，做好标记，放入冰盒，迅速带回实验室。将土壤样品中的石砾、植物残根等杂物去除，研碎混匀，过2 mm 筛，将土样分为2 部分装入对应编号的无菌样品袋内，1 份放至4℃冰箱内冷藏保存用于测定土壤酶活性和有效氮（铵态氮（NH4+-N）、硝态氮（NO3－-N））；另1 份土样风干测C、N、P养分含量。

2.3 样品分析

2.3.1 土壤理化性质 土壤pH 值（土水比1:2.5）用pH 计测定；土壤含水量（SWC）用烘干法测定；有机碳(SOC) 用高锰酸钾-硫酸亚铁滴定法测定；全氮(TN) 和全磷(TP) 采用HClO4-H2SO4消煮法，然后在连续流动分析仪上测定；铵态氮(NH4+-N) 和硝态氮(NO3－-N) 采用2 mol·L－1KCl 作为浸提液（土液比1:10），在180 转的振荡器上振荡2 h，静置过滤后采用连续流动分析仪测定[14]。

2.3.2 土壤酶活性 用荧光标记底物法[15]测定5 种土壤酶（C 水解酶：β-1,4-葡萄糖苷酶（BG）、半纤维素酶（CB）；N 水解酶：β-1,4-N-乙酰基-葡糖胺糖苷酶（NAG）、亮氨酸氨基肽酶（LAP）、P 水解酶:酸性磷酸酶（AP））活性。

2.4 数据分析

利用SPSS 19.0 软件进行配对样本t 检验，分析林下植被保留和去除处理对土壤酶活性、酶化学计量比以及土壤理化性质的差异，显著性水平为0.05；采用Pearson 相关分析土壤酶活性及化学计量比与土壤理化性质之间的关系；利用Origin 2018 作图。

3 结果与分析

3.1 林下植被对杉木人工林土壤理化性质及养分化学计量比的影响

表1 表明：在不同土层中，UR 处理的SWC、SOC、 TN、 TP、 NH4+-N、 NO3－-N 均 显 著 低 于UP 处理，不同处理的pH 值差异不显著；除pH 值外，UR和UP 处理的理化性质均表现为腐殖质层含量最高，5～10 cm 土层最低；UR 处理的SOC、NO3－-N、TN、NH4+-N 含量和UP 处理的SOC、NO3－-N 含量均在3 个土层间差异显著；UP 处理腐殖质层的SWC、TN、NH4+-N 含量显著高于0～5、5～10 cm 土层，而0～5 cm 与5～10 cm 土层间差异不显著。

图1 表明：在不同土层中，UR 处理的C/N 和N/P 与UP 处理的差异不显著，而UR 处理的C/P在腐殖质层和0～5 cm 土层显著低于UP 处理；UR 和UP 处理的C/P 和N/P 均表现为腐殖质层含量最高，5～10 cm 土层最低，但C/N 却与之相反；UR 和UP 处理的N/P 在3 个土层间差异显著，C/N 仅在UR 处理的3 个土层间差异显著，C/P 仅在UP 处理的3 个土层间差异显著。

3.2 林下植被对杉木人工林土壤酶活性及其化学计量比的影响

图2 表明：在不同土层中，UR 处理的5 种酶活性均低于UP 处理。在腐殖质层中，UP 处理的LAP 活性显著高于UR 处理，增幅为48%；在0～5 cm 土层中，UP 处理的CB、NAG 和LAP 活性显著高于UR 处理，增幅分别为38%、53% 和32%；在5～10 cm 土层中，UR 处理的CB 和LAP 活性显著低于UP 处理，降幅分别为45% 和18%。在UR 处理中，腐殖质层的AP、BG、NAG 活性显著高于0～5 cm 和5～10 cm 土层，但0～5 cm 与5～10 cm 间差异不显著，而CB 和LAP 活性在3 个土层间均差异显著；在UP 处理中，腐殖质层的AP、CB、NAG 活性显著高于0～5 cm 和5～10 cm土层，但0～5 cm 与5～10 cm 之间差异不显著，而BG 和LAP 活性在3 个土层间均差异显著。

表1 不同林下植被处理对杉木人工林土壤理化性质的影响 Table 1 Effects of different understory vegetation treatments on physical and chemical properties of Chinese fir plantation soil

图1 不同林下植被处理对杉木人工林土壤养分化学计量比的影响Fig. 1 Influences of different understory vegetation treatments on soil nutrient stoichiometric ratio of Chinese fir plantation

图2 不同林下植被处理对杉木人工林土壤酶活性的影响Fig. 2 Effects of different understory vegetation treatments on soil enzyme activity of Chinese fir plantation

图3 表明：在3 个土层中，UP 处理与UR 处理的ln(BG+CB)/ln(AP) 和ln(NAG+LAP)/ln(AP) 差异不显著；UP 处理的ln(BG+CB)/ln(NAG+LAP) 仅在0～5 cm 土层中显著低于UR 处理，降幅为16%；UR 处理的腐殖质层的ln(BG+CB)/ln(AP) 和ln(NAG+LAP)/ln(AP) 显著高于0～5 cm 和5～10 cm土层，但0～5 cm 与5～10 cm 间差异不显著；而UP 处理的ln(NAG+LAP)/ln(AP) 在3 个土层间均差异显著，腐殖质层的ln(BG+CB)/ln(AP) 显著高于0～5 cm 和5～10 cm 土层，但UR 和UP 处理的ln(BG+CB)/ln(NAG+LAP) 在腐殖质层显著低于5～10 cm 土层，降幅分别为20% 和23%。

图3 不同林下植被处理对杉木人工林土壤酶化学计量比的影响Fig. 3 Influences of different understory vegetation treatments on soil enzyme stoichiometric ratio of Chinese fir plantation

3.3 酶活性及其化学计量比与土壤理化性质的相关分析

图4 表明：杉木人工林土壤的C、N、P 水解酶活性之间呈极显著正相关。土壤酶活性及化学计量比与土壤理化性质之间的相关分析（表2）表明：AP、BG、CB、NAG、LAP 5 种酶活性及ln(BG+CB)/ln(AP)、ln(NAG+LAP)/ln(AP) 与SWC、SOC、 TN、 TP、 NH4+-N、 NO3－-N、 C/P、 N/P 呈显著或极显著正相关，与C/N 呈极显著负相关；土壤ln(BG+CB)/ln(NAG+LAP) 与C/N 呈极显著正相关，但与SWC、SOC、TN、TP、NH4+-N、NO3－-N、C/P、N/P 呈极显著负相关(P＜0.05)；土壤酶活性及其化学计量比与pH 相关不显著。

图4 不同林下植被处理的杉木人工林土壤酶化学计量关系Fig. 4 Stoichiometric relationship of soil enzyme in Chinese fir plantation under different understory vegetation treatments

表2 土壤酶活性和酶化学计量比与土壤理化性质之间的 Pearson 相关分析 Table 2 Pearson correlation analysis between soil enzyme activity and enzyme stoichiometric ratio and soil physical and chemical properties

4 讨论

4.1 林下植被去除对土壤特性和酶活性的影响

本研究表明，林下植被去除降低土壤0～5、5～10 cm 土层的SWC、SOC、TN、TP、有效氮（NH4+-N、NO3－-N）养分含量和CB、LAP 活性，该结果与Yang 等[10]的研究结果一致。就土壤C 水解酶而言，可能是植被去除后，增加了土壤水分的散失量[16]和降低外源C 养分输入量[17]，进而降低酶底物的扩散速率和酶合成，最终导致土壤CB 活性降低[18]。另外，有研究发现，植物根系释放丰富的有效养分，促进微生物分泌更多N 水解酶参与到N 循环中[19-20]，因此，本研究中林下植被去除显著降低了0～5、5～10 cm 土层的有效氮含量和LAP 水解酶活性；但林下植被去除并未对AP 活性产生显著影响，不同土层也产生一致响应，这可能与该试验区受到P 限制有关。我国南方土壤P 通常以有机形式存在，容易被铁铝物质固定成络合物，难以利用[21]，土壤微生物通过产生更多的AP 来矿化有机磷，以满足其对P 的需求[22]，而且杉木根系与真菌共生形成的菌根真菌可以产生土壤AP 分解土壤有机P[23-24]，因此，林下植被去除与否并不影响AP 活性。在缺乏外源养分对土壤进行补充情况下，微生物会持续消耗土壤中原有的磷含量，而AP 活性保持较高水平，最终导致林下植被去除后土壤中TP 含量降低。总之，林下植被去除通过降低土壤水分和养分含量，最终使土壤CB、LAP 活性降低。

4.2 土壤酶化学计量比与P 限制的关系

本研究发现，在林下植被保留和去除处理下，0～5 cm 和5～10 cm 土层中，ln(BG+CB)/ln(AP) 和ln(NAG+LAP)/ln(AP) 均显著低于全球尺度上土壤酶活性比ln(BG+CB)/ln(AP) (0.62) 和ln(NAG+LAP)/ln(AP) (0.44)[25]，这表明P 限制发生在亚热带杉木人工林，该结果与乔航等[26]研究结果一致。这是因为土壤微生物分泌酶的种类需要在C、N、P 的获取上做出权衡，以响应土壤基质中养分资源的变化[27]。土壤养分含量可以通过影响土壤有效养分浓度和土壤C、N、P 化学计量，最终对土壤酶化学计量产生较大影响，说明土壤酶化学计量与土壤养分化学计量之间存在密切关系[8]。在本研究中，林下植被保留和去除的ln(BG+CB)/ln(NAG+LAP) 与C/N 均表现为腐殖质层显著低于5～10 cm，而ln(BG+CB)/ln(AP) 与C/P 及ln(NAG+LAP)/ln(AP) 与N/P 表现为腐殖质层显著高于5～10 cm，这可能是5～10 cm 土层中的土壤有机质比腐殖质层的含量低，腐殖化更高，P 含量在3 个土层中变化程度不大所致。因此，土壤酶化学计量比可作为表征当前土壤养分有效性状况的重要指标。

4.3 影响土壤酶活性及其化学计量比的因子

土壤酶活性及其化学计量的变化与土壤水分、养分及其化学计量比变化密切相关。本研究发现，SWC与C、N、P 酶活性及ln(BG+CB)/ln(AP)、ln(NAG+LAP)/ln(AP) 呈极显著正相关，但与ln(BG+CB)/ln(NAG+LAP) 呈极显著负相关，可能是较高的SWC促进微生物的活性和功能，直接刺激C、N、P 酶的分泌[28-29]，但AP 酶活性变化程度低于C 和N 酶活性的变化程度[30]。研究发现，土壤酶活性及其化学计量比可作为微生物对土壤养分有效性和养分需求的检测指标[31]。在本研究中，C、N、P 水解酶活性及其化学计量比与土壤C、N、P 和有效氮养分含量呈显著或极显著相关，这是由于土壤易分解有机物（如蛋白质、核酸和纤维素）主要被土壤酶分解成有效养分以供给微生物和植物利用[32-33]，进而导致C、N、P 及有效氮养分与土壤酶活性产生极密切关系。另外，ln(BG+CB)/ln(NAG+LAP) 与C/N、 ln(BG+CB)/ln(AP) 与C/P 及ln(NAG+LAP)/ln(AP) 与N/P 均呈极显著正相关。基于此，本研究的土壤水分、养分及其化学计量比可能是影响土壤酶活性及其化学计量比的重要因子。

5 结论

林下植被去除显著降低土壤SWC、SOC、TN、TP、有效氮（NH4+-N、NO3－-N）养分含量和抑制了CB 和LAP 活性，但对AP 活性的影响不显著。通过分析土壤理化性质与酶活性及其化学计量比的相关性，土壤SWC、SOC、TN、TP、有效氮含量和养分化学计量比与C、N、P 酶活性及其化学计量比关系极为密切。因此，土壤水分、养分及其化学计量比可能是影响土壤酶活性及其化学计量比的重要因子。