Search for: All records

Abstract. Tropospheric ozone is a major air pollutant and greenhouse gas. It is also the primary precursor of OH, the main tropospheric oxidant. Global atmospheric chemistry models show large differences in their simulations of tropospheric ozone budgets. Here we implement the widely used GEOS-Chem atmospheric chemistry module as an alternative to CAM-chem within the Community Earth System Model version 2 (CESM2). We compare the resulting GEOS-Chem and CAM-chem simulations of tropospheric ozone and related species within CESM2 to observations from ozonesondes, surface sites, the ATom-1 aircraft campaign over the Pacific and Atlantic, and the KORUS-AQ aircraft campaign over the Seoul Metropolitan Area. We find that GEOS-Chem and CAM-chem within CESM2 have similar tropospheric ozone budgets and concentrations usually within 5 ppb but important differences in the underlying processes including (1) photolysis scheme (no aerosol effects in CAM-chem), (2) aerosol nitrate photolysis, (3) N2O5 cloud uptake, (4) tropospheric halogen chemistry, and (5) ozone deposition to the oceans. Global tropospheric OH concentrations are the same in both models, but there are large regional differences reflecting the above processes. Carbon monoxide is lower in CAM-chem (and lower than observations), at least in part because of higher OH concentrations in the Northern Hemisphere and insufficient production from isoprene oxidation in the Southern Hemisphere. CESM2 does not scavenge water-soluble gases in convective updrafts, leading to some upper-tropospheric biases. Comparison to KORUS-AQ observations shows an overestimate of ozone above 4 km altitude in both models, which at least in GEOS-Chem is due to inadequate scavenging of particulate nitrate in convective updrafts in CESM2, leading to excessive NO production from nitrate photolysis. The KORUS-AQ comparison also suggests insufficient boundary layer mixing in CESM2. This implementation and evaluation of GEOS-Chem in CESM2 contribute to the MUSICA vision of modularizing tropospheric chemistry in Earth system models.

 

Abstract. We quantify future changes in wildfire burned area and carbon emissions inthe 21st century under four Shared Socioeconomic Pathways (SSPs) scenariosand two SSP5-8.5-based solar geoengineering scenarios with a target surfacetemperature defined by SSP2-4.5 – solar irradiance reduction (G6solar) andstratospheric sulfate aerosol injections (G6sulfur) – and explore themechanisms that drive solar geoengineering impacts on fires. This study isbased on fully coupled climate–chemistry simulations with simulatedoccurrence of fires (burned area and carbon emissions) using the WholeAtmosphere Community Climate Model version 6 (WACCM6) as the atmosphericcomponent of the Community Earth System Model version 2 (CESM2). Globally,total wildfire burned area is projected to increase over the 21st centuryunder scenarios without geoengineering and decrease under the twogeoengineering scenarios. By the end of the century, the two geoengineeringscenarios have lower burned area and fire carbon emissions than not onlytheir base-climate scenario SSP5-8.5 but also the targeted-climate scenarioSSP2-4.5. Geoengineering reduces wildfire occurrence by decreasing surfacetemperature and wind speed and increasing relative humidity and soil water,with the exception of boreal regions where geoengineering increases theoccurrence of wildfires due to a decrease in relative humidity and soilwater compared with the present day. This leads to a global reduction in burnedarea and fire carbon emissions by the end of the century relative to theirbase-climate scenario SSP5-8.5. However, geoengineering also yieldsreductions in precipitation compared with a warming climate, which offsetssome of the fire reduction. Overall, the impacts of the different drivingfactors are larger on burned area than fire carbon emissions. In general,the stratospheric sulfate aerosol approach has a stronger fire-reducingeffect than the solar irradiance reduction approach.

 

Abstract. We present the Fire Inventory from National Center for Atmospheric Research (NCAR) version 2.5 (FINNv2.5), a fire emissions inventory that provides publicly available emissions of trace gases and aerosols for various applications, including use in global and regional atmospheric chemistry modeling. FINNv2.5 includes numerous updates to the FINN version 1 framework to better represent burned area, vegetation burned, and chemicals emitted. Major changes include the use of active fire detections from the Visible Infrared Imaging Radiometer Suite (VIIRS) at 375 m spatial resolution, which allows smaller fires to be included in the emissions processing. The calculation of burned area has been updated such that a more rigorous approach is used to aggregate fire detections, which better accounts for larger fires and enables using multiple satellite products simultaneously for emissions estimates. Fuel characterization and emissions factors have also been updated in FINNv2.5. Daily fire emissions for many trace gases and aerosols are determined for 2002–2019 (Moderate Resolution Imaging Spectroradiometer (MODIS)-only fire detections) and 2012–2019 (MODIS + VIIRS fire detections). The non-methane organic gas emissions are allocated to the species of several commonly used chemical mechanisms. We compare FINNv2.5 emissions against other widely used fire emissions inventories. The performance of FINNv2.5 emissions as inputs to a chemical transport model is assessed with satellite observations. Uncertainties in the emissions estimates remain, particularly in Africa and South America during August–October and in southeast and equatorial Asia in March and April. Recommendations for future evaluation and use are given.

 

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. 

