Using the global 3.25-km Simple Cloud Resolving E3SM Atmosphere Model (SCREAM 3 km), a pair of 13-month Cess-Potter simulations are performed to quantify the radiative feedbacks and the hydrologic and circulation responses to warming. Large-scale aspects of SCREAM 3 km’s top-of-atmosphere radiative fluxes, precipitation rates, and circulations are in good agreement with observations and reanalysis, with notable differences, including a drier lower free-troposphere in the Tropics, reduced precipitation and humidity over the Tropical West Pacific, and poleward shifted Southern Hemisphere midlatitude jet. In response to warming, SCREAM 3 km predicts a total radiative feedback within the top 15% of the CMIP5 and CMIP6 models, which puts it substantially higher than the feedback reported by other kilometer-scale models. SCREAM 3 km’s high radiative feedback stems from a strongly positive shortwave cloud feedback, most prominent over the mid- and high-latitudes. SCREAM 3 km’s high precipitation response also puts it among the highest of CMIP models, whereas its circulation response are unremarkable compared to CMIP models. An ensemble of five perturbed initial condition Cess-Potter simulations with a 12 km version of SCREAM (SCREAM 12 km) is performed to characterize uncertainty and resolution sensitivity. It suggests that the uncertainty from analyzing a pair of one-year simulations is small compared to the inter-model spread in feedbacks and precipitation response. SCREAM 12 km also produces a strong precipitation response to warming but a much lower cloud feedback and total radiative feedback. The results from these experiments suggest that the spread in climate feedbacks will likely persist in the next generation of kilometer-scale models.

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.

One of the largest sources of uncertainty in climate models comes from convective cloudssystems parametrizations, which are necessary at current horizontal resolutions, of about 25km.In order to fully resolve the effects of such systems, horizontal resolutions need to reach the 1-3km range, which entail a huge increase in computational costs. In recent years, the performance of the world largest supercomputers has steadilyimproved, to the point that we can now think of running a climate modelat cloud resolving resolution, with a Simulated Years Per Day (SYPD) throughputthat is suitable for decade or century long simulations.However, the architectures of the newest supercomputers are becoming more andmore heterogeneous, which makes the task of keeping a performant code base moreand more challenging. Here, we present our effort to upgrade the new non-hydrostatic atmosphere dycore ofE3SM to a version that is highly performant on a variety of architectures, includingGPUs, conventional CPUs, and many-core CPUs. When using GPUs, our implementation isone of the fastest, achieving 0.97 SYPD on the NGPPS benchmark, when running on the full Summit supercomputer.