Articles | Volume 14, issue 9
Geosci. Model Dev., 14, 5789–5823, 2021
Geosci. Model Dev., 14, 5789–5823, 2021

Model description paper 24 Sep 2021

Model description paper | 24 Sep 2021

GCAP 2.0: a global 3-D chemical-transport model framework for past, present, and future climate scenarios

GCAP 2.0: a global 3-D chemical-transport model framework for past, present, and future climate scenarios
Lee T. Murray1,2, Eric M. Leibensperger3, Clara Orbe4, Loretta J. Mickley5, and Melissa Sulprizio5 Lee T. Murray et al.
  • 1Dept. of Earth and Environmental Sciences, University of Rochester, Rochester, NY, USA
  • 2Dept. of Physics and Astronomy, University of Rochester, Rochester, NY, USA
  • 3Dept. of Physics and Astronomy, Ithaca College, Ithaca, NY, USA
  • 4NASA Goddard Institute for Space Studies, New York, NY, USA
  • 5School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA

Correspondence: Lee T. Murray (


This paper describes version 2.0 of the Global Change and Air Pollution (GCAP 2.0) model framework, a one-way offline coupling between version E2.1 of the NASA Goddard Institute for Space Studies (GISS) general circulation model (GCM) and the GEOS-Chem global 3-D chemical-transport model (CTM). Meteorology for driving GEOS-Chem has been archived from the E2.1 contributions to phase 6 of the Coupled Model Intercomparison Project (CMIP6) for the pre-industrial era and the recent past. In addition, meteorology is available for the near future and end of the century for seven future scenarios ranging from extreme mitigation to extreme warming. Emissions and boundary conditions have been prepared for input to GEOS-Chem that are consistent with the CMIP6 experimental design. The model meteorology, emissions, transport, and chemistry are evaluated in the recent past and found to be largely consistent with GEOS-Chem driven by the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) product and with observational constraints.

1 Introduction

How atmospheric composition and chemistry has changed in the past and will change in the future is of tremendous societal importance. Surface air pollution is the leading cause of preventable death worldwide (GBD 2019 Risk Factor Collaborators2020) and threatens global food security and ecosystem health (e.g., Tai et al.2021). Short-lived climate forcers like methane, tropospheric ozone, and aerosol particles influence global and regional climate (e.g., Fiore et al.2015). The emission of ozone-depleting substances and climate change threatens the overlying stratospheric ozone layer (WMO2018). And solar radiation modification through the purposeful injection of chemical species into the atmosphere is seriously being considered to combat anthropogenic climate change. Yet, it remains highly uncertain in its efficacy or risks (e.g., Eastham et al.2021).

To study and address these issues, scientists, regulators, and policymakers frequently use 3-D chemical-transport models (CTMs). CTMs use archived meteorology to drive the spatial and temporal evolution of trace gases and particles in the atmosphere. By not needing to resolve all the equations of motion, as in a general circulation model (GCM), CTMs can expend additional computational power to resolve more complex chemistry or perform additional simulations. Furthermore, the chain of cause and effect between meteorology and composition is much easier to establish in a CTM than a fully coupled chemistry–climate model (CCM) or Earth system model (ESM) running online chemistry. And because meteorological reanalyses usually drive CTMs, they may easily be matched in space and time to observations, unlike CCMs or ESMs that generate their own winds. However, the reliance of CTMs on existing driving meteorology means that CTMs have traditionally been largely excluded from international assessments aimed at forecasting future or past changes, such as those of the ongoing phase (phase 6) of the Coupled Model Intercomparison Project (CMIP6, Eyring et al.2016) that is set to inform the upcoming Intergovernmental Panel on Climate Change (IPCC) Sixth Assessment Report.

Here, we introduce, describe, and evaluate version 2.0 of the Global Change and Air Pollution (hereafter, “GCAP 2.0”) chemical-transport model framework. GCAP 2.0 represents a major update and expansion of the original GCAP described by Wu et al. (2007) and Murray et al. (2014). Meteorology necessary for driving the grassroots-community GEOS-Chem 3-D chemical-transport model (, last access: 9 September 2021) has been archived for the pre-industrial, recent past, and several future scenarios of the CMIP6 experiment using version E2.1 of the NASA Goddard Institute for Space Studies (GISS) GCM (Kelley et al.2020; Miller et al.2021). In addition, the CMIP6 emissions and surface boundary conditions have been processed for use within GEOS-Chem for consistency with the driving meteorology and to enable GEOS-Chem to perform and contribute to the CMIP6 experiments.

Section 2 summarizes the history of the GCAP framework. Section 3 describes the climate and chemistry models used and their interface. Section 4 summarizes and evaluates the meteorology products in the recent past versus reanalyses. Section 5 describes the emissions and boundary conditions and evaluates the climate-sensitive emissions. Section 6 evaluates the model in the recent past by comparing it to observations. We conclude with a summary section.

2 History

The GISS GCM and GEOS-Chem CTM have a long history of collaborative development. The immediate predecessor to GEOS-Chem was a gas-phase CTM of tropospheric ozone–NOx–CO–hydrocarbon chemistry (Wang et al.1998a, b, c; Wang and Jacob1998) driven by present-day meteorology archived from version II' of the GISS GCM at 4 latitude by 5 longitude horizontal resolution with seven vertical layers extending from the surface to 150 hPa (Hansen et al.1983; Rind and Lerner1996). GEOS-Chem was born when Bey et al. (2001) updated this model to include the tropospheric non-methane hydrocarbon oxidation mechanism of Horowitz et al. (1998) and allowed it to be driven instead by the Goddard Earth Observing System (GEOS) assimilated meteorological reanalyses produced by the NASA Global Modeling and Assimilation Office (GMAO) (Schubert et al.1993). This early version of GEOS-Chem incorporated the GEOS dynamical core (Lin and Rood1996). Shortly afterward, a bulk sulfur–nitrate–ammonium aerosol mechanism was included by Park et al. (2004).

Since these origins, GEOS-Chem has developed a large user and active developer base of hundreds of individuals at more than 150 institutions in over 30 countries (, last access: 9 September 2021). The model is extensively versioned, documented, and benchmarked. Other subsequent major developments include the development of one- and two-way coupled nested regional simulations (Wang et al.2004; Yan et al.2016; Bindle et al.2020), an adjoint for inverse model applications (Henze et al.2007; Kopacz et al.2009), a unified chemical mechanism from the surface to the mesopause (Eastham et al.2014), a flexible emissions pre-processor (Keller et al.2014), and a massively parallel distributed computing framework enabling global simulations down to resolutions of 0.25 latitude by 0.3125 longitude with 72 vertical layers extending to 0.01 hPa (Eastham et al.2018).

The Global Change and Air Pollution (GCAP) framework developed by Wu et al. (2007) re-enabled version 7-02-04 of GEOS-Chem to be driven by GISS meteorology in order to explore how changes in future climate and precursor emissions may influence surface air quality (e.g., Wu et al.2008a, b; Pye et al.2009; Pye and Seinfeld2010; Selin et al.2009; Hui and Hong2013; Zhu et al.2017). GCAP utilized meteorology archived from version III of the GCM (Rind et al.2007) at 4 latitude by 5 longitude resolution with 23 vertical layers extending to 0.002 hPa for the present-day and the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES) “A1B” scenario for 2050 CE (Nakicenovic and Swart2000).

In the subsequent ICE age Chemistry And Proxies (ICECAP) project, Murray et al. (2014) updated the GCAP implementation to enable version 9-01-03 of GEOS-Chem to be driven by paleometeorology archived from version E of the GISS GCM (Schmidt et al.2006) and consistent land cover simulated using terrestrial vegetation models (Kaplan et al.2006; Pfeiffer et al.2013) to explore chemistry–climate changes at and since the Last Glacial Maximum (LGM;  21 kyr before present; Murray et al.2014; Achakulwisut et al.2015; Geng et al.2015, 2017).

However, the original GISS-driven variants of GEOS-Chem suffered from several issues. Most notably, the stratosphere-to-troposphere mass flux was always too large, complicating the tropospheric ozone budget and the interpretation of polar ice-core records once GEOS-Chem developed online interactive stratospheric chemistry. The use of the more accurate, but computationally expensive, GISS dynamical core within GEOS-Chem to improve transport yielded severe performance issues in the CTM. At the time, both GEOS-Chem and the GISS GCM used their own in-house binary formats for file input and output that required translation (versus the standard NetCDF file format used today by both models). Lastly, the different horizontal and vertical resolutions required extensive offline processing of input fields. GCAP was eventually deprecated and removed from the GEOS-Chem codebase in version 11-02d.

However, subsequent developments to both models increased flexibility and capabilities, motivating the development of GCAP 2.0, as described in the subsequent section.

3 Model description

GCAP 2.0 is the second generation of a one-way offline coupling between the NASA GISS GCM and the GEOS-Chem CTM. Meteorology archived from version E2.1 of the GCM for any period of Earth history or its future may be used to drive the GEOS-Chem CTM. The following subsections describe the salient components and edits to the GCM and CTM relevant for GCAP 2.0 simulations.

3.1 NASA GISS ModelE2.1

The version of the GISS GCM frozen and applied to the initial CMIP6 experiments (ModelE2.1, hereafter “E2.1”) is described in detail by Kelley et al. (2020) and Miller et al. (2021). In brief, the standard E2.1 configuration resolves the equations of mass, momentum, and energy in Earth's atmosphere at a horizontal resolution of 2 latitude by 2.5 longitude and with 40 vertical layers extending from the surface to 0.1 hPa ( 28 in the tropical troposphere; see Fig. 1). The model employs a quadratic-upstream scheme for advection that yields finer effective spatial resolutions by transporting higher-order moments of the subgrid distributions (Prather1986). Gravity-wave momentum fluxes resulting from flow over topography and fronts are parameterized as stratospheric drag processes (Rind et al.1988). Moist convection underwent substantial updates relative to E2.1's predecessor version, E2 (Kim et al.2012; Del Genio et al.2012, 2015). Radiation physics includes calculations for major shortwave and longwave absorbers (water vapor, carbon dioxide, ozone, methane, nitrous oxide, chlorofluorocarbons) and aerosol particles (Hansen et al.1983), any of which may be either prescribed or calculated online as a function of emissions, chemistry, and physical losses. The influence of aerosol particles on cloud microphysics and albedo may be explicitly represented or parameterized (Bauer et al.2020). The model may be coupled to a fully interactive ocean model or be applied in atmosphere-only mode through prescribed sea-surface temperatures.

Figure 1Comparison of the vertical resolutions of GEOS-Chem driven by the MERRA-2 or GEOS-FP reanalyses (full and the reduced stratosphere; orange), the original GCAP driven by Model III or ICECAP driven by ModelE (yellow), and GCAP 2.0 driven by E2.1 (red).


GISS contributed several configurations of E2.1 to CMIP6. Here, we use the atmosphere-only configuration with composition prescribed from earlier runs using online and interactive chemistry for computational expediency. At the time of publication, GISS also contributed up to 11 ensemble members per historical and future emission scenario initialized from different moments of the pre-industrial control simulation. We focus on the atmosphere-only ensemble member that contributed to the largest number of Tier 1 and Tier 2 scenarios of the CMIP and Scenario Model Intercomparison Project (ScenarioMIP) experiments, corresponding to variant label “r1i1p1f2” in the CMIP6 data repository (, last access: 9 December 2020).

On top of the E2.1 codebase used for the CMIP6 simulations, we implemented new subdaily diagnostics that archive the same fields used to drive GEOS-Chem as generated by the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) meteorological reanalysis (Gelaro et al.2017). We then re-performed the “r1i1p1f2” variant of the E2.1 contributions to the CMIP and ScenarioMIP experiments using initial and intermediate restart files archived during the original simulations and archiving the subdaily diagnostics necessary for driving GEOS-Chem. We discuss and evaluate the meteorology in Sect. 4. Three-dimensional fields were archived at 3 h temporal resolution and two-dimensional fields were archived at hourly temporal resolution, consistent with the MERRA-2 product. In addition, we archived hourly lightning flash densities and convective cloud depths. The only difference in the repeat simulation configurations with respect to their original runs was a need to call the radiation code every dynamic time step instead of every five to obtain the hourly radiation fields necessary for driving GEOS-Chem; the consequences of this are discussed in Sect. 4.

Table 1GCAP 2.0 meteorology available at the time of publication from the GCAP 2.0 data repository (, last access: 9 September 2021).

All meteorology fields available at 2 latitude by 2.5 longitude with 40 vertical layers from the surface to 0.1 hPa. Two-dimensional fields are archived at hourly temporal resolution. Three-dimensional fields are archived at 3 h temporal resolution.

Download Print Version | Download XLSX

Table 1 summarizes the 178 years of GCAP 2.0 meteorology archived and publicly available at publication time. For comparison, MERRA-2 presently has 41 years of complete meteorology available. The data are publicly served from the new GCAP data server hosted by the University of Rochester Atmospheric Chemistry and Climate Group at (last access: 9 September 2021). Users can point to this repository analogously to the existing GEOS-Chem data servers hosted by Harvard University (, last access: 9 September 2021) or Compute Canada (, last access: 9 September 2021). Historical meteorology has been archived for the pre-industrial era (1851–1860 CE) and recent past (2001–2014 CE). In addition, we archive near-future (2040–2049 CE) and end-of-the-century (2090–2099 CE) meteorology for seven future scenarios ranging from extreme mitigation to extreme warming (see Sect. 5 for a description of the emission scenarios). In addition, to facilitate comparison of GCAP 2.0 meteorology and composition with observations and traditional GEOS-Chem, we have also performed a recent past simulation in which the E2.1 horizontal winds of the r1i1p1f2 variant were “nudged” to match those of the MERRA-2 reanalysis for 2001–2014 CE (Menon et al.2008). Note that we only nudge the winds and not temperature, humidity, or surface pressure as may be done in other models. We urge users to be aware of the challenges involved when interpreting the impact of nudged meteorology on atmospheric composition, especially in the stratosphere (e.g., see Orbe et al.2020a).

3.2 GEOS-Chem

GEOS-Chem (, last access: 9 September 2021) is a global or regional 3-D chemical transport model traditionally driven by assimilated meteorology products produced by the NASA Global Modeling and Assimilation Office (GMAO) Goddard Earth Observing System Data Assimilation System (GEOS-DAS). These include the MERRA-2 science product, which is generated at 0.5 latitude by 0.625 longitude and 72 vertical layers extending from the surface to 0.01 hPa ( 38 layers in the tropical troposphere) and available from 1980 CE to the present (Gelaro et al.2017). There is also a near-real-time product (GEOS-FP) available at 0.25 latitude by 0.3125 longitude horizontal resolution and available from 2012 CE, although with periodic changes to the underlying code. Both products are provided at hourly temporal resolution for two-dimensional fields and at 3 h resolution for three-dimensional fields1. Most GEOS-Chem users make use of these fields that have been pre-processed to coarser horizontal resolutions (4 latitude by 5 longitude or 2 latitude by 2.5 longitude)2 for computational expediency and to minimize storage requirements. Users may also select to run with reduced vertical resolution in the stratosphere (see Fig. 1).

Emissions are the subject of Sect. 5. The original description of the tropospheric chemical mechanism is by Bey et al. (2001), which was expanded to include a stratospheric mechanism by Eastham et al. (2014), the latter of which did not exist in the earlier versions of GCAP. The coupled sulfur–nitrate–ammonium aerosol simulation is described by Park et al. (2004) with aerosol thermodynamics computed via the ISORROPIA II model (Fountoukis and Nenes2007). Advection is handled by a flux-form and partially semi-Lagrangian transport scheme (Lin and Rood1996). Convective transport is parameterized as a single plume acting under the mean upward convective, entrainment, and detrainment mass fluxes for each level of a model column as archived from the GCM.

GEOS-Chem developed the capability to be driven by any horizontal resolution beginning with the “FlexGrid” update in version 12.4.0 (, The International GEOS-Chem User Community2019​​​​​​​). The model can now define any horizontal resolution at initialization and automatically regrid input meteorology upon file read from its archived resolution to the runtime resolution. Because our strategy was to archive from E2.1 the same fields as in the MERRA-2 reanalysis, very few modifications were necessary to the GEOS-Chem source code to allow GEOS-Chem to use E2.1 output as a meteorological driver. The primary additional code required is the inclusion of the specification of the E2.1 vertical resolution. Otherwise, the only other GEOS-Chem changes were removing some hard-coded limitations, e.g., those that prevented the model from running on dates before 1 January 1900. These updates entered the standard GEOS-Chem code in version 13.1.0 (, The International GEOS-Chem User Community2021a). Version 13.0.0 of GEOS-Chem (, The International GEOS-Chem User Community2021c) introduced the ability of the source code to generate run directories, and version 13.1.0 has been updated to do so for GCAP 2.0. We have regridded all restart files to the new 40-layer vertical resolution.

Because of the relatively few required changes, GCAP 2.0 meteorology is easily compatible with any existing GEOS-Chem capability. There are three primary methods by which GEOS-Chem may be used. The first and most common method due to its ease of installation and application is GEOS-Chem Classic or “GCClassic”. Therefore, we have guaranteed that all GCClassic configurations work with GCAP 2.0 by including run directories and regridding all input files that have a vertical dimension. In addition to full-chemistry simulations, these include the speciality simulations (e.g., offline aerosol, methane, carbon dioxide, tagged CO, tagged methane, tagged ozone). Only the existing GCClassic mercury simulation will require some modifications; in the interest of storage, we did not archive the 10 extra sea-ice fields only used by that simulation since they may be determined online from the fraction of sea-ice field that was archived. FlexGrid also enables one to perform a global simulation at a relatively coarse resolution to archive boundary conditions for driving nested regional simulations at higher spatial resolution. Although the underlying meteorology would still be the 2 latitude by 2.5 longitude of the GCAP 2.0 meteorology, one can benefit from the finer spatial resolution of the emission inventories.

The second method of running GEOS-Chem is the Message Passing Interface (MPI) parallelized variant utilizing a cubed-sphere dynamical core known as GEOS-Chem High-Performance (GCHP, Eastham et al.2018). GCAP 2.0 meteorology is fully compatible with GCHP, although we refer the reader to Sect. 5.2 about the necessary pre-processing of emissions for GCAP 2.0 runs using GCHP.

Lastly, there exists an adjoint of GCClassic used for inverse modeling and sensitivity applications (Henze et al.2007). Since the adjoint presently works with MERRA-2 meteorology, the GCAP 2.0 meteorology is also compatible with the adjoint code once the vertical resolution is added, allowing for inverse modeling applications in past and future climates.

4 Meteorology

