Environmental, green and energy chemistry

Atmospheric chemistry models

How box and three-dimensional models combine atmospheric transport, emissions, chemical mechanisms and observations to explain pollutant and climate-relevant species.

1 Lower secondary 2 Upper secondary 3 Undergraduate 4 Advanced 5 Research

IntuitionIntuition: atmosphere as a moving chemical ledger

A model tracks where molecules go, how air carries and mixes them, and how chemistry makes or removes them. A simple box is a ledger; a three-dimensional model is that ledger repeated across a weather-driven grid.

Select a mode. In the box view, trace emissions, chemistry, deposition and exchange; in the 3-D view, follow surface emissions and transport among grid cells and vertical layers.

SchoolSchool level: sources, sinks and lifetime

Definition: Chemical lifetime

For first-order removal at rate constant $k$ , the e-folding lifetime is $au=1/k$ . It is a characteristic time, not necessarily the age of every molecule.

\frac{dC}{dt}=P-kC\qquad C_{ss}=\frac{P}{k}=P\tau

At steady state, production balances loss. If a well-mixed box contains 20 ppb of a species and its first-order loss is 0.10 h⁻¹, its lifetime is 10 h and the production needed to maintain that abundance is 2 ppb h⁻¹.

Example: A two-species box budget

Species A is emitted at 4 ppb h⁻¹ and removed at 0.20 h⁻¹. A reacts to form B at 0.10 h⁻¹; B is removed at 0.25 h⁻¹. Find steady-state A and B.

Solution

A balance gives $A_{ss}=4/(0.20+0.10)=13.3$ ppb. B is produced at $0.10A=1.33$ ppb h⁻¹, so $B_{ss}=1.33/0.25=5.33$ ppb. Conversion is a sink for A and a source for B; omitting it would give the wrong A.

UndergraduateUniversity: equations and photochemistry

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

Here $C$ is a mixing ratio or concentration, $\mathbf u$ wind, $K$ an effective turbulent diffusivity, $P-L$ chemical production minus loss, and $E-D$ emissions minus deposition. Units and coordinate conventions must be consistent; boundary conditions close the transport problem.

A photochemical box model solves coupled reaction rates under prescribed sunlight, temperature and pressure. The Master Chemical Mechanism (MCM) represents explicit near-explicit VOC oxidation chemistry; its large reaction network makes computation demanding.

Example: Leighton photostationary state

For $j_{NO2}=0.008$ s⁻¹, $k_{NO+O3}=1.8\times10^{-14}$ cm³ molecule⁻¹ s⁻¹, $[NO2]=5$ ppb and $[NO]=2$ ppb, find ozone under the Leighton relation.

Solution

At photostationarity, $j[NO2]=k[NO][O3]$ , so $[O3]/[NO2]=j/(k[NO])$ . At 298 K and 1 atm, air has ≈ $2.46\times10^{19}$ molecules cm⁻³; 2 ppb NO is $4.92\times10^{10}$ molecules cm⁻³. Thus $[O3]/[NO2]=9.0$ , and ozone ≈45 ppb. Peroxy-radical conversion of NO to NO2 breaks this simple relationship.

AdvancedAdvanced: coupled solvers and global models

Atmospheric mechanisms are stiff: radical reactions can be extremely fast while reservoir species persist for years. Explicit time steps restricted by the fastest reactions are inefficient; implicit Euler is stable for linear decay, while Gear/BDF and Rosenbrock methods handle stiff systems using Jacobians. QSSA eliminates selected short-lived intermediates when its assumptions hold.

Operator splitting advances transport, chemistry, radiation and deposition in separate substeps. It is efficient but introduces splitting error; symmetric Strang splitting improves formal order when its operators and time scales permit. GEOS-Chem, CMAQ, WRF-Chem and CAM-chem combine meteorology, emissions, chemistry and removal with different architectures and scientific purposes.

Representative chemistry mechanisms
Name	Use
CB05	Condensed regional air-quality chemistry
SAPRC	Photochemical smog mechanisms; lumped VOC classes
MOZART	Global tropospheric and stratospheric chemistry
MCM	Detailed VOC oxidation, often box/chamber applications

Inventories provide gridded or sector totals for anthropogenic and natural emissions; temporal profiles, chemical speciation and plume injection add uncertainty. Dry deposition depends on surface exchange, while wet removal depends on cloud and precipitation processes. Model evaluation compares multiple observables, not just a single concentration.

ResearchResearch: inverse models and evaluation

Adjoint sensitivities and Bayesian inversion infer emissions or chemical parameters from observations. Satellite columns from OMI and TROPOMI constrain pollutants and methane, but retrieval averaging kernels, clouds, transport errors and prior assumptions must be represented. Data assimilation updates model states; inversion estimates uncertain causes or parameters.

Ozone production efficiency (OPE) relates ozone formed to NOx consumed over a defined air-mass history. It is diagnostic, not a universal constant: VOC regime, dilution, deposition and the chosen accounting interval matter. A model–measurement comparison should pair surface networks, aircraft or satellite data with matched sampling and uncertainty analysis.

Methane oxidation consumes OH; a methane increase can lower OH and lengthen methane lifetime, a positive chemical feedback. For the lifetime formula, the OH field is therefore not independent of methane. Coupled chemistry–climate model intercomparison (CCMI) experiments also help assess stratospheric ozone-hole recovery under scenarios; projections remain conditional on emissions and climate forcing. Numerically, with $k_{CH4+OH}\approx6.3\times10^{-15}$ cm³ molecule⁻¹ s⁻¹ and a mass-weighted mean $[OH]\approx10^6$ cm⁻³, $1/(k[OH])\approx5$ y; global budgets give ≈9 y because OH varies and soils and stratosphere also consume CH4.

References

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