## 1. Introduction

In small-amplitude calculations of moist baroclinic instability over a wide range of climates in an idealized GCM, O’Gorman et al. (2018) found that the most unstable mode transitions from a quasi-periodic wave to an isolated vortex at a midlatitude surface air temperature of roughly 292 K.^{1} The structure of the vortex mode that emerged in warm climates (Fig. 1a) consists of a dipole of interlocking potential vorticity (PV) anomalies above the boundary layer: cyclonic in the lower free troposphere and anticyclonic in the upper troposphere. Warm-air advection to the east of the cyclonic anomaly and to the west of the anticyclonic anomaly leads to ascent and diabatic PV generation from latent heat release in the form of a dipole. We note that, because the mode is found for warm-climate simulations, the diabatic generation extends higher in the atmosphere than it would in the current climate.

The constellation of PV anomalies and diabatic PV generation is such that the anomalies are amplified and maintained against the background shear flow. The resulting vortex mode bears the principal characteristic of a diabatic Rossby vortex (DRV) that has been found in a range of different simulations and observation. Its emergence as the fastest-growing mode within the moist baroclinic instability calculations of O’Gorman et al. (2018) points to the profound modifying influence that latent heating has on the structure of fast-growing disturbances in a warming climate.

DRVs first emerged as an alternative mode of instability in idealized studies of moist baroclinic instability. The presence of moisture greatly enriches the dynamics of unstable modes because condensation and precipitation are strongly associated with ascending but not descending motion, and hence an additional nonlinearity is introduced into the thermodynamic equation (O’Gorman 2011). Emanuel et al. (1987) represented condensational heating in Eady and two-level semigeostrophic models by assuming saturated moist-adiabatic ascent. This assumption leads to a nonlinear factor *r*(*w*) that is a function of the vertical velocity *w* and reduces the potential vorticity [or static stability in quasigeostrophic (QG) models] by a factor *r* < 1 in updrafts while leaving it unchanged in downdrafts (*r* = 1). While this parameterization is a simplification, it captures the essential asymmetry that is introduced through the irreversible fallout of condensate during precipitation, and *r* may be calculated from the observed temperature and pressure distributions in the atmosphere. The effect of this condensational heating was to increase the growth rates and decrease the area of ascent of growing modes with respect to dry waves, results that are borne out well by moist baroclinic life cycle studies with shallow water models (Lambaerts et al. 2012) or more comprehensive forecasting models (Booth et al. 2015). When moist instability calculations were done with a more realistic reduction factor *r*(*z*) that varied vertically, the short wavelength cutoff of the Eady model disappeared (Whitaker and Davis 1994; Moore and Montgomery 2004). A new mode of instability emerged at shorter wavelengths that could intensify without the presence of upper-level forcing (Montgomery and Farrell 1991, 1992; Whitaker and Davis 1994; Moore and Montgomery 2004). This mode now grew through the interaction of a surface potential temperature anomaly and an interior PV anomaly, rather than primarily through an interaction of anomalies of potential temperature at the surface and lid, and the budget of eddy available potential energy was dominated by diabatic rather than baroclinic generation.

Subsequent three-dimensional simulations with a mesoscale model by Moore and Montgomery (2005) showed that this alternative diabatic growth mechanism could generate isolated coherent PV-dipole structures consisting of a phase locked low-level cyclonic anomaly and a midtropospheric anticyclonic anomaly starting from an initial moist baroclinic environment without upper-level forcing. The isolated and diabatic character of such a growing disturbance without upper-level forcing led Moore and Montgomery (2004) to classify it as “diabatic Rossby vortex”—a term that we adopt in this paper. More recently, idealized channel simulations of cyclone development using a weather forecasting model (Tierney et al. 2018) showed signs of breakup into “jagged diabatic” PV structures reminiscent of a DRV at sufficiently warm temperatures, in line with the results of O’Gorman et al. (2018). Going beyond initial-value problems, turbulent simulations on a beta plane using moist two-layer QG or shallow water equations showed a transition from a smooth large-scale jet flow to a jet flow disrupted by the presence of small-scale vortices that rapidly intensify through moist-dynamical feedbacks in the strongly precipitating regime of the simulation (Lapeyre and Held 2004; Bembenek et al. 2020). In Bembenek et al. (2020) these vortices were explicitly likened to DRVs. We have also found growing DRVs in the fully nonlinear warm-climate simulations of O’Gorman et al. (2018). An example of such a DRV is shown in Fig. 1c, and it exhibits considerable similarity with the DRV mode calculated by repeated rescaling to small amplitude (Fig. 1a), although the positive PV anomaly extends higher into the atmosphere and the negative PV anomaly and diabatic PV generation are more concentrated at the upper tropopause.

DRVs have also been invoked to account for the initial phase prior to explosive growth of certain cyclones in in operational analyses and realistic simulations. Both the European extreme storm “Lothar” in 1999 (Wernli et al. 2002) and the explosive East Coast winter storm in 2005 (Moore et al. 2008) were shown to propagate and intensify moderately through diabatic effects without upper-level forcing before intensifying explosively through upper-level interactions in a secondary growth process. In this paper, we will remain focused on the initial phase of diabatic self-amplification/propagation without considering upper interactions. While the isolated and vortical structure of the East Coast winter storm led Moore et al. (2008) to classify it as a DRV following the terminology of Moore and Montgomery (2004), the rapid propagational character of Lothar, faster than the ambient winds, led Wernli et al. (2002) to classify it as a “diabatic Rossby wave” (DRW) with the positive diabatic PV generation to the east of the low-level cyclonic PV anomaly playing the role of meridional PV advection in a classic dry Rossby wave as discussed in Parker and Thorpe (1995). Both DRV and DRW refer to the same phenomenon, but neither name is fully satisfactory since such storms are isolated like vortices but propagate through PV generation like a wave (Boettcher and Wernli 2013). The upper-level negative PV anomaly is found to be relatively weak in observed storms, and some uncertainty exists in the literature as to when latent heating leads to growth through interaction of the positive low-level PV anomaly with a self-induced negative upper PV anomaly, or rather just leads to propagation of the low-level PV anomaly. The importance of diabatic effects in individual case studies of rapid cyclogenesis, led Boettcher and Wernli (2013) to study DRVs more systematically by compiling a 10-yr (2001–10) climatology of DRV tracks for the North Pacific and North Atlantic Oceans. DRVs occurred at an average rate of 81 systems per year over the North Pacific and 43 system per year over the North Atlantic. In line with the case studies of Lothar and the East Coast winter storm, DRVs in the current climate were found to propagate with moderate intensification before interacting in a second phase with a preexisting upper-level PV anomaly or jet stream.

It is clear from the literature that DRVs constitute an alternative diabatic growth mechanism that relies both on sufficient baroclinicity and moisture and that produces relatively small-scale modes that can self-amplify exponentially even without the presence of upper-level forcing. Currently, there is no theory for the growth rate and length scale of DRVs or the wave–vortex transition that occurs at higher temperatures in moist baroclinic instability simulations. Analytically tractable models of dry and moist baroclinic instability (Eady 1949; Charney 1947; Phillips 1954; Emanuel et al. 1987; Zurita-Gotor 2005) form much of the basis of our theoretical understanding of cyclones due to their ability to isolate the mechanism of cyclone formation in a conceptually simple model and to relate growth rate and length scale of cyclones to atmospheric parameters in a quantitative way. Given the importance of diabatic effects in cyclogenesis both in the current and future climate, it seems desirable to develop an equivalent conceptually simple model for a DRV.

To this end, we introduce in this paper a minimal moist two-layer QG model with tilted upper and lower boundaries and show that it is capable of producing a DRV mode. Latent heating is represented by an assumption of saturated ascent in updrafts following previous work (Emanuel et al. 1987; Fantini 1995; Zurita-Gotor 2005). We tilt the model boundaries at a slope equal to that of the mean isentropes to make the two-layer model an analog of the interior of the Eady model in which dry-baroclinic instability has been shut off but any moist instability is retained. This allows us to transition to a pure DRV solution within a conceptually simple model. We note that this model is similar in spirit to the unbounded balanced shear flow studied by Snyder and Lindzen (1991) to demonstrate the possibility of growth through diabatically generated interior anomalies in a setup that is dry modally stable. However, Snyder and Lindzen (1991) allowed for negative latent heating in descent regions and so obtained periodic wave solutions rather than an isolated DRV.

We begin in section 2 by formulating the tilted two-layer model and showing that it produces a DRV mode. We then study its PV budget and derive the dispersion relation of the DRV mode analytically, a significant novelty of this paper. Asymptotic solutions for the growth rate and ascent length of the DRV are found in the limit of small *r*. We also solve the dispersion relation for the infinite domain numerically by root finding for the whole range of *r*. In section 3, we study the emergence of DRV modes in the more general case that includes nonzero meridional PV gradients. We first introduce a simple PV argument to explain the wave–vortex transition observed to occur as latent heating becomes dominant in the moist baroclinic instability simulations of O’Gorman et al. (2018). We then generate a phase diagram for when the most unstable mode in a partially tilted two-layer model is a periodic wave versus a DRV as a function of the PV gradients and *r*. In section 4, we compare the warm-climate DRV mode and DRV storm from the idealized GCM simulations of O’Gorman et al. (2018) with two storms in the present climate that have previously been found to have the characteristics of DRVs. Last, in section 5 we summarize our results and discuss their implications.

