基于贝叶斯统计的气溶胶油酸环境氧化过程建模

《ACS Environmental Au》：Modeling Ambient Oxidation of Aerosol Oleic Acid with Bayesian Statistics

　　

有机气溶胶（OA）的非均相氧化是关键的气象过程，但预测建模受限于缺乏环境衰减动力学数据及与环境变量定量关系未定义。研究人员以不饱和脂肪酸（uFAs）为模型OA化合物，基于中国上海两站点长期观测确定了油酸、反油酸与亚油酸的环境衰减速率。尽管温度、臭氧及初始uFA浓度和OA基质发生显著季节变化，衰减速率季节差异极小。分析揭示uFA衰减速率具有区别于以往观察的温度依赖性，且受其他环境参数复合影响。通过应用贝叶斯统计校准，研究人员推导出各温度区间专属的油酸衰减速率定量参数化方案。有机碳（OC）对反应速率具抑制作用，低温下尤甚；而反应物浓度影响随温度升高而增强。结果推进了uFA非均相氧化机理认知，并为真实大气条件下改进OA老化模型参数化提供了鲁棒且受观测约束的框架。

研究背景方面，有机气溶胶（OA）占大气细颗粒物（PM_2.5）相当比例，其非均相氧化主导化学演化与归宿。当前化学传输模式常将OA简化为均匀混合液体，忽视相态与黏度对氧化剂扩散及反应摄取动力学的影响，且缺乏高氧有机分子氧化显式机制，导致模拟偏差。以不饱和脂肪酸（uFAs）如油酸（C18:1）为模型体系虽已有实验室臭氧化动力学，但真实大气衰减受气溶胶基质与环境变量显著调制，实验室外推存在不确定性。为此，研究人员于2019年12月至2020年9月及2022年10月至2023年2月在上海城市站点（SAES）与郊外淀山湖（DSL）开展野外观测，利用贝叶斯框架建立油酸衰减速率与环境变量的定量关系。

主要关键技术方法包括：基于TAG-GC/MS双小时分辨率测定C16:1、C18:1、t-C18:1、C18:2等uFAs及硬脂酸（C18:0）参照物浓度；以晚餐后夜间时段C_uFA/C_18:0比值指数拟合获取拟一级衰减速率k；采用扩展阿伦尼乌斯方程 k = A[O₃]^α[uFA]^β[OC]^γexp(?E_a/RT)+ε 描述温度与浓度依赖；以贝叶斯推断结合MCMC（emcee包，50000迭代）校准A、α、β、γ及σ，先验按物理合理性设定，后验由似然与先验乘积采样获得。

研究结果部分，3.1节气溶胶脂肪酸总体特征显示uFAs占脂肪酸总量约29.2%，餐后夜间归一化浓度持续下降指示连续化学降解。3.2节uFAs衰减动力学表明C18:1平均k=0.21±0.07 h^?1，C18:2（0.27±0.11 h^?1）高于C18:1，t-C18:1（0.14±0.06 h^?1）低于之；初级阿伦尼乌斯图按温度区间（TR₁:1–4℃、TR₂:4–11℃、TR₃:11–21℃、TR₄:21–26℃、TR₅:26–31℃）分段呈不同斜率，TR₁无温度依赖，TR₂与TR₅线性最佳，证实相态转变调控动力学。3.3节扩展阿伦尼乌斯贝叶斯参数化中，各温区独立拟合优度优于全数据混合；TR₂表观活化能E_a,p=50.0 kJ·mol^?1，TR₅为100.9 kJ·mol^?1，中间温区插值获得；后验均值显示A随温区增近9个量级，γ在TR₂为?0.206（OC抑制最强），α、β在TR₃与TR₅绝对值最大，TR₄各指数近零；扩展模型R_p较初级阿伦尼乌斯由0.35升至0.68–0.81，BIC降低，TR₅预测最佳（R_p=0.81）。

讨论部分指出，低温区（TR₂）油酸偏固态，分子扩散限制使OC抑制显著而反应物级次趋零；中温区（TR₃）液态转变令氧化剂级次α与反应物级次β超越γ；高温区（TR₅）黏度下降、氧化剂驱动为主，α主导；TR₄处于扩散—反应控制过渡，残差σ偏高暗示相态与水含量等未控因素。模型虽未涵盖黏度与液态水，但已捕获主要变量协同效应。

结论部分翻译：研究人员基于上海两站10个月野外测量确定三种uFAs环境衰减速率，将贝叶斯建模用于显式纳入臭氧浓度、温度、初始uFA浓度与OC浓度的扩展阿伦尼乌斯方程，推导油酸衰减速率与环境变量定量关系；该方法通过量化参数不确定性与交互作用，灵活描述uFA衰减动态，较初级阿伦尼乌斯方程显著提升预测性能，尤其低高温区；校准参数揭示各温区反应受不同环境变量主导，OC在低温抑制显著、反应物浓度影响随温升增强，从而推进uFA非均相氧化机理理解，为真实大气OA老化模型参数化提供观测约束框架。论文发表于《ACS Environmental Au》。