This section evaluates the GCAP 2.0 meteorology by comparing it to its original CMIP6 simulation, the CMIP6 E2.1 ensemble, and the MERRA-2 reanalysis. Model output contributed to the CMIP6 experiment is archived by an international distributed data repository powered by the Earth System Grid Federation (ESGF) and available online at (last access: 9 September 2021).

Figure 2Temporal evolution of global mean surface temperature in C. Each panel from left to right shows the seven future scenarios arranged by increasing radiative forcing. The black line shows variable tas from the E2.1 ensemble member (r1i1p1f2) that our simulations are based upon obtained from the ESGF repository. The gray shading represents the annual mean and 2σ spread of all E2.1 ensemble members archived on ESGF. The red line shows the value from our reruns of r1i1p1f2 to generate the GCAP 2.0 meteorology. The blue line shows the same but with winds nudged to the MERRA-2 reanalysis. The orange line shows the global mean surface temperature from the MERRA-2 reanalysis. The horizontal lines show the global pre-industrial climatological mean in our simulations (dotted black) and 1.5 C (dashed black) and 2.0 C (solid black) increases on top of the pre-industrial values.


Figure 2 shows the temporal evolution of annual mean surface air temperature. The black line shows the E2.1 r1i1pif2 variant and the gray shading shows the mean and 2σ spread of the submitted E2.1 ensemble from ESGF (ranging from 11 members in the historical to 1 member in some future scenarios). Our repeat simulations of the r1i1p1f2 variant from its archived restart files are shown in red and the same variant nudged to MERRA-2 is shown in blue. The MERRA-2 historical record is shown in orange.

First, we note that the repeat simulations have slightly warmer surface air temperatures than the original simulations. A small portion of this difference results from numerical noise associated with the original and repeat simulations' different computer architectures (NASA Center for Climate Simulations versus the University of Rochester Center for Integrated Research Computing). However, tests revealed that the increased frequency in calls to the radiation code necessary to archive hourly radiation fluxes for driving GEOS-Chem explains most of the difference in surface temperature. Although almost no locations show significant changes with respect to local interannual variability (see Fig. S1 in the Supplement), this does lead to a weakly but statistically different annual mean temperature with respect to the original simulation (+0.6 C for 2001–2014 CE, p value = 0.002). However, the change is largely a linear offset, with temporal correlation remaining high (R=0.87 for 2001–2014 CE), providing confidence in the ability of GCAP 2.0 to produce changes in composition associated with changes in climate accurately. Nudging the winds to MERRA-2 reduces this offset by influencing the rate of mixing of air between the high latitudes and midlatitudes, and the warmer Arctic surface temperatures in our repeat free-running E2.1 simulations are in greater agreement with MERRA-2.

Second, we note that for researchers interested in studying a “Paris Agreement”-like future world, in which future warming is limited to 2 C over pre-industrial levels, one may use the SSP1-1.9, SSP1-2.6, or SSP4-3.4 scenarios. However, if one wishes to study a future world with the more aggressive goal of limiting warming to 1.5 C over pre-industrial levels, then the only scenario that may be used is SSP1-1.9.

Figure 3The same as Fig. 2 but showing the temporal evolution of the global mean precipitation rate in mm d−1 (variable pr on ESGF).


Figure 3 shows the temporal evolution of the global annual mean precipitation rate in the simulations. Global precipitation rates are forecast to increase in the coming century due to the temperature-driven increase in surface evaporation and saturation vapor pressures. Unlike surface air temperature, the repeat simulations closely follow the original values and are statistically identical except in the recent past historical simulation, where they are globally higher by 0.9 %. The MERRA-2 reanalysis shows substantially more interannual variability in its global precipitation rates. Nudging E2.1 decreases the global precipitation rate.

Figures S2 to S51 of the Supplement include detailed comparisons of the seasonal climatologies for all fields in the three meteorology products that may be used to drive GEOS-Chem for 2005–2014 CE (MERRA-2, E2.1 nudged to MERRA-2, and the free-running E2.1). In general, most fields show excellent agreement with high pattern (i.e., spatial) correlation and small mean difference. However, a few fields differ between E2.1 and MERRA-2 that are of interest for chemical-transport modeling, which we now summarize. The primary difference between MERRA-2 and E2.1 is the relative importance of stratiform versus convective precipitation. Whereas both models agree on the total precipitation flux to the surface (Fig. S16), MERRA-2 has a higher rate of stratiform condensation (Fig. S33) balanced by a higher rate of stratiform re-evaporation (Fig. S45). In contrast, E2.1 has a higher rate of convective condensation (Fig. S34) balanced by a higher rate of convective re-evaporation (Fig. S44). Furthermore, E2.1 has consistently smaller surface roughness heights over the ocean and vegetated regions than MERRA-2; in contrast, MERRA-2 does not appear to include an orographic component in its surface roughness calculation and therefore has lower surface roughness heights over non-vegetated land surfaces (Fig. S29). Consequently, the E2.1 simulations have lower planetary boundary layer heights over oceans and heavily vegetated regions (globally  200 m lower; Fig. S12) relative to MERRA-2. The E2.1 simulations also have a lower tropopause pressure by approximately 40 hPa (Fig. S23); this has been corrected in version E2.2 of the GCM by moving to a higher vertical resolution (Orbe et al.2020b). Lastly, E2.1 has a higher fraction of photosynthetically active radiation (PAR) present as diffuse (Fig. S10) rather than direct (Fig. S11) radiation, which will promote higher levels of biogenic emissions (see Sect. 5.2). Note also that MERRA-2 sets PAR fluxes to zero over water. Therefore, coastal and island cells will underestimate the radiation flux in MERRA-2-driven GEOS-Chem simulations, again with consequences for biogenic emissions.

Figure 4Comparison of present-day surface meteorology. Columns from left to right show annual climatological means for 2005–2014 CE in the MERRA-2 reanalysis, our E2.1 simulations with winds nudged to MERRA-2, and our free-running E2.1 simulations. From top to bottom, rows show surface air temperature in K, total precipitation rate in mm d−1, the zonal component of the surface wind in m s−1, and the meridional component of the surface wind in m s−1. Gray dots show where the two E2.1 simulations are statistically different from their MERRA-2 counterparts with respect to interannual variability (p value < 0.05; n= 10 years). The number in the bottom left shows the global mean value for each panel. The top right number shows the pattern correlation, and the number in the bottom right shows the global mean difference in the E2.1 simulations with respect to MERRA-2.

Figure 4 compares the spatial distribution of key surface meteorological variables generated for GCAP 2.0 with their MERRA-2 counterparts. Surface air temperature shows near-perfect agreement in spatial distribution between the E2.1 products and MERRA-2. The E2.1 temperatures are slightly higher than MERRA-2, especially over the Northern Hemisphere's oceans; nudging reduces this difference as previously discussed. Total precipitation in the E2.1 fields has a weaker pattern correlation with MERRA-2 since the free-running model produces a split Intertropical Convergence Zone (ITCZ) in the eastern Pacific (a common issue in free-running GCMs; e.g., see Samanta et al.2019); nudging the winds corrects the spatial patterns but brings the total precipitation rate out of agreement. The surface zonal wind component shows excellent agreement in their spatial patterns, although the magnitudes are greater in the E2.1 simulations relative to MERRA-2. The free-running GCM underestimates the extent of flow towards the Equator over the eastern ocean basins, which is corrected in the nudging simulation (and may be responsible for the improved ITCZ). The E2.1 simulations also lack the relatively large spatial heterogeneity seen in surface winds over the land ice sheets.

Figure 5Comparison of present-day zonal meteorology. Columns from left to right show annual climatological zonal means for 2005–2014 CE in the MERRA-2 reanalysis, our E2.1 simulations with winds nudged to MERRA-2, and our free-running E2.1 simulations. From top to bottom, rows show air temperature in K, specific humidity in kg kg−1, and the zonal wind component in m s−1. Gray dots show where the two E2.1 simulations are statistically different from their MERRA-2 counterparts with respect to interannual variability (p value < 0.05; n= 10 years).


Figure 5 compares key zonal mean meteorological variables generated for GCAP 2.0 with their MERRA-2 counterparts. Lower- and free-tropospheric air temperatures are in agreement between E2.1 and MERRA-2, but the higher tropopause leads to colder temperatures in the upper troposphere with respect to MERRA-2. Nudging the winds removes some of the temperature difference in the extratropical upper troposphere but introduces differences in the free troposphere. Specific humidity agrees between E2.1 and MERRA-2 except for a drier polar free troposphere and northern extratropical surface; nudging leads to an additional drying of the stratosphere. The zonal winds agree well between all simulations, particularly between the MERRA-2 reanalysis and the nudged simulation as to be expected.

5 Emissions and boundary conditions

This section describes the anthropogenic emission inventories and surface boundary conditions from the CMIP6 experiment that have been processed for use by GEOS-Chem, either driven by E2.1 or MERRA-2 meteorology (Sect. 5.1). It then evaluates and compares the emission fluxes that are sensitive to meteorology between MERRA-2 and E2.1-driven GEOS-Chem simulations in the recent past (Sect. 5.2).

5.1 Anthropogenic emissions and surface boundary conditions

Figure 6 shows the time series of annual mean anthropogenic emissions and Fig. 7 shows the time series of annual mean surface boundary conditions developed for the CMIP6 experiments and processed for use in GEOS-Chem. Emissions are used for short-lived climate forcers and air pollution precursors. Surface boundary conditions are used to prescribe long-lived species like chlorofluorocarbons that are well mixed in the troposphere but may advect to and react within the stratosphere. The emissions and boundary conditions developed for the CMIP6 experiments include a historical reconstruction and several future scenarios with different target radiative forcings for the end of this century. They are hosted on ESGF under the input data sets for Model Intercomparison Projects (input4MIPs) project. They are available from (last access: 30 September 2020​​​​​​​). These emissions and boundary conditions have been processed for input to GEOS-Chem/GCAP 2.0. They are consistent with those that influenced the climate of the E2.1 simulations used to generate the respective GCAP 2.0 meteorology. When users generate a GCAP 2.0 run directory, the respective CMIP6 emissions and boundary conditions for a historical or future scenario are enabled by default. However, users may always modify the Harmonized EMissions COmponent (HEMCO, Keller et al.2014) configuration to use any emissions they desire (e.g., to use the alternative inventories, regional overwrites and/or scaling factors from the default GEOS-Chem configuration).

Figure 6Time series of annual mean CMIP6 anthropogenic emissions for 1850–2100 CE for (a) reactive oxides of nitrogen (NOx NO + NO2) in Tg N yr−1, (b) carbon monoxide (CO) in Tg yr−1, (c) non-methane hydrocarbons (NMHCs) in Tg yr−1, (d) sulfur dioxide (SO2) in Tg yr−1, (e) ammonia (NH3) in Tg yr−1, (f) black carbon (BC) in Tg C yr−1, and (g) organic carbon (OC) in Tg C yr−1. Historical emissions for 1850–2014 CE from Hoesly et al. (2018) are shown as black lines. Future scenarios for 2015–2100 CE from Gidden et al. (2019) are shown as colored lines: SSP1-1.9 (green); SSP1-2.6 (blue); SSP4-3.4 (purple); SSP2-4.5 (brown); SSP4-6.0 (orange); SSP3-7.0 (pink); SSP5-8.5 (red). The shaded blue rectangles indicate periods for which GCAP 2.0 input meteorology is available at the time of publication.


Figure 7Time series of annual mean CMIP6 surface boundary conditions for 1850–2100 CE. Individual panels show the temporal evolution of (a) carbon dioxide (CO2) in ppmv (µmol mol−1), (b) methane (CH4) in ppbv ( nmol mol−1), (c) total chlorofluorocarbons (CFCs  CFC11 + CFC12 + CFC113 + CFC114 + CFC115) in pptv ( pmol mol−1), (d) nitrous oxide (N2O) in pptv, (e) methyl chloride (CH3Cl) in pptv, (f) total hydrochlorofluorocarbons (HCFCs  HCFC141b + HCFC142b + HCFC22) in pptv, (g) methyl chloroform (CH3CCl3) in pptv, (h) carbon tetrachloride (CCl4) in pptv, (i) dichloromethane (CH2Cl2) in pptv, (j) chloroform (CHCl3) in pptv, (k) methyl bromide (CH3Br) in pptv, and (l) total halons ( Halon1211 + Halon1301 + Halon2402) in pptv. Historical boundary conditions are from Meinshausen et al. (2017) and future boundary conditions are described by Riahi et al. (2017). The shaded blue rectangles indicate periods for which GCAP 2.0 input meteorology is available at the time of publication.


5.1.1 Historical

Historical anthropogenic emissions in CMIP6 are from the Community Emissions Data System (CEDS, Hoesly et al.2018). For surface emissions, we processed version 2017-05-18 of the CEDS inventory available at monthly temporal and 0.5 spatial resolution for the eight surface sectors listed in Table 2 and the 1850–2014 CE period. CEDS is presently the default global surface anthropogenic emissions inventory used by GEOS-Chem, although only species for the full-chemistry simulation have been processed, and these emissions are overwritten for many locations by regional inventories. Here, we additionally processed the methane and CO2 fluxes for those GEOS-Chem specialty simulations. All surface emissions increased exponentially over the historical period except for sulfur dioxide (SO2), whose emissions peaked in the 1980s (Fig. S52). We also processed the three-dimensional CEDS aircraft emissions for input into GEOS-Chem as a CMIP6-compliant alternative to the Aviation Emissions Inventory Code (AEIC, Stettler et al.2011) source that is the default in GEOS-Chem. We processed version 2017-08-30 of the CEDS aircraft inventory, available for NO, CO, black and organic carbon, SO2, and ammonia (NH3) at monthly temporal and 0.5 spatial resolution and 25 vertical levels of equal thickness from the surface to 15 km. We have vertically regridded to the native 40-level E2.1 resolution and the reduced stratospheric 47-level MERRA-2 resolution. Global aircraft emissions increased mostly linearly across the historical period beginning in approximately 1950 CE (Fig. S53).

Table 2CEDS surface emission sectors.

Download Print Version | Download XLSX

In addition, we processed the biomass burning emissions from version 1.2 of the Biomass Burning for Model Intercomparison Projects (BB4MIPs) inventory (van Marle et al.2017). The BB4MIPs reconstruction combines version 4 of the satellite-based Global Fire Emissions Database (GFED4; van der Werf et al.2017, the default biomass burning inventory in GEOS-Chem and available since 1997 CE) with observational proxies and model simulations from the Fire Model Intercomparison Project (FireMIP) for earlier periods. BB4MIPs provides the total mass flux per species at monthly temporal and 0.25 spatial resolution. We have regridded to 0.5 spatial resolution and speciated for input to GEOS-Chem using a consistent hydrocarbon speciation scheme with CEDS. Fire emissions from BB4MIPs increase slightly across the historical period. However, interannual variability greatly increases in the second half of the 20th century (Fig. S54), driving the large interannual variability observed in the total emissions during this period (e.g., Fig. 6b).

We have prepared version 1.2.0 of the CMIP6 surface boundary conditions derived from historical observations (Meinshausen et al.2017) and available at monthly temporal resolution and as 0.5 latitude bands. These have been regridded to 0.5 global spatial resolution for input into GEOS-Chem. Carbon dioxide, methane, and nitrous oxide have monotonically increased since the pre-industrial era due to anthropogenic activity, whereas shorter-lived stratospheric ozone-depleting substances peaked around the turn of the last century following their ban under the Montreal Protocol (Fig. 7).

5.1.2 Future scenarios

The future anthropogenic emissions and boundary conditions used by CMIP6 and processed here for GCAP 2.0 are summarized by Riahi et al. (2017). In brief, the so-named Shared Socioeconomic Pathways (SSPs) are determined from integrated assessment modeling (IAM) of five future societal narratives that may be followed to limit future warming to a target radiative forcing. The nomenclature for characterizing the SSP scenarios is SSPx-y.z (or SSPxyz), where x is the number of the future narrative (1 to 5; Table 3), and y.z represents the target radiative forcing in W m−2 at the end of the 21st century ranging from 1.9 (extreme mitigation; low warming) to 8.5 (low mitigation; extreme warming). The names of the five narratives are listed in Table 3, and each employs different assumptions about how global society would achieve the target radiative forcing. SSP1 assumes low challenges to mitigation and adaptation (van Vuuren et al.2017), SSP2 assumes medium challenges to mitigation and adaptation (Fricko et al.2017), SSP3 assumes high challenges to mitigation and adaptation (Fujimori et al.2017), SSP4 assumes low challenges to mitigation and high challenges to adaptation (Calvin et al.2017), and SSP5 assumes high challenges to mitigation and low challenges to adaptation (Kriegler et al.2017).

Table 3Shared Socioeconomic Pathway (SSP) narratives.

Download Print Version | Download XLSX

For GCAP 2.0, we focus on seven scenarios corresponding to Tiers 1 and 2 of the ScenarioMIP experiment: SSP1-1.9, SSP1-2.6, SSP4-3.4, SSP2-4.5, SSP4-6.0, SSP3-7.0, and SSP5-8.5 (O'Neill et al.2016). Assumptions about future population growth, urbanization, gross domestic production (GDP), energy, land use, and air pollution trends of the SSP IAM scenarios are described in a series of papers (Crespo Cuaresma2017; Dellink et al.2017; Jiang and O'Neill2017; Leimbach et al.2017; Samir and Lutz2017; Bauer et al.2017; Popp et al.2017; Rao et al.2017). Version 1.1 of these scenarios (Gidden et al.2019) was obtained from input4MIPs on ESGF and available at 0.5 spatial and monthly resolution for 2015 CE and then every 10 years beginning with 2020 CE and ending with 2100 CE. These future emissions were processed for input to GEOS-Chem in the same way as the historical emissions and boundary conditions. Linear interpolation was used to develop individual yearly emissions between the available decadal values. In addition to the sectors of Table 2, future CO2 emissions include a “neg” sector that considers negative emission (i.e., carbon capture).

From the perspective of the simulated climate, the forcing that dominates the end-of-the-century response is the CO2 abundance. Therefore, CO2 is the only gas that monotonically increases with future forcing target values (Fig. 7a). The trajectories of the remainder of the well-mixed greenhouse gases, stratospheric ozone-depleting substances, short-lived climate forcers, and air-pollution precursors vary between the SSP scenarios and target forcings (Figs. 6 and 7b–l). For example, methane in 2100 CE is highest in the SSP3-7.0 scenario, not SSP5-8.5 (Fig. 7b). Because ammonia is primarily produced from agriculture and the world population will continue to grow, it is the only emission expected to remain constant or grow into the future. Otherwise, the narrative assumptions strongly influence the global and regional emission changes. Furthermore, we note that the historical emissions inventories have large amounts of interannual variability (primarily due to biomass burning), whereas the SSP scenarios have very low interannual variability. This highlights the necessity for CTM studies such as those that may be accomplished by GCAP 2.0 that can explore the impact of a wider range of future emission trajectories on air quality, short-lived climate forces, and stratospheric ozone in future warmer climates whose meteorology is primarily driven by changes in CO2. It also highlights the necessity for performing simulations long enough to establish robust statistics for chemistry–climate interactions (e.g., Garcia-Menendez et al.2017).

