scholarly journals Evaluation of the interactive stratospheric ozone (O3v2 module) for the E3SM version 2 Earth System Model

2020 ◽  
Author(s):  
Qi Tang ◽  
Michael J. Prather ◽  
Juno Hsu ◽  
Daniel J. Ruiz ◽  
Philip J. Cameron-Smith ◽  
...  

Abstract. Stratospheric ozone affects climate directly as the predominant heat source in the stratosphere and indirectly through chemical feedbacks controlling other greenhouse gases. The U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1) implemented a new ozone chemistry module that improves the simulation of the sharp tropopause gradients, replacing a version based partly on long-term average climatologies that poorly represented heating rates in the lowermost stratosphere. The new O3v2 module extends seamlessly into the troposphere and preserves the naturally sharp cross-tropopause gradient, with 20-40% less ozone in this region. Additionally, O3v2 enables the diagnosis of stratosphere-troposphere exchange flux of ozone, a key budget term lacking in E3SMv1. Here, we evaluate key features in ozone abundance and other closely related quantities in atmosphere-only E3SMv1 simulations driven by observed sea surface temperatures (SSTs, years 1990-2014), comparing with satellite observations and the University of California, Irvine chemistry transport model (UCI CTM) using the same stratospheric chemistry scheme but driven by European Centre forecast fields for the same period. In terms of stratospheric column ozone, O3v2 shows improved mean bias and northern mid-latitude variability, but not quite as good as the UCI CTM. As expected, SST forcing does not match the observed quasi-biennial oscillation, which is mostly matched with the UCI CTM. This new O3v2 E3SM model retains mostly the same climate state and climate sensitivity as the previous version, and we recommend its use for other climate models that still use ozone climatologies.

2021 ◽  
Vol 14 (3) ◽  
pp. 1219-1236
Author(s):  
Qi Tang ◽  
Michael J. Prather ◽  
Juno Hsu ◽  
Daniel J. Ruiz ◽  
Philip J. Cameron-Smith ◽  
...  

Abstract. Stratospheric ozone affects climate directly as the predominant heat source in the stratosphere and indirectly through chemical reactions controlling other greenhouse gases. The U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1) implemented a new ozone chemistry module that improves the simulation of the sharp tropopause gradients, replacing a version based partly on long-term average climatologies that poorly represented heating rates in the lowermost stratosphere. The new O3v2 module extends seamlessly into the troposphere and preserves the naturally sharp cross-tropopause gradient, with 20 %–40 % less ozone in this region. Additionally, O3v2 enables the diagnosis of stratosphere–troposphere exchange flux of ozone, a key budget term lacking in E3SMv1. Here, we evaluate key features in ozone abundance and other closely related quantities in atmosphere-only E3SMv1 simulations driven by observed sea surface temperatures (SSTs, years 1990–2014), comparing them with satellite observations of ozone and also with the University of California, Irvine chemistry transport model (UCI CTM) using the same stratospheric chemistry scheme but driven by European Centre forecast fields for the same period. In terms of stratospheric column ozone, O3v2 shows reduced mean bias and improved northern midlatitude variability, but it is not quite as good as the UCI CTM. As expected, SST-forced E3SMv1 simulations cannot synchronize with observed quasi-biennial oscillations (QBOs), but they do show the typical QBO pattern seen in column ozone. This new O3v2 E3SMv1 model mostly retains the same climate state and climate sensitivity as the previous version, and we recommend its use for other climate models that still use ozone climatologies.


2021 ◽  
Vol 8 ◽  
Author(s):  
James Keeble ◽  
Yu Yeung Scott Yiu ◽  
Alexander T. Archibald ◽  
Fiona O’Connor ◽  
Alistair Sellar ◽  
...  

