## 1. Introduction

Stratospheric polar vortex in the Arctic winter is significantly more volatile than its Antarctic counterpart, the culmination of which being stratospheric sudden warmings (SSWs). Recurring every 1.5–2 winters on average, SSWs involve a rapid transformation of the polar vortex, reversal of the circumpolar winds and a sudden increase in the polar temperatures (Scherhag 1952; Quiroz 1975; McIntyre and Palmer 1983, 1984; Limpasuvan et al. 2004; Charlton and Polvani 2007; Butler et al. 2015). SSWs not only disrupt stratospheric circulation but also affect tropospheric circulation of the extratropics in weeks to follow (Kodera and Chiba 1995; Baldwin and Dunkerton 1999, 2001; Kushner and Polvani 2004; Polvani and Waugh 2004; Gerber et al. 2009; Hitchcock and Simpson 2014).

SSWs are arguably the most dramatic manifestation of wave–mean flow interaction in Earth’s atmosphere, for which transience in the planetary-scale Rossby waves plays a central role (McIntyre 1982; Limpasuvan et al. 2004). In the first successful modeling of SSW, Matsuno (1971) assumes that wave transience is a forced response of the stratosphere to tropospheric processes such as atmospheric blocking, a paradigm adopted by many subsequent studies (e.g., Holton 1976; Polvani and Saravanan 2000; Scott and Haynes 2000; Sjoberg and Birner 2012), with substantial observational support (Randel and Boville 1987; Polvani and Waugh 2004; Martius et al. 2009; Woollings et al. 2010). An impulsive wave forcing at the tropopause produces a vertically propagating packet of quasi-stationary Rossby waves, and the converging Eliassen–Palm (E–P) flux at the leading edge of the wave packet decelerates the zonal-mean zonal wind (Matsuno 1971).

An alternative paradigm of SSW asserts that wave transience arises from the internal dynamics of the stratosphere. Holton and Mass (1976) demonstrate with an idealized model that even a steady wave forcing at the tropopause generates stratospheric vacillation in wave activity and the zonal-mean zonal wind when a critical forcing amplitude is exceeded. Plumb (1981) proposes instability due to interaction between the steady wave forcing and near-resonant free waves in the stratosphere as a cause of wave transience (see also Geisler 1974; Tung and Lindzen 1979; Smith 1989). Unlike the transient forcing paradigm, the resonance theory predicts the amplification of the stationary wave as part of the solution. Resonance and wave–mean flow interaction also produce multiple equilibria (cf. Charney and DeVore 1979; Held 1983), prompting some to characterize SSW as a transition from a “high-index” (nearly zonal) state to a “low-index” (perpetual resonance) state of the polar vortex (Chao 1985; Yoden 1987, 1990; Ruzmaikin et al. 2003; Birner and Williams 2008; Yasuda et al. 2017). The resonance idea has been applied to nearly barotropic, split-type SSWs (Esler and Scott 2005; Matthewman and Esler 2011; Liu and Scott 2015).

Whether one adopts the first paradigm or the second, modeling studies have long shown that the properties of wave transience are sensitive to the initial profile of the stratospheric winds and their vertical shear (Holton and Dunkerton 1978; Chen and Robinson 1992; Scott and Polvani 2006). This suggests that the condition of the stratosphere has a strong filtering effect on the forcing from below (Hitchcock and Haynes 2016). Generally, weak (but not too weak) zonal winds are conducive to SSWs, but such “preconditioning” of stratospheric winds in turn depends on prior wave events (or lack thereof) (Polvani and Waugh 2004; Jucker 2016; de la Cámara et al. 2017).

Although the short-term predictability of SSWs has improved over time (Mukougawa and Hirooka 2004; Stan and Straus 2009; Tripathi et al. 2016; Taguchi 2016; Rao et al. 2018), there remains considerable uncertainty in GCMs’ ability to predict the frequency of SSWs and its trend under a changing climate (Labitzke and Naujokat 2000; Charlton-Perez et al. 2008; McLandress and Shepherd 2009; Ayarzagüena et al. 2018). This may be due partly to variations in the definition of SSW (Butler et al. 2015), but more fundamentally, our understanding of SSW’s onset is still incomplete, limiting our ability to interpret the disparate model results. Stationary waves in the stratosphere sometimes attain large amplitude and yet do not produce a full-blown SSW; other times they do (Solomon 2014). When an SSW does occur, deceleration and reversal of the zonal-mean zonal wind proceed swiftly, exhibiting a sense of “suddenness.” What separates SSW and non-SSW conditions and why is vortex breakdown abrupt? The resonance paradigm partially addresses this question but supportive evidence from observation is still circumstantial at best (e.g., Smith 1989; Esler et al. 2006).

In what follows, we propose a view that the timing and suddenness of SSWs are determined, on average, by a threshold behavior of Rossby waves after they attain finite-amplitude wave activity flux in the vertical, regardless of how they attain it. Based on the observed finite-amplitude wave activity [FAWA; Nakamura and Zhu (2010), hereafter NZ10, Nakamura and Solomon (2010), hereafter NS10] and its vertical flux, we present evidence that there is an upper bound on the upward E–P flux (“transmission capacity”) of a stationary Rossby wave for a given altitude and flow condition. A rapid, spontaneous vortex breakdown occurs once the incident E–P flux from below reaches this capacity. As we will see, this happens when the zonal-mean zonal wind *u*_{REF} (to be defined below). We thus propose the ratio

The dynamics underlying the threshold behavior is akin to that discussed by Wang and Fyfe (2000), who use a Boussinesq contour dynamics model of the polar vortex to show that the onset of wave breaking occurs once the zonal-mean zonal wind drops below about one-half of the initial value. However, unlike Wang and Fyfe, our diagnostic applies to instantaneous data without a need for solving an initial-value problem. We will also construct an idealized 1D model of SSWs that encapsulates the essence of the threshold behavior. Mathematically similar to the “traffic jam” model of Nakamura and Huang (2018) for atmospheric blocking, this model predicts salient features of SSWs including a rapid vortex breakdown and a downward migration of wave activity/zonal wind anomalies that follows.

In the next section we will review the diagnostic formalism and describe wave–mean flow interaction observed during the life cycles of SSWs in terms of FAWA, *u*_{REF}. Section 3 explains how wave–mean flow interaction makes the vertical component of E–P flux a nonlinear function of FAWA and how this leads to a threshold behavior. We will demonstrate the threshold behavior in response to transient wave forcing with an idealized, semiempirical 1D model of SSW based on the WKB theory. In section 4 we will examine the effects of radiative damping and mixing on the threshold behavior. In section 5 we present an observational support for the theory and demonstrate the utility of

## 2. Wave–mean flow interaction during SSW life cycles

### a. FAWA and reference state

*ϕ*

_{e}[e.g., Eq. (4) of NS10]:

*q*(

*λ*,

*ϕ*,

*z*,

*t*) is PV (see Table 1 for the list of symbols) and

*dS*=

*a*

^{2}cos

*ϕdλdϕ*. Domains

*D*

_{1}and

*D*

_{2}and equivalent latitude

*ϕ*

_{e}are defined by

*Q*(

*ϕ*

_{e},

*z*,

*t*) specifies the PV contour that encloses the same area as the polar cap bounded by

*ϕ*=

*ϕ*

_{e}, that is, the areas of

*D*

_{1}and

*D*

_{2}are the same [see Fig. 1 of NS10]. By definition FAWA

*A*(

*ϕ*

_{e},

*z*,

*t*) is nonnegative and vanishes only when the PV contour is zonally symmetric.

List of symbols.

*A*(

*ϕ*

_{e},

*z*,

*t*) are related through the following relation, assuming

*ϕ*

_{e}=

*ϕ*[e.g., Eq. (10) of NS10 and Eq. (3) of Nakamura and Solomon (2011), hereafter NS11]:

*μ*≡ sin

*ϕ*and

*S*,

*A*,

*A*/∂

*t*= 0),

*reference-state wind u*

_{REF}(

*ϕ*,

*z*,

*t*) that satisfies the following:

*u*

_{REF}is related to the reference-state PV

*q*

_{REF}(

*ϕ*,

*z*,

*t*), which is obtained by zonalizing the instantaneous PV (

*q*) through area-preserving map on each

*z*surface (NZ10, NS10):

^{1}Since

*q*

_{REF}changes only in response to nonconservative processes,

*u*

