Search

Earth system and environmental impact studies need high quality and up-to-date estimates of atmospheric deposition. This study demonstrates the methodological benefits of multimodel ensemble and measurement-model fusion mapping approaches for atmospheric deposition focusing on 2010, a year for which several studies were conducted. Global model-only deposition assessment can be further improved by integrating new model-measurement techniques, including expanded capabilities of satellite observations of atmospheric composition. We identify research and implementation priorities for timely estimates of deposition globally as implemented by the World Meteorological Organization.

We describe the baseline model configuration and simulation characteristics of the Geophysical Fluid Dynamics Laboratory (GFDL)'s Atmosphere Model version 4.1 (AM4.1), which builds on developments at GFDL over 2013–2018 for coupled carbon‐chemistry‐climate simulation as part of the sixth phase of the Coupled Model Intercomparison Project. In contrast with GFDL's AM4.0 development effort, which focused on physical and aerosol interactions and which is used as the atmospheric component of CM4.0, AM4.1 focuses on comprehensiveness of Earth system interactions. Key features of this model include doubled horizontal resolution of the atmosphere (~200 to ~100 km) with revised dynamics and physics from GFDL's previous‐generation AM3 atmospheric chemistry‐climate model. AM4.1 features improved representation of atmospheric chemical composition, including aerosol and aerosol precursor emissions, key land‐atmosphere interactions, comprehensive land‐atmosphere‐ocean cycling of dust and iron, and interactive ocean‐atmosphere cycling of reactive nitrogen. AM4.1 provides vast improvements in fidelity over AM3, captures most of AM4.0's baseline simulations characteristics, and notably improves on AM4.0 in the representation of aerosols over the Southern Ocean, India, and China—even with its interactive chemistry representation—and in its manifestation of sudden stratospheric warmings in the coldest months. Distributions of reactive nitrogen and sulfur species, carbon monoxide, and ozone are all substantially improved over AM3. Fidelity concerns include degradation of upper atmosphere equatorial winds and of aerosols in some regions.

Identifying the contributions of chemistry and transport to observed ozone pollution using regional-to-global models relies on accurate representation of ozone dry deposition. We use a recently developed configuration of the NOAA GFDL chemistry-climate model - in which the atmosphere and land are coupled through dry deposition-to investigate the influence of ozone dry deposition on ozone pollution over northern midlatitudes. In our model, deposition pathways are tied to dynamic terrestrial processes, such as photosynthesis and water cycling through the canopy and soil. Small increases in winter deposition due to more process-based representation of snow and deposition to surfaces reduce hemispheric-scale ozone throughout the lower troposphere by 5-12 ppb, improving agreement with observations relative to a simulation with the standard configuration for ozone dry deposition. Declining snow cover by the end of the 21st-century tempers the previously identified influence of rising methane on winter ozone. Dynamic dry deposition changes summer surface ozone by -4 to +7 ppb. While previous studies emphasize the importance of uptake by plant stomata, new diagnostic tracking of depositional pathways reveals a widespread impact of nonstomatal deposition on ozone pollution. Daily variability in both stomatal and nonstomatal deposition contribute to daily variability in ozone pollution. Twenty-first century changes in summer deposition result from a balance among changes in individual pathways, reflecting differing responses to both high carbon dioxide (through plant physiology versus biomass accumulation) and water availability. Our findings highlight a need for constraints on the processes driving ozone dry deposition to test representation in regional-to-global models.

