酸性茶园土壤羊粪氮替代30 d培养：累积N2O、CO2排放与氮循环功能基因响应

《Agriculture》：Sheep Manure Nitrogen Substitution in Acidic Tea Plantation Soil: Gas Emissions and Nitrogen-Cycling Gene Responses During a 30-Day Incubation

　　

回收粪肥氮（N）可减少矿物肥料投入，但其对土壤排放的影响需要评估。一项为期30 d的无植物培养试验比较了单施尿素（N1）、15%（N2）和30%（N3）腐熟羊粪氮替代以及不施肥对照（N0）。酸性茶园土壤（pH 4.24；有机质49.4 g kg^?1）与粪肥（pH 7.6）在80%充水孔隙度（WFPS）下培养。第1–30天累积分子N₂O在N0–N3中分别为894.1 ± 15.7、911.5 ± 18.4、960.4 ± 6.0和900.6 ± 13.4 ng g^?1干土（平均值 ± 标准误；n = 3）。与N1的差值：N2为48.9 ng g^?1（未校正95%置信区间，2.8–95.1），N3为?10.9 ng g^?1（?57.0–35.3）；Holm校正后均不显著。N3累积CO₂（1478.0 μg g^?1）高于N1（950.6 μg g^?1）。两种替代处理均在N₂O通量下降后的第30天提高了DNA归一化的氨氧化细菌（AOB）amoA和nirS。氨氧化古菌（AOA）amoA和nirK无显著处理×时间交互，nosZ无显著预设日期特异性对比。探索性结构方程模型（SEM）和Mantel分析未确立微生物驱动因素：分类日期校正后无N₂O预测因子保持显著，且无Mantel关联通过错误发现率校正。30%替代值得田间评估，但实验室不显著不能确立排放等效性或农艺效益。

