This study evaluates the performance of a global storm-resolving model (GSRM), the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM). We analyze marine boundary layer clouds in a cold air outbreak over the Norwegian Sea in a 40-day simulation, and compare them to observations from satellite and a field campaign of the Atmospheric Radiation Measurement program (ARM). SCREAM qualitatively captures the cold air outbreak cloud transition in terms of the boundary layer growth, cloud mesoscale structure, and phase partitioning. SCREAM also correctly locates the greatest ice and liquid in the mesoscale updraft. However, the study finds that SCREAM might underestimate cloud supercooled liquid water in the cumulus cloud regime. This study showcases the promise of employing high-resolution and high-frequency observations under similar large-scale conditions for evaluating GSRMs. This approach can help identify model features for future process-level studies before allocating extra resources for a time-matched model intercomparison of a specific case.

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.

Can the current successes of global machine learning-based weather simulators be generalized beyond two-week forecasts to stable and accurate multiyear runs? The recently developed AI2 Climate Emulator (ACE) suggests this is feasible, based upon 10-year simulations trained on a realistic global atmosphere model using a grid spacing of approximately 110~km and forced by a repeating annual cycle of sea-surface temperature. Here we show that ACE, without modification, can be trained to emulate another major atmospheric model, EAMv2, run at a comparable grid spacing for at least ten years with similarly small climate biases. ACE accurately reproduces EAMv2’s frequency distribution of daily-mean precipitation, its time-mean spatial pattern of precipitation, and its space-time structure of tropical precipitation, including the Madden-Julian Oscillation. Moreover, ACE’s climate biases with respect to EAMv2 are substantially smaller than EAMv2’s own biases compared to the observed historical average surface precipitation rate and top-of-atmosphere radiative fluxes.

We develop a doubly periodic version of the Simple Convection-Permitting E3SM Atmosphere Model (SCREAM) to provide an “efficient” configuration for this global storm resolving model (GSRM), akin to a single column model (SCM) often found in conventional general circulation models (GCMs). The design details are explained, in addition to the extensive case library associated with the doubly periodic SCREAM (DP-SCREAM) configuration. We demonstrate that doubly periodic cloud resolving models are useful tools to explore the scale awareness and scale sensitivity of GSRMs, in addition to replicating biases seen in the global models. Using DP-SCREAM, we show that SCREAM is a scale aware model as it is able to realistically partition between sub-grid scale (SGS) and resolved vertical transport across the gray zone of turbulence. We show that SCREAM is reasonably scale insensitive when run at resolutions from 1 to 5 km, but can exhibit sensitivity, particularly for the shallow convective regime, when run at resolutions approaching that of large eddy simulations. We conclude that SGS parameterization improvements are likely needed to reduce this scale sensitivity.

The water cycle is an important component of the earth system and it plays a key role in many facets of society, including energy production, agriculture, and human health and safety. In this study, the Energy Exascale Earth System Model version 1 (E3SMv1) is run with low-resolution (roughly 110 km) and high-resolution (roughly 25 km) configurations — as established by the High Resolution Model Intercomparison Project protocol — to evaluate the atmospheric and terrestrial water budgets over the conterminous United States (CONUS) at the large watershed scale. The water cycle slows down in the HR experiment relative to the LR, with decreasing fluxes of precipitation, evapotranspiration, atmospheric moisture convergence, and runoff. The reductions in these terms exacerbate biases for some watersheds, while reducing them in others. For example, precipitation biases are exacerbated at HR over the Eastern and Central CONUS watersheds, while precipitation biases are reduced at HR over the Western CONUS watersheds. The most pronounced changes to the water cycle come from reductions in precipitation and evapotranspiration, the latter of which results from decreases in evaporative fraction. While the HR simulation is warmer than the LR, moisture convergence decreases despite the increased atmospheric water vapor, suggesting circulation biases are an important factor. Additional exploratory metrics show improvements to water cycle extremes (both in precipitation and streamflow), fractional contributions of different storm types to total precipitation, and mountain snowpack.

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.

