1. Introduction
Since the beginning of the satellite era, observations of the horizontal distribution of tropical cumulus clouds and their time evolution have shown that convective clouds organize into a wide variety of spatial and temporal scales (Chang 1970; Yanai and Murakami 1970b,a; Wallace and Chang 1972; Reed and Recker 1971). Examples of organized convection include mesoscale convective systems (Houze 2004), convectively coupled equatorial waves (Kiladis et al. 2009), and the Madden–Julian oscillation (MJO) (Madden and Julian 1971, 1972). Understanding the mechanisms through which convective clouds organize at different scales and propagate with respect to the mean flow has been a challenging task due to the complex interactions among convection, wind, moisture, and radiation in the tropics, where in situ observations are sparse.
Before satellite-based observations of tropical convection were made, Matsuno (1966) obtained a set of wave solutions to a system of “dry” shallow-water equations on an equatorial beta plane linearized with respect to a resting basic state. The waves identified by Matsuno include the equatorial Rossby wave, the equatorial Kelvin wave, the mixed Rossby–gravity wave, and the eastward- and westward-propagating inertio-gravity waves.
Spectral analysis of satellite brightness temperature data (Takayabu 1994; Wheeler and Kiladis 1999; Wheeler et al. 2000) showed that a majority of the peaks in the spectral power over the wavenumber–frequency space follow the dispersion curves of the waves in the shallow-water system, albeit with a smaller “equivalent depth” than that predicted by purely dry theory. Numerous studies have accounted for this smaller depth in terms of latent heating that is incorporated diagnostically by way of the large-scale vertical velocity. Prescribing convection this way simply slows down the waves by lowering the effective static stability without affecting the wave characteristics in any other ways [see Eq. (19) of Kiladis et al. (2009)]. While this description of convection is simple, the framework does not lead to solutions where a fastest growing mode can be identified. Therefore, it is incomplete since it does not explain the observed horizontal scale of the waves. This shortcoming has led many authors to focus on moist thermodynamics as a mechanism of convective coupling and instability in these waves as well as a means of describing the scale preference (Mapes 2000; Majda et al. 2004; Raymond and Fuchs 2007; Kuang 2008).
Recently, there has been considerable effort to reveal the role of moisture as a dynamically active variable in the tropical atmosphere, especially after a tight coupling between tropical convection and environmental moisture was reported from observations (Sobel 2002; Bretherton et al. 2004; Holloway and Neelin 2009; Rushley et al. 2018). Including moisture as a prognostic variable in simple models of the tropical atmosphere yields new modes that do not exist in a dry system (Fuchs and Raymond 2002, 2005, 2017; Sugiyama 2009; Raymond and Fuchs 2009; Sobel and Maloney 2012, 2013; Adames and Kim 2016), which are collectively called “moisture modes.” The existence of moisture modes suggests that the dry assumption used in the Matsuno model might not describe all of the convectively coupled waves that exist in the tropical atmosphere.
Observational and modeling studies indicate that moisture modes may exist, and the MJO may be a physical manifestation of one (Andersen and Kuang 2012; Pritchard and Bretherton 2014; Sobel et al. 2014; Benedict et al. 2014; Wolding et al. 2016; Jiang et al. 2018). The MJO’s spectral signature exhibits a high coherence between water vapor and precipitation [see Figs. 1 and 4 in Yasunaga and Mapes (2012)]. Furthermore, the central role that water vapor plays in the convective organization of the MJO has been documented by numerous studies (Grabowski and Moncrieff 2004; Del Genio et al. 2012; Andersen and Kuang 2012; Shi et al. 2018).
Recent theoretical models indicate that moisture may also play a central role in the thermodynamics of equatorial Rossby waves and tropical depression (TD)-like disturbances such as some easterly waves and monsoon low pressure systems (Fuchs-Stone et al. 2019; Gonzalez and Jiang 2019; Rydbeck and Maloney 2015; Adames and Ming 2018). These tropical motion systems are also characterized by a balanced wind field and slow propagation (Raymond et al. 2015). Along with the MJO, equatorial Rossby waves and TD-like disturbances can be considered to comprise the moisture mode group of tropical motion systems.
Beside the moisture mode group, another group of tropical motion systems exists. These exhibit fast propagation speeds of ~15 m s−1 or higher. Temperature (buoyancy) fluctuations are considered central to their thermodynamics (Mapes 2000; Mapes et al. 2006; Khouider and Majda 2006; Raymond and Fuchs 2009; Raymond et al. 2015; Herman et al. 2016; Khouider 2019). This category includes Kelvin waves, mixed Rossby–gravity waves and inertio-gravity waves. Collectively, these waves will be referred to as the gravity wave group. This categorical view of convectively coupled tropical motion systems is shown schematically in Fig. 1a.

Two ways in which convectively coupled tropical motion systems can be organized. (a) In a category-based view, the thermodynamics of tropical motion systems correspond to either the moisture mode group or the gravity wave group. (b) In a continuum-based view, the two categories described in (a) form the ends of a spectrum. All tropical motion systems fall somewhere in this spectrum. The nondimensional scale described in this study, that is, Nmode, specifies where in the spectrum a phenomenon may be located.
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1

