多环芳烃（PAHs）在CHASER化学-气候模式中的全球模拟

《Atmosphere》：Global Simulation of Polycyclic Aromatic Hydrocarbons (PAHs) with the CHASER Chemistry–Climate Model

　　

为厘清全球尺度气-粒分配热力学不确定性与持久性有机污染物传输机制，研究人员将多环芳烃（PAHs）耦合进CHASER V4.0化学-气候模式，首次在该框架内实现萘（NAP）、菲（PHE）、芘（PYR）与苯并[a]芘（BaP）的在线全球模拟以研究其大气归趋。在扩展的Dachs–Eisenreich双模方案中，烟炱吸附（KSA(T)）与有机质吸收（KOA(T)）相耦合。耦合模式在城市、背景与北极站点评估中较好捕捉全球空间梯度（对数10尺度合并空间R=0.78；剔除严重子网格热点的BaP空间FAC2=80%）与季节循环（R=0.62–0.93）。KSA(T)方案修正了基线KOA(T)方案在冷高纬对颗粒态分数（φ）的低估（极地φBaP最高提升约0.6）。敏感性试验表明，元素碳质量分数（fEC）升高增强烟炱吸附与颗粒相对氧化的屏蔽，使BaP平均寿命延至1.72天（轻同系物为0.57–0.72天）并维持跨洋流出。该框架为评估跨界有毒物传输、气候敏感健康风险及污染减缓政策提供重要工具。

研究背景方面，多环芳烃（PAHs，半挥发性持久性有机污染物）源于含碳物质不完全燃烧，具致癌致畸致突变效应，可经“全球蒸馏”输运至极地，而变暖导致的冰雪消融会引发二次排放，使全球循环进一步复杂化。气-粒分配是决定PAH环境归趋的核心过程，直接调控氧化降解与沉降效率；颗粒相组分受有机气溶胶“屏蔽效应”保护而延长大气寿命。既有研究存在三方面问题：不同模式对吸收/吸附温度响应的刻画离散，低温下分配系数（Kp）可差一个量级，造成北极地表浓度低估达量级水平；全耦合模拟中气溶胶背景与目标信号共变，难以解耦有机质（f_OM）与元素碳（f_EC）的独立贡献；氧化次级产物毒性更高但全球显式传输困难。因此研究人员在CHASER V4.0（Chemical Atmospheric General Circulation Model for Study of Atmospheric Environment and Radiative Forcing，耦合MIROC5的大气化学-气候模式，T42水平分辨率约2.8°×2.8°，36层σ-p混合垂直坐标，在线耦合SPRINTARS气溶胶模块）中新建PAH分支，选取NAP、PHE、PYR、BaP四种代表性物种，采用扩展Dachs–Eisenreich双模分配（有机质吸收K_OA(T)叠加烟炱吸附K_SA(T)），驱动2004年GEMS-PAH底向上排放清单（0.1°重采样至模型网格，分生物质燃烧、化石燃料与工业、交通三类），气相氧化考虑OH、O₃、NO₂二级反应，光解以甲基丙烯醛（MACR）为代理并乘环数相关修正因子，干沉降用阻力模型、湿沉降按有机碳类体扫并引入低温（<268 K）冰云增强因子，借此开展在线全球模拟与多站点评估。论文发表于《Atmosphere》。

关键技术方法上，研究人员以CHASER V4.0自由运行模式为宿主，水平T42、垂直36层，在线耦合气象与化学；排放采用2004年GEMS-PAH全球清单并按三大宏观部门聚合；化学损失用OH/O₃/NO₂速率常数与MACR代理光解；气-粒分配用Dachs–Eisenreich双模公式，K_OA(T)依van ’t Hoff关系、K_SA(T)取同温度敏感系数（b=1）并做离线b=0.5–1.5敏感性；设置BASE/A1/B1/A2/B2五组f_OM/f_EC固定场景解耦组分贡献；验证数据来自城市（北京、广州）、区域背景（Zingst、Ko?etice、Aspvreten、Lahemaa）与北极（Alert）站点2003–2008年观测。

研究结果部分，3.1节大气浓度、季节变化与验证显示，2004年全球年均表面浓度NAP 2.00、PHE 0.42、PYR 0.17、BaP 0.05 ng m^?3，北低纬高于南半球，东亚、印度、东南亚、非洲为高峰；BaP寿命1.72天、轻同系物0.57–0.72天，颗粒相屏蔽抵消快速物理移除从而维持跨太平洋与大西洋羽流。季节上冬高夏低，NH中纬受采暖与交通排放及夏季OH光化学清除主导，南亚早春受残旱季生物质燃烧与弱季风清除控制，亚马逊9月单峰、非洲多月中峰对应火季-雨季交替。与观测对比，欧洲背景站BaP月序列R=0.62–0.74、RMSE极低；北京站R=0.93但冬季峰值因2.8°网格稀释低估（NMB?46%），5–10月偏高反映城市强OH损耗与季风湿清除表征不足；全局空间log₁₀R=0.78，剔广州、Agra后BaP空间FAC2=80%。3.2节气-粒分配表明，仅K_OA(T)方案在极地φBaP偏低，加入K_SA(T)后70°N/S出现正差峰，极地φBaP提升达0.6；半挥发PHE、PYR的φ对K_SA(T)相对增幅最大，NAP仍气态主导；五场景敏感性中f_EC升高放大φ对组分响应，东亚/欧洲/北美中纬φ跳跃显著，证明K_SA(T)引入使模式对f_OM/f_EC偏差更敏感并强化向北极的颗粒载带传输。3.3节不确定性指出：2004年排放仅作高排放基准，近年中国减排使东亚绝对值偏低；瞬时平衡假设无法复现广州夏季新鲜烟炱非平衡包裹致BaP高估挥发；b=1假设对BaP冷区φ无影响、对PHE/PYR存参数不确定；北京f_EC=0.10敏感试验改善匹配属网格稀释补偿而非组分真实值；观测跨年拼接与相位定义差异限制超级城市评估。

讨论与结论翻译：研究人员将专用PAH模块植入CHASER V4.0，建立涵盖全球排放、气相氧化、沉降与温度依赖双模分配的在线系统，得出（1）K_SA(T)方案修正高纬φ低估，烟炱吸附主导冷区低挥发PAH分配；（2）敏感试验证实烟炱吸附增强φ对f_OM/f_EC的响应；（3）高污染羽流中增强的屏蔽效应使BaP寿命延至1.72天，支撑跨洋传输；（4）尽管2.8°分辨率致城市峰值代表误差，模式仍再现全球空间（R=0.78，BaP FAC2=80%）与季节（R=0.62–0.93）规律。该双模框架为气候变化下跨界有毒物传输、人群暴露与POPs削减政策评估提供预报工具，后续将扩展含氧PAH次级衍生物在线传输。