本文档属于 类型a，即报告了一项单一的原创性研究。
本研究由中国浙江大学环境与资源科学学院顾宝晶教授团队牵头，联合了浙江大学政策模拟实验室、英国生态与水文中心、爱丁堡大学、荷兰环境评估署等多家机构的国际研究团队共同完成，主要作者包括Cui Jinglan, Zhang Xiuming, Reis Stefan, Wang Chen, Wang Sitong, He Peiying, Chen Hongyi, van Grinsven Hans J. M. 和 Gu Baojing。该项重要研究成果于2023年6月22日在线发表于国际顶级学术期刊《Nature Sustainability》（第6卷，2023年10月刊，页码1166-1176）。
本研究聚焦于全球变化生态学与农业环境科学的交叉领域，核心是探究气候变化关键驱动因子——大气二氧化碳浓度升高（elevated CO2, eCO2）对全球农田氮循环的深刻影响及其连锁效应。
农田是人类粮食安全的基础，也是地球上氮素流动规模最大的生态系统。自工业革命以来，大气CO2浓度已上升了约47%，预计到本世纪末将达到600-1000 ppm。CO2作为主要的温室气体，同时也是一种“气体肥料”，通过促进植物光合作用和提高初级生产力（即CO2施肥效应）来刺激陆地生物圈的碳固存。然而，生物圈的碳固存能力在很大程度上受到氮素有效性的限制。因此，在CO2浓度升高的背景下，氮循环如何响应并反馈于碳循环，是决定未来农田生态系统是碳汇还是碳源的关键，也是评估气候变化下粮食生产和环境可持续性的核心科学问题。
尽管已有大量研究关注气候变化（如变暖、极端天气）对作物产量的影响，以及CO2施肥效应对碳循环（如净初级生产力、土壤有机碳）的促进作用，但关于CO2浓度升高对复杂农田氮循环的整体量化研究却严重缺乏。先前研究多集中于单一或少数氮通量变量（如N2O排放、氮素利用效率），缺乏从系统视角综合考虑氮循环各环节的联动响应。此外，这种响应可能存在显著的空间异质性，从而影响区域发展和国家间的不平等。因此，量化eCO2对全球农田氮循环的综合影响、预测其未来变化并评估其社会经济价值，对于完善地球系统模型、指导农业可持续发展和气候政策制定具有迫切的科学意义和现实需求。
本研究旨在填补这一知识空白，具体目标包括：1) 基于全球eCO2模拟实验数据，系统评估农田关键碳、氮循环变量对eCO2的响应；2) 利用耦合模型，在0.5°×0.5°的空间分辨率下，预测多种未来情景下（2050年）全球农田的氮收支变化；3) 开展影响评估，对eCO2引起的生态系统和人类健康影响进行货币化估值。
本研究包含三个主要且环环相扣的步骤：全球整合分析（Meta-analysis）、情景模拟与预测、货币化影响评估。整个工作流程逻辑严密，从实验观测到全球模拟，再到社会经济影响量化，构成了一个完整的研究链条。
第一流程：建立全球数据库与整合分析
此流程旨在从已有实验研究中提炼eCO2对农田碳氮循环影响的普适性规律。 1. 数据收集与处理：研究团队系统性地检索并筛选了全球范围内通过自由大气CO2富集实验（FACE）、开顶式气室（OTC） 和温室/生长箱（GC） 等方法模拟CO2浓度升高的田间和室内实验研究。最终，建立了一个包含全球众多试验站点的数据库，生成了总计1003个响应比，覆盖小麦、水稻、玉米、大豆等多种作物类型。 2. 变量提取与计算：从已发表文献中提取了eCO2处理组和对照组（环境CO2浓度）中关键变量的均值，变量涵盖碳循环（作物产量、土壤呼吸、土壤有机碳）、氮循环（生物固氮、氮素矿化、氮素利用效率、反硝化、氨挥发、氮氧化物排放、N2O排放、硝酸盐淋溶径流）以及植物组织氮含量和碳氮比等。使用响应比作为效应大小的度量，其计算公式为：RR = ln(X_eCO2 / X_aCO2)，其中X代表变量均值。 3. 统计分析与调节因子检验：采用自助法（Bootstrapping） 计算了各变量响应比的总体均值及其95%置信区间。通过元回归分析，检验了CO2升高幅度、年均温、年降水量以及作物类型等潜在调节因子对响应模式的影响。所有统计分析在MetaWin和R语言的metafor包中完成。
第二流程：全球农田氮收支情景模拟
此流程旨在将第一流程得到的定量响应关系“上推”到全球尺度，预测未来变化。 1. 模型构建与数据融合：研究采用耦合人类与自然系统模型（Coupled Human and Natural Systems model, CHANS）作为核心模拟工具。为了获得更可靠、空间分辨率更高（0.5°×0.5°）的全球农田氮收支数据，研究创新性地将全球营养模型（Global Nutrient Model, GNUM，属于IMAGE综合评估模型框架的一部分）的输出结果与CHANS模型相结合。具体而言，将GNUM生成的网格化氮收支数据输入CHANS模型，并利用CHANS模型中嵌入的国家级历史数据进行验证和优化，从而最小化预测的不确定性。 2. 情景设计：研究设计了两大类情景进行对比：基线情景（无气候变化，假设CO2浓度自2020年起保持固定）和eCO2情景（仅考虑CO2浓度升高这一单一气候变化因子，不考虑伴随的升温和降水变化）。每类情景下又包含三个子情景，分别对应不同的共享社会经济路径（SSP） 和典型浓度路径（RCP） 组合，以涵盖不同的未来发展可能性： * 可持续社会路径：SSP1（基线）与 SSP1-RCP1.9（eCO2） * 中间路径（一切照旧）：SSP2（基线）与 SSP2-RCP4.5（eCO2，作为主要分析情景） * 分层社会路径：SSP4（基线）与 SSP4-RCP6.0（eCO2） 3. 模型参数化与模拟：将整合分析得到的eCO2对作物产量、籽粒氮含量、氮素利用效率、生物固氮、各种氮损失通量等关键参数的响应比，作为调节因子整合到CHANS模型中。例如，未来作物的氮收获量计算中，既考虑了eCO2带来的增产效应，也计入了籽粒氮含量下降的效应。模型考虑了氮收支各组分（化肥、粪肥、生物固氮、大气沉降等输入；作物收获、氨挥发、N2O排放、硝酸盐流失等输出）之间的质量平衡关系，并允许人为氮投入（如化肥）作为一个灵活组分，随未来土壤养分状况和农业生产适应性调整而变动。利用蒙特卡洛模拟进行了1000次迭代，以评估预测结果的不确定性。
第三流程：货币化影响评估
此流程旨在量化eCO2变化带来的社会经济收益或损失。 1. 影响类别界定：评估涵盖五个方面的影响：生态系统效益（因活性氮损失减少而避免的生态系统服务损害）、人类健康效益（因空气污染物前体物减少而避免的健康损害）、作物增产收益、化肥节省收益以及气候影响（N2O的增温效应与NH3/NOx形成气溶胶的冷却效应的净值）。 2. 货币化方法：采用效益转移法和损害成本法进行估值。例如，生态系统效益以欧盟的氮损害成本评估为基础，根据各国的支付意愿和购买力平价进行调整；人类健康效益通过量化氮排放对PM2.5浓度的贡献及其健康损害来计算；作物收益和化肥节省则分别基于作物价格和化肥价格计算。评估在国家尺度进行，再汇总到区域和全球。
整合分析结果：eCO2会引发全球农田碳氮循环的协同强化。碳循环方面，作物产量平均提高21%，土壤呼吸增加25%，土壤有机碳小幅增加6%。氮循环方面，响应更为显著和复杂：生物固氮大幅增加55%，氮素矿化增加22%；氮素利用效率提高19%，这直接导致活性氮损失全面减少，其中硝酸盐淋溶径流减少45%，氨挥发减少21%，氮氧化物排放减少33%。然而，N2O排放增加了29%，主要归因于反硝化过程增强（+24%）。植物组织呈现“氮稀释”现象，籽粒、叶片和茎秆的氮含量分别下降7%、15%和10%，叶片碳氮比提高19%。这些结果共同表明，eCO2下碳可用性的增加刺激了微生物活动和氮循环，而增强的氮循环又反过来缓解了碳同化的氮限制，形成了一个正反馈回路，维持了CO2施肥效应对作物产量的促进作用。
情景模拟结果：以SSP2-RCP4.5（“一切照旧”）情景下的2050年预测为例，与无气候变化的基线情景相比，eCO2将显著改变全球农田氮收支格局： * 氮输入总量减少27 Tg N/年：主要得益于化肥施用量大幅减少34 Tg N/年，粪肥投入也减少5 Tg N/年。虽然生物固氮增加15 Tg N/年，但不足以抵消化肥的减量。大气氮沉降因NH3和NOx排放减少而降低3 Tg N/年。 * 氮收获量增加12 Tg N/年：尽管籽粒氮含量下降，但CO2施肥带来的增产效应更强，导致作物带走的氮总量增加。增产热点区域主要分布在东亚、南亚、北美五大湖地区等粮食主产区。 * 氮盈余（损失）减少39 Tg N/年：具体表现为活性氮损失总量减少46 Tg N/年，其中硝酸盐淋溶径流减少32 Tg N/年，氨挥发减少12 Tg N/年，N2O减少2 Tg N/年，NOx减少0.9 Tg N/年。但同时，非活性的N2排放因反硝化增强而增加了7 Tg N/年。 * 空间异质性显著：氮输入和损失减少的“热点”区域主要在中国、印度、北美和欧洲等农业集约化程度高的地区。氮素利用效率的增幅在空间上不均，高基础效率地区的相对提升更大，这可能加剧区域不平等。
货币化影响评估结果：仅考虑eCO2这一单一因子（未计入升温等负面效应），其在205年将为全球带来6680亿美元的社会效益。其中，生态系统效益最大（3590亿美元），其次为人类健康效益（1280亿美元） 和作物增产收益（1240亿美元），化肥节省贡献了约570亿美元的收益。气候影响的净成本很小（约-10亿美元）。中国、印度、北美和欧洲将获得总效益的65%，其中中国和印度在获得巨大健康与生态效益的同时，也需承担因N2O排放增加带来的少量气候成本。撒哈拉以南非洲将主要从作物增产中获益（206亿美元）。
本研究得出核心结论：未来大气CO2浓度升高，作为一个独立的全球变化驱动因子，将通过协同强化碳氮循环，显著提高全球农田的氮素利用效率和生物固氮能力。这不仅能在2050年前增加粮食产量、大幅减少化肥需求和活性氮污染，还可产生巨大的社会经济净效益。这为在气候变化背景下，协同实现粮食安全、环境可持续和人类健康目标提供了新的机遇和科学依据。
研究的科学价值在于：首次从系统性、全球性和空间显式的角度，全面量化了eCO2对农田氮循环各关键通量的影响，并预测了其未来变化格局。研究将实验生态学的整合分析、复杂系统模型模拟和货币化评估有机结合，构建了一个完整的研究范式。其应用价值在于：研究结果明确指出，未来农业管理实践必须考虑并适应CO2升高带来的氮循环变化，例如利用生物固氮增加和氮需求下降的契机，优化化肥-有机肥配施，选育适应高CO2的作物品种。同时，研究警示籽粒氮含量下降可能影响膳食蛋白质供应，需调整未来营养策略。
文章在讨论部分特别指出，本研究的范围限定于eCO2的直接效应。在现实世界中，CO2升高总是伴随着全球变暖、降水格局改变和极端天气事件。这些因素可能产生负面冲击，其损害成本未在本评估中计入。因此，eCO2的净效益能否完全抵消气候变化的总体负面影响，仍需更深入的综合分析。未来的研究需要采用更复杂的集成评估方法，全面解析氮循环对多重气候压力因子的响应机制，这是设计有效气候政策和可持续农业系统的关键科学基础。