Stratospheric ozone projections in the tropics, modeled using the UKESM1 Earth system model, are explored under different Shared Socioeconomic Pathways (SSPs). Consistent with other studies, it is found that tropical stratospheric column ozone does not return to 1980s values by the end of the 21st century under any SSP scenario as increased ozone mixing ratios in the tropical upper stratosphere are offset by continued ozone decreases in the tropical lower stratosphere. Stratospheric column ozone is projected to be largest under SSP scenarios with the smallest change in radiative forcing, and smallest for SSP scenarios with larger radiative forcing, consistent with a faster Brewer-Dobson circulation at high greenhouse gas loadings. This study explores the use of machine learning (ML) techniques to make accurate, computationally inexpensive projections of tropical stratospheric column ozone. Four ML techniques are investigated: Ridge regression, Lasso regression, Random Forests and Extra Trees. All four techniques investigated here are able to make projections of future tropical stratospheric column ozone which agree well with those made by the UKESM1 Earth system model, often falling within the ensemble spread of UKESM1 simulations for a broad range of SSPs. However, all techniques struggle to make accurate projects for the final decades of the SSP5-8.5 scenario. Accurate projections can only be achieved when the ML methods are trained on sufficient data, including both historical and future simulations. When trained only on historical data, the projections made using models based on ML techniques fail to accurately predict tropical stratospheric ozone changes. Results presented here indicate that, when sufficiently trained, ML models have the potential to make accurate, computationally inexpensive projections of tropical stratospheric column ozone. Further development of these models may reduce the computational burden placed on fully coupled chemistry-climate and Earth system models and enable the exploration of tropical stratospheric column ozone recovery under a much broader range of future emissions scenarios.


2012 ◽  
Vol 5 (3) ◽  
pp. 2811-2842 ◽  
Author(s):  
M. A. Chandler ◽  
L. E. Sohl ◽  
J. A. Jonas ◽  
H. J. Dowsett

Abstract. Climate reconstructions of the mid-Pliocene Warm Period (mPWP) bear many similarities to aspects of future global warming as projected by the Intergovernmental Panel on Climate Change. In particular, marine and terrestrial paleoclimate data point to high latitude temperature amplification, with associated decreases in sea ice and land ice and altered vegetation distributions that show expansion of warmer climate biomes into higher latitudes. NASA GISS climate models have been used to study the Pliocene climate since the USGS PRISM project first identified that the mid-Pliocene North Atlantic sea surface temperatures were anomalously warm. Here we present the most recent simulations of the Pliocene using the AR5/CMIP5 version of the GISS Earth System Model known as ModelE2-R. These simulations constitute the NASA contribution to the Pliocene Model Intercomparison Project (PlioMIP) Experiment 2. Many findings presented here corroborate results from other PlioMIP multi-model ensemble papers, but we also emphasize features in the ModelE2-R simulations that are unlike the ensemble means. We provide discussion of features that show considerable improvement compared with simulations from previous versions of the NASA GISS models, improvement defined here as simulation results that more closely resemble the ocean core data as well as the PRISM3D reconstructions of the mid-Pliocene climate. In some regions even qualitative agreement between model results and paleodata are an improvement over past studies, but the dramatic warming in the North Atlantic and Greenland-Iceland-Norwegian Sea in these new simulations is by far the most accurate portrayal ever of this key geographic region by the GISS climate model. Our belief is that continued development of key physical routines in the atmospheric model, along with higher resolution and recent corrections to mixing parameterizations in the ocean model, have led to an Earth System Model that will produce more accurate projections of future climate.


2018 ◽  
Vol 18 (19) ◽  
pp. 14715-14735 ◽  
Author(s):  
Simon Chabrillat ◽  
Corinne Vigouroux ◽  
Yves Christophe ◽  
Andreas Engel ◽  
Quentin Errera ◽  
...  