5.2 Natural emissions

Figure 8 shows the climatological mean emission fluxes for key species whose emissions are sensitive to meteorology. Figures S55–S60 in the Supplement provide seasonal details for each species.

Figure 8Annual mean spatial distribution of meteorology-dependent emission fluxes for 2005–2014 CE. Each column from left to right shows emission fluxes calculated using: MERRA-2 meteorology, E2.1 meteorology nudged to MERRA-2, and the free-running E2.1 meteorology, respectively. Each row from top to bottom shows emission fluxes for: (a–c) isoprene (2-methyl-1,3-butadiene; CH2=C(CH3)CH=CH2) from terrestrial plants in 10−9 kg m−2 s−1, (d–f) dimethylsulfide ((CH3)2S; DMS) from marine organisms, (g–i) aeolian mineral dust, (j–l) aeolian sea salt, (m–o) the vertically integrated source of NO from lightning, and (p–r) NO from soil microbial activity, respectively. Gray dots indicate locations where the two E2.1-driven simulations show statistically significant differences (p value < 0.05; n= 10 years) with respect to the MERRA-2-driven simulation. The value in the lower left of each panel gives the globally integrated source in (a–l) Tg yr−1 or (m–r) Tg N yr−1. The number in the lower (upper) right of each panel gives the total difference (pattern correlation) of the E2.1-driven simulations with respect to their respective MERRA-2-driven values.

Emission fluxes sensitive to meteorology and thereby grid resolution began to be pre-processed offline at high spatial resolution in version 12.4.0 of GEOS-Chem. This was to facilitate the calculation of consistent emissions between the various cubed-sphere geometries of the GCHP variant of the model (Eastham et al.2018, see Sect. 3.2), although the option to calculate these emissions online was maintained through online extensions in the HEMCO processing code (Keller et al.2014). The default behavior for GCAP 2.0 is to use online calculations to respond to the underlying climate. Users who wish to run GCHP with GCAP 2.0 may quickly pre-process offline natural emission fluxes using the “standalone” version of HEMCO with a HEMCO_Config.rc file from a GCClassic version of GCAP 2.0. Here, we compare the online emission fluxes in our MERRA-2-driven simulation to those of our E2.1-driven simulations.

Panels a–c of Fig. 8 show the spatial distribution of isoprene from terrestrial plants. Emissions of non-methane hydrocarbons (NMHCs) from terrestrial plants follow version 2.1 of the Model of Emissions from Gases and Aerosols from Nature (MEGAN), which responds positively to changes in diffuse photosynthetically active radiation (PAR), recent surface air temperature, and soil root wetness (Guenther et al.2012). There is an option for emissions to respond to CO2 abundance as well (Tai et al.2013). Because E2.1 has a greater proportion of PAR present as diffuse radiation, isoprene emissions in E2.1-driven simulations are about 40 % higher than in the MERRA-2-driven simulation.

Panels d–f of Fig. 8 show a very tight agreement in the spatial pattern and magnitude of the flux of dimethylsulfide (DMS; (CH3)2S) produced by marine phytoplankton. Emissions of NMHCs from marine environments are represented as the product of prescribed seawater concentration distributions and sea-to-air transfer velocities calculated via the parameterization of Nightingale et al. (2000a, b). The latter respond to sea-surface temperatures and surface wind velocities.

Panels g–i of Fig. 8 compare the source of mineral dust between the simulations. We use the Dust Entrainment and Deposition (DEAD) scheme for mineral dust evasion (Zender et al.2003), which responds to changes in surface friction velocity (u*), roughness height, snow/ice cover and depth, soil wetness, pressure, specific humidity, and temperature. Dust mobilization was found in our tests to be extremely sensitive to the meteorology product used, with poor spatial correlation (R≤0.26) and with each meteorology product yielding a different order of magnitude in its global total. Therefore, we have determined respective scaling factors for the E2.1 simulations that bring the present-day global total into agreement with the MERRA-2-driven value. These are included by default in the GCAP 2.0 run directories for the DEAD dust scheme.

Panels j–l of Fig. 8 show the source of sea-salt aerosol. The sea-salt mobilization scheme is described by Jaeglé et al. (2011) and responds to sea-surface temperature, surface wind velocity, and the fraction of sea-ice coverage. There is an excellent agreement between each source's spatial distributions, although the stronger surface winds in E2.1 lead to 15 %–20 % higher emissions in the E2.1 simulations, especially over the Southern Ocean.

Panels m–o of Fig. 8 show the column-integrated source of NO from lightning. In MERRA-2, lightning flash densities (flashes km−2 s−1) are pre-calculated offline from MERRA-2 convective cloud depths at 0.5 latitude by 0.625 longitude spatial and 3 h temporal resolution. These densities are then input as a meteorological parameter to GEOS-Chem, from which vertical profiles of NO production are determined following Murray et al. (2012). For the GCAP 2.0 meteorology, flash rates are calculated online in the E2.1 moist convection code following the description in Kelley et al. (2020) and archived at E2.1 native spatial and hourly temporal resolution for input as a meteorological parameter to GEOS-Chem (see Sect. 4). Both lightning flash density calculations are ultimately based on the same cloud-top height scheme of Price and Rind (1992) and global mean lightning flash rates are tuned to climatology in the recent past (Cecil et al.2014). However, because the spatial and seasonal climatology in the MERRA-2-driven simulations is constrained by satellite observations (Murray et al.2012), which is not appropriate for a free-running GCM, the spatial patterns differ between the simulations. Therefore, E2.1 overestimates the fraction of lightning in the tropics with respect to the extratropics and puts too much lightning over South America and Oceania and not enough over Africa. It is also worth emphasizing that we do not know how lightning has changed since the pre-industrial era or will change in a warming world (e.g., Williams2005; Price2013; Murray2016, 2018; Finney et al.2018).

Panels p–r of Fig. 8 show the source of NO from soil microbial activity. The parameterization is described by Hudman et al. (2012) and responds to surface air temperature, wind speed, soil wetness, cloud fraction, downwelling shortwave radiation, and snow/ice cover. To a lesser degree, lightning can also influence the soil NO source through its impact on nitrate deposition to the soils. The spatial correlation is excellent between the different sources, although the E2.1 magnitude is higher by about 40 %.

Lastly, we note the important and variable geologic source of SO2. Volcanic emissions of SO2 in GEOS-Chem are normally prescribed from the Aerosol Comparisons between Observations and Models (AeroCom) point-source inventory (Carn et al.2015), with data available since 1978 CE. The CMIP6 experiment did not provide historical or future emission fluxes for volcanism. Instead, input4MIPs provided time series of stratospheric aerosol surface area densities and effective radii with which to force the GCMs. Therefore, when users generate a GCAP 2.0 run directory, they are given the option to select a fixed historical AeroCom year from which to prescribe their volcanic emissions.

6 Model evaluation

This section evaluates the performance of GCAP 2.0 driven by E.21 versus MERRA-2 meteorology for the recent past through comparison with observations. We first evaluate model physics and transport using the “TransportTracers” variant of GEOS-Chem (Sect. 6.1). We then evaluate the standard full-chemistry mechanism (Sect. 6.2).

All simulations were performed at 4 latitude by 5 longitude horizontal resolution for the period 2001–2014 CE, with meteorology respectively prescribed from MERRA-2, E2.1 nudged to MERRA-2, and the free-running E2.1 simulation. All simulations used version 12.9.3 of GEOS-Chem (, The International GEOS-Chem User Community2020) with modifications as described throughout the paper. The E2.1 meteorology fields were regridded from their native resolution upon input to GEOS-Chem by FlexGrid. The E2.1 simulations used the native 40-layer resolution and the MERRA-2 simulations used the 47-layer reduced-stratospheric resolution (see Fig. 1). Identical initial conditions were regridded to each model's respective vertical resolution. The first 4 years of each simulation were discarded as initialization, with the remaining 10 years used for evaluation and statistics. All prescribed emissions were identical between each simulation.

6.1 Transport

Model transport and physical processes may be evaluated against observations using the “TransportTracers” variant of GEOS-Chem. We focus on four tracers of particular utility: sulfur hexafluoride (SF6), radon-222 (222Rn), lead-210 (210Pb), and beryllium-7 (7Be).

Sulfur hexafluoride is a trace gas of anthropogenic origin that is chemically and physically inert on human timescales (lifetime of 3200 years). It is primarily emitted at the surface in the Northern Hemisphere (Maiss and Brenninkmeijer1998). Its meridional gradient may be used to test the rate of interhemispheric mixing (Rigby et al.2010; Hall et al.2011) and its vertical gradients may be used to infer the age of air in the stratosphere (Waugh and Hall2002; Waugh2009). Its meridional gradients may also be used to infer the tropospheric age of air (Waugh et al.2013). Here, we use emissions from version 4.2 of the Emissions Database for Global Atmospheric Research (EDGAR), available at 0.1 global resolution for 1970–2008.

Terrigenic 222Rn is an inert, insoluble, short-lived (half-life of 3.8 d) noble gas produced from the slow decay of 226Ra (half-life of 1600 years) found in uranium ores. Its evasion from surface soils is relatively uniform and constant and is as described by Jacob et al. (1997). Its insolubility and timescale of decay make it a useful tracer for diagnosing quick vertical mixing within atmospheric models from boundary layer processes and moist convection (e.g., Allen et al.1996; Brost and Chatfield1989; Considine et al.2005; Feichter and Crutzen1990; Hauglustaine et al.2004; Jacob and Prather1990; Jacob et al.1997; Lambert et al.1982; Mahowald et al.1995; Stockwell et al.1998).

Radiogenic 210Pb is the chemically inert decay product of 222Rn. It is readily taken up by submicron aerosol particles and subsequently removed from the atmosphere by deposition or decay (Bondietti et al.1988; Maenhaut et al.1979; Preiss et al.1996; Sanak et al.1981). Because of its relatively long lifetime (half-life of 22.2 years), nearly all 210Pb is removed via deposition. As its source from 222Rn is relatively well known, and there is a global and long-term surface deposition flux inventory (Preiss et al.1996), it is the standard test for model deposition.

Cosmogenic 7Be is produced by cosmic-ray spallation of N2 and O2, predominantly in the polar upper troposphere and lower stratosphere (Lal et al.1958). The source of 7Be is updated for this work to use the parameterization of Usoskin and Kovaltsov (2008). Mean solar activity is assumed (solar modulation potential Φ= 670 MV), leading to an average production rate of 0.065 atoms cm−2 s−1; about 60 % in the stratosphere and 40 % in the troposphere. Like 210Pb, 7Be is rapidly taken up by submicron aerosol particles (Bondietti et al.1988; Maenhaut et al.1979; Papastefanou2009; Papastefanou and Ioannidou1996; Sanak et al.1981). It is subsequently transported until removal by deposition or radioactive decay (half-life of 53.3 d). 7Be has been used to constrain vertical transport, wet deposition fluxes, and stratosphere–troposphere exchange in models (e.g., Allen et al.2003; Brost et al.1991; Koch et al.1996; Liu et al.2001, 2016; Barrett et al.2012).

Table 4Radionuclide budgets for 2005–2014 CE.

Download Print Version | Download XLSX

Table 4 gives the atmospheric budget of the three radionuclides driven by the three meteorological products.

6.1.1 Horizontal mixing

Figure 9 shows observed meridional and vertical gradients of SF6 with respect to Cape Matatula, American Samoa (SMO), in the remote tropical southern Pacific. The observations are version 2.1.1 of the NOAA Carbon Cycle Group SF6 ObsPack (, NOAA Carbon Cycle Group ObsPack Team2018) and represent a mixture of surface in situ, flask, tower, and aircraft sources from 2005–2014 CE. Observations were aggregated at model spatial and monthly temporal resolution, compared to that month's SMO value, from which zonal climatologies were determined. Also shown is the value of each simulation sampled and processed as in the observations. In all simulations, GEOS-Chem underestimates the cross-equatorial meridional gradient by 17 %–26 %, an improvement over GEOS-Chem driven by earlier meteorology products (see Supplement of Murray et al.2014). However, this suggests that the interhemispheric mixing rate in the model is too fast and/or that the EDGAR inventory underestimates the emission growth rate in the Northern Hemisphere relative to the Southern Hemisphere. Meanwhile, meridional mixing rates in the Southern Hemisphere are consistent with the observations in all simulations. The E2.1 simulation slightly better matches the cross-equatorial gradient than the MERRA-2 and nudged simulations, but E2.1 greatly underestimates the lower stratospheric gradient (see Sect. 6.1.3).

Figure 9Lower atmospheric gradients of SF6 shown as zonal mean percent difference with respect to Cape Matatula, American Samoa (SMO; 14.2 S, 170.6 W; 42 m a.s.l.), for 2005–2014 CE. The left panel shows observed values from version 2.1.1 of the NOAA SF6 ObsPack aggregated at E2.1 vertical resolution. The three panels on the right show the values from the model driven by the three meteorology products sampled at each observation month and location. Gray dots indicate locations where the simulated gradient is statistically different from the observations with respect to interannual variability (p value < 0.05; n=10 years). The number in each model panel's top left shows the pattern correlation (R) of that simulation to the observations. The lower left number shows the mean absolute percent bias of each simulation relative to the observations. The dashed black line shows each simulation's climatological zonal mean tropopause height.


6.1.2 Vertical mixing – troposphere

We assess vertical mixing within the troposphere using vertical profiles of 222Rn and the ratio of 7Be to 210Pb in surface air.

Figure 10 shows simulated climatological 222Rn profiles sampled at the month and location of the available observations, also plotted. Observations are scarce and available only at northern extratropical continental locations (Bradley and Pearson1970; Nazarov et al.1970; Wilkening1970; Moore et al.1973; Kritz et al.1998). In an overly convective atmosphere, the vertical gradient of 222Rn would disappear. There is a slight overestimate of 222Rn abundance within the boundary layer and underestimate above in all of our simulations, implying a small underestimate in boundary layer ventilation. Our results are comparable to or better than other atmospheric models (e.g., see Fig. 5 of Considine et al.2005).

Figure 10Mean observed vertical profile of 222Rn compared to the model sampled at month and location of observations. The units are mBq per standard cubic meter at 0 C and 1 atm, equivalent to a linear transformation of the molar mixing ratio (5.637 mBq SCM−1= 1.0 ×10-22 mol 222Rn (mol air)−1).


Figure 11Surface mixing ratios of 7Be (top row; in mBq SCM−1; 4.05 mBq SCM−1= 1.0 ×10-21 mol 7Be (mol air)−1), 210Pb (middle row; in mBq SCM−1; 2.66 mBq SCM−1= 1.0 ×10-19 mol 210Pb (mol air)−1), and the log10-transformed ratio of the 7Be to 210Pb activities in surface air (bottom row; unitless). The left column shows long-term mean observations from the DOE Surface Air Sampling Program (SASP) program for 1969–1999 CE. The 7Be observations have been selected for periods of average solar activity (Φ=670±50 MV from Usoskin et al.2005). The three columns on the right show (from left to right) values simulated by GEOS-Chem driven by meteorology archived for 2005–2014 CE from MERRA-2, E2.1 nudged to MERRA-2, and the free-running E2.1. The number in each panel's lower left shows the mean value of the observations or models sampled at the observed locations. The top right number shows the pattern correlation (R) of the simulated values with the respective observations. The lower right number shows the mean bias of the simulated values with respect to the observations.

Figure 11 shows the annual mean surface mixing ratios of 7Be, 210Pb and their ratio in our simulations. Since 7Be is produced in the upper troposphere/lower stratosphere and 210Pb is produced near the surface, and because the ratio of 7Be to 210Pb is unaffected by deposition, the ratio serves as a useful measure for tropospheric vertical mixing (Koch et al.1996). A persistent high bias would indicate excessive downward transport and/or insufficient upward transport, assuming no bias in either source. The left column shows the climatological long-term data from the surface monitoring stations of the DOE Environmental Measurements Laboratory (EML) Surface Air Sampling Program (SASP) (obtained from, last access: 15 January 2021​​​​​​​). SASP recorded the spatial and temporal distribution of various radionuclides in surface ambient air from 1957 until 1999 CE, including 7Be and 210Pb. For comparison with our simulations, we select data from periods of average solar activity (solar modulation potential Φ= 670 ± 50 MV from the Usoskin et al.2005 reconstruction). All simulations show only minor biases in surface abundance in either 7Be and 210Pb. The 7Be to 210Pb ratio shows tight agreement between the different model simulations, with the E2.1 simulation slightly outperforming the MERRA-2-driven simulation relative to the observations.

6.1.3 Vertical mixing – stratosphere

Figure 12 shows the simulated zonal climatology of the age of air in the stratosphere, which is defined as the mean time since an air mass at a given location was last in the troposphere (Hall and Waugh2000). Age of air increases away from the equatorial tropopause where most tropospheric air enters the stratosphere (Holton et al.1995). We determine age of air by using SF6 as a chronological tracer (e.g., Waugh and Hall2002) to determine the average temporal lag between a mixing ratio at a given location in the stratosphere relative to the tropical tropopause for the period 2005–2014 CE.

Figure 12Average age of air in the stratosphere in GEOS-Chem simulations driven by MERRA-2, E2.1 nudged to MERRA-2, and E2.1, using the temporal lag in the simulated 2005–2014 CE time series of SF6 relative to the tropical tropopause as a chronological tracer. The dashed black line shows the simulated zonal mean tropopause altitude.


Models traditionally underestimate the stratospheric age of air implied by observations, which increases up to 7 years by 35 km in the poles (Waugh and Hall2002). All three of our simulations, including the MERRA-2 reanalysis, underestimate the age of air. E2.1 has the youngest air, with values consistent with those determined online within the GCM and reported by Orbe et al. (2020b). E2.1 nudged to MERRA-2 has the oldest air but is still only about two-thirds of what observational constraints suggest they should be and remain difficult to interpret given the challenges of nudged simulations (e.g., Orbe et al.2020a). The young age in E2.1 results from too strong an ascent in the tropical pipe and a relatively leaky lower branch of the Brewer–Dobson circulation (e.g., Ray et al.2010); updates to version E2.2 of the GISS GCM greatly improve the stratospheric circulation (Orbe et al.2020b) and will be included in future GCAP 2.0 meteorology products and scenarios.