论文主体解读：农业土壤是N₂O的重要来源，集约施肥茶园在强酸性条件下尤其受关注。粪肥N部分替代矿物N可减少矿物肥料投入并回流养分，但替代会改变速效N与有机态N平衡并引入有机质，排放响应取决于材料、速率和受体土壤。成熟羊粪呈碱性，与强酸性茶园土组合需要评估；已有报道在作物系统、材料、速率和土壤间不一致，且近期茶土培养多采用猪粪或稻草，尚需明确强酸性茶土中较低替代率下羊粪N对N有效性和排放的影响，并将矿物N、微生物响应时序与同期累积排放比较。研究人员在四川西部强酸性茶园土中开展30 d无植物密闭培养，设置不施肥对照（N0）、单施尿素（N1）、15%（N2）和30%（N3）腐熟羊粪N替代，田间等效总N 300 kg N ha^?1，在80%充水孔隙度（WFPS，表示土壤孔隙充满水的比例）下培养，比较累积N₂O、矿质N轨迹、N循环基因丰度及CO₂和CH₄响应。研究表明，在Holm校正标准下两种替代均未显著改变累积N₂O，N3增加累积CO₂，两替代处理在N₂O通量下降后的第30天提高DNA归一化AOB amoA和nirS；探索性关联未确立微生物驱动。其意义在于为强酸性茶土羊粪N循环提供短期室内生化筛选证据，并提示30%替代需田间评估，但实验室不显著不能证明排放等效、最佳速率或农艺效益。
关键技术方法（≤250字）：供试土壤采自中国四川西部茶园0–20 cm层（30°11′ N，103°22′ E），商业腐熟羊粪购自安徽富民生物有机肥有限公司（Bozhou, China）。采用完全随机设计，四处理三重复；气体监测采用独立密闭顶空瓶，土壤采样采用另设瓶，密闭顶空思路参考Hu et al.并调整。于25 °C、80% WFPS培养30 d，按日期采集气体，用气相色谱测定N₂O、CO₂和CH₄通量并以梯形积分计算累积排放；同步测定溶解性有机碳（DOC）、NH₄⁺-N和NO₃^?-N；提取DNA并用实时定量聚合酶链反应（qPCR）测定AOA amoA、AOB amoA、nirS、nirK和nosZ丰度；统计采用重复测量方差分析、Holm校正、Pearson/Mantel检验及探索性结构方程模型（SEM）。
3.1 温室气体通量的时间变化：通过30 d周期气体监测，N₂O通量随时间下降且处理×时间交互显著（p<0.001），Holm校正后N2与N3及N1的差异随日期改变；CO₂通量第1至3–5天升高后总体下降，N3第20天再升高（交互p=0.00242），第5天N2和N3高于N1、第20天N3高于二者；CH₄在净吸收与释放间波动，交互不显著（GG p=0.0691）。
3.2 累积气体排放与表观排放因子：通过30 d梯形积分，N0–N3累积N₂O分别为894.1±15.7、911.5±18.4、960.4±6.0和900.6±13.4 ng g^?1干土，六组两两比较Holm校正后均不显著，N2与N1差48.9 ng g^?1（+5.4%；未校正95% CI 2.8–95.1），N3差?10.9 ng g^?1（?1.2%；?57.0–35.3）；N3累积CO₂ 1478.0±53.8 μg g^?1高于N1 950.6±64.8和N2 1137.4±64.8，Holm校正后显著；累积CH₄仅N0与N3校正后差异显著（p=0.0402），施肥处理间不显著；表观N₂O排放因子N1–N3为0.0085%、0.0324%和0.0032%。
3.3 DOC与矿质氮动态：通过定期土壤采样，DOC第1至5天下降后较稳定（交互p<0.001），第7天N3低于N2；NH₄⁺-N两替代处理第1和30天低于N1、第5–20天高于N1（交互p<0.001）；NO₃^?-N第3天下降后向第30天增加（交互p<0.001），第30天N3>N1>N2（N3均值81.6 mg N kg^?1）；总提取矿质氮第3天N1和N2为435.7和525.9 mg N kg^?1而N0为28.7，第5天降至219.3和287.1，N0第20天达168.2 mg N kg^?1。
3.4 功能基因响应：通过第1、5、30天qPCR，DNA归一化AOB amoA和nirS在第30天两替代处理均高于N1，且发生在N₂O通量降至处理均值<700 ng kg^?1 h^?1后；AOA amoA和nirK无显著处理×时间交互（GG p=0.0501和0.345）；nosZ交互显著（GG p=0.0182）但N1–N3无显著预设日期特异性对比。
3.5 探索性关联与时间混杂：通过12个处理×天均值，N₂O通量与DOC正相关（r=0.905）、与NO₃^?-N负相关（r=?0.739），DOC与CO₂通量负相关（r=?0.741），三项Pearson关联BH-FDR后p<0.05；30项Mantel关联均未通过BH-FDR；探索性SEM分类日期校正后无N₂O预测因子保持显著，nirS–N₂O关联仅在49.7%的可估计重复标签映射中p<0.05，显示对气瓶与土壤瓶任意对齐敏感。
讨论部分总结：讨论认为两种替代均无累积N₂O显著变化，也无速率依赖增加证据；N0与N1不显著，湿培养背景排放可能涉及原生土壤N，但对照扣除不能分配肥料源与原生N。通量对比方向随采样日变化，单日比较不完整；与猪粪或稻草茶土培养减排报道不同，材料、土壤和培养条件差异未独立检验。矿质N方面，第1天24 h已有大量NH₄⁺-N积累，高NH₄⁺不能仅归因肥料N，原土总N 2.87 g kg^?1且N0矿质N池变化大，第3天大池不证明正激发；尿素快速水解、高湿限氧、强酸改变氨氧化器生态位及粪肥有机N逐步释放可作解释但未量化；第3–5天总提取矿质氮下降不能仅由NH₄⁺转NO₃^?解释，轨迹为净提取池而非总过程速率；第30天N3较高NO₃^?提示潜在淋失风险而非实测田间损失。DOC早期下降与N₂O下降同步，但日校正模型中关联不显著，不能区分底物效应与共同时间；80% WFPS不复制田间干湿循环，CH₄交换不能决定氧化还原状态。第30天AOB amoA和nirS增加发生在N₂O下降后，显示DNA基因池与气体交换时间解耦；DNA可积累或持久，瞬时活性依赖底物、氧和转录调控，不能证明底物耗竭或这些基因驱动早期排放峰；碱性粪肥可改变氨有效性和酸胁迫，但粪肥pH本身不量化中和能力或处理水平土壤pH变化；nosZ DNA丰度不测N₂O还原。管理上，15%和30%替代分别替代45和90 kg N ha^?1矿物N；30%值得田间评估，在测试条件下无统计检测累积N₂O增加；15%同样不显著，数据不能确定最佳或更安全速率。N3累积CO₂比N1多55.5%，粪肥C分解和原生C周转改变均可贡献但未分配，CO₂响应不证明粪肥C矿化、长期土壤C保留或净气候效益。实验受一个复合土、一种粪肥、三重复、30 d无植物固定温湿度培养限制，无根、排水和田间波动限制外推。
研究结论翻译：在此酸性土壤培养中，按Holm校正标准，15%和30%羊粪N替代均未显著改变相对于单施尿素的累积N₂O。30%处理增加累积CO₂，而两种替代处理在N₂O通量下降后均提高第30天DNA归一化AOB amoA和nirS。这些基因响应和探索性关联未确立因果排放机制。30%替代值得田间评估，但实验室不显著不能证明排放等效、最佳速率或农艺效益。田间评估应包括茶叶N吸收、产量、季节性排放以及高残留硝酸盐相关淋失风险。