_{REF}is invariant under conservative dynamics [Eq. (6), see Methven and Berrisford (2015) for a related discussion]. In what follows,

*u*

_{REF}is inverted hemispherically from

*q*

_{REF}with Eq. (7) with no-slip lower boundary condition [see supplementary material of Nakamura and Huang (2018)].

### b. Reversible and irreversible components of wave–mean flow interaction

*A*= 0) one obtains

*u*

_{REF},

*A*. Their relationship at 60°N and

*z*= 32 km (10.34 hPa) during SSW’s life cycle is demonstrated in Fig. 1a for the 2009 event (Harada et al. 2010) [data source: ERA-Interim (Dee et al. 2011)]. Vortex conditions at four stages of this splitting event, labeled A–D, are illustrated in Fig. 2. As the vortex is distorted,

*A*increases and Δ

*u*decreases, following a nearly linear relationship (Fig. 1a, red dots). Note that Δ

*u*is always negative, meaning that

*u*

_{REF}and it becomes even more so as the event matures. As the flow recovers after the event (blue dots),

*A*decreases and Δ

*u*increases, tracing closely the linear path of the lead. Both

*A*and Δ

*u*almost return to their original values. Figure 1b shows the same diagram for a composite of 18 SSW events between 1979 and 2016. (To obtain a typical picture of SSWs, we have excluded a few late-season events and those that are preceded by another significant event during the 30-day lead time. See Table 2 for the list of the events.) Despite the broad type variations from one event to another (Table 2), and despite the visually evident irreversibility in Fig. 2, the reversible, linear relationship between

*A*and Δ

*u*is robust through the life cycles of SSWs.

Potential vorticity (color shading; PVU; 1 PVU = 10^{−6} K kg^{−1} m^{2} s^{−1}) and Montgomery streamfunction (contours; m^{2} s^{−2}) on 850-K isentropic surface (approximately the same altitude as in Fig. 1) during the 2009 SSW: (a) 8 Jan; (b) 19 Jan; (c) 28 Jan; (d) 5 Feb. The four panels correspond to the four stages labeled A–D in Figs. 1a and 1c.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Potential vorticity (color shading; PVU; 1 PVU = 10^{−6} K kg^{−1} m^{2} s^{−1}) and Montgomery streamfunction (contours; m^{2} s^{−2}) on 850-K isentropic surface (approximately the same altitude as in Fig. 1) during the 2009 SSW: (a) 8 Jan; (b) 19 Jan; (c) 28 Jan; (d) 5 Feb. The four panels correspond to the four stages labeled A–D in Figs. 1a and 1c.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Potential vorticity (color shading; PVU; 1 PVU = 10^{−6} K kg^{−1} m^{2} s^{−1}) and Montgomery streamfunction (contours; m^{2} s^{−2}) on 850-K isentropic surface (approximately the same altitude as in Fig. 1) during the 2009 SSW: (a) 8 Jan; (b) 19 Jan; (c) 28 Jan; (d) 5 Feb. The four panels correspond to the four stages labeled A–D in Figs. 1a and 1c.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

List of SSWs included and excluded in Figs. 1a, 1b, and 11. SSW type categories are based on Charlton and Polvani (2007) and Kodera et al. (2016).

*A*/

*u*

_{REF}during 60-day leads of SSW events on different pressure levels, using both ERA-Interim and MERRA-2 (Gelaro et al. 2017). Data points in both panels cluster around a linear relationship expressed approximately

^{2}

*α*≈ 0.4 means that about 40% of FAWA is used to drive

*u*

_{REF}, whereas the remainder is used to drive the residual circulation (NS10). Although the local relation [Eq. (9)] cannot be expected directly from Eq. (8) (except for barotropic case in which

*α*= 1), given the robust observational support, we will use Eq. (9) as a key empirical constraint in the subsequent theoretical development.

Relationship between *A*/*u*_{REF} (abscissa) and *ϕ* = 60°N during 60 days leading to SSW events. Red: 10 hPa. Orange: 20 hPa. Purple: 30 hPa. Green: 50 hPa. To emphasize the average structure, error bars are limited to ±0.25 standard deviations. (a) ERA-Interim (Dee et al. 2011) (daily composite of 26 cases). (b): MERRA-2 (Gelaro et al. 2017) (daily composite of 25 cases). The event generally progresses from upper left to lower right. See Table 1 of Lubis et al. (2018a) for the dates of SSWs.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Relationship between *A*/*u*_{REF} (abscissa) and *ϕ* = 60°N during 60 days leading to SSW events. Red: 10 hPa. Orange: 20 hPa. Purple: 30 hPa. Green: 50 hPa. To emphasize the average structure, error bars are limited to ±0.25 standard deviations. (a) ERA-Interim (Dee et al. 2011) (daily composite of 26 cases). (b): MERRA-2 (Gelaro et al. 2017) (daily composite of 25 cases). The event generally progresses from upper left to lower right. See Table 1 of Lubis et al. (2018a) for the dates of SSWs.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Relationship between *A*/*u*_{REF} (abscissa) and *ϕ* = 60°N during 60 days leading to SSW events. Red: 10 hPa. Orange: 20 hPa. Purple: 30 hPa. Green: 50 hPa. To emphasize the average structure, error bars are limited to ±0.25 standard deviations. (a) ERA-Interim (Dee et al. 2011) (daily composite of 26 cases). (b): MERRA-2 (Gelaro et al. 2017) (daily composite of 25 cases). The event generally progresses from upper left to lower right. See Table 1 of Lubis et al. (2018a) for the dates of SSWs.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

In contrast to the reversible nature of Δ*u*, the behavior of *u*_{REF} during SSW is distinctly irreversible. Figures 1c and 1d plot *u*_{REF} cos*ϕ* against *A* cos*ϕ* for the corresponding periods in Figs. 1a and 1b at 60°N and *z* = 32 km. During the 30-day lead, *u*_{REF} remains more or less steady for small values of *A*. For larger values of *A*, *u*_{REF} starts to decrease steadily and it does not recover even after *A* reverses its course. The steady decline of *u*_{REF} is driven by nonconservative processes [Eq. (6)], and although part of this driving is seasonal radiative forcing, the majority of it comes from mixing (Lubis et al. 2018a,b; Martineau and Son 2015). The irreversible loss of *u*_{REF} through mixing delays the vortex recovery in some SSW events (Fig. 9 of Lubis et al. 2018a).

## 3. Theory of threshold behavior

In the previous section we described how finite-amplitude Rossby waves modify the zonal-mean zonal wind during the life cycles of SSWs. In this section we discuss how this affects the propagation of Rossby waves and gives rise to a threshold behavior, using a semiempirical theory.^{3}

### a. Vertical group velocity of a stationary Rossby wave

^{4}The same theory also governs FAWA at small amplitude under conservative dynamics (NZ10). To maximize the affinity of quasigeostrophic theory for data, we formulate the dispersion relation for a near plane wave in the Mercator coordinate (e.g., section 5b of Hoskins and Karoly 1981). This leads to the following expression for the vertical group velocity for the stationary Rossby wave (see appendix for derivation):

*k*=

*n*/

*a*is equatorial wavenumber and

*n*> 0 is the integer zonal wavenumber, whereas

*m*is vertical wavenumber given by Eq. (A9). We assume a positive root of

*m*and

*N*nearly constant with height in the stratosphere (≈2.1 × 10

^{−2}s

^{−1}). The factor

*N*, this is the only other empirical constraint we use to simplify the theory.

Relationship between *ϕ* = 60°N during 60 days leading up to SSW events. Red: 10 hPa. Orange: 20 hPa. Purple: 30 hPa. Green: 50 hPa. Error bars indicate ±0.25 standard deviation. (a) ERA-Interim (daily composite of 26 cases). (b) MERRA2 (daily composite of 25 cases). The event generally progresses from upper right to lower left.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Relationship between *ϕ* = 60°N during 60 days leading up to SSW events. Red: 10 hPa. Orange: 20 hPa. Purple: 30 hPa. Green: 50 hPa. Error bars indicate ±0.25 standard deviation. (a) ERA-Interim (daily composite of 26 cases). (b) MERRA2 (daily composite of 25 cases). The event generally progresses from upper right to lower left.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Relationship between *ϕ* = 60°N during 60 days leading up to SSW events. Red: 10 hPa. Orange: 20 hPa. Purple: 30 hPa. Green: 50 hPa. Error bars indicate ±0.25 standard deviation. (a) ERA-Interim (daily composite of 26 cases). (b) MERRA2 (daily composite of 25 cases). The event generally progresses from upper right to lower left.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