Two ways in which convectively coupled tropical motion systems can be organized. (a) In a category-based view, the thermodynamics of tropical motion systems correspond to either the moisture mode group or the gravity wave group. (b) In a continuum-based view, the two categories described in (a) form the ends of a spectrum. All tropical motion systems fall somewhere in this spectrum. The nondimensional scale described in this study, that is, Nmode, specifies where in the spectrum a phenomenon may be located.
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
Two ways in which convectively coupled tropical motion systems can be organized. (a) In a category-based view, the thermodynamics of tropical motion systems correspond to either the moisture mode group or the gravity wave group. (b) In a continuum-based view, the two categories described in (a) form the ends of a spectrum. All tropical motion systems fall somewhere in this spectrum. The nondimensional scale described in this study, that is, Nmode, specifies where in the spectrum a phenomenon may be located.
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
While the separation of waves into two categories may provide a useful framework to analyze these waves, several studies indicate that the boundary between the two groups is not as distinct as previously thought. For example, Roundy (2012a) found that the horizontal and vertical structure of convectively coupled Kelvin waves becomes more similar to that of the MJO as its equivalent depth becomes smaller. Another study by Roundy (2012b) showed that the spectral signatures of Kelvin waves and the MJO are continuous over regions of mean westerly winds. These two studies indicate that the distinction between the MJO and Kelvin waves may not be clear. A separate study by Sobel and Kim (2012) analyzed the MJO’s circulation features as they decouple from convection in the equatorial central Pacific. They found that a transition between the slowly propagating MJO convection and the faster planetary-scale Kelvin wave may exist during this decoupling process. A recent study by Powell (2017) indicates that MJO convection over the Indo-Pacific warm pool can be interpreted as a moisture mode. Away from the warm pool, the propagation of convection occurs at a phase speed that is consistent with a Kelvin wave that is weakly coupled to convection. During the transition period that was described by Sobel and Kim (2012), Powell (2017) suggests that the MJO may exhibit characteristics that are a blend of a moisture mode and a Kelvin wave.
In this study, we will show that the moist enthalpy budget offers a unified view of the two groups of tropical motion systems. A nondimensional scale defined as Nmode describes the conditions that lead to moisture mode behavior versus gravity wave behavior. Moisture mode behavior predominates when moisture is the dominant contributor to moist enthalpy while gravity wave behavior is dominant when temperature governs the distribution of moist enthalpy. The two groups are bridged by an intermediate “mixed moisture–gravity” wave group: waves whose thermodynamics exhibit properties of both moisture modes and gravity waves. Our results indicate that moisture modes and gravity waves are the bookends of a thermodynamic spectrum of waves, as depicted in Fig. 1b. While the framework presented here may not adequately describe the structure and instability of observed waves, it does offer a framework in which their thermodynamics can be qualitatively understood.
This study is structured as follows. The following section discusses the theoretical framework including the basic equations, scale analysis and linear wave solutions. Section 3 discusses the conditions that lead to gravity wave and moisture mode behavior in observations and reanalysis data. Our results are synthesized and discussed in further detail in section 4.
2. Theory
a. Basic equations
We will perform our analysis on a set of shallow-water basic equations that correspond to a convectively coupled wave with no meridional wind (υ = 0). The coupling with convection is dictated through a prognostic moisture equation under moist convective adjustment (Betts and Miller 1986; Betts 1986). This is a similar framework to that analyzed by Neelin and Yu (1994), Fuchs and Raymond (2002), and Fuchs and Raymond (2017), who analyzed υ = 0 moist waves under moist convective adjustment. This analysis has numerous limitations, since observed Kelvin waves exhibit significant tilts in their fields (Straub and Kiladis 2003; Kiladis et al. 2009; Herman et al. 2016) and other waves, such as the MJO, have a meridional wind field that play an important role in their dynamics (Kim et al. 2014; Adames et al. 2016). However, we posit that the scaling obtained from this system of equations is also valid for waves that exhibit more complex features, including the Kelvin wave and the MJO, as shown in appendixes B–D. Our choice of focusing our discussion on the υ = 0 wave stems from its simple wave solutions, which will be used to support the scale analysis.
Mean state variables used in this study. In the units column a long dash indicates nondimensional.


b. Scale analysis
Important insight into the dynamics of moist equatorial waves can be obtained by considering a scale analysis for motions that exhibit no meridional wind variation (υ = 0). While this analysis is heavily idealized, appendix B shows that the scaling shown here is the same for most waves in which υ ≠ 0. The main scaling parameters we obtain from this analysis are summarized in Table 2. The scaling presented in this section differs from previous studies (Charney 1963; Yano and Bonazzola 2009) in that the time scale of motion is determined by wave propagation (i.e., a wave time scale) rather than by advection by the mean flow (i.e., an advective time scale). Since we are dealing with equatorially trapped motions, the zonal and meridional length scales are assumed to be different.
Description of the meaning of the nondimensional numbers described in section 2 and what their large and small limits physically describe. The Rossby number (Ro) is included for reference. For the Froude number, gravity waves are the fastest waves that can be obtained from the system of equations shown in Eq. (6), and thus Fr cannot be larger than unity.


Equation (14) indicates that WTG balance is achieved in waves that are substantially slower than dry gravity waves. More quantitatively, one can show that any phase speed less than 15 m s−1 satisfies WTG balance since it will yield Fr2 < 0.1. The equation also indicates that in a dry atmosphere,
c. Linear wave solutions
The dispersion is composed of a growing mode and two damped modes. The solution to Eq. (22) that corresponds to a growing mode is shown in Fig. 2 for different values of τc and

(left) Frequency and (right) growth rate for the growing solution in Eq. (22) (solid) and the approximate solution in Eq. (23a) (dashed). The solution for
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1

(left) Frequency and (right) growth rate for the growing solution in Eq. (22) (solid) and the approximate solution in Eq. (23a) (dashed). The solution for
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
(left) Frequency and (right) growth rate for the growing solution in Eq. (22) (solid) and the approximate solution in Eq. (23a) (dashed). The solution for
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
In the next subsections, we will show that for large
1) Convectively coupled gravity wave
2) υ = 0 moisture mode
3. Insights from observations and reanalysis
a. Data and methods
In this section, we will apply the scaling in section 2b to observational and reanalysis data. We make use of the daily mean column-integrated water vapor ⟨q⟩ in the ERA-Interim dataset (Dee et al. 2011). The data has a 1.5° longitude × 1.5° latitude horizontal resolution. Our precipitation measurements come from the Tropical Rainfall Measurement Mission, product 3b42 (Huffman et al. 2007). We have also verified our results with data from the Global Precipitation Climatology Project (Huffman et al. 2001). Finally, daily outgoing longwave radiation (OLR) observations are taken from the dataset compiled by Liebmann and Smith (1996). All fields shown in this study are anomalies obtained by using a 100-day high-pass Lanczos filter (Duchon 1979).
To determine the properties of observed convectively coupled equatorial waves we make use of space–time spectral analysis (Wheeler and Kiladis 1999; Hendon and Wheeler 2008). To extract the signal of convectively coupled equatorial waves, the time series of OLR is divided into 128-day segments that overlap by 64 days. The segments are tapered to zero through the use of a Hanning window. After tapering, complex fast Fourier transforms (FFTs) are computed in longitude and then in time. Finally, the power spectrum is averaged over all segments and over the 10°N–10°S latitude belt. The number of degrees of freedom is calculated to be ~136 [2 (amplitude and phase) × 34 (years) × 365 (days in a year)/128 (segment length)]. We calculate the signal as Sxx = (Pxx − Pred)/Pxx, where Pred is the red spectrum, following the method of (Masunaga 2007), and a value of 0.2 is considered to be statistically significant in this study. The results shown herein are not sensitive to the choice of window as long as it is of 96 days or longer, which is sufficiently long to capture MJO variability in each segment of time.