Abstract. We present a consistent intercomparison of the mean age of air (AoA) according to five modern reanalyses: the European Centre for Medium-Range Weather Forecasts Interim Reanalysis (ERA-Interim), the Japanese Meteorological Agency's Japanese 55-year Reanalysis (JRA-55), the National Centers for Environmental Prediction Climate Forecast System Reanalysis (CFSR) and the National Aeronautics and Space Administration's Modern Era Retrospective analysis for Research and Applications version 1 (MERRA) and version 2 (MERRA-2). The modeling tool is a kinematic transport model driven only by the surface pressure and wind fields. It is validated for ERA-I through a comparison with the AoA computed by another transport model. The five reanalyses deliver AoA which differs in the worst case by 1 year in the tropical lower stratosphere and more than 2 years in the upper stratosphere. At all latitudes and altitudes, MERRA-2 and MERRA provide the oldest values (∼5–6 years in midstratosphere at midlatitudes), while JRA-55 and CFSR provide the youngest values (∼4 years) and ERA-I delivers intermediate results. The spread of AoA at 50 hPa is as large as the spread obtained in a comparison of chemistry–climate models. The differences between tropical and midlatitude AoA are in better agreement except for MERRA-2. Compared with in situ observations, they indicate that the upwelling is too fast in the tropical lower stratosphere. The spread between the five simulations in the northern midlatitudes is as large as the observational uncertainties in a multidecadal time series of balloon observations, i.e., approximately 2 years. No global impact of the Pinatubo eruption can be found in our simulations of AoA, contrary to a recent study which used a diabatic transport model driven by ERA-I and JRA-55 winds and heating rates. The time variations are also analyzed through multiple linear regression analyses taking into account the seasonal cycles, the quasi-biennial oscillation and the linear trends over four time periods. The amplitudes of AoA seasonal variations in the lower stratosphere are significantly larger when using MERRA and MERRA-2 than with the other reanalyses. The linear trends of AoA using ERA-I confirm those found by earlier model studies, especially for the period 2002–2012, where the dipole structure of the latitude–height distribution (positive in the northern midstratosphere and negative in the southern midstratosphere) also matches trends derived from satellite observations of SF6. Yet the linear trends vary substantially depending on the considered period. Over 2002–2015, the ERA-I results still show a dipole structure with positive trends in the Northern Hemisphere reaching up to 0.3 yr dec−1. No reanalysis other than ERA-I finds any dipole structure of AoA trends. The signs of the trends depend strongly on the input reanalysis and on the considered period, with values above 10 hPa varying between approximately −0.4 and 0.4 yr dec−1. Using ERA-I and CFSR, the 2002–2015 trends are negative above 10 hPa, but using the three other reanalyses these trends are positive. Over the whole period (1989–2015) each reanalysis delivers opposite trends; i.e., AoA is mostly increasing with CFSR and ERA-I but mostly decreasing with MERRA, JRA-55 and MERRA-2. In view of this large disagreement, we urge great caution for studies aiming to assess AoA trends derived only from reanalysis winds. We briefly discuss some possible causes for the dependency of AoA on the input reanalysis and highlight the need for complementary intercomparisons using diabatic transport models.


2017 ◽  
Vol 10 (1) ◽  
pp. 271-319 ◽  
Author(s):  
Thomas Gasser ◽  
Philippe Ciais ◽  
Olivier Boucher ◽  
Yann Quilcaille ◽  
Maxime Tortora ◽  
...  

Abstract. This paper provides a comprehensive description of OSCAR v2.2, a simple Earth system model. The general philosophy of development is first explained, followed by a complete description of the model's drivers and various modules. All components of the Earth system necessary to simulate future climate change are represented in the model: the oceanic and terrestrial carbon cycles – including a book-keeping module to endogenously estimate land-use change emissions – so as to simulate the change in atmospheric carbon dioxide; the tropospheric chemistry and the natural wetlands, to simulate that of methane; the stratospheric chemistry, for nitrous oxide; 37 halogenated compounds; changing tropospheric and stratospheric ozone; the direct and indirect effects of aerosols; changes in surface albedo caused by black carbon deposition on snow and land-cover change; and the global and regional response of climate – in terms of temperature and precipitation – to all these climate forcers. Following the probabilistic framework of the model, an ensemble of simulations is made over the historical period (1750–2010). We show that the model performs well in reproducing observed past changes in the Earth system such as increased atmospheric concentration of greenhouse gases or increased global mean surface temperature.


2020 ◽  
Vol 13 (2) ◽  
pp. 717-734 ◽  
Author(s):  
Nicholas A. Davis ◽  
Sean M. Davis ◽  
Robert W. Portmann ◽  
Eric Ray ◽  
Karen H. Rosenlof ◽  
...  