*l*[and through Eq. (A9),

*m*] as follows. Invariance of Eq. (11) through wave–mean flow interaction suggests

*aμ*),

^{5}the lhs of Eq. (12) may be rewritten as

*A*with respect to latitude:

*α*= 0.4 in the last approximation. This gives the meridional scale of FAWA,

*L*≈ 1.45 × 10

^{6}m, at 60°N. Then assuming that

*q*′ is modal in latitude and proportional to sin

*ly*[

*y*is Mercator latitude; Eq. (A1)]:

^{6}

*m*= 2.36 × 10

^{−4}m

^{−1}for

*n*= 1 and 2.18 × 10

^{−4}m

^{−1}for

*n*= 2. Substituting Eq. (15) into Eq. (10), we obtain the vertical group velocity of the stationary Rossby wave as

*C*is the ratio of

*c*

_{gz}to

*C*is nearly independent of height and

*C*≈ 2.23 × 10

^{−3}for

*n*= 1 and ≈4.11 × 10

^{−3}for

*n*= 2 at 60°N. Equation (17) highlights the proportionality between the group velocity and

*A*according to Eq. (9),

### b. Wave activity budget and threshold

^{7}

*c*

_{gz}

*Ae*

^{−z/H}≡

*F*

_{z}is the vertical component of wave activity (E–P) flux density. Here the convergence of the meridional wave activity flux and nonconservative terms are ignored. These omissions would make Eq. (19) a poor approximation for describing the full life cycle of SSW (e.g., O’Neill and Taylor 1979; Lubis et al. 2018a,b). Yet there is mounting evidence that early development of SSW is dictated by the upward E–P flux and its convergence (Polvani and Waugh 2004; Jucker 2016), so we take Eq. (19) as a model for the SSW’s onset. [Alternatively, Eq. (19) may be thought of as the budget of meridionally integrated wave activity (e.g., Plumb and Semeniuk 2003).] With Eq. (18) we have

*c*

_{gz}is a linear function of

*A*,

*F*

_{z}is a quadratic function of

*A*. We will compare this theoretical prediction of

*F*

_{z}with observation in section 5. For now we demonstrate how the nonlinear flux of FAWA brings about a threshold behavior.

*F*

_{0}> 0 for upward propagation):

*F*

_{0}may be determined, for example, by a steady wave forcing at the lower boundary. Solving Eq. (22) for a steady-state

*ξ*:

*ξ*is unstable and hence unrealizable (we will see why shortly). The discriminant condition in Eq. (23) means that there is a height-dependent upper bound on the wave activity flux that a steady stationary wave can transmit. This places a “transmission capacity” for finite-amplitude Rossby waves, whose vertical transmission is guaranteed at small amplitude (Charney and Drazin 1961). For a given

*F*

_{0}and Γ(

*z*), Eq. (23) predicts the altitude at which this upper bound is first reached. Because of the exponential factor, Γ(

*z*) is a decreasing function of

*z*at sufficiently high altitude [Eq. (21)], and therefore even if

*F*

_{0}< Γ(

*z*)/4 at low altitudes,

*F*

_{0}eventually equals Γ/4 at some

*z*. This altitude

*z*

_{c}is lower for larger

*F*

_{0}. At this altitude, the following conditions are met:

*u*

_{REF}at

*z*=

*z*

_{c}, and this criterion does not depend on the value of

*α*as long as it is positive. The above conditions correspond to the vertex of the quadratic curve in Fig. 5 (point “c”) touching the line

*F*

_{z}=

*F*

_{0}. For

*z*>

*z*

_{c}, the vertex will be lower than

*F*

_{0}, so Eq. (22) has no real root. To ensure a real value for

*ξ*,

*F*

_{z}must be capped at the transmission capacity:

*F*

_{z}/∂

*z*< 0 above

*z*

_{c}(see the left panel of Fig. 5b), which is incompatible with the steady-state assumption:

*A*(and hence

*ξ*) must increase with time according to Eq. (19). From Eqs. (19), (25), and (21), the growth rate of wave activity is estimated to be

*e*-folding time of 1.5–3 days, consistent with the suddenness of SSW. The growth of

*ξ*beyond

*ξ*

_{c}causes

*F*

_{z}to decrease further, since

*F*

_{z}is a decreasing function of

*ξ*for

*ξ*>

*ξ*

_{c}[Eq. (21), Fig. 5a]. This sets up a positive feedback between wave activity and its flux: the decreasing flux above the threshold altitude

*z*

_{c}reinforces the flux convergence and further increase in

*A*, which then promotes further reduction in

*F*

_{z}, provided that there is a continued supply of incident flux from below. Therefore, above

*z*=

*z*

_{c}where the threshold (24) is first met, wave activity increases rapidly until it reaches

*ξ*>

*ξ*

_{c}and corresponds to “vortex breakdown.”

^{8}The threshold behavior and spontaneous breakdown of the vortex arises only because the E–P flux is a nonlinear function of

*ξ*and unimodal.

(a) Schematic diagram for the relationship between FAWA density (*Ae*^{−z/H}; horizontal axis) and E–P flux (*F*_{z}; vertical axis) at a given altitude. The corresponding values of *ξ* are also labeled on the horizontal axis. For *ξ* < *ξ*_{c} = 0.5 (“a” and “b”), a stable steady state is possible. Once *ξ* exceeds *ξ*_{c}, the state moves quickly from “c” to “d” (i.e., the state of vortex breakdown), through a positive feedback. The slopes of the line segments provide the migration speeds of the resulting shocks. (b) Schematic diagrams for the vertical profile of *F*_{z}. (left) At the onset of SSW. *F*_{z} above *z*_{c} is in state “c” in (a). (right) After SSW. *F*_{z} above *z*_{b} is in state “d” in (a). Because of the gap in the flux, *z*_{b} gradually descends. See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

(a) Schematic diagram for the relationship between FAWA density (*Ae*^{−z/H}; horizontal axis) and E–P flux (*F*_{z}; vertical axis) at a given altitude. The corresponding values of *ξ* are also labeled on the horizontal axis. For *ξ* < *ξ*_{c} = 0.5 (“a” and “b”), a stable steady state is possible. Once *ξ* exceeds *ξ*_{c}, the state moves quickly from “c” to “d” (i.e., the state of vortex breakdown), through a positive feedback. The slopes of the line segments provide the migration speeds of the resulting shocks. (b) Schematic diagrams for the vertical profile of *F*_{z}. (left) At the onset of SSW. *F*_{z} above *z*_{c} is in state “c” in (a). (right) After SSW. *F*_{z} above *z*_{b} is in state “d” in (a). Because of the gap in the flux, *z*_{b} gradually descends. See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

(a) Schematic diagram for the relationship between FAWA density (*Ae*^{−z/H}; horizontal axis) and E–P flux (*F*_{z}; vertical axis) at a given altitude. The corresponding values of *ξ* are also labeled on the horizontal axis. For *ξ* < *ξ*_{c} = 0.5 (“a” and “b”), a stable steady state is possible. Once *ξ* exceeds *ξ*_{c}, the state moves quickly from “c” to “d” (i.e., the state of vortex breakdown), through a positive feedback. The slopes of the line segments provide the migration speeds of the resulting shocks. (b) Schematic diagrams for the vertical profile of *F*_{z}. (left) At the onset of SSW. *F*_{z} above *z*_{c} is in state “c” in (a). (right) After SSW. *F*_{z} above *z*_{b} is in state “d” in (a). Because of the gap in the flux, *z*_{b} gradually descends. See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Note that vortex breakdown starts when

### c. Numerical experiments

*A*and

*z*≤ 110 km and discretized with grid spacing of 100 m. A small numerical diffusion with a diffusion coefficient of 10 m

^{2}s

^{−1}is added to the rhs to smooth out noise. We assume

*C*= 3 × 10

^{−3},

*H*= 7 km, and

*A*as

*α*= 0.4. The initial condition is

*u*

_{REF}increases linearly from 20 m s

^{−1}at

*z*= 10 km to 80 m s

^{−1}at

*z*= 50 km and remains 80 m s

^{−1}above 50 km. This is broadly consistent with the climatology of

*u*

_{REF}at the beginning of the winter (see Fig. 12 below). At the lower boundary (

*z*= 10 km), wave forcing is prescribed as

