气候变化与氮氧化物减排：地表臭氧响应未来的不确定性评估
一项由麻省理工学院（Massachusetts Institute of Technology）Emmie J. Le Roy、Anthony Y.H. Wong、Sebastian D. Eastham、Arlene M. Fiore和Noelle E. Selin共同完成的研究，于2025年5月22日发表在Environmental Science & Technology期刊上（Volume 59, Issue 15）。该研究聚焦于气候变化背景下，减少人为氮氧化物（NOx）排放对地表臭氧（O3）浓度的影响如何演变，旨在为未来不确定气候条件下的空气质量政策评估提供关键科学依据。
学术背景与研究目标 对流层臭氧是一种重要的温室气体，对人体健康和农作物生产力具有显著危害。尽管全球已实施多种臭氧前体物控制措施，但由于臭氧水平上升、人口变化等因素，臭氧相关的死亡率仍在持续增加。臭氧的形成与其前体物（主要为NOx和挥发性有机化合物，VOCs）之间的关系具有高度非线性。更为复杂的是，气候变化通过改变生物源排放（如土壤NOx、异戊二烯）、大气化学过程和输送条件，可能进一步影响这种非线性关系。因此，理解气候变化如何改变臭氧对前体物减排的响应，对于评估当前减排策略在未来气候下的持续有效性至关重要。
以往的研究多关注气候变化对臭氧绝对浓度的影响，或在使用未来排放和未来气象场的情景下，探究臭氧对排放变化的响应。然而，这些研究通常受限于单一的未来气候实现，未能充分考量气候系统内部变率（Internal Climate Variability）带来的不确定性。鉴于空气质量标准（如美国环保署基于三年内第四高日最大8小时平均臭氧浓度）关注的是臭氧分布的极端高值部分，忽略内部变率可能导致对政策效果的误判。
本研究旨在填补这一空白。其核心目标是：量化在历史与两种不同变暖程度（高变暖与低变暖）的未来气候情景下，地表臭氧对全球人为NOx排放减少10%的响应（文中定义为δO3-m10%enox）有何不同，并特别关注这种响应的概率分布变化，而不仅仅是平均值。研究选取了北半球中纬度三个具有高人为NOx排放和成熟监测网络的区域：北美东部（ENA）、西欧中部（WCE）和东北亚（NEA）。通过使用全球大气化学模型集合（ensemble），研究试图回答：未来气候变化是会增强还是会削弱NOx减排对改善臭氧污染的效果？这种影响在不同区域有何差异？其背后的驱动机制是什么？
详细研究流程与方法 本研究采用了集合模拟与气候-化学传输模型相结合的复杂工作流程，主要包括以下几个步骤：
气候情景与气象场生成：研究采用了社区大气模型（Community Atmosphere Model, CAM）v3.1模拟的历史（2000-2015年）及两种未来（2080-2095年）气候情景下的气象场。两种未来情景分别代表“高变暖”（无约束温室气体排放，2100年辐射强迫达10 W m⁻²）和“低变暖”（严格减排，2100年辐射强迫达3.7 W m⁻²）。为了捕捉气候内部变率，每个情景均包含5个不同的模型实现，这些实现通过对海表风应力施加微小扰动生成。气象场水平分辨率约为2°×2.5°，垂直26层。
大气化学传输模拟：研究者使用高性能版GEOS-Chem化学传输模型，并将其驱动气象替换为上述CAM输出，构建了“GCHP-CAM”模拟框架。模型运行在c48立方球网格（约200 km分辨率）和26个垂直层次上。为了评估水平分辨率的影响，还进行了一套c90分辨率（约100 km）的敏感性模拟。
排放与实验设计：为了隔离气候本身的影响，研究将基准人为排放（基于2014年社区排放数据系统，CEDS）在所有模拟时段内固定不变。而生物源排放（异戊二烯和土壤NOx）则根据模型中的气象条件（温度、光照、土壤湿度等）动态计算。特别考虑了二氧化碳（CO2）浓度升高对异戊二烯排放的抑制作用。闪电和野火排放则基于2014年气候态固定不变。
扰动实验与核心指标：对于每个气候情景的每个模型实现，研究者运行了两组模拟：“基准”模拟和“低NOx”模拟。“低NOx”模拟将全球人为NOx排放统一减少10%。核心响应指标δO3-m10%enox定义为两者日最大8小时平均臭氧浓度的差值。若δO3-m10%enox为负值，表示“NOx敏感”条件（减排降低臭氧）；若为正值或零，则表示“NOx不敏感”条件（减排导致臭氧不变或升高）。此外，研究还统计了每年4月至8月期间出现“臭氧负效益”（即δO3-m10%enox > 0）的天数比例，以及“高臭氧日”（定义为超过历史模拟期间区域特定95百分位阈值的天数）的频率变化。
模型评估与不确定性分析：研究通过与2014年6月和8月北美、欧洲和中国的实地观测数据对比，评估了GCHP-CAM模拟地表臭氧分布的能力。由于使用气候模式气象会导致系统性高估夜间臭氧（垂直混合过程模拟差异），研究者明确表示主要进行相对值分析（即比较不同情景下的δO3-m10%enox），而非关注臭氧绝对浓度。研究还对比了GCHP-CAM与使用再分析气象驱动的GCHP-M2模型在计算δO3-m10%enox上的差异，并评估了模型水平和垂直分辨率对结果的影响，确保了主要结论对模型设置的稳健性。
驱动机制诊断：为了解释未来δO3-m10%enox变化的原因，研究者分析了在基准排放下，未来气候情景中几个关键变量的变化：土壤NOx排放、生物源异戊二烯排放、过氧乙酰硝酸酯（PAN）浓度以及HOx自由基浓度。这些变量分别代表了NOx和VOCs源汇、自由基化学以及NOx输送和储存过程的变化，是影响臭氧生成化学机制的关键因素。
主要研究结果 研究通过集合模拟，获得了关于δO3-m10%enox在未来气候下演变规律的详细且具有统计意义的结果。
首先，在历史情景（2000-2015年） 下，模型模拟的区域平均δO3-m10%enox在北美东部和西欧中部均为负值（分别为-1.1 ppbv和-0.7 ppbv），表明这两个区域整体上处于NOx敏感状态，这与2000年代以来卫星和地面观测到的光化学状态一致。在东北亚，区域平均δO3-m10%enox也为负值（-0.4 ppbv），但其空间分布不均，部分网格呈现正值，反映了该区域在2014年左右高排放下部分城市地区存在NOx饱和光化学现象。
其次，未来气候下的响应变化呈现显著的区域差异。统计分析（Kolmogorov-Smirnov检验）表明，与历史情景相比： * 在北美东部和西欧中部，高变暖情景下δO3-m10%enox的整个概率分布向负值更小的方向偏移（即绝对值变小）。这意味着，虽然NOx减排仍能降低臭氧（响应符号未变），但臭氧浓度下降的幅度在高变暖的未来将小于历史时期。 * 在东北亚，情况截然相反。在高变暖和低变暖情景下，δO3-m10%enox的分布均向更负的方向偏移。这表明，未来在这两个区域，同样的NOx减排将带来比历史时期更大的臭氧浓度降幅。
第三，臭氧负效益天数频率的变化进一步印证了区域差异。在东北亚，高变暖情景下，区域平均的出现臭氧负效益（减排反而增臭氧）的天数比例，从历史时期的27.1%显著下降到23.7%。这意味着未来在该地区进行NOx减排时，出现适得其反效果的风险降低。而在北美东部和西欧中部，区域尺度的臭氧负效益天数比例没有发生统计学上的显著变化。
第四，结合基线臭氧水平变化的分析显示，在高变暖情景下，西欧中部和东北亚的“高臭氧日”频率相对于历史时期显著增加（分别超过三倍和两倍）。尽管如此，在所有区域和情景下，10%的NOx减排都能有效降低高臭氧日的频率。然而，减排的相对效力因区域和气候而异：在北美东部，减排降低高臭氧日频率的百分比在历史和高变暖情景下相近（~40%）；在西欧中部，高变暖下该百分比（-18%）小于历史时期（-25%）；而在东北亚，高变暖下的百分比（-16%）则大于历史时期（-13%）。这综合表明，在东北亚，未来NOx减排的“效力”可能相对增强，但在西欧中部，其效力可能相对减弱。
第五，潜在驱动机制的分析揭示了造成上述区域差异的关键竞争过程： * 北美东部和西欧中部：由于基准人为NOx排放较低，土壤NOx排放占总NOx的比例相对可观。在高变暖下，土壤NOx排放大幅增加（西欧中部增加65%），稀释了人为NOx排放的占比，并改变了NOx与ROx自由基的平衡，从而削弱了臭氧对人为NOx减排的敏感性。这是导致这两个区域δO3-m10%enox绝对值减小的主导因素。 * 东北亚：由于基准人为NOx排放很高，土壤NOx排放的占比很小（从5%升至7%），其相对影响有限。相反，高温高湿条件下HOx自由基产量的显著增加成为主导机制。增加的HOx会拓宽NOx促进臭氧生成的光化学区间（即提高NOx敏感区的上限），并提升单位NOx的臭氧生成效率，尤其在原先接近NOx饱和的区域效果更明显。这导致该区域臭氧对NOx减排的敏感性增强，δO3-m10%enox负值更大。 * 其他因素：PAN热分解加速、生物源异戊二烯排放的变化（受温度升高和CO2抑制的共同作用）等过程也参与竞争，使得响应变化并非总是随变暖程度单调变化。
结论与意义 本研究得出核心结论：气候变化对NOx减排臭氧效益的影响具有强烈的区域依赖性，这主要取决于各区域基准的人为NOx排放水平及由此决定的主导化学机制。 * 在北美东部和西欧中部这类已过渡到区域尺度NOx敏感状态的地区，高温驱动的土壤NOx排放增加是主要因素，它可能削弱未来NOx减排的边际效益，部分抵消因人为NOx长期下降而预期的臭氧控制收益。 * 在东北亚这类目前仍存在NOx不敏感区域的排放密集区，气候变化（特别是HOx产量增加）可能使臭氧化学向更有利于NOx减排的方向转变，降低减排导致臭氧升高的风险，并可能增强减排效力。然而，必须注意到，高变暖下该区域的基线臭氧水平本身是上升的，控制挑战依然严峻。
本研究的科学价值在于：首次系统性地运用集合模拟方法，量化了气候内部变率及不同变暖程度下，臭氧对NOx减排响应的不确定性范围，并揭示了其背后的关键化学-气候反馈机制。它强调了在进行未来空气质量影响评估时，必须超越对平均变化的关注，充分考虑概率分布、极端事件频率以及内部变率的遮蔽效应。
其应用价值在于：为空气质量管理者和政策制定者提供了更全面的未来情景分析工具和视角。研究结果表明，未来NOx控制策略的有效性评估需要结合具体区域的气候变化路径和排放基线。特别是在制定长期空气质量达标规划时，需考虑气候变化可能带来的“气候惩罚”（如西欧中部）或“气候协同”（如东北亚）效应。
研究亮点 1. 方法创新：采用多成员集合的气候-化学耦合模拟框架，首次将气候内部变率的定量影响纳入臭氧减排响应研究，显著提升了未来预估的稳健性和可靠性。 2. 视角全面：不仅分析平均响应，更深入探讨了响应概率分布、臭氧负效益天数、高臭氧日频率等与空气质量标准密切相关的统计指标。 3. 机制深入：通过诊断土壤NOx、HOx、PAN等多个竞争性过程，清晰阐明了不同区域响应差异背后的主导驱动因素，建立了“排放基线-主导机制-响应变化”的逻辑链条。 4. 区域对比鲜明：研究揭示了北美/欧洲与东北亚在未来气候变化影响下截然不同的臭氧化学演变轨迹，为针对性的区域空气质量管理提供了重要科学启示。
其他有价值的内容 研究还探讨了模型不确定性的影响，包括水平与垂直分辨率对模拟NOx敏感性的影响（尤其是在污染地区对NO滴定作用的模拟），以及CO2对异戊二烯排放抑制作用参数化带来的不确定性。这些分析增强了研究结论的透明度。作者指出，未来工作需进一步改进生物圈-大气相互作用过程（如土壤NOx排放的温度依赖性）的表述，并深入理解PAN在不同化学环境下对臭氧生成的影响。此外，除了气象变率，未来人为排放、土地利用变化等因素的不确定性同样重要，特别是甲烷变化的影响需要单独评估。