This study highlights a previously under-appreciated “climate penalty” feedback mechanism - namely, substantial reductions of ozone uptake by water stressed vegetation – as a missing piece to the puzzle of why European ozone pollution episodes have not decreased as expected in recent decades, despite marked reductions in regional emissions of ozone precursors due to regulatory changes. The most extreme ozone pollution episodes are linked to heatwaves and droughts, which are increasing in frequency and intensity over Europe, with severe impacts on natural and human systems. Under drought stress, plants close their stomata to reduce water loss, consequently limiting the ozone uptake by vegetation (a component of dry deposition), leading to increased surface ozone concentrations. Such land-biosphere feedbacks are often overlooked in prior air quality projections, owing to a lack of process-based model formulations. Here, we use six decades of observations and Earth system model simulations (1960-2018) with an interactive dry deposition scheme to show that declining ozone removal by water-stressed vegetation in the warming climate exacerbate ozone air pollution over Europe. Incorporated into a dynamic vegetation land – atmospheric chemistry – climate model, the dry deposition scheme mechanistically describes the response of ozone deposition to atmospheric CO2 concentration, canopy air vapor pressure deficit, and soil water availability. Our observational and modeling analyses reveal drought stress causing as much as 70% reductions in ozone removal by forests. Reduced ozone removal by water-stressed vegetation worsens peak ozone episodes during European mega-droughts, such as the 2003 event, offsetting much of the air quality improvements gained from regional emission controls. Accounting for vegetation feedbacks leads to a three-fold increase in high surface ozone events above 80 ppbv (8-hour average) and a 20% increase in the sensitivity of ozone pollution extremes (95th percentile) to increasing temperature. As the frequency of hot and dry summers is expected to increase in the coming decades, this ozone climate penalty could be severe and therefore needs to be considered when designing clean air policy in the European Union. Notes: This study is currently under review for possible publication in Nature Climate Change.

The top-of-atmosphere (TOA) outgoing longwave flux over the 9.6 µm ozone band is a fundamental quantity for understanding chemistry–climate coupling. However, observed TOA fluxes are hard to estimate as they exhibit considerable variability in space and time that depend on the distributions of clouds, ozone (O3), water vapor (H2O), air temperature (Ta), and surface temperature (Ts). Benchmarking present-day fluxes and quantifying the relative influence of their drivers is the first step for estimating climate feedbacks from ozone radiative forcing and predicting radiative forcing evolution. To that end, we constructed observational instantaneous radiative kernels (IRKs) under clear-sky conditions, representing the sensitivities of the TOA flux in the 9.6 µm ozone band to the vertical distribution of geophysical variables, including O3, H2O, Ta, and Ts based upon the Aura Tropospheric Emission Spectrometer (TES) measurements. Applying these kernels to present-day simulations from the Chemistry-Climate Model Initiative (CCMI) project as compared to a 2006 reanalysis assimilating satellite observations, we show that the models have large differences in TOA flux, attributable to different geophysical variables. In particular, model simulations continue to diverge from observations in the tropics, as reported in previous studies of the Atmospheric Chemistry Climate Model Intercomparison Project (ACCMIP) simulations. The principal culprits are tropical middle and upper tropospheric ozone followed by tropical lower tropospheric H2O. Five models out of the eight studied here have TOA flux biases exceeding 100 mW m−2 attributable to tropospheric ozone bias. Another set of five models have flux biases over 50 mW m−2 due to H2O. On the other hand, Ta radiative bias is negligible in all models (no more than 30 mW m−2). We found that the atmospheric component (AM3) of the Geophysical Fluid Dynamics Laboratory (GFDL) general circulation model and Canadian Middle Atmosphere Model (CMAM) have the lowest TOA flux biases globally but are a result of cancellation of opposite biases due to different processes. Overall, the multi-model ensemble mean bias is −133±98 mW m−2, indicating that they are too atmospherically opaque due to trapping too much radiation in the atmosphere by overestimated tropical tropospheric O3 and H2O. Having too much O3 and H2O in the troposphere would have different impacts on the sensitivity of TOA flux to O3 and these competing effects add more uncertainties on the ozone radiative forcing. We find that the inter-model TOA outgoing longwave radiation (OLR) difference is well anti-correlated with their ozone band flux bias. This suggests that there is significant radiative compensation in the calculation of model outgoing longwave radiation., Atmospheric Chemistry and Physics, 20 (1), ISSN:1680-7375, ISSN:1680-7367