*z*= 110 km) a radiation condition is assumed.

In the first experiment, we turn on the wave forcing around day 10 (*t*_{on} = 10 days) relatively quickly (*t*_{0} = 6 h) and leave it on. The result is shown in the top row of Fig. 6 in height–time cross sections for *z* ≤ 50 km. After forcing is turned on, a trace of wave activity moves up and quickly amplifies toward the top of the stratosphere (Fig. 6a). During days 12–15, *A* undergoes a rapid increase in the upper stratosphere, whereas *A* and zero

Numerical solutions of Eq. (19) with transient wave forcing applied at the lower boundary. (a),(c),(e) *A*, with contours plotted for 2, 10, 20, … m s^{−1}. (b),(d),(f) ^{−1}. (a),(b) *t*_{0} = 6 h, *t*_{on} = 10 days (“quick spinup”). (c),(d) *t*_{0} = 96 h, *t*_{on} =10 days (“slow spinup”). (e),(f) *t*_{0} = 6 h, *t*_{on} = 10 days, *t*_{off} =15 days (“short pulse”). (g),(h) As in (e) and (f), but for *α* = 0. Contour interval above 200 m s^{−1} in (g) is variable. The thin gray contour in (b), (d), and (f) marks *z* = *z*_{c}. See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Numerical solutions of Eq. (19) with transient wave forcing applied at the lower boundary. (a),(c),(e) *A*, with contours plotted for 2, 10, 20, … m s^{−1}. (b),(d),(f) ^{−1}. (a),(b) *t*_{0} = 6 h, *t*_{on} = 10 days (“quick spinup”). (c),(d) *t*_{0} = 96 h, *t*_{on} =10 days (“slow spinup”). (e),(f) *t*_{0} = 6 h, *t*_{on} = 10 days, *t*_{off} =15 days (“short pulse”). (g),(h) As in (e) and (f), but for *α* = 0. Contour interval above 200 m s^{−1} in (g) is variable. The thin gray contour in (b), (d), and (f) marks *z* = *z*_{c}. See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Numerical solutions of Eq. (19) with transient wave forcing applied at the lower boundary. (a),(c),(e) *A*, with contours plotted for 2, 10, 20, … m s^{−1}. (b),(d),(f) ^{−1}. (a),(b) *t*_{0} = 6 h, *t*_{on} = 10 days (“quick spinup”). (c),(d) *t*_{0} = 96 h, *t*_{on} =10 days (“slow spinup”). (e),(f) *t*_{0} = 6 h, *t*_{on} = 10 days, *t*_{off} =15 days (“short pulse”). (g),(h) As in (e) and (f), but for *α* = 0. Contour interval above 200 m s^{−1} in (g) is variable. The thin gray contour in (b), (d), and (f) marks *z* = *z*_{c}. See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

The progression of vortex breakdown is understood as follows. Soon after the wave forcing is turned on, *F*_{z} reaches a steady value of *F*_{0} = 7.87 × 10^{−2} m^{2} s^{−2} in the lower stratosphere. For this value of *F*_{0}, the threshold (24) is reached at *z*_{c} = 27.5 km. Therefore, a runaway flux convergence and vortex breakdown occur above this altitude first, shifting the vortex state from “c” to “d” in Fig. 5a. Subsequently the region of vortex breakdown (a stagnant reservoir of FAWA like critical layer) expands downward by absorbing the incident wave activity flux from below, which causes the vortex terminus to move down (right panel of Fig. 5b). This is akin to a traffic jam on a highway, where the jam expands backward by absorbing the traffic from behind (Lighthill and Whitham 1955; Richards 1956). Nakamura and Huang (2017, 2018) and Paradise et al. (2019) demonstrate that a similar dynamics is at play in the formation of atmospheric blocking.

*A*and define its altitude as

*z*

_{b}(

*t*), we can estimate its migration speed

*dz*

_{b}(

*t*)/

*dt*using the Rankine–Hugoniot relation [LeVeque 2002, chapter 3]. Just above

*z*

_{b}, from Eq. (27):

*z*

_{b}, from Eq. (23):

*z*

_{b}=

*z*

_{c}where the shock first forms,

*dz*

_{b}/

*dt*= −

*Cu*

_{REF}(

*z*

_{c})/2. This corresponds to the slope of the line segment

*z*

_{c}. We have computed

*z*

_{b}(

*t*) by numerically integrating Eq. (33) and show it in Fig. 6b. It captures the descent of the vortex boundary very well below

*z*

_{c}= 27.5 km.

When the wave forcing is switched on more slowly (*t*_{0} = 96 h, *t*_{on} = 10 days), the result is qualitatively similar but the vortex breakdown occurs earlier and the descent of the vortex terminus is noticeably slower in the upper stratosphere (Figs. 6c,d). In this case, a weak wave activity flux during the early stage of spinup manages to penetrate higher and initiate vortex breakdown from *z*_{c} > 50 km. The vortex terminus is very sharp throughout the stratosphere.

*z*>

*z*

_{b}. In reality, we expect some wave reflection to occur at

*z*=

*z*

_{b}. For example, Killworth and McIntyre (1985) show that a Rossby wave critical layer changes from an absorber to a perfect reflector of waves as it matures. We can crudely model the effects of wave reflection by making the wave forcing shorter, since the reflected (downward) wave packet partially cancels the upward wave activity flux and reduces the persistence of the net wave forcing (Harnik 2009; Kodera et al. 2016). In Figs. 6e and 6f, we apply wave forcing only for 5 days. We accomplish this by adding

*z*= 22 km, below which

*A*and

*z*

_{b}, therefore

*z*

_{b}remains steady. Greater and longer forcings cause the vortex terminus to reach lower altitudes (not shown).

To highlight the role of wave–mean flow interaction in the threshold behavior, in Figs. 6g and 6h we repeat the same experiment with *α* = 0. Without wave–mean flow interaction, the pulse of wave activity propagates upward without modifying

## 4. Effects of damping

*S*is nonconservative sources-sinks of

*A*, whereas

*c*

_{gz}is given by Eq. (18).

^{9}Note that

*u*

_{REF}now varies with time due to nonconservative processes. Note also that Eqs. (35)–(37) ensure

### a. Radiative damping

*τ*(

*z*):

*τ*

^{−1}increases linearly from (30 days)

^{−1}at

*z*= 10 km to (5 days)

^{−1}at

*z*= 50 km and remains (5 days)

^{−1}for

*z*> 50 km. Similar values are used by Holton and Mass (1976) and Plumb and Semeniuk (2003). With (39) and the same diffusion coefficient as before, we solve Eqs. (35)–(37) numerically.

*t*

_{on}= 2 days,

*t*

_{0}= 6 h). The subsequent evolution is qualitatively similar to the undamped case, except the descent of the vortex terminus is much slower: 38 days after the forcing is switched on, it is still 4–5 km above the lower boundary, whereas it reached the boundary in less than 20 days in the adiabatic case (Figs. 6a,b). Also, after the vortex breakdown,

*τ*

^{−1}varies in

*z*much more slowly than exp(−

*z*/

*H*), the approximate solution to Eq. (40) is

*F*

_{1}is independent of height and

*F*

_{1}= 0 for

*z*>

*z*

_{b}since the rhs of Eq. (41) vanishes as

*z*→ ∞. This leads to

*u*

_{REF}with

*u*

_{REF}−

*H*/(

*Cτ*),

*F*

_{0}with

*F*

_{1}, and

*F*

_{z}with

*u*

_{REF}−

*H*/(

*Cτ*) in Eq. (33), that is, radiative damping slows down the descent.

Numerical solutions of Eqs. (35) and (36) with radiative damping. Transient wave forcing is applied at the lower boundary. (a),(c) *A*. (b),(d) *t*_{0} = 6 h, *t*_{on} = 2 days (quick switch-on). (c),(d) *t*_{0} = 6 h, *t*_{on} = 2 days, *t*_{off} = 7 days (short pulse). In (b) and (d), the minimum contour level is 4 m s^{−1} and the minimum value of ^{−1}. See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Numerical solutions of Eqs. (35) and (36) with radiative damping. Transient wave forcing is applied at the lower boundary. (a),(c) *A*. (b),(d) *t*_{0} = 6 h, *t*_{on} = 2 days (quick switch-on). (c),(d) *t*_{0} = 6 h, *t*_{on} = 2 days, *t*_{off} = 7 days (short pulse). In (b) and (d), the minimum contour level is 4 m s^{−1} and the minimum value of ^{−1}. See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Numerical solutions of Eqs. (35) and (36) with radiative damping. Transient wave forcing is applied at the lower boundary. (a),(c) *A*. (b),(d) *t*_{0} = 6 h, *t*_{on} = 2 days (quick switch-on). (c),(d) *t*_{0} = 6 h, *t*_{on} = 2 days, *t*_{off} = 7 days (short pulse). In (b) and (d), the minimum contour level is 4 m s^{−1} and the minimum value of ^{−1}. See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

