Impacts of updated reaction kinetics on the global GEOS-Chem simulation of atmospheric chemistry

We updated the chemical mechanism of the GEOS-Chem global 3-D model of atmospheric chemistry to include new recommendations from the NASA Jet Propulsion Laboratory (JPL) chemical kinetics Data Evaluation 19-5 and from the International Union of Pure and Applied Chemistry (IUPAC) and to balance carbon and nitrogen. We examined the impact of these updates on the GEOS-Chem version 14.0.1 simulation. Notable changes include 11 updates to reactions of reactive nitrogen species, resulting in a 7% net increase in the stratospheric NOx (NO + NO2) burden; an updated CO + OH rate formula leading to a 2.7% reduction in total tropospheric CO; adjustments to the rate coefficient and branching ratios of propane + OH, leading to reduced tropospheric propane (-17%) and increased acetone (+3.5%) burdens; a 41% increase in the tropospheric burden of peroxyacetic acid due to a decrease in the rate coefficient for its reaction with OH, further contributing to reductions in peroxyacetyl nitrate (PAN; -3.8%) and acetic acid (-3.4%); and a number of minor adjustments to halogen radical cycling. Changes to the global tropospheric burdens of other species include -0.7% for ozone, +0.3% for OH (-0.4% for methane lifetime against oxidation by tropospheric OH), +0.8% for formaldehyde, and -1.7% for NOx. The updated mechanism reflects the current state of the science, including complex chemical dependencies of key atmospheric species on temperature, pressure, and concentrations of other compounds. The improved conservation of carbon and nitrogen will facilitate future studies of their overall atmospheric budgets.

Atmospheric Gas-Phase Formation of Methanesulfonic Acid

Despite its impact on the climate, the mechanism of methanesulfonic acid (MSA) formation in the oxidation of dimethyl sulfide (DMS) remains unclear. The DMS + OH reaction is known to form methanesulfinic acid (MSIA), methane sulfenic acid (MSEA), the methylthio radical (CH3S), and hydroperoxymethyl thioformate (HPMTF). Among them, HPMTF reacts further to form SO2 and OCS, while the other three form the CH3SO2 radical. Based on theoretical calculations, we find that the CH3SO2 radical can add O2 to form CH3S(O)2OO, which can react further to form MSA. The branching ratio is highly temperature sensitive, and the MSA yield increases with decreasing temperature. In warmer regions, SO2 is the dominant product of DMS oxidation, while in colder regions, large amounts of MSA can form. Global modeling indicates that the proposed temperature-sensitive MSA formation mechanism leads to a substantial increase in the simulated global atmospheric MSA formation and burden.

Overview of ICARUS-A Curated, Open Access, Online Repository for Atmospheric Simulation Chamber Data

