Moist Static Potential Vorticity Budget in Tropical Motion Systems

Ángel F. Adames aDepartment of Atmospheric and Oceanic Sciences, University of Wisconsin–Madison, Madison, Wisconsin

Search for other papers by Ángel F. Adames in
Current site
Google Scholar
PubMed
Close
https://orcid.org/0000-0003-3822-5347
,
Rosa M. Vargas Martes aDepartment of Atmospheric and Oceanic Sciences, University of Wisconsin–Madison, Madison, Wisconsin

Search for other papers by Rosa M. Vargas Martes in
Current site
Google Scholar
PubMed
Close
,
Haochang Luo bDepartment of Climate and Space Science and Engineering, University of Michigan, Ann Arbor, Michigan

Search for other papers by Haochang Luo in
Current site
Google Scholar
PubMed
Close
, and
Richard B. Rood bDepartment of Climate and Space Science and Engineering, University of Michigan, Ann Arbor, Michigan

Search for other papers by Richard B. Rood in
Current site
Google Scholar
PubMed
Close
Open access

Abstract

Analyses of simple models of moist tropical motion systems reveal that the column-mean moist static potential vorticity (MSPV) can explain their propagation and growth. The MSPV is akin to the equivalent PV except it uses moist static energy (MSE) instead of the equivalent potential temperature. Examination of an MSPV budget that is scaled for moist off-equatorial synoptic-scale systems reveals that α, the ratio between the vertical gradients of latent and dry static energies, describes the relative contribution of dry and moist advective processes to the evolution of MSPV. Horizontal advection of the moist component of MSPV, a process akin to horizontal MSE advection, governs the evolution of synoptic-scale systems in regions of high humidity. On the other hand, horizontal advection of dry PV predominates in a dry atmosphere. Derivation of a “moist static” wave activity density budget reveals that α also describes the relative importance of moist and dry processes to wave activity amplification and decay. Linear regression analysis of the MSPV budget in eastern Pacific easterly waves shows that the MSPV anomalies originate over the eastern Caribbean and propagate westward due to dry PV advection. They are amplified by the fluxes of the moist component of MSPV over the Caribbean Sea and over the eastern Pacific from 105° to 130°W, underscoring the importance of moist processes in these waves. On the other hand, dry PV convergence amplifies the waves from 90° to 100°W, likely as a result of the barotropic energy conversions that occur in this region.

© 2022 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Ángel F. Adames, angel.adamescorraliza@wisc.edu

Abstract

Analyses of simple models of moist tropical motion systems reveal that the column-mean moist static potential vorticity (MSPV) can explain their propagation and growth. The MSPV is akin to the equivalent PV except it uses moist static energy (MSE) instead of the equivalent potential temperature. Examination of an MSPV budget that is scaled for moist off-equatorial synoptic-scale systems reveals that α, the ratio between the vertical gradients of latent and dry static energies, describes the relative contribution of dry and moist advective processes to the evolution of MSPV. Horizontal advection of the moist component of MSPV, a process akin to horizontal MSE advection, governs the evolution of synoptic-scale systems in regions of high humidity. On the other hand, horizontal advection of dry PV predominates in a dry atmosphere. Derivation of a “moist static” wave activity density budget reveals that α also describes the relative importance of moist and dry processes to wave activity amplification and decay. Linear regression analysis of the MSPV budget in eastern Pacific easterly waves shows that the MSPV anomalies originate over the eastern Caribbean and propagate westward due to dry PV advection. They are amplified by the fluxes of the moist component of MSPV over the Caribbean Sea and over the eastern Pacific from 105° to 130°W, underscoring the importance of moist processes in these waves. On the other hand, dry PV convergence amplifies the waves from 90° to 100°W, likely as a result of the barotropic energy conversions that occur in this region.

© 2022 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Ángel F. Adames, angel.adamescorraliza@wisc.edu

1. Introduction

It is well known that large-scale (≥1000 km across) tropical motion systems organize and are themselves modulated by deep convection (Kiladis et al. 2006, 2009; Feng et al. 2020). These so-called convectively coupled tropical motion systems play an important role in the hydrologic cycle of the tropics, occasionally leading to floods (Lubis and Jacobi 2015; Hunt and Fletcher 2019; Dominguez et al. 2020). They can also lead to tropical cyclone formation (Maloney and Hartmann 2000; Gall et al. 2010; Blake and Kimberlain 2013; Núñez Ocasio et al. 2020). Some of these systems exhibit slow propagation and strong vorticity anomalies. In this group we find equatorial Rossby waves, easterly waves, monsoon low pressure systems, the Madden–Julian oscillation (MJO), as well as other tropical depression–like disturbances (Yasunaga and Mapes 2012; Inoue et al. 2020).

Numerous studies into these slowly propagating systems have revealed that tropospheric water vapor plays an important role in their evolution (Raymond and Fuchs 2009; Rydbeck and Maloney 2015; Fuchs-Stone et al. 2019; Adames and Maloney 2021). There is increasing evidence that the MJO has properties of a moisture mode, a type of tropical disturbance whose thermodynamics are largely governed by water vapor [Raymond and Fuchs 2009; see Adames and Maloney (2021) for further review]. New studies also suggest that convectively coupled equatorial Rossby waves may exhibit properties of moisture modes (Gonzalez and Jiang 2019; Fuchs-Stone et al. 2019). Water vapor also explains the majority of variance in rainfall in east Pacific easterly waves (Wolding et al. 2020). It may also play an important role in monsoon low pressure systems (Diaz and Boos 2019; Clark et al. 2020; Adames 2021), and in African easterly waves (Russell et al. 2020).

In spite of the growing evidence of water vapor’s importance in slowly evolving tropical motion systems, our understanding on how moisture and convection interact with the large-scale circulation remains incomplete (Grabowski and Moncrieff 2004; Ahmed and Neelin 2018). Nonetheless, recent studies based on simple moist models of the atmosphere suggest that we may gain insight onto these interactions by examining a quantity defined as “gross” potential vorticity (GPV) (Adames and Ming 2018; Adames 2021). The GPV combines moist static energy (MSE) and PV into a single quantity that is quasi-conserved under moist adiabatic, inviscid, and hydrostatic processes. Furthermore, the processes that lead to the propagation and growth of waves in their simple models were succinctly described by the terms in the GPV equation.

More recently, Adames (2021) hypothesized that GPV may be physically related to the equivalent PV (Emanuel 1983; Martin et al. 1992) or, more specifically, a quantity that resembles equivalent PV but defined with MSE. While the former is applied to shallow-water and two-layer models of the atmosphere, the latter could be used to diagnose the processes that lead to the evolution of systems in the real atmosphere. The goal of this study is to examine this hypothesis. We will show that a quantity referred to as the moist static PV (MSPV) does indeed capture many of the features of GPV and quantifies the relative role of dry and moist processes in the evolution of off-equatorial tropical motion systems. We will also demonstrate the usefulness of the MSPV budget by applying it to eastern Pacific easterly waves (PEWs).

The remainder of this study is structured as follows. Section 2 reviews the shallow-water model of Adames and Ming (2018, hereafter AM18) and discusses the physical relevance of GPV. Section 3 defines the MSPV as the three-dimensional analog of GPV and discusses the processes that lead to its evolution. Scale analysis of the MSPV budget is shown in section 4. The MSPV budget is linearized and interpreted in section 5. An example of the potential applications of MSPV and its budget is shown in section 6. The conclusions are offered in section 7.

2. Evolution of tropical motion systems in a linear shallow-water model with prognostic moisture

a. Review of the AM18 model

Using a set of shallow-water equations on a beta plane with prognostic moisture, AM18 derived a quantity they defined as GPV. We can write the GPV as follows:
qG=M˜ζ+f0s^p(CpT+Lq),
where ζ′ is the anomalous relative vorticity, T′ is the anomalous temperature, q′ is the specific humidity anomaly, L = 2.5 × 106 J kg−1 is the latent energy of vaporization, Cp = 1005 J kg−1 K−1 is the specific heat at constant pressure, f0 is the planetary vorticity at a reference latitude, s^p can be thought of as a shallow-water dry static stability, and M˜ is the dimensionless moist static stability (referred to as the gross moist stability). The variable M˜ is defined as
M˜=s^pLq^ps^p=1α^,
where q^p can be thought of as a shallow-water equivalent of a vertical moisture gradient, referred to as the gross moisture stratification (Yu et al. 1998). The hats denote that they are shallow-water mean fields, and their units are different from the equivalent fields that will be shown in the next section. The nomenclature of the shallow-water model has been slightly modified from AM18 to facilitate comparison with the subsequent sections. For the same purpose we have also defined α^=Lq^p/s^p as the shallow-water climatological-mean Chikira (2014) parameter. For reference, the relevant variables discussed in this section are summarized in Table 1.
Table 1

The variables and definitions used in section 2.