Abstract. Specified dynamics (SD) schemes relax the circulation in climate models toward a reference meteorology to simulate historical variability. These simulations are widely used to isolate the dynamical contributions to variability and trends in trace gas species. However, it is not clear if trends in the stratospheric overturning circulation are properly reproduced by SD schemes. This study assesses numerous SD schemes and modeling choices in the Community Earth System Model (CESM) Whole Atmosphere Chemistry Climate Model (WACCM) to determine a set of best practices for reproducing interannual variability and trends in tropical stratospheric upwelling estimated by reanalyses. Nudging toward the reanalysis meteorology as is typically done in SD simulations does not accurately reproduce lower-stratospheric upwelling trends present in the underlying reanalysis. In contrast, nudging to anomalies from the climatological winds or anomalies from the zonal-mean winds and temperatures better reproduces trends in lower-stratospheric upwelling, possibly because these schemes do not disrupt WACCM's climatology. None of the schemes substantially alter the structure of upwelling trends – instead, they make the trends more or less AMIP-like. An SD scheme's performance in simulating the acceleration of the shallow branch of the mean meridional circulation from 1980 to 2017 hinges on its ability to simulate the downward shift of subtropical lower-stratospheric wave momentum forcing. Key to this is not nudging the zonal-mean temperature field. Gravity wave momentum forcing, which drives a substantial fraction of the upwelling in WACCM, cannot be constrained by nudging and presents an upper limit on the performance of these schemes.


2016 ◽  
Author(s):  
Thomas Gasser ◽  
Philippe Ciais ◽  
Olivier Boucher ◽  
Yann Quilcaille ◽  
Maxime Tortora ◽  
...  

Abstract. This paper provides a comprehensive description of OSCAR v2.2, a simple Earth system model. The general philosophy of development is first explained, it is then followed by a complete description of the model's drivers and various modules. All components of the Earth system necessary to simulate future climate change are represented in the model: the oceanic and terrestrial carbon-cycles – including a book-keeping module to endogenously estimate land-use change emissions – so as to simulate the change in atmospheric carbon dioxide; the tropospheric OH chemistry and the natural wetlands, to simulate that of methane; the stratospheric chemistry, for nitrous oxide; thirty-seven halogenated compounds; changing tropospheric and stratospheric ozone; the direct and indirect effects of aerosols; changes in surface albedo caused by black carbon deposition on snow and land-cover change; and the global and regional response of climate – in terms of temperatures and precipitations – to all these climate forcers. Following the probabilistic framework of the model, an ensemble of simulations is made over the historical period (1750–2010). We show that the model performs well in reproducing observed past changes in the Earth system such as increased atmospheric concentration of greenhouse gases or increased global mean surface temperature.


2020 ◽  
Vol 20 (11) ◽  
pp. 6521-6539 ◽  
Author(s):  
Hans Brenna ◽  
Steffen Kutterolf ◽  
Michael J. Mills ◽  
Kirstin Krüger

