针对上述问题，中国科学院沈阳应用生态研究所方运霆研究员团队与中国科学院地理科学与资源研究所张扬建团队等建立了15N标记量化土壤N2O和N2产生速率的技术，对中国18个森林土壤反硝化作用N2O和N2产生速率的温度敏感性进行了研究，探讨了主要驱动因素，并利用生态系统过程模型（DyN-LPJ）结合室内培养的温度敏感性结果模拟了未来变暖情形下全球森林土壤反硝化作用气态氮释放的响应。研究的18个森林站点从南到北相距4000km，纬度相差33 （图1），覆盖了广泛的气候梯度（年均温相差25 C，年降水量相差2000mm）。
该图基于自然资源部标准地图服务网站下载的审图号为GS(2020)4619号的标准地图制作，底图边界无修改
研究发现，中国不同气候带森林土壤的反硝化作用 N2O 和 N2 产生速率随着温度的升高呈指数增长并呈现一定的地理格局（图 2 ）。然而，不同森林土壤之间 N2O 和 N2 释放的 Q10 值与土壤理化和微生物指标及气候因素无相关性，表明不同气候带下的森林土壤 N2O 和 N2 释放具有一致的温度敏感性， Q10 值分别为 2.1 0.5 和 2.6 0.6 （图 3 ）。该研究的陆地生态系统反硝化作用 Q10 值（ 1.4-3.2 ）与已报道的水生生态系统反硝化作用 Q10 值（ 1.6-3.5 ）非常相近（图 4 ）。该研究表明不管是土壤还是海洋沉积物环境，反硝化作用的温度敏感性相对来说比较一致，有利于模型模拟和预测全球变暖情况下反硝化作用的响应。 N2O 是一种强效的温室气体。气候变暖将促进更多的 N2O 排放，更多 N2O 释放将进一步加剧变暖，对气候变化形成正反馈。该研究还发现， N2 产生的温度敏感性高于 N2O ，表明气候变暖会使反硝化作用更加彻底，使土壤气态氮损失更多以 N2 的方式进行。
图2 18个森林土壤反硝化作用N2O、N2、N2O+N2产生速率对温度的响应
目前，大多数森林仍然处于氮限制状态，气候变暖引起的气态氮损失和由此产生的氮限制很有可能进一步限制陆地生态系统的初级生产力及其碳汇功能。为进一步揭示全球尺度上气候变暖对森林土壤气态氮损失的影响，该研究将森林土壤反硝化作用 N2O 和 N2 释放的 Q10 值与生态系统过程模型 DyN-LPJ 结合，预测了在全球变暖和其他全球变化因子（比如氮沉降、降水、大气 CO2 ）改变的背景下未来全球森林土壤反硝化作用气态氮的释放。模型的结果显示，在 1991-2000 年全球森林生态系统土壤反硝化作用 N2O 和 N2 的释放速率分别为 3.8 和 81.9 Tg N yr-1 ，气态氮释放以氮气为主。在 SSP2-4.5 和 SSP5-8.5 情景下， 2100 年森林土壤反硝化作用 N2O 的释放速率分别升高至 6.3 和 7.8 Tg N yr-1 ， N2 的释放速率将分别为 121.1 和 133.0 Tg N yr-1 （图 5 ）。
图5 DyN-LPJ模型模拟不同温室气体排放情景模式下2100年全球森林生态系统反硝化作用N2O和N2释放速率与2000年相比的年增量
a-c和d-f图分别为SSP2-4.5和SSP5-8.5情景模式下的气态氮释放模拟结果
综上所述，该研究在大陆尺度上量化了森林土壤反硝化作用产生的N2O和N2的温度敏感性，为模型模拟提供重要的数据参数，有助于加深对未来变暖背景下森林生态系统复杂的碳氮耦合过程及其对变暖反馈的理解。
这项研究成果以“Universal temperature sensitivity of denitrification nitrogen losses in forest soils”为题于2023年6月22日发表在Nature Climate Change上（IF2021=28.66）。中国科学院沈阳应用生态研究所于浩明博士生为第一作者，方运霆研究员和中国科学院地理科学与资源研究所张扬建研究员为共同通讯作者。合作者还有挪威生命科学大学Jan Mulder教授和Peter Dorsch教授、美国纽约州立大学朱伟兴教授、中国科学院青藏高原研究所旭日研究员等。该研究受到国家重点研发计划（2016YFA0600802)和中国科学院卢嘉锡国际团队项目（GJTD-2018-07）等项目资助。
