While internal climate variability is known to affect climate projections, its influence is often underappreciated and confused with model error. Why? In general, modeling centers contribute a small number of realizations to international climate model assessments [e.g., phase 5 of the Coupled Model Intercomparison Project (CMIP5)]. As a result, model error and internal climate variability are difficult, and at times impossible, to disentangle. In response, the Community Earth System Model (CESM) community designed the CESM Large Ensemble (CESM-LE) with the explicit goal of enabling assessment of climate change in the presence of internal climate variability. All CESM-LE simulations use a single CMIP5 model (CESM with the Community Atmosphere Model, version 5). The core simulations replay the twenty to twenty-first century (1920–2100) 30 times under historical and representative concentration pathway 8.5 external forcing with small initial condition differences. Two companion 1000+-yr-long preindustrial control simulations (fully coupled, prognostic atmosphere and land only) allow assessment of internal climate variability in the absence of climate change. Comprehensive outputs, including many daily fields, are available as single-variable time series on the Earth System Grid for anyone to use. Early results demonstrate the substantial influence of internal climate variability on twentieth- to twenty-first-century climate trajectories. Global warming hiatus decades occur, similar to those recently observed. Internal climate variability alone can produce projection spread comparable to that in CMIP5. Scientists and stakeholders can use CESM-LE outputs to help interpret the observational record, to understand projection spread and to plan for a range of possible futures influenced by both internal climate variability and forced climate change.

Stratospheric sudden warmings are the clearest and strongest manifestation of dynamical coupling in the stratosphere–troposphere system. While many sudden warmings have been individually documented in the literature, this study aims at constructing a comprehensive climatology: all major midwinter warming events are identified and classified, in both the NCEP–NCAR and 40-yr ECMWF Re-Analysis (ERA-40) datasets. To accomplish this a new, objective identification algorithm is developed. This algorithm identifies sudden warmings based on the zonal mean zonal wind at 60°N and 10 hPa, and classifies them into events that do and do not split the stratospheric polar vortex. Major midwinter stratospheric sudden warmings are found to occur with a frequency of approximately six events per decade, and 46% of warming events lead to a splitting of the stratospheric polar vortex. The dynamics of vortex splitting events is contrasted to that of events where the vortex is merely displaced off the pole. In the stratosphere, the two types of events are found to be dynamically distinct: vortex splitting events occur after a clear preconditioning of the polar vortex, and their influence on middle-stratospheric temperatures lasts for up to 20 days longer than vortex displacement events. In contrast, the influence of sudden warmings on the tropospheric state is found to be largely insensitive to the event type. Finally, a table of dynamical benchmarks for major stratospheric sudden warming events is compiled. These benchmarks are used in a companion study to evaluate current numerical model simulations of the stratosphere.

An overview of the Community Earth System Model Version 2 (CESM2) is provided, including a discussion of the challenges encountered during its development and how they were addressed. In addition, an evaluation of a pair of CESM2 long preindustrial control and historical ensemble simulations is presented. These simulations were performed using the nominal 1° horizontal resolution configuration of the coupled model with both the “low‐top” (40 km, with limited chemistry) and “high‐top” (130 km, with comprehensive chemistry) versions of the atmospheric component. CESM2 contains many substantial science and infrastructure improvements and new capabilities since its previous major release, CESM1, resulting in improved historical simulations in comparison to CESM1 and available observations. These include major reductions in low‐latitude precipitation and shortwave cloud forcing biases; better representation of the Madden‐Julian Oscillation; better El Niño‐Southern Oscillation‐related teleconnections; and a global land carbon accumulation trend that agrees well with observationally based estimates. Most tropospheric and surface features of the low‐ and high‐top simulations are very similar to each other, so these improvements are present in both configurations. CESM2 has an equilibrium climate sensitivity of 5.1–5.3 °C, larger than in CESM1, primarily due to a combination of relatively small changes to cloud microphysics and boundary layer parameters. In contrast, CESM2's transient climate response of 1.9–2.0 °C is comparable to that of CESM1. The model outputs from these and many other simulations are available to the research community, and they represent CESM2's contributions to the Coupled Model Intercomparison Project Phase 6.

The NCAR Community Earth System Model (CESM) now includes an atmospheric component that extends in altitude to the lower thermosphere. This atmospheric model, known as the Whole Atmosphere Community Climate Model (WACCM), includes fully interactive chemistry, allowing, for example, a self-consistent representation of the development and recovery of the stratospheric ozone hole and its effect on the troposphere. This paper focuses on analysis of an ensemble of transient simulations using CESM1(WACCM), covering the period from the preindustrial era to present day, conducted as part of phase 5 of the Coupled Model Intercomparison Project. Variability in the stratosphere, such as that associated with stratospheric sudden warmings and the development of the ozone hole, is in good agreement with observations. The signals of these phenomena propagate into the troposphere, influencing near-surface winds, precipitation rates, and the extent of sea ice. In comparison of tropospheric climate change predictions with those from a version of CESM that does not fully resolve the stratosphere, the global-mean temperature trends are indistinguishable. However, systematic differences do exist in other climate variables, particularly in the extratropics. The magnitude of the difference can be as large as the climate change response itself. This indicates that the representation of stratospheretroposphere coupling could be a major source of uncertainty in climate change projections in CESM.

