Are the winters 2010 and 2012 archetypes exhibiting extreme opposite behavior of the North Atlantic jet stream?

The atmospheric circulation over the North Atlantic-European sector experienced exceptional but highly contrasting conditions in the recent 2010 and 2012 winters (November-March; with the year dated by the relevant January). Evidence is given for the remarkably different locations of the eddy-driven westerly jet over the North Atlantic. In the 2010 winter the maximum of the jet stream was systematically between 30ºN and 40ºN (in the ‘south jet regime’), while in the 2012 winter it was predominantly located around 55ºN (north jet regime). These jet features underline the occurrence of either weak flow (2010) or strong and persistent ridges throughout the troposphere (2012). This is confirmed by the very different occurrence of blocking systems over the North Atlantic, associated with episodes of strong cyclonic (anticyclonic) Rossby wave breaking in 2010 (2012) winters. These dynamical features underlie strong precipitation and temperature anomalies over parts of Europe, with detrimental impacts on many socioeconomic sectors. Despite the highly contrasting atmospheric states, mid and high-latitude boundary conditions do not reveal strong differences in these two winters. The two winters were associated with opposite ENSO phases, but there is no causal evidence of a remote forcing from the Pacific sea surface temperatures. Finally, the exceptionality of the two winters is demonstrated in relation to the last 140 years. It is suggested that these winters may be seen as archetypes of North Atlantic jet variability under current climate conditions.

Dynamical evolution of North Atlantic ridges and poleward Jet Stream displacements

The development of a particular wintertime atmospheric circulation regime over the North Atlantic, comprising a northward shift of the North Atlantic eddy-driven jet stream and an associated strong and persistent ridge in the subtropics, is investigated. Several different methods of analysis are combined to describe the temporal evolution of the events and relate it to shifts in the phase of the North Atlantic Oscillation and East Atlantic pattern. First, the authors identify a close relationship between northward shifts of the eddy-driven jet, the establishment and maintenance of strong and persistent ridges in the subtropics, and the occurrence of upper-tropospheric anticyclonic Rossby wave breaking over Iberia. Clear tropospheric precursors are evident prior to the development of the regime, suggesting a preconditioning of the Atlantic jet stream and an upstream influence via a large-scale Rossby wave train from the North Pacific. Transient (2–6 days) eddy forcing plays a dual role, contributing to both the initiation and then the maintenance of the circulation anomalies. During the regime there is enhanced occurrence of anticyclonic Rossby wave breaking, which may be described as low-latitude blocking-like events over the southeastern North Atlantic. A strong ridge is already established at the time of wave-breaking onset, suggesting that the role of wave-breaking events is to amplify the circulation anomalies rather than to initiate them. Wave breaking also seems to enhance the persistence, since it is unlikely that a persistent ridge event occurs without being also accompanied by wave breaking.

Changing Northern Hemisphere storm tracks in an ensemble of IPCC climate change simulations

Winter storm-track activity over the Northern Hemisphere and its changes in a greenhouse gas scenario (the Special Report on Emission Scenarios A1B forcing) are computed from an ensemble of 23 single runs from 16 coupled global climate models (CGCMs). All models reproduce the general structures of the observed climatological storm-track pattern under present-day forcing conditions. Ensemble mean changes resulting from anthropogenic forcing include an increase of baroclinic wave activity over the eastern North Atlantic, amounting to 5%–8% by the end of the twenty-first century. Enhanced activity is also found over the Asian continent and over the North Pacific near the Aleutian Islands. At high latitudes and over parts of the subtropics, activity is reduced. Variations of the individual models around the ensemble average signal are not small, with a median of the pattern correlation near r = 0.5. There is, however, no evidence for a link between deviations in present-day climatology and deviations with respect to climate change.

The GCM–Reality Intercomparison Project for SPARC (GRIPS): scientific issues and initial results

