skip to main content


Title: Implementation and evaluation of the GEOS-Chem chemistry module version 13.1.2 within the Community Earth System Model v2.1
Abstract. We implement the GEOS-Chem chemistry module as a chemical mechanism in version 2 of the Community Earth System Model (CESM). Our implementation allowsthe state-of-the-science GEOS-Chem chemistry module to be used with identical emissions, meteorology, and climate feedbacks as the CAM-chemchemistry module within CESM. We use coupling interfaces to allow GEOS-Chem to operate almost unchanged within CESM. Aerosols are converted at eachtime step between the GEOS-Chem bulk representation and the size-resolved representation of CESM's Modal Aerosol Model (MAM4). Land-type informationneeded for dry-deposition calculations in GEOS-Chem is communicated through a coupler, allowing online land–atmosphere interactions. Wet scavengingin GEOS-Chem is replaced with the Neu and Prather scheme, and a common emissions approach is developed for both CAM-chem and GEOS-Chem in CESM. We compare how GEOS-Chem embedded in CESM (C-GC) compares to the existing CAM-chem chemistry option (C-CC) when used to simulate atmosphericchemistry in 2016, with identical meteorology and emissions. We compare the atmospheric composition and deposition tendencies between the twosimulations and evaluate the residual differences between C-GC and its use as a stand-alone chemistry transport model in the GEOS-Chem HighPerformance configuration (S-GC). We find that stratospheric ozone agrees well between the three models, with differences of less than 10 % inthe core of the ozone layer, but that ozone in the troposphere is generally lower in C-GC than in either C-CC or S-GC. This is likely due to greatertropospheric concentrations of bromine, although other factors such as water vapor may contribute to lesser or greater extents depending on theregion. This difference in tropospheric ozone is not uniform, with tropospheric ozone in C-GC being 30 % lower in the Southern Hemisphere whencompared with S-GC but within 10 % in the Northern Hemisphere. This suggests differences in the effects of anthropogenic emissions. Aerosolconcentrations in C-GC agree with those in S-GC at low altitudes in the tropics but are over 100 % greater in the upper troposphere due todifferences in the representation of convective scavenging. We also find that water vapor concentrations vary substantially between the stand-aloneand CESM-implemented version of GEOS-Chem, as the simulated hydrological cycle in CESM diverges from that represented in the source NASA Modern-Era Retrospective analysis for Research and Applications (Version 2; MERRA-2)reanalysis meteorology which is used directly in the GEOS-Chem chemistrytransport model (CTM). Our implementation of GEOS-Chem as a chemistry option in CESM (including full chemistry–climate feedback) is publicly available and is beingconsidered for inclusion in the CESM main code repository. This work is a significant step in the MUlti-Scale Infrastructure for Chemistry andAerosols (MUSICA) project, enabling two communities of atmospheric researchers (CESM and GEOS-Chem) to share expertise through a common modelingframework, thereby accelerating progress in atmospheric science.  more » « less
Award ID(s):
1914903 1914920
NSF-PAR ID:
10417639
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Geoscientific Model Development
Volume:
15
Issue:
23
ISSN:
1991-9603
Page Range / eLocation ID:
8669 to 8704
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Bromine released from the decomposition of short-lived brominated source gases contributes as a sink of ozone in the lower stratosphere.The two major contributors are CH2Br2 and CHBr3.In this study, we investigate the global seasonal distribution of these two substances, based on four High Altitude and Long Range Research Aircraft (HALO) missions, the HIAPER Pole-to-Pole Observations (HIPPO) mission, and the Atmospheric Tomography (ATom) mission.Observations of CH2Br2 in the free and upper troposphere indicate a pronounced seasonality in both hemispheres, with slightly larger mixing ratios in the Northern Hemisphere (NH).Compared to CH2Br2, CHBr3 in these regions shows larger variability and less clear seasonality, presenting larger mixing ratios in winter and autumn in NH midlatitudes to high latitudes.The lowermost stratosphere of SH and NH shows a very similar distribution of CH2Br2 in hemispheric spring with differences well below 0.1 ppt, while the differences in hemispheric autumn are much larger with substantially smaller values in the SH than in the NH.This suggests that transport processes may be different in both hemispheric autumn seasons, which implies that the influx of tropospheric air (“flushing”) into the NH lowermost stratosphere is more efficient than in the SH.The observations of CHBr3 support the suggestion, with a steeper vertical gradient in the upper troposphere and lower stratosphere in SH autumn than in NH autumn.However, the SH database is insufficient to quantify this difference.We further compare the observations to model estimates of TOMCAT (Toulouse Off-line Model of Chemistry And Transport) and CAM-Chem (Community Atmosphere Model with Chemistry, version 4), both using the same emission inventory of Ordóñez et al. (2012).The pronounced tropospheric seasonality of CH2Br2 in the SH is not reproduced by the models,presumably due to erroneous seasonal emissions or atmospheric photochemical decomposition efficiencies.In contrast, model simulations of CHBr3 show a pronounced seasonality in both hemispheres, which is not confirmed by observations.The distributions of both species in the lowermost stratosphere of the Northern and Southern hemispheres are overall well captured by the models with the exception of southern hemispheric autumn,where both models present a bias that maximizes in the lowest 40 K above the tropopause, with considerably lower mixing ratios in the observations.Thus, both models reproduce equivalent flushing in both hemispheres, which is not confirmed by the limited available observations.Our study emphasizes the need for more extensive observations in the SH to fully understand the impact of CH2Br2 and CHBr3 on lowermost-stratospheric ozone loss and to help constrain emissions.

     
    more » « less
  2. Abstract. In this paper, we present a new version of the chemistry–climate model SOCOL-AERv2 supplemented by an iodine chemistry module. We perform three 20-year ensemble experiments to assess the validity of the modeled iodine and to quantify the effects of iodine on ozone. The iodine distributions obtained with SOCOL-AERv2-I agree well with AMAX-DOAS observations and with CAM-chem model simulations. For the present-day atmosphere, the model suggests that the iodine-induced chemistry leads to a 3 %–4 % reduction in the ozone column, which is greatest at high latitudes. The model indicates the strongest influence of iodine in the lower stratosphere with 30 ppbv less ozone at low latitudes and up to 100 ppbv less at high latitudes. In the troposphere, the account of the iodine chemistry reduces the tropospheric ozone concentration by 5 %–10 % depending on geographical location. In the lower troposphere, 75 % of the modeled ozone reduction originates from inorganic sources of iodine, 25 % from organic sources of iodine. At 50 hPa, the results show that the impacts of iodine from both sources are comparable. Finally, we determine the sensitivity of ozone to iodine by applying a 2-fold increase in iodine emissions, as it might be representative for iodine by the end of this century. This reduces the ozone column globally by an additional 1.5 %–2.5 %. Our results demonstrate the sensitivity of atmospheric ozone to iodine chemistry for present and future conditions, but uncertainties remain high due to the paucity of observational data of iodine species. 
    more » « less
  3. Abstract. We present an updated mechanism for tropospheric halogen (Cl + Br + I) chemistry in the GEOS-Chem global atmospheric chemical transportmodel and apply it to investigate halogen radical cycling and implications for tropospheric oxidants. Improved representation of HOBr heterogeneouschemistry and its pH dependence in our simulation leads to less efficient recycling and mobilization of bromine radicals and enables the model toinclude mechanistic sea salt aerosol debromination without generating excessive BrO. The resulting global mean tropospheric BrO mixingratio is 0.19 ppt (parts per trillion), lower than previous versions of GEOS-Chem. Model BrO shows variable consistency and biases in comparison tosurface and aircraft observations in marine air, which are often near or below the detection limit. The model underestimates the daytimemeasurements of Cl2 and BrCl from the ATom aircraft campaign over the Pacific and Atlantic, which if correct would imply a very largemissing primary source of chlorine radicals. Model IO is highest in the marine boundary layer and uniform in the free troposphere, with a globalmean tropospheric mixing ratio of 0.08 ppt, and shows consistency with surface and aircraft observations. The modeled global meantropospheric concentration of Cl atoms is 630 cm−3, contributing 0.8 % of the global oxidation of methane, 14 % of ethane,8 % of propane, and 7 % of higher alkanes. Halogen chemistry decreases the global tropospheric burden of ozone by 11 %,NOx by 6 %, and OH by 4 %. Most of the ozone decrease is driven by iodine-catalyzed loss. The resulting GEOS-Chem ozonesimulation is unbiased in the Southern Hemisphere but too low in the Northern Hemisphere. 
    more » « less
  4. Abstract. Mercury (Hg) is emitted to the atmosphere mainly as volatile elemental Hg0. Oxidation to water-soluble HgII plays a major role in Hg deposition to ecosystems. Here, we implement a new mechanism for atmospheric Hg0HgII redox chemistry in the GEOS-Chem global model and examine the implications for the global atmospheric Hg budget and deposition patterns. Our simulation includes a new coupling of GEOS-Chem to an ocean general circulation model (MITgcm), enabling a global 3-D representation of atmosphere–ocean Hg0HgII cycling. We find that atomic bromine (Br) of marine organobromine origin is the main atmospheric Hg0 oxidant and that second-stage HgBr oxidation is mainly by the NO2 and HO2 radicals. The resulting chemical lifetime of tropospheric Hg0 against oxidation is 2.7 months, shorter than in previous models. Fast HgII atmospheric reduction must occur in order to match the  ∼ 6-month lifetime of Hg against deposition implied by the observed atmospheric variability of total gaseous mercury (TGM ≡ Hg0+HgII(g)). We implement this reduction in GEOS-Chem as photolysis of aqueous-phase HgII–organic complexes in aerosols and clouds, resulting in a TGM lifetime of 5.2 months against deposition and matching both mean observed TGM and its variability. Model sensitivity analysis shows that the interhemispheric gradient of TGM, previously used to infer a longer Hg lifetime against deposition, is misleading because Southern Hemisphere Hg mainly originates from oceanic emissions rather than transport from the Northern Hemisphere. The model reproduces the observed seasonal TGM variation at northern midlatitudes (maximum in February, minimum in September) driven by chemistry and oceanic evasion, but it does not reproduce the lack of seasonality observed at southern hemispheric marine sites. Aircraft observations in the lowermost stratosphere show a strong TGM–ozone relationship indicative of fast Hg0 oxidation, but we show that this relationship provides only a weak test of Hg chemistry because it is also influenced by mixing. The model reproduces observed Hg wet deposition fluxes over North America, Europe, and China with little bias (0–30%). It reproduces qualitatively the observed maximum in US deposition around the Gulf of Mexico, reflecting a combination of deep convection and availability of NO2 and HO2 radicals for second-stage HgBr oxidation. However, the magnitude of this maximum is underestimated. The relatively low observed Hg wet deposition over rural China is attributed to fast HgII reduction in the presence of high organic aerosol concentrations. We find that 80% of HgII deposition is to the global oceans, reflecting the marine origin of Br and low concentrations of organic aerosols for HgII reduction. Most of that deposition takes place to the tropical oceans due to the availability of HO2 and NO2 for second-stage HgBr oxidation.

     
    more » « less
  5. Abstract

    Many Chemistry‐Climate Models (CCMs) include a simplified treatment of brominated very short‐lived (VSLBr) species by assuming CH3Br as a surrogate for VSLBr. However, neglecting a comprehensive treatment of VSLBrin CCMs may yield an unrealistic representation of the associated impacts. Here, we use the Community Atmospheric Model with Chemistry (CAM‐Chem) CCM to quantify the tropospheric and stratospheric changes between various VSLBrchemical approaches with increasing degrees of complexity (i.e., surrogate, explicit, and full). Our CAM‐Chem results highlight the improved accuracy achieved by considering a detailed treatment of VSLBrphotochemistry, including sea‐salt aerosol dehalogenation and heterogeneous recycling on ice‐crystals. Differences between the full and surrogate schemes maximize in the lowermost stratosphere and midlatitude free troposphere, resulting in a latitudinally dependent reduction of ∼1–7 DU in total ozone column and a ∼5%–15% decrease of the OH/HO2ratio. We encourage all CCMs to include a complete chemical treatment of VSLBrin the troposphere and stratosphere.

     
    more » « less