2060年代气候与排放情景下全球尺度O3污染的动力学研究

《Atmospheric Environment》：Kinetic study of global-scale O 3 pollution under 2060s climate and emission scenarios

　　

开展了全球尺度的空气质量模拟，以从大气化学角度阐明未来臭氧（O3）污染状况。开发了一个零维动力学模型，以阐明HOX循环增强动力学的重要性。温度升高增强了HOX循环，部分抵消了O3前体物减少对O3生成动力学的影响。

论文解读

一、研究背景与问题提出

臭氧（O₃）污染是全球性环境问题，对生态环境、人类和动物健康以及农作物生产均产生负面影响。世界卫生组织报告指出，地面O₃每年导致约50万人过早死亡，凸显了O₃污染治理的紧迫性。同时，O₃作为一种气候强迫因子，贡献了0.4 W m^-2的辐射强迫，约占全球增温总量的10%。

在对流层中，O₃通过二氧化氮（NO₂）的光解反应生成，并受到挥发性有机化合物（VOCs）和一氧化碳（CO）等前体物的增强。因此，减少氮氧化物（NO_X）和VOCs等前体物排放被认为是缓解O₃污染的有效途径。政策制定者正致力于建立净零碳社会以应对全球变暖，可再生能源、碳捕集技术和电动汽车的引入预计将同时减少CO₂和一次空气污染物的排放。

然而，已有研究多基于气候模型输出对O₃前体物与气候变化之间的关系进行总体分析，忽略了O₃生成过程中涉及的具体反应动力学。反应速率系数依赖于温度，气候驱动的温度变化预计会影响HO_X化学中基元反应的动力学。例如，OH + VOC反应通常随温度升高而加速，而HO₂ + HO₂等自由基-自由基反应因无能垒活化能而呈现负温度依赖性。因此，温度升高可能对O₃生成的反应动力学产生正反两方面的影响，但以往研究尚未深入探讨这些动力学效应。基于此，本研究旨在通过聚焦反应动力学，揭示气候变化情景下O₃生成化学的内在机制。

二、研究内容与主要结论

研究人员利用全球化学传输模型（CTM）GEOS-Chem v14.4.3，在2° × 2.5°空间分辨率和72个垂直层设置下，模拟了2017–2019年（当前）和2059–2061年（未来）两个时期的O₃及相关大气物种浓度。研究选取了17个分析站点，包括15个城市站点（东京、北京、莫斯科、纽约、巴黎、德里、达卡、雅加达、卡拉奇、拉各斯、墨西哥城、珠江三角洲、布宜诺斯艾利斯、约翰内斯堡、圣保罗）和2个偏远站点（亚马逊、大维多利亚沙漠）。未来气候情景采用SSP1-1.9和SSP2-4.5，分别对应实现1.5 °C目标的强CO₂限制情景和中等气候行动情景。通过设置8种CTM模拟情景（EcMc、EfMf、EfMc、EcMf），分别隔离排放变化和气候变化对空气质量的影响。

研究结果表明，在SSP1-1.9和SSP2-4.5情景下，北半球大部分地区和赤道站点在气候变化和减排共同作用下，地面O₃浓度呈现正增长。排放减少情景下O₃的贡献较低，大部分O₃增加归因于气候变化。但在美国东部和中国东部等站点呈现不同特征：美国东部在SSP1-1.9情景下人为排放大幅减少，超过气候变化效应，O₃浓度下降；中国东部在两种情景下O₃增加可归因于NO_X减少引起的NO滴定效应增强。北极地区O₃大幅增加主要归因于从大陆地区输送的O₃增加。

因子分析显示，在SSP1-1.9和SSP2-4.5情景下，大多数站点的O₃增加与绝对湿度（ΔAH）、短波辐射（ΔSW）和温度（ΔT）中至少一个参数的增加相关。排放减少情景下，NO_X、VOC和CO浓度的降低在大多数城市和偏远站点均有效。但在北京站点，NO_X减少速率远高于其他城市，通过抑制NO + O₃ → NO₂ + O₂反应导致O₃增加。

三、关键技术方法

本研究采用以下关键技术方法：其一，全球化学传输模型GEOS-Chem v14.4.3（含KPP fullchem标准机制，涵盖300余种物种和1500余个化学反应），用于模拟O₃及相关物种的时空分布；其二，气象输入采用MERRA2再分析资料（当前时期）和MRI-ESM2.0地球系统模式输出（未来时期），并通过季节平均偏差校正方法消除GCM与再分析资料之间的系统偏差；其三，排放清单采用CMIP6的SSP1-1.9和SSP2-4.5人为排放情景，生物质燃烧排放取自GFED数据库，BVOC排放采用MEGAN v2.1在线模型计算；其四，构建了基于HO_X化学（包含R1–R9九个关键反应的零维动力学模型，通过稳态近似和牛顿-拉夫森方法求解OH、HO₂和RO₂浓度，进而解析计算O₃生成速率P_O3。

四、研究结果