To investigate the effects of the middle atmosphere on climate, the World Climate Research Programme is supporting the project "Stratospheric Processes and their Role in Climate" (SPARC). A central theme of SPARC, to examine model simulations of the coupled troposphere—middle atmosphere system, is being performed through the initiative called GRIPS (GCM—Reality Intercomparison Project for SPARC). In this paper, an overview of the objectives of GRIPS is given. Initial activities include an assessment of the performance of middle atmosphere climate models, and preliminary results from this evaluation are presented here. It is shown that although all 13 models evaluated represent most major features of the mean atmospheric state, there are deficiencies in the magnitude and location of the features, which cannot easily be traced to the formulation (resolution or the parameterizations included) of the models. Most models show a cold bias in all locations, apart from the tropical tropopause region where they can be either too warm or too cold. The strengths and locations of the major jets are often misrepresented in the models. Looking at three—dimensional fields reveals, for some models, more severe deficiencies in the magnitude and positioning of the dominant structures (such as the Aleutian high in the stratosphere), although undersampling might explain some of these differences from observations. All the models have shortcomings in their simulations of the present—day climate, which might limit the accuracy of predictions of the climate response to ozone change and other anomalous forcing.

Mechanisms for tropical upwelling in the stratosphere

The dynamics of the tropical upwelling branch of the stratospheric Brewer–Dobson circulation are examined, with a particular focus on the role of the middle-atmosphere Hadley circulation. Upwelling is examined in terms of both the diabatic circulation and Lagrangian trajectories using a zonally symmetric balance model. The behavior of the wave-driven circulation in the presence of angular momentum redistribution by the Hadley circulation is also considered. The results of the zonally symmetric model are compared with fields from a middle-atmosphere GCM. It is found that the Hadley circulation makes a significant contribution to annual mean tropical upwelling at the upwelling maximum in the vicinity of the stratopause, and can account for most of the annual mean upwelling seen in the GCM. In the mid- to lower stratosphere, the role of the Hadley circulation is much weaker and wave drag appears to be required to explain the observed upwelling, although the Hadley circulation makes a nonnegligible contribution to the annual cycle of the upwelling. Subtropical wave drag can produce annual mean upwelling through a nonlinear mechanism; viscosity is not required. However, the magnitude of the observed upwelling suggests that wave drag must penetrate quite close to the equator.

The middle-atmosphere Hadley circulation and equatorial inertial adjustment

In the tropical middle atmosphere the climatological radiative equilibrium temperature is inconsistent with gradient-wind balance and the available angular momentum, especially during solstice seasons. Adjustment toward a balanced state results in a type of Hadley circulation that lies outside the “downward control” view of zonally averaged dynamics. This middle-atmosphere Hadley circulation is reexamined here using a zonally symmetric balance model driven through an annual cycle. It is found that the inclusion of a realistic radiation scheme leads to a concentration of the circulation near the stratopause and to its closing off in the mesosphere, with no need for relaxational damping or a rigid lid. The evolving zonal flow is inertially unstable, leading to a rapid process of inertial adjustment, which becomes significant in the mesosphere. This short-circuits the slower process of angular momentum homogenization by the Hadley circulation itself, thereby weakening the latter. The effect of the meridional circulation associated with extratropical wave drag on the Hadley circulation is considered. It is shown that the two circulations are independent for linear (quasigeostrophic) zonal-mean dynamics, and interact primarily through the advection of temperature and angular momentum. There appears to be no significant coupling in the deep Tropics via temperature advection since the wave-driven circulation is unable to alter meridional temperature gradients in this region. However, the wave-driven circulation can affect the Hadley circulation by advecting angular momentum out of the Tropics. The validity of the zonally symmetric balance model with parameterized inertial adjustment is tested by comparison with a three-dimensional primitive equations model. Fields from a middle-atmosphere GCM are also examined for evidence of these processes. While many aspects of the GCM circulation are indicative of the middle-atmosphere Hadley circulation, particularly in the upper stratosphere, it appears that the circulation is obscured in the mesosphere and lower stratosphere by other processes.