Widespread efforts to abate ozone (O3) smog have significantly reduced nitrogen oxides (NOx) emissions over the past two decades in the Southeast U.S. (SEUS), a place heavily influenced by both anthropogenic and biogenic emissions. How reactive nitrogen speciation responds to the reduction in NOx emissions in this region remains to be elucidated. Here we exploit aircraft measurements from ICARTT (July–August, 2004), SENEX (June–July, 2013), and SEAC4RS (August–September, 2013) and long-term ground measurement networks alongside a global chemistry-climate model to examine decadal changes in summertime reactive nitrogen species and ozone over the Southeast U.S. We find that most reactive nitrogen species, including NOx, peroxyacetyl nitrate (PAN) and nitric acid (HNO3) decline proportionally with decreasing NOx emissions in this region, leading to a similar decline in exported NOy. This linear response is in part due to the nearly constant summertime supply of biogenic VOC emissions in this region. Our model captures the observed relative change of reactive nitrogen species and surface ozone from 2004 to 2013. Model sensitivity tests indicate that further reductions of NOx emissions will lead to a continued decline in surface ozone and less frequent extreme ozone events.

Peroxyacetyl nitrate (PAN) formed in the atmospheric oxidation of non-methane volatile organic compounds (NMVOCs) is the principal tropospheric reservoir for nitrogen oxide radicals (NOx = NO + NO2). PAN enables the transport and release of NOx to the remote troposphere with major implications for the global distributions of ozone and OH, the main tropospheric oxidants. Simulation of PAN is a challenge for global models because of the dependence of PAN on vertical transport as well as complex and uncertain NMVOC sources and chemistry. Here we use an improved representation of NMVOCs in a global 3-D chemical transport model (GEOS-Chem) and show that it can simulate PAN observations from aircraft campaigns worldwide. The immediate carbonyl precursors for PAN formation include acetaldehyde (44% of the global source), methylglyoxal (30%), acetone (7%), and a suite of other isoprene and terpene oxidation products (19%). A diversity of NMVOC emissions is responsible for PAN formation globally including isoprene (37%) and alkanes (14%). Anthropogenic sources are dominant in the extratropical Northern Hemisphere outside the growing season. Open fires appear to play little role except at high northern latitudes in spring, although results are very sensitive to plume chemistry and plume rise. Lightning NOx is the dominant contributor to the observed PAN maximum in the free troposphere over the South Atlantic.

Alkyl and multifunctional organic nitrates, molecules of the chemical form RONO2, are products of chain terminating reactions in the tropospheric HOx and NOx catalytic cycles and thereby impact ozone formation locally. Many of the molecules in the class have lifetimes that are long enough that they can be transported over large distances. If the RONO2 then decompose to deliver NOx to remote regions they affect ozone production rates in locations distant from the original NOx source. While measurements of total RONO2 (ΣANs) and small straight chain alkyl nitrates are routine, measurements of the specific multifunctional RONO2 molecules that are believed to dominate the total have rarely been reported and never reported in coincidence with ambient ΣANs measurements. Here we describe observations obtained during the BEARPEX 2009 experiment including ΣANs and a suite of multifunctional nitrates including isoprene derived hydroxynitrates, oxidation products of those nitrates, 2-methyl-3-buten-2-ol (MBO) derived hydroxynitrates, and monoterpene nitrates. At the BEARPEX field site, the sum of the individual biogenically derived nitrates account for two-thirds of the ΣANs, confirming predictions of the importance of biogenic nitrates to the NOy budget. Isoprene derived nitrates, transported to the site, are a much larger fraction of the ΣANs at the site than the nitrates derived from the locally emitted MBO. Evidence for additional nitrates, possibly from nocturnal chemistry of isoprene and α-pinene, is presented.

Tropical tropospheric ozone affects Earth's radiative forcing and the oxidative capacity of the atmosphere. Considerable work has been devoted to the study of the processes controlling its budget. Yet, large discrepancies between simulated and observed tropical tropospheric ozone remain. Here, we characterize some of the mechanisms by which the photochemistry of isoprene impacts the budget of tropical ozone. At the regional scale, we use forward sensitivity simulation to explore the sensitivity to the representation of isoprene nitrates. We find that isoprene nitrates can account for up to 70% of the local NOx = NO+NO2 sink. The resulting modulation of ozone can be well characterized by their net modulation of NOx. We use adjoint sensitivity simulations to demonstrate that the oxidation of isoprene can affect ozone outside of continental regions through the transport of NOx over near-shore regions (e.g., South Atlantic) and the oxidation of isoprene outside of the boundary layer far from its emissions regions. The latter mechanism is promoted by the simulated low boundary-layer oxidative conditions. In our simulation, ~20% of the isoprene is oxidized above the boundary layer in the tropics. Changes in the interplay between regional and global effect are discussed in light of the forecasted increase in anthropogenic emissions in tropical regions.