## 2. A simple model for a DRV

### a. Model formulation

*f*plane:

*ψ*is the streamfunction,

*w*is the vertical velocity,

*N*

^{2}is the constant static stability,

*f*is the Coriolis parameter,

*J*(

*A*,

*B*) =

*A*−

_{x}B_{y}*B*is the Jacobian, and

_{x}A_{y}*r*in regions of ascent. Under an assumption of saturated moist-adiabatic ascent,

*θ*and

*θ** are the potential and saturated equivalent potential temperature, respectively, and Γ

*and Γ*

_{d}*are the dry-adiabatic and moist-adiabatic lapse rates, respectively. The reduction factor in the ascent region varies strongly in the vertical direction. In cyclones with strong diabatic heating,*

_{m}*r*can go all the way to zero in the interior and tend toward 1 as the tropopause is reached. Averaged in the vertical direction,

*r*= 0.1 is a typical value for the current climate and

*r*= 0.01 for the warm-climate GCM simulations in O’Gorman et al. (2018). In physical terms, the nonlinear factor

*r*(

*w*) represents the fact that while moist ascending air releases latent heat upon condensation and feels a locally reduced static stability, the descending air is subsaturated (after irreversible fallout of condensate by precipitation) and thus feels the full static stability. Moist thermodynamics thus introduces an additional nonlinearity into the equations that greatly enriches the dynamics.

We simplify the dynamics further by discretizing the equations in the vertical direction into two equal layers of height (Fig. 2), anticipating that the two layers will be sufficient to represent the PV-dipole structure of the DRV. We introduce a barotropic streamfunction *ϕ* = (*ψ*_{1} + *ψ*_{2})/2 and a baroclinic streamfunction *τ* = (*ψ*_{1} − *ψ*_{2})/2, where 1 refers to the upper layer and 2 refers to the lower layer. The layer interface height is *η* = −(*f*/*g*′)(*ψ*_{1} − *ψ*_{2}), with *g*′ = *g*(*θ*_{1} − *θ*_{2})/*θ*_{0} where *g* is the gravitational constant, *θ*_{1} and *θ*_{2} are potential temperatures in each layer, and *θ*_{0} is a reference potential temperature.

Schematic of the (a) untilted and (b) tilted two-layer model with basic-state interface height *η*, boundary slopes *h*_{1}* _{y}* =

*h*

_{2}

*=*

_{y}*η*, and basic-state PV gradients

_{y}*q*

_{1}

*and*

_{y}*q*

_{2}

*. Also shown is the basic-state zonal mean wind profile (the dot represents flow out of the page, and the cross represents flow into the page), which is the same for the tilted and untilted models.*

_{y}Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Schematic of the (a) untilted and (b) tilted two-layer model with basic-state interface height *η*, boundary slopes *h*_{1}* _{y}* =

*h*

_{2}

*=*

_{y}*η*, and basic-state PV gradients

_{y}*q*

_{1}

*and*

_{y}*q*

_{2}

*. Also shown is the basic-state zonal mean wind profile (the dot represents flow out of the page, and the cross represents flow into the page), which is the same for the tilted and untilted models.*

_{y}Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Schematic of the (a) untilted and (b) tilted two-layer model with basic-state interface height *η*, boundary slopes *h*_{1}* _{y}* =

*h*

_{2}

*=*

_{y}*η*, and basic-state PV gradients

_{y}*q*

_{1}

*and*

_{y}*q*

_{2}

*. Also shown is the basic-state zonal mean wind profile (the dot represents flow out of the page, and the cross represents flow into the page), which is the same for the tilted and untilted models.*

_{y}Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

We assume small perturbations about a basic state *τ*_{0} = −*Uy* corresponding to a flow *u*_{1} = −*ψ*_{1}* _{y}* =

*U*in the upper layer and

*u*

_{2}= −

*ψ*

_{2}

*= −*

_{y}*U*in the lower layer. The small amplitude of the perturbations allows us to linearize the advection terms, but the thermodynamic equation remains nonlinear because of the latent heating term.

*h*

_{1}(

*y*) and

*h*

_{2}(

*y*), respectively, to have slopes in the meridional direction of

*h*

_{1}

*=*

_{y}*h*

_{2}

*=*

_{y}*η*= −(2

_{y}*f*/

*g*′)

*τ*

_{0}

*so as to match the slope of the basic-state layer interface*

_{y}*η*(Fig. 2b) in contrast to the standard untilted two-layer model (Fig. 2a). This makes our two-layer model an analog of the interior of the Eady model with zero meridional PV gradients

*q*

_{1}

*=*

_{y}*q*

_{2}

*= 0 in the basic state. The dry modal instability through interlocking Rossby waves is thus shut off, but any instability solely due to the moist processes is retained. The equations for the perturbations about the basic state are derived in appendix A in the limit of small-amplitude perturbations and are given here in nondimensional form:*

_{y}*τ*,

*ϕ*, and

*w*are now perturbations about the basic state.

^{2}We note that a dry two-layer model with sloping upper and lower boundaries was already introduced by Bretherton (1966) to allow independent variation of PV gradients and vertical shear in an investigation of the short-wavelength cutoff of baroclinic instability. Evidently, the reason for introducing the tilted boundaries here is different, as we are interested in eliminating the basic-state PV gradients all together.

### b. Numerical simulation

We first solve the tilted-model equations numerically to isolate and study the fastest-growing mode for a given static-stability reduction factor 0 ≤ *r* ≤ 1. To this end, we discretize the equations using second-order central finite differences in a periodic domain in *x*.

We integrate the barotropic and baroclinic vorticity equations [Eqs. (4) and (5)] forward in time for the variables **Φ** = *ϕ _{xx}* and

**T**=

*τ*. Time stepping is performed with MATLAB’S “ode45” function, which is based on an explicit Runge–Kutta (4, 5) formula with an adaptive time step.

_{xx}*w*at each time step from the nonlinear omega equation,

*r*(

*w*) and requires an iterative approach to finding the solution. We solve it iteratively at each time step as [

*r*(

*w*)

^{n}*w*

^{n}^{+1}]

*−*

_{xx}*w*

^{n}^{+1}= rhs, where

*n*is the iteration step. We start the iteration from a random guess for

*w*to define the initial

*r*(

*w*), and we iterate until the root-mean-square (rms) of (

*w*

^{n}^{+1}−

*w*) is smaller than 10

^{n}^{−12}.

We start the time stepping from random initial conditions for **Φ** and **T**. At each time step, we invert **Φ** = *ϕ _{xx}* and

**T**=

*τ*to obtain

_{xx}*ϕ*and

*τ*by imposing that

*ϕ*and

*τ*have zero mean. We then solve for

*w*using the iterative approach to the omega equation described above, and we then update

**Φ**and

**T**using Eqs. (4) and (5). We rescale the amplitudes of the vectors

**Φ**and

**T**by a factor of 100 each time rms (

**x**) > 10, where

**x**= (

**T**,

**Φ**), to avoid large numbers that could cause problems with the numerical representation. We integrate until the nondimensional time is

*t*= 200 when we find that the solution has converged to a normal mode.

The vertical velocity (at time *t* = 200) for *r* = 0.01 is shown in Fig. 3a where we have used a grid spacing of Δ*x* = 0.025 and a domain size of 8*π*. Remarkably, the solution evolves into a DRV with a single spatially localized peak in vertical velocity just like in the warm limit of the idealized GCM calculations of O’Gorman et al. (2018, their Figs. 1f and 2f). The isolated solution is in stark contrast to the spatially periodic structure of moist baroclinic waves. The solution is exponentially growing and fixed in space because the basic-state zonal wind is equal and opposite in each layer, but the DRV would propagate zonally with a more realistic vertical wind profile. We have repeated the calculations using a linear drag on the relative vorticity in the lower layer with a damping time scale of either 10 days (weak drag) or 2.5 days (strong drag). Isolated DRV solutions persist even with drag included, with similar length scale but reduced growth rate relative to the default case with no drag (not shown).

(a) Vertical velocity, (b) PV anomalies (solid) and diabatic PV generation rate (dashed), and (c) total meridional winds (solid) and meridional winds induced by PV anomalies in the same layer (dashed). All quantities are plotted vs *x* for the fastest-growing mode of the tilted-model equations at *r* = 0.01, which is a DRV. In (b) and (c), quantities in the upper layer are shown in blue and quantities in the lower layer are shown in red. The domain size is *L* = 8*π*, and the grid spacing is Δ*x* = 0.025. All quantities are nondimensional, and the overall magnitude of the DRV is arbitrary. The results in (b) and (c) have been zoomed in around the location of ascending motion to better show the structure of the fields since the DRV occupies only a small fraction of the domain.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

(a) Vertical velocity, (b) PV anomalies (solid) and diabatic PV generation rate (dashed), and (c) total meridional winds (solid) and meridional winds induced by PV anomalies in the same layer (dashed). All quantities are plotted vs *x* for the fastest-growing mode of the tilted-model equations at *r* = 0.01, which is a DRV. In (b) and (c), quantities in the upper layer are shown in blue and quantities in the lower layer are shown in red. The domain size is *L* = 8*π*, and the grid spacing is Δ*x* = 0.025. All quantities are nondimensional, and the overall magnitude of the DRV is arbitrary. The results in (b) and (c) have been zoomed in around the location of ascending motion to better show the structure of the fields since the DRV occupies only a small fraction of the domain.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

