The Whole Atmosphere Community Climate Model version 6 (WACCM6) is a major update of the whole atmosphere modeling capability in the Community Earth System Model (CESM), featuring enhanced physical, chemical and aerosol parameterizations. This work describes WACCM6 and some of the important features of the model. WACCM6 can reproduce many modes of variability and trends in the middle atmosphere, including the quasi‐biennial oscillation, stratospheric sudden warmings, and the evolution of Southern Hemisphere springtime ozone depletion over the twentieth century. WACCM6 can also reproduce the climate and temperature trends of the 20th century throughout the atmospheric column. The representation of the climate has improved in WACCM6, relative to WACCM4. In addition, there are improvements in high‐latitude climate variability at the surface and sea ice extent in WACCM6 over the lower top version of the model (CAM6) that comes from the extended vertical domain and expanded aerosol chemistry in WACCM6, highlighting the importance of the stratosphere and tropospheric chemistry for high‐latitude climate variability.

The NCAR Community Climate System Model, version 3 (CCSM3) exhibits persistent errors in its simulation of the El Niño–Southern Oscillation (ENSO) mode of coupled variability. The amplitude of the oscillation is too strong, the dominant 2-yr period too regular, and the width of the sea surface temperature response in the Pacific too narrow, with positive anomalies extending too far into the western Pacific. Two changes in the parameterization of deep convection result in a significant improvement to many aspects of the ENSO simulation. The inclusion of convective momentum transport (CMT) and a dilution approximation for the calculation of convective available potential energy (CAPE) are used in development integrations, and a striking improvement in ENSO characteristics is seen. An increase in the periodicity of ENSO is achieved by a reduction in the strength of the existing “short-circuited” delayed-oscillator mode. The off-equatorial response is weaker and less tropically confined, largely as a result of the CMT and an associated redistribution of zonal momentum. The Pacific east–west structure is improved in response to the presence of convective dilution and cooling provided by increased surface fluxes. The initiation of El Niño events is fundamentally different. Enhanced intraseasonal surface stress variability leads to absolute surface westerlies and a cooling–warming dipole between the Philippine Sea and western Pacific. Lag-regression analysis shows that intraseasonal variability may play a significant role in event initiation and maintenance as opposed to being a benign response to increased SSTs. Recent observational evidence appears to support such a leading relationship.

Middle atmospheric general circulation models (GCMs) must employ a parameterization for small-scale gravity waves (GWs). Such parameterizations typically make very simple assumptions about gravity wave sources, such as uniform distribution in space and time or an arbitrarily specified GW source function. The authors present a configuration of the Whole Atmosphere Community Climate Model (WACCM) that replaces the arbitrarily specified GW source spectrum with GW source parameterizations. For the nonorographic wave sources, a frontal system and convective GW source parameterization are used. These parameterizations link GW generation to tropospheric quantities calculated by the GCM and provide a model-consistent GW representation. With the new GW source parameterization, a reasonable middle atmospheric circulation can be obtained and the middle atmospheric circulation is better in several respects than that generated by a typical GW source specification. In particular, the interannual NH stratospheric variability is significantly improved as a result of the source-oriented GW parameterization. It is also shown that the addition of a parameterization to estimate mountain stress due to unresolved orography has a large effect on the frequency of stratospheric sudden warmings in the NH stratosphere by changing the propagation of stationary planetary waves into the polar vortex.

This work documents version two of the Department of Energy's Energy Exascale Earth System Model (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.

The Energy Exascale Earth System Model Atmosphere Model version 1, the atmospheric component of the Department of Energy's Energy Exascale Earth System Model is described. The model began as a fork of the well‐known Community Atmosphere Model, but it has evolved in new ways, and coding, performance, resolution, physical processes (primarily cloud and aerosols formulations), testing and development procedures now differ significantly. Vertical resolution was increased (from 30 to 72 layers), and the model top extended to 60 km (~0.1 hPa). A simple ozone photochemistry predicts stratospheric ozone, and the model now supports increased and more realistic variability in the upper troposphere and stratosphere. An optional improved treatment of light‐absorbing particle deposition to snowpack and ice is available, and stronger connections with Earth system biogeochemistry can be used for some science problems. Satellite and ground‐based cloud and aerosol simulators were implemented to facilitate evaluation of clouds, aerosols, and aerosol‐cloud interactions. Higher horizontal and vertical resolution, increased complexity, and more predicted and transported variables have increased the model computational cost and changed the simulations considerably. These changes required development of alternate strategies for tuning and evaluation as it was not feasible to “brute force” tune the high‐resolution configurations, so short‐term hindcasts, perturbed parameter ensemble simulations, and regionally refined simulations provided guidance on tuning and parameterization sensitivity to higher resolution. A brief overview of the model and model climate is provided. Model fidelity has generally improved compared to its predecessors and the CMIP5 generation of climate models.

Transport of momentum by convection is an important process affecting global circulation. Owing to the lack of global observations, the quantification of the impact of this process on the tropospheric climate is difficult. Here an implementation of two convective momentum transport parameterizations, presented by Schneider and Lindzen and Gregory et al., in the Community Atmosphere Model, version 3 (CAM3) is presented, and their effect on global climate is examined in detail. An analysis of the tropospheric zonal momentum budget reveals that convective momentum transport affects tropospheric climate mainly through changes to the Coriolis torque. These changes result in improvement of the representation of the Hadley circulation: in December–February, the upward branch of the circulation is weakened in the Northern Hemisphere and strengthened in the Southern Hemisphere, and the lower northerly branch is weakened. In June–August, similar improvements are noted. The inclusion of convective momentum transport in CAM3 reduces many of the model’s biases in the representation of surface winds, as well as in the representation of tropical convection. In an annual mean, the tropical easterly bias, subtropical westerly bias, and the bias in the 60°S jet are improved. Representation of convection is improved along the equatorial belt with decreased precipitation in the Indian Ocean and increased precipitation in the western Pacific. The improvements of the representation of tropospheric climate are greater with the implementation of the Schneider and Lindzen parameterization.