We report fluxes and dry deposition velocities for 16 atmospheric compounds above a southeastern United States forest, including: hydrogen peroxide (H_2O_2), nitric acid (HNO_3), hydrogen cyanide (HCN), hydroxymethyl hydroperoxide, peroxyacetic acid, organic hydroxy nitrates, and other multifunctional species derived from the oxidation of isoprene and monoterpenes. The data suggest that dry deposition is the dominant daytime sink for small, saturated oxygenates. Greater than 6 wt %C emitted as isoprene by the forest was returned by dry deposition of its oxidized products. Peroxides account for a large fraction of the oxidant flux, possibly eclipsing ozone in more pristine regions. The measured organic nitrates comprise a sizable portion (15%) of the oxidized nitrogen input into the canopy, with HNO_3 making up the balance. We observe that water-soluble compounds (e.g., strong acids and hydroperoxides) deposit with low surface resistance whereas compounds with moderate solubility (e.g., organic nitrates and hydroxycarbonyls) or poor solubility (e.g., HCN) exhibited reduced uptake at the surface of plants. To first order, the relative deposition velocities of water-soluble compounds are constrained by their molecular diffusivity. From resistance modeling, we infer a substantial emission flux of formic acid at the canopy level (∼1 nmol m^(−2)⋅s^(−1)). GEOS−Chem, a widely used atmospheric chemical transport model, currently underestimates dry deposition for most molecules studied in this work. Reconciling GEOS−Chem deposition velocities with observations resulted in up to a 45% decrease in the simulated surface concentration of trace gases.

The CMAQ (Community Multiscale Air Quality) us model in combination with observations for INTEX-NA/ICARTT (Intercontinental Chemical Transport Experiment–North America/International Consortium for Atmospheric Research on Transport and Transformation) 2004 are used to evaluate recent advances in isoprene oxidation chemistry and provide constraints on isoprene nitrate yields, isoprene nitrate lifetimes, and NOx recycling rates. We incorporate recent advances in isoprene oxidation chemistry into the SAPRC-07 chemical mechanism within the US EPA (United States Environmental Protection Agency) CMAQ model. The results show improved model performance for a range of species compared against aircraft observations from the INTEX-NA/ICARTT 2004 field campaign. We further investigate the key processes in isoprene nitrate chemistry and evaluate the impact of uncertainties in the isoprene nitrate yield, NOx (NOx = NO + NO2) recycling efficiency, dry deposition velocity, and RO2 + HO2 reaction rates. We focus our examination on the southeastern United States, which is impacted by both abundant isoprene emissions and high levels of anthropogenic pollutants. We find that NOx concentrations increase by 4–9% as a result of reduced removal by isoprene nitrate chemistry. O3 increases by 2 ppbv as a result of changes in NOx. OH concentrations increase by 30%, which can be primarily attributed to greater HOx production. We find that the model can capture observed total alkyl and multifunctional nitrates (∑ANs) and their relationship with O3 by assuming either an isoprene nitrate yield of 6% and daytime lifetime of 6 hours or a yield of 12% and lifetime of 4 h. Uncertainties in the isoprene nitrates can impact ozone production by 10% and OH concentrations by 6%. The uncertainties in NOx recycling efficiency appear to have larger effects than uncertainties in isoprene nitrate yield and dry deposition velocity. Further progress depends on improved understanding of isoprene oxidation pathways, the rate of NOx recycling from isoprene nitrates, and the fate of the secondary, tertiary, and further oxidation products of isoprene.