(a) Vertical velocity, (b) PV anomalies (solid) and diabatic PV generation rate (dashed), and (c) total meridional winds (solid) and meridional winds induced by PV anomalies in the same layer (dashed). All quantities are plotted vs *x* for the fastest-growing mode of the tilted-model equations at *r* = 0.01, which is a DRV. In (b) and (c), quantities in the upper layer are shown in blue and quantities in the lower layer are shown in red. The domain size is *L* = 8*π*, and the grid spacing is Δ*x* = 0.025. All quantities are nondimensional, and the overall magnitude of the DRV is arbitrary. The results in (b) and (c) have been zoomed in around the location of ascending motion to better show the structure of the fields since the DRV occupies only a small fraction of the domain.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

*r*(

*w*)

*w*term as

*r*(

*w*)

*w*=

*w*− [1 −

*r*(

*w*)]

*w*in the thermodynamic equation [Eq. (6)], eliminating the

*w*term using the baroclinic vorticity equation [Eq. (5)], and adding the barotropic vorticity equation [Eq. (4)] to give

*q*

_{2}=

*ϕ*−

_{xx}*τ*+

_{xx}*τ*is the lower-layer PV,

*σ*is the growth rate,

*q*

_{2}

*is zonal advection, and*

_{x}*q*

_{1}=

*ϕ*+

_{xx}*τ*−

_{xx}*τ*, which has diabatic generation given by

*r*shows that the basic PV structure remains similar for all values of 0 ≤

*r*< 1, although the growth rate and horizontal length scale of the ascent region do change when

*r*is varied. At

*r*= 1 the system is dry and stable because there are no contributions from latent heating, and by construction there are no meridional PV gradients to otherwise support baroclinic instability.

We calculate the growth rate of the mode by assuming exponential growth of the rms of **x** = [**Φ**, **T**] over each time step Δ*t* to give *t* can vary. We then average *σ* over the end period of the calculation (*t* = 195–200).

From Fig. 4 we see that within the narrow region of ascent the growth of the positive PV anomaly is due to diabatic PV generation through latent heating that is partially offset by zonal advection. In the region of descending motion to the west, the PV generation due to latent heating is zero and the growth of the PV anomaly is due to zonal advection over a more extended spatial scale. The PV tendency from radiative cooling is spatially constant with a value of

Terms in the lower-layer PV budget vs *x* for the fastest-growing mode of the tilted-model equations at *r* = 0.01, which is a DRV. The domain size is *L* = 8*π*, and the grid spacing is Δ*x* = 0.025. The terms in the PV budget that are shown are the total tendency (blue) and the contributions from zonal advection (red dashed) and latent heating (black dashed–dotted). The PV tendency from radiative cooling is a small constant with a value of

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Terms in the lower-layer PV budget vs *x* for the fastest-growing mode of the tilted-model equations at *r* = 0.01, which is a DRV. The domain size is *L* = 8*π*, and the grid spacing is Δ*x* = 0.025. The terms in the PV budget that are shown are the total tendency (blue) and the contributions from zonal advection (red dashed) and latent heating (black dashed–dotted). The PV tendency from radiative cooling is a small constant with a value of

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Terms in the lower-layer PV budget vs *x* for the fastest-growing mode of the tilted-model equations at *r* = 0.01, which is a DRV. The domain size is *L* = 8*π*, and the grid spacing is Δ*x* = 0.025. The terms in the PV budget that are shown are the total tendency (blue) and the contributions from zonal advection (red dashed) and latent heating (black dashed–dotted). The PV tendency from radiative cooling is a small constant with a value of

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Now that we have isolated the DRV solution within a simplified model, it is possible to develop analytical solutions for its characteristics.

### c. Analytic theory

*w*:

*w*derived for moist baroclinic modes in an untilted two-layer model [cf. with Eq. (12) in Zurita-Gotor 2005] except for the two extra terms −2(

*rw*)

*and (*

_{xx}*r*− 1)

*w*on the left-hand side of the equation and the radiative cooling term

We look for symmetric solutions about the peak in *w* since that is what was obtained in the numerical solutions and since the equation for *w* is symmetric under *x* → −*x*. We put the peak in *w* at *x* = 0, and by symmetry we need only consider the half of the domain *x* ≥ 0, where *w* is ascending between 0 ≤ *x* ≤ *b* and descending for *x* > *b* (see Fig. 5). Here *b* is the location of the boundary between ascent and descent that must be found as part of the solution.

Schematic of the solution of the *w* equation [see Eq. (12)] for the DRV mode. For 0 ≤ *x* ≤ *b* we have ascending motion *w* = *w*_{↑} > 0 and *r* < 1, and for *x* > *b* we have descending motion *w* = *w*_{↓} < 0 and *r* = 1. Here, *x* = *b* is the location of the boundary between ascent and descent that must be found as part of the solution.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Schematic of the solution of the *w* equation [see Eq. (12)] for the DRV mode. For 0 ≤ *x* ≤ *b* we have ascending motion *w* = *w*_{↑} > 0 and *r* < 1, and for *x* > *b* we have descending motion *w* = *w*_{↓} < 0 and *r* = 1. Here, *x* = *b* is the location of the boundary between ascent and descent that must be found as part of the solution.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Schematic of the solution of the *w* equation [see Eq. (12)] for the DRV mode. For 0 ≤ *x* ≤ *b* we have ascending motion *w* = *w*_{↑} > 0 and *r* < 1, and for *x* > *b* we have descending motion *w* = *w*_{↓} < 0 and *r* = 1. Here, *x* = *b* is the location of the boundary between ascent and descent that must be found as part of the solution.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

*x*→ ∞ assuming growing solutions

*σ*≥ 0. In the ascending region, the solution is symmetric about

*x*= 0 and is given by

*r*and

*σ*.

*L*and then take the limit of an infinite domain

*L*→ ∞. Mass conservation expressed as

*a*→ 0 as

*L*→ ∞ if we want solutions for which

*b*and

*w*

_{↑}remain bounded. We need additional constraints to determine the constants

*c*

_{1},

*c*

_{2},

*d*

_{1},

*d*

_{2},

*σ*, and

*b*. We impose that

*w*

_{↑}=

*w*

_{↓}= 0 at

*x*=

*b*, continuity of (

*rw*)

*and (*

_{x}*rw*)

*at*

_{xx}*x*=

*b*and mass conservation. The continuity conditions follow from the continuity of

*ϕ*,

*τ*and the existence of the derivatives in the governing Eqs. (9)–(11). The constraint that

*w*

_{↓}= 0 at

*x*=

*b*gives

*w*solution is arbitrary, which allows us to fix one of the amplitudes without loss of generality. We choose

*d*

_{2}= 1 when

*σ*> 1 and

*d*

_{2}= −1 when

*σ*< 1 to ensure that

*w*

_{↓}< 0. In the limit of

*L*→ ∞ and

*a*→ 0, the resulting equations are

*w*

_{↑}= 0 at

*x*=

*b*, continuity of [

*r*(

*w*)

*w*]

*, continuity of [*

_{x}*r*(

*w*)

*w*]

*, and mass conservation, respectively. The limit of*

_{xx}*L*→ ∞ and

*a*→ 0 in the mass conservation equation must be taken carefully—a subtle point that is discussed in appendix C.

*c*

_{1}and

*c*

_{2}, which also gets rid of the arbitrary constant

*d*

_{2}(see appendix D), we obtain two equations:

*k*

_{1}and

*k*

_{2}[Eqs. (15) and (16)] yield the dispersion relationship for the growth rate

*σ*and half-ascent length

*b*as a function of the static-stability reduction factor

*r*, a key novel result of this paper.

*r*→ 0 it is possible to show analytically that at leading order the growth rate is

^{3}For comparison, the growth rate

*σ*and half-ascent length (equal to

*b*for the DRV) for the fastest-growing modes in an untilted dry two-layer model (Phillips 1954) and moist two-layer model (Emanuel et al. 1987) are given in Table 1. Note that the half-ascent length is just one-quarter of the wavelength for the dry mode and that the results in Emanuel et al. (1987) need to be rescaled by a factor

*r*limit, consistent with the fact that the DRV emerges as the fastest-growing solution in the warm-climate simulations of moist baroclinic instability in O’Gorman et al. (2018). The ascent length of the DRV and moist wave are the same in the small

*r*limit and are smaller than that of a dry wave by about a factor of 15 for

*r*= 0.01.