*t*

_{on}= 2,

*t*

_{off}= 7 days in Eqs. (30) and (34)], not only is the vortex breakdown interrupted halfway, but the vortex recovery proceeds in a discontinuous fashion (Figs. 7c,d). In a stark contrast with Figs. 6c and 6d,

*A*and

*z*=

*z*

_{b}) that moves up at a near constant rate. The upward movement of this “recovery line” is due to a gap in the wave activity flux across it. After the wave forcing is shut off,

*z*

_{b}since

*A*≈ 0. However,

*F*

_{z}is positive for

*z*>

*z*

_{b}, since

*A*≈ [

*u*

_{REF}−

*H*/(

*Cτ*)]/

*α*> 0. Due to this gap in the flux, the recovery line

*z*=

*z*

_{b}must move up, and the rate of the movement is given by

*H*= 7 km and

*τ*≈ 5 days, the ascent rate is 1.47 km day

^{−1}, roughly consistent with Fig. 7c (≈1.4 km day

^{−1}).

### b. Mixing

*A*in Eq. (39) with a flow-dependent damping through

*S*:

*τ*(

*z*) is the radiative damping time introduced earlier. When

*A*= 0),

*τ*

_{mix}recovers the radiative damping time

*τ*. When

*τ*

_{mix}minimizes to 0.2

*τ*. The above formula is ad hoc and meant as a crude parameterization of mixing inhibition by the mean flow (e.g., Ferrari and Nikurashin 2010).

When Eq. (43) is used instead of *S*_{rad} in Eq. (39), the numerical results change substantially. (In this set of experiments, we also triple the diffusion coefficient to 30 m^{2} s^{−1}). Figures 8a and 8d correspond to Figs. 7a and 7d, the only difference being *S* used in Eq. (35). In both the switch-on and the short pulse forcings, FAWA is heavily damped by mixing after it reaches maximum in the upper stratosphere (Figs. 8a,d). In fact, in both cases *A* survives only for a short duration in a very similar way, suggesting an important role of mixing in determining the duration of a wave event. The main difference between the two forcing types is the response of *promotes* vortex recovery when the wave forcing is short. White contours in Fig. 8b identify regions in which

As in Fig. 7, but including mixing parameterization for *A*. (a),(c) *A*. (b),(d) *t*_{0} = 6 h, *t*_{on} = 2 days (quick switch-on). (c),(d) *t*_{0} = 6 h, *t*_{on} = 2 days, *t*_{off} = 7 days (short pulse). White contours in (b) and (d) show

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

As in Fig. 7, but including mixing parameterization for *A*. (a),(c) *A*. (b),(d) *t*_{0} = 6 h, *t*_{on} = 2 days (quick switch-on). (c),(d) *t*_{0} = 6 h, *t*_{on} = 2 days, *t*_{off} = 7 days (short pulse). White contours in (b) and (d) show

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

As in Fig. 7, but including mixing parameterization for *A*. (a),(c) *A*. (b),(d) *t*_{0} = 6 h, *t*_{on} = 2 days (quick switch-on). (c),(d) *t*_{0} = 6 h, *t*_{on} = 2 days, *t*_{off} = 7 days (short pulse). White contours in (b) and (d) show

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Figure 9 shows relationships between *A* and *A* and *u*_{REF} at *z* = 32 km during the two experiments shown in Fig. 8. In both cases *A* and Δ*u* form a linear relation as expected (Figs. 9a,c), whereas *u*_{REF} monotonically decreases as *A* grows and decays. Only in the recovery stage of Fig. 9d, *u*_{REF} is radiatively restored and the hysteresis loop is complete.^{10} The qualitative resemblance between Figs. 9 and 1 (except that Fig. 1d does not form a closed loop partly due to the seasonal demise of the polar vortex) is another evidence that the model is capturing the rudimentary dynamics of the observed SSWs.

(a) *A* (horizontal axis) vs *z* = 32 km during the experiment shown in Figs. 8a and 8b. (b) As in (a), but the vertical axis is *u*_{REF}. (c),(d) As in (a) and (b), but for the experiment shown in Figs. 8c and 8d.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

(a) *A* (horizontal axis) vs *z* = 32 km during the experiment shown in Figs. 8a and 8b. (b) As in (a), but the vertical axis is *u*_{REF}. (c),(d) As in (a) and (b), but for the experiment shown in Figs. 8c and 8d.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

(a) *A* (horizontal axis) vs *z* = 32 km during the experiment shown in Figs. 8a and 8b. (b) As in (a), but the vertical axis is *u*_{REF}. (c),(d) As in (a) and (b), but for the experiment shown in Figs. 8c and 8d.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

## 5. Observational evidence for threshold behavior

*F*

_{z}is a nonlinear (unimodal) function of FAWA due to wave–mean flow interaction and (ii) the mode of

*F*

_{z}defines the threshold for a spontaneous growth of FAWA. The WKB theory based on a vertically propagating stationary Rossby wave predicts

*F*

_{z}to be a quadratic function of

*A*[Eq. (20)], but the foregoing conclusions hold for a more generalized form of Eqs. (19) and (21):

*M*(

*ξ*) is any smooth, nonnegative, unimodal function that does not depend on

*z*, and Γ(

*z*) is defined in Eq. (21). For the WKB theory

*M*(

*ξ*) = (1 −

*ξ*)

*ξ*. In the steady state, the lhs of Eq. (44) vanishes and Eq. (23) is replaced by

*M*(

*ξ*

_{c}) =

*M*

_{max}. The threshold behavior occurs at

*ξ*=

*ξ*

_{c}. If the observed flux fits the above form and if we can demonstrate the modality of

*M*with respect to

*ξ*and exceedance of

*ξ*

_{c}during SSWs, it will be a strong evidence for a threshold behavior, irrespective of the WKB theory or how the wave activity flux is generated. Since [e.g., Andrews et al. (1987), their Eq. (3.5.6); NS10, their Eq. (7)]

*M*at multiple altitudes using reanalysis data as

*ξ*. If the data matches the WKB theory, we should observe

*M*(

*ξ*) = (1 −

*ξ*)

*ξ*at all altitudes.

Figure 10 shows *M*(*ξ*) during the growing stage of SSW (60-day lead plus 2-day lag) for 5, 7, 10, and 20 hPa at 60°N (daily composites of 26 events in ERA-Interim and 25 events in MERRA-2). In Eq. (47), we have assumed *n* = 2 for all events to evaluate *C* [Eq. (17)] but used instantaneous *u*_{REF} since it varies significantly from one season to another. In doing so, we are effectively treating Γ(*z*) as a function of time also. This brings about some scatter in our analysis, but for both reanalysis products, data points for 5-, 7-, and 10-hPa collapse on average on a single curve (dashed red) with a clearly defined peak around *ξ* ≈ 0.7(=*ξ*_{c}). The 20-hPa data trace a similar curve (dashed green) but its peak is lower and shifted to the left of the other curve. On both curves, the last few days around the event are on the right-hand side of the peak (*ξ* > *ξ*_{c}). This is highly suggestive of a threshold behavior, and at least for 5–10 hPa, the height-independent *M*(*ξ*) is consistent with Eq. (44). On the other hand, although the diagnosed *M*_{max} is close to theory (0.25), its location *ξ*_{c} is shifted from the theoretical prediction (0.7 as opposed to 0.5). It may be that the weak dependence of *u*_{REF} on *A* due to mixing toward the event (Figs. 1c and 1d) modifies the form of *M*(*ξ*) from theory. Some discrepancies are undoubtedly related to the limitation of the WKB approximation employed in section 4. For example, *M*(*ξ*) remains positive even after the zonal-mean zonal wind is reversed; this is likely due to a memory effect of rapidly evolving Rossby waves – an aspect that cannot be handled by the quasi-steady assumption in the WKB theory.