Table 1
It is worth noting that, to a reasonable approximation, CpT′ can be replaced by s′ in Eq. (1), where s′ is the dry static energy (DSE) anomaly. By doing this, the gross PV can be written as qGM˜ζ+f0s^p1m, where m′ = s′ + Lq′ is the anomalous MSE. We can define the shallow-water moist static stability perturbation as M˜=ms^p1, which leads to
qGM˜ζ+f0M˜.
This definition of GPV will be particularly helpful in establishing the similarities between GPV and the moist static PV discussed in the following section.
One of the useful features of GPV comes from its decomposition into dry and moist components:
qG=M˜qd+(1M˜)qm,
where qd and qm are the dry and moist quasigeostrophic (QG) potential vorticity anomalies, respectively, written as follows:
qd=ζ+f0CpTs^p,
qm=f0(qq^p+CpTs^p).
While qd is simply the shallow-water equivalent of the PV in QG theory (Holton and Hakim 2012), the physical meaning of qm warrants further discussion. Inspection of Eq. (4c) shows that qm is reminiscent to moist enthalpy or MSE, except that q′ and CpT′ are divided by their mean vertical gradients. AM18 showed that the budget equation for qm is similar to that of MSE. Thus, we can interpret qm as describing the evolution of the moist thermodynamics of a wave. A similar quantity to qm was derived by Lapeyre and Held (2004) [see their Eq. (2)].

AM18 found that when M˜0 (or α^1), which can be thought to represent a saturated troposphere, moist enthalpy dominates the evolution of convectively coupled waves. Conversely, when M˜1 (or α^0), the atmosphere is dry and the GPV equation simplifies to the dry QG PV equation. Thus, the GPV can be thought of as a generalization of QG PV that includes the effects of moist thermodynamics.

The usefulness of using qG′ comes from evaluation of its budget. In AM18, the GPV budget was written as
qGt=υβG,
where υ′ is the meridional wind anomaly, βG is the GPV gradient, the weighted sum of the mean meridional dry and moist PV gradients:
βGM˜βd+(1M˜)βm.
The dispersion relation of moist Rossby waves can be obtained from Eq. (5) if the moisture anomalies are neglected, in which case it can be written as follows:
ω=βGkM˜(k2+l2)+kd2,
where kd is the inverse of the Rossby radius of deformation, k is the zonal wavenumber, l is the meridional wavenumber and ω is the wave’s angular frequency. Examining Eq. (7) reveals that the numerator on the right-hand side (rhs) of Eq. (7) is the rhs term of Eq. (5). Meanwhile, the denominator term on the rhs term of Eq. (7) is the left-hand side of Eq. (5) (AM18). When moisture is included and precipitation is parameterized in terms of it, the dispersion relation will contain an imaginary term that is positive when both the moist and dry PV tendency are in the same direction, a mechanism of growth AM18 termed as “moisture-vortex” instability. AM18 also found that systems that grow from moisture-vortex instability exhibit an in-phase component between qG and the advection of mean GPV by the anomalous meridional wind.

An example of how moist and dry PV advection contribute to the evolution of a convectively coupled Rossby wave is shown in Fig. 1. For M˜=0.2, moist PV advection contributes nearly twice as much to GPV evolution as dry PV advection does, even though βm is less than half the magnitude of βd in this figure. Thus, the value of M˜ is critical to determining the importance of moist processes in the evolution of a convectively coupled system.

Fig. 1.
Fig. 1.

Contribution of (top) dry PV advection M˜υβd and (middle) moist PV advection (1M˜)υβm to the (bottom) GPV advection υβG. The contoured fields are (top) M˜qd, (middle) (1M˜)qm, and (bottom) qG. Values of M˜=0.2, βd = 2.2 × 10−11 m−1 s−1, and βm = 0.8 × 10−11 m−1 s−1 were used. In the three panels the advection has been multiplied by 1011 to facilitate interpretation. The contour interval for all panels is 2 × 10−6 s−1.

Citation: Journal of the Atmospheric Sciences 79, 3; 10.1175/JAS-D-21-0161.1

b. Wave activity evolution

The aforementioned results indicate that qG can provide information about the dynamics of convectively coupled systems that cannot be obtained by examining dry PV alone. To examine this possibility, we can turn Eq. (5) into an equation for the evolution of wave activity akin to those that describe the Eliassen–Palm (EP) flux (see chapter 10 in Vallis 2017). We will multiply all terms in Eq. (5) by qG/βG and define the shallow wave activity density (or more accurately, its pseudo-momentum) as follows:
AqG22βG,
which leads to the following equation for the evolution of wave activity:
At=υqG.
Equation (9) indicates that wave activity increases where there is a flux of anomalous GPV by the anomalous meridional winds. We can write this equation in terms of the convergence of a wave activity flux, which we break down into dry and moist components:
At=M˜Cd+(1M˜)Cm,
Cdυqd,
Cm −υqm.
Equation (10) reveals that the evolution of A is due to the weighted sum of the convergences of dry and moist wave activity (Cd and Cm, respectively). In a moist atmosphere, where M˜0, the evolution of wave activity is determined by qm. Conversely, in a dry atmosphere, M˜1, and dry wave activity governs the evolution of wave activity.

By examining Fig. 1 we observe that qm is shifted westward toward meridional wind anomalies of the opposite sign. This phasing indicates that there is a net positive flux of qm when averaged over the horizontal domain of Fig. 1, indicating that Cm > 0. In contrast, qd is in quadrature with υ′. Thus these waves are growing from moist wave energy convergence alone. This result indicates that the moisture-vortex instability mechanism described in AM18 can be understood in terms of Cm.

3. Moist static potential vorticity as an analog of GPV

The results of the previous section suggest that GPV can summarize the dynamics of convectively coupled tropical motion systems in a shallow-water model. However, the GPV equation cannot be readily applied to observations, reanalysis and GCM output because of the underlying shallow-water and linearity assumptions. In this section, we seek to obtain a quantity that is reminiscent to the GPV but applicable to the stably stratified atmosphere. Table 2 summarizes the main variables that will be used from this section onward.

Table 2

The main variables and definitions used in this study from section 3 onward. Note that the PVU is defined as 10−3 J m2 kg−2 s−1.

Table 2

A candidate field variable that exhibits many similarities to GPV is the moist static Ertel PV (MSPV). In isobaric coordinates, the MSPV can be written as follows:
Pms=gωam,
where g = 9.81 m s−2 is the gravitational acceleration;
ωa=ηi+ξj+ζak
is the absolute vorticity vector, where ζa = ζ + f is the vertical component of the absolute vorticity;
=ix+jy+kp
is the gradient operator in isobaric coordinates; and
m=CpT+Φ+Lq
is the MSE. The MSPV in Eq. (11) is very similar to the equivalent PV defined in Martin et al. (1992); Moore and Lambert (1993) and McCann (1995), except that the equivalent potential temperature (θe) is replaced by MSE. Both MSE and θe are approximately conserved in moist adiabatic motions. Specifically, MSE is conserved under hydrostatic, moist adiabatic processes. In deriving the MSE, the temporal tendency and horizontal advection of pressure, as well as nonhydrostatic pressure perturbations are neglected (Madden and Robitaille 1970; Betts 1974). For large-scale tropical motions these terms are negligible, but can be large in systems with strong winds and pressure perturbations such as mature tropical cyclones. In spite of these differences, we choose to use MSE because it allows us to make a more direct comparison with the GPV defined in the previous section. Furthermore, the fact that MSE is defined as a sum of dry static energy and latent energy will allow for a straightforward decomposition of the MSPV into its dry and moist components, as was done for the GPV. Because m is used instead of θe, the units of MSPV differ from those of equivalent PV. We address this issue by defining the PV unit (PVU) as 10−3 J m2 kg−2 s−1. This will cause the MSPV to have the same PVUs as equivalent PV.

It is worth noting that in isobaric coordinates k points toward the surface. It is also worth pointing out that the vertical components of both and ωa have different units than their horizontal components. While this notation may cause some confusion, it makes Eq. (12) more compact and easily comparable to PV in height and isentropic coordinates. The same notation has been used by previous studies that analyzed PV in isobaric coordinates (Zhang and Ling 2012; Russell et al. 2020).

a. MSPV budget

The conservation equation for MSPV can be written as follows:
DPmsDt=G+BFr,
where
DDt=t+ux+υy+ωp
is the material derivative in isobaric coordinates;
B=gfvgphm
is the baroclinic (solenoidal) generation of MSPV, written as the horizontal advection of MSE by the thermal wind as in Martin et al. (1992), where vg = ugi + υgj is the geostrophic wind vector;
Fr=gm(×Fr)
is the frictional dissipation of MSPV, where Fr is the frictional dissipation of momentum; and
G=gωam˙
is the diabatic generation of MSPV, where
m˙DmDt=Qrωm¯p
describes the sources and sinks of MSE, where Qr is the radiative heating rate and pωm¯ is the vertical divergence of the turbulent flux of MSE.

Equation (12a) is similar to the one shown by Martin et al. (1992), except it includes the frictional dissipation term and we replace θe with MSE. An alternate form of Eq. (12a), in which all the rhs terms are expressed in the form of a forcing vector, is shown in the first section of the appendix.

b. Interpretation of MSPV

By definition, the MSPV is the absolute vorticity that is normal to the MSE surfaces. The definition of MSPV differs from conventional Ertel PV in that m replaces the potential temperature (θ). However, this replacement leads to differences in the properties and physical interpretation of the two PV quantities. For inviscid, adiabatic processes, dry PV is conserved following the flow. A simplified, linear version of PV behaves as the wave equation for Rossby waves. Thus, dry PV not only serves as a tracer, it also summarizes the dynamics of inviscid, adiabatic geophysical flows.

