《Agronomy for Sustainable Development》：Sustainable crop production under future climates: national mapping of cover-crop rotations and soil organic carbon–N2O trade-offs
编辑推荐：
可持续作物生产需要既能维持产量又能限制氮损失和温室气体排放的管理措施。然而，国家尺度的评估很少在未来气候条件下比较覆盖作物轮作与现行种植制度，同时联合评价生产力、氮淋溶、土壤有机碳（SOC）和氧化亚氮（N2O）。这限制了对区域权衡的识别和因地制宜轮作制度的设计
可持续作物生产需要既能维持产量又能限制氮损失和温室气体排放的管理措施。然而，国家尺度的评估很少在未来气候条件下比较覆盖作物轮作与现行种植制度，同时联合评价生产力、氮淋溶、土壤有机碳（SOC）和氧化亚氮（N2O）。这限制了对区域权衡的识别和因地制宜轮作制度的设计。为填补这一空白，研究人员使用DSSAT模型，在六种气候模式和RCP4.5与RCP8.5情景下，比较了含豆科和非豆科覆盖作物的多样化4年轮作与常规（business-as-usual）轮作在德国各农业生态区的表现。评估指标包括谷物单位产量、氮利用效率、氮淋溶、SOC变化和净温室气体平衡。直接和间接N2O排放采用依赖氮投入和淋溶的IPCC Tier 2方法估算；净土壤CO2交换由模拟的SOC变化推算。覆盖作物轮作在所有情景和区域均提高了谷物单位产量和氮利用效率，并减少了氮淋溶。豆科覆盖作物的生产力增益最强，谷物产量增加约16–33%（取决于区域和情景），同时氮淋溶下降约40%。SOC增加，但增益随时间放缓。N2O响应取决于区域、气候和氮管理。到本世纪末，SOC变化与N2O排放合计使净田块尺度温室气体平衡较常规轮作降低约35–45%。该研究的新颖贡献在于开展了全德国范围、相对于基线（baseline-relative）的评估，联合绘制了未来气候下不同覆盖作物轮作的生产力、氮损失和SOC–N2O权衡。这些基于模型的田块尺度通量比较，在相对差异和空间格局方面比绝对量级更为稳健。
研究背景方面，农田生产粮食、纤维和燃料，同时也是温室气体（GHG）排放的重要来源。2021年德国农业活动排放56.3 Mt CO2-eq，占总GHG排放的7.4%，其中CO2约占80%，甲烷（CH4）占11%，N2O占6%。在作物生产中，氮（N）肥管理是排放的主要驱动因素，约50%的施氮未被植物吸收，导致N2O排放升高；德国谷类作物系统中氮肥使用约占GHG排放的80%，因此提高氮利用效率至关重要。除减排外，土壤有机碳（SOC）固存被广泛提出作为缓解途径（如“4 per 1000”倡议），但其潜力因饱和效应、可逆性以及与氮循环的相互作用（可能通过增加N2O排放抵消增益）而存在争议。多样化轮作和整合覆盖作物（CCs）可能有助于改善SOC、提高氮利用效率并减少环境损失，但这些效应并非普遍正面且存在内在权衡：增强SOC的措施在潮湿或排水不良条件下可能增加N2O排放，而CCs在水分受限环境中可能因土壤水分竞争而降低产量。实证证据显示既有增产（如豆科CCs下可达14%）也有减产（干旱条件下），结果强烈依赖于气候、土壤性质和管理实践。尽管研究广泛，多数研究孤立或在站点尺度评价这些过程，限制了大尺度决策的适用性；很少有研究在国家尺度上同时评估当前和未来气候条件下的产量、氮动态（包括淋溶）、SOC固存和净GHG平衡。此外，CC系统相对于常规（BAU）轮作的优势往往未被明确量化，难以评估其实际效益；相对于基线系统表述结果可提高可解释性并突出管理优势。同时，结果的区域差异常被报道，但其驱动因素（如气候变异性、土壤特征、水分可用性）往往未被明确讨论，限制了结果的可迁移性。因此，本研究相对于BAU系统评估改进轮作在德国各农业生态区的表现，并将空间差异与气候、土壤和水分相关约束联系起来。研究人员使用DSSAT模型评估在當前和未来气候情景（RCP4.5和RCP8.5）下含不同CC组合的多样化4年轮作。研究目标为：（i）确定将豆科和非豆科CCs纳入4年轮作的不同情景对作物产量、SOC发展、氮相关变量和GHG排放相对于无CC和BAU情景的影响；（ii）探索这些建议管理选项在未来气候条件下大空间尺度的潜在结果。
关键技术方法方面，研究使用基于过程的农业生态系统模型DSSAT 4.8.0，以日步长模拟作物生长、土壤碳氮动态和水分平衡；SOC采用基于CENTURY的模块模拟。模型基于德国西南部6个站点2008–2018年十个生长季（收获年2009–2018）的高质量数据集（60站点-年）进行校准，并用德国联邦植物品种局的官方品种试验数据在四个农业生态区进行独立验证。研究设计了六种改进CC轮作处理（基于两种4年结构，CC类型分豆科/非豆科、频率为0、2、3次）和两种BAU基线处理。气候数据采用CMIP5/CORDEX的六种GCM（RCP4.5三种、RCP8.5三种），降尺度至约12.5 km分辨率；土壤数据来自WISE数据库（聚合至10 km分辨率）。模拟在德国15 km分辨率点图层上于2006–2099年执行八种处理，共1,340,064次模拟，并通过R脚本自动化；CO2施肥效应按各RCP情景动态考虑。采用德国联邦农业部门的“谷物单位”（GU）概念比较轮作，NUE以每kg施用合成氮的GU计；N2O排放按IPCC（2006）Tier 2指南结合德国区域特定排放因子估算直接和间接排放，GHG通量为N2O（GWP 100年为273）与净CO2排放之和。结果采用kriging插值至1 km分辨率并用农田矢量图层掩膜，统计上用单因素方差分析和Fisher's LSD（α=0.05）比较。
研究结果方面，（1）谷物单位及其氮利用效率的空间格局：与BAU相比，改进轮作情景在所有区域和RCP下均一致提高GU，无论是否包含CCs；GU增幅从北部16%到南部33%。含豆科CCs的情景（如N1–L2和N0–L3）在各地理位置、时期或气候情景下均表现有利，平均比N0–L0高23%；而无豆科CCs情景在2070s RCP8.5下GU最低，比无CC情景在北部、东部、南部和西部区分别降低约5.1%、4.0%、6.3%和8.1%。空间上，北部区BAU下约83%耕地GU为4.7–5.5 Mg GU ha?1，改进轮作无CC和含CC（N0–L3）分别使其提高13.7%和34.7%，约75%区域在N0–L3下达6.4–7.3 Mg GU ha?1；东部和西部区BAU主要为4.7–5.5，N0–L0和N1–L1为5.5–6.4，N0–L3扩展至6.4–8.2；南部区各情景GU最高，N0–L3下超过70%面积达7.3–8.2，而BAU峰值仅4.7–5.5。NUE随豆科CCs纳入而逐渐提高；改进轮作无CC和含非豆科CCs的NUE均显著高于BAU，2070s略高于2030s。BAU下全德耕地NUE为25–42 kg GU kg?1施氮，而豆科CC情景下北部约70%区域平均提高100%、东部约一半区域提高130%、南部90%区域提高130%、西部约一半区域提高100%；无CC情景NUE介于两者之间，四区大部分为42–51 kg GU kg?1。N0–L0相对BAU的较高NUE反映了改进轮作结构本身的效应。含豆科CCs情景比BAU或无CC情景GU产量平均提高30–60%（取决于区域和时期），但干旱年份CCs的水分竞争可能削弱或抵消水分受限区的产量增益，因此响应应视为区域和气候条件依赖的，需适应性管理（如提早终止CC或采用低耗水CC混合）。（2）土壤有机氮、氮淋溶和净SOC平衡的空间格局：土壤有机氮在各情景间差异轻微但显著，CC纳入优于BAU和无CC，N1–L1和N0–L3更明显，2070s更突出；2030s平均约5.7 Mg ha?1，2070s约6.8 Mg ha?1，北部和西部较低、南部较高。CCs使北部和东部土壤有机氮增加（多数区域7–8 Mg ha?1）。氮淋溶方面，CCs明显减少损失，北部和东部效应最强：BAU和无CC下平均约40 kg N ha?1，含CC轮作降至约20–25 kg N ha?1；南部和西部趋势相似但差异较小，因BAU淋溶已较低。时间变化总体较小，但RCP8.5下南部和西部2070s淋溶略高。淋溶减少是轮作设计、氮管理和CC纳入的综合效应，不能仅归因于CCs。（3）土壤通量的时间格局：RCP4.5下北部BAU和无CC情景土壤碳固存可忽略（中位数接近零），而累积N2O排放持续上升至2046年33.8和2096年76.5 Mg CO2-eq ha?1；RCP8.5下BAU和N0–L0到2096年土壤CO2释放分别相当于12.4和14.7 Mg CO2-eq ha?1，含CC情景平均固碳相当于23.2 Mg CO2-eq ha?1，使GHG排放降至49.5对89.0 Mg CO2-eq ha?1（RCP4.5），RCP8.5下程度较小。CC情景累积SOC轨迹的曲率显示本世纪后期十年增量递减（趋于平台），RCP8.5下增温加速分解使净累积相对世纪中叶减少。东部和西部BAU与无CC的累积SOC增益相近（到2096年约17和24 Mg CO2-eq ha?1），无CC的N2O排放较BAU低7 Mg CO2-eq ha?1，净GHG略低。含CC情景SOC增益更高且N2O排放更低。净GHG在两种RCP下到世纪末接近中性或略正；南部区效益最明显，N1–L1和N0–L3较BAU的N2O排放分别减少10.5和23.1 Mg CO2-eq ha?1。全德范围，RCP4.5下BAU和无CC的土壤碳固存到世纪中叶平台化于约23 Mg CO2-eq ha?1，而CC情景继续增至世纪末51；RCP8.5下所有情景SOC固存到世纪中叶达最大后下降。N2O排放在RCP间变化不大（BAU约87、CC情景约65 Mg CO2-eq ha?1）。净GHG差异主要反映SOC固存差异：BAU在RCP4.5和RCP8.5下平均约69和75，无CC为58和72，含CC情景为18和28 Mg CO2-eq ha?1。含CC改进轮作使全德净GHG较BAU降低约40%，豆科CC存在时更明显；含两或三个CCs的改进轮作N2O排放较BAU和无CC分别低18%和13%。
讨论部分总结：本研究量化了直接田块尺度GHG通量（N2O直接+间接和SOC变化净CO2），未包括上游/农场外过程（如肥料和种子生产运输、田间作业燃料），这些通常由生命周期评价处理；纳入这些组分可能调节CCs的表观缓解效果，尽管豆科CCs可降低矿质氮投入从而可能增强净气候效益。SOC累积趋于平台且响应随区域和气候而异；干旱年份水分竞争可降低产量和SOC增益，因此在水分受限或排水不良条件下建议适应性管理。未来研究方向包括：（i）将过程模型与LCA耦合以捕捉农场能源使用和上游投入；（ii）多模型集合以界定SOC–N2O–产量相互作用的结构不确定性；（iii）在数据允许处靶向验证水分平衡和N2O通量；（iv）测试将CCs与其他保护性措施结合的因地制宜方案。全德模拟需考虑土壤类型、气候和耕作实践的区域差异；本研究中播种和终止日期标准化以保证情景可比，但可能未充分捕捉区域差异。CC轮作的大规模采用取决于社会经济和管理约束（额外劳动力、成本、知识要求和农民感知风险），需考虑政策激励、咨询支持和区域特定管理约束。研究结果对德国以外其他温带谷物产区也有相关性，但受益量级应视为区域特定并依赖土壤、气候和管理相互作用。建模假设方面，因DSSAT中并非所有模拟轮作作物都有作物模块，部分作物用功能相似的替代作物表示；模型评价nRMSE一般约20–30%，因此情景间小差异不应过度解读，结果最好按稳健的方向性响应、相对差异和空间格局解释，而非绝对预测。SOC（0–30 cm）和土壤CO2通量性能此前已在德国西南部10年6站点数据集上评价，而全德范围适合直接验证N2O通量和氮淋溶的长期观测仍有限；水分平衡组分未在全德尺度验证，干旱年份结果应谨慎解释。
研究结论部分翻译：为评估含不同覆盖作物组成和氮管理的改进轮作在未来气候下如何影响全德生产力、氮损失、SOC和净温室气体平衡，本研究比较了改进4年轮作与常规系统。其新颖贡献在于开展了全德国范围、相对于基线的评估，联合绘制了未来气候下不同覆盖作物轮作的生产力、氮损失和SOC–N2O权衡。结果突出了将豆科CCs纳入改进轮作并配以考虑豆科氮 credit 给后茬作物的建议氮管理的显著效益，可增强生产力和长期可持续性。这些集成系统提高了生产力和NUE，减少了对合成和有机氮投入的依赖，改善了土壤有机氮动态同时减少氮淋溶。轮作设计、CC纳入和改进氮管理的综合效应不仅增加了SOC固存还减少了N2O排放，使净田块尺度GHG平衡较BAU实践降低约40%，豆科与非豆科CCs交替时降低尤为明显。总体而言，结果强调了集成轮作和氮管理策略在德国各农业生态区促进气候变化缓解和可持续农业生产的潜力。然而，这些效益在长时间尺度上实现并随区域土壤、气候和管理条件而变化。此外，结果应结合模型和数据局限性解读，特别是氮淋溶估算和标准化管理假设。未来工作应进一步细化区域特定管理实践并评估采用约束以支持这些系统的实际实施。
