The Community Earth System Model version 2 (CESM2) includes a detailed representation of chemistry throughout the atmosphere in the Community Atmosphere Model with chemistry and Whole Atmosphere ...Community Climate Model configurations. These model configurations use the Model for Ozone and Related chemical Tracers (MOZART) family of chemical mechanisms, covering the troposphere, stratosphere, mesosphere, and lower thermosphere. The new MOZART tropospheric chemistry scheme (T1) has a number of updates over the previous version (MOZART‐4) in CESM, including improvements to the oxidation of isoprene and terpenes, organic nitrate speciation, and aromatic speciation and oxidation and thus improved representation of ozone and secondary organic aerosol precursors. An evaluation of the present‐day simulations of CESM2 being provided for Climate Model Intercomparison Project round 6 (CMIP6) is presented. These simulations, using the anthropogenic and biomass burning emissions from the inventories specified for CMIP6, as well as online calculation of emissions of biogenic compounds, lightning NO, dust, and sea salt, indicate an underestimate of anthropogenic emissions of a variety of compounds, including carbon monoxide and hydrocarbons. The simulation of surface ozone in the southeast United States is improved over previous model versions, largely due to the improved representation of reactive nitrogen and organic nitrate compounds resulting in a lower ozone production rate than in CESM1 but still overestimates observations in summer. The simulation of tropospheric ozone agrees well with ozonesonde observations in many parts of the globe. The comparison of NOx and PAN to aircraft observations indicates the model simulates the nitrogen budget well.
Plain Language Summary
The set of chemical reactions for tropospheric chemistry used in the Community Earth System Model version 2 (CESM2) has been updated significantly over CESM1 in the Community Atmosphere Model with chemistry (CAM‐chem) and Whole Atmosphere Community Climate Model (WACCM) configurations. The emissions used for the CESM2 simulations are documented here, with anthropogenic and biomass burning emissions based on the specified inventories for Climate Model Intercomparison Project 6 (CMIP6), and emissions of biogenic compounds, lightning NO, dust, and sea salt are calculated online and dependent on the simulated meteorology. Evaluation of the CAM‐chem and WACCM configurations of CESM2 with observations indicate an underestimate of anthropogenic emissions of a variety of compounds, including carbon monoxide and hydrocarbons. The updated chemistry leads to an improvement in the simulation of tropospheric ozone.
Key Points
This paper fully documents the significant updates to the chemistry mechanisms in version 2 of the Community Earth System Model
The new tropospheric chemistry scheme improves representation of isoprene oxidation as well as other ozone precursors over earlier versions
The simulation of tropospheric ozone is improved in comparison to observations
Hydroxyl radical (OH) oxidation of toluene produces ring-retaining products: cresol and benzaldehyde, and ring-opening products: bicyclic intermediate compounds and epoxides. Here, first- and ...later-generation OH oxidation products from cresol and benzaldehyde are identified in laboratory chamber experiments. For benzaldehyde, first-generation ring-retaining products are identified, but later-generation products are not detected. For cresol, low-volatility (saturation mass concentration, C* ∼ 3.5 × 104 − 7.7 × 10−3 µg m−3), first- and later-generation ring-retaining products are identified. Subsequent OH addition to the aromatic ring of o-cresol leads to compounds such as hydroxy, dihydroxy, and trihydroxy methyl benzoquinones and dihydroxy, trihydroxy, tetrahydroxy, and pentahydroxy toluenes. These products are detected in the gas phase by chemical ionization mass spectrometry (CIMS) and in the particle phase using offline direct analysis in real-time mass spectrometry (DART-MS). Our data suggest that the yield of trihydroxy toluene from dihydroxy toluene is substantial. While an exact yield cannot be reported as authentic standards are unavailable, we find that a yield for trihydroxy toluene from dihydroxy toluene of ∼ 0.7 (equal to the reported yield of dihydroxy toluene from o-cresol; Olariu et al., 2002) is consistent with experimental results for o-cresol oxidation under low-NO conditions. These results suggest that even though the cresol pathway accounts for only ∼ 20 % of the oxidation products of toluene, it is the source of a significant fraction (∼ 20–40 %) of toluene secondary organic aerosol (SOA) due to the formation of low-volatility products.
Ozone is a greenhouse gas and air pollutant that is harmful to human health and plants. During the summer in the southeastern US, many regional and global models are biased high for surface ozone ...compared to observations. Past studies have suggested different solutions including the need for updates to model representation of clouds, chemistry, ozone deposition, and emissions of nitrogen oxides (NOx) or biogenic hydrocarbons. Here, due to the high biogenic emissions in the southeastern US, more comprehensive and updated isoprene and terpene chemistry is added into CESM/CAM-chem (Community Earth System Model/Community Atmosphere Model with full chemistry) to evaluate the impact of chemistry on simulated ozone. Comparisons of the model results with data collected during the Studies of Emissions Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys (SEAC4RS) field campaign and from the US EPA (Environmental Protection Agency) CASTNET (Clean Air Status and Trends Network) monitoring stations confirm the updated chemistry improves simulated surface ozone, ozone precursors, and NOx reservoir compounds. The isoprene and terpene chemistry updates reduce the bias in the daily maximum 8 h average (MDA8) surface ozone by up to 7 ppb. In the past, terpene oxidation in particular has been ignored or heavily reduced in chemical schemes used in many regional and global models, and this study demonstrates that comprehensive isoprene and terpene chemistry is needed to reduce surface ozone model biases. Sensitivity tests were performed in order to evaluate the impact of lingering uncertainties in isoprene and terpene oxidation on ozone. Results suggest that even though isoprene emissions are higher than terpene emissions in the southeastern US, remaining uncertainties in isoprene and terpene oxidation have similar impacts on ozone due to lower uncertainties in isoprene oxidation. Additionally, this study identifies the need for further constraints on the aerosol uptake of organic nitrates derived from isoprene and terpenes in order to reduce uncertainty in simulated ozone. Although the updates to isoprene and terpene chemistry greatly reduce the ozone bias in CAM-chem, a large bias remains. Evaluation against SEAC4RS field campaign results suggests future improvements to horizontal resolution and cloud parameterizations in CAM-chem may be particularly important for further reducing this bias.
We use a large laboratory, modeling, and field dataset to investigate the isoprene + O3 reaction, with the goal of better understanding the fates of the C1 and C4 Criegee intermediates in the ...atmosphere. Although ozonolysis can produce several distinct Criegee intermediates, the C1 stabilized Criegee (CH2OO, 61 ± 9%) is the only one observed to react bimolecularly. We suggest that the C4 Criegees have a low stabilization fraction and propose pathways for their decomposition. Both prompt and non-prompt reactions are important in the production of OH (28% ± 5%) and formaldehyde (81% ± 16%). The yields of unimolecular products (OH, formaldehyde, methacrolein (42 ± 6%) and methyl vinyl ketone (18 ± 6%)) are fairly insensitive to water, i.e., changes in yields in response to water vapor (≤4% absolute) are within the error of the analysis. We propose a comprehensive reaction mechanism that can be incorporated into atmospheric models, which reproduces laboratory data over a wide range of relative humidities. The mechanism proposes that CH2OO + H2O (k(H2O)∼ 1 × 10(-15) cm(3) molec(-1) s(-1)) yields 73% hydroxymethyl hydroperoxide (HMHP), 6% formaldehyde + H2O2, and 21% formic acid + H2O; and CH2OO + (H2O)2 (k(H2O)2∼ 1 × 10(-12) cm(3) molec(-1) s(-1)) yields 40% HMHP, 6% formaldehyde + H2O2, and 54% formic acid + H2O. Competitive rate determinations (kSO2/k(H2O)n=1,2∼ 2.2 (±0.3) × 10(4)) and field observations suggest that water vapor is a sink for greater than 98% of CH2OO in a Southeastern US forest, even during pollution episodes (SO2 ∼ 10 ppb). The importance of the CH2OO + (H2O)n reaction is demonstrated by high HMHP mixing ratios observed over the forest canopy. We find that CH2OO does not substantially affect the lifetime of SO2 or HCOOH in the Southeast US, e.g., CH2OO + SO2 reaction is a minor contribution (<6%) to sulfate formation. Extrapolating, these results imply that sulfate production by stabilized Criegees is likely unimportant in regions dominated by the reactivity of ozone with isoprene. In contrast, hydroperoxide, organic acid, and formaldehyde formation from isoprene ozonolysis in those areas may be significant.
Recent advances in our knowledge of the gas-phase oxidation of isoprene, the impact of chamber walls on secondary organic aerosol (SOA) mass yields, and aerosol measurement analysis techniques ...warrant reevaluating SOA yields from isoprene. In particular, SOA from isoprene oxidation under high-NOx conditions forms via two major pathways: (1) low-volatility nitrates and dinitrates (LV pathway) and (2) hydroxymethyl-methyl-α-lactone (HMML) reaction on a surface or the condensed phase of particles to form 2-methyl glyceric acid and its oligomers (2MGA pathway). These SOA production pathways respond differently to reaction conditions. Past chamber experiments generated SOA with varying contributions from these two unique pathways, leading to results that are difficult to interpret. This study examines the SOA yields from these two pathways independently, which improves the interpretation of previous results and provides further understanding of the relevance of chamber SOA yields to the atmosphere and regional or global modeling. Results suggest that low-volatility nitrates and dinitrates produce significantly more aerosol than previously thought; the experimentally measured SOA mass yield from the LV pathway is∼0.15. Sufficient seed surface area at the start of the reaction is needed to limit the effects of vapor wall losses of low-volatility compounds and accurately measure the complete SOA mass yield. Under dry conditions, substantial amounts of SOA are formed from HMML ring-opening reactions with inorganic ions and HMML organic oligomerization processes. However, the lactone organic oligomerization reactions are suppressed under more atmospherically relevant humidity levels, where hydration of the lactone is more competitive. This limits the SOA formation potential from the 2MGA pathway to HMML ring-opening reactions with water or inorganic ions under typical atmospheric conditions. The isoprene SOA mass yield from the LV pathway measured in this work is significantly higher than previous studies have reported, suggesting that low-volatility compounds such as organic nitrates and dinitrates may contribute to isoprene SOA under high-NOx conditions significantly more than previously thought and thus deserve continued study.
Understanding the efficiency and variability of photochemical ozone (O3) production from western wildfire plumes is important to accurately estimate their influence on North American air quality. A ...set of photochemical measurements were made from the NOAA Twin Otter research aircraft as a part of the Fire Influence on Regional to Global Environments and Air Quality (FIREX-AQ) experiment. We use a zero-dimensional (0-D) box model to investigate the chemistry driving O3 production in modeled plumes. Modeled afternoon plumes reached a maximum O3 mixing ratio of 140 ± 50 ppbv (average ± standard deviation) within 20 ± 10 min of emission compared to 76 ± 12 ppbv in 60 ± 30 min in evening plumes. Afternoon and evening maximum O3 isopleths indicate that plumes were near their peak in NO x efficiency. A radical budget describes the NO x volatile - organic compound (VOC) sensitivities of these plumes. Afternoon plumes displayed a rapid transition from VOC-sensitive to NO x -sensitive chemistry, driven by HO x (=OH + HO2) production from photolysis of nitrous acid (HONO) (48 ± 20% of primary HO x ) and formaldehyde (HCHO) (26 ± 9%) emitted directly from the fire. Evening plumes exhibit a slower transition from peak NO x efficiency to VOC-sensitive O3 production caused by a reduction in photolysis rates and fire emissions. HO x production in evening plumes is controlled by HONO photolysis (53 ± 7%), HCHO photolysis (18 ± 9%), and alkene ozonolysis (17 ± 9%).
Methacryloyl peroxynitrate (MPAN), the acyl peroxynitrate of methacrolein, has been suggested to be an important secondary organic aerosol (SOA) precursor from isoprene oxidation. Yet, the mechanism ...by which MPAN produces SOA through reaction with the hydroxyl radical (OH) is unclear. We systematically evaluate three proposed mechanisms in controlled chamber experiments and provide the first experimental support for the theoretically-predicted lactone formation pathway from the MPAN + OH reaction, producing hydroxymethyl-methyl-α-lactone (HMML). The decomposition of the MPAN-OH adduct yields HMML + NO3 (∼75%) and hydroxyacetone + CO + NO3 (∼25%), out-competing its reaction with atmospheric oxygen. The production of other proposed SOA precursors, e.g., methacrylic acid epoxide (MAE), from MPAN and methacrolein are negligible (<2%). Furthermore, we show that the beta-alkenyl moiety of MPAN is critical for lactone formation. Alkyl radicals formed cold via H-abstraction by OH do not decompose to HMML, even if they are structurally identical to the MPAN-OH adduct. The SOA formation from HMML, from polyaddition of the lactone to organic compounds at the particle interface or in the condensed phase, is close to unity under dry conditions. However, the SOA yield is sensitive to particle liquid water and solvated ions. In hydrated inorganic particles, HMML reacts primarily with H2O to produce the monomeric 2-methylglyceric acid (2MGA) or with aqueous sulfate and nitrate to produce the associated organosulfate and organonitrate, respectively. 2MGA, a tracer for isoprene SOA, is semivolatile and its accommodation in aerosol water decreases with decreasing pH. Conditions that enhance the production of neutral 2MGA suppress SOA mass from the HMML channel. Considering the liquid water content and pH ranges of ambient particles, 2MGA will exist largely as a gaseous compound in some parts of the atmosphere.
Decades of air quality improvements have substantially reduced the motor vehicle emissions of volatile organic compounds (VOCs). Today, volatile chemical products (VCPs) are responsible for half of ...the petrochemical VOCs emitted in major urban areas. We show that VCP emissions are ubiquitous in US and European cities and scale with population density. We report significant VCP emissions for New York City (NYC), including a monoterpene flux of 14.7 to 24.4 kg ⋅ d−1 ⋅ km−2 from fragranced VCPs and other anthropogenic sources, which is comparable to that of a summertime forest. Photochemical modeling of an extreme heat event, with ozone well in excess of US standards, illustrates the significant impact of VCPs on air quality. In the most populated regions of NYC, ozone was sensitive to anthropogenic VOCs (AVOCs), even in the presence of biogenic sources. Within this VOC-sensitive regime, AVOCs contributed upwards of ∼20 ppb to maximum 8-h average ozone. VCPs accounted for more than 50% of this total AVOC contribution. Emissions from fragranced VCPs, including personal care and cleaning products, account for at least 50% of the ozone attributed to VCPs. We show that model simulations of ozone depend foremost on the magnitude of VCP emissions and that the addition of oxygenated VCP chemistry impacts simulations of key atmospheric oxidation products. NYC is a case study for developed megacities, and the impacts of VCPs on local ozone are likely similar for other major urban regions across North America or Europe.
The extent to which emission control technologies and policies have reduced anthropogenic NO x emissions from motor vehicles is large but uncertain. We evaluate a fuel-based emission inventory for ...southern California during the June 2021 period, coinciding with the Re-Evaluating the Chemistry of Air Pollutants in CAlifornia (RECAP-CA) field campaign. A modified version of the Fuel-based Inventory of Vehicle Emissions (FIVE) is presented, incorporating 1.3 km resolution gridding and a new light-/medium-duty diesel vehicle category. NO x concentrations and weekday–weekend differences were predicted using the WRF-Chem model and evaluated using satellite and aircraft observations. Model performance was similar on weekdays and weekends, indicating appropriate day-of-week scaling of NO x emissions and a reasonable distribution of emissions by sector. Large observed weekend decreases in NO x are mainly due to changes in on-road vehicle emissions. The inventory presented in this study suggests that on-road vehicles were responsible for 55–72% of the NO x emissions in the South Coast Air Basin, compared to the corresponding fraction (43%) in the planning inventory from the South Coast Air Quality Management District. This fuel-based inventory suggests on-road NO x emissions that are 1.5 ± 0.4, 2.8 ± 0.6, and 1.3 ± 0.7 times the reference EMFAC model estimates for on-road gasoline, light- and medium-duty diesel, and heavy-duty diesel, respectively.
Chemical mechanisms describe the atmospheric
transformations of organic and inorganic species and connect air emissions
to secondary species such as ozone, fine particles, and hazardous air
...pollutants (HAPs) like formaldehyde. Recent advances in our understanding of
several chemical systems and shifts in the drivers of atmospheric chemistry
warrant updates to mechanisms used in chemical transport models such as the
Community Multiscale Air Quality (CMAQ) modeling system. This work builds on
the Regional Atmospheric Chemistry Mechanism version 2 (RACM2) and develops
the Community Regional Atmospheric Chemistry Multiphase Mechanism (CRACMM)
version 1.0, which demonstrates a fully coupled representation of chemistry
leading to ozone and secondary organic aerosol (SOA) with consideration of
HAPs. CRACMMv1.0 includes 178 gas-phase species, 51 particulate species,
and 508 reactions spanning gas-phase and heterogeneous pathways. To support
estimation of health risks associated with HAPs, nine species in CRACMM
cover 50 % of the total cancer and 60 % of the total non-cancer
emission-weighted toxicity estimated for primary HAPs from anthropogenic and
biomass burning sources in the US, with the coverage of toxicity higher
(> 80 %) when secondary formaldehyde and acrolein are
considered. In addition, new mechanism species were added based on the
importance of their emissions for the ozone, organic aerosol, or atmospheric
burden of total reactive organic carbon (ROC): sesquiterpenes, furans,
propylene glycol, alkane-like low- to intermediate-volatility organic
compounds (9 species), low- to intermediate-volatility oxygenated species (16 species), intermediate-volatility aromatic hydrocarbons (2 species), and
slowly reacting organic carbon. Intermediate- and lower-volatility organic
compounds were estimated to increase the coverage of anthropogenic and
biomass burning ROC emissions by 40 % compared to current operational
mechanisms. Autoxidation, a gas-phase reaction particularly effective in
producing SOA, was added for C10 and larger alkanes, aromatic
hydrocarbons, sesquiterpenes, and monoterpene systems including second-generation aldehydes. Integrating the radical and SOA chemistry put
additional constraints on both systems and enabled the implementation of
previously unconsidered SOA pathways from phenolic and furanone compounds,
which were predicted to account for ∼ 30 % of total aromatic
hydrocarbon SOA under typical atmospheric conditions. CRACMM organic aerosol
species were found to span the atmospherically relevant range of species carbon
number, number of oxygens per carbon, and oxidation state with a slight high
bias in the number of hydrogens per carbon. In total, 11 new emitted species
were implemented as precursors to SOA compared to current CMAQv5.3.3
representations, resulting in a bottom-up prediction of SOA, which is
required for accurate source attribution and the design of control strategies.
CRACMMv1.0 is available in CMAQv5.4.