Higher Vertical Resolution for Select Physical Processes in the Energy Exascale Earth System Model (E3SM)

Author(s):  
Hsiang-He Lee ◽  
Peter Bogenschutz ◽  
Takanobu Yamaguchi

<p>The low cloud bias in atmospheric models for climate and weather remains an unsolved problem. Coarse vertical resolution in the current global climate models (GCM) may be a significant cause of low cloud bias because planetary boundary layer (PBL) parameterizations, including higher-order turbulence closure (HOC), cannot resolve sharp temperature and moisture gradients often found at the top of subtropical stratocumulus layers. The aim of this work is to implement a new computational method, the Framework for Improvement by Vertical Enhancement (FIVE) into the Energy Exascale Earth System Model (E3SM) and its single column model. Three physics schemes are interfaced to vertically enhanced physics (VEP), which allows for these schemes to be computed on a higher vertical resolution grid compared to rest of the E3SM model.  In this presentation we use VEP for turbulence, microphysics, and radiation parameterizations and demonstrate better representation of subtropical boundary layer clouds while limiting additional computational cost from the increased number of levels.  We will also briefly discuss future plans for an adaptive vertical grid for VEP, which will allow for additional layers to be added only when/where they are needed.</p>

2020 ◽  
Author(s):  
Chih-Chieh Chen ◽  
Changhai Liu ◽  
Mitch Moncrieff ◽  
Yaga Richter

<p>The importance of convective organization on the global circulation has been recognized for a long time, but parameterizations of the associated processes are missing in global climate models. Contemporary convective parameterizations commonly use a convective plume model (or a spectrum of plumes). This is perhaps appropriate for unorganized convection but the assumption of a gap between the small cumulus scale and the large-scale motion fails to recognize mesoscale dynamics manifested in mesoscale convective systems (MCSs) and multi-scale cloud systems associated with the MJO. Organized convection is abundant in environments featuring vertical wind shear, and significantly modulates the life cycle of moist convection, the transport of heat and momentum, and accounts for a large percentage of precipitation in the tropics. Mesoscale convective organization is typically associated with counter-gradient momentum transport, and distinct heating profiles between the convective and stratiform regions.</p><p>Moncrieff, Liu and Bogenschutz (2017) recently developed a dynamical based parameterization of organized moisture convection, referred to as multiscale coherent structure parameterization (MCSP), for global climate models. A prototype version of MCSP has been implemented in the NCAR Community Earth System Model (CESM) and the Energy Exascale Earth System Model (E3SM), positively affecting the distribution of tropical precipitation, convectively coupled tropical waves, and the Madden-Julian oscillation. We will show the further development of the MCSP and its impact on the simulation of mean precipitation and variability in the two global climate models.</p>


2018 ◽  
Vol 2018 ◽  
pp. 1-24 ◽  
Author(s):  
Jacob Agyekum ◽  
Thompson Annor ◽  
Benjamin Lamptey ◽  
Emmannuel Quansah ◽  
Richard Yao Kuma Agyeman

A selected number of global climate models (GCMs) from the fifth Coupled Model Intercomparison Project (CMIP5) were evaluated over the Volta Basin for precipitation. Biases in models were computed by taking the differences between the averages over the period (1950–2004) of the models and the observation, normalized by the average of the observed for the annual and seasonal timescales. The Community Earth System Model, version 1-Biogeochemistry (CESM1-BGC), the Community Climate System Model Version 4 (CCSM4), the Max Planck Institute Earth System Model, Medium Range (MPI-ESM-MR), the Norwegian Earth System Model (NorESM1-M), and the multimodel ensemble mean were able to simulate the observed climatological mean of the annual total precipitation well (average biases of 1.9% to 7.5%) and hence were selected for the seasonal and monthly timescales. Overall, all the models (CESM1-BGC, CCSM4, MPI-ESM-MR, and NorESM1-M) scored relatively low for correlation (<0.5) but simulated the observed temporal variability differently ranging from 1.0 to 3.0 for the seasonal total. For the annual cycle of the monthly total, the CESM1-BGC, the MPI-ESM-MR, and the NorESM1-M were able to simulate the peak of the observed rainy season well in the Soudano-Sahel, the Sahel, and the entire basin, respectively, while all the models had difficulty in simulating the bimodal pattern of the Guinea Coast. The ensemble mean shows high performance compared to the individual models in various timescales.


2020 ◽  
Author(s):  
Charlotte Lang ◽  
Charles Amory ◽  
Alison Delhasse ◽  
Stefan Hofer ◽  
Christoph Kittel ◽  
...  