3.1 2017–2019与2059–2061年间地面O₃的变化

在SSP1-1.9和SSP2-4.5情景下，北半球大部分地区和赤道站点在气候变化和排放减少共同作用下，地面O₃呈现正增长。排放减少情景（EfMc）和气候变化情景（EcMf）的对比分析表明，大部分O₃增加归因于气候变化，而非排放变化。美国东部在SSP1-1.9情景下因人为排放大幅减少导致O₃下降；中国东部因NO_X减少引起的NO滴定效应增强导致O₃增加。北极地区O₃增加主要归因于从大陆地区输送的O₃增加。

3.2 地面O₃变化的因子分析

在SSP1-1.9和SSP2-4.5情景下，大多数站点O₃增加与ΔAH、ΔSW和ΔT中至少一个参数的增加相关，表明气候变化引起的气象参数变化影响O₃变化。莫斯科站点在两种情景下均显示O₃下降，可能与风速、风向和边界层高度等其他气象因素变化有关。排放减少情景下，NO_X、VOC和CO浓度降低在大多数站点有效，但北京站点因NO_X减少速率远高于其他城市，通过NO滴定效应导致O₃增加。研究结果与Turnock等（2018，2020）和Allen等（2020）的预测在O₃减少量级上基本一致。

3.3 地面O₃变化的动力学分析

基于HO_X化学的零维动力学模型，研究人员推导了O₃生成速率P_O3的解析表达式。模型验证显示，动力学模型估算的P_O3日变化特征与CTM模拟结果高度吻合，尽管存在约1.7倍的高估，这主要源于简化机制中R7反应（RO₂ + RO₂）对HO₂产率的假设。

3.3.3 主要空气污染物、环境温度和HO_X生成速率对地面O₃生成速率的相对敏感性

相对敏感性分析表明，城市站点的P_O3对VOC和NO_X浓度均表现出较高敏感性，接近过渡区特征。P_O3对CO的相对敏感性与VOC和NO_X相当，表明相同比例的CO减排与VOC和NO_X减排在缓解光化学O₃生成方面效果相当。P_O3对温度、P_OH和P_XO2的敏感性高于对VOC和NO_X的敏感性，表明HO_X循环增强动力学是O₃生成增加的重要贡献因素。

3.3.4 O₃前体物浓度和温度变化对P_O3的影响

绝对敏感性分析显示，在SSP1-1.9情景下，大多数站点的P_O3对AVOC呈负敏感性，表明AVOC减排对缓解光化学O₃生成有效。在SSP2-4.5情景下，纽约、莫斯科、拉各斯、德里、达卡和珠江三角洲等站点的P_O3随AVOC和BVOC增加而增加，这归因于BVOC排放及其反应产物的增加。P_O3对NO_X在除北京外的大多数站点呈负敏感性。温度升高对P_O3的贡献与前体物浓度变化相当，表明HO_X化学速率系数变化的重要性与前体物浓度变化相当。

3.3.5 P_O3对HO_X化学速率系数的敏感性分析

对R1–R9九个反应的速率系数敏感性分析表明，SSP2-4.5情景下纽约、莫斯科、北京、东京、德里、珠江三角洲和圣保罗等站点对HO_X化学所有速率系数均表现出强烈敏感性。k₁（OH + VOC）、k₃（HO₂ + NO）、k₄（RO₂ + NO）和k₉（OH + CO）具有正活化能，随温度升高而增大，增强O₃生成；k₂（OH + NO₂）、k₆（HO₂ + RO₂）、k₇（RO₂ + RO₂）和k₈（OH + OH）为无能垒反应，随温度升高而减小，削弱终止反应，整体上对P_O3产生正响应。

3.3.6 研究局限性

敏感性分析聚焦于P_O3而非实际O₃浓度。地面O₃浓度不仅受HO_X化学控制，还受沉降、输送等物理过程影响。尽管如此，研究结果提供了前体物、温度、P_OH和P_XO2如何基于HO_X化学动力学影响O₃生成效率的基础认识。

五、讨论与结论

研究结果表明，虽然人为排放减少在SSP1-1.9情景下的大多数站点和SSP2-4.5情景的部分站点导致O₃下降，但气候变化效应（特别是BVOC排放增加）抵消了这些减少，导致全球尺度O₃增加。敏感性分析显示，环境温度升高对P_O3的驱动作用与VOC和NO_X排放变化相当，其机制在于温度增强HO_X循环动力学。具体而言，OH + VOC和RO₂ + NO反应具有正能垒，温度升高增强这些反应和HO_X循环；而HO₂ + HO₂、HO₂ + RO₂和RO₂ + RO₂等终止反应具有负能垒，温度升高削弱这些反应路径，同样贡献于P_O3增加。

本研究证明了将动力学模型与CTM结果耦合以阐明气候变化对O₃污染影响的有效性，从物理和化学角度深化了对未来空气污染问题的理解，为O₃污染治理策略的制定提供了动力学层面的科学依据。