Comparison of the nondimensional and dimensional growth rate *σ* and half-ascent length in two-layer models for the most unstable dry wave (Phillips 1954), moist wave (Emanuel et al. 1987), and DRV. The half-ascent length is given by one-quarter of the wavelength for the dry wave, and by one-half of the length of the region of ascent for the DRV and moist wave (for the DRV it is *b*). Dimensional values for the growth rate [*σ*(*U*/*L _{D}*)] and ascent length (

*bL*) are calculated using typical scales

_{D}*L*=

_{D}*NH*/(

*f*) = 1000 km/

*U*= 10 m s

^{−1}. The factor of

*L*follows from our choice of nondimensionalization. Growth rates for the moist wave and DRV are presented in the limit of a convectively neutral stratification (

_{D}*r*→ 0). Half-ascent lengths for the moist wave and DRV are presented as the small-

*r*asymptotic expressions for the nondimensional results and are evaluated at

*r*= 0.01 (representative of a warm climate) for the dimensional results, since the ascent length would be zero for

*r*→ 0.

To obtain *σ* and *b* for the full range of *r*, we solve the dispersion equations numerically. Equations (22) and (23) are solved using MATLAB’S fsolve. We start by solving at *r* = 10^{−3} with initial guess *σ* = 1.53 and *b* = 0.06 for the first two values of *r*, and we use linear extrapolation for the initial guesses at each subsequent value of *r*. The results are compared with the time-marching solutions of Eqs. (4), (5), and (7) for a finite domain with periodic boundary conditions in Fig. 6, where we use a larger domain *L* = 32*π* with Δ*x* = 0.084 for all values of *r* to resolve the large and weakly growing solutions as *r* → 1. Note that the time-marching solution at *r* = 1 is not growing and is not shown in Fig. 6. The growth rates from the time-marching and dispersion-relation approaches are in good agreement for values below a critical value of *r* = 0.38, whereas the ascent lengths are in good agreement only for values of *r* less than roughly 0.2. A sample vertical-velocity profile at *r* = 0.01 (Fig. 7) confirms that the *w* profiles from the time marching and dispersion relations are in very good agreement at small *r*. As the critical value of *r* = 0.38 is approached, the root-finding solution for *b* tends to large numbers. For values of *r* > 0.38, we only find solutions to the dispersion relation for which the half-ascent length *b* < 0. These solutions are unphysical and can be discarded.

(a) Growth rate *σ* and (b) half-ascent length *b* from the time-marching solutions of the tilted-model equations [Eqs. (4), (5), and (7)] in a finite periodic domain *L* = 32*π* with Δ*x* = 0.084 (solid red) and from the root-finding of the dispersion Eqs. (22) and (23) for an infinite domain (dashed blue).

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

(a) Growth rate *σ* and (b) half-ascent length *b* from the time-marching solutions of the tilted-model equations [Eqs. (4), (5), and (7)] in a finite periodic domain *L* = 32*π* with Δ*x* = 0.084 (solid red) and from the root-finding of the dispersion Eqs. (22) and (23) for an infinite domain (dashed blue).

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

(a) Growth rate *σ* and (b) half-ascent length *b* from the time-marching solutions of the tilted-model equations [Eqs. (4), (5), and (7)] in a finite periodic domain *L* = 32*π* with Δ*x* = 0.084 (solid red) and from the root-finding of the dispersion Eqs. (22) and (23) for an infinite domain (dashed blue).

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Comparison of the vertical velocity vs *x* for *r* = 0.01 from the time-marching solution [Eqs. (4), (5), and (7)] in a finite periodic domain (solid red) and from the root-finding of the dispersion relation [Eqs. (22) and (23)] for an infinite domain (dashed blue).

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Comparison of the vertical velocity vs *x* for *r* = 0.01 from the time-marching solution [Eqs. (4), (5), and (7)] in a finite periodic domain (solid red) and from the root-finding of the dispersion relation [Eqs. (22) and (23)] for an infinite domain (dashed blue).

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Comparison of the vertical velocity vs *x* for *r* = 0.01 from the time-marching solution [Eqs. (4), (5), and (7)] in a finite periodic domain (solid red) and from the root-finding of the dispersion relation [Eqs. (22) and (23)] for an infinite domain (dashed blue).

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

In mathematical terms, the breakdown of the solution on an infinite domain at *r* = 0.38 can be traced to the point at which *σ*^{2} + *r* – 1 = 0 and the right-hand sides of Eqs. (22) and (23) diverge to infinity. Empirically we find that as this point is approached both *k*_{1} and *k*_{2} also go to zero, such that *b* → ∞ is needed to balance a diverging right-hand side. Setting *k*_{1} = *k*_{2} = 0 we obtain *r* = (3 − *σ*^{2} + *r* − 1 < 0, which implies from Eqs. (15) and (16) that *k*_{2} becomes imaginary while *k*_{1} remains real. Hence, “tan” in Eq. (23) switches to “tanh” while the right-hand side of Eq. (23) becomes negative. To satisfy the equation for growing modes, this requires *b* < 0, which is unphysical.

Thus, while DRV solutions continue to exist for *r* > 0.38 on a finite domain with periodic boundary conditions, isolated DRV solutions on an infinite domain cease to exist.

## 3. The role of meridional PV gradients and the wave–vortex transition

So far we have discussed the emergence of DRVs in a moist two-layer model with zero meridional PV gradients in which PV is generated purely from diabatic effects. While such a setup is a useful idealization for a DRV mode, we are interested in studying how the occurrence of DRVs generalizes to a more realistic situation with PV gradients. We start by considering a qualitative PV argument for how the dynamics changes as diabatic effects become dominant over meridional PV advection, and we then extend our tilted two-layer model to include meridional PV gradients.

### a. PV dynamical perspective on the wave–vortex transition

The transition from a periodic wave to an isolated vortex mode in simulations of moist baroclinic instability in which diabatic effects become dominant can be interpreted elegantly within the PV framework of a moist two-layer model (see Fig. 8). Focusing on the lower layer (the upper layer is analogous), we start from the initial condition of a PV wave train consisting of positive and negative anomalies (Fig. 8a).

A PV perspective on the transition from a periodic wave to an isolated vortex in simulations of moist baroclinic instability. Cyclonic PV anomalies are shown as red plus signs, and anticylonic PV anomalies are shown as blue minus signs. Blue arrows illustrate the direction of horizontal winds induced by the anomalies. Starting from an initial condition consisting of (a) a wave train of PV anomalies in the lower layer, (b) only positive anomalies survive because of the asymmetry in diabatic generation. However, a series of positive anomalies have weaker meridional flow between them that leads to weaker diabatic PV generation as compared with an isolated anomaly, and thus (c) one anomaly is a faster-growing mode of the system.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

A PV perspective on the transition from a periodic wave to an isolated vortex in simulations of moist baroclinic instability. Cyclonic PV anomalies are shown as red plus signs, and anticylonic PV anomalies are shown as blue minus signs. Blue arrows illustrate the direction of horizontal winds induced by the anomalies. Starting from an initial condition consisting of (a) a wave train of PV anomalies in the lower layer, (b) only positive anomalies survive because of the asymmetry in diabatic generation. However, a series of positive anomalies have weaker meridional flow between them that leads to weaker diabatic PV generation as compared with an isolated anomaly, and thus (c) one anomaly is a faster-growing mode of the system.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

A PV perspective on the transition from a periodic wave to an isolated vortex in simulations of moist baroclinic instability. Cyclonic PV anomalies are shown as red plus signs, and anticylonic PV anomalies are shown as blue minus signs. Blue arrows illustrate the direction of horizontal winds induced by the anomalies. Starting from an initial condition consisting of (a) a wave train of PV anomalies in the lower layer, (b) only positive anomalies survive because of the asymmetry in diabatic generation. However, a series of positive anomalies have weaker meridional flow between them that leads to weaker diabatic PV generation as compared with an isolated anomaly, and thus (c) one anomaly is a faster-growing mode of the system.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

In a regime in which PV gradients are dominant, meridional advection of the negative background PV gradient in the lower layer would generate a positive PV tendency to the east of a positive PV anomaly, and a negative PV tendency to its west. These advective PV tendencies give rise to a eastward propagating Rossby wave that can phase lock and grow by interacting with a counterpropagating Rossby wave in the upper layer.

If instead we are in a regime in which the diabatic generation of PV is dominant, ascent to the east of a positive PV anomaly causes latent heat release that generates a positive PV tendency in the lower layer, whereas descent to its west does not generate latent heating and thus there is no negative PV tendency in the lower layer. Thus, only positive PV anomalies survive in the lower layer (with repeated amplitude rescaling to calculate the growing mode) due to the asymmetry in diabatic PV generation (Fig. 8b).

However, a series of positive PV anomalies have weaker meridional flow between them (because of cancellation of the induced flow from neighboring PV anomalies; Fig. 8b) as compared with the meridional flow surrounding an isolated single PV anomaly. Thus, one PV anomaly is a faster-growing mode of the system when diabatic PV generation is dominant over meridional PV advection (Fig. 8c).

### b. Including PV gradients in the two-layer tilted model

*h*

_{1}

*=*

_{y}*α*

_{1}and

*h*

_{2}

*=*

_{y}*α*

_{2}. We retain the basic state

*τ*

_{0}= −

*y*(nondimensionalized) and

*ϕ*

_{0}= 0 corresponding to a shear flow in thermal wind balance. The perturbation equations around this basic state are derived in appendix A and take the form

*q*

_{1}

*= 1 −*

_{y}*α*

_{1}and

*q*

_{2}

*= −1 +*

_{y}*α*

_{2}.

We solve the equations numerically for given values of *r*, *α*_{1}, and *α*_{2} using time marching as in section 2b. We use a domain size of *L* = 8*π*, but to calculate solutions for a large parameter space we use a larger grid spacing (Δ*x* = 0.13) than before. We classify the most unstable solution as either an isolated DRV (single local maximum in *w*) or a periodic wave (multiple local maxima in *w*) or stable. Local maxima in *w* for which *w* < 0 are not counted. Cases that are stable or very weakly growing (*σ* < 0.09) are counted as stable. The results of this classification are shown in Fig. 9 along with the growth rate *σ* in Figs. 9a–c and half-ascent length *b* in Figs. 9d–f.

Phase diagram showing whether the most unstable mode of the tilted two-layer model is a periodic wave (blue) or a DRV (red) or stable (white) for a setup where (a),(d) the PV gradients are equal and opposite in strength *q*_{1}* _{y}* = −

*q*

_{2}

*for multiple values of*

_{y}*r*and for a setup where the PV gradients vary independently from each other and (b),(e)

*r*= 0.1 or (c),(f)

*r*= 0.01. Dashed lines in (a)–(c) show the growth rate, and dashed lines in (d)–(f) show the half-ascent length

*b*. The domain length is

*L*= 8

*π*, and the grid spacing is Δ

*x*= 0.13. Note that the standard two-layer configuration has positive upper-layer PV gradient (

*q*

_{1}

*> 0) and negative lower-layer PV gradient (*

_{y}*q*

_{2}

*< 0), which corresponds to the upper-left quadrants in (b), (c), (e), and (f).*

_{y}Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Phase diagram showing whether the most unstable mode of the tilted two-layer model is a periodic wave (blue) or a DRV (red) or stable (white) for a setup where (a),(d) the PV gradients are equal and opposite in strength *q*_{1}* _{y}* = −

*q*

_{2}

*for multiple values of*

_{y}*r*and for a setup where the PV gradients vary independently from each other and (b),(e)

*r*= 0.1 or (c),(f)

*r*= 0.01. Dashed lines in (a)–(c) show the growth rate, and dashed lines in (d)–(f) show the half-ascent length

*b*. The domain length is

*L*= 8

*π*, and the grid spacing is Δ

*x*= 0.13. Note that the standard two-layer configuration has positive upper-layer PV gradient (

*q*

_{1}

*> 0) and negative lower-layer PV gradient (*

_{y}*q*

_{2}

*< 0), which corresponds to the upper-left quadrants in (b), (c), (e), and (f).*

_{y}Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Phase diagram showing whether the most unstable mode of the tilted two-layer model is a periodic wave (blue) or a DRV (red) or stable (white) for a setup where (a),(d) the PV gradients are equal and opposite in strength *q*_{1}* _{y}* = −

*q*

_{2}

*for multiple values of*

_{y}*r*and for a setup where the PV gradients vary independently from each other and (b),(e)

*r*= 0.1 or (c),(f)

*r*= 0.01. Dashed lines in (a)–(c) show the growth rate, and dashed lines in (d)–(f) show the half-ascent length

*b*. The domain length is

*L*= 8

*π*, and the grid spacing is Δ

*x*= 0.13. Note that the standard two-layer configuration has positive upper-layer PV gradient (

*q*

_{1}

*> 0) and negative lower-layer PV gradient (*

_{y}*q*

_{2}

*< 0), which corresponds to the upper-left quadrants in (b), (c), (e), and (f).*

_{y}Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

We begin by focusing on the equally tilted case *α*_{1} = *α*_{2} = *α* for which *q*_{1}* _{y}* = −

*q*

_{2}

*= 1 −*

_{y}*α*. We let 0 ≤

*α*≤ 2, which includes the classic untilted regime with

*q*

_{1}

*= 1 and*

_{y}*q*

_{2}

*= −1 at*

_{y}*α*= 0, the no-PV-gradient regime

*q*

_{1}

*=*

_{y}*q*

_{2}

*= 0 from the previous section at*

_{y}*α*= 1, and a reversed-PV-gradient regime

*q*

_{1}

*= −1 and*

_{y}*q*

_{2}

*= 1 at*

_{y}*α*= 2. For

*α*< 1, a DRV emerges as the fastest-growing solution when the magnitudes of the upper and lower PV gradients are weaker than a threshhold

*q*

_{y}_{crit}of roughly 0.7 such that

*q*

_{1}

*= −*

_{y}*q*

_{2}

*<*

_{y}*q*

_{y}_{crit}as shown in Fig. 9a.

^{4}When instead the magnitudes of the upper and lower PV gradients are greater than the critical threshold (

*q*

_{1}

*= −*

_{y}*q*

_{2}

*>*

_{y}*q*

_{y}_{crit}), the periodic wave emerges as the fastest-growing solution for all values of

*r*. For the cases with

*α*> 1 in which the PV gradients are reversed from their usual directions, a DRV is the fastest-growing solution, but it is necessary for

*r*to be sufficiently low for the solution to be unstable, consistent with the dry solutions (

*r*= 1) being Fjørtoft stable (Fjørtoft 1950; Pedlosky 1964) because of the mismatch between the directions of the shear and PV gradients. Evidently, this mismatch is also sufficient to inhibit the growth of the moist modes unless

*r*is sufficiently low such that the dynamics is dominated by latent heating. Overall, these results confirm that the emergence of isolated DRV modes in the tilted model is not an artifact of having exactly zero PV gradients—although this is a useful limit to consider for theory—but rather generalizes to more realistic configurations that do include PV gradients. Again, we have repeated the calculations using a linear drag acting on the relative vorticity in the lower layer with a damping time scale of either 10 days (weak drag) or 2.5 days (strong drag). The wave–vortex transition persists, with the precise boundary largely unaffected by drag. Length scales are similar, but the growth rates are reduced (not shown).

Surprisingly, Fig. 9a suggests that the transition from wave to DRV regime in the two-layer model with equal and opposite PV gradients is independent of *r* (the boundary at *q _{y}*

_{crit}is entirely horizontal) but does rely on weakening or reversing the PV gradients relative to the classic untilted two-layer model. We investigate this result further by repeating the calculations with PV gradients that are not equal and opposite, but rather allowed to vary independently from each other, for two example values of

*r*(Figs. 9b,c). We observe that lowering of

*r*, as expected in a warmer climate whose stratification is closer to moist adiabatic, does make a difference since it produces DRVs as the fastest-growing solution for a larger range of PV gradients, particularly away from the diagonal line where the PV gradients are exactly equal and opposite. We also recall from the analytical solutions the existence of an upper bound on

*r*for DRVs to occur in an infinite domain when the PV gradients are zero (Fig. 6). Overall, we find that both weak PV gradients and weak moist static stability (small

*r*) can favor DRVs.

## 4. Comparison with storms at finite amplitude

Our two-layer theory for DRV modes and the warm-climate DRV mode in the idealized GCM calculation of O’Gorman et al. (2018) are both based on an assumption of small-amplitude disturbances. In this section, we analyze DRV storms in reanalysis and a warm-climate simulation of the idealized GCM to see how finite amplitude affects storm structure. We are particularly interested in the question of whether finite-amplitude effects can weaken the upper-level anticyclonic PV anomaly and possibly lead to DRVs that propagate but do not grow strongly because of the lower PV anomaly not having a strong enough upper PV anomaly with which to interact. Differences in the vertical structure of PV anomalies are expected at finite amplitude for two reasons. First, diabatic PV generation is weaker at finite amplitude in anticyclonic regions, an effect that is captured in semigeostrophic models (Hoskins 1975; Emanuel et al. 1987). Second, vertical advection of PV at finite amplitude can significantly offset negative diabatic generation at upper levels because positive PV anomalies at lower levels are advected upward.

We compare the structure of the warm-climate DRV mode (Fig. 1a) with a finite-amplitude DRV storm in the same idealized GCM and climate (Fig. 1c) and with two finite-amplitude storms from reanalysis. The first storm from reanalysis is the East Coast winter storm (Fig. 10a) that was identified as a propagating and moderately growing DRV by Moore et al. (2008). It later experienced explosive growth through interaction with a prexisting upper PV anomaly, but we consider the earlier diabatic phase. The second storm from reanalysis is an example midlatitude summer cyclone (Fig. 10c) from an updated version of the DRV climatology of Boettcher and Wernli (2013) that is based on ERA5 reanalysis (Hersbach et al. 2020). It was identified as a DRV by a tracking algorithm, selecting for substantial baroclinicity, sufficient moisture, fast propagation, and weak upper-level forcing. These example storms are meant to illustrate some of the variations in the constellation of PV anomalies and diabatic generation in observed DRVs.

The PV structure and PV generation rate due to latent heating of the (a),(b) 2005 winter storm (latitude 37° at 0000 UTC 25 Feb 2005) that was identified as a DRV in Moore et al. (2008), and (c),(d) an example summer cyclone (latitude 41.25° at 0500 UTC 10 Jul 2009) that was identified as a DRV in the climatology of Boettcher and Wernli (2013). Shading shows PV anomalies with respect to a 4-day moving average (using 6-hourly fields). Contours show (a),(c) the diabatic PV generation [the first term on the right-hand side of Eq. (30)] and (b),(d) the generalized diabatic PV generation including both diabatic PV generation and diabatic vertical advection as in Eq. (29). In all cases only diabatic effects from latent heating are included. Red contours are positive, and blue contours are negative, and the contour interval is 0.44 pvu h^{−1} for (a) and (b) and 0.10 pvu h^{−1} for (c) and (d). The zero contour line is not shown. All fields are calculated from ERA5 reanalysis, and PV is calculated using the hydrostatic approximation to Ertel’s PV.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

The PV structure and PV generation rate due to latent heating of the (a),(b) 2005 winter storm (latitude 37° at 0000 UTC 25 Feb 2005) that was identified as a DRV in Moore et al. (2008), and (c),(d) an example summer cyclone (latitude 41.25° at 0500 UTC 10 Jul 2009) that was identified as a DRV in the climatology of Boettcher and Wernli (2013). Shading shows PV anomalies with respect to a 4-day moving average (using 6-hourly fields). Contours show (a),(c) the diabatic PV generation [the first term on the right-hand side of Eq. (30)] and (b),(d) the generalized diabatic PV generation including both diabatic PV generation and diabatic vertical advection as in Eq. (29). In all cases only diabatic effects from latent heating are included. Red contours are positive, and blue contours are negative, and the contour interval is 0.44 pvu h^{−1} for (a) and (b) and 0.10 pvu h^{−1} for (c) and (d). The zero contour line is not shown. All fields are calculated from ERA5 reanalysis, and PV is calculated using the hydrostatic approximation to Ertel’s PV.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

The PV structure and PV generation rate due to latent heating of the (a),(b) 2005 winter storm (latitude 37° at 0000 UTC 25 Feb 2005) that was identified as a DRV in Moore et al. (2008), and (c),(d) an example summer cyclone (latitude 41.25° at 0500 UTC 10 Jul 2009) that was identified as a DRV in the climatology of Boettcher and Wernli (2013). Shading shows PV anomalies with respect to a 4-day moving average (using 6-hourly fields). Contours show (a),(c) the diabatic PV generation [the first term on the right-hand side of Eq. (30)] and (b),(d) the generalized diabatic PV generation including both diabatic PV generation and diabatic vertical advection as in Eq. (29). In all cases only diabatic effects from latent heating are included. Red contours are positive, and blue contours are negative, and the contour interval is 0.44 pvu h^{−1} for (a) and (b) and 0.10 pvu h^{−1} for (c) and (d). The zero contour line is not shown. All fields are calculated from ERA5 reanalysis, and PV is calculated using the hydrostatic approximation to Ertel’s PV.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Ertel PV anomalies are defined with respect to a zonal mean for the idealized GCM (using once-daily fields for the mode and 6 hourly for the macroturbulent state) and with respect to a 4-day moving average for the 6-hourly fields from ERA5 reanalysis forecasts. The forecast mode is chosen because it provides the temperature tendencies necessary for the calculation of latent heating. The fields are first interpolated from model to pressure levels, for both GCM and reanalysis fields, prior to calculating the PV and PV generation rates.

*Q*is the potential vorticity,

*θ*is the potential temperature,

^{5}Including vertical advection of PV is important because it can strongly offset diabatic PV generation (Büeler and Pfahl 2017; Lackmann 2002; Martínez-Alvarado et al. 2016; Stoelinga 1996; Persson 1995; Pfahl et al. 2015; Wernli and Davies 1997), and its inclusion clearly improves the match between the vertical structures of PV anomalies and PV generation for the small-amplitude DRV mode in the idealized GCM (cf. Figs. 1a,b). In addition, considering diabatic vertical advection as part of the diabatic source of PV makes a closer connection to our QG theory in which the pseudo-PV is not advected in the vertical.

We only consider diabatic effects due to latent heating. For the idealized GCM, we have confirmed that there is no convective precipitation in the region of the finite-amplitude DRV storm, and

The finite-amplitude DRVs from reanalysis (Fig. 10) do not extend as high in the atmosphere as the DRV in the warm climate of the idealized GCM, and this is as expected given that they occur in the current climate in which tropopause is lower and latent heating occurs lower in the troposphere. The generalized diabatic PV generation (contours in Figs. 10b,d) is noticeably smaller in magnitude for the upper-level negative generation rates as compared with lower-level positive generation rates. From the point of view of Eq. (30), the magnitude of diabatic PV generation is reduced in the upper anticylonic region as compared with the lower cyclonic region in which *Q* is larger, an effect that has been seen before in the context of warm conveyor belts (Joos and Wernli 2012), and the negative diabatic PV generation is also offset by upward diabatic advection of positive PV from the positive PV anomaly lower in the atmosphere. Alternatively, from the point of view of Eq. (29), the factor of *Q*^{2} tends to be much smaller in magnitude in anticyclonic regions as compared with cyclonic regions of a finite-amplitude storm. As a result, the upper-level negative PV anomaly is weaker in magnitude than the lower-level positive PV anomaly, especially in the case of the winter storm in which it is difficult to clearly identify a negative upper-level PV anomaly that is diabatically generated.

The reason the winter storm has a greater asymmetry between lower and upper PV anomalies as compared with the summer storm seems to be because it is a stronger storm (which emphasizes the finite-amplitude effects) but also because of more subtle effects related to its vertical profile of

The finite-amplitude DRV in the warm-climate simulation of the idealized GCM (Figs. 1c,d) shows some similarities to the corresponding small-amplitude mode (Figs. 1a,b), although the upper-level negative PV anomaly and the generalized diabatic PV generation are considerably reduced in the storm relative to the mode for vertical levels between 300 and 500 hPa and the positive PV anomaly extends higher, both as a result of diabatic vertical advection of PV. Negative generation of PV is nonetheless strong in the upper troposphere near the tropopause, and this seems to be because of diabatic vertical advection up the mean vertical PV gradient at those levels.

Overall, our analysis of finite-amplitude DRV storms shows that finite-amplitude effects must be taken into account to relate the structure of PV anomalies and diabatic generation in observed DRVs particularly for the upper-level PV anomalies. Our results also show the value of combining diabatic PV generation and diabatic vertical advection in a generalized diabatic PV generation diagnostic (following Hoskins et al. 1985), especially when trying to connect to simpler QG models and modal solutions.

## 5. Conclusions

We have analyzed a moist two-layer QG model with conditional latent heating and tilted boundaries and shown that it is capable of producing a DRV mode. The emergence of a DRV solution in a minimal model retaining the essential physics of baroclinicity and moisture clarifies the physical mechanisms involved and allows us to derive the first analytical expressions for the growth rate and horizontal length scale of DRVs.

A key step in our approach is the tilting of the model boundaries at a slope equal to the mean isentropes, which makes the two-layer model an analog of the interior of the Eady model in which dry-baroclinic instability has been shut off but moist instabilities are still possible. This allowed us to obtain a pure DRV solution within a conceptually simple two-layer model. PV-budget analysis revealed two distinct dynamical regimes. In the ascending branch, growth of the anomalies was maintained by diabatic heating partly offset by zonal advection, while in the descending branch growth was maintained solely by zonal advection.

We went on to derive the analytical dispersion relation for the growth rate and horizontal length scale of a DRV on an infinite domain, a significant novelty of this paper. The governing equation for the vertical velocity in the DRV is similar to the equation for the vertical velocity of moist baroclinic waves (Emanuel et al. 1987; Zurita-Gotor 2005) except for the presence of two extra terms that lead to isolated rather than periodic solutions—a distinctive characteristic of the DRV. Analytic solutions to the dispersion equations were found in the limit of small static-stability reduction factor (i.e., in the limit in which the stratification is neutral to moist convection). While the ascent length remains the same for the DRV as for the moist wave solutions of Emanuel et al. (1987) in this limit, the DRV grows faster by 54% as compared with the moist wave. This faster growth is consistent with the fact that the DRV emerged as the fastest-growing solution in the moist baroclinic instability simulations of O’Gorman et al. (2018) in a warm climate with small moist static stability. Root solving of the dispersion equations for a larger range of *r* values showed that physical solutions cease to exist when *r* > 0.38. This is an indication that isolated DRV disturbances cannot exist on an infinite domain when the moist static stability is not small enough.

Including nonzero meridional PV gradients in the tilted two-layer model and varying their strengths and varying the moist static stability (as represented by the reduction factor *r*), we showed that isolated DRV solutions emerge even in more realistic model setups and are not an artifact of the assumption of zero PV gradients in our simplest version of the tilted two-layer model. The most unstable mode transitions from periodic waves to isolated DRVs when the magnitude of the PV gradients is weakened or entirely reversed relative to the standard two-layer setup. This suggests that the vertical structure of meridional PV gradients may be an important additional factor that helps to determine DRV genesis zones in addition to small moist static stability. Weak QG PV gradients can be found particularly at polar latitudes in the current climate, which could help strengthen the links that have been previously established between the growth mechanism of DRVs and polar lows (Montgomery and Farrell 1991, 1992; Moore and Montgomery 2005; Moreno-Ibáñez et al. 2021).

The stark transition from periodic wave solutions to isolated DRV disturbances when diabatic heating becomes more important than meridional PV advection was also explained qualitatively using “PV thinking”: in a diabatically dominated regime, the asymmetry of the diabatic heating regenerates only positive PV anomalies in the lower layer and negative PV anomalies in the upper layer. However, a series of like-signed PV anomalies in each layer leads to weaker meriodional flow between the PV anomalies than occurs for a single PV dipole consisting of one anomaly in each layer. Thus the single PV dipole has stronger ascent and latent heating and emerges as the fastest-growing mode of the system.

Last, we compared the structure of small-amplitude DRV modes with finite-amplitude storms from reanalysis in winter and summer and from a warm-climate simulation in an idealized GCM. The finite-amplitude storms have similarities with the small-amplitude DRV modes but also some differences. In the storms from reanalysis, the upper-level negative PV anomaly is substantially weaker than the lower-level positive PV anomaly. This asymmetry arises because diabatic PV generation is weaker in anticyclonic regions at finite amplitude, and also because upward PV advection from the positive PV anomaly at lower levels can offset the upper-level negative PV generation. For the finite-amplitude DRV in the warm climate of the idealized GCM, vertical advection of the mean vertical PV gradient near the tropopause meant there was still a strong upper-level negative PV anomaly. In the case of the winter storm, the upper-level PV anomaly was sufficiently weak that it was difficult to identify, and we hypothesize that a weak upper-level PV anomaly may explain why some DRV storms in the current climate propagate but do not grow strongly. This hypothesis could be tested in future work with a semigeostrophic model that has sufficient vertical levels to accurately resolve vertical PV advection, and by tracking DRVs across a range of climates in idealized GCM simulations to study the relationship between growth rates and the structure of the PV anomalies.

Future work could also investigate whether a DRV solution and wave–vortex transition can also be isolated within a continuous Eady model in which dry baroclinic instability is eliminated by removing the upper lid. This setup is likely no longer tractable analytically (because solutions are no longer separable in the presence of nonlinear heating; Zurita-Gotor 2005), but a numerical analysis would make for a useful extension of this work, in which realistic features such as near surface temperature advection, vertically dependent drag and vertically dependent static stability reduction factor could be more readily incorporated.

The most unstable modes were calculated in O’Gorman et al. (2018) through repeated rescaling of perturbations to small amplitude, assuming upward motion to be saturated, and using a basic state equal to the zonal and time-mean of a fully nonlinear simulation for that climate.

Equations (4) and (5) are identical to the two-layer moist QG equations (Zurita-Gotor 2005) except for the addition of the terms −*τ _{x}* and −

*ϕ*in Eq. (4) and Eq. (5), respectively, which arise because of the tilted boundaries, and except for the presence of the mean radiative cooling term

_{x}*w*in which case any mean radiative cooling term drops out for an untilted model. We will see shortly, however, that the mean radiative cooling does not drop out when forming the

*w*equation for the tilted model.

The dimensional growth rate for *r* → 0 is 1.62*Uf* /(*NH*), and the dimensional half-ascent length is *π**NH*/*f*, where *H* is the depth of one layer and the vertical shear is 2*U*/*H*.

Note that the solutions classified as DRVs at *r* = 1 could instead be considered to be waves. From the dry dispersion relation, the wavelength of these most unstable dry modes becomes infinite (not shown), and thus one maximum in *w* is found numerically in the domain no matter how large of a domain is chosen. Note also that the pure DRV solution with *q*_{1}* _{y}* =

*q*

_{2}

*= 0 at*

_{y}*r*= 0.9 is shown to be stable but would grow weakly on a larger domain.

An alternative approach of including vertical advection of PV in pressure coordinates (rather than diabatic vertical advection) gives similar results except that there can be additional vertical advection of PV in the upper troposphere and stratosphere in regions where latent heating is small.

## Acknowledgments.

We acknowledge helpful discussions with Santiago Benavides, Andre da Silva, Minmin Fu, Martín Vélez-Pardo, Stephan Pfahl, Kerry Emanuel, Glenn Flierl, Maxi Boettcher, Heini Wernli, and Raymond Pierrehumbert. We thank Maxi Boettcher and Heini Wernli for providing us with the database of tracked DRVs from reanalysis. We are grateful to an anonymous reviewer, Heini Wernli, and Michael T. Montgomery for their helpful comments. We acknowledge support from NSF AGS 2031472 and the mTerra Catalyst Fund.

## Data availability statement.

ERA5 data are available from ECMWF (https://cds.climate.copernicus.eu), and code for the tilted two-layer model calculations is available on GitHub (https://github.com/matthieukohl/DRV_Paper).

## APPENDIX A

### Derivation of the Tilted Two-Layer Model

*f*plane [Eqs. (1) and (2)] in the vertical taking into account the tilted boundaries

*h*

_{1}(

*y*) and

*h*

_{2}(

*y*) through a modified boundary condition on

*w*at the top and bottom (Fig. A1). For the vorticity equations in the two layers we obtain

*ϕ*and baroclinic streamfunction

*τ*:

*x*,

*y*∼

*L*, where

_{D}*L*=

_{D}*NH*/(

*f*) is the deformation radius with individual layer height

*H*,

*z*∼

*H*,

*u*and

*υ*∼

*U*,

*ϕ*=

*τ*∼

*UL*,

_{D}*t*∼

*L*/

_{D}*U*, and

*h*

_{1}and

*h*

_{2}∼

*UH*/(

*fL*) to obtain

_{D}*τ*

_{0}= −

*y*and

*ϕ*

_{0}= 0 corresponding to a shear flow

*u*

_{1}= 1 and

*u*

_{2}= −1 such that the advection terms are linearized. Note, however, that the thermodynamic equation remains nonlinear because of the nonlinear dependence of

*r*(

*w*) on

*w*. Using the tilted boundary conditions

*h*

_{1}=

*h*

_{2}=

*y*and assuming that the perturbations are independent of

*y*gives the perturbation Eqs. (4)–(6). If instead we assume

*h*

_{1}=

*α*

_{1}

*y*and

*h*

_{2}=

*α*

_{2}

*y*, we obtain the perturbation Eqs. (26)–(28).

Schematic of the tilted two-layer model in the *x*–*z* plane.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Schematic of the tilted two-layer model in the *x*–*z* plane.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

Schematic of the tilted two-layer model in the *x*–*z* plane.

Citation: Journal of the Atmospheric Sciences 79, 10; 10.1175/JAS-D-22-0022.1

## APPENDIX B

### Derivation of the Equation for *w*

*σ*and taking one derivative, and substituting for

*τ*and

_{xx}*τ*using Eq. (10), we find

_{xxxx}*ϕ*−

_{xxx}*ϕ*can be substituted using Eq. (B1) to give

_{x}*σ*

^{2}from Eq. (B5), we get

*ϕ*terms in Eq. (B6) to finally obtain Eq. (12), which is an equation in

*w*only.

## APPENDIX C

### Mass Conservation for a DRV on an Infinite Domain

*L*→ ∞. It is important that we took the integral in

*x*to obtain Eq. (C1) prior to taking the limit

*L*→ ∞ (i.e., taking the limit of the statement of mass conservation) because the order of taking the limit and integral affects the result. As shown in section 2c, the definition of

*a*together with mass conservation implies that

*a*→ 0 if we want solutions for which

*b*and

*w*

_{↑}remain bounded. This leaves us with the indeterminate term

*aL*in the mass conservation equation [Eq. (C1)]. We can eliminate this term by using

*aL*:

*L*→ ∞ and

*a*→ 0 to arrive at Eq. (21) in the main text. For completeness we note that by substituting Eq. (21) into Eq. (C2), and taking the limit

*L*→ ∞ and

*a*→ 0, we obtain the expression

## APPENDIX D

### Derivation of the Dispersion Relation

*c*

_{1}and

*c*

_{2}. Writing Eqs. (18) and (20) as

*c*sin

_{i}*k*and

_{i}b*c*cos

_{i}*k*, with

_{i}b*i*= 1, 2, by each other, we obtain the two “tan” equations that form the dispersion relation [Eqs. (22) and (23)] in section 2c.

## APPENDIX E

### Asymptotic Expressions for Growth Rate and Half-Ascent Length

*r*≪ 1, the wavenumbers can be simplified to

*k*

_{1}= 1/

*O*(

*k*

_{2}=

*O*(

*r*), where we assume that

*σ*=

*σ*

_{0}is an

*O*(1) quantity to be determined. Plugging these expressions into the two tangent Eqs. (22) and (23) that form the dispersion relation we obtain to leading order:

*r*becomes small. To balance it, we use the

*ansatz b*/

*π*/2 +

*ϵ*, where

*ϵ*→ 0 as

*r*→ 0, from which we obtain

*ϵ*=

*σ*

_{0}) and

*b*= (

*π*/2)

*r*(1 +

*σ*

_{0}) such that to leading order

*b*is

*O*(

*O*(1), the right-hand side must vanish. Thus 1 −

*σ*

_{0}(

*σ*

_{0}− 1) = 0, which implies for growing solutions

*σ*

_{0}> 0 that

## REFERENCES

Bembenek, E., D. N. Straub, and T. M. Merlis, 2020: Effects of moisture in a two-layer model of the midlatitude jet stream.

,*J. Atmos. Sci.***77**, 131–147, https://doi.org/10.1175/JAS-D-19-0021.1.Boettcher, M., and H. Wernli, 2013: A 10-yr climatology of diabatic Rossby waves in the Northern Hemisphere.

,*Mon. Wea. Rev.***141**, 1139–1154, https://doi.org/10.1175/MWR-D-12-00012.1.Booth, J. F., L. Polvani, P. A. O’Gorman, and S. Wang, 2015: Effective stability in a moist baroclinic wave.

,*Atmos. Sci. Lett.***16**, 56–62, https://doi.org/10.1002/asl2.520.Bretherton, F. P., 1966: Baroclinic instability and the short wavelength cut-off in terms of potential vorticity.

,*Quart. J. Roy. Meteor. Soc.***92**, 335–345, https://doi.org/10.1002/qj.49709239303.Büeler, D., and S. Pfahl, 2017: Potential vorticity diagnostics to quantify effects of latent heating in extratropical cyclones. Part I: Methodology.

,*J. Atmos. Sci.***74**, 3567–3590, https://doi.org/10.1175/JAS-D-17-0041.1.Charney, J. G., 1947: The dynamics of long waves in a baroclinic westerly current.

,*J. Atmos. Sci.***4**, 136–162, https://doi.org/10.1175/1520-0469(1947)004<0136:tdolwi>2.0.co;2.Eady, E. T., 1949: Long waves and cyclone waves.

,*Tellus***1**(3), 33–52, https://doi.org/10.3402/tellusa.v1i3.8507.Emanuel, K. A., M. Fantini, and A. J. Thorpe, 1987: Baroclinic instability in an environment of small stability to slantwise moist convection. Part I: Two-dimensional models.

,*J. Atmos. Sci.***44**, 1559–1573, https://doi.org/10.1175/1520-0469(1987)044<1559:BIIAEO>2.0.CO;2.Fantini, M., 1995: Moist Eady waves in a quasigeostrophic three-dimensional model.

,*J. Atmos. Sci.***52**, 2473–2485, https://doi.org/10.1175/1520-0469(1995)052<2473:MEWIAQ>2.0.CO;2.Fjørtoft, R., 1950:

*Application of Integral Theorems in Deriving Criteria of Stability for Laminar Flows and for the Baroclinic Circular Vortex*. Vol. 17, Grøndahl and Søns Boktr, 52 pp.Hersbach, H., and Coauthors, 2020: The ERA5 global reanalysis.

,*Quart. J. Roy. Meteor. Soc.***146**, 1999–2049, https://doi.org/10.1002/qj.3803.Hoskins, B. J., 1975: Geostrophic momentum approximation and the semi-geostrophic equations.

,*J. Atmos. Sci.***32**, 233–242, https://doi.org/10.1175/1520-0469(1975)032<0233:TGMAAT>2.0.CO;2.Hoskins, B. J., M. E. McIntyre, and A. W. Robertson, 1985: On the use and significance of isentropic potential vorticity maps.

,*Quart. J. Roy. Meteor. Soc.***111**, 877–946, https://doi.org/10.1002/qj.49711147002.Joos, H., and H. Wernli, 2012: Influence of microphysical processes on the potential vorticity development in a warm conveyor belt: A case-study with the limited-area model COSMO.

,*Quart. J. Roy. Meteor. Soc.***138**, 407–418, https://doi.org/10.1002/qj.934.Lackmann, G. M., 2002: Cold-frontal potential vorticity maxima, the low-level jet, and moisture transport in extratropical cyclones.

,*Mon. Wea. Rev.***130**, 59–74, https://doi.org/10.1175/1520-0493(2002)130<0059:CFPVMT>2.0.CO;2.Lambaerts, J., G. Lapeyre, and V. Zeitlin, 2012: Moist versus dry baroclinic instability in a simplified two-layer atmospheric model with condensation and latent heat release.

,*J. Atmos. Sci.***69**, 1405–1426, https://doi.org/10.1175/JAS-D-11-0205.1.Lapeyre, G., and I. M. Held, 2004: The role of moisture in the dynamics and energetics of turbulent baroclinic eddies.

,*J. Atmos. Sci.***61**, 1693–1710, https://doi.org/10.1175/1520-0469(2004)061<1693:tromit>2.0.co;2.Martínez-Alvarado, O., S. L. Gray, and J. Methven, 2016: Diabatic processes and the evolution of two contrasting summer extratropical cyclones.

,*Mon. Wea. Rev.***144**, 3251–3276, https://doi.org/10.1175/MWR-D-15-0395.1.Montgomery, M. T., and B. F. Farrell, 1991: Moist surface frontogenesis associated with interior potential vorticity anomalies in a semigeostrophic model.

,*J. Atmos. Sci.***48**, 343–368, https://doi.org/10.1175/1520-0469(1991)048<0343:msfawi>2.0.co;2.Montgomery, M. T., and B. F. Farrell, 1992: Polar low dynamics.

,*J. Atmos. Sci.***49**, 2484–2505, https://doi.org/10.1175/1520-0469(1992)049<2484:PLD>2.0.CO;2.Moore, R. W., and M. T. Montgomery, 2004: Reexamining the dynamics of short-scale, diabatic Rossby waves and their role in midlatitude moist cyclogenesis.

,*J. Atmos. Sci.***61**, 754–768, https://doi.org/10.1175/1520-0469(2004)061<0754:RTDOSD>2.0.CO;2.Moore, R. W., and M. T. Montgomery, 2005: Analysis of an idealized, three-dimensional diabatic Rossby vortex: A coherent structure of the moist baroclinic atmosphere.

,*J. Atmos. Sci.***62**, 2703–2725, https://doi.org/10.1175/JAS3472.1.Moore, R. W., M. T. Montgomery, and H. C. Davies, 2008: The integral role of a diabatic Rossby vortex in a heavy snowfall event.

,*Mon. Wea. Rev.***136**, 1878–1897, https://doi.org/10.1175/2007MWR2257.1.Moreno-Ibáñez, M., R. Laprise, and P. Gachon, 2021: Recent advances in polar low research: Current knowledge, challenges and future perspectives.

,*Tellus***73A**, 1890412, https://doi.org/10.1080/16000870.2021.1890412.O’Gorman, P. A., 2011: The effective static stability experienced by eddies in a moist atmosphere.

,*J. Atmos. Sci.***68**, 75–90, https://doi.org/10.1175/2010JAS3537.1.O’Gorman, P. A., T. M. Merlis, and M. S. Singh, 2018: Increase in the skewness of extratropical vertical velocities with climate warming: Fully nonlinear simulations versus moist baroclinic instability.

,*Quart. J. Roy. Meteor. Soc.***144**, 208–217, https://doi.org/10.1002/qj.3195.Parker, D. J., and A. J. Thorpe, 1995: Conditional convective heating in a baroclinic atmosphere: A model of convective frontogenesis.

,*J. Atmos. Sci.***52**, 1699–1711, https://doi.org/10.1175/1520-0469(1995)052<1699:CCHIAB>2.0.CO;2.Pedlosky, J., 1964: The stability of currents in the atmosphere and the ocean: Part I.

,*J. Atmos. Sci.***21**, 201–219, https://doi.org/10.1175/1520-0469(1964)021<0201:TSOCIT>2.0.CO;2.Persson, P. O. G., 1995: Simulations of the potential vorticity structure and budget of fronts 87 IOP8.

,*Quart. J. Roy. Meteor. Soc.***121**, 1041–1081, https://doi.org/10.1002/qj.49712152506.Pfahl, S., P. A. O’Gorman, and M. S. Singh, 2015: Extratropical cyclones in idealized simulations of changed climates.

,*J. Climate***28**, 9373–9392, https://doi.org/10.1175/JCLI-D-14-00816.1.Phillips, N. A., 1954: Energy transformations and meridional circulations associated with simple baroclinic waves in a two-level, quasi-geostrophic model.

,*Tellus***6**, 273–286, https://doi.org/10.1111/j.2153-3490.1954.tb01123.x.Snyder, C., and R. S. Lindzen, 1991: Quasi-geostrophic wave-CISK in an unbounded baroclinic shear.

,*J. Atmos. Sci.***48**, 76–86, https://doi.org/10.1175/1520-0469(1991)048<0076:QGWCIA>2.0.CO;2.Stoelinga, M. T., 1996: A potential vorticity-based study of the role of diabatic heating and friction in a numerically simulated baroclinic cyclone.

,*Mon. Wea. Rev.***124**, 849–874, https://doi.org/10.1175/1520-0493(1996)124<0849:APVBSO>2.0.CO;2.Tierney, G., D. J. Posselt, and J. F. Booth, 2018: An examination of extratropical cyclone response to changes in baroclinicity and temperature in an idealized environment.

,*Climate Dyn.***51**, 3829–3846, https://doi.org/10.1007/s00382-018-4115-5.Wernli, H., and H. C. Davies, 1997: A Lagrangian-based analysis of extratropical cyclones. I: The method and some applications.

,*Quart. J. Roy. Meteor. Soc.***123**, 467–489, https://doi.org/10.1002/qj.49712353811.Wernli, H., S. Dirren, M. A. Liniger, and M. Zillig, 2002: Dynamical aspects of the life cycle of the winter storm ‘Lothar’ (26–26 December 1999).

,*Quart. J. Roy. Meteor. Soc.***128**, 405–429, https://doi.org/10.1256/003590002321042036.Whitaker, J. S., and C. A. Davis, 1994: Cyclogenesis in a saturated environment.

,*J. Atmos. Sci.***51**, 889–908, https://doi.org/10.1175/1520-0469(1994)051<0889:ciase>2.0.co;2.Zurita-Gotor, P., 2005: Updraft/downdraft constraints for moist baroclinic modes and their implications for the short-wave cutoff and maximum growth rate.

,*J. Atmos. Sci.***62**, 4450–4458, https://doi.org/10.1175/JAS3630.1.