A closer look at chaotic advection in the stratosphere: part I: geometric structure

The relevance of chaotic advection to stratospheric mixing and transport is addressed in the context of (i) a numerical model of forced shallow-water flow on the sphere, and (ii) a middle-atmosphere general circulation model. It is argued that chaotic advection applies to both these models if there is suitable large-scale spatial structure in the velocity field and if the velocity field is temporally quasi-regular. This spatial structure is manifested in the form of “cat’s eyes” in the surf zone, such as are commonly seen in numerical simulations of Rossby wave critical layers; by analogy with the heteroclinic structure of a temporally aperiodic chaotic system the cat’s eyes may be thought of as an “organizing structure” for mixing and transport in the surf zone. When this organizing structure exists, Eulerian and Lagrangian autocorrelations of the velocity derivatives indicate that velocity derivatives decorrelate more rapidly along particle trajectories than at fixed spatial locations (i.e., the velocity field is temporally quasi-regular). This phenomenon is referred to as Lagrangian random strain.

A closer look at chaotic advection in the stratosphere: part II: statistical diagnostics

Statistical diagnostics of mixing and transport are computed for a numerical model of forced shallow-water flow on the sphere and a middle-atmosphere general circulation model. In particular, particle dispersion statistics, transport fluxes, Liapunov exponents (probability density functions and ensemble averages), and tracer concentration statistics are considered. It is shown that the behavior of the diagnostics is in accord with that of kinematic chaotic advection models so long as stochasticity is sufficiently weak. Comparisons with random-strain theory are made.

Coupled Kelvin-wave and mirage-wave instabilities in semigeostrophic dynamics

A weak instability mode, associated with phase-locked counterpropagating coastal Kelvin waves in horizontal anticyclonic shear, is found in the semigeostrophic (SG) equations for stratified flow in a channel. This SG instability mode approximates a similar mode found in the Euler equations in the limit in which particle-trajectory slopes are much smaller than f/N, where f is the Coriolis frequency and N > f the buoyancy frequency. Though weak under normal parameter conditions, this instability mode is of theoretical interest because its existence accounts for the failure of an Arnol’d-type stability theorem for the SG equations. In the opposite limit, in which the particle motion is purely vertical, the Euler equations allow only buoyancy oscillations with no horizontal coupling. The SG equations, on the other hand, allow a physically spurious coastal “mirage wave,” so called because its velocity field vanishes despite a nonvanishing disturbance pressure field. Counterpropagating pairs of these waves can phase-lock to form a spurious “mirage-wave instability.” Closer examination shows that the mirage wave arises from failure of the SG approximations to be self-consistent for trajectory slopes f/N.

On the group-velocity property for wave-activity conservation laws

The density and the flux of wave-activity conservation laws are generally required to satisfy the group-velocity property: under the WKB approximation (i.e., for nearly monochromatic small-amplitude waves in a slowly varying medium), the flux divided by the density equals the group velocity. It is shown that this property is automatically satisfied if, under the WKB approximation, the only source of rapid variations in the density and the flux lies in the wave phase. A particular form of the density, based on a self-adjoint operator, is proposed as a systematic choice for a density verifying this condition.

On nonlinear symmetric stability and the nonlinear saturation of symmetric instability

A nonlinear symmetric stability theorem is derived in the context of the f-plane Boussinesq equations, recovering an earlier result of Xu within a more general framework. The theorem applies to symmetric disturbances to a baroclinic basic flow, the disturbances having arbitrary structure and magnitude. The criteria for nonlinear stability are virtually identical to those for linear stability. As in Xu, the nonlinear stability theorem can be used to obtain rigorous upper bounds on the saturation amplitude of symmetric instabilities. In a simple example, the bounds are found to compare favorably with heuristic parcel-based estimates in both the hydrostatic and non-hydrostatic limits.