(a) Convective adjustment time scale τc estimated using Eq. (26). (b) Nc calculated using the values of τc from (a). (c) Nmode estimated using Eq. (19) and utilizing the values of τc obtained from (a). (d) As in (c), but using the simplified equation for Nmode [Eq. (27)]. In all panels solid lines correspond to the dispersion curves of Kelvin waves, equatorial Rossby waves, and n = 1 inertio-gravity waves with equivalent depths of 12 and 90 m. The dashed box encloses eastward-propagating wavenumbers 1–10 and time scales of 20–100 days, which correspond to the traditional MJO band. The dashed lines in westward-propagating zonal wavenumbers correspond to constant phase speeds of −7 and −11 m s−1, which delineate the spectral region where TD-like disturbances are found (Yasunaga and Mapes 2012). The frequency is in units of cycles per day and the zonal wavenumber is expressed as cycles around the circumference of the equator.
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1

(a) Convective adjustment time scale τc estimated using Eq. (26). (b) Nc calculated using the values of τc from (a). (c) Nmode estimated using Eq. (19) and utilizing the values of τc obtained from (a). (d) As in (c), but using the simplified equation for Nmode [Eq. (27)]. In all panels solid lines correspond to the dispersion curves of Kelvin waves, equatorial Rossby waves, and n = 1 inertio-gravity waves with equivalent depths of 12 and 90 m. The dashed box encloses eastward-propagating wavenumbers 1–10 and time scales of 20–100 days, which correspond to the traditional MJO band. The dashed lines in westward-propagating zonal wavenumbers correspond to constant phase speeds of −7 and −11 m s−1, which delineate the spectral region where TD-like disturbances are found (Yasunaga and Mapes 2012). The frequency is in units of cycles per day and the zonal wavenumber is expressed as cycles around the circumference of the equator.
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
(a) Convective adjustment time scale τc estimated using Eq. (26). (b) Nc calculated using the values of τc from (a). (c) Nmode estimated using Eq. (19) and utilizing the values of τc obtained from (a). (d) As in (c), but using the simplified equation for Nmode [Eq. (27)]. In all panels solid lines correspond to the dispersion curves of Kelvin waves, equatorial Rossby waves, and n = 1 inertio-gravity waves with equivalent depths of 12 and 90 m. The dashed box encloses eastward-propagating wavenumbers 1–10 and time scales of 20–100 days, which correspond to the traditional MJO band. The dashed lines in westward-propagating zonal wavenumbers correspond to constant phase speeds of −7 and −11 m s−1, which delineate the spectral region where TD-like disturbances are found (Yasunaga and Mapes 2012). The frequency is in units of cycles per day and the zonal wavenumber is expressed as cycles around the circumference of the equator.
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
b. Estimation of Nmode from data
Figure 3a shows the wavenumber–frequency distribution of τc, estimated through Eq. (26). Even though τc exhibits a distribution reminiscent to a red spectrum, it is clear that τc is smallest over the Kelvin wave band, with τc values of 3 h or less. It is largest over the band corresponding to large-scale westward-propagating disturbances, where it is on the order of a day. Values corresponding to ~17 h are observed over the MJO band. It is worth noting that the spectrum of τc broadly resembles the cospectrum between precipitation and moist static energy shown in Fig. 2a of Yasunaga et al. (2019), revealing the dominance of moisture at low frequencies. The only significant difference between Fig. 3a and Fig. 2a of Yasunaga et al. (2019) lies in the local maximum in the cospectrum of P′ and MSE found along the dispersion curves of Kelvin waves, which further indicates the nonnegligible role that dry static energy plays in those waves.
Figure 3b shows the wavenumber–frequency distribution of Nmode. The contours approximately follow lines of constant slope, indicating that the phase speed of a wave is its dominant contributor. As a result, Nmode is largest at the highest frequencies and horizontal scales, which corresponds to the fastest waves (i.e., cp = ω/k). The lowest values of Nmode are located at the lowest frequencies and smaller spatial scales, which corresponds to the slowest waves.
The dispersion curves for inertio-gravity are well embedded in the region where Nmode ≫ 1, indicating that the thermodynamics of these waves are predominantly driven by thermal fluctuations, as indicated by previous studies (Mapes et al. 2006; Raymond et al. 2007; Raymond and Fuchs 2009). Although not shown, the spectral signature of mixed Rossby–gravity waves is also well within the gravity wave region of the spectrum (Nmode ≫ 1). The Kelvin wave band lies mostly along the gravity wave end of the spectrum, although smaller values of Nmode of roughly 3 are also observed. This indicates that gravity wave dynamics are, on average, dominant in Kelvin waves, although the role of moisture may not be negligible, as posited by Kuang (2008).
Along the MJO band (dashed box in Fig. 3) Nmode exhibits values ~0.1 near zonal wavenumbers 3–10 and time scales longer than 50 days. Larger values are observed near zonal wavenumber 1 and a time scale of 20 days. This result indicates that the largest-scale components of the MJO might not be explained as a moisture mode, but may also exhibit gravity wave properties (Nmode ~ 1). It has been shown by numerous studies that the part of the MJO cycle when the anomalous circulation travels through the Western Hemisphere can be understood as a circumnavigating Kelvin wave (Matthews 2008; Haertel et al. 2015; Powell 2017).
That the wavenumber–frequency distribution of Nmode in Fig. 3b approximately follows lines of constant phase speed indicates that Nc contributes little to Nmode [see Eq. (19)]. We can further investigate this by analyzing the wavenumber–frequency distribution of Nc, shown in Fig. 3c. Indeed, there is less variation in Nc than in Nmode, although this distribution is not constant. The region corresponding to westward-propagating TD-like disturbances exhibits the largest value of Nc of ~0.1 while spectral regions corresponding to the highest and lowest frequencies exhibit values of ~0.02. Similar values of ~0.03 are observed at most other wavenumbers and frequencies. When Nc is averaged across all wavenumbers and frequencies, a mean value of
We can obtain a simplified spectrum of Nmode using Eq. (21) but using

The approximate phase speed of a wave as a function of Nmode as indicated by Eq. (27). A constant value of
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1

The approximate phase speed of a wave as a function of Nmode as indicated by Eq. (27). A constant value of
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
The approximate phase speed of a wave as a function of Nmode as indicated by Eq. (27). A constant value of
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
To synthesize our results, we can use Eq. (27) to modify the wavenumber–frequency spectrum (Fig. 5a) to show it in terms of the simplified Nmode. This alternate way of showing the signal strength is shown in Fig. 5b. By showing the spectrum this way, the gravity wave group is located near the top of the figure, while the moisture mode group is located near the bottom. As discussed previously, Kelvin waves have among the largest values of Nmode in the symmetric spectrum. Westward-propagating TD-type waves exhibit Nmode values of near unity. Equatorial Rossby waves exhibit Nmode values of ~0.2, although the analysis shown here may not be fully applicable to these waves (see appendix B). The MJO signal shows the most range in its Nmode value. The MJO-related signal at zonal wavenumber 1 is centered near the Nmode = 1. The higher wavenumber components of the MJO fall closer to the moisture mode limit of the spectrum.