Abstract. The supereruption of Los Chocoyos (14.6∘ N, 91.2∘ W) in Guatemala ∼84 kyr ago was one of the largest volcanic events of the past 100 000 years. Recent petrologic data show that the eruption released very large amounts of climate-relevant sulfur and ozone-destroying chlorine and bromine gases (523±94 Mt sulfur, 1200±156 Mt chlorine, and 2±0.46 Mt bromine). Using the Earth system model (ESM) of the Community Earth System Model version 2 (CESM2) coupled with the Whole Atmosphere Community Climate Model version 6 (WACCM6), we simulated the impacts of the sulfur- and halogen-rich Los Chocoyos eruption on the preindustrial Earth system. Our simulations show that elevated sulfate burden and aerosol optical depth (AOD) persists for 5 years in the model, while the volcanic halogens stay elevated for nearly 15 years. As a consequence, the eruption leads to a collapse of the ozone layer with global mean column ozone values dropping to 50 DU (80 % decrease) and leading to a 550 % increase in surface UV over the first 5 years, with potential impacts on the biosphere. The volcanic eruption shows an asymmetric-hemispheric response with enhanced aerosol, ozone, UV, and climate signals over the Northern Hemisphere. Surface climate is impacted globally due to peak AOD of >6, which leads to a maximum surface cooling of >6 K, precipitation and terrestrial net primary production decrease of >25 %, and sea ice area increases of 40 % in the first 3 years. Locally, a wetting (>100 %) and strong increase in net primary production (NPP) (>700 %) over northern Africa is simulated in the first 5 years and related to a southward shift of the Intertropical Convergence Zone (ITCZ) to the southern tropics. The ocean responds with pronounced El Niño conditions in the first 3 years that shift to the southern tropics and are coherent with the ITCZ change. Recovery to pre-eruption ozone levels and climate takes 15 years and 30 years, respectively. The long-lasting surface cooling is sustained by an immediate increase in the Arctic sea ice area, followed by a decrease in poleward ocean heat transport at 60∘ N which lasts up to 20 years. In contrast, when simulating Los Chocoyos conventionally by including sulfur and neglecting halogens, we simulate a larger sulfate burden and AOD, more pronounced surface climate changes, and an increase in column ozone. By comparing our aerosol chemistry ESM results to other supereruption simulations with aerosol climate models, we find a higher surface climate impact per injected sulfur amount than previous studies for our different sets of model experiments, since the CESM2(WACCM6) creates smaller aerosols with a longer lifetime, partly due to the interactive aerosol chemistry. As the model uncertainties for the climate response to supereruptions are very large, observational evidence from paleo archives and a coordinated model intercomparison would help to improve our understanding of the climate and environment response.


2018 ◽  
Vol 9 (3) ◽  
pp. 1045-1062 ◽  
Author(s):  
Andrés Navarro ◽  
Raúl Moreno ◽  
Francisco J. Tapiador

Abstract. ESMs (Earth system models) are important tools that help scientists understand the complexities of the Earth's climate. Advances in computing power have permitted the development of increasingly complex ESMs and the introduction of better, more accurate parameterizations of processes that are too complex to be described in detail. One of the least well-controlled parameterizations involves human activities and their direct impact at local and regional scales. In order to improve the direct representation of human activities and climate, we have developed a simple, scalable approach that we have named the POPEM module (POpulation Parameterization for Earth Models). This module computes monthly fossil fuel emissions at grid-point scale using the modeled population projections. This paper shows how integrating POPEM parameterization into the CESM (Community Earth System Model) enhances the realism of global climate modeling, improving this beyond simpler approaches. The results show that it is indeed advantageous to model CO2 emissions and pollutants directly at model grid points rather than using the same mean value globally. A major bonus of this approach is the increased capacity to understand the potential effects of localized pollutant emissions on long-term global climate statistics, thus assisting adaptation and mitigation policies.


2021 ◽  
Vol 31 (06) ◽  
pp. 2130017
Author(s):  
Thomas E. Mulder ◽  
Heiko Goelzer ◽  
Fred W. Wubs ◽  
Henk A. Dijkstra

There is now much geological evidence that the Earth was fully glaciated during several periods in the geological past (about 700[Formula: see text]Myr ago) and attained a so-called Snowball Earth (SBE) state. Additional support for this idea has come from climate models of varying complexity that show transitions to SBE states and undergo hysteresis under changes in solar radiation. In this paper, we apply large-scale bifurcation analyses to a novel, fully-implicit Earth System Model of Intermediate Complexity (I-EMIC) to study SBE transitions. The I-EMIC contains a primitive equation ocean model, a model for atmospheric heat and moisture transport, a sea ice component and formulations for the adjustment of albedo over snow and ice. With the I-EMIC, high-dimensional branches of the SBE bifurcation diagram are obtained through parameter continuation. We are able to identify stable and unstable equilibria and uncover an intricate bifurcation structure associated with the ice-albedo feedback. Moreover, large-scale linear stability analyses are performed near major bifurcations, revealing the spatial nature of destabilizing perturbations.


Sign in / Sign up

Export Citation Format

Share Document