In contrast, MSPV is not conserved under inviscid, moist adiabatic processes, as in the equivalent PV discussed by Schubert et al. (2001). Furthermore, because pm > 0 in the lower troposphere, it follows that MSPV tends to be negative in the lower troposphere when the absolute vorticity is positive (Fig. 2). In contrast, dry PV is positive in this layer for the same vorticity anomaly since pθ < 0 in the free troposphere. Furthermore, the fact that the tropical lower troposphere is potentially unstable (pm > 0, Fig. 3a) can lead to the development of deep moist convection. Undiluted rising parcels within deep convection approximately conserve MSE. If we assume that the mass circulation associated with the undiluted updraft is bounded by two MSE surfaces of the same value in the upper and lower troposphere (Riehl and Malkus 1958), it follows that the updraft will be associated with a constant MSE surface that connects the lower and upper troposphere, as shown in Fig. 3b.1 Because MSE surfaces are impermeable to MSPV in the same way that isentropes are impermeable to Ertel PV (Haynes and McIntyre 1987, 1990), it follows that the isobaric MSPV concentration of rising parcels is constrained to be enclosed within the MSE surfaces described in Fig. 3b. Thus, in the absence of large sources of MSPV, deep tropical convection acts to redistribute MSPV in the troposphere. Furthermore, because the MSE gradient reverses sign in the mid troposphere (Fig. 3a), it follows that parcels must rotate anticyclonically as they exit the updraft in the upper troposphere. The transport of MSPV from the lower to the upper troposphere occurs at time scales that are much shorter than the time scale of synoptic-scale motion systems (Fig. 3).

Fig. 2.
Fig. 2.

Schematic describing MSPV quasi-conservation. For inviscid, adiabatic processes, a vortex filament will remain enclosed between two moist isentropic surfaces as it follows the flow.

Citation: Journal of the Atmospheric Sciences 79, 3; 10.1175/JAS-D-21-0161.1

Fig. 3.
Fig. 3.

Schematic describing how MSPV is redistributed in a synoptic-scale system that exhibits cyclonic rotation in the lower troposphere in a moist atmosphere with no planetary vorticity, so that Pmsgζpm. (a) Because the lower troposphere is potentially unstable, deep moist convection can develop that rises at the moist adiabatic lapse rate (approximately conserving m). (b) The convection will transport MSPV from the lower to the upper troposphere as vortex filaments rise following the surfaces of constant MSE. (c) Because the MSE gradient reverses sign in the midtroposphere, it follows that parcels must rotate anticyclonically in order to quasi-conserve MSPV.

Citation: Journal of the Atmospheric Sciences 79, 3; 10.1175/JAS-D-21-0161.1

Because of the vertical redistribution of MSPV by convection, it follows that the vertically averaged MSPV ( Pms) is the actual quasi-conserved variable, as is also the case with MSE (Andersen and Kuang 2012). We define the column-averaged MSPV as
Pms=1ΔpptpsPmsdp,
where pt = 100 hPa is the tropopause pressure, ps = 1000 hPa is the surface pressure, and Δp = 900 hPa is the depth of the troposphere. We choose to perform column-averaging rather than column integrating so that the units of the MSPV tendency are in the form of PVU s−1.

4. Scaling of the MSPV budget for off-equatorial synoptic-scale tropical motion systems

A comparison of Eq. (12) with Eq. (5) reveals that the MSPV and GPV budgets are quite different. This is a result of the assumptions made in the formulation of the AM18 model as well as the linearization. A budget that is more reminiscent of Eq. (5) can be obtained from scaling an linearization of Eq. (12a). In this section we will scale Eq. (12a) for convectively coupled synoptic-scale tropical motion systems with a vertical scale that spans the entire troposphere. While MSPV can be applied even near the equator, we will emphasize off-equatorial tropical motion systems in order to compare the MSPV budget more directly to the GPV discussed in section 2. Magnitudes for the different field variables are shown in Table 3. The scales were obtained using ERA5 data and closely follow the scale analysis done by Adames et al. (2021). Scaling of the MSPV reveals the following magnitudes:
Pms= −gηmx106gξmy106gζamp105,
which indicates that the vertical component of MSPV governs its distribution. Thus, the MSPV can be expressed to leading order as follows:
Pms −gζamp.
We can use the simplified expression for MSPV in Eq. (15) to obtain the scales of the material derivative of MSPV:
DPmsDt=Pmst1010+uPmsx1010+υPmsy1010+ωPmsp1011,
which indicates that the evolution of MSPV is governed by the Eulerian derivative and the horizontal advection terms.
Table 3

Scales used for the different field variables.

Table 3
We can scale the baroclinic term B by assuming that the geostrophic wind is of similar magnitude to the horizontal winds (∼10 m s−1). Through the thermal wind equation, such a wind could be obtained by assuming a horizontal temperature gradient of 1 K (1000 km)−1. The resulting scaling is written as follows:
B=gfvgphm1011,
which shows that it is much smaller than the MSPV tendency.
Last, scale analysis of the MSPV generation term yields the following magnitudes:
G= −gηm˙x1011gξm˙y1011gζam˙p1011.
By neglecting G, Fr, and B, and dropping the vertical advection term yields the following approximate MSPV budget:
DhPmsDt0,
where
DhDt=t+vh
is the horizontal material derivative and v = ui + υj is the horizontal wind vector. Equation (19a) reveals that MSPV is quasi-conserved following the horizontal flow in synoptic-scale tropical motion systems.

It is important to note that the scaling presented here only applies to synoptic-scale systems. The leading-order terms are likely to be different if the analysis is performed on disturbances whose spatial and/or temporal scale differ from those assumed in this section. Furthermore, because the leading-order terms are only an order of magnitude larger than the other terms, it follows that small departures from the scaling assumptions of this section can cause the smaller terms to become nonnegligible. Nonetheless, we will use the scaled MSPV budget as a starting point to understand the evolution of moist tropical motion systems.

5. Linearized MSPV

a. Linearized MSPV and its relation to GPV

We can linearize MSPV with respect to a slowly varying background state. Anomalies of MSPV will be denoted by primes, while overbars depict the background state. By expanding Pms into its mean and anomaly fields, and dropping the nonlinear terms yields the following:
Pms= −g(ζm¯p+ζ¯amp),
where
m¯p=Lq¯ps¯pmp=Lqpsp
is the mean and anomalous vertical MSE gradient, respectively, where
  s¯p=s¯psp= −sp
q¯p=q¯p  qp=qp
are the mean and anomalous static stability and moisture stratification, where s = CpT + Φ is the dry static energy. With these definitions we can also define the dry and moist components of the MSPV:
Pms=(1α¯)Pd+α¯Pm,
where the dry PV anomaly is written as
Pd=g(ζs¯p+ζ¯asp),
and the moist component of the MSPV is
Pm=P¯d(sps¯pqpq¯p)
and
α¯=Lq¯ps¯p
The term α can be thought of as a measure of potential instability since a layer is potentially unstable when α > 1 (see the second section of the appendix). Note that in a dry atmosphere α¯=0, so that MSPV approximately becomes the Ertel PV. During moist adiabatic processes, where parcels conserve m as they ascend through the troposphere α = 1. In this case Pd plays no role in the evolution of convectively coupled waves, and it is Pm that drives the evolution of the system. As Eq. (22c) indicates, Pm describes the relative change in the stratification of DSE and latent energy. If m¯p<0, increasing the moisture stratification will act to further decrease mp, and thus leading to a more negative MSPV assuming that ζ¯a>0. Similarly, Pm can also be thought of as a measure of relative changes in the stability of the column since Pm=P¯dα/α¯. Similar variables to Pm have been obtained in other studies that analyze variants of MSPV (Stechmann and Majda 2006; Smith and Stechmann 2017; Wetzel et al. 2020).
Equations (22) and (22) can be thought as the full atmosphere equivalents of Eqs. (3) and (4a). We can show that
Pmsgs¯pqG.
The main difference between MSPV and GPV are in their dimensional units. In spite of having different dimensions, the two variables have the same physical interpretation.

b. Anomalous MSPV budget

After scaling and linearization, the equation for the evolution of the MSPV anomalies can be written as
D¯hPmsDt = vhP¯ms,
where
D¯hDt=t+v¯h.
We can decompose P¯ms into the dry and moist PV gradient by noting that P¯ms can be rewritten as
P¯ms=(1α¯)P¯d.
By multiplying Eq. (26) by h and after some rearrangement of the terms, we can express the horizontal MSPV gradient as
hP¯ms=(1α¯)βd+α¯βm,
where
βdhP¯d
βmP¯d(hs¯ps¯phq¯pq¯p)
are the dry and moist PV gradients, respectively. We use a similar notation to the moist and dry PV gradients in section 2 to suggest their physical connection, although the PV gradients in this section are vectors with different units (cf. Tables 1 and 2).
With the decomposition of hP¯ms, we can write the anomalous MSPV budget as
D¯hPmsDt=(1α¯)vβdα¯vβm,
where the horizontal advection terms are now decomposed into the advection of dry and moist PV. Note that Eq. (30) can be thought as equivalent to Eq. (5) but more general. A similar equation can be obtained without linearization, as shown in the third section of the appendix.
The column-averaged budget takes the following form:
Pmst=vhPms=v¯hPms(1α¯)vβdα¯vβm,
where the contribution of the background flow to the advection of MSPV is now shown explicitly on the rhs of Eq. (31).