(a) Signal strength of equatorially symmetric OLR overlaid by the dispersion curves shown in Fig. 3. (b) As in (a), but using the simplified Nmode from Eq. (27) as the y axis. The limits for the moisture mode and gravity wave are shown as dashed horizontal gray lines. Dotted horizontal lines correspond to Nmode values of 0.2, 1, and 5.
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1

(a) Signal strength of equatorially symmetric OLR overlaid by the dispersion curves shown in Fig. 3. (b) As in (a), but using the simplified Nmode from Eq. (27) as the y axis. The limits for the moisture mode and gravity wave are shown as dashed horizontal gray lines. Dotted horizontal lines correspond to Nmode values of 0.2, 1, and 5.
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
(a) Signal strength of equatorially symmetric OLR overlaid by the dispersion curves shown in Fig. 3. (b) As in (a), but using the simplified Nmode from Eq. (27) as the y axis. The limits for the moisture mode and gravity wave are shown as dashed horizontal gray lines. Dotted horizontal lines correspond to Nmode values of 0.2, 1, and 5.
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
4. Discussion and conclusions
The purpose of this study is to analyze the processes that lead to two seemingly distinct groups of convectively coupled waves. These are waves in which their thermodynamics are determined by the evolution of moisture (moisture modes), and waves in which their thermodynamics are determined by temperature (gravity waves) (see Fig. 1a). Our study is motivated by the suggestion that these families of waves may not be distinct, but comprise part of a spectrum of observed waves, as indicated by Fig. 1b and suggested by previous studies (Roundy 2012a,b; Sobel and Kim 2012; Powell 2017). Our study provides simple scaling arguments that elucidate the physical processes that lead to this spectrum of waves, with the two main groups at the ends of the spectrum.
a. Nmode and the spectrum between moisture modes and gravity waves
Whether a convectively coupled wave exhibits moisture mode or gravity wave behavior can be understood in terms of a nondimensional mode number Nmode. A schematic description of how this ratio affects the distribution of h′ is shown in Fig. 6. When τc is very large relative to τg (small Nmode), precipitation is sufficiently long-lasting that gravity waves can adjust the temperature anomalies to satisfy WTG balance. This results in moisture anomalies being dominant and the behavior of the wave will be that of a moisture mode. When Nmode is large, precipitation removes moisture from the column at a faster rate than gravity waves can relax the temperature anomalies from latent heating toward WTG balance. In such a case, moisture is eliminated quickly and temperature dominates h′, resulting in convectively coupled gravity waves. For Nmode ~ 1 the wave will exhibit behavior of both moisture mode and gravity wave, which we refer to as a mixed moisture–gravity (MMG) wave.

Schematic describing the mechanism in which a wave could exhibit either gravity wave or moisture mode behavior. (left) In the case when Nmode ≪ 1 dry gravity waves propagating away from the region of convection homogenize the anomalous temperature (pink shading), resulting in WTG balance. This occurs more quickly than the removal of moisture by convection. As a result, moisture dominates the distribution of h′ and the dynamics are those of a moisture mode. (right) In the case when Nmode ≫ 1 convection removes moisture from the column more quickly than what dry gravity waves can adjust the column to WTG balance. As a result, temperature dominates the distribution of h′ and the resulting dynamics resemble those of a convectively coupled gravity wave.
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1