6.1.4 Stratosphere–troposphere exchange

Beryllium-7 has often been used as a tracer of downward transport from the stratosphere (Dibb et al.1992, 1994; Husain et al.1977; Rehfeld and Heimann1995; Sanak et al.1985; Viezee and Singh1980) and as an indicator of stratosphere–troposphere exchange (STE) performance within global atmospheric models (Allen et al.2003; Liu et al.2001, 2016; Barrett et al.2012; Murray et al.2014).

Figure 13Zonal mean fraction of 7Be that was produced in the stratosphere in GEOS-Chem driven by MERRA-2, E2.1 nudged to MERRA-2 and E2.1. Zonal mean tropopause height in each simulation is shown as a dashed line.


Figure 13 shows the annual zonal fraction of 7Be of stratospheric origin in each simulation. Using E2.1 (MERRA-2) meteorology, we find that 23 % (30 %) of annual average surface 7Be abundance from 38–51 N is of stratospheric origin, slightly lower (higher) than the observational constraint of 25 % reported by Dutkiewicz and Husain (1985). This is a dramatic improvement over the earlier GCAP/ICECAP studies in which the stratospheric downwelling source in the E2.1-driven simulations was greatly overestimated (Murray et al.2014) and was also seen in the NASA Global Modeling Initiative (GMI) CTM driven by GISS Model II' meteorology (Liu et al.2016). This possibly reflects improvements in downward mass flux in the GCM associated with the increase in vertical resolution (Fig. 1) and updates to the GEOS-Chem dynamical core code that occurred since the original GCAP/ICECAP.

6.1.5 Deposition

Figure 14 (top row) compares the simulated wet deposition flux of 210Pb in each simulation to the observational values reported by Preiss et al. (1996) aggregated to model resolution. Wet deposition of 210Pb in all simulations is biased low by  20 % in all simulations, consistent with earlier versions of GEOS-Chem (Murray et al.2014). These results imply a low bias in the 210Pb source (and therefore 222Rn emission) or a high bias in the dry deposition flux in all simulations.

Figure 14Radionuclide wet deposition fluxes. Annual average wet deposition flux in mBq m−2 s−1 for 210Pb (top row) and 7Be (bottom row). The left column shows observed values aggregated to the model resolution for 210Pb and 7Be. The three columns on the right show (from left to right) values simulated by GEOS-Chem driven by meteorology archived for 2005–2014 CE from MERRA-2, E2.1 nudged to MERRA-2, and free-running E2.1. The number in each panel's lower left shows the mean value of the observations or models sampled at the observed locations. The top right number shows the pattern correlation (R) of the simulated values with the respective observations. The lower right number shows the mean bias of the simulated values with respect to the observations.

Figure 14 (bottom row) compares the simulated wet deposition fluxes of 7Be to the few observations that exist (Baskaran et al.1993; Bleichrodt1978; Brown et al.1989; Dibb1989; Du et al.2008; Harvey and Matthews1989; Hasebe et al.1981; Hirose et al.2004; Igarashi et al.1998; Narazaki and Fujitaka2010; Nijampurkar and Rao1993; Olsen et al.1985; Papastefanou et al.1995; Schuler et al.1991; Turekian et al.1983; Wallbrink and Murray1994, and references therein). The wet deposition flux is increased in the E2.1 simulation with respect to the MERRA-2 simulation, particularly over the midlatitude oceans, and is the closest simulation to matching the observations.

We find a tropospheric residence time for 210Pb-containing aerosols against deposition of 8.3 d in MERRA-2 versus 8.8 d in the E2.1 simulations (Table 4), an improvement in consistency with respect to GCAP/ICECAP, and within the range of previous estimates of 6.5–12.5 d (Turekian et al.1977; Lambert et al.1982; Balkanski et al.1993; Koch et al.1996; Guelle et al.1998b, a; Liu et al.2001). We find a similar relative increase in the lifetime of 7Be-containing aerosols of 19.7 d versus 21.7 d, also consistent with earlier findings of 23 d (Koch et al.1996) and 21 d (Liu et al.2001). Whereas wet deposition in MERRA-2-driven simulations is primarily due to stratiform clouds, it is primarily due to convective clouds in E2.1-driven simulations, consistent with the meteorological fields (see Sect. 4) and leading to the longer residence times in the E2.1 simulations.

6.2 Chemistry

Next, we evaluate the “standard” full-chemistry simulation of version 12.9.3 of GEOS-Chem driven by the three meteorology products for 2005–2014 CE. The standard mechanism contains a unified chemical mechanism of ozone–NOx–hydrocarbon–halogen–aerosol chemistry from the surface to the mesopause (Bey et al.2001; Park et al.2004; Eastham et al.2014; Wang et al.2019). As of version 12.9.3, this includes 255 interactive species, 597 gas-phase reactions, 101 heterogeneous reactions, and 153 photolysis reactions. All our simulations use identical prescribed emissions from the CMIP6 experiments as described in Sect. 5.1 (as opposed to the default GEOS-Chem anthropogenic inventories) and natural emissions vary with each model meteorology as described in Sect. 5.2. Volcanic emissions of SO2 from 2005 were used in all simulation years.

6.2.1 Hydroxyl radical

Table 5 assesses hydroxyl radical (OH) concentrations in the model by comparing the simulated lifetime of relatively long-lived molecules whose main sink is tropospheric OH with observational constraints. It is common for atmospheric models to be biased high with respect to OH in these observational constraints (i.e., low with respect to lifetime, e.g., Naik et al.2013; Voulgarakis et al.2013), which is true as well in our three simulations. However, the lifetime of methyl chloroform (CH3CCl3) and methane in the E2.1 simulations is statistically consistent with the low end of the observed constraint of 6.0-0.4+0.5 years and 10.2-0.7+0.9 years from Prinn et al. (2005), respectively. All simulations fall within the multi-model estimates of a tropospheric methane lifetime of 10.2 ± 1.7 years (Fiore et al.2009) and 9.8 ± 1.6 years (Voulgarakis et al.2013) but fall short of the observationally derived estimates of 11.2 ± 1.3 years from Prather et al. (2012). In all simulations, the E2.1 simulations better match the observational constraints. The E2.1-driven simulations in the tropics have thinner overhead ozone columns (see Sect. 6.2.4), greater free-tropospheric water vapor abundances, and greater lightning NO emissions (see Fig. 8), all of which promote increased OH (e.g., Murray et al.2014), so their lower OH may reflect their higher biogenic NMHC emissions (Fig. 8).

Table 5Lifetime against oxidation by tropospheric OH (years).

a Prinn et al. (2005) for CH3CCl3 and Prather et al. (2012) for CH4. b 2005–2014 CE annual mean and standard deviation.

Download Print Version | Download XLSX

6.2.2 Oxidized nitrogen

Table 6 gives the tropospheric budget for total reactive nitrogen (NOy) in all three simulations, which includes NOx ( NO + NO2) and its longer-lived reservoir species such as nitric acid (HONO2) and peroxyacetyl nitrate (CH3C(O)OONO2; PAN). The E2.1 simulations have greater total sources of NOy because of the larger natural sources from lightning and soils (see Sect. 5.2) and a greater flux of NOy transported from the stratosphere from the products of nitrous oxide (N2O) oxidation. The NOy speciation between family members is largely consistent between the simulations. The lifetime of NOy in the nudged E2.1 simulation is longer than in the free-running E2.1 simulation, reflecting the reduction in that simulation's global mean precipitation flux.

Table 6Tropospheric total reactive nitrogen (NOy)a family budget for 2005–2014 CE.


Download Print Version | Download XLSX

Figure 15Annual average tropospheric columns of nitrogen dioxide (in 1016 molec cm−2; top row), total columns of formaldehyde (in 1016 molec cm−2; middle row), and carbon monoxide mixing ratio at 500 hPa (in ppbv; bottom row) for 2005–2014 CE. The left column from top to bottom shows the respective observations from OMI (Krotkov et al.2017; González Abad et al.2015) and AIRS (Tian et al.2020). The three columns on the right show equivalent values determined from GEOS-Chem driven by MERRA-2, E2.1 nudged to MERRA-2, and free-running E2.1 meteorology. Gray dots show locations where the simulated values are statistically different from the observations with respect to interannual variability (p value < 0.05; n= 10 years). The number in the lower left of each panel shows the global mean value. The number in each model panel's top right shows the pattern correlation (R) between the simulated and observed values. The number in the lower right shows the mean bias of the model with respect to the observations. The tropopause was determined in the simulations using the thermal lapse rate for comparison with satellite tropospheric products.

Panels a–d of Fig. 15 compare the spatial distribution of tropospheric columns of nitrogen dioxide (NO2) for 2005–2014 CE from version 3 of the OMNO2d product from the Ozone Monitoring Instrument (OMI) on the Aura satellite (, Krotkov et al.2019, 2017) to the three simulations. Figure S61 in the Supplement provides seasonal details. The simulations have been sampled at the satellite's overpass time and the tropopause was determined online within the model following the thermal definition is used to calculate the partial columns. All simulations well reproduce the spatial distribution of tropospheric NO2 (all R≥0.94) and have small global mean low biases but statistically disagree with the satellite product in the subtropical latitudes. The subtropical disagreement is potentially due to uncertainties in the tropopause height between the various products. Tropospheric columns match over East Asia during this period but are underestimated over North America and Europe.

6.2.3 Reduced carbon

Table 7 gives the tropospheric emissions and lifetimes for key reduced carbon species in all three simulations. Species primarily emitted from terrestrial plants such as isoprene and monoterpenes have higher emission rates in E2.1 due to the higher diffuse radiation fluxes than in MERRA-2 (see Sects. 4 and 5.2). Species that are primarily lost via oxidation by OH, such as isoprene, have shorter lifetimes in E2.1 due to the higher OH abundances (see Table 5). In contrast, soluble species such as methanol have longer atmospheric lifetimes in E2.1 due to the lower large-scale stratiform precipitation rates in E2.1 relative to MERRA-2.

Table 7Biogenic emissions and tropospheric lifetimes of select NMHC species for 2005–2014 CE.

Annual mean and standard deviation. The percentage of total emission is given per biogenic source.

Download Print Version | Download XLSX

Panels e–h of Fig. 15 compare total columns of formaldehyde (HCHO) from version 3 of the OMHCHOd product from OMI on the Aura satellite (, Chance2019; González Abad et al.2015) to the three simulations for 2005–2014 CE. Figure S62 in the Supplement provides seasonal details. The simulations have been sampled at the overpass time of the satellite. Formaldehyde is a common product of hydrocarbon oxidation. Despite higher mean values over the continents, all three simulations are statistically consistent with respect to the large amount of interannual variability in the simulated and observed HCHO columns. There is strong spatial correlation between the simulations and the satellite product (all R≥0.8). Terrestrial columns simulated by the E2.1 simulations are higher than the MERRA-2 simulations, reflecting the higher biogenic emissions. All simulations underestimate HCHO columns over the remote ocean aside from the continental outflows of North America and Asia and over the Arctic.

Table 8Tropospheric carbon monoxide (CO) budget for 2005–2014 CE.

Annual mean and standard deviation. The percentage of the total is given per source and sink.

Download Print Version | Download XLSX

Table 8 gives the tropospheric budget for carbon monoxide (CO) in all three simulations. The direct emissions of CO from anthropogenic and biomass burning sources between the three simulations are identical by experimental design. However, the chemical source of CO from methane and non-methane hydrocarbon oxidation is higher in the E2.1 simulations due to the higher OH abundances (Table 5). This is balanced by the increased chemical loss of CO by OH. The influence of OH on CO production from short-lived non-methane hydrocarbon species seems to outweigh the influence on CO loss, and consequently there are slightly longer CO tropospheric lifetimes in the E2.1 simulations.

Panels i–l of Fig. 15 compare CO mixing ratios at 500 hPa from version 7 of the AIRX3STD product from the Atmospheric Infrared Sounder (AIRS) on the Aqua satellite (, AIRS project2019; Tian et al.2020) to the three simulations for 2005–2014 CE. Figure S63 in the Supplement provides seasonal details. The simulations have been sampled at the overpass time of the satellite. All simulations have strong spatial correlation with the observations (R=0.93). CO in the free troposphere is higher in the E2.1 simulations, consistent with its longer lifetime (Table 8). The E2.1 simulations are globally biased high by 15 % compared to the AIRS values. However, they are statistically consistent almost everywhere with respect to the large amount of interannual variability in the observations and simulations. In contrast, the MERRA-2 simulations are globally biased low by 5 % and significantly so throughout the tropics and Southern Hemisphere. All models underestimate the curious enhancement of CO seen in the AIRS climatology over Antarctica.

6.2.4 Ozone

Table 9 gives the tropospheric budget for the odd-oxygen family (Ox O3+ rapid cycling species) for all three simulations. The individual budget terms are all are consistent with the range of reported values from the Tropospheric Ozone Assessment Report (TOAR) multi-model assessment (see Fig. 3 of Young et al.2018), the CMIP6 models that performed interactive tropospheric chemistry (Griffiths et al.2021), as well as the last extensive tropospheric ozone budget evaluation within the standard GEOS-Chem model (Hu et al.2017). The E2.1 simulations are on the high end of the previously reported values due to the higher tropopause height in those simulations; the upper troposphere and lower stratosphere regions contribute strongly to each Ox budget term due to the rapidly increasing abundances of ozone with altitude there. We point out that the stratosphere-to-troposphere flux calculated using the “residual method” of the other budget terms yields consistent results when we track the mass flux of ozone across the dynamic tropopause in the model. Relative to the earlier GCAP and ICECAP simulations (Murray et al.2014), the transport of ozone from the stratosphere is dramatically improved.

Table 9Tropospheric odd-oxygen (Ox) family budget for 2005–2014 CE.

A molar mass of 48 g is assumed for Ox.

Download Print Version | Download XLSX

Figure 16Comparison of the annual cycle of ozone for 2005–2014 CE between ozonesonde observations (black circles) and the MERRA-2 (solid orange line), E2.1 nudged to MERRA-2 (blue line), and E2.1 (red line) simulations. Model and observational data were grouped into four latitude bands (90 to 30 S, 30 S to 0, 0 S to 30 N, and 30 S to 90 N) and sampled at three altitudes (700, 500, and 250 hPa), with the models sampled at locations and months of the ozonesonde measurements before averaging together. Error bars on the observations indicate the average interannual standard deviation for each group of observations. The correlation (r) and mean normalized bias error (mnbe) for the MERRA-2 (orange), E2.1 nudged to MERRA-2 (blue), and E2.1 (red) means versus the observations are also indicated in each panel.


Figure 16 evaluates the zonal and seasonal distribution of ozone versus in situ measurements. We use the ozonesonde measurements archived by the World Ozone and Ultraviolet Radiation Data Centre (WOUDC) of the World Meteorological Organization/Global Atmosphere Watch Program (WMO/GAW). The data were accessed on 4 November 2019 from (WMO/GAW2019). All models fall within the variability of the measurements, except the free troposphere of the northern extratropics, where the simulations are biased low, especially during the summer months. The MERRA-2 simulations better reproduce the seasonality of ozone in the tropical free and upper troposphere, likely reflecting the climatologically constrained lightning NO source in that version (see Sect. 5.2).

Figure 17Annual average surface ozone mixing ratio in ppbv ( nmol mol−1; top row), tropospheric columns of ozone (TCOs; in Dobson units; middle row), and total ozone columns (TOC; in Dobson units; bottom row) for 2005–2014 CE. The left column from top to bottom shows observations from TOAR (Schultz et al.2017a), OMI/MLS (Ziemke et al.2006) and OMI (Dobber et al.2006), respectively. The three columns on the right show equivalent values determined from GEOS-Chem driven by MERRA-2, E2.1 nudged to MERRA-2, and free-running E2.1 meteorology. Gray dots show locations where the simulated values are statistically different from the observations with respect to interannual variability (p value < 0.05; n= 10 years). The number in the lower left of each panel shows the global mean value. The number in each model panel's top right shows the pattern correlation (R) between the simulated and observed values. The number in the lower right shows the mean bias of the model with respect to the observations. The tropopause was determined in the simulations using the thermal lapse rate for comparison with satellite tropospheric products.

Figure 17 evaluates spatial distributions of ozone in the three simulations against surface in situ and satellite observations. The top row shows the annual average surface ozone mixing ratio in ppbv ( nmol mol−1) from our simulations versus the gridded mean non-urban surface ozone product from the Tropospheric Ozone Assessment Report (TOAR) (, Schultz et al.2017b, a). The middle row shows tropospheric columns of ozone (TCO; in Dobson units) versus the joint Ozone Monitoring Instrument (OMI) and Microwave Limb Sounder (MLS) product from the Aura satellite (Ziemke et al.2006). The bottom row shows total ozone columns (TOCs; in Dobson units) versus the OMDOAO3e product from OMI (, Veefkind2012; Dobber et al.2006). All observational products have been aggregated to model resolution for comparison. In the case of the satellite products, we sampled the model at the overpasses' time and location. We use the online thermal lapse rate tropopause to calculate TCO. Seasonal details are given in Figs. S64–S66 of the Supplement.

All simulations are biased high by 15 %–17 % with respect to surface ozone, mostly driven by the eastern North American data. The E2.1 simulations are especially higher over the Amazon than either MERRA-2 or the observations. Comparisons of TCO to OMI/MLS are sensitive to uncertainties in the tropopause location in the satellite product versus the models (Griffiths et al.2021). The E2.1 TCOs are globally higher than the OMI/MLS product by 12 %, reflecting the lower tropopause pressures. However, the only locations in which it is statistically different regarding interannual variability are over the tropical oceans, where it is biased low. This likely reflects the more vigorous convection in E2.1 that leads to ozone destruction (e.g., Murray et al.2013). Tropospheric columns in MERRA-2 match the global mean from OMI/MLS but also underestimate the western Pacific and additionally underestimate northern extratropical ozone. All simulations underestimate tropospheric columns in the southern extratropics with respect to OMI/MLS. The models all show excellent agreement with respect to total ozone columns with small positive mean global biases of 4 % and high pattern correlation (all R≥0.96). Total ozone in the tropics is higher in MERRA-2 than E2.1, consistent with differences in the rate of vertical ascent in the tropical pipe implied by the stratospheric age of air comparison (see Sect. 6.1.3). The E2.1 simulation overestimates Antarctic ozone relative to MERRA-2 or the nudged simulation, although not significantly compared to interannual variability.

6.2.5 Particulate matter

