梯度有机肥替代重塑茶园土壤微生态与微生物碳利用特征

《Agronomy》：Gradient Organic Fertilizer Substitution Reshapes Soil Microecology and Microbial Carbon Utilization in Tea Plantations

　　

有机肥替代被广泛用于改善酸化的茶园土壤，但替代梯度对土壤团聚体、微生物碳利用以及微生物C-N-P限制的耦合效应仍未得到充分量化。本研究旨在确定梯度有机肥替代如何调控亚热带茶园土壤理化性质、团聚体粒径分布、胞外酶化学计量以及微生物碳源利用，并在实验条件下确定适宜的替代范围。研究人员在中国福建省武夷山酸性红壤茶园开展了一项长期田间试验。设置了五个处理：纯化肥（CF）以及有机肥分别替代33%（OF33）、50%（OF50）、67%（OF67）和100%（OF100）的化肥氮素。测定了茶树生长指标、土壤理化性质、团聚体机械稳定性、微生物生物量碳、氮、磷（MBC、MBN、MBP）、胞外酶活性以及BIOLOG-ECO碳源利用图谱。结果表明，相对于CF，OF33和OF67使土壤pH提高了0.46–0.60个单位。OF33总体上改善了速效氮和钾含量以及团聚体稳定性指标，而OF50增加了百芽重和微生物碳源丰富度。OF100增加了MBC，但也增加了向量长度和向量角度，表明相对更强的碳和磷限制。BIOLOG-ECO图谱将各处理区分开来，适度替代（OF33–OF50）拓宽了碳源利用生态位。然而，不同指标在不同处理下达到峰值；因此，在本实验条件下，33–50%的替代应被视为一个适宜的范围，而非普遍最优比例。这些发现为红壤茶园减少化肥施用提供了特定地点和品种的依据。

研究背景与问题

亚热带红壤区茶园长期过量单一施用氮磷钾化肥，导致土壤严重酸化（超过65%的山地茶园土壤pH低于4.5）、活性铝离子大量淋溶抑制茶树侧根发育与养分转运蛋白表达、活性有机碳持续耗竭、碳氮磷元素比例失衡、根际微生物群落多样性下降及养分转化功能减弱，最终表现为茶叶百芽重和叶绿素含量降低，同时面源氮磷流失风险上升。有机肥部分替代化肥是缓解茶园土壤酸化、重建土壤碳库和恢复根际微生态的核心生态改土技术，但现有研究多仅设置两个施肥水平（纯化肥vs.全有机肥），缺乏对有机肥替代梯度下土壤理化性质、团聚体结构与微生物碳代谢差异调控效应的系统性整合分析，也极少将土壤团聚体结构、理化养分、生态酶学指标和微生物碳源利用纳入同一分析框架。此外，多数研究仅关注土壤层面变化，缺乏农艺指标与土壤微生态系统关系的定量分析。

研究内容与意义

研究人员利用福建武夷山酸性红壤茶园长周期定位施肥试验平台，设置了纯化肥（CF）以及有机肥替代33%（OF33）、50%（OF50）、67%（OF67）和100%（OF100）化肥氮素的五个梯度处理，系统测定了茶树生长指标、土壤理化性质、团聚体机械稳定性、微生物生物量碳氮磷（MBC、MBN、MBP）、胞外酶活性及BIOLOG-ECO碳源利用图谱。研究揭示了梯度有机氮替代对茶园土壤微生态系统的差异化调控规律，确定了实验条件下的适宜替代范围为33–50%，为红壤茶园减施化肥和发展微生物驱动的绿色栽培体系提供了理论依据。

主要关键技术方法

1. 1.
   土壤团聚体分级与稳定性计算：采用干筛法对风干土壤进行粒径分级（2 mm、1 mm、0.5 mm、0.25 mm），计算>0.25 mm团聚体质量占比（R_0.25）、平均重量直径（MWD）和几何平均直径（GMD）。
2. 2.
   胞外酶活性测定与生态酶化学计量：测定β-葡萄糖苷酶（BG）、N-乙酰-β-D-氨基葡萄糖苷酶（NAG）、亮氨酸氨肽酶（LAP）和酸性磷酸酶（ACP）活性，计算lnBG/lnACP和lnBG/ln(NAG+LAP)比值，以及微生物营养限制向量的长度（Vector L）和角度（Vector A）。