Schematic describing the mechanism in which a wave could exhibit either gravity wave or moisture mode behavior. (left) In the case when Nmode ≪ 1 dry gravity waves propagating away from the region of convection homogenize the anomalous temperature (pink shading), resulting in WTG balance. This occurs more quickly than the removal of moisture by convection. As a result, moisture dominates the distribution of h′ and the dynamics are those of a moisture mode. (right) In the case when Nmode ≫ 1 convection removes moisture from the column more quickly than what dry gravity waves can adjust the column to WTG balance. As a result, temperature dominates the distribution of h′ and the resulting dynamics resemble those of a convectively coupled gravity wave.
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
Schematic describing the mechanism in which a wave could exhibit either gravity wave or moisture mode behavior. (left) In the case when Nmode ≪ 1 dry gravity waves propagating away from the region of convection homogenize the anomalous temperature (pink shading), resulting in WTG balance. This occurs more quickly than the removal of moisture by convection. As a result, moisture dominates the distribution of h′ and the dynamics are those of a moisture mode. (right) In the case when Nmode ≫ 1 convection removes moisture from the column more quickly than what dry gravity waves can adjust the column to WTG balance. As a result, temperature dominates the distribution of h′ and the resulting dynamics resemble those of a convectively coupled gravity wave.
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
Further examination of Nmode reveals that its temporal and spatial variations are predominantly determined by the wave’s phase speed. Waves whose phase speed is on the order of ~3 m s−1 are likely to be moisture modes while fast waves with phase speeds of ~30 m s−1 are likely to be gravity waves. Waves with intermediate phase speeds will exhibit mixed behavior. These results indicate that the moisture budget better describes the processes that lead to the evolution of slow waves while the internal energy budget describes the evolution of fast waves. This finding is consistent with those of previous studies (Raymond et al. 2007; Raymond and Fuchs 2009; Herman et al. 2016; Mapes et al. 2006). It is worth mentioning that the description of waves using Nmode is purely diagnostic. It can only describe the properties of a wave if its phase speed is already known. It cannot be used to explain why different phenomena exhibit different propagation speeds.
b. Analogy to geostrophic adjustment
The adjustment process that is described by Nmode is analogous to geostrophic adjustment [see section 3.9 in Vallis (2017)]. In geostrophic adjustment, gravity waves induce acceleration in a region of anomalous geopotential that starts at rest. As the gravity waves propagate away from the geopotential anomaly, they become deflected by the Coriolis force, eventually becoming trapped. This trapping results in the temperature anomalies not being fully smoothed out, which results in the wind field adjusting to geostrophic balance. Because it takes a certain amount of time for geostrophic adjustment to occur, not all observed phenomena are in geostrophic balance. Classic scaling arguments for the midlatitudes reveal that whether a phenomenon is able to reach geostrophic balance or not can be understood via the Rossby number (Ro). As described in the previous subsection, Nmode describes an analogous adjustment process. In the same sense that Ro describes the interplay between gravity wave acceleration and the Coriolis force, Nmode describes the interplay between gravity wave adjustment toward WTG balance and convective moisture adjustment.
c. On the spatial and temporal variations of τc
A salient feature of the wavenumber–frequency distribution of τc in Fig. 3a is that τc is larger at the largest scales and lowest frequencies, while exhibiting smaller values at higher frequencies and smaller scales. These are also the spatial and temporal scales in which precipitation and water vapor exhibit the highest coherence (see Fig. 4 in Yasunaga and Mapes 2012). One possible explanation is that precipitation for faster waves is more tightly coupled to fluctuations in CAPE and CIN, as described by Khouider and Majda (2006); Raymond et al. (2007) and Herman et al. (2016), rather than to fluctuations in moisture alone.
Such a modulation by temperature would increase the amplitude of precipitation without changing the magnitude of the anomalies in specific humidity, thus effectively reducing τc. Ongoing and future work will seek to answer this question using alternative parameterizations of precipitation that are not just limited to moisture fluctuations (F. Ahmed 2019, personal communication).
d. Relevance to the MJO cycle
Our results can be used to reinterpret the MJO cycle. Many studies have suggested that a circumnavigating Kelvin wave initiated subsequent MJO events in the Indian Ocean (Knutson and Weickmann 1987; Matthews 2008; Haertel et al. 2015; Powell and Houze 2015). Our results indicate that the circumnavigating signal could be a Kelvin wave that acquires characteristics of a moisture mode as it approaches the Indian Ocean. This transformation may be due to the MJO exhibiting a lower effective GMS value over the warm pool, which results in a smaller Nmode. A similar conclusion was obtained by Powell (2017), who found that a parameter μ, which can be thought of as
We can investigate the results of Powell (2017) within the context of this study. If we use Eq. (27) along with the phase speed values that Powell (2017) found, we obtain that Nmode ~ 0.2 over the warm pool and Nmode ~ 10 over the Western Hemisphere. Thus, our results agree with those of Powell (2017) in that the MJO likely exhibits moisture mode behavior in the warm pool, while exhibiting gravity wave behavior over the Western Hemisphere. Thus, it seems like MJO initiation could be a circumnavigating Kelvin wave that obtains moisture mode characteristics as it reaches the warm pool. Such a hypothesis should be examined in the future.
e. Relevance to modeling studies
Our results may also shed some light as to why some models are able to simulate Kelvin waves but not MJO-like variability, or vice versa. We hypothesize that models that do the former exhibit variability whose GMS is too large, or their convective schemes remove water vapor from the column too quickly (a too-large Nmode) (Benedict et al. 2014). Additionally, it may also explain why when entrainment is increased in some models MJO activity is enhanced (Del Genio et al. 2012; Kim et al. 2012; Klingaman and Woolnough 2014; Hannah and Maloney 2014; Zhu and Hendon 2015). When the convective scheme is modified to inhibit convection, Nmode is reduced, which increases the likelihood of simulating moisture mode variability in these models.
f. Caveats and future directions
While the results presented in this study may provide useful insights of the thermodynamic properties of moist waves, there are nonetheless multiple caveats. While some of these have already been discussed, we summarize them here for convenience:
A single vertical structure composed of a first baroclinic mode is assumed. Observed waves usually exhibit multiple vertical structures and significant vertical tilts. While this does not necessarily invalidate the scale analysis (section 2b), it does limit the interpretation of the wave solutions (section 2c).
Precipitation is parameterized only as a function of column-integrated water vapor. Rainfall is also dependent on other parameters such as CAPE and CIN.
While appendix B shows that the scaling can be extended to waves in which υ ≠ 0, the corresponding wave solutions were not examined in detail.
The analysis is performed in a linear framework in a precipitating atmosphere. Observed phenomena often exhibit significant nonlinearities and the assumption of a precipitating atmosphere may not always be adequate.
Future work should seek to address these caveats. For example, it is possible to extend the analysis performed here to a framework where no assumptions about the wave’s vertical structure are made.
Acknowledgments
ÁFA was supported by the University of Michigan’s startup package and by the National Science Foundation Grant AGS-1841559. DK was supported by the National Aeronautics and Space Administration Grant 80NSSC17K0227, the U.S. Department of Energy’s Regional and Global Model Analysis program under Grant DE-SC0019495, and Korea Meteorological Administration Research and Development Program under Grant KMI2018-03110. SKC was supported by a National Defense Science and Engineering Graduate Fellowship. KI was supported by an appointment to the NASA Postdoctoral Program at the NASA Goddard Institute for Space Studies, administered by Universities Space Research Association under contract with NASA.
APPENDIX A
Vertical Truncation of the Basic Equations

Basis functions of vertical velocity Ω (solid), horizontal winds and geopotential Λ (dashed), and (c) temperature a (dot–dashed). The value of Ω has been normalized by dividing by
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1

