This work describes the implementation and evaluation of the Slab Ocean Model component of the Energy Exascale Earth System Model version 2 (E3SMv2-SOM), and its application to understanding the climate sensitivity to ocean heat transport (OHT) strength and CO$_{2}$ forcing. E3SMv2-SOM reproduces the baseline climate and Equilibrium Climate Sensitivity (ECS) of the E3SMv2 fully coupled experiments, reasonably well, with a pattern correlation close to 1 and global mean bias that is less than 1$\%$ of the fully coupled surface temperature, precipitation and sea ice extent and volume. Similar to other model behaviour, the ECS estimated from the SOM (4.5$^\circ$C) is greater than the estimate from fully coupled model (4.0$^\circ$C; from 150 years regression). The E3SMv2 baseline climate is also very sensitive to the strength of the OHT from which the prescribed ocean heat convergence (OHC) for the SOM is derived, with a surface temperature difference of about 4.0$^\circ$C between high- and low-OHT SOM experiments. The surface temperature response in the high/low-OHT experiments occur through a positive/negative Shortwave cloud radiative effect, caused by a decrease/increase in marine low-level clouds over subpolar regions. This surface temperature sensitivity to prescribed OHCs is particularly large in the Southern hemisphere and is associated with an overcompensation of between prescribed OHC/OHT by atmosphere heat transports. This large sensitivity indicates stronger low-level cloud feedbacks in E3SM. The SOMâ\euro™s ECS estimate is also sensitive to the baseline climate it is initialized from, with an ECS difference of 0.5$^\circ$C between the high- and low- OHT CO$_2$ increase experiments.

The new generation of heterogeneous CPU/GPU computer systems offer much greater computational performance but are not yet widely used for climate modeling. One reason for this is that traditional climate models were written before GPUs were available and would require an extensive overhaul to run on these new machines. In addition, even conventional “high–resolution’ simulations don’t provide enough parallel work to keep GPUs busy, so the benefits of such overhaul would be limited for the types of simulations climate scientists are accustomed to. The vision of the Simple Cloud-Resolving Energy Exascale Earth System (E3SM) Atmosphere Model (SCREAM) project is to create a global atmospheric model with the architecture to efficiently use GPUs and horizontal resolution sufficient to fully take advantage of GPU parallelism. After 5 yrs of model development, SCREAM is finally ready for use. In this paper, we describe the design of this new code, its performance on both CPU and heterogeneous machines, and its ability to simulate real-world climate via a set of four 40 day simulations covering all 4 seasons of the year.

This paper describes the first implementation of the d x=3.25 km version of the Energy Exascale Earth System Model (E3SM) global atmosphere model and its behavior in a 40 day prescribed-sea-surface-temperature simulation (Jan 20-Feb 28, 2020). This simulation was performed as part of the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) phase 2 model intercomparison. Effective resolution is found to be $\sim 6x the horizontal grid resolution despite using a coarser grid for physical parameterizations. Despite this new model being in an immature and untuned state, moving to 3.25 km grid spacing solves several long-standing problems with the E3SM model. In particular, Amazon precipitation is much more realistic, the frequency of light and heavy precipitation is improved, agreement between the simulated and observed diurnal cycle of tropical precipitation is excellent, and the vertical structure of tropical convection and coastal stratocumulus look good. In addition, the new model is able to capture the frequency and structure of important weather events (e.g. hurricanes, midlatitude storms including atmospheric rivers, and cold air outbreaks). Interestingly, this model does not get rid of the erroneous southern branch of the intertropical convergence zone nor the tendency for strongest convection to occur over the Maritime Continent rather than the West Pacific, both of which are classic climate model biases. Several other problems with the simulation are identified, underscoring the fact that this model is a work in progress.