c. The moist static wave activity equation

One feature of PV is that it can describe the propagation of wave energy. In the midlatitudes, the propagation of wave energy can be described in terms of the EP flux (Andrews and McIntyre 1976) and its more generalized forms that have been derived thereafter (Plumb 1985; Takaya and Nakamura 2001). We can obtain an analogous equation but for MSPV by following Takaya and Nakamura (1997). We multiply Eq. (24a) by Pms/|hP¯ms|. Assuming that the horizontal MSPV gradient varies in space more slowly than Pms2, we can define the equivalent wave activity density (or more accurately, the equivalent pseudo-momentum) as follows:
A=12Pms2|hP¯ms|,
and turn Eq. (24a) into an equation for wave activity evolution:
At=Ct,
where Ct is the total wave activity convergence:
Ct −PmsvhP¯ms|hP¯ms|hv¯A.
Note that for a dry atmosphere, Ct can be written as the convergence of a wave activity flux as in Takaya and Nakamura (2001), Ct = −F. It also describes the EP flux after zonal averaging. The EP flux can be extended to include the effects of water vapor for a near-saturated atmosphere (Stone and Salustri 1984; Dwyer and O’Gorman 2017). However, such an assumption would lead to large inaccuracies at the time scales that we are interested in since water vapor fluctuates significantly in the tropics while temperature does not. Thus, it would be more appropriate to derive an EP flux that is applicable for an unsaturated atmosphere. However, we were unable to obtain such an expression. As a result, we keep the notation in terms of a wave activity convergence.
As was done in the previous section, we can break down Ct into contributions from dry and moist processes:
Ct=(1α¯)Cd+α¯Cmhv¯A,
where we define the convergence of dry and moist wave activity fluxes as follows:
CdPdvhP¯ms|hP¯ms|,
CmPmvhP¯ms|hP¯ms|,
noting that vPd=Fd, where Fd is the dry wave activity flux, or equivalently the extended EP flux (Plumb 1985; Takaya and Nakamura 2001). With these definitions and after column-averaging the wave activity tendency equation becomes:
At=hv¯A+(1α¯)Cd+α¯Cm
Equation (38) shows the relative importance of moist and dry wave energy convergence (Cd and Cm, respectively) on the evolution of wave activity. The dry wave energy convergence is given by the flux of dry PV, as in the EP flux formulation, except that the zonal and meridional fluxes are weighted by xP¯ms/|hP¯ms| and yP¯ms/|hP¯ms|, respectively. The moist wave energy convergence can be thought to be similar to a weighted MSE flux. It is worth noting that the wave activity equation largely mirrors the MSPV budget in Eq. (31). Both equations contain a term that is weighted by (1α¯), a term that is weighted by α¯ and a term involving the mean flow. In spite of the similarities between the two equations, a process that is dominant in Eq. (31) might not contribute significantly to Eq. (38). For example, it is possible for (1α¯)vβd to dominate the evolution of MSPV while α¯Cm to dominate the evolution of wave activity.

6. MSPV evolution in east Pacific easterly waves

a. Data and methods

To show an example of the usefulness of MSPV, we will investigate its variability in PEWs. We will make use of data from the fifth reanalysis from the European Centre for Medium-Range Weather Forecasts (ECMWF) (ERA5; Hersbach et al. 2019). We used instantaneous fields with a time interval of 12 h, spanning the 40-yr interval of 1979–2018. The ERA5 data used have horizontal resolution of 0.5° × 0.5° and 27 vertical levels from 1000 to 100 hPa. We make use of the following ERA5 pressure level fields: specific humidity (q), temperature (T), geopotential (Φ), vorticity (ζ), and the zonal and meridional winds (u and υ, respectively). Additionally, we make use of ERA5 column-integrated water vapor (CWV), precipitation rate (P), and outgoing longwave radiation (OLR).

Anomaly fields are obtained by removing the seasonal mean at each grid point for the months of July–September (JAS), when easterly wave activity in the Northern Hemisphere is the strongest (e.g., Mekonnen et al. 2006; Cheng et al. 2019). A Lanczos filter (Duchon 1979) is performed on the anomalous data to retain easterly wave-related variability on the order of 2–6 days (Mekonnen et al. 2006; Leroux et al. 2010; Cheng et al. 2019). From here on, primed variables (e.g., Pms) will describe 2–6-day filtered anomalies while fields with an overline (e.g., P¯ms) denote JAS averages. While it would be more accurate to calculate the overlined fields using a 6-day low-pass filter, we use a simple climatological average since our main interest is to showcase the potential usefulness of the MSPV budget. A more rigorous calculation of the MSPV budget applied to easterly waves will be shown in a future paper.

The variability associated with PEWs is extracted by performing an EOF analysis on the filtered OLR data for the months of JAS, over a domain that ranges from 10°S–30°N and from 150° to 75°W. The first two EOFs of OLR variability in the eastern Pacific explain 3.16% and 3.01% of the filtered OLR variance, respectively. This procedure mirrors that of Rydbeck and Maloney (2014, 2015), except that we choose to perform the analysis on OLR rather than 700-hPa vorticity and our range of months is narrower. In spite of these differences, our resulting patterns are similar to those found by Rydbeck and Maloney (2014, 2015). The leading EOF patterns that correspond to PEWs are insensitive to the choice of domain or field variable, although the PEWs explain a larger fraction of the total variance if a smaller domain is chosen. They also explain a larger fraction of total variance if a spatially smoother field such as geopotential is used. However, we opt to use OLR since we are specifically interested in the moist processes associated with PEWs and how it impacts the MSPV budget.

The figures shown in this section are obtained through regression analysis of the first two principal components (PCs) of the aforementioned EOF analysis. We will refer to these as the PCs of PEW activity. Following Adames and Wallace (2014), the regression maps are obtained via the following equation:
D=N1SPT,
where D is the linear regression pattern in dimensional units, S is a two-dimensional matrix representing a field variable (e.g., Pms), P is the nondimensional PC, and N is the sample size. Statistical significance of the regression maps is obtained by applying a two-tailed t-test at the 95% level under the null hypothesis that the regression maps are uncorrelated with the PCs. We have defined the contour interval in such a way that all the contoured and shaded fields presented are statistically significant. In plots where vector arrows are shown, only arrows that are statistically significant are shown.

b. MSPV evolution in PEWs

Maps of anomalous geopotential height z, horizontal winds v, precipitation P, and Pms regressed upon the first PC of anomalous OLR are shown in Fig. 4a. The maps are lag regressions that show the evolution of PEWs from lag day −2 to lag day 2. Negative (positive) lags represent the days before (after) PC1 attains a maximum amplitude. At lag day −2, an anomalous anticyclone collocated with suppressed rainfall can be seen southwest of Mexico. Both the z and P anomalies associated with these features exhibit a southwest to northeast tilt, a common feature of easterly waves in this region (Serra et al. 2008; Rydbeck and Maloney 2015). The anticyclone is flanked by regions of enhanced rainfall.

Fig. 4.
Fig. 4.

(a)–(e) Anomalous precipitation (shading), 700-hPa geopotential height anomalies (contours) and the horizontal wind field (arrows) regressed onto the leading PC of 2–6-day OLR variance over the eastern Pacific basin (10°S–30°N, 150°–75°W). Each panel is a lag regression with the lag days indicated in the title. Contour interval is every 0.5 m starting at 1 m. Positive z anomalies are shown as solid contours while negative anomalies are shown as dashed contours. (f)–(j) As in (a)–(e), but showing the column-integrated water vapor (CWV′) as shading and the column-averaged MSPV anomalies as contours. Contour interval is every 1.5 × 10−2 PVU starting at 3 × 10−2 PVU. (k)–(o) As in (a)–(e), but showing anomalous horizontal MSPV advection as shading and the MSPV tendency as contours. Contour interval is every 0.0025 PVU day−1 starting at 0.005 PVU day−1.

Citation: Journal of the Atmospheric Sciences 79, 3; 10.1175/JAS-D-21-0161.1

Another anomalous cyclone is seen over the Caribbean Sea near the Dominican Republic. From examination of the z anomalies, it is unclear whether this cyclone is related to the PEWs. It is clearly separated from the PEWs and exhibits a broader horizontal structure that does not tilt with longitude. A region of weak suppressed rainfall is also seen to the east of the Yucatan Peninsula in association with an anticyclone that dissipated over this region.

At lag day −1 (Fig. 4b), the anomalous anticyclone propagated slightly to the northwest. It has strengthened slightly, although the suppressed precipitation anomalies that are associated with it have weakened. An anomalous cyclone has developed to the east of the anticyclone nearly in phase with the enhanced rainfall. Suppressed rainfall anomalies have now appeared south of Guatemala, possibly in association to the suppressed rainfall anomalies that were previously east of the Yucatan Peninsula. The anomalous cyclone over the Caribbean has propagated westward, and has reduced its distance to the PEWs. It has also weakened substantially. To the east of the Caribbean cyclone, an anomalous anticyclone has appeared.