Figure 18 evaluates the spatial distribution of particulate matter in the three simulations against satellite observations. The top row shows the simulated concentration of fine particulate matter under 2.5 µm (PM2.5) in µg m−3 versus the historical reconstruction from Hammer et al. (2020). The middle row shows the total aerosol optical thickness (AOT) at 550 nm (unitless) in the simulations versus the combined Dark Target and Deep Blue AOT at 0.55 µm from Collection 6.1 of the MODerate resolution Imaging Spectroradiometer (MODIS) MYD08 product from the Aqua satellite (, Platnick et al.2015). The bottom row shows the total column of sulfur dioxide (SO2) in Dobson units in the simulations versus the second public release of version 3 of the OMSO2e product from OMI on the Aura satellite (, Li et al.2020b, a). The model has been sampled at the time of the Aura and Aqua overpasses for comparison to the satellite products. Figures S67–S68 of the Supplement provide seasonal details for AOT and the SO2 columns.

Figure 18Annual average surface concentration of fine particulate matter (PM2.5; in µg m−3; top row), aerosol optical thickness at 550 nm (unitless; middle row), and total column of sulfur dioxide (in Dobson units; bottom row) for 2005–2014 CE. The left column from top to bottom shows observations from Hammer et al. (2020), Aqua MODIS, and OMI (Li et al.2020a), respectively. The three columns on the right show equivalent values determined from GEOS-Chem driven by MERRA-2, E2.1 nudged to MERRA-2, and free-running E2.1 meteorology. Gray dots show locations where the simulated values are statistically different from the observations with respect to interannual variability (p value < 0.05; n= 10 years). The number in the lower left of each panel shows the global mean value. The number in each model panel's top right shows the pattern correlation (R) between the simulated and observed values. The number in the lower right shows the mean bias of the model with respect to the observations.

Surface PM2.5 agrees between the simulations and Hammer et al. (2020) over heavily industrialized regions. However, we note that the Hammer et al. (2020) product used ratios of surface PM2.5 to AOT from GEOS-Chem to generate their proxy reconstruction from satellite AOT measurements, so it is not an entirely independent comparison. Surface concentrations are underestimated almost everywhere else in all simulations, especially in regions heavily influenced by mineral dust and biomass burning. The same story is reflected in the MODIS AOT, except the stronger sea-salt aerosol emission in the E2.1 simulations yielding better comparison in the Southern Ocean and portions of the Pacific. Simulated total columns of sulfur dioxide are underestimated everywhere in all simulations except over East Asia and the locations of the major volcanic eruptions of 2005 (Sierra Negra in the Galápagos and Anatahan of the Northern Mariana Islands).

7 Summary

This paper described and evaluated version 2.0 of the Global Change and Air Pollution (GCAP) chemical-transport model framework.

GCAP 2.0 is a one-way offline coupling between the E2.1 version of the NASA GISS GCM frozen for the CMIP6 experiments (Kelley et al.2020; Miller et al.2021) and the GEOS-Chem 3-D chemical-transport model (, last access: 9 September 2021; Bey et al.2001). Additional subdaily diagnostics were added to E2.1 to archive the same fields as the MERRA-2 reanalysis product (Gelaro et al.2017) that is normally used to drive GEOS-Chem. We then re-performed one of the atmosphere-only members of the E2.1 contributions to the CMIP6 ensemble, archiving the meteorology necessary for driving GEOS-Chem. The E2.1 meteorology is available at 2 latitude by 2.5 longitude and for 40 vertical layers ranging from the surface to 0.1 hPa. At publication time, meteorology is available for the pre-industrial era (1851–1860 CE) and recent past (2001–2014 CE), including a recent-past simulation nudged to MERRA-2 to assist users in comparing with observations. Also available is meteorology for the near future (2040–2049 CE) and the end of the century (2090–2099 CE) for seven future SSP scenarios ranging from extreme mitigation (SSP1-1.9) to extreme warming (SSP5-8.5). In addition, the CMIP6 emissions and surface boundary conditions (Hoesly et al.2018; van Marle et al.2017; Riahi et al.2017; Meinshausen et al.2017; Gidden et al.2019) have been processed for input into GEOS-Chem. GCAP 2.0 is operational in all current variants of the GEOS-Chem model, with all GCClassic run directories and input files provided. All GCAP 2.0 input data are publicly served at (last access: 9 September 2021).

The meteorology was evaluated by comparing to both the original simulation and the MERRA-2 reanalysis for the recent past. Surface air in the repeat simulation is slightly warmer than the original run due to increased calls to the radiation code necessary for archiving shortwave fluxes for input to GEOS-Chem. The E2.1 climatology in the recent past largely agrees with the MERRA-2 climatology for that period, with the primary difference being in the relative amount of precipitation in stratiform versus convective clouds as well as a higher tropopause height in E2.1. Emissions that respond to meteorology in GEOS-Chem are slightly higher in the E2.1-driven simulations, including biogenic emissions from terrestrial plants, the lightning and soil microbial sources of reactive nitrogen, and sea-salt evasion. The dust mobilization parameterization was found to be extremely sensitive to resolution and meteorology, and scaling factors have been determined to constrain the global source.

Model physics and transport were evaluated using simulations and observations of sulfur hexafluoride and radionuclides. In all cases, transport is substantially improved over the original GCAP, and the E2.1-driven simulations perform comparably to the MERRA-2-driven simulations. Most importantly, whereas age of air remains too young in both E2.1 and MERRA-2, the stratosphere-to-troposphere mass flux now yields comparable values, with consistent stratosphere-to-troposphere fluxes of ozone with multi-model means. This is a major improvement over the previous versions of GCAP (Wu et al.2007; Murray et al.2014). However, we urge users to be cautious when using these fields for stratospheric chemistry–climate applications. Future simulations will provide CMIP6 meteorology from the 102-layer version of the GISS GCM (E2.2, Rind et al.2020), which will be better suited for studies of the middle atmosphere; E2.2 includes an interactive quasi-biennial oscillation and improved polar vortex variability including sudden warmings, which could contribute additional dynamical variability that GEOS-Chem may otherwise not see (Orbe et al.2020b).

Lastly, the chemistry of the model using the CMIP6 emissions and the different meteorology was evaluated against a suite of satellite products and in situ observations. The E2.1-driven simulations have lower OH and therefore more accurate methyl chloroform and methane lifetimes. Greater biogenic fluxes yield greater abundances of oxidation products (e.g., CO) in the E2.1-driven simulations, improving comparison with observations in the Southern Hemisphere. However, the MERRA-2-driven simulations have a superior representation of free-tropospheric ozone, likely due to the constrained lightning seasonality and distribution as well as the more realistic tropopause pressure in those simulations. All simulations underestimate particulate matter abundances outside of industrialized areas and AOT in most places. Otherwise, model performance was very similar in all simulations.

Code and data availability

GCAP 2.0 entered the public GEOS-Chem repository in version 13.1.0 (, The International GEOS-Chem User Community2021a; ​​, The International GEOS-Chem User Community2021b;, Lundgren et al.2021). Source code for generating E2.1 output to drive GEOS-Chem is available at (Murray et al.2021). Archived meteorology, emissions, and boundary conditions for GCAP 2.0 are hosted online at (last access: 9 September 2021, Murray2021). LTM can generate additional CMIP6 scenarios and time periods upon request.


The supplement related to this article is available online at:

Author contributions

LTM conceived the project, performed the model development, simulation, and evaluation of the E2.1 and GEOS-Chem interface, and wrote the manuscript. EML and LJM provided code edits from ModelE that were updated to E2.1 by LTM. CO provided feedback on the implementation of the new subdaily diagnostics. MS was a developer of the FlexGrid code in GEOS-Chem. All authors contributed feedback to the manuscript.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We thank the entire GISS and GEOS-Chem developer and user communities.

Financial support

This research has been supported by the National Science Foundation (grant nos. AGS-1702106 and AGS-2002414 to Lee T. Murray) and the National Aeronautics and Space Administration, Earth Sciences Division (grant no. NNX13AO08G to Loretta J. Mickley).

Review statement

This paper was edited by Fiona O'Connor and reviewed by two anonymous referees.


Achakulwisut, P., Mickley, L. J., Murray, L. T., Tai, A. P. K., Kaplan, J. O., and Alexander, B.: Uncertainties in isoprene photochemistry and emissions: implications for the oxidative capacity of past and present atmospheres and for climate forcing agents, Atmos. Chem. Phys., 15, 7977–7998,, 2015. a

AIRS project: Aqua/AIRS L3 Daily Standard Physical Retrieval (AIRS+AMSU) 1 degree x 1 degree V7.0, Greenbelt, MD, USA, Goddard Earth Sciences Data and Information Services Center (GES DISC),, 2019. a

Allen, D. J., Rood, R., Thompson, A. M., and Hudson, R.: Three-dimensional radon 222 calculations using assimilated meteorological data and a convective mixing algorithm, J. Geophys. Res.-Atmos., 101, 6871–6881,, 1996. a

Allen, D. J., Dibb, J., Ridley, B., Pickering, K., and Talbot, R.: An estimate of the stratospheric contribution to springtime tropospheric ozone maxima using TOPSE measurements and beryllium-7 simulations, J. Geophys. Res.-Atmos., 108, 8355,, 2003. a, b

Balkanski, Y., Jacob, D. J., Gardner, G., Graustein, W., and Turekian, K.: Transport and Residence Times of Tropospheric Aerosols Inferred from a Global 3-Dimensional Simulation of Pb-210, J. Geophys. Res.-Atmos., 98, 20573–20586,, 1993. a

Barrett, S. R. H., Yim, S. H. L., Gilmore, C. K., Murray, L. T., Kuhn, S. R., Tai, A. P. K., Yantosca, R. M., Byun, D. W., Ngan, F., Li, X., Levy, J. I., Ashok, A., Koo, J., Wong, H. M., Dessens, O., Balasubramanian, S., Fleming, G. G., Pearlson, M. N., Wollersheim, C., Malina, R., Arunachalam, S., Binkowski, F. S., Leibensperger, E. M., Jacob, D. J., Hileman, J. I., and Waitz, I. A.: Public Health, Climate, and Economic Impacts of Desulfurizing Jet Fuel, Environ. Sci. Technol., 46, 4275–4282,, 2012. a, b

Baskaran, M., Coleman, C., and Santschi, P.: Atmospheric Depositional Fluxes of Be-7 and Pb-210 at Galveston and College-Station, Texas, J. Geophys. Res.-Atmos., 98, 20555–20571,, 1993. a

Bauer, N., Calvin, K., Emmerling, J., Fricko, O., Fujimori, S., Hilaire, J., Eom, J., Krey, V., Kriegler, E., Mouratiadou, I., Sytze de Boer, H., van den Berg, M., Carrara, S., Daioglou, V., Drouet, L., Edmonds, J. E., Gernaat, D., Havlik, P., Johnson, N., Klein, D., Kyle, P., Marangoni, G., Masui, T., Pietzcker, R. C., Strubegger, M., Wise, M., Riahi, K., and van Vuuren, D. P.: Shared Socio-Economic Pathways of the Energy Sector – Quantifying the Narratives, Glob. Environ. Change, 42, 316–330,, 2017. a

Bauer, S. E., Tsigaridis, K., Faluvegi, G., Kelley, M., Lo, K. K., Miller, R. L., Nazarenko, L., Schmidt, G. A., and Wu, J.: Historical (1850–2014) Aerosol Evolution and Role on Climate Forcing Using the GISS ModelE2.1 Contribution to CMIP6, J. Adv. Model. Earth Syst., 12, e2019MS001978,, 2020. a

Bey, I., Jacob, D. J., Yantosca, R. M., Logan, J. A., Field, B. D., Fiore, A. M., Li, Q., Liu, H. Y., Mickley, L. J., and Schultz, M. G.: Global modeling of tropospheric chemistry with assimilated meteorology: Model description and evaluation, J. Geophys. Res.-Atmos., 106, 23073–23095,, 2001. a, b, c, d

Bindle, L., Martin, R. V., Cooper, M. J., Lundgren, E. W., Eastham, S. D., Auer, B. M., Clune, T. L., Weng, H., Lin, J., Murray, L. T., Meng, J., Keller, C. A., Pawson, S., and Jacob, D. J.: Grid-Stretching Capability for the GEOS-Chem 13.0.0 Atmospheric Chemistry Model, Geosci. Model Dev. Discuss. [preprint],, in review, 2020. a

Bleichrodt, J.: Mean Tropospheric Residence Time of Cosmic-Ray-Produced Beryllium-7 at North Temperate Latitudes, J. Geophys. Res.-Oc. Atm., 83, 3058–3062,, 1978. a

Bondietti, E., Brantley, J., and Rangarajan, C.: Size Distributions and Growth of Natural and Chernobyl-Derived Sub-Micron Aerosols in Tennessee, J. Environ. Radioactiv., 6, 99–120,, 1988. a, b

Bradley, W. and Pearson, J.: Aircraft Measurements of Vertical Distribution of Radon in Lower Atmosphere, J. Geophys. Res., 75, 5890,, 1970. a

Brost, R. and Chatfield, R.: Transport of Radon in a 3-Dimensional, Subhemispheric Model, J. Geophys. Res.-Atmos., 94, 5095–5119,, 1989. a

Brost, R., Feichter, J., and Heimann, M.: 3-Dimensional Simulation of Be-7 in a Global Climate Model, J. Geophys. Res.-Atmos., 96, 22423–22445,, 1991. a

Brown, L., Stensland, G., Klein, J., and Middleton, R.: Atmospheric Deposition of Be-7 and Be-10, Geochim. Cosmochim. Ac., 53, 135–142,, 1989. a

Calvin, K., Bond-Lamberty, B., Clarke, L., Edmonds, J., Eom, J., Hartin, C., Kim, S., Kyle, P., Link, R., Moss, R., McJeon, H., Patel, P., Smith, S., Waldhoff, S., and Wise, M.: The SSP4: A world of deepening inequality, Glob. Environ. Change, 42, 284–296,, 2017. a

Carn, S. A., Yang, K., Prata, A. J., and Krotkov, N. A.: Extending the long-term record of volcanic SO2 emissions with the Ozone Mapping and Profiler Suite nadir mapper, Geophys. Res. Lett., 42, 925–932,, 2015. a

Cecil, D. J., Buechler, D. E., and Blakeslee, R. J.: Gridded lightning climatology from TRMM-LIS and OTD: Dataset description, Atmos. Res., 135–136, 404–414,, 2014. a

Chance, K.: OMI/Aura Formaldehyde (HCHO) Total Column Daily L3 Weighted Mean Global 0.1deg Lat/Lon Grid V003, Greenbelt, MD, USA, Goddard Earth Sciences Data and Information Services Center (GES DISC),, 2019. a

Considine, D. B., Bergmann, D. J., and Liu, H.: Sensitivity of Global Modeling Initiative chemistry and transport model simulations of radon-222 and lead-210 to input meteorological data, Atmos. Chem. Phys., 5, 3389–3406,, 2005. a, b

Crespo Cuaresma, J.: Income projections for climate change research: A framework based on human capital dynamics, Glob. Environ. Change, 42, 226–236,, 2017. a

Del Genio, A. D., Chen, Y., Kim, D., and Yao, M.-S.: The MJO Transition from Shallow to Deep Convection in CloudSat/CALIPSO Data and GISS GCM Simulations, J. Climate, 25, 3755–3770,, 2012. a

Del Genio, A. D., Wu, J., Wolf, A. B., Chen, Y., Yao, M.-S., and Kim, D.: Constraints on Cumulus Parameterization from Simulations of Observed MJO Events, J. Climate, 28, 6419–6442,, 2015. a

Dellink, R., Chateau, J., Lanzi, E., and Magné, B.: Long-term economic growth projections in the Shared Socioeconomic Pathways, Glob. Environ. Change, 42, 200–214,, 2017. a

Dibb, J.: Atmospheric Deposition of Beryllium-7 in the Chesapeake Bay Region, J. Geophys. Res.-Atmos., 94, 2261–2265,, 1989. a

Dibb, J., Talbot, R., and Gregory, G.: Beryllium-7 and Pb-210 in the Western-Hemisphere Arctic Atmosphere - Observations from Three Recent Aircraft-Based Sampling Programs, J. Geophys. Res.-Atmos., 97, 16709–16715,, 1992. a

Dibb, J., Meeker, L., Finkel, R., Southon, J., Caffee, M., and Barrie, L.: Estimation of Stratospheric Input to the Arctic Troposphere - Be-7 and Be-10 in Aerosols at Alert, Canada, J. Geophys. Res.-Atmos., 99, 12855–12864,, 1994. a

Dobber, M., Dirksen, R., Levelt, P., van den Oord, G., Voors, R., Kleipool, Q., Jaross, G., Kowalewski, M., Hilsenrath, E., Leppelmeier, G., Johan, d. V., Dierssen, W., and Rozemeijer, N.: Ozone monitoring instrument calibration, IEEE T. Geosci. Remote, 44, 1209–1238,, 2006. a, b

Du, J., Zhang, J., Zhang, J., and Wu, Y.: Deposition patterns of atmospheric Be-7 and Pb-210 in coast of East China Sea, Shanghai, China, Atmos. Environ., 42, 5101–5109,, 2008. a

Dutkiewicz, V. and Husain, L.: Stratospheric and Tropospheric Components of Be-7 in Surface Air, J. Geophys. Res.-Atmos., 90, 5783–5788,, 1985. a

Eastham, S., Doherty, S., Keith, D., Richter, J. H., and Xia, L.: Air quality and climate connections., Eos, 102,, 2021. a

Eastham, S. D., Weisenstein, D. K., and Barrett, S. R.: Development and evaluation of the unified tropospheric–stratospheric chemistry extension (UCX) for the global chemistry-transport model GEOS-Chem, Atmos. Environ., 89, 52–63,, 2014. a, b, c

Eastham, S. D., Long, M. S., Keller, C. A., Lundgren, E., Yantosca, R. M., Zhuang, J., Li, C., Lee, C. J., Yannetti, M., Auer, B. M., Clune, T. L., Kouatchou, J., Putman, W. M., Thompson, M. A., Trayanov, A. L., Molod, A. M., Martin, R. V., and Jacob, D. J.: GEOS-Chem High Performance (GCHP v11-02c): a next-generation implementation of the GEOS-Chem chemical transport model for massively parallel applications, Geosci. Model Dev., 11, 2941–2953,, 2018. a, b, c

Eyring, V., Bony, S., Meehl, G. A., Senior, C. A., Stevens, B., Stouffer, R. J., and Taylor, K. E.: Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization, Geosci. Model Dev., 9, 1937–1958,, 2016. a

