Chemistry Labs

环境、绿色与能源化学

大气化学模型

箱体与三维模型如何结合大气输送、排放、化学机理和观测,解释污染物与气候相关物种。

直觉直觉:流动的化学收支账本

模型追踪分子的位置、空气输送与混合,以及化学反应的生成和去除。简单箱体是收支账本;三维模型则在受天气驱动的网格上重复计算。

选择一种模式。在箱体视图中追踪排放、化学、沉降和交换;在三维视图中观察地表排放及网格和垂直层之间的输送。

中学阶段中学阶段:来源、汇与寿命

定义: 化学寿命

对于速率常数为 kk 的一级去除过程,e 倍衰减寿命为 au=1/k au=1/k。这是特征时间,不一定是每个分子的实际年龄。

dCdt=P−kCCss=Pk=Pτ\frac{dC}{dt}=P-kC\qquad C_{ss}=\frac{P}{k}=P\tau

稳态时生成与损失相平衡。若充分混合的箱体中某物质浓度为 20 ppb,一级损失率为 0.10 h⁻¹,则寿命为 10 h,维持该浓度所需生成速率为 2 ppb h⁻¹。

例题: 双物种箱体收支

A 以 4 ppb h⁻¹ 排放,以 0.20 h⁻¹ 去除;A 以 0.10 h⁻¹ 转化为 B,B 以 0.25 h⁻¹ 去除。求稳态 A、B 浓度。

解答

由收支得 Ass=4/(0.20+0.10)=13.3A_{ss}=4/(0.20+0.10)=13.3 ppb。B 的生成速率为 0.10A=1.330.10A=1.33 ppb h⁻¹,因此 Bss=1.33/0.25=5.33B_{ss}=1.33/0.25=5.33 ppb。转化是 A 的汇、B 的源;忽略它会算错 A。

大学大学:方程与光化学

∂C∂t=−∇⋅(uC)+∇⋅(K∇C)+P−L+E−D\frac{\partial C}{\partial t}=-\nabla\cdot(\mathbf u C)+\nabla\cdot(K\nabla C)+P-L+E-D

其中 CC 为混合比或浓度,u\mathbf u 为风,KK 为有效湍流扩散系数,P−LP-L 为化学生成减损失,E−DE-D 为排放减沉降。单位和坐标约定必须一致,边界条件使输送问题闭合。

光化学箱体模型在给定日照、温度和压力下求解耦合反应速率。Master Chemical Mechanism(MCM)详细表示 VOC 氧化化学;庞大的反应网络计算量很大。

例题: Leighton 光稳态

当 jNO2=0.008j_{NO2}=0.008 s⁻¹、kNO+O3=1.8×10−14k_{NO+O3}=1.8\times10^{-14} cm³ molecule⁻¹ s⁻¹、[NO2]=5[NO2]=5 ppb、[NO]=2[NO]=2 ppb 时,用 Leighton 关系求臭氧。

解答

光稳态下 j[NO2]=k[NO][O3]j[NO2]=k[NO][O3],故 [O3]/[NO2]=j/(k[NO])[O3]/[NO2]=j/(k[NO])。298 K、1 atm 的空气约含 2.46×10192.46\times10^{19} molecules cm⁻³,2 ppb 的 NO 为 4.92×10104.92\times10^{10} molecules cm⁻³。比值因此为 9.0,臭氧约 45 ppb。过氧自由基将 NO 转化为 NO2 会破坏这一简单关系。

进阶进阶:耦合求解器与全球模型

大气反应机理是刚性系统:自由基反应可能极快,而储库物种可持续多年。显式步长受最快反应限制,效率低;隐式 Euler 对线性衰减稳定,Gear/BDF 和 Rosenbrock 方法利用雅可比矩阵处理刚性系统。满足假设时,QSSA 可消去选定的短寿命中间体。

算子分裂将输送、化学、辐射和沉降分成独立子步推进,效率高但会产生分裂误差;条件适合时,对称 Strang 分裂可提高形式阶数。GEOS-Chem、CMAQ、WRF-Chem 和 CAM-chem 以不同架构和科学目标耦合气象、排放、化学与去除。

代表性化学机理
名称用途
CB05简化的区域空气质量化学
SAPRC光化学烟雾机理;VOC 分类归并
MOZART全球对流层和平流层化学
MCM详细 VOC 氧化,常用于箱体或烟雾箱

排放清单提供人为和自然源的网格或部门总量;时间分布、化学物种分配和烟羽注入高度都会带来不确定性。干沉降取决于地表交换,湿清除取决于云和降水。模型评估应比较多种观测量,而非仅看一个浓度。

科研研究:反演与评估

伴随灵敏度和贝叶斯反演可由观测推断排放或化学参数。OMI、TROPOMI 卫星柱浓度约束污染物和甲烷,但必须考虑反演平均核、云、输送误差和先验假设。数据同化更新模型状态;反演估计不确定的成因或参数。

臭氧生成效率(OPE)描述特定气团历史中生成臭氧与消耗 NOx 的关系,是诊断量而非普适常数;VOC 条件、稀释、沉降及统计时段都会影响它。模型—观测比较应匹配地面网络、飞机或卫星采样,并分析不确定性。

甲烷氧化消耗 OH;甲烷增加可能降低 OH、延长甲烷寿命,形成正化学反馈。因此寿命公式中的 OH 场并非独立于甲烷。化学—气候模型相互比较(CCMI)也用于评估情景下平流层臭氧洞恢复,但预测取决于排放和气候强迫。 数值上,取 kCH4+OH≈6.3×10−15k_{CH4+OH}\approx6.3\times10^{-15} cm³ molecule⁻¹ s⁻¹、质量加权平均 [OH]≈106[OH]\approx10^6 cm⁻³,则 1/(k[OH])≈51/(k[OH])\approx5 年;全球收支给出约 9 年,因为 OH 时空变化且土壤和平流层也消耗 CH4。

参考文献

  • Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, 3rd ed. · J. H. Seinfeld, S. N. Pandis, 2016
  • Introduction to Atmospheric Chemistry · D. J. Jacob, 1999
  • Global modeling of tropospheric chemistry with assimilated meteorology: Model description and evaluation · H. Bey et al., 2001
  • Protocol for the development of the Master Chemical Mechanism, MCM v3 (Part A): tropospheric degradation of non-aromatic volatile organic compounds · S. M. Saunders, M. E. Jenkin, R. G. Derwent, M. J. Pilling, 2003
  • Modeling of Atmospheric Chemistry · G. Brasseur, D. J. Jacob, 2017