(a) *ξ* = *αA*/*u*_{REF} (horizontal axis) vs *M*(*ξ*) (vertical axis) at 60°N during 60 days before and 2 days after SSW events. Composite of 26 events in ERA-Interim between 1979 and 2016. Each dot shows a daily composite on a pressure level. Dots generally progress from low *ξ* to high *ξ* during this period. Error bars indicate ±0.25 standard deviation. Red: 5 hPa. Orange: 7 hPa. Purple: 10 hPa. Green: 20 hPa. Third-order polynomial fits are shown for 5, 7, and 10 hPa (red dashed curve) and for 20 hPa (green dashed curve). Black curve is the theoretical prediction from section 4. (b) As in (a), but based on MERRA-2 (1979–2016, composite of 25 events). See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

(a) *ξ* = *αA*/*u*_{REF} (horizontal axis) vs *M*(*ξ*) (vertical axis) at 60°N during 60 days before and 2 days after SSW events. Composite of 26 events in ERA-Interim between 1979 and 2016. Each dot shows a daily composite on a pressure level. Dots generally progress from low *ξ* to high *ξ* during this period. Error bars indicate ±0.25 standard deviation. Red: 5 hPa. Orange: 7 hPa. Purple: 10 hPa. Green: 20 hPa. Third-order polynomial fits are shown for 5, 7, and 10 hPa (red dashed curve) and for 20 hPa (green dashed curve). Black curve is the theoretical prediction from section 4. (b) As in (a), but based on MERRA-2 (1979–2016, composite of 25 events). See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

(a) *ξ* = *αA*/*u*_{REF} (horizontal axis) vs *M*(*ξ*) (vertical axis) at 60°N during 60 days before and 2 days after SSW events. Composite of 26 events in ERA-Interim between 1979 and 2016. Each dot shows a daily composite on a pressure level. Dots generally progress from low *ξ* to high *ξ* during this period. Error bars indicate ±0.25 standard deviation. Red: 5 hPa. Orange: 7 hPa. Purple: 10 hPa. Green: 20 hPa. Third-order polynomial fits are shown for 5, 7, and 10 hPa (red dashed curve) and for 20 hPa (green dashed curve). Black curve is the theoretical prediction from section 4. (b) As in (a), but based on MERRA-2 (1979–2016, composite of 25 events). See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Since *ξ*_{c} ≈ 0.7, a threshold behavior is expected to occur when *r* at 60°N and *z* = 32 km are compared for 18 winters with SSWs and non-SSW winters using ERA-Interim. While *r* stays above 0.3 for the majority of the non-SSW winters, the winters with SSW experience rapid decrease in *r* after it reaches 0.3, which marks the onset of the event. There is a wide spread after the event, as some events recover quickly (“reflecting” type in Table 2) while others recover more slowly (“absorbing” type), but the threshold behavior during the onset period is shared by all events analyzed. There were a few winters in which *r* dropped below 0.3 temporarily but a full SSW did not follow, presumably because the wave forcing from below was terminated prematurely.

The red line shows a composite time series of *z* = 32 km (10.34 hPa). The same 18 events as in Fig. 1 (listed in Table 1) are used. Lag 0 corresponds to the first reversal of *r* = 0.3. See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

The red line shows a composite time series of *z* = 32 km (10.34 hPa). The same 18 events as in Fig. 1 (listed in Table 1) are used. Lag 0 corresponds to the first reversal of *r* = 0.3. See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

The red line shows a composite time series of *z* = 32 km (10.34 hPa). The same 18 events as in Fig. 1 (listed in Table 1) are used. Lag 0 corresponds to the first reversal of *r* = 0.3. See text for details.

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

Given the ability of *r* to identify an SSW threshold, we propose it as a useful instantaneous diagnostic for the proximity of the polar vortex to breakdown and a tool to understand the intermittency of SSWs. Figure 12 shows height–time cross sections of *u*_{REF} (color shading) and *r* (contours, plotted for 0.4 or less) at 60°N for strong-vortex winters and weak-vortex winters. To differentiate vortex strengths, we have computed the January–March NAM index based on the first EOF of the 50-hPa geopotential height anomaly north of 30°N. We have chosen 15 years in which the index exceeded 0.5 standard deviation as strong-vortex winters, and 11 years in which it exceeded −0.5 standard deviation as weak-vortex winters (Perlwitz and Graf 2001; Baldwin and Dunkerton 1999). (See caption for the list of winters in each category.) Composite for the strong-vortex winters shows a strong and unperturbed *u*_{REF} in the stratosphere through January–February, and *r* drops below threshold only at the end of March (final warming, Fig. 12a). SSWs are characteristically absent, as exemplified by the streak of SSW-free winters during the mid-1990s. Whereas during the weak-vortex winters, both *u*_{REF} and *r* are significantly smaller already at the beginning of January (Fig. 12b). Then *r* drops below 0.3 from late January through February, starting at the top of the stratosphere and subsequently at lower altitudes, as SSWs form around this time of winter. Mixing associated with SSWs sharply decreases *u*_{REF} from the mid- to lower stratosphere, and it remains weak through the rest of the winter. However, persistence of *u*_{REF} in the upper stratosphere increases after SSW, and this significantly delays the timing of final warming and easterly transition compared to the strong-vortex winters (Fig. 12b), consistent with the findings of previous studies (Hu et al. 2014; Lubis et al. 2017).

(a) Composite of height–time cross sections of *u*_{REF} (color shading) and

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

(a) Composite of height–time cross sections of *u*_{REF} (color shading) and

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

(a) Composite of height–time cross sections of *u*_{REF} (color shading) and

Citation: Journal of the Atmospheric Sciences 77, 3; 10.1175/JAS-D-19-0249.1

A weaker *u*_{REF} is more conducive to SSWs because it lowers Γ [Eq. (21)], allowing the threshold to be reached in the mid- to lower stratosphere with less wave forcing [Eq. (45)]. On the other hand, if *u*_{REF} is too weak, it reduces the growth rate, and hence suddenness, of SSWs [Eq. (26)]. This may be the reason why final warming is more pronounced during the strong-vortex winters (Fig. 12a). Since *u*_{REF} depends on mixing and radiative damping associated with preceding wave events, early season troposphere–stratosphere coupling likely plays a key role in the preconditioning of *u*_{REF} for SSWs. Interestingly, during the weak-vortex winters there are pockets of below-threshold *r*s in the lower troposphere (Fig. 12b). Their connection to SSWs is unclear from this figure alone, but it strongly suggests that wave coupling between the troposphere and the stratosphere is more active during the weak-vortex winters.

## 6. Summary and discussion

In the spirit of Jucker (2016), we have sought canonical properties of SSWs that transcend type variations. In particular, we have highlighted the role of wave–mean flow interaction [Eq. (9)] giving rise to suddenness and intermittency of SSWs through a threshold behavior. The threshold arises from a competition between an increasing wave activity *A* and a decreasing zonal-mean zonal wind *F*_{z}. This renders *F*_{z} a quadratic function of *A* [Eq. (20)] and the mode of *F*_{z} defines the threshold [Eq. (24)]. Once *A* exceeds the threshold value, vortex breakdown proceeds spontaneously through a positive feedback between *A* and *F*_{z}, until *F*_{z} between the two, the boundary migrates downward (Fig. 6). Mathematically this process is analogous to a traffic jam on a highway, where the jam expands backward by absorbing the incoming traffic. The resultant 1D model is akin to the traffic jam model of atmospheric blocking considered by Nakamura and Huang (2017, 2018) and Paradise et al. (2019). Further simplified from previous theoretical models (e.g., Holton and Mass 1976; Plumb and Semeniuk 2003), this model permits certain analytical predictions besides the threshold, including the growth rate of SSW [Eq. (26)], migration speeds of the vortex terminus [Eq. (33)] and the recovery line [Eq. (42)]. The model’s response has been tested for a variety of wave forcing and damping and shown to capture the salient features of observed SSWs (Fig. 8) despite the obvious limitations of the WKB theory it is built on.

Examination of reanalysis data lends strong support for the threshold behavior during the onset of SSWs. In the mid- to upper stratosphere, the observed wave activity follows a canonical unimodal function of FAWA and the onset of SSW lies on the right-hand side of its mode (Fig. 10). Consistent with this, the time series of *r* suggests *r* ≈ 0.3 as a “tipping point” for SSW’s onset (Fig. 11). The utility of *r* as an instantaneous diagnostic for the proximity to vortex breakdown has been illustrated (Fig. 12).

Given that the threshold is reached on average 4–5 days prior to an SSW event (Fig. 11), *r* should provide a forecast skill of SSW with a few days of lead time without a need for running a numerical weather prediction model. More importantly, the threshold may be used to better constrain the effects of climate change on the frequency of SSWs. Currently climate models’ confidence in the future projection of SSW frequency is low (e.g., Ayarzagüena et al. 2018). The lack of convergence among the model projections is likely due to model-to-model variations in (i) the threshold itself and (ii) the fashion in which the threshold is reached. On the first point, analyses like Figs. 10 and 11 should reveal model biases. Models with *ξ*_{c} > 0.7 likely underestimate SSW frequencies for the same wave forcing, and vice versa. On the second point, one must ultimately understand how *u*_{REF} and *A* are determined and how they respond to climate perturbations. Particularly important is the early season behavior (preconditioning) of *u*_{REF} as it determines the subsequent wave forcing required to reach the threshold [Eq. (47)]. Sensitivity of the model prediction of *u*_{REF} and its impact on the SSW frequency to the early season troposphere–stratosphere wave coupling is a worthy topic for further investigation.

## Acknowledgments

This work has been supported by NSF Grants AGS1563307 and AGS1909522. The authors thank two anonymous reviewers for constructive critiques on the early version of the manuscript.

## APPENDIX

### Dispersion Relation in the Mercator Coordinate

^{2}

*ϕ*, reads

^{A1}

*ψ*′ is the perturbation geostrophic streamfunction and the following coefficients are considered “slowly varying”:

*k*=

*n*/