&lt;p&gt;We have compared the surface mass (SMB) and energy balance of the Earth System model (ESM) CESM (Community Earth System Model) with those of the regional climate model (RCM) MAR (Mod&amp;#232;le Atmosph&amp;#233;rique R&amp;#233;gional) forced by CESM over the present era (1981 &amp;#8212; 2010) and the future (2011 &amp;#8212; 2100 with SSP585 scenario).&lt;/p&gt;&lt;p&gt;Until now, global climate models (GCM) and ESMs forcing RCMs such as MAR didn&amp;#8217;t include a module able to simulate snow and energy balance at the surface of a snow pack like the SISVAT module of MAR and were therefore not able to simulate the SMB of an ice sheet. Evaluating the added value of an RCM compared to a GCM could only be done by comparing atmospheric outputs (temperature, wind, precipitation &amp;#8230;) in both models. CESM is the first ESM including a land model capable of simulating the surface of an ice sheet and thus to directly compare the SMB of an RCM and an ESM the first time.&lt;/p&gt;&lt;p&gt;Our results show that, if the SMB and is components are very similar in CESM and MAR over the present era, they quickly start to diverge in our future projection, the SMB of MAR decreasing more than that of CESM. This difference in SMB evolution is almost exclusively explained by a much larger increase of the melter runoff in MAR compared to CESM whereas the temporal evolution of snowfall, rainfall and sublimation is comparable in both runs.&lt;/p&gt;


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 11 (10) ◽  
pp. 4155-4174 ◽  
Author(s):  
Benjamin Brown-Steiner ◽  
Noelle E. Selin ◽  
Ronald Prinn ◽  
Simone Tilmes ◽  
Louisa Emmons ◽  
...  

Abstract. While state-of-the-art complex chemical mechanisms expand our understanding of atmospheric chemistry, their sheer size and computational requirements often limit simulations to short lengths or ensembles to only a few members. Here we present and compare three 25-year present-day offline simulations with chemical mechanisms of different levels of complexity using the Community Earth System Model (CESM) Version 1.2 CAM-chem (CAM4): the Model for Ozone and Related Chemical Tracers, version 4 (MOZART-4) mechanism, the Reduced Hydrocarbon mechanism, and the Super-Fast mechanism. We show that, for most regions and time periods, differences in simulated ozone chemistry between these three mechanisms are smaller than the model–observation differences themselves. The MOZART-4 mechanism and the Reduced Hydrocarbon are in close agreement in their representation of ozone throughout the troposphere during all time periods (annual, seasonal, and diurnal). While the Super-Fast mechanism tends to have higher simulated ozone variability and differs from the MOZART-4 mechanism over regions of high biogenic emissions, it is surprisingly capable of simulating ozone adequately given its simplicity. We explore the trade-offs between chemical mechanism complexity and computational cost by identifying regions where the simpler mechanisms are comparable to the MOZART-4 mechanism and regions where they are not. The Super-Fast mechanism is 3 times as fast as the MOZART-4 mechanism, which allows for longer simulations or ensembles with more members that may not be feasible with the MOZART-4 mechanism given limited computational resources.


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.


2019 ◽  
Vol 32 (9) ◽  
pp. 2627-2649 ◽  
Author(s):  
Yassir A. Eddebbar ◽  
Keith B. Rodgers ◽  
Matthew C. Long ◽  
Aneesh C. Subramanian ◽  
Shang-Ping Xie ◽  
...  

AbstractThe oceanic response to recent tropical eruptions is examined in Large Ensemble (LE) experiments from two fully coupled global climate models, the Community Earth System Model (CESM) and the Geophysical Fluid Dynamics Laboratory Earth System Model (ESM2M), each forced by a distinct volcanic forcing dataset. Following the simulated eruptions of Agung, El Chichón, and Pinatubo, the ocean loses heat and gains oxygen and carbon, in general agreement with available observations. In both models, substantial global surface cooling is accompanied by El Niño–like equatorial Pacific surface warming a year after the volcanic forcing peaks. A mechanistic analysis of the CESM and ESM2M responses to Pinatubo identifies remote wind forcing from the western Pacific as a major driver of this El Niño–like response. Following eruption, faster cooling over the Maritime Continent than adjacent oceans suppresses convection and leads to persistent westerly wind anomalies over the western tropical Pacific. These wind anomalies excite equatorial downwelling Kelvin waves and the upwelling of warm subsurface anomalies in the eastern Pacific, promoting the development of El Niño conditions through Bjerknes feedbacks a year after eruption. This El Niño–like response drives further ocean heat loss through enhanced equatorial cloud albedo, and dominates global carbon uptake as upwelling of carbon-rich waters is suppressed in the tropical Pacific. Oxygen uptake occurs primarily at high latitudes, where surface cooling intensifies the ventilation of subtropical thermocline waters. These volcanically forced ocean responses are large enough to contribute to the observed decadal variability in oceanic heat, carbon, and oxygen.


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.


2019 ◽  
Author(s):  
Takasumi Kurahashi-Nakamura ◽  
André Paul ◽  
Guy Munhoven ◽  
Ute Merkel ◽  
Michael Schulz

Abstract. We developed a coupling scheme for the Community Earth System Model version 1.2 (CESM1.2) and the Model of Early Diagenesis in the Upper Sediment of Adjustable complexity (MEDUSA), and explored the effects of the coupling on solid components in the upper sediment and on bottom seawater chemistry by comparing the coupled model's behaviour with that of the uncoupled CESM having a simplified treatment of sediment processes. CESM is a fully-coupled atmosphere-ocean-sea ice-land model and its ocean component (the Parallel Ocean Program version 2, POP2) includes a biogeochemical component (BEC). MEDUSA was coupled to POP2 in an off-line manner so that each of the models ran separately and sequentially with regular exchanges of necessary boundary condition fields. This development was done with the ambitious aim of a future application for long-term (spanning a full glacial cycle; i.e., ~ 105 years) climate simulations with a state-of-the-art comprehensive climate model including the carbon cycle, and was motivated by the fact that until now such simulations have been done only with less-complex climate models. We found that the sediment-model coupling already had non-negligible immediate advantages for ocean biogeochemistry in millennial-time-scale simulations. First, the MEDUSA-coupled CESM outperformed the uncoupled CESM in reproducing an observation-based global distribution of sediment properties, especially for organic carbon and opal. Thus, the coupled model is expected to act as a better bridge between climate dynamics and sedimentary data, which will provide another measure of model performance. Second, in our experiments, the MEDUSA-coupled model and the uncoupled model had a difference of 0.2‰ or larger in terms of δ13C of bottom water over large areas, which implied potential significant model biases for bottom seawater chemical composition due to a different way of sediment treatment. Such a model bias would be a fundamental issue for paleo model–data comparison often relying on data derived from benthic foraminifera.


Sign in / Sign up

Export Citation Format

Share Document