3. 3.
   BIOLOG-ECO碳源利用指纹分析：将根际土壤悬浮液接种至BIOLOG-ECO微孔板，于28℃恒温黑暗培养7天，测定590 nm吸光度，计算平均颜色变化率（AWCD）、碳源丰富度、香农多样性指数（Shannon H'）、均匀度和Levins生态位宽度（Levins B），并通过主坐标分析（PCoA）和随机森林模型识别特征性差异碳源。
4. 4.
   冗余分析（RDA）与相关性分析：使用RDA分析茶树生长指标、土壤理化性质和胞外酶活性、土壤团聚体稳定性指标与微生物碳源利用模式之间的关联，结合Mantel检验和Pearson/Spearman相关性矩阵揭示多变量间的相互关系。

研究结果

3.1 不同施肥方式对茶树生长的影响

有机肥替代比例对茶树叶片叶绿素含量（SPAD值）和百芽重具有显著的差异化调控作用。OF67处理的叶绿素含量最高，其次为OF33，两者均显著高于OF50、CF和OF100处理。百芽重在OF50处理中最高，其次为OF33，均显著高于CF处理，而高替代比例处理（OF67和OF100）的百芽重显著低于CF。芽密度在CF处理中最高，OF33最低。总体而言，中等有机替代（33–50%）与最高的叶绿素含量和百芽重相关联。

3.2 不同施肥方式对茶树根际土壤理化性质的影响

所有有机替代处理的土壤pH均显著高于纯化肥对照，其中OF33和OF67效果最佳，pH提升0.46–0.60个单位。土壤有机质（OM）随替代比例增加呈稳步上升趋势，OF100较CF提高35.2%。全氮（TN）在OF33处理中最高，较CF增加57.7%。全钾（TK）和全磷（TP）在有机替代处理中总体下降，其中TP降幅显著（OF50较CF降低72.8%）。速效氮（AN）和速效钾（AK）在OF33处理中最高，分别较CF提高17.9%和17.1%；速效磷（AP）在CF中最高，OF67最低（较CF降低32.9%）。67%替代比例导致所有速效养分指标全面下降。

3.3 不同施肥方式对茶树土壤酶及其化学计量比的影响

碳获取酶BG活性在OF100中最高，显著高于CF和OF67。磷获取酶ACP活性在所有有机替代处理中均显著高于CF，OF33和OF67活性最高。氮获取酶（NAG+LAP）活性在CF中最高，随替代比例增加总体呈下降趋势。生态酶化学计量显示，向量长度（Vector L）在OF100中显著高于其他处理，表明微生物碳限制最严重；OF67的向量长度最低。所有处理的向量角度（Vector A）均大于45°，指示相对磷限制；OF100和OF33的向量角度最高，CF最低。线性回归显示向量长度与向量角度呈显著正相关，表明碳和磷限制协同加剧。

3.4 不同施肥方式下茶树根际土壤微生物生物量C、N、P化学计量特征

MBC随有机肥施用量增加而持续上升，OF100处理MBC是CF的3.5倍。MBN在CF中最高，有机替代后显著下降，OF33最低。MBP在OF33中最高，显著高于其他处理，CF最低。微生物生物量碳氮比（MBC:MBN）和碳磷比（MBC:MBP）均随替代比例增加而持续升高，表明有机肥输入主要扩张了微生物碳库，而磷库增长滞后。

3.5 不同施肥方式下茶树根际土壤BIOLOG分析

AWCD在OF50和OF100中最高，OF67最低。碳源丰富度在OF50中显著最高。香农多样性和均匀度在OF33、OF50和OF100中较高，OF67最低。Levins生态位宽度证实中等有机替代（OF33、OF50）显著拓宽了微生物碳源利用生态位。PCoA分析显示各处理碳源利用模式分离显著，PERMANOVA确认组间差异显著。随机森林筛选出L-天冬酰胺、L-精氨酸、D-苹果酸、γ-羟基丁酸、吐温40、吐温80、D-甘露醇和N-乙酰-D-葡萄糖胺等为核心差异碳源。中等替代处理微生物偏好利用氨基酸和小分子有机酸，OF100偏好多糖底物，CF偏好单糖和酚酸。

3.6 团聚体机械稳定性

OF33处理的大团聚体（>2 mm）质量分数、R_0.25、MWD和GMD均为最优，其次为OF50。所有有机改良处理中微团聚体（<0.25 mm）的平均比例均低于纯化肥处理。稳定性指标排序一致为OF33 > OF50 > OF100 > OF67 > CF，表明低比例有机肥（33%）最有利于促进大团聚体形成和提高土壤机械稳定性。