*a*(

*n*is the integer zonal wavenumber) is the equatorial wavenumber and (

*Y*,

*Z*,

*T*) ≡ (

*γy*,

*γz*,

*γt*), (

*γ*≪ 1) is the slow coordinate in the meridional plane and time. Substitution in Eq. (A3) yields the following dispersion relation:

*ω*= 0) Eq. (A6) gives

## REFERENCES

Andrews, D., and M. McIntyre, 1976: Planetary waves in horizontal and vertical shear: The generalized Eliassen–Palm relation and the mean zonal acceleration.

,*J. Atmos. Sci.***33**, 2031–2048, https://doi.org/10.1175/1520-0469(1976)033<2031:PWIHAV>2.0.CO;2.Andrews, D., J. Holton, and C. Leovy, 1987:

. Academic Press, 489 pp.*Middle Atmosphere Dynamics*Ayarzagüena, B., and Coauthors, 2018: No robust evidence of future changes in major stratospheric sudden warmings: A multi-model assessment from CCMI.

,*Atmos. Chem. Phys.***18**, 11 277–11 287, https://doi.org/10.5194/acp-18-11277-2018.Baldwin, M. P., and T. J. Dunkerton, 1999: Propagation of the Arctic oscillation from the stratosphere to the troposphere.

,*J. Geophys. Res.***104**, 30 937–30 946, https://doi.org/10.1029/1999JD900445.Baldwin, M. P., and T. J. Dunkerton, 2001: Stratospheric harbingers of anomalous weather regimes.

,*Science***294**, 581–584, https://doi.org/10.1126/science.1063315.Birner, T., and P. D. Williams, 2008: Sudden stratospheric warmings as noise-induced transitions.

,*J. Atmos. Sci.***65**, 3337–3343, https://doi.org/10.1175/2008JAS2770.1.Butler, A. H., D. J. Seidel, S. C. Hardiman, N. Butchart, T. Birner, and A. Match, 2015: Defining sudden stratospheric warmings.

,*Bull. Amer. Meteor. Soc.***96**, 1913–1928, https://doi.org/10.1175/BAMS-D-13-00173.1.Chao, W. C., 1985: Sudden stratospheric warmings as catastrophes.

,*J. Atmos. Sci.***42**, 1631–1646, https://doi.org/10.1175/1520-0469(1985)042<1631:SSWAC>2.0.CO;2.Charlton, A. J., and L. M. Polvani, 2007: A new look at stratospheric sudden warmings. Part I: Climatology and modeling benchmarks.

,*J. Climate***20**, 449–469, https://doi.org/10.1175/JCLI3996.1.Charlton-Perez, A. J., L. M. Polvani, J. Austin, and F. Li, 2008: The frequency and dynamics of stratospheric sudden warmings in the 21st century.

,*J. Geophys. Res.***113**, D16116, https://doi.org/10.1029/2007JD009571.Charney, J., and P. Drazin, 1961: Propagation of planetary-scale disturbances from the lower into the upper atmosphere.

,*J. Geophys. Res.***66**, 83–109, https://doi.org/10.1029/JZ066I001P00083.Charney, J., and J. DeVore, 1979: Multiple flow equilibria in the atmosphere and blocking.

,*J. Atmos. Sci.***36**, 1205–1216, https://doi.org/10.1175/1520-0469(1979)036<1205:MFEITA>2.0.CO;2.Chen, P., and W. A. Robinson, 1992: Propagation of planetary waves between the troposphere and stratosphere.

,*J. Atmos. Sci.***49**, 2533–2545, https://doi.org/10.1175/1520-0469(1992)049<2533:POPWBT>2.0.CO;2.de la Cámara, A., J. R. Albers, T. Birner, R. R. Garcia, P. Hitchcock, D. E. Kinnison, and A. K. Smith, 2017: Sensitivity of sudden stratospheric warmings to previous stratospheric conditions.

,*J. Atmos. Sci.***74**, 2857–2877, https://doi.org/10.1175/JAS-D-17-0136.1.Dee, D. P., and Coauthors, 2011: The Era-Interim reanalysis: Configuration and performance of the data assimilation system.

,*Quart. J. Roy. Meteor. Soc.***137**, 553–597, https://doi.org/10.1002/qj.828.Dickinson, R. E., 1969: Theory of planetary wave-zonal flow interaction.

,*J. Atmos. Sci.***26**, 73–81, https://doi.org/10.1175/1520-0469(1969)026<0073:TOPWZF>2.0.CO;2.Esler, J. G., and R. K. Scott, 2005: Excitation of transient Rossby waves on the stratospheric polar vortex and the barotropic sudden warming.

,*J. Atmos. Sci.***62**, 3661–3682, https://doi.org/10.1175/JAS3557.1.Esler, J. G., L. M. Polvani, and R. K. Scott, 2006: The Antarctic stratospheric sudden warming of 2002: A self-tuned resonance?

,*Geophys. Res. Lett.***33**, L12804, https://doi.org/10.1029/2006GL026034.Ferrari, R., and M. Nikurashin, 2010: Suppression of eddy diffusivity across jets in the Southern Ocean.

,*J. Phys. Oceanogr.***40**, 1501–1519, https://doi.org/10.1175/2010JPO4278.1.Geisler, J. E., 1974: A numerical model of sudden stratospheric warming mechanism.

,*J. Geophys. Res.***79**, 4989–4999, https://doi.org/10.1029/JC079I033P04989.Gelaro, R., and Coauthors, 2017: The Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2).

,*J. Climate***30**, 5419–5454, https://doi.org/10.1175/JCLI-D-16-0758.1.Gerber, E. P., C. Orbe, and L. M. Polvani, 2009: Stratospheric influence on the tropospheric circulation revealed by idealized ensemble forecasts.

,*Geophys. Res. Lett.***36**, L24801, https://doi.org/10.1029/2009GL040913.Harada, Y., A. Goto, H. Hasegawa, N. Fujikawa, H. Naoe, and T. Hirooka, 2010: A major stratospheric sudden warming event in January 2009.

,*J. Atmos. Sci.***67**, 2052–2069, https://doi.org/10.1175/2009JAS3320.1.Harnik, N., 2009: Observed stratospheric downward reflection and its relation to upward pulses of wave activity.

,*J. Geophys. Res.***114**, D08120, https://doi.org/10.1029/2008JD010493.Held, I., 1983: Stationary and quasi-stationary eddies in the extratropical troposphere: Theory.

*Large-Scale Dynamical Processes in the Atmosphere*, B. Hoskins and R. Pearce, Eds., Academic Press, 127–168.Hitchcock, P., and I. R. Simpson, 2014: The downward influence of stratospheric sudden warmings.

,*J. Atmos. Sci.***71**, 3856–3876, https://doi.org/10.1175/JAS-D-14-0012.1.Hitchcock, P., and P. H. Haynes, 2016: Stratospheric control of planetary waves.