Basis functions of vertical velocity Ω (solid), horizontal winds and geopotential Λ (dashed), and (c) temperature a (dot–dashed). The value of Ω has been normalized by dividing by
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
Basis functions of vertical velocity Ω (solid), horizontal winds and geopotential Λ (dashed), and (c) temperature a (dot–dashed). The value of Ω has been normalized by dividing by
Citation: Journal of the Atmospheric Sciences 76, 12; 10.1175/JAS-D-19-0121.1
APPENDIX B
Nmode Scaling for Motion with υ ≠ 0
a. Scaling for nonnegligible zonal wind acceleration
Scaling of the zonal momentum equation reveals that when Roe is of magnitude 100 or greater, then the scaling for geopotential is [ϕ] ~ Ucp. It follows that the scaling for this case identical to that presented in section 2. Given that equatorially trapped waves exhibit a meridional scale of Ly ~ 106 m, and given that β ~ 10−11 m−1 s−1, it follows that this scaling applies to waves that exhibit a phase speed on the order of 101 m s−1.
b. Scaling when zonal momentum equation satisfies geostrophic balance
APPENDIX C
Nmode Scaling for the MJO
In this section we will show that the Nmode scaling discussed in the main text can also be applied to the MJO, even though the MJO exhibits significant meridional flow (υ ≠ 0). To show this, we invoke Eq. (B2) and perform scale analysis with scaling values that qualitatively correspond to observations of the MJO. Since the MJO is planetary-scale phenomenon, its zonal scale is on the order of Lx ~ 107 m while equatorial trapping indicates that Ly ~ 106 m. Given that the MJO’s phase speed has been estimated to be between 4.7 and 6.6 cp ~ 101 m s−1 (Adames and Kim 2016; Powell 2017), we will use a scale for the phase speed of cp ~ 101 m s−1, which yields Roe ~ 100. As a result, the two terms in square brackets in Eq. (B2) are of comparable magnitude, which indicates that φ scales following Eq. (8). It follows that the resulting definition of Nmode is identical to that shown in the main text. If we used more realistic values of cp and β we would obtain values of 1/Roe that range from 3 to 4.5, which still qualitatively satisfies the scaling in Eq. (8).
APPENDIX D
The Second Baroclinic Mode and the Observed Structure of Convectively Coupled Kelvin Waves
It is important to note that the convectively coupled gravity wave solution obtained in section 2c(1) exhibits only a first baroclinic structure in vertical motion, temperature and wind. Observed convectively coupled Kelvin waves exhibit more complex structures that include a second baroclinic mode (Wheeler et al. 2000; Straub and Kiladis 2003; Kiladis et al. 2009; Herman et al. 2016). While the linear wave solutions discussed do not explain the vertical structure and instability mechanism of observed Kelvin waves, the scale analysis is still be applicable. If we assume that convectively coupled Kelvin waves exhibit a predominantly second baroclinic structure, then c = 27 m s−1. This slower phase speed will cause τg to be nearly twice as long as the case for the first baroclinic mode, resulting in a larger value of Nmode. Waves with phase speeds ~15–20 m s−1, as observed in Kelvin waves exhibit Nmode ~ 101. Thus, second baroclinic convectively coupled Kelvin would correspond to the gravity wave group of waves.
APPENDIX E
Meridional Structure of the υ = 0 Wave Solution
REFERENCES
Adames, Á. F., 2017: Precipitation budget of the Madden–Julian oscillation. J. Atmos. Sci., 74, 1799–1817, https://doi.org/10.1175/JAS-D-16-0242.1.
Adames, Á. F., and J. M. Wallace, 2015: Three-dimensional structure and evolution of the moisture field in the MJO. J. Atmos. Sci., 72, 3733–3754, https://doi.org/10.1175/JAS-D-15-0003.1.
Adames, Á. F., and D. Kim, 2016: The MJO as a dispersive, convectively coupled moisture wave: Theory and observations. J. Atmos. Sci., 73, 913–941, https://doi.org/10.1175/JAS-D-15-0170.1.
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, 2083–2106, https://doi.org/10.1175/JAS-D-17-0310.1.
Adames, Á. F., J. M. Wallace, and J. M. Monteiro, 2016: Seasonality of the structure and propagation characteristics of the MJO. J. Atmos. Sci., 73, 3511–3526, https://doi.org/10.1175/JAS-D-15-0232.1.
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, 2782–2804, https://doi.org/10.1175/JCLI-D-11-00168.1.
Benedict, J. J., and D. A. Randall, 2007: Observed characteristics of the MJO relative to maximum rainfall. J. Atmos. Sci., 64, 2332–2354, https://doi.org/10.1175/JAS3968.1.
Benedict, J. J., E. D. Maloney, A. H. Sobel, and D. M. W. Frierson, 2014: Gross moist stability and MJO simulation skill in three full-physics GCMs. J. Atmos. Sci., 71, 3327–3349, https://doi.org/10.1175/JAS-D-13-0240.1.
Betts, A. K., 1986: A new convective adjustment scheme. Part I: Observational and theoretical basis. Quart. J. Roy. Meteor. Soc., 112, 677–691, https://doi.org/10.1002/QJ.49711247307.
Betts, A. K., and M. Miller, 1986: A new convective adjustment scheme. Part II: Single column tests using GATE wave, BOMEX, ATEX and Arctic air-mass data sets. Quart. J. Roy. Meteor. Soc., 112, 693–709, https://doi.org/10.1002/QJ.49711247308.
Bretherton, C. S., and P. K. Smolarkiewicz, 1989: Gravity waves, compensating subsidence and detrainment around cumulus clouds. J. Atmos. Sci., 46, 740–759, https://doi.org/10.1175/1520-0469(1989)046<0740:GWCSAD>2.0.CO;2.
Bretherton, C. S., M. E. Peters, and L. E. Back, 2004: Relationships between water vapor path and precipitation over the tropical oceans. J. Climate, 17, 1517–1528, https://doi.org/10.1175/1520-0442(2004)017<1517:RBWVPA>2.0.CO;2.
Chang, C.-P., 1970: Westward propagating cloud patterns in the tropical Pacific as seen from time-composite satellite photographs. J. Atmos. Sci., 27, 133–138, https://doi.org/10.1175/1520-0469(1970)027<0133:WPCPIT>2.0.CO;2.
Charney, J. G., 1963: A note on large-scale motions in the tropics. J. Atmos. Sci., 20, 607–609, https://doi.org/10.1175/1520-0469(1963)020<0607:ANOLSM>2.0.CO;2.
Dee, D. P., and Coauthors, 2011: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Quart. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828.
Del Genio, A. D., Y. Chen, D. Kim, and M.-S. Yao, 2012: The MJO transition from shallow to deep convection in CloudSat/CALIPSO data and GISS GCM simulations. J. Climate, 25, 3755–3770, https://doi.org/10.1175/JCLI-D-11-00384.1.
Duchon, C. E., 1979: Lanczos filtering in one and two dimensions. J. Appl. Meteor., 18, 1016–1022, https://doi.org/10.1175/1520-0450(1979)018<1016:LFIOAT>2.0.CO;2.
Emanuel, K. A., 1987: An air–sea interaction model of intraseasonal oscillations in the tropics. J. Atmos. Sci., 44, 2324–2340, https://doi.org/10.1175/1520-0469(1987)044<2324:AASIMO>2.0.CO;2.
Emanuel, K. A., 1994: Atmospheric Convection. Oxford University, 580 pp.
Frierson, D. M. W., 2007: The dynamics of idealized convection schemes and their effect on the zonally averaged tropical circulation. J. Atmos. Sci., 64, 1959–1976, https://doi.org/10.1175/JAS3935.1.
Fuchs, Ž., and D. J. Raymond, 2002: Large-scale modes of a nonrotating atmosphere with water vapor and cloud–radiation feedbacks. J. Atmos. Sci., 59, 1669–1679, https://doi.org/10.1175/1520-0469(2002)059<1669:LSMOAN>2.0.CO;2.
Fuchs, Ž., and D. J. Raymond, 2005: Large-scale modes in a rotating atmosphere with radiative-convective instability and WISHE. J. Atmos. Sci., 62, 4084–4094, https://doi.org/10.1175/JAS3582.1.
Fuchs, Ž., and D. J. Raymond, 2017: A simple model of intraseasonal oscillations. J. Adv. Model. Earth Syst., 9, 1195–1211, https://doi.org/10.1002/2017MS000963.
Fuchs-Stone, Ž., D. J. Raymond, and S. Senti, 2019: A simple model of convectively coupled equatorial Rossby waves. J. Adv. Model. Earth Syst., 11, 173–184, https://doi.org/10.1029/2018MS001433.
Gill, A. E., 1980: Some simple solutions for heat-induced tropical circulation. Quart. J. Roy. Meteor. Soc., 106, 447–462, https://doi.org/10.1002/qj.49710644905.
Gonzalez, A. O., and X. Jiang, 2019: Distinct propagation characteristics of intraseasonal variability over the tropical west Pacific. J. Geophys. Res. Atmos., 124, 5332–5351, https://doi.org/10.1029/2018JD029884.
Grabowski, W. W., and M. W. Moncrieff, 2004: Moisture–convection feedback in the tropics. Quart. J. Roy. Meteor. Soc., 130, 3081–3104, https://doi.org/10.1256/qj.03.135.
Haertel, P., K. Straub, and A. Budsock, 2015: Transforming circumnavigating Kelvin waves that initiate and dissipate the Madden–Julian oscillation. Quart. J. Roy. Meteor. Soc, 141, 1586–1602, https://doi.org/10.1002/QJ.2461.
Hannah, W. M., and E. D. Maloney, 2014: The moist static energy budget in NCAR CAM5 hindcasts during DYNAMO. J. Adv. Model. Earth Syst., 6, 420–440, https://doi.org/10.1002/2013MS000272.
Hayashi, Y., 1979: A generalized method of resolving transient disturbances into standing and traveling waves by space–time spectral analysis. J. Atmos. Sci., 36, 1017–1029, https://doi.org/10.1175/1520-0469(1979)036<1017:AGMORT>2.0.CO;2.
Hendon, H. H., and M. C. Wheeler, 2008: Some space–time spectral analyses of tropical convection and planetary-scale waves. J. Atmos. Sci., 65, 2936–2948, https://doi.org/10.1175/2008JAS2675.1.
Herman, M. J., Z. Fuchs, D. J. Raymond, and P. Bechtold, 2016: Convectively coupled Kelvin waves: From linear theory to global models. J. Atmos. Sci., 73, 407–428, https://doi.org/10.1175/JAS-D-15-0153.1.
Holloway, C. E., and J. D. Neelin, 2009: Moisture vertical structure, column water vapor, and tropical deep convection. J. Atmos. Sci., 66, 1665–1683, https://doi.org/10.1175/2008JAS2806.1.
Houze, R. A., Jr., 2004: Mesoscale convective systems. Rev. Geophys., 42, RG4003, https://doi.org/10.1029/2004RG000150.
Huffman, G. J., R. F. Adler, M. M. Morrissey, D. T. Bolvin, S. Curtis, R. Joyce, B. McGavock, and J. Susskind, 2001: Global precipitation at one-degree daily resolution from multisatellite observations. J. Hydrometeor., 2, 36–50, https://doi.org/10.1175/1525-7541(2001)002<0036:GPAODD>2.0.CO;2.
Huffman, G. J., and Coauthors, 2007: The TRMM Multisatellite Precipitation Analysis (TMPA): Quasi-global, multiyear, combined-sensor precipitation estimates at fine scales. J. Hydrometeor., 8, 38–55, https://doi.org/10.1175/JHM560.1.
Inoue, K., and L. E. Back, 2015: Gross moist stability assessment during TOGA COARE: Various interpretations of gross moist stability. J. Atmos. Sci., 72, 4148–4166, https://doi.org/10.1175/JAS-D-15-0092.1.
Inoue, K., and L. E. Back, 2017: Gross moist stability analysis: Assessment of satellite-based products in the GMS plane. J. Atmos. Sci., 74, 1819–1837, https://doi.org/10.1175/JAS-D-16-0218.1.
Jiang, X., M. Zhao, E. D. Maloney, and D. E. Waliser, 2016: Convective moisture adjustment time scale as a key factor in regulating model amplitude of the Madden-Julian oscillation. Geophys. Res. Lett., 43, 10 412–10 419, https://doi.org/10.1002/2016GL070898.
Jiang, X., A. F. Adames, M. Zhao, D. Waliser, and E. Maloney, 2018: A unified moisture mode framework for seasonality of the Madden–Julian oscillation. J. Climate, 31, 4215–4224, https://doi.org/10.1175/JCLI-D-17-0671.1.
Khouider, B., 2019: Simple models for moist gravity waves. Models for Tropical Climate Dynamics: Waves, Clouds, and Precipitation, Springer, 75–98, https://doi.org/10.1007/978-3-030-17775-1_5.
Khouider, B., and A. J. Majda, 2006: A simple multicloud parameterization for convectively coupled tropical waves. Part I: Linear analysis. J. Atmos. Sci., 63, 1308–1323, https://doi.org/10.1175/JAS3677.1.
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, https://doi.org/10.1029/2008RG000266.
Kim, D., A. H. Sobel, A. D. Del Genio, Y. Chen, S. J. Camargo, M.-S. Yao, M. Kelley, and L. Nazarenko, 2012: The tropical subseasonal variability simulated in the NASA GISS general circulation model. J. Climate, 25, 4641–4659, https://doi.org/10.1175/JCLI-D-11-00447.1.
Kim, D., J.-S. Kug, and A. H. Sobel, 2014: Propagating versus nonpropagating Madden–Julian oscillation events. J. Climate, 27, 111–125, https://doi.org/10.1175/JCLI-D-13-00084.1.
Kim, D., M.-S. Ahn, I.-S. Kang, and A. D. Del Genio, 2015: Role of longwave cloud–radiation feedback in the simulation of the Madden–Julian oscillation. J. Climate, 28, 6979–6994, https://doi.org/10.1175/JCLI-D-14-00767.1.
Klingaman, N. P., and S. J. Woolnough, 2014: Using a case-study approach to improve the Madden–Julian oscillation in the Hadley Centre model. Quart. J. Roy. Meteor. Soc., 140, 2491–2505, https://doi.org/10.1002/qj.2314.
Knutson, T. R., and K. M. Weickmann, 1987: 30–60 day atmospheric oscillations: Composite life cycles of convection and circulation anomalies. Mon. Wea. Rev., 115, 1407–1436, https://doi.org/10.1175/1520-0493(1987)115<1407:DAOCLC>2.0.CO;2.
Kuang, Z., 2008: A moisture-stratiform instability for convectively coupled waves. J. Atmos. Sci., 65, 834–854, https://doi.org/10.1175/2007JAS2444.1.
Liebmann, B., and C. A. Smith, 1996: Description of a complete (interpolated) outgoing longwave radiation dataset. Bull. Amer. Meteor. Soc., 77, 1275–1277, https://doi.org/10.1175/1520-0477-77.6.1274.
Madden, R., and P. Julian, 1971: Detection of a 40–50 day oscillation in the zonal wind in the tropical Pacific. J. Atmos. Sci., 28, 702–708, https://doi.org/10.1175/1520-0469(1971)028<0702:DOADOI>2.0.CO;2.
Madden, R., and P. Julian, 1972: Description of global scale circulation cells in the tropics with a 40–50 day period. J. Atmos. Sci., 29, 1109–1123, https://doi.org/10.1175/1520-0469(1972)029<1109:DOGSCC>2.0.CO;2.
Majda, A. J., B. Khouider, G. N. Kiladis, K. H. Straub, and M. G. Shefter, 2004: A model for convectively coupled tropical waves: Nonlinearity, rotation, and comparison with observations. J. Atmos. Sci., 61, 2188–2205, https://doi.org/10.1175/1520-0469(2004)061<2188:AMFCCT>2.0.CO;2.
Maloney, E. D., 2009: The moist static energy budget of a composite tropical intraseasonal oscillation in a climate model. J. Climate, 22, 711–729, https://doi.org/10.1175/2008JCLI2542.1.
Mapes, B. E., 2000: Convective inhibition, subgrid-scale triggering energy, and stratiform instability in a toy tropical wave model. J. Atmos. Sci., 57, 1515–1535, https://doi.org/10.1175/1520-0469(2000)057<1515:CISSTE>2.0.CO;2.
Mapes, B. E., S. Tulich, J. Lin, and P. Zuidema, 2006: The mesoscale convection life cycle: Building block or prototype for large-scale tropical waves? Dyn. Atmos. Oceans, 42, 3–29, https://doi.org/10.1016/j.dynatmoce.2006.03.003.
Masunaga, H., 2007: Seasonality and regionality of the Madden–Julian oscillation, Kelvin wave and equatorial Rossby wave. J. Atmos. Sci., 64, 4400–4416, https://doi.org/10.1175/2007JAS2179.1.
Matsuno, T., 1966: Quasi-geostrophic motions in the equatorial area. J. Meteor. Soc. Japan, 44, 25–43, https://doi.org/10.2151/jmsj1965.44.1_25.
Matthews, A. J., 2008: Primary and successive events in the Madden-Julian oscillation. Quart. J. Roy. Meteor. Soc., 134, 439–453, https://doi.org/10.1002/qj.224.
Neelin, J. D., and I. M. Held, 1987: Modeling tropical convergence based on the moist static energy budget. Mon. Wea. Rev., 115, 3–12, https://doi.org/10.1175/1520-0493(1987)115<0003:MTCBOT>2.0.CO;2.
Neelin, J. D., and J.-Y. Yu, 1994: Modes of tropical variability under convective adjustment and the Madden–Julian oscillation. Part I: Analytical theory. J. Atmos. Sci., 51, 1876–1894, https://doi.org/10.1175/1520-0469(1994)051<1876:MOTVUC>2.0.CO;2.
Neelin, J. D., and N. Zeng, 2000: A quasi-equilibrium tropical circulation model—Formulation. J. Atmos. Sci., 57, 1741–1766, https://doi.org/10.1175/1520-0469(2000)057<1741:AQETCM>2.0.CO;2.
Neelin, J. D., I. M. Held, and K. H. Cook, 1987: Evaporation–wind feedback and low-frequency variability in the tropical atmosphere. J. Atmos. Sci., 44, 2341–2348, https://doi.org/10.1175/1520-0469(1987)044<2341:EWFALF>2.0.CO;2.
Peters, M. E., and C. S. Bretherton, 2005: A simplified model of the Walker circulation with an interactive ocean mixed layer and cloud-radiative feedbacks. J. Climate, 18, 4216–4234, https://doi.org/10.1175/JCLI3534.1.
Powell, S. W., 2017: Successive MJO propagation in MERRA-2 reanalysis. Geophys. Res. Lett., 44, 5178–5186, https://doi.org/10.1002/2017GL073399.
Powell, S. W., and R. A. Houze, 2015: Effect of dry large-scale vertical motions on initial MJO convective onset. J. Geophys. Res. Atmos., 120, 4783–4805, https://doi.org/10.1002/2014JD022961.
Pritchard, M. S., and C. S. Bretherton, 2014: Causal evidence that rotational moisture advection is critical to the superparameterized Madden–Julian oscillation. J. Atmos. Sci., 71, 800–815, https://doi.org/10.1175/JAS-D-13-0119.1.
Raymond, D. J., and Ž. Fuchs, 2007: Convectively coupled gravity and moisture modes in a simple atmospheric model. Tellus, 59, 627–640, https://doi.org/10.1111/j.1600-0870.2007.00268.x.
Raymond, D. J., and Ž. Fuchs, 2009: Moisture modes and the Madden–Julian oscillation. J. Climate, 22, 3031–3046, https://doi.org/10.1175/2008JCLI2739.1.
Raymond, D. J., S. L. Sessions, and Ž. Fuchs, 2007: A theory for the spinup of tropical depressions. Quart. J. Roy. Meteor. Soc., 133, 1743–1754, https://doi.org/10.1002/QJ.125.
Raymond, D. J., S. L. Sessions, A. H. Sobel, and Ž. Fuchs, 2009: The mechanics of gross moist stability. J. Adv. Model. Earth Syst., 1, 9, https://doi.org/10.3894/JAMES.2009.1.9.
Raymond, D., Ž. Fuchs, S. Gjorgjievska, and S. Sessions, 2015: Balanced dynamics and convection in the tropical troposphere. J. Adv. Model. Earth Syst., 7, 1093–1116, https://doi.org/10.1002/2015MS000467.
Reed, R. J., and E. E. Recker, 1971: Structure and properties of synoptic-scale wave disturbances in the equatorial western Pacific. J. Atmos. Sci., 28, 1117–1133, https://doi.org/10.1175/1520-0469(1971)028<1117:SAPOSS>2.0.CO;2.
Roundy, P. E., 2012a: Observed structure of convectively coupled waves as a function of equivalent depth: Kelvin waves and the Madden–Julian oscillation. J. Atmos. Sci., 69, 2097–2106, https://doi.org/10.1175/JAS-D-12-03.1.
Roundy, P. E., 2012b: The spectrum of convectively coupled Kelvin waves and the Madden–Julian oscillation in regions of low-level easterly and westerly background flow. J. Atmos. Sci., 69, 2107–2111, https://doi.org/10.1175/JAS-D-12-060.1.
Ruppert, J. H., and C. Hohenegger, 2018: Diurnal circulation adjustment and organized deep convection. J. Climate, 31, 4899–4916, https://doi.org/10.1175/JCLI-D-17-0693.1.
Rushley, S. S., D. Kim, C. S. Bretherton, and M.-S. Ahn, 2018: Reexamining the nonlinear moisture-precipitation relationship over the tropical oceans. Geophys. Res. Lett., 45, 1133–1140, https://doi.org/10.1002/2017GL076296.
Rydbeck, A. V., and E. D. Maloney, 2015: On the convective coupling and moisture organization of east Pacific easterly waves. J. Atmos. Sci., 72, 3850–3870, https://doi.org/10.1175/JAS-D-15-0056.1.