At lag day 0 (Fig. 4c), a new anticyclone is seen developing over the eastern Pacific south of Oaxaca. At this time the anomalous cyclone over the Caribbean has dissipated. Only a weak cyclonic flow is seen in association with it to the east of the Yucatan Peninsula as well as weak precipitation anomalies. However, the anticyclone to the west has amplified. At lag day 1, the weak cyclonic flow is centered over and to the South of the Yucatan Peninsula, and its features are reminiscent of a Central American Gyre (Papin et al. 2017). The regression maps at lag day 2 mirror those of lag day −1 with reversed polarity. The only notable difference is that circulation anomalies over the Caribbean are no longer evident.

The results of Figs. 4a–e suggest that westward-propagating anomalies exist both in the eastern Pacific and over the Caribbean. However, it is not clear from these maps whether the anomalies at the two locations are physically related. Examination of Pms (Figs. 4f–j) reveals what appears to be a single wave train of Pms that extends from 130° to 60°W. Anomalies in Pms develop over the eastern Caribbean and propagate westward. While the z anomalies substantially weaken as they approach the Yucatan Peninsula, the Pms anomalies do not. Instead, the Pms anomalies approach and cross central America over the course of two days. The emergence of the Pms anomalies over the eastern Pacific appears to be related to the development of the geopotential and precipitation anomalies in this region, although the latter are shifted to the southeast of the Pms anomalies.

Positive anomalies in CWV′ precede the development of the Pms anomalies and are initially associated with anomalous southerly winds (Fig. 4f). Once the Pms anomalies develop they closely follow the evolution of the CWV. The CWV anomalies exhibit a similar horizontal structure to the Pms anomalies over the Caribbean, although they are shifted to the east of them. Such a zonal shift is expected and is also seen in the moist shallow-water waves shown in Fig. 1. That the Pms anomalies resemble the CWV anomalies over the Caribbean, combined with the weak signature in geopotential height observed in that region suggests that the Pms anomalies are largely the result of anomalies of the moist component of the MSPV (Pm). Over the eastern Pacific, the Pms anomalies match the location of the CWV anomalies near the coast of Mexico, but not farther south. This is likely a result of the Pms anomalies being scaled by the mean ambient absolute vorticity, which becomes small near the equator.

The propagation of the MSPV anomalies is largely the result of the anomalous horizontal advection of MSPV ( vhPms), as shown by the strong correspondence of the contours and shading in Figs. 4k–o. The shading is in phase with the contours at nearly all stages of the PEW life cycle, with only two exceptions. One exception is the time when the anomalies are developing over the eastern Caribbean near lag days −2 and 0. The other exception is during the decay stage of the wave over the eastern Pacific, when the contours have nearly vanished, indicating a weak to nonexistent MSPV tendency. The mismatch between the MSPV tendency and horizontal MSPV advection suggests that a nonnegligible amount of MSPV generation is occurring at this time. Nonetheless, the MSPV advection closely matches the MSPV tendency during the rest of the PEW life cycle, consistent with the scaling discussed in section 4.

Since horizontal MSPV advection is the dominant contributor to the propagation of the MSPV anomalies in PEWs, it is instructive to consider the processes that drive this advection. In Fig. 5, we show the three contributions to the anomalous horizontal MSPV advection. Dry PV advection (Figs. 5a,b) is largely in phase with the MSPV tendency throughout the domain. This phasing is robust at the times when both PC1 and PC2 become a maximum. The only time the correspondence is slightly weaker is during the initial stage of PEW development in the eastern Caribbean at the time when PC1 is a maximum, consistent with the overall weaker correspondence of total MSPV advection found in Fig. 4.

Fig. 5.
Fig. 5.

As in Fig. 4, but showing the anomalous horizontal advection of MSPV as contours and (a),(b) the horizontal advection of dry PV by the anomalous horizontal winds [ (1α)vβd]; (c),(d) the horizontal advection of the moist PV by the anomalous horizontal winds [ αvβd] and (e),(f) the total advection of mean MSPV by the anomalous horizontal winds [ vhP¯ms] as shading. Linear regressions onto (a),(c),(e) PC1 and (b),(d),(f) PC2 of eastern Pacific 2–6-day filtered OLR variability. Contour interval is every 0.0025 PVU day−1 starting at 0.005 PVU day−1.

Citation: Journal of the Atmospheric Sciences 79, 3; 10.1175/JAS-D-21-0161.1

Moist PV advection (Figs. 5a,b), on the other hand, is stronger to the south of the MSPV tendency, and does not align with the MSPV tendency, except over the far western portion of the domain at the time when PC1 is a maximum. Instead, it acts to slow down propagation to the east of 100°W. Last, advection of MSPV by the mean flow contributes to the propagation of the MSPV anomalies over the Caribbean Sea, but does not contribute significantly to MSPV propagation over the eastern Pacific.

We can summarize the results from Fig. 5 by projecting the individual advective terms in Eq. (31) to the Pms tendency. This projection is obtained as in Clark et al. (2020) by calculating the integral of the advective term with tPms over a region with an area A. The result is then divided by the integral of the square of tPms:
Proj(X,tPms)=AXtPmsdAAtPmstPmsdA,
where A is the region covering 5°–30°N, 130°–60°W, which encompasses the Caribbean Sea and eastern Pacific, and X is a term in Eq. (31). We perform this projection for both PC1 and PC2 of easterly wave activity, and average the result. The result of this calculation reveals that dry PV advection is indeed the largest contributor to MSPV propagation, explaining up ∼90% of the tendency (Fig. 6). Advection by the background winds also favors westward propagation, but its contribution is roughly 1/3 that of the dry PV advection. The overall contribution from advection of the moist component is small and acts against the westward propagation of the anomalies. A small residual in the budget is found, possibly a result of neglecting variations in the background state, or from nonlinear MSPV advection.
Fig. 6.
Fig. 6.

Fractional contribution of the terms in Eq. (31) to the propagation of the MSPV anomalies over the Caribbean and eastern Pacific region (5°–30°N, 130°–60°W). The bars are obtained by calculating each term for PC1 and PC2 and then averaging the two PCs together.

Citation: Journal of the Atmospheric Sciences 79, 3; 10.1175/JAS-D-21-0161.1

To understand how the PEWs grow over the Caribbean and eastern Pacific region, Fig. 7 shows the contribution of (1α¯)Cd and α¯Cm to the evolution of A over the region. Because we are interested in wave activity convergence over the PEW life cycle, Fig. 7 is constructed as an average of the two PCs that represent PEW activity. Examination of Fig. 7 reveals that dry wave activity convergence (1α¯)Cd is dominant near 95°W. In this region the meridional shear in the climatological-mean zonal wind is a maximum. Hence the dry wave energy convergence is possibly the result of barotropic energy conversion, as suggested by previous studies (Rydbeck and Maloney 2014; Molinari and Vollaro 2000). Over the Caribbean, moist energy conversion α¯Cm explains nearly all of the wave energy convergence. Convergence of A by the mean flow is comparable to the contributions from (1α¯)Cd and α¯Cm, and are largest over the eastern Pacific. Together with α¯Cm, it explains nearly all of the wave energy convergence from 105° to 130°W. When the three processes are considered together, it is clear that the amplitude of the total wave activity convergence Ct is consistent with the amplitude of MSPV advection shown in Figs. 4k–o. Thus, while the wave energy convergence explains the amplification of the MSPV anomalies as they propagate across the Caribbean Sea and eastern Pacific, it does not explain their appearance over the eastern Caribbean Sea.

Fig. 7.
Fig. 7.

Wave activity convergence from the anomalous winds during the PEW life cycle averaged over 10°–20°N. The red line shows the convergence of dry wave energy (1α¯)Cd, the blue line shows the convergence of moist wave energy α¯Cm, the green line shows hv¯A, and the black line shows the sum of the three terms, Ct. The lines are constructed by calculating each term for PC1 and PC2 and then averaging the two PCs together. For increased clarity, the data were smoothed using a 10° running mean.

Citation: Journal of the Atmospheric Sciences 79, 3; 10.1175/JAS-D-21-0161.1

7. Discussion and conclusions

Shallow-water models that include moist processes and prognostic moisture indicate that a quantity referred to as “gross potential vorticity” (GPV) can summarize the dynamics of moist waves (Adames and Ming 2018; Adames 2021). Guided by the shallow-water equations, we propose that the moist static potential vorticity (MSPV) may serve as an analog to GPV in a stratified atmosphere. The MSPV is similar to the equivalent PV, which has been used in studies of supercell thunderstorms (Rotunno and Klemp 1985), heavy precipitation bands (Clark et al. 2002), moist symmetric instability (Emanuel 1983; Martin et al. 1992) and in extratropical cyclones (Emanuel et al. 1987; Cao and Cho 1995). Furthermore, the budget equation for MSPV largely mirrors that of equivalent PV (Martin et al. 1992). Quantities similar to GPV and MSPV have also been derived in other idealized studies (Lapeyre and Held 2004; Monteiro and Sukhatme 2016). However, to the authors’ best knowledge, MSPV has not been extensively applied to the study of moist tropical motion systems.