Feichter, J. and Crutzen, P.: Parameterization of vertical tracer transport due to deep cumulus convection in a global transport model and its evaluation with 222Radon measurements, Tellus B, 42, 100–117,, 1990. a

Finney, D. L., Doherty, R. M., Wild, O., Stevenson, D. S., MacKenzie, I. A., and Blyth, A. M.: A projected decrease in lightning under climate change, Nat. Clim. Change, 8, 210–213,, 2018. a

Fiore, A., Naik, V., and Leibensperger, E.: Air quality and climate connections, J. Air Waste Manag. Assoc., 65, 645–685,, 2015. a

Fiore, A. M., Dentener, F. J., Wild, O., Cuvelier, C., Schultz, M. G., Hess, P., Textor, C., Schulz, M., Doherty, R. M., Horowitz, L. W., MacKenzie, I. A., Sanderson, M. G., Shindell, D. T., Stevenson, D. S., Szopa, S., Van Dingenen, R., Zeng, G., Atherton, C., Bergmann, D., Bey, I., Carmichael, G., Collins, W. J., Duncan, B. N., Faluvegi, G., Folberth, G., Gauss, M., Gong, S., Hauglustaine, D., Holloway, T., Isaksen, I. S. A., Jacob, D. J., Jonson, J. E., Kaminski, J. W., Keating, T. J., Lupu, A., Marmer, E., Montanaro, V., Park, R. J., Pitari, G., Pringle, K. J., Pyle, J. A., Schroeder, S., Vivanco, M. G., Wind, P., Wojcik, G., Wu, S., and Zuber, A.: Multimodel estimates of intercontinental source-receptor relationships for ozone pollution, J. Geophys. Res.-Atmos., 114, D04301,, 2009. a

Fountoukis, C. and Nenes, A.: ISORROPIA II: a computationally efficient thermodynamic equilibrium model for K+Ca2+Mg2+NH4+Na+SO42-NO3-ClH2O aerosols, Atmos. Chem. Phys., 7, 4639–4659,, 2007. a

Fricko, O., Havlik, P., Rogelj, J., Klimont, Z., Gusti, M., Johnson, N., Kolp, P., Strubegger, M., Valin, H., Amann, M., Ermolieva, T., Forsell, N., Herrero, M., Heyes, C., Kindermann, G., Krey, V., McCollum, D. L., Obersteiner, M., Pachauri, S., Rao, S., Schmid, E., Schoepp, W., and Riahi, K.: The marker quantification of the Shared Socioeconomic Pathway 2: A middle-of-the-road scenario for the 21st century, Glob. Environ. Change, 42, 251–267,, 2017. a

Fujimori, S., Hasegawa, T., Masui, T., Takahashi, K., Herran, D. S., Dai, H., Hijioka, Y., and Kainuma, M.: SSP3: AIM implementation of Shared Socioeconomic Pathways, Glob. Environ. Change, 42, 268–283,, 2017. a

Garcia-Menendez, F., Monier, E., and Selin, N. E.: The role of natural variability in projections of climate change impacts on U.S. ozone pollution, Geophys. Res. Lett., 44, 2911–2921,, 2017. a

GBD 2019 Risk Factor Collaborators: Global burden of 87 risk factors in 204 countries and territories, 1990-2019: a systematic analysis for the Global Burden of Disease Study 2019, Lancet, 396, 1223–1249,, 2020. a

Gelaro, R., McCarty, W., Suárez, M., Todling, R., Molod, A., Takacs, L., Randles, C., Darmenov, A., Bosilovich, M., Reichle, R., Wargan, K., Coy, L., Cullather, R., Draper, C., Akella, S., Buchard, V., Conaty, A., da Silva, A., Gu, W., Kim, G., Koster, R., Lucchesi, R., Merkova, D., Nielsen, J., Partyka, G., Pawson, S., Putman, W., Rienecker, M., Schubert, S., Sienkiewicz, M., and Zhao, B.: The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), J. Climate, 30, 5419–5454,, 2017. a, b, c

Geng, L., Zatko, M. C., Alexander, B., Fudge, T. J., Schauer, A. J., Murray, L. T., and Mickley, L. J.: Effects of postdepositional processing on nitrogen isotopes of nitrate in the Greenland Ice Sheet Project 2 ice core, Geophys. Res. Lett., 42, 5346–5354,, 2015. a

Geng, L., Murray, L. T., Mickley, L. J., Lin, P., Fu, Q., Schauer, A. J., and Alexander, B.: Isotopic evidence of multiple controls on atmospheric oxidants over climate transitions, Nature, 546, 133–136,, 2017. a

Gidden, M. J., Riahi, K., Smith, S. J., Fujimori, S., Luderer, G., Kriegler, E., van Vuuren, D. P., van den Berg, M., Feng, L., Klein, D., Calvin, K., Doelman, J. C., Frank, S., Fricko, O., Harmsen, M., Hasegawa, T., Havlik, P., Hilaire, J., Hoesly, R., Horing, J., Popp, A., Stehfest, E., and Takahashi, K.: Global emissions pathways under different socioeconomic scenarios for use in CMIP6: a dataset of harmonized emissions trajectories through the end of the century, Geosci. Model Dev., 12, 1443–1475,, 2019. a, b, c

González Abad, G., Liu, X., Chance, K., Wang, H., Kurosu, T. P., and Suleiman, R.: Updated Smithsonian Astrophysical Observatory Ozone Monitoring Instrument (SAO OMI) formaldehyde retrieval, Atmos. Meas. Tech., 8, 19–32,, 2015. a, b

Griffiths, P. T., Murray, L. T., Zeng, G., Shin, Y. M., Abraham, N. L., Archibald, A. T., Deushi, M., Emmons, L. K., Galbally, I. E., Hassler, B., Horowitz, L. W., Keeble, J., Liu, J., Moeini, O., Naik, V., O'Connor, F. M., Oshima, N., Tarasick, D., Tilmes, S., Turnock, S. T., Wild, O., Young, P. J., and Zanis, P.: Tropospheric ozone in CMIP6 simulations, Atmos. Chem. Phys., 21, 4187–4218,, 2021. a, b

Guelle, W., Balkanski, Y., Dibb, J., Schulz, M., and Dulac, F.: Wet deposition in a global size-dependent aerosol transport model 2. Influence of the scavenging scheme on Pb-210 vertical profiles, surface concentrations, and deposition, J. Geophys. Res.-Atmos., 103, 28875–28891,, 1998a. a

Guelle, W., Balkanski, Y., Schulz, M., Dulac, F., and Monfray, P.: Wet deposition in a global size-dependent aerosol transport model – 1. Comparison of a 1 year Pb-210 simulation with ground measurements, J. Geophys. Res.-Atmos., 103, 11429–11445,, 1998b. a

Guenther, A. B., Jiang, X., Heald, C. L., Sakulyanontvittaya, T., Duhl, T., Emmons, L. K., and Wang, X.: The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1): an extended and updated framework for modeling biogenic emissions, Geosci. Model Dev., 5, 1471–1492,, 2012. a

Hall, B. D., Dutton, G. S., Mondeel, D. J., Nance, J. D., Rigby, M., Butler, J. H., Moore, F. L., Hurst, D. F., and Elkins, J. W.: Improving measurements of SF6 for the study of atmospheric transport and emissions, Atmos. Meas. Tech., 4, 2441–2451,, 2011. a

Hall, T. M. and Waugh, D. W.: Stratospheric residence time and its relationship to mean age, J. Geophys. Res., 105, 6773,, 2000. a

Hammer, M., van Donkelaar, A., Li, C., Lyapustin, A., Sayer, A., Hsu, N., Levy, R., Garay, M., Kalashnikova, O., Kahn, R., Brauer, M., Apte, J., Henze, D., Zhang, L., Zhang, Q., Ford, B., Pierce, J., and Martin, R.: Global Estimates and Long-Term Trends of Fine Particulate Matter Concentrations (1998–2018), Environ. Sci. Technol., 54, 7879–7890,, 2020. a, b, c, d

Hansen, J., Russell, G., Rind, D., Stone, P., Lacis, A., Lebedeff, S., Ruedy, R., and Travis, L.: Efficient Three-Dimensional Global Models for Climate Studies: Models I and II, Mon. Weather Rev., 111, 609–662,<0609:ETDGMF>2.0.CO;2, 1983. a, b

Harvey, M. and Matthews, K.: Be-7 Deposition in a High-Rainfall Area of New-Zealand, J. Atmos. Chem., 8, 299–306,, 1989. a

Hasebe, N., Doke, T., Kikuchi, J., Takeuchi, Y., And Sugiyama, T.: Observation of Fallout Rates of Atmospheric Be-7 and Na-22 Produced by Cosmic-Rays Concerning Estimation of the Fallout Rate of Atmospheric Al-26, J. Geophys. Res.-Space, 86, 520–524,, 1981. a

Hauglustaine, D., Hourdin, F., Jourdain, L., Filiberti, M., Walters, S., Lamarque, J., and Holland, E.: Interactive chemistry in the Laboratoire de Météorologie Dynamique general circulation model: Description and background tropospheric chemistry evaluation, J. Geophys. Res.-Atmos., 109, D04314,, 2004. a

Henze, D. K., Hakami, A., and Seinfeld, J. H.: Development of the adjoint of GEOS-Chem, Atmos. Chem. Phys., 7, 2413–2433,, 2007. a, b

Hirose, K., Honda, T., Yagishita, S., Igarashi, Y., and Aoyama, M.: Deposition behaviors of Pb-210, Be-7 and thorium isotopes observed in Tsukuba and Nagasaki, Japan, Atmos. Environ., 38, 6601–6608,, 2004. a

Hoesly, R. M., Smith, S. J., Feng, L., Klimont, Z., Janssens-Maenhout, G., Pitkanen, T., Seibert, J. J., Vu, L., Andres, R. J., Bolt, R. M., Bond, T. C., Dawidowski, L., Kholod, N., Kurokawa, J.-I., Li, M., Liu, L., Lu, Z., Moura, M. C. P., O'Rourke, P. R., and Zhang, Q.: Historical (1750–2014) anthropogenic emissions of reactive gases and aerosols from the Community Emissions Data System (CEDS), Geosci. Model Dev., 11, 369–408,, 2018. a, b, c

Holton, J., Haynes, P., Mcintyre, M., Douglass, A., Rood, R., and Pfister, L.: Stratosphere-Troposphere Exchange, Rev. Geophys., 33, 403–439,, 1995. a

Horowitz, L. W., Liang, J., Gardner, G. M., and Jacob, D. J.: Export of reactive nitrogen from North America during summertime: Sensitivity to hydrocarbon chemistry, J. Geophys. Res.-Atmos., 103, 13451–13476,, 1998. a

Hu, L., Jacob, D. J., Liu, X., Zhang, Y., Zhang, L., Kim, P. S., Sulprizio, M. P., and Yantosca, R. M.: Global budget of tropospheric ozone: Evaluating recent model advances with satellite (OMI), aircraft (IAGOS), and ozonesonde observations, Atmos. Environ., 167, 323–334,, 2017. a

Hudman, R. C., Moore, N. E., Mebust, A. K., Martin, R. V., Russell, A. R., Valin, L. C., and Cohen, R. C.: Steps towards a mechanistic model of global soil nitric oxide emissions: implementation and space based-constraints, Atmos. Chem. Phys., 12, 7779–7795,, 2012. a

Hui, J. and Hong, L.: Projected Changes in NOx Emissions from Lightning as a Result of 2000–2050 Climate Change, Atmospheric and Oceanic Science Letters, 6, 284–289,, 2013. a

Husain, L., Coffey, P., Meyers, R., and Cederwall, R.: Ozone Transport from Stratosphere to Troposphere, Geophys. Res. Lett., 4, 363–365,, 1977. a

Igarashi, Y., Hirose, I., and Otsuji-Hatori, M.: Beryllium-7 deposition and its relation to sulfate deposition, J. Atmos. Chem., 29, 217–231,, 1998. a

Jacob, D. J. and Prather, M. J.: Radon-222 as a test of convective transport in a general circulation model, Tellus B, 42, 118–134,, 1990. a

Jacob, D. J., Prather, M., Rasch, P., Shia, R., Balkanski, Y., Beagley, S., Bergmann, D., Blackshear, W., Brown, M., Chiba, M., Chipperfield, M., deGrandpre, J., Dignon, J., Feichter, J., Genthon, C., Grose, W., Kasibhatla, P., Kohler, I., Kritz, M., Law, K., PENNER, J., Ramonet, M., Reeves, C., Rotman, D., Stockwell, D., VanVelthoven, P., Verver, G., Wild, O., Yang, H., and Zimmermann, P.: Evaluation and intercomparison of global atmospheric transport models using Rn-222 and other short-lived tracers, J. Geophys. Res.-Atmos, 102, 5953–5970,, 1997. a, b

Jaeglé, L., Quinn, P. K., Bates, T. S., Alexander, B., and Lin, J.-T.: Global distribution of sea salt aerosols: new constraints from in situ and remote sensing observations, Atmos. Chem. Phys., 11, 3137–3157,, 2011. a

Jiang, L. and O'Neill, B. C.: Global urbanization projections for the Shared Socioeconomic Pathways, Glob. Environ. Change, 42, 193–199,, 2017. a

Kaplan, J. O., Folberth, G., and Hauglustaine, D. A.: Role of methane and biogenic volatile organic compound sources in late glacial and Holocene fluctuations of atmospheric methane concentrations, Global Biogeochem. Cy., 20, GB2016,, 2006. a

Keller, C. A., Long, M. S., Yantosca, R. M., Da Silva, A. M., Pawson, S., and Jacob, D. J.: HEMCO v1.0: a versatile, ESMF-compliant component for calculating emissions in atmospheric models, Geosci. Model Dev., 7, 1409–1417,, 2014. a, b, c

Kelley, M., Schmidt, G. A., Nazarenko, L. S., Bauer, S. E., Ruedy, R., Russell, G. L., Ackerman, A. S., Aleinov, I., Bauer, M., Bleck, R., Canuto, V., Cesana, G., Cheng, Y., Clune, T. L., Cook, B. I., Cruz, C. A., Del Genio, A. D., Elsaesser, G. S., Faluvegi, G., Kiang, N. Y., Kim, D., Lacis, A. A., Leboissetier, A., LeGrande, A. N., Lo, K. K., Marshall, J., Matthews, E. E., McDermid, S., Mezuman, K., Miller, R. L., Murray, L. T., Oinas, V., Orbe, C., García-Pando, C. P., Perlwitz, J. P., Puma, M. J., Rind, D., Romanou, A., Shindell, D. T., Sun, S., Tausnev, N., Tsigaridis, K., Tselioudis, G., Weng, E., Wu, J., and Yao, M.-S.: GISS-E2.1: Configurations and Climatology, J. Adv. Model. Earth Syst., 12, e2019MS002025,, 2020. a, b, c, d

Kim, D., Sobel, A. H., Del Genio, A. D., Chen, Y., Camargo, S. J., Yao, M.-S., Kelley, M., and Nazarenko, L.: The Tropical Subseasonal Variability Simulated in the NASA GISS General Circulation Model, J. Climate, 25, 4641–4659,, 2012. a

Koch, D., Jacob, D. J., and Graustein, W.: Vertical transport of tropospheric aerosols as indicated by Be-7 and Pb-210 in a chemical tracer model, J. Geophys. Res.-Atmos., 101, 18651–18666,, 1996. a, b, c, d

Kopacz, M., Jacob, D. J., Henze, D. K., Heald, C. L., Streets, D. G., and Zhang, Q.: Comparison of adjoint and analytical Bayesian inversion methods for constraining Asian sources of carbon monoxide using satellite (MOPITT) measurements of CO columns, J. Geophys. Res., 114, D04305,, 2009. a

Kriegler, E., Bauer, N., Popp, A., Humpenöder, F., Leimbach, M., Strefler, J., Baumstark, L., Bodirsky, B. L., Hilaire, J., Klein, D., Mouratiadou, I., Weindl, I., Bertram, C., Dietrich, J.-P., Luderer, G., Pehl, M., Pietzcker, R., Piontek, F., Lotze-Campen, H., Biewald, A., Bonsch, M., Giannousakis, A., Kreidenweis, U., Müller, C., Rolinski, S., Schultes, A., Schwanitz, J., Stevanovic, M., Calvin, K., Emmerling, J., Fujimori, S., and Edenhofer, O.: Fossil-fueled development (SSP5): An energy and resource intensive scenario for the 21st century, Glob. Environ. Change, 42, 297–315,, 2017. a

Kritz, M., Rosner, S., and Stockwell, D.: Validation of an off-line three-dimensional chemical transport model using observed radon profiles – 1. Observations, J. Geophys. Res.-Atmos., 103, 8425–8432,, 1998. a

Krotkov, N. A., Lamsal, L. N., Celarier, E. A., Swartz, W. H., Marchenko, S. V., Bucsela, E. J., Chan, K. L., Wenig, M., and Zara, M.: The version 3 OMI NO2 standard product, Atmos. Meas. Tech., 10, 3133–3149,, 2017. a, b

Krotkov, N. A., Lamsal, L. N., Marchenko, S. V., Celarier, E. A., Bucsela, E. J., Swartz, W. H., Joiner, J., and the OMI core team: OMI/Aura NO2 Cloud-Screened Total and Tropospheric Column L3 Global Gridded 0.25 degree x 0.25 degree V3, NASA Goddard Space Flight Center, Goddard Earth Sciences Data and Information Services Center (GES DISC),, 2019. a

Lal, D., Malhotra, P., and Peters, B.: On the Production of Radioisotopes in the Atmosphere by Cosmic Radiation and Their Application to Meteorology, J. Atmos. Terr. Phys., 12, 306–328,, 1958. a

Lambert, G., Polian, G., Sanak, J., Ardouin, B., Buisson, A., Jegou, A., and Leroulley, J.: Cycle du radon et de ses descendants: application à l'étude des échanges troposphère-stratosphère, Ann Geophys, 38, 497–531, 1982. a, b

Leimbach, M., Kriegler, E., Roming, N., and Schwanitz, J.: Future growth patterns of world regions – A GDP scenario approach, Glob. Environ. Change, 42, 215–225,, 2017. a

Li, C., Krotkov, N. A., Leonard, P. J. T., Carn, S., Joiner, J., Spurr, R. J. D., and Vasilkov, A.: Version 2 Ozone Monitoring Instrument SO2 product (OMSO2 V2): new anthropogenic SO2 vertical column density dataset, Atmos. Meas. Tech., 13, 6175–6191,, 2020a. a, b