,*Geophys. Res. Lett.***43**, 11 884–11 892, https://doi.org/10.1002/2016GL071372.Holton, J. R., 1976: A semi-spectral numerical model for wave-mean flow interactions in the stratosphere: Application to sudden stratospheric warmings.

,*J. Atmos. Sci.***33**, 1639–1649, https://doi.org/10.1175/1520-0469(1976)033<1639:ASSNMF>2.0.CO;2.Holton, J. R., and C. Mass, 1976: Stratospheric vacillation cycles.

,*J. Atmos. Sci.***33**, 2218–2225, https://doi.org/10.1175/1520-0469(1976)033<2218:SVC>2.0.CO;2.Holton, J. R., and T. Dunkerton, 1978: On the role of wave transience and dissipation in stratospheric mean flow vacillations.

,*J. Atmos. Sci.***35**, 740–744, https://doi.org/10.1175/1520-0469(1978)035<0740:OTROWT>2.0.CO;2.Hoskins, B. J., and D. J. Karoly, 1981: The steady linear response of a spherical atmosphere to thermal and orographic forcing.

,*J. Atmos. Sci.***38**, 1179–1196, https://doi.org/10.1175/1520-0469(1981)038<1179:TSLROA>2.0.CO;2.Hu, J., R. Ren, and H. Xu, 2014: Occurrence of winter stratospheric sudden warming events and the seasonal timing of spring stratospheric final warming.

,*J. Atmos. Sci.***71**, 2319–2334, https://doi.org/10.1175/JAS-D-13-0349.1.Jucker, M., 2016: Are sudden stratospheric warmings generic? Insights from an idealized GCM.

,*J. Atmos. Sci.***73**, 5061–5080, https://doi.org/10.1175/JAS-D-15-0353.1.Killworth, P., and M. McIntyre, 1985: Do Rossby-wave critical layers absorb, reflect, or over-reflect?

,*J. Fluid Mech.***161**, 449–492, https://doi.org/10.1017/S0022112085003019.Kodera, K., and M. Chiba, 1995: Tropospheric circulation changes associated with stratospheric sudden warmings—A case study.

,*J. Geophys. Res.***100**, 11 055–11 068, https://doi.org/10.1029/95JD00771.Kodera, K., H. Mukougawa, P. Maury, M. Ueda, and C. Claud, 2016: Absorbing and reflecting sudden stratospheric warming events and their relationship with tropospheric circulation.

,*J. Geophys. Res. Atmos.***121**, 80–94, https://doi.org/10.1002/2015JD023359.Kushner, P. J., and L. M. Polvani, 2004: Stratosphere–troposphere coupling in a relatively simple AGCM: The role of eddies.

,*J. Climate***17**, 629–639, https://doi.org/10.1175/1520-0442(2004)017<0629:SCIARS>2.0.CO;2.Labitzke, K., and B. Naujokat, 2000: The lower Arctic stratosphere in winter since 1952.

,*SPARC Newsl.***15**, 11–14.LeVeque, R., 2002:

. Cambridge University Press, 558 pp.*Finite Volume Methods for Hyperbolic Problems*Lighthill, M. J., and G. B. Whitham, 1955: On kinematic waves II. A theory of traffic flow on long crowded roads.

,*Proc. Roy. Soc. London***229**, 317–345, https://doi.org/10.1098/rspa.1955.0089.Limpasuvan, V., D. W. J. Thompson, and D. L. Hartmann, 2004: The life cycle of the Northern Hemisphere sudden stratospheric warmings.

,*J. Climate***17**, 2584–2596, https://doi.org/10.1175/1520-0442(2004)017<2584:TLCOTN>2.0.CO;2.Liu, Y. S., and R. K. Scott, 2015: The onset of the barotropic sudden warming in a global model.

,*Quart. J. Roy. Meteor. Soc.***141**, 2944–2955, https://doi.org/10.1002/QJ.2580.Lubis, S. W., V. Silverman, K. Matthes, N. Harnik, N.-E. Omrani, and S. Wahl, 2017: How does downward planetary wave coupling affect polar stratospheric ozone in the Arctic winter stratosphere?

,*Atmos. Chem. Phys.***17**, 2437–2458, https://doi.org/10.5194/acp-17-2437-2017.Lubis, S. W., C. S. Y. Huang, and N. Nakamura, 2018a: Role of finite-amplitude eddies and mixing in the life cycle of stratospheric sudden warmings.

,*J. Atmos. Sci.***75**, 3987–4003, https://doi.org/10.1175/JAS-D-18-0138.1.Lubis, S. W., C. S. Y. Huang, N. Nakamura, N. E. Omrani, and M. Jucker, 2018b: Role of finite-amplitude Rossby waves and nonconservative processes in downward migration of extratropical flow anomalies.

,*J. Atmos. Sci.***75**, 1385–1401, https://doi.org/10.1175/JAS-D-17-0376.1.Martineau, P., and S.-W. Son, 2015: Onset of circulation anomalies during stratospheric vortex weakening events: The role of planetary-scale waves.

,*J. Climate***28**, 7347–7370, https://doi.org/10.1175/JCLI-D-14-00478.1.Martius, O., L. M. Polvani, and H. C. Davies, 2009: Blocking precursors to stratospheric sudden warming events.

,*Geophys. Res. Lett.***36**, L14806, https://doi.org/10.1029/2009GL038776.Matsuno, T., 1971: A dynamical model of the stratospheric sudden warming.

,*J. Atmos. Sci.***28**, 1479–1494, https://doi.org/10.1175/1520-0469(1971)028<1479:ADMOTS>2.0.CO;2.Matthewman, N. J., and J. G. Esler, 2011: Stratospheric sudden warmings as self-tuning resonances. Part I: Vortex splitting events.

,*J. Atmos. Sci.***68**, 2481–2504, https://doi.org/10.1175/JAS-D-11-07.1.McIntyre, M., 1982: How well do we understand the dynamics of stratospheric warmings?

,*J. Meteor. Soc. Japan***60**, 37–65, https://doi.org/10.2151/jmsj1965.60.1_37.McIntyre, M., and T. Palmer, 1983: Breaking planetary waves in the stratosphere.

,*Nature***305**, 593–600, https://doi.org/10.1038/305593a0.McIntyre, M., and T. Palmer, 1984: The “surf zone” in the stratosphere.

,*J. Atmos. Terr. Phys.***46**, 825–849, https://doi.org/10.1016/0021-9169(84)90063-1.McLandress, C., and T. G. Shepherd, 2009: Impact of climate change on stratospheric sudden warmings as simulated by the Canadian middle atmosphere model.

,*J. Climate***22**, 5449–5463, https://doi.org/10.1175/2009JCLI3069.1.Methven, J., and P. Berrisford, 2015: The slowly evolving background state of the atmosphere.

,*Quart. J. Roy. Meteor. Soc.***141**, 2237–2258, https://doi.org/10.1002/QJ.2518.Mukougawa, H., and T. Hirooka, 2004: Predictability of stratospheric sudden warming: A case study for 1998/99 winter.

,*Mon. Wea. Rev.***132**, 1764–1776, https://doi.org/10.1175/1520-0493(2004)132<1764:POSSWA>2.0.CO;2.Nakamura, N., and A. Solomon, 2010: Finite-amplitude wave activity and mean flow adjustments in the atmospheric general circulation. Part I: Quasigeostrophic theory and analysis.

,*J. Atmos. Sci.***67**, 3967–3983, https://doi.org/10.1175/2010JAS3503.1.Nakamura, N., and D. Zhu, 2010: Finite-amplitude wave activity and diffusive flux of potential vorticity in eddy–mean flow interaction.

,*J. Atmos. Sci.***67**, 2701–2716, https://doi.org/10.1175/2010JAS3432.1.Nakamura, N., and A. Solomon, 2011: Finite-amplitude wave activity and mean flow adjustments in the atmospheric general circulation. Part II: Analysis in the isentropic coordinate.

,*J. Atmos. Sci.***68**, 2783–2799, https://doi.org/10.1175/2011JAS3685.1.Nakamura, N., and C. S. Y. Huang, 2017: Local wave activity and the onset of blocking along a potential vorticity front.

,*J. Atmos. Sci.***74**, 2341–2362, https://doi.org/10.1175/JAS-D-17-0029.1.Nakamura, N., and C. S. Y. Huang, 2018: Atmospheric blocking as a traffic jam in the jet stream.

,*Science***361**, 42–47, https://doi.org/10.1126/science.aat0721.