A scale analysis of the MSPV budget was performed using scales that are typical for off-equatorial synoptic-scale tropical motion systems. The scaling indicates that horizontal advection of MSPV is the largest contributor to the MSPV tendency. This advection can be broken down into contributions from dry and moist processes, as in the shallow-water GPV. The relative contributions of these processes to the evolution of MSPV depends on α, the ratio between the vertical gradients of latent energy and dry static energy. In moist regions where the tropospheric-mean value of α is near unity, moist processes are expected to govern the evolution of MSPV, while dry processes are expected to play a minor role. In a dry troposphere (α → 0), the MSPV equation approximately becomes the traditional Ertel PV (referred to as the dry PV). The scale analysis indicates that other processes in the MSPV budget are substantially smaller than the horizontal advection. While these processes are neglected in this study, we acknowledge that they could be important when seeking a more complete understanding of moist systems.

The MSPV budget can also be used to obtain an equation for the evolution of wave activity. In a similar vein to the MSPV budget, the convergence of wave activity can be broken into dry and moist contributions, with α once again determining their relative contributions.

That α appears in the MSPV and equivalent wave activity budgets may be an indication of its physical importance to moist tropical motion systems. Previous studies have shown that α determines the sign and magnitude of diabatic processes in the moisture equation (Chikira 2014; Wolding and Maloney 2015; Wolding et al. 2016; Janiga and Zhang 2016) (see the second section of the appendix). In simple models, α is related to 1M˜, where M˜ is the normalized gross moist stability (section 2). Studies based on these models have shown that moisture governs the thermodynamics of tropical motion systems when M˜0 (Adames et al. 2019; Ahmed et al. 2021). Other studies have found that a realistic value of α is needed to simulate a robust MJO (Chikira 2014). The vertical profile of α is also important in understanding the MJO’s response to climate change (Bui and Maloney 2019; Maloney et al. 2019). Our results indicate that a realistic value of α is also needed in order to properly simulate moist tropical motion systems in general. More work is needed in order to fully understand the role of α in the tropics, and its implication for MSPV and moisture evolution.

Application of the MSPV budget to PEWs shows that MSPV advection governs the evolution of MSPV, consistent with the scaling presented in section 4. It also showed that the MSPV anomalies originate over the eastern Caribbean Sea. Horizontal MSPV advection does not explain the MSPV tendency at this point, indicating that the scaling assumptions in section 4 break down at this time. We hypothesize that at this time, MSPV is generated by the baroclinic term B. Over this region, vertical wind shear in association with the Caribbean low-level jet (Amador 1998; Muñoz et al. 2008) can zonally advect the moisture anomalies that are seen at this time (Fig. 4f), generating negative MSPV. The enhanced moisture at this time is likely due to anomalous southerly flow advecting humid air from northern South America into the Caribbean.

Once developed, the MSPV anomalies continue to amplify as they propagate westward across the Caribbean. Examination of the wave activity equation suggests that amplification of the waves at this time is due to the flux of the moist component of MSPV ( Pm), which can be interpreted as a convergence of moist wave activity. Once the MSPV anomalies arrive at the eastern Pacific, moist wave activity conversion diminishes while dry wave energy convergence significantly amplifies. At this time, the larger geopotential height anomalies associated with PEWs is observed. We attribute the large dry wave energy convergence from 105° to 85°W as the result of barotropic energy conversions, consistent with the findings of previous studies (Molinari et al. 1997; Molinari and Vollaro 2000; Rydbeck and Maloney 2014). Once the MSPV anomalies propagate to the west of ∼100°W moist convergence becomes the dominant source of wave activity once again.

Our results suggest that the life cycle of PEWs is at least partly tied to waves that arrive to the eastern Pacific from the Caribbean, in agreement with previous studies (Molinari and Vollaro 2000; Serra et al. 2010). However, we do not reject the possibility that some PEWs are excited in the eastern Pacific, as was shown by Rydbeck et al. (2017). Our results also indicate that both dry and moist processes play a key role in these waves, and both need to be considered in order to gain a complete understanding of their behavior.

In addition to PEWs, the column MSPV budget can be used to gain new insights onto the dynamics of moist tropical motion systems, specifically in understanding the relative role of moist and dry processes in their evolution. Application of the MSPV budget to monsoon low pressure systems (Hurley and Boos 2015) could lead to new insights to the dynamics of these systems. Recent studies suggest that these systems grow from a combination of barotropic instability and processes involving moist convection, possibly involving moisture-vortex instability (AM18; Adames 2021; Diaz and Boos 2019, 2021). A combination of moist convection and dry processes may also explain the dynamics of African easterly waves (Berry and Thorncroft 2012; Janiga and Thorncroft 2013; Russell et al. 2020).

Analysis of the MSPV budget could also provide insights onto the processes that lead to tropical cyclone (TC) genesis. A recent study by Murthy and Boos (2018) showed that the radial gradient of surface latent heat fluxes was important for the formation of tropical depressions. By inspecting Eq. (12a), we can see that horizontal gradients in surface heat fluxes are included in the MSPV generation term G. Furthermore, studies by Raymond et al. (2007) and Tang (2017a,b), suggest that a horizontal gradient in the normalized gross moist stability, i.e., the normalized divergence of column MSE fluxes, is also important for TC genesis. We hypothesize that this contribution is implicit in the flux convergence of MSPV. Thus, the MSPV budget elucidates the physical processes that lead to TC genesis. However, we note that there is evidence that the MSPV in a TC evolves from being negative for weak TCs to being approximately zero in mature TCs (Peng et al. 2019), consistent with the notion of slantwise moist neutrality in these systems (Emanuel 1986). Future work should examine how MSPV evolves both during cyclogenesis and during intensification.

There are some caveats to the use of MSPV. First, MSPV is not conserved for moist adiabatic processes since it contains a baroclinic generation term (Schubert et al. 2001). Furthermore, it cannot be inverted without some additional information about the moisture field. For example, Smith and Stechmann (2017) show that a similar quantity to MSPV, referred to precipitating QGPV can be jointly inverted with a quantity M that includes additional information about the water vapor field. While PV definitions that use the virtual potential temperature may still be preferable when looking for an invertible quantity, we posit that the evolution of MSPV offers a more complete summary of the dynamics and thermodynamics of moist tropical motion systems. Nonetheless, we acknowledge that MSPV is a complex quantity whose behavior differs significantly from Ertel PV. It is unclear whether the benefit of using a variable that summarizes the dynamics of moist systems outweighs the potential hindrances that its complexity brings.

In conclusion, the column-averaged MSPV budget has the potential to further our understanding of tropical motion systems. The budget could be applied to systems that exhibit such characteristics such as African easterly waves, monsoon low pressure systems, tropical depressions, among others. It can also be used in GCM intercomparison studies to analyze how well models capture the relative contribution of dry and moist processes to the evolution of these systems. All of these have the potential to be fruitful directions for future research.

Acknowledgments.

ÁFA and HL were supported by the National Science Foundation’s Grant AGS-1841559. RVM was supported by National Science Foundation Graduate Research Fellowship Program under Grant DGE-1747503. Support was also provided by the Graduate School and the Office of the Vice Chancellor for Research and Graduate Education at the University of Wisconsin–Madison with funding from the Wisconsin Alumni Research Foundation. ÁFA would like to thank Kuniaki Inoue and Hannah Zanowski for conversations that helped in the interpretation of the MSPV budget.

APPENDIX

Additional Details on the Interpretation of the MSPV Budget

a. Alternate form of the MSPV budget

An interesting feature of Eq. (12a) is that all the terms on the rhs can be written in terms of the divergence of a forcing vector as in Haynes and McIntyre (1987):
Pmst=Jms,
where
Jms=uPms+g(ωam˙mfvgpm×Fr)
is the MSPV forcing. The forcing Jms consists of the MSPV flux, MSPV generation along vortex lines, MSE flux by the thermal wind, and frictional dissipation. In steady state, Jms describes the circulation of MSPV.

b. Further discussion on α

The term α was originally described by Chikira (2014), who used it to understand the evolution of moisture in the MJO cycle. Through the application of the WTG approximation, the equation for moisture conservation can be written as
DhLqDt(1α)Qc+αQrLωq¯p,
where Qc is the convective heating, Qr is the radiative heating, and ωq¯ is the vertical eddy flux of moisture. Examination of the WTG moisture equation reveals that radiative heating becomes more effective at moistening the troposphere as α increases. Furthermore, increasing α reduces the net drying by latent heat release via compensating vertical moisture advection (Adames and Maloney 2021). When α > 1 latent heat release moistens the atmosphere instead of drying it. With all other processes held constant, a value of α near unity implies that convection is less efficient at drying the troposphere. This occurrence tends to be related with a near zero “gross moist stability” (Neelin and Held 1987; Raymond and Fuchs 2009; section 2), a condition that leads to the existence of moisture modes.
The physical relevance of α can be considered more generally by considering the criterion for potential instability. In height coordinates, a layer is potentially unstable if zθe < 0. In isobaric coordinates, a layer is potentially unstable if mp > 0, where mp = Lqpsp. By rearranging the terms we can show that the moist static stability can be written as
mp=sp(α1)

Examination of Eq. (A4) reveals that a layer is potentially unstable when α > 1 and it is stable when α < 1.

c. Scaled MSPV equation without linearization

The horizontal MSPV can be decomposed into dry and moist contributions even without linearization. Noting that Pms(1α)Pd, and expanding the horizontal gradient in Pms yields the following expression:
Pmst(1α)vβdαvβm,
where βd and βm are as in Eq. (30), but for the full q and s fields, not their background-mean components.