3.7 相关性分析与冗余分析

Mantel检验和Pearson相关矩阵显示，茶树生长指标与MBP、TN、TP、AN、AK、BG、ACP、NAG、LAP显著正相关。团聚体稳定性指标（R_0.25、MWD、GMD）与MBP、AN、AK、BG、NAG显著正相关。RDA分析中，以茶树生长指标为约束变量时前两轴解释84.29%总方差；以上壤理化性质和胞外酶活性为约束变量时前两轴解释51.15%总方差，TK、AP、MBP、BG、TN与OF50显著关联，ACP、OM、MBC、LAP与OF100显著关联，AN与OF33和CF正相关；以团聚体稳定性指标为约束变量时前两轴解释82.14%总方差。

讨论总结

研究人员指出，长期单施化肥驱动了亚热带红壤茶园酸化和养分失衡，形成了低pH、易得氮、低碳活性的极端资源基质。梯度有机肥施用通过腐殖质官能团提供酸碱缓冲，提升了土壤pH，有机质含量随替代比例线性增加。土壤养分呈现显著阈值分化：33%替代处理的全氮、碱解氮和速效钾含量最高，体现了有机肥缓释养分与无机速效养分的时序互补供应；而全磷和速效磷随替代比例持续下降，归因于有机磷形态释放缓慢、植物与微生物吸收增强以及pH介导的磷吸附平衡变化。

生态酶化学计量分析表明，OF100处理中碳获取酶BG活性和营养限制向量长度均达最大值，表明全有机肥处理下微生物面临严重碳限制，迫使群落偏向A型矿化策略。所有处理的向量角度均大于45°，指示相对磷限制贯穿整个施肥梯度，且向量长度与角度呈显著正相关，揭示碳磷限制协同加剧。基于微生物碳泵理论可推断，OF33和OF50处理下养分供给匹配微生物代谢需求，群落转向Y型生长策略，同化碳更多转化为微生物残体形成稳定有机碳。

BIOLOG-ECO分析证实，OF33和OF50处理下的AWCD、碳源丰富度、香农多样性和Levins生态位宽度均达到峰值，表明该替代范围能最大程度拓宽微生物碳利用光谱。随机森林鉴定出氨基酸、小分子羧酸和多聚物等为区分各处理的特征碳源。中等替代处理微生物对氨基酸和低分子有机酸的利用能力较强，全有机肥处理偏好多糖底物，纯化肥处理优先利用单糖和酚酸。

团聚体方面，33%替代处理在R_0.25、MWD和GMD上表现最优，表明适量有机胶结物质促进了稳定大团聚体的形成。相关性分析显示团聚体稳定性与OM、MBP、ACP高度显著正相关。研究人员提出了一条合理的调控路径：适量有机肥补充有机胶结剂优化团聚体结构，构建分层多样的微生物微生境，多样化的孔隙结构支持广谱碳利用微生物群落的定殖，碳磷获取酶的平衡分泌同步增强有机质固持和有机磷活化能力，协同缓解茶园多重养分限制。

研究结论

基于对土壤理化性质、胞外酶化学计量、BIOLOG微生物碳代谢特征、土壤团聚体组成及多维冗余分析的综合分析，本长期田间试验表明，有机氮替代通过土壤理化性质、团聚体结构和微生物碳代谢的协同变化，而非单一途径，调控亚热带茶园土壤的碳、氮、磷养分循环。在本长期田间试验条件下，基于所测量的指标（土壤pH、速效氮钾、微生物生物量磷、团聚体机械稳定性、胞外酶化学计量和微生物碳源利用），33–50%的有机替代比例可被推荐为适宜的应用范围。该施用方案能在多个方面产生协同效益：同时缓解茶园土壤酸化，拓宽微生物碳利用生态位，改变碳磷获取酶的分泌模式，改善土壤团聚体结构，同时促进茶树叶片叶绿素积累和更高的百芽重。相比之下，单独施用化肥与微生物碳源利用谱的收窄相关，而100%有机肥则诱导了严重的微生物碳限制（以最高向量长度为标志）。由于实验在单一地点、使用单一茶树品种和单一有机肥类型进行，这些发现对其他茶园条件的普适性有待进一步验证。对于与本实验条件相似的酸性红壤茶园，推荐使用有机肥提供总氮需求的33–50%；具体施用量应根据所用有机肥的含氮量进行计算。本研究为红壤茶园减施化肥和可持续绿色栽培提供了全面的土壤微生态学理论基础。