Atmospheric simulation chambers continue to be indispensable tools for research in the atmospheric sciences. Insights from chamber studies are integrated into atmospheric chemical transport models, which are used for science-informed policy decisions. However, a centralized data management and access infrastructure for their scientific products had not been available in the United States and many parts of the world. ICARUS (Integrated Chamber Atmospheric data Repository for Unified Science) is an open access, searchable, web-based infrastructure for storing, sharing, discovering, and utilizing atmospheric chamber data [https://icarus.ucdavis.edu]. ICARUS has two parts: a data intake portal and a search and discovery portal. Data in ICARUS are curated, uniform, interactive, indexed on popular search engines, mirrored by other repositories, version-tracked, vocabulary-controlled, and citable. ICARUS hosts both legacy data and new data in compliance with open access data mandates. Targeted data discovery is available based on key experimental parameters, including organic reactants and mixtures that are managed using the PubChem chemical database, oxidant information, nitrogen oxide (NOx) content, alkylperoxy radical (RO₂) fate, seed particle information, environmental conditions, and reaction categories. A discipline-specific repository such as ICARUS with high amounts of metadata works to support the evaluation and revision of atmospheric model mechanisms, intercomparison of data and models, and the development of new model frameworks that can have more predictive power in the current and future atmosphere. The open accessibility and interactive nature of ICARUS data may also be useful for teaching, data mining, and training machine learning models.

An Inversion Framework for Optimizing Non-Methane VOC Emissions Using Remote Sensing and Airborne Observations in Northeast Asia During the KORUS-AQ Field Campaign

We aim to reduce uncertainties in CH₂O and other volatile organic carbon (VOC) emissions through assimilation of remote sensing data. We first update a three-dimensional (3D) chemical transport model, GEOS-Chem with the KORUSv5 anthropogenic emission inventory and inclusion of chemistry for aromatics and C₂H₄, leading to modest improvements in simulation of CH₂O (normalized mean bias (NMB): -0.57 to -0.51) and O₃ (NMB: -0.25 to -0.19) compared against DC-8 aircraft measurements during KORUS-AQ; the mixing ratio of most VOC species are still underestimated. We next constrain VOC emissions using CH₂O observations from two satellites (OMI and OMPS) and the DC-8 aircraft during KORUS-AQ. To utilize data from multiple platforms in a consistent manner, we develop a two-step Hybrid Iterative Finite Difference Mass Balance and four-dimensional variational inversion system (Hybrid IFDMB-4DVar). The total VOC emissions throughout the domain increase by 47%. The a posteriori simulation reduces the low biases of simulated CH₂O (NMB: -0.51 to -0.15), O₃ (NMB: -0.19 to -0.06), and VOCs. Alterations to the VOC speciation from the 4D-Var inversion include increases of biogenic isoprene emissions in Korea and anthropogenic emissions in Eastern China. We find that the IFDMB method alone is adequate for reducing the low biases of VOCs in general; however, 4D-Var provides additional refinement of high-resolution emissions and their speciation. Defining reasonable emission errors and choosing optimal regularization parameters are crucial parts of the inversion system. Our new hybrid inversion framework can be applied for future air quality campaigns, maximizing the value of integrating measurements from current and upcoming geostationary satellite instruments.

Gas-Phase Oxidation Rates and Products of 1,2-Dihydroxy Isoprene

1,2-Dihydroxy isoprene (1,2-DHI), a product of isoprene oxidation from multiple chemical pathways, is produced in the atmosphere in large quantities; however, its chemical fate has not been comprehensively studied. Here, we perform chamber experiments to investigate its gas-phase reactions. We find that the reactions of 1,2-DHI with OH radicals and ozone are rapid (k_OH = 8.0 (±1.3) × 10^-11 cm³ molecule^-1 s^-1; k_O3 = 7.2 (±1.1) × 10^-18 cm³ molecule^-1 s^-1). Reaction with OH, which dominates 1,2-DHI loss, leads primarily to fragmentation and radical recycling; major products under both high- and low-NO conditions include hydroxyacetone, glycolaldehyde, and 2,3-dihydroxy-2-methyl-propanal (DHMP). Radical-terminating hydroperoxide formation from the peroxy radical (RO₂) reaction with HO₂ and organonitrate formation from RO₂ + NO are not observed in the gas phase, possibly due to low volatility; constraints for their branching ratios are instead derived by mass balance. We also measure secondary organic aerosol mass yields from 1,2-DHI (0-23%) and show that oxidation in the presence of aqueous particles leads to formic and acetic acid production. Finally, we incorporate results into GEOS-Chem, a global chemical transport model, to compute the global production (25.3 Tg a^-1) and gas-phase loss (20.2 Tg a^-1) of 1,2-DHI and show that its oxidation provides non-negligible contributions to the atmospheric budgets of hydroxyacetone, glycolaldehyde, hydroxymethyl hydroperoxide, formic acid, and DHMP.

Ozone pollution in the North China Plain spreading into the late-winter haze season

Surface ozone is a severe air pollution problem in the North China Plain, which is home to 300 million people. Ozone concentrations are highest in summer, driven by fast photochemical production of hydrogen oxide radicals (HOx) that can overcome the radical titration caused by high emissions of nitrogen oxides (NOx) from fuel combustion. Ozone has been very low during winter haze (particulate) pollution episodes. However, the abrupt decrease of NOx emissions following the COVID-19 lockdown in January 2020 reveals a switch to fast ozone production during winter haze episodes with maximum daily 8-h average (MDA8) ozone concentrations of 60 to 70 parts per billion. We reproduce this switch with the GEOS-Chem model, where the fast production of ozone is driven by HOx radicals from photolysis of formaldehyde, overcoming radical titration from the decreased NOx emissions. Formaldehyde is produced by oxidation of reactive volatile organic compounds (VOCs), which have very high emissions in the North China Plain. This remarkable switch to an ozone-producing regime in January-February following the lockdown illustrates a more general tendency from 2013 to 2019 of increasing winter-spring ozone in the North China Plain and increasing association of high ozone with winter haze events, as pollution control efforts have targeted NOx emissions (30% decrease) while VOC emissions have remained constant. Decreasing VOC emissions would avoid further spreading of severe ozone pollution events into the winter-spring season.

Satellite isoprene retrievals constrain emissions and atmospheric oxidation

Isoprene is the dominant non-methane organic compound emitted to the atmosphere1-3. It drives ozone and aerosol production, modulates atmospheric oxidation and interacts with the global nitrogen cycle4-8. Isoprene emissions are highly uncertain1,9, as is the nonlinear chemistry coupling isoprene and the hydroxyl radical, OH-its primary sink10-13. Here we present global isoprene measurements taken from space using the Cross-track Infrared Sounder. Together with observations of formaldehyde, an isoprene oxidation product, these measurements provide constraints on isoprene emissions and atmospheric oxidation. We find that the isoprene-formaldehyde relationships measured from space are broadly consistent with the current understanding of isoprene-OH chemistry, with no indication of missing OH recycling at low nitrogen oxide concentrations. We analyse these datasets over four global isoprene hotspots in relation to model predictions, and present a quantification of isoprene emissions based directly on satellite measurements of isoprene itself. A major discrepancy emerges over Amazonia, where current underestimates of natural nitrogen oxide emissions bias modelled OH and hence isoprene. Over southern Africa, we find that a prominent isoprene hotspot is missing from bottom-up predictions. A multi-year analysis sheds light on interannual isoprene variability, and suggests the influence of the El Niño/Southern Oscillation.

Sulfate Formation via Cloud Processing from Isoprene Hydroxyl Hydroperoxides (ISOPOOH)

The oxidation of sulfur dioxide (SO₂) by peroxides leads to the formation of sulfate in cloudwater, contributing to particulate matter (PM) formation. The reaction with hydrogen peroxide (H₂O₂) is considered to be the main cloud oxidation pathway. Previous studies have examined the oxidation of SO₂ in cloudwater by small organic peroxides with one functional group; however, oxidation by multifunctional organic hydroperoxides, which are expected to have higher water solubility and reactivity, has not been examined. We investigate the aqueous oxidation of SO₂ by the two main isomers of isoprene hydroxyl hydroperoxide (ISOPOOH), the primary low-NO_x isoprene oxidation products in the atmosphere. Having large Henry's law constants and being among the most abundant multifunctional hydroperoxides, they are among the most important organic hydroperoxides present in clouds. The pH dependence of the reactions was investigated at cloud relevant pH of 3-6, and the results reveal their importance compared to the oxidation of SO₂ via H₂O₂. Model simulations in GEOS-Chem, updated with the chemistry described herein, highlight the importance of these pathways for sulfate formation in regions with high isoprene emissions and low-NO_x atmospheric conditions, especially if they maintain significant SO₂ emissions.

Anthropogenic drivers of 2013-2017 trends in summer surface ozone in China

Observations of surface ozone available from ∼1,000 sites across China for the past 5 years (2013-2017) show severe summertime pollution and regionally variable trends. We resolve the effect of meteorological variability on the ozone trends by using a multiple linear regression model. The residual of this regression shows increasing ozone trends of 1-3 ppbv a-1 in megacity clusters of eastern China that we attribute to changes in anthropogenic emissions. By contrast, ozone decreased in some areas of southern China. Anthropogenic NOx emissions in China are estimated to have decreased by 21% during 2013-2017, whereas volatile organic compounds (VOCs) emissions changed little. Decreasing NOx would increase ozone under the VOC-limited conditions thought to prevail in urban China while decreasing ozone under rural NOx-limited conditions. However, simulations with the Goddard Earth Observing System Chemical Transport Model (GEOS-Chem) indicate that a more important factor for ozone trends in the North China Plain is the ∼40% decrease of fine particulate matter (PM2.5) over the 2013-2017 period, slowing down the aerosol sink of hydroperoxy (HO2) radicals and thus stimulating ozone production.

Gas-Phase Reactions of Isoprene and Its Major Oxidation Products

Isoprene carries approximately half of the flux of non-methane volatile organic carbon emitted to the atmosphere by the biosphere. Accurate representation of its oxidation rate and products is essential for quantifying its influence on the abundance of the hydroxyl radical (OH), nitrogen oxide free radicals (NO _x), ozone (O₃), and, via the formation of highly oxygenated compounds, aerosol. We present a review of recent laboratory and theoretical studies of the oxidation pathways of isoprene initiated by addition of OH, O₃, the nitrate radical (NO₃), and the chlorine atom. From this review, a recommendation for a nearly complete gas-phase oxidation mechanism of isoprene and its major products is developed. The mechanism is compiled with the aims of providing an accurate representation of the flow of carbon while allowing quantification of the impact of isoprene emissions on HO _x and NO _x free radical concentrations and of the yields of products known to be involved in condensed-phase processes. Finally, a simplified (reduced) mechanism is developed for use in chemical transport models that retains the essential chemistry required to accurately simulate isoprene oxidation under conditions where it occurs in the atmosphere-above forested regions remote from large NO _x emissions.

A multi-year data set on aerosol-cloud-precipitation-meteorology interactions for marine stratocumulus clouds

Airborne measurements of meteorological, aerosol, and stratocumulus cloud properties have been harmonized from six field campaigns during July-August months between 2005 and 2016 off the California coast. A consistent set of core instruments was deployed on the Center for Interdisciplinary Remotely-Piloted Aircraft Studies Twin Otter for 113 flight days, amounting to 514 flight hours. A unique aspect of the compiled data set is detailed measurements of aerosol microphysical properties (size distribution, composition, bioaerosol detection, hygroscopicity, optical), cloud water composition, and different sampling inlets to distinguish between clear air aerosol, interstitial in-cloud aerosol, and droplet residual particles in cloud. Measurements and data analysis follow documented methods for quality assurance. The data set is suitable for studies associated with aerosol-cloud-precipitation-meteorology-radiation interactions, especially owing to sharp aerosol perturbations from ship traffic and biomass burning. The data set can be used for model initialization and synergistic application with meteorological models and remote sensing data to improve understanding of the very interactions that comprise the largest uncertainty in the effect of anthropogenic emissions on radiative forcing.

Atmospheric fates of Criegee intermediates in the ozonolysis of isoprene

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.

Kinetics and Products of the Reaction of the First-Generation Isoprene Hydroxy Hydroperoxide (ISOPOOH) with OH

The atmospheric oxidation of isoprene by the OH radical leads to the formation of several isomers of an unsaturated hydroxy hydroperoxide, ISOPOOH. Oxidation of ISOPOOH by OH produces epoxydiols, IEPOX, which have been shown to contribute mass to secondary organic aerosol (SOA). We present kinetic rate constant measurements for OH + ISOPOOH using synthetic standards of the two major isomers: (1,2)- and (4,3)-ISOPOOH. At 297 K, the total OH rate constant is 7.5 ± 1.2 × 10(-11) cm(3) molecule(-1) s(-1) for (1,2)-ISOPOOH and 1.18 ± 0.19 × 10(-10) cm(3) molecule(-1) s(-1) for (4,3)-ISOPOOH. Abstraction of the hydroperoxy hydrogen accounts for approximately 12% and 4% of the reactivity for (1,2)-ISOPOOH and (4,3)-ISOPOOH, respectively. The sum of all H-abstractions account for approximately 15% and 7% of the reactivity for (1,2)-ISOPOOH and (4,3)-ISOPOOH, respectively. The major product observed from both ISOPOOH isomers was IEPOX (cis-β and trans-β isomers), with a ∼ 2:1 preference for trans-β IEPOX and similar total yields from each ISOPOOH isomer (∼ 70-80%). An IEPOX global production rate of more than 100 Tg C each year is estimated from this chemistry using a global 3D chemical transport model, similar to earlier estimates. Finally, following addition of OH to ISOPOOH, approximately 13% of the reactivity proceeds via addition of O2 at 297 K and 745 Torr. In the presence of NO, these peroxy radicals lead to formation of small carbonyl compounds. Under HO2 dominated chemistry, no products are observed from these channels. We suggest that the major products, highly oxygenated organic peroxides, are lost to the chamber walls. In the atmosphere, formation of these compounds may contribute to organic aerosol mass.

Mechanism of the hydroxyl radical oxidation of methacryloyl peroxynitrate (MPAN) and its pathway toward secondary organic aerosol formation in the atmosphere

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.

Atmospheric fate of methyl vinyl ketone: peroxy radical reactions with NO and HO2

First generation product yields from the OH-initiated oxidation of methyl vinyl ketone (3-buten-2-one, MVK) under both low and high NO conditions are reported. In the low NO chemistry, three distinct reaction channels are identified leading to the formation of (1) OH, glycolaldehyde, and acetyl peroxy R2a , (2) a hydroperoxide R2b , and (3) an α-diketone R2c . The α-diketone likely results from HOx-neutral chemistry previously only known to occur in reactions of HO2 with halogenated peroxy radicals. Quantum chemical calculations demonstrate that all channels are kinetically accessible at 298 K. In the high NO chemistry, glycolaldehyde is produced with a yield of 74 ± 6.0%. Two alkyl nitrates are formed with a combined yield of 4.0 ± 0.6%. We revise a three-dimensional chemical transport model to assess what impact these modifications in the MVK mechanism have on simulations of atmospheric oxidative chemistry. The calculated OH mixing ratio over the Amazon increases by 6%, suggesting that the low NO chemistry makes a non-negligible contribution toward sustaining the atmospheric radical pool.

Gas phase production and loss of isoprene epoxydiols

Isoprene epoxydiols (IEPOX) form in high yields from the OH-initiated oxidation of isoprene under low-NO conditions. These compounds contribute significantly to secondary organic aerosol formation. Their gas-phase chemistry has, however, remained largely unexplored. In this study, we characterize the formation of IEPOX isomers from the oxidation of isoprene by OH. We find that cis-β- and trans-β-IEPOX are the dominant isomers produced, and that they are created in an approximate ratio of 1:2 from the low-NO oxidation of isoprene. Three isomers of IEPOX, including cis-β- and trans-β, were synthesized and oxidized by OH in environmental chambers under high- and low-NO conditions. We find that IEPOX reacts with OH at 299 K with rate coefficients of (0.84 ± 0.07) × 10(-11), (1.52 ± 0.07) × 10(-11), and (0.98 ± 0.05) × 10(-11) cm(3) molecule(-1) s(-1) for the δ1, cis-β, and trans-β isomers. Finally, yields of the first-generation products of IEPOX + OH oxidation were measured, and a new mechanism of IEPOX oxidation is proposed here to account for the observed products. The substantial yield of glyoxal and methylglyoxal from IEPOX oxidation may help explain elevated levels of those compounds observed in low-NO environments with high isoprene emissions.