Abstract. Emissions are a central component of atmosphericchemistry models. The Harmonized Emissions Component (HEMCO) is a softwarecomponent for computing emissions from a user-selected ensemble of emissioninventories and algorithms. It allows users to re-grid, combine, overwrite,subset, and scale emissions from different inventories through aconfiguration file and with no change to the model source code. Theconfiguration file also maps emissions to model species with appropriateunits. HEMCO can operate in offline stand-alone mode, but more importantlyit provides an online facility for models to compute emissions at runtime.HEMCO complies with the Earth System Modeling Framework (ESMF) forportability across models. We present a new version here, HEMCO 3.0, thatfeatures an improved three-layer architecture to facilitate implementationinto any atmospheric model and improved capability for calculatingemissions at any model resolution including multiscale and unstructuredgrids. The three-layer architecture of HEMCO 3.0 includes (1) the Data InputLayer that reads the configuration file and accesses the HEMCO library ofemission inventories and other environmental data, (2) the HEMCO Core thatcomputes emissions on the user-selected HEMCO grid, and (3) the ModelInterface Layer that re-grids (if needed) and serves the data to theatmospheric model and also serves model data to the HEMCO Core forcomputing emissions dependent on model state (such as from dust or vegetation). The HEMCO Core is common to the implementation in all models, whilethe Data Input Layer and the Model Interface Layer are adaptable to themodel environment. Default versions of the Data Input Layer and ModelInterface Layer enable straightforward implementation of HEMCO in any simplemodel architecture, and options are available to disable features such asre-gridding that may be done by independent couplers in more complexarchitectures. The HEMCO library of emission inventories and algorithms iscontinuously enriched through user contributions so that new inventoriescan be immediately shared across models. HEMCO can also serve as a generaldata broker for models to process input data not only for emissions but forany gridded environmental datasets. We describe existing implementations ofHEMCO 3.0 in (1) the GEOS-Chem “Classic” chemical transport model withshared-memory infrastructure, (2) the high-performance GEOS-Chem (GCHP)model with distributed-memory architecture, (3) the NASA GEOS Earth SystemModel (GEOS ESM), (4) the Weather Research and Forecasting model withGEOS-Chem (WRF-GC), (5) the Community Earth System Model Version 2 (CESM2),and (6) the NOAA Global Ensemble Forecast System – Aerosols(GEFS-Aerosols), as well as the planned implementation in the NOAA Unified ForecastSystem (UFS). Implementation of HEMCO in CESM2 contributes to theMulti-Scale Infrastructure for Chemistry and Aerosols (MUSICA) by providinga common emissions infrastructure to support different simulations ofatmospheric chemistry across scales. 

Abstract. Secondary organic aerosol (SOA) is a dominant contributor of fine particulate matter in the atmosphere, but the complexity of SOA formation chemistry hinders the accurate representation of SOA in models. Volatility-based SOA parameterizations have been adopted in many recent chemistry modeling studies and have shown a reasonable performance compared to observations. However, assumptions made in these empirical parameterizations can lead to substantial errors when applied to future climatic conditions as they do not include the mechanistic understanding of processes but are rather fitted to laboratory studies of SOA formation. This is particularly the case for SOA derived from isoprene epoxydiols (IEPOX SOA), for which we have a higher level of understanding of the fundamental processes than is currently parameterized in most models. We predict future SOA concentrations using an explicit mechanism and compare the predictions with the empirical parameterization based on the volatility basis set (VBS) approach. We then use the Community Earth System Model 2 (CESM2.1.0) with detailed isoprene chemistry and reactive uptake processes for the middle and end of the 21st century under four Shared Socioeconomic Pathways (SSPs): SSP1–2.6, SSP2–4.5, SSP3–7.0, and SSP5–8.5. With the explicit chemical mechanism, we find that IEPOX SOA is predicted to increase on average under all future SSP scenarios but with some variability in the results depending on regions and the scenario chosen. Isoprene emissions are the main driver of IEPOX SOA changes in the future climate, but the IEPOX SOA yield from isoprene emissions also changes by up to 50 % depending on the SSP scenario, in particular due to different sulfur emissions. We conduct sensitivity simulations with and without CO2 inhibition of isoprene emissions that is highly uncertain, which results in factor of 2 differences in the predicted IEPOX SOA global burden, especially for the high-CO2 scenarios (SSP3–7.0 and SSP5–8.5). Aerosol pH also plays a critical role in the IEPOX SOA formation rate, requiring accurate calculation of aerosol pH in chemistry models. On the other hand, isoprene SOA calculated with the VBS scheme predicts a nearly constant SOA yield from isoprene emissions across all SSP scenarios; as a result, it mostly follows isoprene emissions regardless of region and scenario. This is because the VBS scheme does not consider heterogeneous chemistry; in other words, there is no dependency on aerosol properties. The discrepancy between the explicit mechanism and VBS parameterization in this study is likely to occur for other SOA components as well, which may also have dependencies that cannot be captured by VBS parameterizations. This study highlights the need for more explicit chemistry or for parameterizations that capture the dependence on key physicochemical drivers when predicting SOA concentrations for climate studies. 