Li, C., Krotkov, N. A., Leonard, P.: OMI/Aura Sulfur Dioxide (SO2) Total Column L3 1 day Best Pixel in 0.25 degree x 0.25 degree V3, Greenbelt, MD, USA, Goddard Earth Sciences Data and Information Services Center (GES DISC),, 2020b. a

Lin, S.-J. and Rood, R. B.: Multidimensional Flux-Form Semi-Lagrangian Transport Schemes, Mon. Weather Rev., 124, 2046–2070,<2046:MFFSLT>2.0.CO;2, 1996. a, b

Liu, H., Jacob, D. J., Bey, I., and Yantosca, R.: Constraints from Pb-210 and Be-7 on wet deposition and transport in a global three-dimensional chemical tracer model driven by assimilated meteorological fields, J. Geophys. Res.-Atmos, 106, 12109–12128,, 2001. a, b, c, d

Liu, H., Considine, D. B., Horowitz, L. W., Crawford, J. H., Rodriguez, J. M., Strahan, S. E., Damon, M. R., Steenrod, S. D., Xu, X., Kouatchou, J., Carouge, C., and Yantosca, R. M.: Using beryllium-7 to assess cross-tropopause transport in global models, Atmos. Chem. Phys., 16, 4641–4659,, 2016. a, b, c

Lundgren, L., Bindle, L., Yantosca, B., Sulprizio, M., Downs, W., and Eastham, S. D.: geoschem/GCHP: GCHP 13.1.0, Zenodo [code], ​​​, 2021. a

Maenhaut, W., Zoller, W., and Coles, D.: Radionuclides in the South Pole Atmosphere, J. Geophys. Res.-Oc. Atm., 84, 3131–3138,, 1979. a, b

Mahowald, N., Rasch, P., and Prinn, R.: Cumulus parameterizations in chemical transport models, J. Geophys. Res.-Atmos., 100, 26173–26189,, 1995. a

Maiss, M. and Brenninkmeijer, C. A. M.: Atmospheric SF6: Trends, Sources, and Prospects, Environ. Sci. Technol., 32, 3077–3086,, 1998. a

Meinshausen, M., Vogel, E., Nauels, A., Lorbacher, K., Meinshausen, N., Etheridge, D. M., Fraser, P. J., Montzka, S. A., Rayner, P. J., Trudinger, C. M., Krummel, P. B., Beyerle, U., Canadell, J. G., Daniel, J. S., Enting, I. G., Law, R. M., Lunder, C. R., O'Doherty, S., Prinn, R. G., Reimann, S., Rubino, M., Velders, G. J. M., Vollmer, M. K., Wang, R. H. J., and Weiss, R.: Historical greenhouse gas concentrations for climate modelling (CMIP6), Geosci. Model Dev., 10, 2057–2116,, 2017. a, b, c

Menon, S., Del Genio, A. D., Kaufman, Y., Bennartz, R., Koch, D., Loeb, N., and Orlikowski, D.: Analyzing signatures of aerosol-cloud interactions from satellite retrievals and the GISS GCM to constrain the aerosol indirect effect, J. Geophys. Res., 113, D14S22,, 2008. a

Miller, R. L., Schmidt, G. A., Nazarenko, L. S., Bauer, S. E., Kelley, M., Ruedy, R., Russell, G. L., Ackerman, A. S., Aleinov, I., Bauer, M., Bleck, R., Canuto, V., Cesana, G., Cheng, Y., Clune, T. L., Cook, B. I., Cruz, C. A., Del Genio, A. D., Elsaesser, G. S., Faluvegi, G., Kiang, N. Y., Kim, D., Lacis, A. A., Leboissetier, A., LeGrande, A. N., Lo, K. K., Marshall, J., Matthews, E. E., McDermid, S., Mezuman, K., Murray, L. T., Oinas, V., Orbe, C., Pérez García-Pando, C., Perlwitz, J. P., Puma, M. J., Rind, D., Romanou, A., Shindell, D. T., Sun, S., Tausnev, N., Tsigaridis, K., Tselioudis, G., Weng, E., Wu, J., and Yao, M.-S.: CMIP6 Historical Simulations (1850–2014) With GISS-E2.1, J. Adv. Model. Earth Syst., 13, e2019MS002034,, 2021. a, b, c

Moore, H., Poet, S., and Martell, E.: Rn-222, Pb-210, Bi-210, and Po-210 Profiles and Aerosol Residence Times Versus Altitude, J. Geophys. Res., 78, 7065–7075,, 1973. a

Murray, L. T.: Lightning NOx and Impacts on Air Quality, Curr. Pollut. Rep., 2, 115–133,, 2016. a

Murray, L. T.: An uncertain future for lightning, Nat. Clim. Change, 8, 191–192,, 2018. a

Murray, L. T.: GCAP 2.0 input files, GCAP Data Repository [data set], available at:, last access: 9 September 2021. a

Murray, L. T., Jacob, D. J., Logan, J. A., Hudman, R. C., and Koshak, W. J.: Optimized regional and interannual variability of lightning in a global chemical transport model constrained by LIS/OTD satellite data, J. Geophys. Res.-Atmos., 117, D20307,, 2012. a, b

Murray, L. T., Logan, J. A., and Jacob, D. J.: Interannual variability in tropical tropospheric ozone and OH: The role of lightning, J. Geophys. Res.-Atmos., 118, 11468–11480,, 2013. a

Murray, L. T., Mickley, L. J., Kaplan, J. O., Sofen, E. D., Pfeiffer, M., and Alexander, B.: Factors controlling variability in the oxidative capacity of the troposphere since the Last Glacial Maximum, Atmos. Chem. Phys., 14, 3589–3622,, 2014. a, b, c, d, e, f, g, h, i, j

Murray, L. T., Leibensperger, E. M., and Mickley, L. J.: MERRA-2 like diagnostics for the GISS ModelE2.1 GCM, Zenodo [code],, 2021. a

Naik, V., Voulgarakis, A., Fiore, A. M., Horowitz, L. W., Lamarque, J.-F., Lin, M., Prather, M. J., Young, P. J., Bergmann, D., Cameron-Smith, P. J., Cionni, I., Collins, W. J., Dalsøren, S. B., Doherty, R., Eyring, V., Faluvegi, G., Folberth, G. A., Josse, B., Lee, Y. H., MacKenzie, I. A., Nagashima, T., van Noije, T. P. C., Plummer, D. A., Righi, M., Rumbold, S. T., Skeie, R., Shindell, D. T., Stevenson, D. S., Strode, S., Sudo, K., Szopa, S., and Zeng, G.: Preindustrial to present-day changes in tropospheric hydroxyl radical and methane lifetime from the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), Atmos. Chem. Phys., 13, 5277–5298,, 2013. a

Nakicenovic, N. and Swart, R. (Eds.): IPCC Special Report on Emissions Scenarios, Cambridge University Press, Cambridge, UK, 2000. a

Narazaki, Y. and Fujitaka, K.: The Geographical Distribution and Features of 7Be Deposition in Japan, Japan Health Physics Society, 37, 317–324, 2010. a

Nazarov, L., Kuzenkov, A., Malakhov, S., Volokitina, L., Gaziyev, Y., and Vasilyev, A.: Radioactive Aerosol Distribution in Middle and Upper Troposphere over USSR in 1963-1968, J. Geophys. Res., 75, 3575,, 1970. a

Nightingale, P. D., Liss, P. S., and Schlosser, P.: Measurements of air-sea gas transfer during an open ocean algal bloom, Geophys. Res. Lett., 27, 2117–2120,, 2000a. a

Nightingale, P. D., Malin, G., Law, C. S., Watson, A. J., Liss, P. S., Liddicoat, M. I., Boutin, J., and Upstill-Goddard, R. C.: In situ evaluation of air-sea gas exchange parameterizations using novel conservative and volatile tracers, Global Biogeochem. Cy., 14, 373–387,, 2000b. a

Nijampurkar, V. and Rao, D.: Polar Fallout of Radionuclides Si-32, Be-7 and Pb-210 and Past Accumulation Rate of Ice at Indian Station, Dakshin Gangotri, East Antarctica, J. Environ. Radioactiv., 21, 107–117,, 1993. a

NOAA Carbon Cycle Group ObsPack Team: Multi-laboratory compilation of atmospheric sulfur hexafluoride data for the period 1983–2017; obspack_sf6_1_v2.1.1_2018-08-17, NOAA Earth System Research Laboratory, Global Monitoring Division,, 2018. a

Olsen, C., Larsen, I., Lowry, P., Cutshall, N., Todd, J., Wong, G., and Casey, W.: Atmospheric Fluxes and Marsh-Soil Inventories of Be-7 and Pb-210, J. Geophys. Res.-Atmos., 90, 10487–10495,, 1985. a

O'Neill, B. C., Tebaldi, C., van Vuuren, D. P., Eyring, V., Friedlingstein, P., Hurtt, G., Knutti, R., Kriegler, E., Lamarque, J.-F., Lowe, J., Meehl, G. A., Moss, R., Riahi, K., and Sanderson, B. M.: The Scenario Model Intercomparison Project (ScenarioMIP) for CMIP6, Geosci. Model Dev., 9, 3461–3482,, 2016. a

Orbe, C., Plummer, D. A., Waugh, D. W., Yang, H., Jöckel, P., Kinnison, D. E., Josse, B., Marecal, V., Deushi, M., Abraham, N. L., Archibald, A. T., Chipperfield, M. P., Dhomse, S., Feng, W., and Bekki, S.: Description and Evaluation of the specified-dynamics experiment in the Chemistry-Climate Model Initiative , Atmos. Chem. Phys., 20, 3809–3840,, 2020a. a, b

Orbe, C., Rind, D., Jonas, J., Nazarenko, L., Faluvegi, G., Murray, L. T., Shindell, D. T., Tsigaridis, K., Zhou, T., Kelley, M., and Schmidt, G. A.: GISS Model E2.2: A Climate Model Optimized for the Middle Atmosphere – 2. Validation of Large-Scale Transport and Evaluation of Climate Response, J. Geophys. Res.-Atmos., 125, e2020JD033151,, 2020b. a, b, c, d

Papastefanou, C.: Beryllium-7 Aerosols in Ambient Air, Aerosol Air Qual. Res., 9, 187–197,, 2009. a

Papastefanou, C. and Ioannidou, A.: Beryllium-7 aerosols in ambient air, Environ. Int., 22, S125–S130,, 1996. a

Papastefanou, C., Ioannidou, A., Stoulos, S., and Manolopoulou, M.: Atmospheric Deposition of Cosmogenic Be-7 and Cs-137 from Fallout of the Chernobyl Accident, Sci. Total Environ., 170, 151–156,, 1995. a

Park, R. J., Jacob, D. J., Field, B. D., Yantosca, R., and Chin, M.: Natural and transboundary pollution influences on sulfate-nitrate-ammonium aerosols in the United States: Implications for policy, J. Geophys. Res., 109, D15204,, 2004. a, b, c

Pfeiffer, M., Spessa, A., and Kaplan, J. O.: A model for global biomass burning in preindustrial time: LPJ-LMfire (v1.0), Geosci. Model Dev., 6, 643–685,, 2013. a

Platnick, S., et al.: MODIS Atmosphere L3 Monthly Product. NASA MODIS Adaptive Processing System, Goddard Space Flight Center, USA,, 2015. a

Popp, A., Calvin, K., Fujimori, S., Havlik, P., Humpenöder, F., Stehfest, E., Bodirsky, B. L., Dietrich, J. P., Doelmann, J. C., Gusti, M., Hasegawa, T., Kyle, P., Obersteiner, M., Tabeau, A., Takahashi, K., Valin, H., Waldhoff, S., Weindl, I., Wise, M., Kriegler, E., Lotze-Campen, H., Fricko, O., Riahi, K., and Vuuren, D. P. v.: Land-use futures in the shared socio-economic pathways, Glob. Environ. Change, 42, 331–345,, 2017. a

Prather, M. J.: Numerical advection by conservation of second-order moments, J. Geophys. Res., 91, 6671,, 1986. a

Prather, M. J., Holmes, C. D., and Hsu, J.: Reactive greenhouse gas scenarios: Systematic exploration of uncertainties and the role of atmospheric chemistry, Geophys. Res. Lett., 39, L09803,, 2012. a, b

Preiss, N., Melieres, M., and Pourchet, M.: A compilation of data on lead 210 concentration in surface air and fluxes at the air-surface and water-sediment interfaces, J. Geophys. Res.-Atmos., 101, 28847–28862,, 1996. a, b, c

Price, C. and Rind, D.: A simple lightning parameterization for calculating global lightning distributions, J. Geophys. Res.-Atmos., 97, 9919–9933,, 1992. a

Price, C. G.: Lightning Applications in Weather and Climate Research, Surv. Geophys., 34, 755–767,, 2013. a

Prinn, R., Huang, J., Weiss, R., Cunnold, D., Fraser, P., Simmonds, P., McCulloch, A., Harth, C., Reimann, S., Salameh, P., O'Doherty, S., Wang, R., Porter, L., Miller, B., and Krummel, P.: Evidence for variability of atmospheric hydroxyl radicals over the past quarter century, Geophys. Res. Lett., 32, L07809,, 2005. a, b

Pye, H. O. T. and Seinfeld, J. H.: A global perspective on aerosol from low-volatility organic compounds, Atmos. Chem. Phys., 10, 4377–4401,, 2010. a

Pye, H. O. T., Liao, H., Wu, S., Mickley, L. J., Jacob, D. J., Henze, D. K., and Seinfeld, J. H.: Effect of changes in climate and emissions on future sulfate-nitrate-ammonium aerosol levels in the United States, J. Geophys. Res.-Atmos, 114, D01205,, 2009. a

Rao, S., Klimont, Z., Smith, S. J., Van Dingenen, R., Dentener, F., Bouwman, L., Riahi, K., Amann, M., Bodirsky, B. L., van Vuuren, D. P., Aleluia Reis, L., Calvin, K., Drouet, L., Fricko, O., Fujimori, S., Gernaat, D., Havlik, P., Harmsen, M., Hasegawa, T., Heyes, C., Hilaire, J., Luderer, G., Masui, T., Stehfest, E., Strefler, J., van der Sluis, S., and Tavoni, M.: Future air pollution in the Shared Socio-economic Pathways, Glob. Environ. Change, 42, 346–358,, 2017. a

Ray, E. A., Moore, F. L., Rosenlof, K. H., Davis, S. M., Boenisch, H., Morgenstern, O., Smale, D., Rozanov, E., Hegglin, M., Pitari, G., Mancini, E., Braesicke, P., Butchart, N., Hardiman, S., Li, F., Shibata, K., and Plummer, D. A.: Evidence for changes in stratospheric transport and mixing over the past three decades based on multiple data sets and tropical leaky pipe analysis, J. Geophys. Res., 115, D21304,, 2010. a

Rehfeld, S. and Heimann, M.: Three dimensional atmospheric transport simulation of the radioactive tracers Pb-210, Be-7, Be-10, and Sr-90, J. Geophys. Res.-Atmos., 100, 26141–26161,, 1995. a

Riahi, K., van Vuuren, D. P., Kriegler, E., Edmonds, J., O'Neill, B. C., Fujimori, S., Bauer, N., Calvin, K., Dellink, R., Fricko, O., Lutz, W., Popp, A., Cuaresma, J. C., KC, S., Leimbach, M., Jiang, L., Kram, T., Rao, S., Emmerling, J., Ebi, K., Hasegawa, T., Havlik, P., Humpenöder, F., Da Silva, L. A., Smith, S., Stehfest, E., Bosetti, V., Eom, J., Gernaat, D., Masui, T., Rogelj, J., Strefler, J., Drouet, L., Krey, V., Luderer, G., Harmsen, M., Takahashi, K., Baumstark, L., Doelman, J. C., Kainuma, M., Klimont, Z., Marangoni, G., Lotze-Campen, H., Obersteiner, M., Tabeau, A., and Tavoni, M.: The Shared Socioeconomic Pathways and their energy, land use, and greenhouse gas emissions implications: An overview, Glob. Environ. Change, 42, 153–168,, 2017. a, b, c

Rigby, M., Mühle, J., Miller, B. R., Prinn, R. G., Krummel, P. B., Steele, L. P., Fraser, P. J., Salameh, P. K., Harth, C. M., Weiss, R. F., Greally, B. R., O'Doherty, S., Simmonds, P. G., Vollmer, M. K., Reimann, S., Kim, J., Kim, K.-R., Wang, H. J., Olivier, J. G. J., Dlugokencky, E. J., Dutton, G. S., Hall, B. D., and Elkins, J. W.: History of atmospheric SF6 from 1973 to 2008, Atmos. Chem. Phys., 10, 10305–10320,, 2010. a

Rind, D. and Lerner, J.: Use of on-line tracers as a diagnostic tool in general circulation model development: 1. Horizontal and vertical transport in the troposphere, J. Geophys. Res.-Atmos., 101, 12667–12683,, 1996. a

Rind, D., Suozzo, R., Balachandran, N. K., Lacis, A., and Russell, G.: The GISS Global Climate-Middle Atmosphere Model. Part I: Model Structure and Climatology, J. Atmos. Sci., 45, 329–370,<0329:TGGCMA>2.0.CO;2, 1988. a

Rind, D., Lerner, J., Jonas, J., and McLinden, C.: Effects of resolution and model physics on tracer transports in the NASA Goddard Institute for Space Studies general circulation models, J. Geophys. Res., 112, D09315,, 2007. a

Rind, D., Orbe, C., Jonas, J., Nazarenko, L., Zhou, T., Kelley, M., Lacis, A., Shindell, D., Faluvegi, G., Romanou, A., Russell, G., Tausnev, N., Bauer, M., and Schmidt, G.: GISS Model E2.2: A Climate Model Optimized for the Middle Atmosphere – Model Structure, Climatology, Variability, and Climate Sensitivity, J. Geophys. Res.-Atmos., 125, e2019JD032204,, 2020. a

Samanta, D., Karnauskas, K. B., and Goodkin, N. F.: Tropical Pacific SST and ITCZ Biases in Climate Models: Double Trouble for Future Rainfall Projections, Geophys. Res. Lett., 46, 2242–2252,, 2019. a

Samir, K. and Lutz, W.: The human core of the shared socioeconomic pathways: Population scenarios by age, sex and level of education for all countries to 2100., Glob. Environ. Change, 42, 181–192,, 2017. a

Sanak, J., Gaudry, A., and Lambert, G.: Size Distribution of Pb-210 Aerosols Over Oceans, Geophys. Res. Lett., 8, 1067–1069,, 1981. a, b