On Hamiltonian balanced dynamics and the slowest invariant manifold

The concept of a slowest invariant manifold is investigated for the five-component model of Lorenz under conservative dynamics. It is shown that Lorenz's model is a two-degree-of-freedom canonical Hamiltonian system, consisting of a nonlinear vorticity-triad oscillator coupled to a linear gravity wave oscillator, whose solutions consist of regular and chaotic orbits. When either the Rossby number or the rotational Froude number is small, there is a formal separation of timescales, and one can speak of fast and slow motion. In the same regime, the coupling is weak, and the Kolmogorov–Arnold-Moser theorem is shown to apply. The chaotic orbits are inherently unbalanced and are confined to regions sandwiched between invariant tori consisting of quasi-periodic regular orbits. The regular orbits generally contain free fast motion, but a slowest invariant manifold may be geometrically defined as the set of all slow cores of invariant tori (defined by zero fast action) that are smoothly related to such cores in the uncoupled system. This slowest invariant manifold is not global; in fact, its structure is fractal; but it is of nearly full measure in the limit of weak coupling. It is also nonlinearly stable. As the coupling increases, the slowest invariant manifold shrinks until it disappears altogether. The results clarify previous definitions of a slowest invariant manifold and highlight the ambiguity in the definition of “slowness.” An asymptotic procedure, analogous to standard initialization techniques, is found to yield nonzero free fast motion even when the core solutions contain none. A hierarchy of Hamiltonian balanced models preserving the symmetries in the original low-order model is formulated; these models are compared with classic balanced models, asymptotically initialized solutions of the full system and the slowest invariant manifold defined by the core solutions. The analysis suggests that for sufficiently small Rossby or rotational Froude numbers, a stable slowest invariant manifold can be defined for this system, which has zero free gravity wave activity, but it cannot be defined everywhere. The implications of the results for more complex systems are discussed.

Nonlinear stability of Eady's model

A nonlinear stability theorem is established for Eady's model of baroclinic flow. In particular, the Eady basic state is shown to be nonlinearly stable (for arbitrary shear) provided (Δz)/(Δy) > 2(5)^1/2f/(πN),where Δz is the height of the domain, Δy the channel width, f the Coriolis parameter, and N the buoyancy frequency. When this criterion is satisfied, explicit bounds can be derived on the disturbance potential enstrophy, the disturbance energy, and the disturbance available potential energy on the rigid lids, which are expressed in terms of the initial disturbance fields. The disturbances are completely general (with nonzero potential vorticity) and are not assumed to be of small amplitude. The results may be regarded as an extension of Arnol'd's second nonlinear stability theorem to continuously stratified quasigeostrophic baroclinic flow.

Nonlinear saturation of baroclinic instability: part III: bounds on the energy

Rigorous upper bounds are derived on the saturation amplitude of baroclinic instability in the two-layer model. The bounds apply to the eddy energy and are obtained by appealing to a finite amplitude conservation law for the disturbance pseudoenergy. These bounds are to be distinguished from those derived in Part I of this study, which employed a pseudomomentum conservation law and provided bounds on the eddy potential enstrophy. The bounds apply to conservative (inviscid, unforced) flow, as well as to forced-dissipative flow when the dissipation is proportional to the potential vorticity. Bounds on the eddy energy are worked out for a general class of unstable westerly jets. In the special case of the Phillips model of baroclinic instability, and in the limit of infinitesimal initial eddy amplitude, the bound states that the eddy energy cannot exceed ϵβ2/6F where ϵ = (U − Ucrit)/Ucrit is the relative supercriticality. This bound captures the essential dynamical scalings (i.e., the dependence on ϵ, β, and F) of the saturation amplitudes predicted by weakly nonlinear theory, as well as exhibiting remarkable quantitative agreement with those predictions, and is also consistent with heuristic baroclinic adjustment estimates.