Abstract. Secondary organic aerosol derived from isopreneepoxydiols (IEPOX-SOA) is thought to contribute the dominant fraction oftotal isoprene SOA, but the current volatility-based lumped SOAparameterizations are not appropriate to represent the reactive uptake ofIEPOX onto acidified aerosols. A full explicit modeling of this chemistryis however computationally expensive owing to the many species and reactionstracked, which makes it difficult to include it in chemistry–climate modelsfor long-term studies. Here we present three simplified parameterizations(version 1.0) for IEPOX-SOA simulation, based on an approximateanalytical/fitting solution of the IEPOX-SOA yield and formation timescale.The yield and timescale can then be directly calculated using the globalmodel fields of oxidants, NO, aerosol pH and other key properties, and drydeposition rates. The advantage of the proposed parameterizations is thatthey do not require the simulation of the intermediates while retaining thekey physicochemical dependencies. We have implemented the newparameterizations into the GEOS-Chem v11-02-rc chemical transport model,which has two empirical treatments for isoprene SOA (the volatility-basis-set, VBS, approach and a fixed 3 % yield parameterization), and comparedall of them to the case with detailed fully explicit chemistry. The bestparameterization (PAR3) captures the global tropospheric burden of IEPOX-SOAand its spatiotemporal distribution (R²=0.94) vs. thosesimulated by the full chemistry, while being more computationally efficient(∼5 times faster), and accurately captures the response tochanges in NO_x and SO₂ emissions. On the other hand, the constant3 % yield that is now the default in GEOS-Chem deviates strongly (R²=0.66), as does the VBS (R²=0.47, 49 % underestimation), withneither parameterization capturing the response to emission changes. Withthe advent of new mass spectrometry instrumentation, many detailed SOAmechanisms are being developed, which will challenge global and especiallyclimate models with their computational cost. The methods developed in thisstudy can be applied to other SOA pathways, which can allow includingaccurate SOA simulations in climate and global modeling studies in thefuture.

 

ABSTRACT To explore the various couplings across space and time and between ecosystems in a consistent manner, atmospheric modeling is moving away from the fractured limited-scale modeling strategy of the past toward a unification of the range of scales inherent in the Earth system. This paper describes the forward-looking Multi-Scale Infrastructure for Chemistry and Aerosols (MUSICA), which is intended to become the next-generation community infrastructure for research involving atmospheric chemistry and aerosols. MUSICA will be developed collaboratively by the National Center for Atmospheric Research (NCAR) and university and government researchers, with the goal of serving the international research and applications communities. The capability of unifying various spatiotemporal scales, coupling to other Earth system components, and process-level modularization will allow advances in both fundamental and applied research in atmospheric composition, air quality, and climate and is also envisioned to become a platform that addresses the needs of policy makers and stakeholders. 

Model intercomparison studies often report a large spread in simulation results, but quantifying the causes of these differences is hindered by the fact that several processes contribute to the model spread simultaneously. Here we use the Multi‐Scale Infrastructure for Chemistry and Aerosols (MUSICA) version 0 to investigate the model resolution dependencies of simulated chemical species, with a focus on the differences between global uniform grid and regional refinement grid simulations with the same modeling framework. We construct two global (ne30 [∼112 km] and ne60 [∼56 km]) and two regional refinement grids over Korea (ne30x8 [∼14 km] and ne30x16 [∼7 km]). The grid resolution can change chemical concentrations by an order of magnitude in the boundary layer, and the importance increases as the species' reactivity increases (e.g., up to 50% and 1,000% changes for ethane and xylenes, respectively). The diurnal cycle of oxidants (OH, O₃, and NO₃) also varies with the grid resolution, which leads to different oxidation pathways of volatile organic compounds (e.g., the fraction of monoterpenes reacting with NO₃in Seoul around midnight is 90% for ne30, but 65% for ne30x16). The models with high‐resolution grids usually do a better job at reproducing aircraft observations during the KORUS‐AQ campaign, but not always, implying compensating errors in the coarse grid simulations. For example, ozone is better reproduced by the coarse grid due to the artificial mixing of NO_x. When developing new chemical mechanisms and evaluating models over urban areas, the uncertainties associated with model resolution should be considered.