Sanak, J., Lambert, G., and Ardouin, B.: Measurement of Stratosphere-to-Troposphere Exchange in Antarctica by Using Short-Lived Cosmonuclides, Tellus B, 37, 109–115,, 1985. a

Schmidt, G. A., Ruedy, R., Hansen, J. E., Aleinov, I., Bell, N., Bauer, M., Bauer, S., Cairns, B., Canuto, V., Cheng, Y., Del Genio, A., Faluvegi, G., Friend, A. D., Hall, T. M., Hu, Y., Kelley, M., Kiang, N. Y., Koch, D., Lacis, A. A., Lerner, J., Lo, K. K., Miller, R. L., Nazarenko, L., Oinas, V., Perlwitz, J., Perlwitz, J., Rind, D., Romanou, A., Russell, G. L., Sato, M., Shindell, D. T., Stone, P. H., Sun, S., Tausnev, N., Thresher, D., and Yao, M.-S.: Present-Day Atmospheric Simulations Using GISS ModelE: Comparison to In Situ, Satellite, and Reanalysis Data, J. Climate, 19, 153–192,, 2006. a

Schubert, S., Rood, R. B., and Pfaendtner, J.: An Assimilated Dataset for Earth Science Applications, Bull. Am. Meteorol. Soc., 74, 2331–2342,<2331:AADFES>2.0.CO;2, 1993. a

Schuler, C., Wieland, E., Santschi, P., Sturm, M., Lueck, A., Bollhalder, S., Beer, J., Bonani, G., Hofmann, H., Suter, M., and Wolfli, W.: A Multitracer Studyof Radionuclides in Lake Zurich, Switzerland .1. Comparison of Atmospheric and Sedimentary Fluxes of Be-7, Be-10, Pb-210, Po-210, and Cs-137, J. Geophys. Res.pOceans, 96, 17051–17065,, 1991. a

Schultz, M. G., Schröder, S., Lyapina, O., Cooper, O., Galbally, I., Petropavlovskikh, I., Von Schneidemesser, E., Tanimoto, H., Elshorbany, Y., Naja, M., Seguel, R., Dauert, U., Eckhardt, P., Feigenspahn, S., Fiebig, M., Hjellbrekke, A.-G., Hong, Y.-D., Christian Kjeld, P., Koide, H., Lear, G., Tarasick, D., Ueno, M., Wallasch, M., Baumgardner, D., Chuang, M.-T., Gillett, R., Lee, M., Molloy, S., Moolla, R., Wang, T., Sharps, K., Adame, J. A., Ancellet, G., Apadula, F., Artaxo, P., Barlasina, M., Bogucka, M., Bonasoni, P., Chang, L., Colomb, A., Cuevas, E., Cupeiro, M., Degorska, A., Ding, A., Fröhlich, M., Frolova, M., Gadhavi, H., Gheusi, F., Gilge, S., Gonzalez, M. Y., Gros, V., Hamad, S. H., Helmig, D., Henriques, D., Hermansen, O., Holla, R., Huber, J., Im, U., Jaffe, D. A., Komala, N., Kubistin, D., Lam, K.-S., Laurila, T., Lee, H., Levy, I., Mazzoleni, C., Mazzoleni, L., McClure-Begley, A., Mohamad, M., Murovic, M., Navarro-Comas, M., Nicodim, F., Parrish, D., Read, K. A., Reid, N., Ries, L., Saxena, P., Schwab, J. J., Scorgie, Y., Senik, I., Simmonds, P., Sinha, V., Skorokhod, A., Spain, G., Spangl, W., Spoor, R., Springston, S. R., Steer, K., Steinbacher, M., Suharguniyawan, E., Torre, P., Trickl, T., Weili, L., Weller, R., Xu, X., Xue, L., and Zhiqiang, M.: Tropospheric Ozone Assessment Report: Database and Metrics Data of Global Surface Ozone Observations, Elem. Sci. Anth., 5, 58,, 2017a. a, b

Schultz, M. G., Schröder, S., Lyapina, O., Cooper, O. R., Galbally, I., Petropavlovskikh, I., von Schneidemesser, E., Tanimoto, H., Elshorbany, Y., Naja, M., Seguel, R. J., Dauert, U., Eckhardt, P., Feigenspan, S., Fiebig, M., Hjellbrekke, A.-G., Hong, Y.-D., Kjeld, P. C., Koide, H., Lear, G., Tarasick, D., Ueno, M., Wallasch, M., Baumgardner, D., Chuang, M.-T., Gillett, R., Lee, M., Molloy, S., Moolla, R., Wang, T., Sharps, K., Adame, J. A., Ancellet, G., Apadula, F., Artaxo, P., Barlasina, M. E., Bogucka, M., Bonasoni, P., Chang, L., Colomb, A., Cuevas-Agulló, E., Cupeiro, M., Degorska, A., Ding, A., Fröhlich, M., Frolova, M., Gadhavi, H., Gheusi, F., Gilge, S., Gonzalez, M. Y., Gros, V., Hamad, S. H., Helmig, D., Henriques, D., Hermansen, O., Holla, R., Hueber, J., Im, U., Jaffe, D. A., Komala, N., Kubistin, D., Lam, K.-S., Laurila, T., Lee, H., Levy, I., Mazzoleni, C., Mazzoleni, L. R., McClure-Begley, A., Mohamad, M., Murovec, M., Navarro-Comas, M., Nicodim, F., Parrish, D., Read, K. A., Reid, N., Ries, L., Saxena, P., Schwab, J. J., Scorgie, Y., Senik, I., Simmonds, P., Sinha, V., Skorokhod, A. I., Spain, G., Spangl, W., Spoor, R., Springston, S. R., Steer, K., Steinbacher, M., Suharguniyawan, E., Torre, P., Trickl, T., Weili, L., Weller, R., Xu, X., Xue, L., and Zhiqiang, M.​​​​​​​: Tropospheric Ozone Assessment Report, links to Global surface ozone datasets, PANGAEA [data set],, 2017b. a

Selin, N. E., Wu, S., Nam, K. M., Reilly, J. M., Paltsev, S., Prinn, R. G., and Webster, M. D.: Global health and economic impacts of future ozone pollution, Environ. Res. Lett., 4, 044014,, 2009. a

Stettler, M., Eastham, S., and Barrett, S.: Air quality and public health impacts of UK airports. Part I: Emissions, Atmos. Environ., 45, 5415–5424,, 2011. a

Stockwell, D., Kritz, M., Chipperfield, M., and Pyle, J.: Validation of an off-line three-dimensional chemical transport model using observed radon profiles – 2. Model results, J. Geophys. Res.-Atmos., 103, 8433–8445,, 1998. a

Tai, A. P. K., Mickley, L. J., Heald, C. L., and Wu, S.: Effect of CO2 inhibition on biogenic isoprene emission: Implications for air quality under 2000 to 2050 changes in climate, vegetation, and land use, Geophys. Res. Lett., 40, 3479–3483,, 2013. a

Tai, A. P. K., Sadiq, M., Pang, J. Y. S., Yung, D. H. Y., and Feng, Z.: Impacts of Surface Ozone Pollution on Global Crop Yields: Comparing Different Ozone Exposure Metrics and Incorporating Co-effects of CO2, Frontiers in Sustainable Food Systems, 5, 534616,, 2021. a

The International GEOS-Chem User Community: geoschem/geos-chem: GEOS-Chem 12.4.0, Zenodo [code],, 2019. a

The International GEOS-Chem User Community: geoschem/geos-chem: GEOS-Chem 12.9.3, Zenodo [code],, 2020. a

The International GEOS-Chem User Community​​​​​​​: geoschem/GCClassic: GEOS-Chem 13.1.0, Zenodo [code],, 2021a. a, b

The International GEOS-Chem User Community: geoschem/HEMCO: HEMCO 3.0.0, Zenodo [code], ​​​, 2021b. a

The International GEOS-Chem User Community: geoschem/GCClassic: GEOS-Chem 13.0.0, Zenodo [code], ​​​, 2021c. a

Tian, B., Manning, E., Roman, J., Thrastarson, H., Fetzer, E., and Monarrez, R.: AIRS Version 7 Level 3 Product User Guide, Version 1.1, Tech. rep., Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA, 2020. a, b

Turekian, K., Nozaki, Y., and Benninger, L.: Geochemistry of Atmospheric Radon and Radon Products, Annu. Rev. Earth Planet. Sci., 5, 227–255,, 1977. a

Turekian, K., Benninger, L., and Dion, E.​​​​​​​: 7Be and 210Pb Total Deposition Fluxes at New-Haven, Connecticut and at Bermuda, J. Geophys. Res.-Oc. Atm., 88, 5411–5415,, 1983. a

Usoskin, I., Alanko-Huotari, K., Kovaltsov, G., and Mursula, K.: Heliospheric modulation of cosmic rays: Monthly reconstruction for 1951-2004, J. Geophys. Res.-Space, 110, A12108,, 2005. a, b

Usoskin, I. G. and Kovaltsov, G. A.: Production of cosmogenic Be-7 isotope in the atmosphere: Full 3-D modeling, J. Geophys. Res.-Atmos., 113, D12107,, 2008. a

van der Werf, G. R., Randerson, J. T., Giglio, L., van Leeuwen, T. T., Chen, Y., Rogers, B. M., Mu, M., van Marle, M. J. E., Morton, D. C., Collatz, G. J., Yokelson, R. J., and Kasibhatla, P. S.: Global fire emissions estimates during 1997–2016, Earth Syst. Sci. Data, 9, 697–720,, 2017. a

van Marle, M. J. E., Kloster, S., Magi, B. I., Marlon, J. R., Daniau, A.-L., Field, R. D., Arneth, A., Forrest, M., Hantson, S., Kehrwald, N. M., Knorr, W., Lasslop, G., Li, F., Mangeon, S., Yue, C., Kaiser, J. W., and van der Werf, G. R.: Historic global biomass burning emissions for CMIP6 (BB4CMIP) based on merging satellite observations with proxies and fire models (1750–2015), Geosci. Model Dev., 10, 3329–3357,, 2017. a, b

van Vuuren, D. P., Stehfest, E., Gernaat, D. E., Doelman, J. C., van den Berg, M., Harmsen, M., de Boer, H. S., Bouwman, L. F., Daioglou, V., Edelenbosch, O. Y., Girod, B., Kram, T., Lassaletta, L., Lucas, P. L., van Meijl, H., Müller, C., van Ruijven, B. J., van der Sluis, S., and Tabeau, A.: Energy, land-use and greenhouse gas emissions trajectories under a green growth paradigm, Glob. Environ. Change, 42, 237–250,, 2017. a

Veefkind, P.: OMI/Aura Ozone (O3) DOAS Total Column L3 1 day 0.25 degree x 0.25 degree V3, Greenbelt, MD, USA, Goddard Earth Sciences Data and Information Services Center (GES DISC),, 2012. a

Viezee, W. and Singh, H.: The Distribution of Beryllium-7 in the Troposphere – Implications on Stratospheric-Tropospheric Air Exchange, Geophys. Res. Lett., 7, 805–808,, 1980. a

Voulgarakis, A., Naik, V., Lamarque, J.-F., Shindell, D. T., Young, P. J., Prather, M. J., Wild, O., Field, R. D., Bergmann, D., Cameron-Smith, P., Cionni, I., Collins, W. J., Dalsøren, S. B., Doherty, R. M., Eyring, V., Faluvegi, G., Folberth, G. A., Horowitz, L. W., Josse, B., MacKenzie, I. A., Nagashima, T., Plummer, D. A., Righi, M., Rumbold, S. T., Stevenson, D. S., Strode, S. A., Sudo, K., Szopa, S., and Zeng, G.: Analysis of present day and future OH and methane lifetime in the ACCMIP simulations, Atmos. Chem. Phys., 13, 2563–2587,, 2013. a, b

Wallbrink, P. and Murray, A.: Fallout of Be-7 in South Eastern Australia, J. Environ. Radioactiv., 25, 213–228,, 1994. a

Wang, X., Jacob, D. J., Eastham, S. D., Sulprizio, M. P., Zhu, L., Chen, Q., Alexander, B., Sherwen, T., Evans, M. J., Lee, B. H., Haskins, J. D., Lopez-Hilfiker, F. D., Thornton, J. A., Huey, G. L., and Liao, H.: The role of chlorine in global tropospheric chemistry, Atmos. Chem. Phys., 19, 3981–4003,, 2019. a

Wang, Y. and Jacob, D. J.: Anthropogenic forcing on tropospheric ozone and OH since preindustrial times, J. Geophys. Res.-Atmos., 103, 31123–31135,, 1998. a

Wang, Y., Jacob, D. J., and Logan, J. A.: Global simulation of tropospheric O3-NOx-hydrocarbon chemistry: 1. Model formulation, J. Geophys. Res.-Atmos, 103, 10713–10725,, 1998a. a

Wang, Y., Jacob, D. J., and Logan, J. A.: Global simulation of tropospheric O3-NOx-hydrocarbon chemistry: 3. Origin of tropospheric ozone and effects of nonmethane hydrocarbons, J. Geophys. Res.-Atmos., 103, 10757–10767,, 1998b. a

Wang, Y., Logan, J. A., and Jacob, D. J.: Global simulation of tropospheric O3-NOx-hydrocarbon chemistry: 2. Model evaluation and global ozone budget, J. Geophys. Res.-Atmos., 103, 10727–10755,, 1998c. a

Wang, Y. X., McElroy, M. B., Jacob, D. J., and Yantosca, R. M.: A nested grid formulation for chemical transport over Asia: Applications to CO, J. Geophys. Res.-Atmos., 109, D22307,, 2004. a

Waugh, D.: Atmospheric dynamics: The age of stratospheric air, Nat. Geosci., 2, 14–16,, 2009. a

Waugh, D. W. and Hall, T. M.: Age of stratospheric air: Theory, observations, and models, Rev. Geophys., 40, 1010,, 2002. a, b, c

Waugh, D. W., Crotwell, A. M., Dlugokencky, E. J., Dutton, G. S., Elkins, J. W., Hall, B. D., Hintsa, E. J., Hurst, D. F., Montzka, S. A., Mondeel, D. J., Moore, F. L., Nance, J. D., Ray, E. A., Steenrod, S. D., Strahan, S. E., and Sweeney, C.: Tropospheric SF6: Age of air from the Northern Hemisphere midlatitude surface, J. Geophys. Res.-Atmos., 118, 11429–11441,, 2013. a

Wilkening, M.: Rn-222 Concentrations in Convective Patterns of a Mountain Environment, J. Geophys. Res., 75, 1733,, 1970. a

Williams, E.: Lightning and climate: A review, Atmos. Res., 76, 272–287,, 2005. a

World Meteorological Society (WMO): Scientific Assessment of Ozone Depletion: 2018, Global Ozone Research and Monitoring Project – Report No. 58, Tech. rep.,, 2018. a

WMO/GAW Ozone Monitoring Community, World Meteorological Organization-Global Atmosphere Watch Program (WMO-GAW)/World Ozone and Ultraviolet Radiation Data Centre (WOUDC) [data set], retrieved from: (last access: 4 November 2019), a list of all contributors is available on the website,, 2019. a

Wu, S., Mickley, L. J., Jacob, D. J., Logan, J. A., Yantosca, R. M., and Rind, D.: Why are there large differences between models in global budgets of tropospheric ozone, J. Geophys. Res., 112, D05302,, 2007.  a, b, c

Wu, S., Mickley, L. J., Jacob, D. J., Rind, D., and Streets, D. G.: Effects of 2000–2050 changes in climate and emissions on global tropospheric ozone and the policy-relevant background surface ozone in the United States, J. Geophys. Res., 113, D18312,, 2008a. a

Wu, S., Mickley, L. J., Leibensperger, E. M., Jacob, D. J., Rind, D., and Streets, D. G.: Effects of 2000–2050 global change on ozone air quality in the United States, J. Geophys. Res., 113, D06302,, 2008b. a

Yan, Y., Lin, J., Chen, J., and Hu, L.: Improved simulation of tropospheric ozone by a global-multi-regional two-way coupling model system, Atmos. Chem. Phys., 16, 2381–2400,, 2016. a

Young, P. J., Naik, V., Fiore, A. M., Gaudel, A., Guo, J., Lin, M. Y., Neu, J. L., Parrish, D. D., Rieder, H. E., Schnell, J. L., Tilmes, S., Wild, O., Zhang, L., Ziemke, J. R., Brandt, J., Delcloo, A., Doherty, R. M., Geels, C., Hegglin, M. I., Hu, L., Im, U., Kumar, R., Luhar, A., Murray, L. T., Plummer, D., Rodriguez, J., Saiz-Lopez, A., Schultz, M. G., Woodhouse, M. T., and Zeng, G.: Tropospheric Ozone Assessment Report: Assessment of global-scale model performance for global and regional ozone distributions, variability, and trends, Elem. Sci. Anth., 6, 10,, 2018. a

Zender, C. S., Bian, H., and Newman, D.: Mineral Dust Entrainment and Deposition (DEAD) model: Description and 1990s dust climatology, J. Geophys. Res., 108, 4416,, 2003. a

Zhu, J., Liao, H., Mao, Y., Yang, Y., and Jiang, H.: Interannual variation, decadal trend, and future change in ozone outflow from East Asia, Atmos. Chem. Phys., 17, 3729–3747,, 2017. a

Ziemke, J. R., Chandra, S., Duncan, B. N., Froidevaux, L., Bhartia, P. K., Levelt, P. F., and Waters, J. W.: Tropospheric ozone determined from aura OMI and MLS: Evaluation of measurements and comparison with the Global Modeling Initiative's Chemical Transport Model, J. Geophys. Res.-Atmos., 111, D19303,, 2006. a, b


Like many free-running climate models, E2.1 uses a 365 d calendar, whereas GEOS-Chem includes leap days; the default behavior of GCAP 2.0 is to repeat 28 February meteorology on 29 February. Users alternatively may stop the model at the end of 28 February and apply the restart file to 1 March for leap years to avoid meteorological discontinuities.


Note that the native horizontal grid of E2.1 is offset from that traditionally used by GEOS-Chem at comparable resolutions. The former has the International Date Line as a cell edge, whereas the latter has it as a cell midpoint. In addition, E2.1 does not make use of half-polar cells as does GEOS-DAS or GEOS-Chem, making the total number of latitude bands one fewer in E2.1 as opposed to GEOS-Chem.

Short summary
Chemical-transport models are tools used to study air pollution and inform public policy. However, they are limited by the availability of archived meteorology. Here, we describe how the GEOS-Chem chemical-transport model may now be driven by meteorology archived from a state-of-the-art general circulation model for past and future climates, allowing it to be used to explore the impact of climate change on air pollution and atmospheric composition.