The increasing richness of exoplanet observations has motivated a variety of three-dimensional atmospheric circulation models of these planets. Under strongly irradiated conditions, models of tidally locked, short-period planets (both hot Jupiters and terrestrial planets) tend to exhibit a circulation dominated by a fast eastward, or "superrotating," jet stream at the equator. When the radiative and advection time scales are comparable, this phenomenon can cause the hottest regions to be displaced eastward from the substellar point by tens of degrees longitude. Such an offset has been subsequently observed on HD 189733b, supporting the possibility of equatorial jets on short-period exoplanets. Despite its relevance, however, the dynamical mechanisms responsible for generating the equatorial superrotation in such models have not been identified. Here, we show that the equatorial jet results from the interaction of the mean flow with standing Rossby waves induced by the day-night thermal forcing. The strong longitudinal variations in radiative heating-namely intense dayside heating and nightside cooling-trigger the formation of standing, planetary-scale equatorial Rossby and Kelvin waves. The Rossby waves develop phase tilts that pump eastward momentum from high latitudes to the equator, thereby inducing equatorial superrotation. We present an analytic theory demonstrating this mechanism and explore its properties in a hierarchy of one-layer (shallow-water) calculations and fully 3D models. The wave-mean-flow interaction produces an equatorial jet whose latitudinal width is comparable to that of the Rossby waves, namely the equatorial Rossby deformation radius modified by radiative and frictional effects. For conditions typical of synchronously rotating hot Jupiters, this length is comparable to a planetary radius, explaining the broad scale of the equatorial jet obtained in most hot Jupiter models. Our theory illuminates the dependence of the equatorial jet speed on forcing amplitude, strength of friction, and other parameters, as well as the conditions under which jets can form at all.

The importance of stratospheric ozone depletion on the atmospheric circulation of the troposphere is studied with an atmospheric general circulation model, the Community Atmospheric Model, version 3 (CAM3), for the second half of the twentieth century. In particular, the relative importance of ozone depletion is contrasted with that of increased greenhouse gases and accompanying sea surface temperature changes. By specifying ozone and greenhouse gas forcings independently, and performing long, time-slice integrations, it is shown that the impacts of ozone depletion are roughly 2–3 times larger than those associated with increased greenhouse gases, for the Southern Hemisphere tropospheric summer circulation. The formation of the ozone hole is shown to affect not only the polar tropopause and the latitudinal position of the midlatitude jet; it extends to the entire hemisphere, resulting in a broadening of the Hadley cell and a poleward extension of the subtropical dry zones. The CAM3 results are compared to and found to be in excellent agreement with those of the multimodel means of the recent Coupled Model Intercomparison Project (CMIP3) and Chemistry–Climate Model Validation (CCMVal2) simulations. This study, therefore, strongly suggests that most Southern Hemisphere tropospheric circulation changes, in austral summer over the second half of the twentieth century, have been caused by polar stratospheric ozone depletion.

Horizontal temperature gradients are small in the tropical atmosphere, as a consequence of the smallness of the Coriolis parameter near the equator. This provides a strong constraint on both large-scale fluid dynamics and diabatic processes. This work is a step toward the construction of a balanced dynamical theory for the tropical circulation that is based on this constraint, and in which the diabatic processes are explicit and interactive. The authors first derive the basic fluid-dynamical scaling under the weak temperature gradient (WTG) approximation in a shallow water system with a fixed mass source representing an externally imposed heating. This derivation follows an earlier similar one by Held and Hoskins, but extends the analysis to the nonlinear case (though on an f plane), examines the resulting system in more detail, and presents a solution for an axisymmetric ''top-hat'' forcing. The system is truly balanced, having no gravity waves, but is different from other balance models in that the heating is included a priori in the scaling. The WTG scaling is then applied to a linear moist model in which the convective heating is controlled by a moisture variable that is advected by the flow. This moist model is derived from the Quasi-equilibrium Tropical Circulation Model (QTCM) equations of Neelin and Zeng but can be viewed as somewhat more general. A number of additional approximations are made in order to consider balanced dynamical modes, apparently not studied previously, which owe their existence to interactions of the moisture and flow fields. A particularly interesting mode arises on an f plane with a constant background moisture gradient. In the limit of low frequency and zero meridional wavenumber this mode has a dispersion relation mathematically identical to that of a barotropic Rossby wave, though the phase speed is eastward (for moisture decreasing poleward in the background state) and the propagation mechanism is quite different. This mode also has significant positive growth rate for low wavenumbers. The addition of the ␤ effect complicates matters. For typical parameters, when ␤ is included the direction of phase propagation is ambiguous, and the growth rate reduced, as the effects of the background gradients in moisture and planetary vorticity appear to cancel to a large degree. Possible relevance to intraseasonal variability and easterly wave dynamics is briefly discussed. *Lamont-Doherty Earth Observatory Contribution Number LDEO 6190.

It has recently been shown that extreme stratospheric events (ESEs) are followed by surface weather anomalies (for up to 60 days), suggesting that stratospheric variability might be used to extend weather prediction beyond current time scales. In this paper, attention is drawn away from the stratosphere to demonstrate that the originating point of ESEs is located in the troposphere. First, it is shown that anomalously strong eddy heat fluxes at 100 hPa nearly always precede weak vortex events, and conversely, anomalously weak eddy heat fluxes precede strong vortex events, consistent with wave-mean flow interaction theory. This finding clarifies the dynamical nature of ESEs and suggests that a major source of stratospheric variability (and thus predictability) is located in the troposphere below and not in the stratosphere itself. Second, it is shown that the daily time series of eddy heat flux found at 100 hPa and integrated over the prior 40 days, exhibit a remarkably high anticorrelation (0.8) with the Arctic Oscillation (AO) index at 10 hPa. Following Baldwin and Dunkerton, it is then demonstrated that events with anomalously strong (weak) integrated eddy heat fluxes at 100 hPa are followed by anomalously large (small) surface values of the AO index up to 60 days following each event. This suggests that the stratosphere is unlikely to be the dominant source of the anomalous surface weather regimes discussed in Thompson et al.