REFERENCES

  • Adames, Á. F., 2021: Interactions between water vapor, potential vorticity, and vertical wind shear in quasi-geostrophic motions: Implications for rotational tropical motion systems. J. Atmos. Sci., 78, 903923, https://doi.org/10.1175/JAS-D-20-0205.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Adames, Á. F., and J. M. Wallace, 2014: Three-dimensional structure and evolution of the MJO and its relation to the mean flow. J. Atmos. Sci., 71, 20072026, https://doi.org/10.1175/JAS-D-13-0254.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Adames, Á. F., and Y. Ming, 2018: Interactions between water vapor and potential vorticity in synoptic-scale monsoonal disturbances: Moisture vortex instability. J. Atmos. Sci., 75, 20832106, https://doi.org/10.1175/JAS-D-17-0310.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Adames, Á. F., and E. D. Maloney, 2021: Moisture mode theory’s contribution to advances in our understanding of the Madden-Julian oscillation and other tropical disturbances. Curr. Climate Change Rep., 7, 7285, https://doi.org/10.1007/s40641-021-00172-4.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Adames, Á. F., D. Kim, S. K. Clark, Y. Ming, and K. Inoue, 2019: Scale analysis of moist thermodynamics in a simple model and the relationship between moisture modes and gravity waves. J. Atmos. Sci., 76, 38633881, https://doi.org/10.1175/JAS-D-19-0121.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Adames, Á. F., S. W. Powell, F. Ahmed, V. C. Mayta, and J. D. Neelin, 2021: Tropical precipitation evolution in a buoyancy-budget framework. J. Atmos. Sci., 78, 509528, https://doi.org/10.1175/JAS-D-20-0074.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ahmed, F., and J. D. Neelin, 2018: Reverse engineering the tropical precipitation–buoyancy relationship. J. Atmos. Sci., 75, 15871608, https://doi.org/10.1175/JAS-D-17-0333.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ahmed, F., J. D. Neelin, and A. F. Adames, 2021: Quasi-equilibrium and weak temperature gradient balances in an equatorial beta-plane model. J. Atmos. Sci., 78, 209227, https://doi.org/10.1175/JAS-D-20-0184.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Amador, J. A., 1998: A climatic feature of the tropical Americas: The trade wind easterly jet. Top. Meteor. Oceanogr., 5, 91102.

  • Andersen, J. A., and Z. Kuang, 2012: Moist static energy budget of MJO-like disturbances in the atmosphere of a zonally symmetric aquaplanet. J. Climate, 25, 27822804, https://doi.org/10.1175/JCLI-D-11-00168.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Andrews, D. G., and M. E. McIntyre, 1976: Planetary waves in horizontal and vertical shear: The generalized Eliassen–Palm relation and the mean zonal acceleration. J. Atmos. Sci., 33, 20312048, https://doi.org/10.1175/1520-0469(1976)033<2031:PWIHAV>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Berry, G. J., and C. D. Thorncroft, 2012: African easterly wave dynamics in a mesoscale numerical model: The upscale role of convection. J. Atmos. Sci., 69, 12671283, https://doi.org/10.1175/JAS-D-11-099.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Betts, A. K., 1974: Further comments on “A comparison of the equivalent potential temperature and the static energy.” J. Atmos. Sci., 31, 17131715, https://doi.org/10.1175/1520-0469(1974)031<1713:FCOCOT>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Blake, E. S., and T. B. Kimberlain, 2013: Eastern North Pacific hurricane season of 2011. Mon. Wea. Rev., 141, 13971412, https://doi.org/10.1175/MWR-D-12-00192.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Bui, H. X., and E. D. Maloney, 2019: Mechanisms for global warming impacts on Madden–Julian Oscillation precipitation amplitude. J. Climate, 32, 69616975, https://doi.org/10.1175/JCLI-D-19-0051.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Cao, Z., and H.-R. Cho, 1995: Generation of moist potential vorticity in extratropical cyclones. J. Atmos. Sci., 52, 32633282, https://doi.org/10.1175/1520-0469(1995)052<3263:GOMPVI>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Cheng, Y.-M., C. D. Throncroft, and G. N. Kiladis, 2019: Two contrasting African easterly wave behaviors. J. Atmos. Sci., 76, 17531768, https://doi.org/10.1175/JAS-D-18-0300.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Chikira, M., 2014: Eastward-propagating intraseasonal oscillation represented by Chikira–Sugiyama cumulus parameterization. Part II: Understanding moisture variation under weak temperature gradient balance. J. Atmos. Sci., 71, 615639, https://doi.org/10.1175/JAS-D-13-038.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Clark, J. H. E., R. P. James, and R. H. Grumm, 2002: A reexamination of the mechanisms responsible for banded precipitation. Mon. Wea. Rev., 130, 30743086, https://doi.org/10.1175/1520-0493(2002)130<3074:AROTMR>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Clark, S. K., Y. Ming, and Á. F. Adames, 2020: Monsoon low pressure system–like variability in an idealized moist model. J. Climate, 33, 20512074, https://doi.org/10.1175/JCLI-D-19-0289.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Diaz, M., and W. R. Boos, 2019: Barotropic growth of monsoon depressions. Quart. J. Roy. Meteor. Soc., 145, 824844, https://doi.org/10.1002/qj.3467.

  • Diaz, M., and W. R. Boos, 2021: Evolution of idealized vortices in monsoon-like shears: Application to monsoon depressions. J. Atmos. Sci., 78, 12071225, https://doi.org/10.1175/JAS-D-20-0286.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Dominguez, C., J. M. Done, and C. L. Bruyère, 2020: Easterly wave contributions to seasonal rainfall over the tropical Americas in observations and a regional climate model. Climate Dyn., 54, 191209, https://doi.org/10.1007/s00382-019-04996-7.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Duchon, C. E., 1979: Lanczos filtering in one and two dimensions. J. Appl. Meteor., 18, 10161022, https://doi.org/10.1175/1520-0450(1979)018<1016:LFIOAT>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Dwyer, J. G., and P. A. O’Gorman, 2017: Moist formulations of the Eliassen–Palm flux and their connection to the surface westerlies. J. Atmos. Sci., 74, 513530, https://doi.org/10.1175/JAS-D-16-0111.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Emanuel, K. A., 1983: The Lagrangian parcel dynamics of moist symmetric instability. J. Atmos. Sci., 40, 23682376, https://doi.org/10.1175/1520-0469(1983)040<2368:TLPDOM>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Emanuel, K. A., 1986: An air-sea interaction theory for tropical cyclones. Part I: Steady-state maintenance. J. Atmos. Sci., 43, 585605, https://doi.org/10.1175/1520-0469(1986)043<0585:AASITF>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 15591573, https://doi.org/10.1175/1520-0469(1987)044<1559:BIIAEO>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Feng, T., X.-Q. Yang, J.-Y. Yu, and R. Huang, 2020: Convective coupling in tropical-depression-type waves. Part I: Rainfall characteristics and moisture structure. J. Atmos. Sci., 77, 34073422, https://doi.org/10.1175/JAS-D-19-0172.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fuchs-Stone, Ž., D. J. Raymond, and S. Sentic, 2019: A simple model of convectively coupled equatorial Rossby waves. J. Adv. Model. Earth Syst., 11, 173184, https://doi.org/10.1029/2018MS001433.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gall, J. S., W. M. Frank, and M. C. Wheeler, 2010: The role of equatorial Rossby waves in tropical cyclogenesis. Part I: Idealized numerical simulations in an initially quiescent background environment. Mon. Wea. Rev., 138, 13681382, https://doi.org/10.1175/2009MWR3114.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gonzalez, A. O., and X. Jiang, 2019: Distinct propagation characteristics of intraseasonal variability over the tropical west Pacific. J. Geophys. Res. Atmos., 124, 53325351, https://doi.org/10.1029/2018JD029884.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Grabowski, W. W., and M. W. Moncrieff, 2004: Moisture–convection feedback in the tropics. Quart. J. Roy. Meteor. Soc., 130, 30813104, https://doi.org/10.1256/qj.03.135.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Haynes, P. H., and M. E. McIntyre, 1987: On the evolution of vorticity and potential vorticity in the presence of diabatic heating and frictional or other forces. J. Atmos. Sci., 44, 828841, https://doi.org/10.1175/1520-0469(1987)044<0828:OTEOVA>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Haynes, P. H., and M. McIntyre, 1990: On the conservation and impermeability theorems for potential vorticity. J. Atmos. Sci., 47, 20212031, https://doi.org/10.1175/1520-0469(1990)047<2021:OTCAIT>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hersbach, H., and Coauthors, 2019: Global reanalysis: Goodbye ERA-Interim, hello ERA5. ECMWF Newsletter, No. 147, ECMWF, Reading, United Kingdom, 1724, https://doi.org/10.21957/vf291hehd7.

    • Search Google Scholar
    • Export Citation
  • Holton, J. R., and G. J. Hakim, 2012: An Introduction to Dynamic Meteorology. Academic Press, 552 pp.

  • Hunt, K. M. R., and J. K. Fletcher, 2019: The relationship between Indian monsoon rainfall and low-pressure systems. Climate Dyn., 53, 18591871, https://doi.org/10.1007/s00382-019-04744-x.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hurley, J. V., and W. R. Boos, 2015: A global climatology of monsoon low-pressure systems. Quart. J. Roy. Meteor. Soc., 141, 10491064, https://doi.org/10.1002/qj.2447.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Inoue, K., Á. F. Adames, and K. Yasunaga, 2020: Vertical velocity profiles in convectively coupled equatorial waves and MJO: New diagnoses of vertical velocity profiles in the wavenumber–frequency domain. J. Atmos. Sci., 77, 21392162, https://doi.org/10.1175/JAS-D-19-0209.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Janiga, M. A., and C. D. Thorncroft, 2013: Regional differences in the kinematic and thermodynamic structure of African easterly waves. Quart. J. Roy. Meteor. Soc., 139, 15981614, https://doi.org/10.1002/qj.2047.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Janiga, M. A., and C. Zhang, 2016: MJO moisture budget during dynamo in a cloud-resolving model. J. Atmos. Sci., 73, 22572278, https://doi.org/10.1175/JAS-D-14-0379.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Kiladis, G. N., C. D. Thorncroft, and N. M. J. Hall, 2006: Three-dimensional structure and dynamics of African easterly waves. Part I: Observations. J. Atmos. Sci., 63, 22122230, https://doi.org/10.1175/JAS3741.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Kiladis, G. N., M. C. Wheeler, P. T. Haertel, K. H. Straub, and P. E. Roundy, 2009: Convectively coupled equatorial waves. Rev. Geophys., 47, RG2003, 142, https://doi.org/10.1029/2008RG000266.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lapeyre, G., and I. M. Held, 2004: The role of moisture in the dynamics and energetics of turbulent baroclinic eddies. J. Atmos. Sci., 61, 16931710, https://doi.org/10.1175/1520-0469(2004)061<1693:TROMIT>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Leroux, S., N. M. Hall, and G. N. Kiladis, 2010: A climatological study of transient–mean-flow interactions over West Africa. Quart. J. Roy. Meteor. Soc., 136 (Suppl. 1), 397410, https://doi.org/10.1002/qj.474.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lubis, S. W., and C. Jacobi, 2015: The modulating influence of convectively coupled equatorial waves (CCEWs) on the variability of tropical precipitation. Int. J. Climatol., 35, 14651483, https://doi.org/10.1002/joc.4069.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Madden, R., and F. Robitaille, 1970: A comparison of the equivalent potential temperature and the static energy. J. Atmos. Sci., 27, 327329, https://doi.org/10.1175/1520-0469(1970)027<0327:ACOTEP>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Maloney, E. D., and D. L. Hartmann, 2000: Modulation of eastern North Pacific hurricanes by the Madden–Julian oscillation. J. Climate, 13, 14511460, https://doi.org/10.1175/1520-0442(2000)013<1451:MOENPH>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Maloney, E. D., Á. F. Adames, and H. X. Bui, 2019: Madden–Julian oscillation changes under anthropogenic warming. Nat. Climate Change, 9, 2633, https://doi.org/10.1038/s41558-018-0331-6.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Martin, J. E., J. D. Locatelli, and P. V. Hobbs, 1992: Organization and structure of clouds and precipitation on the mid-Atlantic coast of the United States. Part V: The role of an upper-level front in the generation of a rainband. J. Atmos. Sci., 49, 12931303, https://doi.org/10.1175/1520-0469(1992)049<1293:OASOCA>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • McCann, D. W., 1995: Three-dimensional computations of equivalent potential vorticity. Wea. Forecasting, 10, 798802, https://doi.org/10.1175/1520-0434(1995)010<0798:TDCOEP>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mekonnen, A., C. D. Thorncroft, and A. R. Aiyyer, 2006: Analysis of convection and its association with African easterly waves. J. Climate, 19, 54055421, https://doi.org/10.1175/JCLI3920.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Molinari, J., and D. Vollaro, 2000: Planetary- and synoptic-scale influences on eastern Pacific tropical cyclogenesis. Mon. Wea. Rev., 128, 32963307, https://doi.org/10.1175/1520-0493(2000)128<3296:PASSIO>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Molinari, J., D. Knight, M. Dickinson, D. Vollaro, and S. Skubis, 1997: Potential vorticity, easterly waves, and eastern Pacific tropical cyclogenesis. Mon. Wea. Rev., 125, 26992708, https://doi.org/10.1175/1520-0493(1997)125<2699:PVEWAE>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Monteiro, J. M., and J. Sukhatme, 2016: Quasi-geostrophic dynamics in the presence of moisture gradients. Quart. J. Roy. Meteor. Soc., 142, 187195, https://doi.org/10.1002/qj.2644.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Moore, J. T., and T. E. Lambert, 1993: The use of equivalent potential vorticity to diagnose regions of conditional symmetric instability. Wea. Forecasting, 8, 301308, https://doi.org/10.1175/1520-0434(1993)008<0301:TUOEPV>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Muñoz, E., A. J. Busalacchi, S. Nigam, and A. Ruiz-Barradas, 2008: Winter and summer structure of the Caribbean low-level jet. J. Climate, 21, 12601276, https://doi.org/10.1175/2007JCLI1855.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Murthy, V. S., and W. R. Boos, 2018: Role of surface enthalpy fluxes in idealized simulations of tropical depression spinup. J. Atmos. Sci., 75, 18111831, https://doi.org/10.1175/JAS-D-17-0119.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Neelin, J. D., and I. M. Held, 1987: Modeling tropical convergence based on the moist static energy budget. Mon. Wea. Rev., 115, 312, https://doi.org/10.1175/1520-0493(1987)115<0003:MTCBOT>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Núñez Ocasio, K. M., J. L. Evans, and G. S. Young, 2020: A wave-relative framework analysis of AEW–MCS interactions leading to tropical cyclogenesis. Mon. Wea. Rev., 148, 46574671, https://doi.org/10.1175/MWR-D-20-0152.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Papin, P. P., L. F. Bosart, and R. D. Torn, 2017: A climatology of Central American gyres. Mon. Wea. Rev., 145, 19832000, https://doi.org/10.1175/MWR-D-16-0411.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Peng, K., R. Rotunno, G. H. Bryan, and J. Fang, 2019: Evolution of an axisymmetric tropical cyclone before reaching slantwise moist neutrality. J. Atmos. Sci., 76, 18651884, https://doi.org/10.1175/JAS-D-18-0264.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Plumb, R. A., 1985: On the three-dimensional propagation of stationary waves. J. Atmos. Sci., 42, 217229, https://doi.org/10.1175/1520-0469(1985)042<0217:OTTDPO>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Raymond, D. J., and Ž. Fuchs, 2009: Moisture modes and the Madden–Julian oscillation. J. Climate, 22, 30313046, https://doi.org/10.1175/2008JCLI2739.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Raymond, D. J., S. L. Sessions, and Ž. Fuchs, 2007: A theory for the spinup of tropical depressions. Quart. J. Roy. Meteor. Soc., 133, 17431754, https://doi.org/10.1002/qj.125.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Riehl, H., and J. S. Malkus, 1958: On the heat balance of the equatorial trough zone. Geophysica, 6, 503538.

  • Rotunno, R., and J. Klemp, 1985: On the rotation and propagation of simulated supercell thunderstorms. J. Atmos. Sci., 42, 271292, https://doi.org/10.1175/1520-0469(1985)042<0271:OTRAPO>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Russell, J. O. H., A. Aiyyer, and J. Dylan White, 2020: African easterly wave dynamics in convection-permitting simulations: Rotational stratiform instability as a conceptual model. J. Adv. Model. Earth Syst., 12, e2019MS001706, https://doi.org/10.1029/2019MS001706.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rydbeck, A. V., and E. D. Maloney, 2014: Energetics of east Pacific easterly waves during intraseasonal events. J. Climate, 27, 76037621, https://doi.org/10.1175/JCLI-D-14-00211.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rydbeck, A. V., and E. D. Maloney, 2015: On the convective coupling and moisture organization of east Pacific easterly waves. J. Atmos. Sci., 72, 38503870, https://doi.org/10.1175/JAS-D-15-0056.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rydbeck, A. V., E. D. Maloney, and G. J. Alaka, 2017: In situ initiation of east Pacific easterly waves in a regional model. J. Atmos. Sci., 74, 333351, https://doi.org/10.1175/JAS-D-16-0124.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Schubert, W. H., S. A. Hausman, M. Garcia, K. V. Ooyama, and H.-C. Kuo, 2001: Potential vorticity in a moist atmosphere. J. Atmos. Sci., 58, 31483157, https://doi.org/10.1175/1520-0469(2001)058<3148:PVIAMA>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Serra, Y. L., G. N. Kiladis, and M. F. Cronin, 2008: Horizontal and vertical structure of easterly waves in the Pacific ITCZ. J. Atmos. Sci., 65, 12661284, https://doi.org/10.1175/2007JAS2341.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Serra, Y. L., G. N. Kiladis, and K. I. Hodges, 2010: Tracking and mean structure of easterly waves over the Intra-Americas Sea. J. Climate, 23, 48234840, https://doi.org/10.1175/2010JCLI3223.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Smith, L. M., and S. N. Stechmann, 2017: Precipitating quasigeostrophic equations and potential vorticity inversion with phase changes. J. Atmos. Sci., 74, 32853303, https://doi.org/10.1175/JAS-D-17-0023.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Stechmann, S. N., and A. J. Majda, 2006: