Satellite observations and their corresponding instrument simulators are used to document global cloud biases in the Community Atmosphere Model (CAM) versions 4 and 5. The model-observation comparisons show that, despite having nearly identical cloud radiative forcing, CAM5 has a much more realistic representation of cloud properties than CAM4. In particular, CAM5 exhibits substantial improvement in three long-standing climate model cloud biases: 1) the underestimation of total cloud, 2) the overestimation of optically thick cloud, and 3) the underestimation of midlevel cloud. While the increased total cloud and decreased optically thick cloud in CAM5 result from improved physical process representation, the increased midlevel cloud in CAM5 results from the addition of radiatively active snow. Despite these improvements, both CAM versions have cloud deficiencies. Of particular concern, both models exhibit large but differing biases in the subtropical marine boundary layer cloud regimes that are known to explain intermodel differences in cloud feedbacks and climate sensitivity. More generally, this study demonstrates that simulatorfacilitated evaluation of cloud properties, such as amount by vertical level and optical depth, can robustly expose large and at times radiatively compensating climate model cloud biases.

This paper describes the first implementation of the Δ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 (January 20 through February 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 the horizontal dynamics 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., tropical cyclones, extratropical cyclones 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.

Results from the new Department of Energy super-parameterized (SP) Energy ExascaleEarth System Model (SP-E3SM) are analyzed and compared to the traditionally parameterized E3SMv1 and previous studies using SP models. SP-E3SM is unique in that it utilizes Graphics Processing Unit hardware acceleration, cloud resolving model mean-state acceleration, and reduced radiation to dramatically increase the model throughput and allow decadal experiments at 100-km external resolution. It also differs from other SP models by using a spectral element dynamical core on a cubed-sphere grid and a finer vertical grid with a higher model top. Despite these differences, SP-E3SM generally reproduces the behavior of other SP models. Tropical wave variability is improved relative to E3SM, including the emergence of a Madden-Julian Oscillation and a realistic slowdown of Moist Kelvin Waves. However, the distribution of precipitation exhibits indicates an overly frequent occurrence of rain rates less than 1 mm day −1 , and while the timing of diurnal rainfall shows modest improvements the signal is not as coherent as observations. A notable grid imprinting bias is identified in the precipitation field and attributed to a unique feedback associated with the interactions between the explicit cloud resolving model convection and the spectral element grid structure. Spurious zonal mean column water tendencies due to grid imprinting are quantified-while negligible for the conventionally parameterized E3SM, they become large with super-parameterization, approaching 10% of the physical tendencies. The implication is that finding a remedy to grid imprinting will become especially important as spectral element dynamical cores begin to be combined with explicitly resolved convection.

This work documents version two of the Department of Energy's Energy Exascale Earth SystemModel (E3SM). 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 with Coupled Model Intercomparison Project Phase 6 Diagnosis, Evaluation, and Characterization of Klima simulations augmented with historical simulations as well as simulations to evaluate impacts of different forcing agents. The simulated climate has many realistic features of the climate system, 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 Program assessment. However, a number of important biases remain including a weak Atlantic Meridional Overturning Circulation, deficiencies in the characteristics and spectral distribution of tropical atmospheric variability, and a significant underestimation of the observed warming in the second half of the historical period. 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.

Satellite simulators are often used to account for limitations in satellite retrievals of cloud properties in comparisons between models and satellite observations. The purpose of this framework is to enable more robust evaluation of model cloud properties, so that differences between models and observations can more confidently be attributed to model errors. A critical step in this process is accounting for the difference between the spatial scales at which cloud properties are retrieved with those at which clouds are simulated in global models. In this study, we create a series of sensitivity tests using 4‐km global model output from the multiscale modeling framework to evaluate the sensitivity of simulated satellite retrievals to common assumptions about cloud and precipitation overlap and condensate variability used in climate models whose grid spacing is many tens to hundreds of kilometers. We find the simulated retrievals are sensitive to these assumptions. Using maximum‐random overlap with homogeneous cloud and precipitation condensate leads to errors in Multiangle Imaging Spectroradiometer and International Satellite Cloud Climatology Project‐simulated cloud cover and in CloudSat‐simulated radar reflectivity that are significant compared to typical differences between the model simulations and observations. A more realistic treatment of unresolved clouds and precipitation is shown to substantially reduce these errors. The sensitivity to these assumptions underscores the need for the adoption of more realistic subcolumn treatments in models and the need for consistency among subcolumn assumptions between models and simulators to ensure that simulator‐diagnosed errors are consistent with the model formulation.

No abstract

Classical numerical models for the global atmosphere, as used for numerical weather forecasting or climate research, have been developed for conventional central processing unit (CPU) architectures. This hinders the employment of such models on current top-performing supercomputers, which achieve their computing power with hybrid architectures, mostly using graphics processing units (GPUs). Thus also scientific applications of such models are restricted to the lesser computer power of CPUs. Here we present the development of a GPU-enabled version of the ICON atmosphere model (ICON-A), motivated by a research project on the quasi-biennial oscillation (QBO), a globalscale wind oscillation in the equatorial stratosphere that depends on a broad spectrum of atmospheric waves, which originates from tropical deep convection. Resolving the relevant scales, from a few kilometers to the size of the globe, is a formidable computational problem, which can only be realized now on top-performing supercomputers. This motivated porting ICON-A, in the specific configuration needed for the research project, in a first step to the GPU architecture of the Piz Daint computer at the Swiss National Supercomputing Centre and in a second step to the JUWELS Booster computer at the Forschungszentrum Jülich. On Piz Daint, the ported code achieves a single-node GPU vs. CPU speedup factor of 6.4 and allows for global experiments at a horizontal resolution of 5 km on 1024 computing nodes with 1 GPU per node with a turnover of 48 simulated days per day. On JUWELS Booster, the more modern hardware in combination with an upgraded code base allows for simulations at the same resolution on 128 computing nodes with 4 GPUs per node and a turnover of 133 simulated days per day. Additionally, the code still remains functional on CPUs, as is demonstrated by additional experiments on the Levante compute system at the German Climate Computing Center. While the application shows good weak scaling over the tested 16fold increase in grid size and node count, making also higher resolved global simulations possible, the strong scaling on GPUs is relatively poor, which limits the options to increase Published by Copernicus Publications on behalf of the European Geosciences Union. 6986 M. A. Giorgetta et al.: The ICON-A model on GPUs turnover with more nodes. Initial experiments demonstrate that the ICON-A model can simulate downward-propagating QBO jets, which are driven by wave-mean flow interaction.

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