This work documents version two of the Department of Energy’s Energy Exascale Earth System Model (E3SM). E3SM version 2 (E3SMv2) is a significant evolution from its predecessor E3SMv1, resulting in a model that is nearly twice as fast and with a simulated climate that is improved in many metrics. We describe the physical climate model in its lower horizontal resolution configuration consisting of 110 km atmosphere, 165 km land, 0.5° river routing model, and an ocean and sea ice with mesh spacing varying between 60 km in the mid-latitudes and 30 km at the equator and poles. The model performance is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima (DECK) simulations augmented with historical simulations as well as simulations to evaluate impacts of different forcing agents. The simulated climate is generally realistic, with notable improvements in clouds and precipitation compared to E3SMv1. E3SMv1 suffered from an excessively high equilibrium climate sensitivity (ECS) of 5.3 K. In E3SMv2, ECS is reduced to 4.0 K which is now within the plausible range based on a recent World Climate Research Programme (WCRP) assessment. However, E3SMv2 significantly underestimates the global mean surface temperature in the second half of the historical record. An analysis of single-forcing simulations indicates that correcting the historical temperature bias would require a substantial reduction in the magnitude of the aerosol-related forcing.

Snow and ice albedo reduction due to deposition of absorbing particles (i.e., snow darkening effect (SDE)) warms the Earth system and is largely attributed to black carbon (BC) and dust. Absorbing organic aerosol (BrC) also contributes to SDE but has received less attention due to uncertainty and challenges in model representation. This work incorporates the SDE of absorbing organic aerosol (BrC) from biomass burning and biofuel sources into the Snow Ice and Aerosol Radiative (SNICAR) model within a variant of the Community Earth System Model (CESM). Additionally, 12 different emission regions of BrC and BC from biomass burning and biofuel sources are tagged to quantify the relative contribution to global and regional SDE. BrC global SDE (0.021–0.056 Wm-2) is larger than other model estimates, corresponding to 37%–98% of the SDE from BC. When compared to observations, BrC simulations have a range in median bias (-2.5%–+21%), with better agreement in the simulations that include BrC photochemical bleaching. The largest relative contributions to global BrC SDE are traced to Northern Asia (23%–31%), Southeast Asia (16%–21%), and South Africa (13%–17%). Transport from Southeast Asia contributes nearly half of the regional BrC SDE in Antarctica (0.084–0.3 Wm-2), which is the largest regional input to global BrC SDE. Lower latitude BrC SDE is correlated with snowmelt, in-snow BrC concentrations, and snow cover fraction, while polar BrC SDE is correlated with surface insolation and snowmelt. This indicates the importance of in-snow processes and snow feedbacks on modeled BrC SDE.

Nitrate aerosol plays an important role in affecting regional air quality as well as Earth’s climate. However, it is not well represented or even neglected in many global climate models. In this study, we couple the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) module with the four-mode version of the Modal Aerosol Module (MAM4) in DOE’s Energy Exascale Earth System Model version 2 (E3SMv2) to treat nitrate aerosol and its radiative effects. We find that nitrate aerosol simulated by E3SMv2-MAM4-MOSAIC is sensitive to the treatment of gaseous HNO3 transfer to/from interstitial particles related to accommodation coefficients of HNO3 (αHNO3) on dust and non-dust particles. We compare three different treatments of HNO3 transfer: 1) a treatment (MTC_SLOW) that uses a low αHNO3 in the mass transfer coefficient (MTC) calculation; 2) a dust-weighted MTC treatment (MTC_WGT) that uses a high αHNO3 on non-dust particles; and 3) a dust-weighted MTC treatment that also splits coarse mode aerosols into the coarse dust and sea salt sub-modes in MOSAIC (MTC_SPLC). MTC_WGT and MTC_SPLC increase the global annual mean (2005-2014) nitrate burden from 0.096 (MTC_SLOW) to 0.237 and 0.185 Tg N, respectively, mostly in the coarse mode. They also produce stronger nitrate direct radiative forcing (–0.048 and –0.051 W m–2, respectively) and indirect forcing (–0.33 and –0.35 W m–2, respectively) than MTC_SLOW (–0.021 and –0.24 W m–2). All three treatments overestimate nitrate surface concentrations compared with ground-based observations. MTC_WGT and MTC_SPLC improve the vertical profiles of nitrate concentrations against aircraft measurements below 400 hPa.