In this study, we find significant sensitivity to the choice of time step for the Energy Exascale Earth System Model's atmospheric component, leading to large decreases in the magnitude of cloud forcing when the time step is reduced to 10 seconds. Reducing the time step size for the microphysics increases precipitation, leading to a drying of the atmosphere and an increase in surface evaporation. This effect is amplified when the microphysics is substepped together with other cloud physics processes. Coupling the model's dynamics and physics more frequently reduces cloud fraction at lower altitudes, while producing more cloud liquid at higher altitudes. Reducing the deep convection time step also reduces low cloud mass and cloud fraction. Together, these results suggest that cloud physics in a global circulation model can depend strongly on time step, and in particular on the frequency with which cloud-related processes are coupled with each other and with the model dynamics.

In this paper, we introduce a testbed for evaluating and comparing climate modeling systems at cloud resolving scales using hindcasts of the June 2012 North American derecho. The testbed is applied to two models: the regionally-refined Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) at horizontal resolutions ranging from 6.5 to 1.625 km and the Weather Research and Forecasting (WRF) model with 4 km grid spacing. We find the simulation results to be highly sensitive to the initial conditions, initialization time, and model configurations, with initial conditions from the Rapid Refresh (RAP) producing the best simulation. Significant improvement is identified in the SCREAM simulations as horizontal grid spacing is refined. While a propagation delay of approximately 2 hours is found in both models, SCREAM at 1.625 km simulates the observed bow echo structure of the derecho well and predicts strong surface gusts that exceed 30 m/s. In comparison, WRF hardly produces surface wind over 25 m/s, and the derecho wind gust in WRF is 42-46% lower than in SCREAM. Moreover, WRF has a lower bias in simulating cold clouds but overestimates the precipitation intensity. Both models well reproduce the observed outgoing longwave radiation spatial patterns (Pearson correlation > 0.88) while they simulate larger areas of composite radar reflectivity > 40 dBZ by up to 4 times and underestimate the precipitating area by ~ 70\% in WRF and 47\% in SCREAM compared to observations.

Climate models rely on parameterizations of a variety of processes in the atmospheric physics, but a common concern is that the temporal resolution is too coarse to consistently resolve the behavior that individual parameterizations are designed to capture. This study examines timescales numerically derived from the Morrison-Gettelman (MG2) microphysics as implemented within the Energy Exascale Earth System Model, version 1 (E3SMv1). Numerically-relevant timescales in MG2 are derived by computing the eigenspectrum of its Jacobian. These timescales are found to often be smaller than the default 5 min timestep used for MG2. The fast timescales are then heuristically connected to individual microphysics processes. By substepping a few particular rain processes within MG2, the time discretization error for those processes was significantly reduced with minimal additional expense to the overall microphysics. While this improvement has a substantial effect on the target processes and on the vertical distribution of stratiform-derived rain within E3SMv1, the overall model climate is found to not be sensitive to MG2 time step.

A common goal of next-generation Global Circulation Models (GCMs) is that they should be “scale-aware”, which typically implies that such models should not be excessively sensitive to grid spacings, and that they should in some sense converge monotonically towards a result as grid spacings decrease. While both horizontal and vertical resolution have been treated in this manner, time resolution is typically viewed differently. Specifically, a decrease in time step size is often viewed as a “necessary evil”, being decreased only in cases where spatial resolution is also decreased, requiring a change to the time resolution to satisfy a CFL condition. Our experiments with the E3SM Atmosphere Model suggest that cloud physics and precipitation in GCMs is in fact quite sensitive to process coupling time step size, and that the biases affected by time integration error are independent from (and of comparable size to) biases due to other common sources of error, such as grid spacing and choice of sub-grid-scale physics parameterizations. This suggests that process coupling frequency is a key feature that should be adjusted for future models.

The spread in global mean precipitation among climate models is explored in two ensembles using the complementary perspectives of surface evaporation and energy budgets. Models with higher global-mean precipitation have stronger oceanic evaporation, driven by drier near-surface air. The drier surface conditions occur alongside increases in near-surface temperature and moisture at 925 hPa, which point to stronger boundary layer mixing. Correlations suggest that the degree of lower tropospheric mixing explains 18%-49% of the intermodel precipitation variance. To test this hypothesis, the degree of mixing is varied in a single-model experiment by adjusting the relative humidity threshold that controls low-cloud fraction. Indeed, increasing lower tropospheric mixing results in more global mean precipitation. Energetically, increased precipitation rates are associated with more downwelling longwave radiation to the surface and weaker sensible heat fluxes. These results highlight how lower-tropospheric processes must be better constrained to reduce the precipitation discrepancy among climate models.