The C(5)-hydroperoxyenals (C(5)-HPALDs) are a newly-recognized class of multi-functional hydrocarbons produced during the hydroxyl radical (OH)-initiated oxidation of isoprene. Recent theoretical calculations suggest that fast photolysis of these compounds may be an important OH source in high-isoprene, low-NO regions. We report experimental constraints for key parameters of photolysis, OH reaction and ozone reaction of these compounds as derived from a closely-related, custom-synthesized C(6)-HPALD. The photolysis quantum yield is 1.0 ± 0.4 over the range 300-400 nm, assuming an absorption cross section equal to the average of those measured for several analogous enals. The yield of OH from photolysis was determined as 1.0 ± 0.8. The OH reaction rate constant is (5.1 ± 1.8) × 10(-11) cm(3) molecule(-1) s(-1) at 296 K. The ozone reaction rate constant is (1.2 ± 0.2) × 10(-18) cm(3) molecule(-1) s(-1) at 296 K. These results are consistent with previous first-principles estimates, though the nature and fate of secondary oxidation products remains uncertain. Incorporation of C(5)-HPALD chemistry with the above parameters in a 0-D box model, along with experimentally-constrained rates for C(5)-HPALD production from isomerization of first-generation isoprene hydroxyperoxy radicals, is found to enhance modeled OH concentrations by 5-16% relative to the traditional isoprene oxidation mechanism for the chemical regimes of recent observational studies in rural and remote regions. This enhancement in OH will increase if C(5)-HPALD photo-oxidation products also photolyze to yield additional OH or if the C(5)-HPALD production rate is faster than has been observed.

We implement a new isoprene oxidation mechanism in a global 3-D chemical transport model (GEOS-Chem). Model results are evaluated with observations for ozone, isoprene oxidation products, and related species from the International Consortium for Atmospheric Research on Transport and Transformation aircraft campaign over the eastern United States in summer 2004. The model achieves an unbiased simulation of ozone in the boundary layer and the free troposphere, reflecting canceling effects from recent model updates for isoprene chemistry, bromine chemistry, and HO2 loss to aerosols. Simulation of the ozone-CO correlation is improved relative to previous versions of the model, and this is attributed to a lower and reversible yield of isoprene nitrates, increasing the ozone production efficiency per unit of nitrogen oxides (NO_x ≡ NO + NO_2). The model successfully reproduces the observed concentrations of organic nitrates (∑ANs) and their correlations with HCHO and ozone. ∑ANs in the model is principally composed of secondary isoprene nitrates, including a major contribution from nighttime isoprene oxidation. The correlations of ∑ANs with HCHO and ozone then provide sensitive tests of isoprene chemistry and argue in particular against a fast isomerization channel for isoprene peroxy radicals. ∑ANs can provide an important reservoir for exporting NO_x from the U.S. boundary layer. We find that the dependence of surface ozone on isoprene emission is positive throughout the U.S., even if NO_x emissions are reduced by a factor of 4. Previous models showed negative dependences that we attribute to erroneous titration of OH by isoprene.

The response of ozone (O3) dry deposition to ecosystem‐atmosphere interactions is poorly understood but is central to determining the potential for extreme pollution events under current and future climate conditions. Using observations and an interactive dry deposition scheme within two dynamic vegetation land models (GFDL LM3.0/LM4.0) driven by observation‐based meteorological forcings over 1948‐2014, we investigate the factors controlling seasonal and interannual variability (IAV) in O3 deposition velocities (Vd,O3). Stomatal activity in this scheme is determined mechanistically, depending on phenology, soil moisture, vapor pressure deficit, and CO2 concentration. Soil moisture plays a key role in modulating the observed and simulated Vd,O3 seasonal changes over evergreen forests in Mediterranean Europe, South Asia, and the Amazon. Analysis of multi‐year observations at forest sites in Europe and North America reveals drought stress to reduce Vd,O3 by ~50%. Both LM3.0 and LM4.0 capture the observed Vd,O3 decreases due to drought; however, IAV is weaker by a factor of two in LM3.0 coupled to atmospheric models, particularly in regions with large precipitation biases. IAV in summertime Vd,O3 to forests, driven primarily by the stomatal pathway, is largest (15‐35%) in semi‐arid regions of western Europe, eastern North America, and northeastern China. Monthly mean Vd,O3 for the highest year is two to four times that of the lowest, with significant implications for surface O3 variability and extreme events. Using Vd,O3 from LM4.0 in an atmospheric chemistry model improves the simulation of surface O3 abundance and spatial variability (reduces mean biases by ~10 ppb) relative to the widely‐used Wesely scheme.