1. Introduction
Deep convection plays a key role in the meteorology and climatology of the tropics. It controls the lapse rate of the troposphere and influences the climatologicalmean transport of energy (Hartmann 2015). Furthermore, a multitude of tropical phenomena are tightly coupled to deep convection: the intertropical convergence zone, monsoon systems, El Niño–Southern Oscillation (ENSO), the Madden–Julian oscillation (MJO) (Madden and Julian 1972), convectively coupled waves (Kiladis et al. 2009), and tropical cyclones. These phenomena can have major impacts on tropical and extratropical weather and are important sources of predictability (Kim et al. 2018; Tang et al. 2018). While the importance of tropical deep convection cannot be understated, the processes that lead to its onset, evolution and coupling with the largescale circulation have remained elusive (Kiladis et al. 2009; Kuo et al. 2017; Schiro and Neelin 2019). Furthermore, an accurate representation of convection and its variability in global climate models (GCMs) remains a significant challenge (Taylor et al. 2012; Stanfield et al. 2016).
Scale analysis of the tropical belt (Charney 1963; Yano and Bonazzola 2009) revealed that a simple thermodynamic balance prevails: Diabatic heating is approximately balanced by vertical advection of dry static energy (or similarly potential temperature or dry entropy). Because of the weak Coriolis force, dry gravity waves are very effective at eliminating temperature fluctuations (Bretherton and Smolarkiewicz 1989; Wolding et al. 2016). This adjustment process results in a relatively homogeneous distribution of temperature, where horizontal and temporal fluctuations in freetropospheric temperatures rarely exceed 1 K (Sobel and Bretherton 2000). Thus, to leading order, the tropical troposphere is in weak temperature gradient (WTG) balance (Sobel et al. 2001). The application of WTG balance has led to many advances in our understanding of tropical phenomena such as the MJO, the Walker circulation, the diurnal cycle of convection and tropical cyclogenesis (Bretherton and Sobel 2002; Raymond and Sessions 2007; Chikira 2014; Ruppert and Hohenegger 2018).
While WTG balance is the dominant thermodynamic balance in the tropical free troposphere, it is inapplicable in the planetary boundary layer (PBL) (Sobel and Bretherton 2000). In this layer, energy and momentum exchanges with the surface result in strong turbulent mixing, a process that dominates over the gravity wave adjustment process that leads to WTG balance. Instead, observations indicate that in convecting regions, a balance exists between the energy input from surface fluxes and the energy sink that arises from convective downdrafts and turbulent entrainment (Emanuel 1993; Raymond 1995; ThayerCalder and Randall 2015; de Szoeke 2018). This balance is often referred to as boundary layer quasi equilibrium ^{ 1 } (BLQE; Raymond 1997 and references therein), as is often used as a simple convective parameterization (Yano and Emanuel 1991; Emanuel 2019).
While the strict application of WTG and BLQE have provided numerous insights about the nature and occurrence of deep convection, they are still incomplete treatments of the tropical atmosphere. Observations reveal that temperature fluctuations in the lower free troposphere (LFT), although small, modulate the convective available potential energy (CAPE) and convective inhibition (CIN) in ways that can enhance or suppress deep convection (Mapes 2000; Raymond et al. 2006; Kuang 2008). Furthermore, observations suggest that water vapor fluctuations in the LFT play a central role in the organization of deep convection (Raymond 2000; Grabowski and Moncrieff 2004; Sahany et al. 2012). This regulating role is achieved through dry air entrainment and dilution reducing the buoyancy of rising cumulus clouds (Lucas et al. 1994; Hannah 2017; Kuo et al. 2017). To further complicate matters, large fluctuations of PBL moist enthalpy (ME) are observed over daily and synoptic time scales, which modulate CAPE (Donner and Phillips 2003). Indeed, observed variations of moisture and temperature in the LFT as well as fluctuations of PBL ME coincide with a multitude of atmospheric conditions that enhance or suppress rainfall (Powell 2019).
In spite of these challenges, recent efforts to elucidate the relationship between the atmospheric thermodynamic environment and rainfall have yielded promising results. Recently, Ahmed and Neelin (2018, henceforth AN18) developed a framework aimed at understanding precipitation within the context of an entraining plume model. In their framework, plumes rising out of the PBL entrained environmental air through a mass inflow profile akin to that observed in organized convection (Kingsmill and Houze 1999; Mechem et al. 2002; Schiro et al. 2018). They found that the mean buoyancy of the plume in the LFT was highly correlated with precipitation. Thus, their precipitation–buoyancy relation may serve as the foundation for a conceptual model of tropical rainfall that includes the aforementioned key thermodynamic variables.
In Fig. 1, leadingorder dynamic balances in the atmosphere (left column) are paired with formulations of secondorder balances (right column) in the samecolored boxes that describe the evolution of the quantities listed by samecolored text. For example, in the midlatitudes, the leadingorder balance in the horizontal momentum equations (green) is geostrophic balance. The secondorder balance describes the evolution of the geostrophic wind as formulated in the quasigeostrophic approximation. Similar balances for the tropics can and have been derived in the thermodynamic equations. In this study, the secondorder balances for the freetropospheric thermodynamic and moisture and PBL ME equations are derived from their leadingorder balances. These are then used to derive a buoyancy tendency equation based on the precipitation–buoyancy relation of AN18, which itself can be derived from the secondorder balance in vertical momentum.
This study is structured as follows. Data and methods are described in section 2. Section 3 discusses the buoyancy of an entraining plume and its relation to tropical rainfall. Section 4 derives the leading and secondorder thermodynamic equations used to understand buoyancy evolution. In section 5 we obtain a prognostic equation for plume buoyancy and analyze the processes that lead to its evolution. A concluding discussion is offered in section 6.
2. Data and methods
a. ERA5 and TRMM 3B42
We utilized output from the fifth reanalysis from the European Centre for MediumRange Weather Forecasts (ECMWF) (ERA5; Hersbach et al. 2019). We used instantaneous fields with a time interval of 3 h, spanning the 40yr interval of 1979–2018. The ERA5 data used have horizontal resolution of 0.25° × 0.25° and 14 vertical levels from 1000 to 600 hPa. We make use of the following ERA5 fields: precipitation rate (P), temperature (T), specific humidity (q), and vertical velocity (ω). Only data within the 15°N/S latitude belt are used.
In section 5 we make use of precipitation data from the 3hourly, 0.25° × 0.25° horizontal resolution Tropical Rainfall Measuring Mission (TRMM) 3B42 version 7A (Huffman et al. 2007) dataset. The data are area averaged over the Dynamics of the Madden–Julian Oscillation (DYNAMO) (Yoneyama et al. 2013) northern array domain (0°–5°N, 73°–80°E) in order to create a time series of rainfall that can be compared to the in situ soundingbased data.
b. Soundings
We complemented the reanalysis data by also using sounding data from Majuro (7.1°N, 171.4°E), Manaus (3.2°S, 59.9°W), and Gan (0.6°S, 73.1°E). For the former two, twicedaily observations from a landbased rawinsonde station network for 2000–13 were obtained from the Integrated Global Radiosonde Archive (IGRA; Durre et al. 2006). The rawinsonde dataset was linearly interpolated in pressure to regular 25hPa intervals. The sounding data are used to obtain the mean values of the mean thermodynamic variables and their variance, which are in turn used in the scale analysis (see Table 3).
Furthermore, evaluation of the plume buoyancy equation in section 5 was performed with the rawinsonde measurements taken from the DYNAMO northern sounding array (NSA), located in the central equatorial Indian Ocean (Ciesielski et al. 2014). Stations in the NSA were Gan Island (0.69°N, 73.51°E), the R/V Revelle (0°, 80.5°E), Colombo (6.91°N, 79.878°E), and Malé (4.91°N, 73.53°E). Ciesielski et al. (2014) describes the details of the sounding data, observation characteristics, and qualitycontrol procedures for DYNAMO soundings. A complete description of the thermodynamic budget terms used in section 5 is provided by Johnson et al. (2015).
c. Linear regression
The relationship between LFTaveraged plume buoyancy B _{ L } and terms describing its tendency was determined with DYNAMO soundings using linear regression following the same procedure as many previous studies (e.g., Straub and Kiladis 2002; Sumi and Masunaga 2016; among others). The results presented in section 5 were obtained by linearly regressing 3hourly plume buoyancy equation terms and TRMM rainfall against the B _{ L } anomaly time series. The B _{ L } anomalies were calculated for 10 October–31 December. The linear regression results were then scaled to a 0.025 m s^{−2} anomaly at zero lag, a typical value that is associated with the occurrence of deep convection. Because of the use of linear regression, the fields that are discussed in section 5 are anomalies that correspond to fluctuations with respect to a B _{ L } time series.
3. Plume buoyancy and tropical precipitation
Constants used in this study with their units and values.
Main variables and definitions used in this study.
The measure of plume buoyancy in Eq. (1a) is a variant of the standard potential temperaturebased buoyancy [see Eq. (8.1) in Petty 2008]. It is rewritten so that the plume’s contribution to B is expressed in terms of θ _{ e }, a variable that is approximately conserved during moist adiabatic processes. Equation (1a) defines buoyancy in the same way as AN18, except that it is divided by κ. The inclusion of κ results in buoyancy values that are ~3 times smaller than those reported in AN18. However, the salient results of both papers are nearly unaffected by this factor. We will use the definition of buoyancy from Eq. (1a) since it can be obtained directly from the traditional temperaturebased definition of B, as discussed in appendix A.
Previous studies have shown that water vapor in the LFT as well as fluctuations in CAPE and CIN play an important role in the occurrence and organization of tropical convection (Mapes 2000; Grabowski and Moncrieff 2004; Raymond et al. 2006; Tulich and Mapes 2010; Kuang 2010). Based on results from previous studies, AN18 posited that buoyancy averaged over the LFT (B
_{
L
}) can robustly diagnose precipitation over the tropics. This measure of plume buoyancy combines the known sensitivity of convection to lowertropospheric temperature and moisture into a single variable. AAN simplified the framework by showing that B
_{
L
} can be described in terms of the θ
_{
e
} and
Schematic describing the entraining plume buoyancy described by AN18 and AAN. By applying a deep inflow profile to prescribe entrainment, the mean LFT plume buoyancy B
_{
L
} can be described by Eq. (4) assuming that
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
Schematic describing the entraining plume buoyancy described by AN18 and AAN. By applying a deep inflow profile to prescribe entrainment, the mean LFT plume buoyancy B
_{
L
} can be described by Eq. (4) assuming that
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
Schematic describing the entraining plume buoyancy described by AN18 and AAN. By applying a deep inflow profile to prescribe entrainment, the mean LFT plume buoyancy B
_{
L
} can be described by Eq. (4) assuming that
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
The first term in Eq. (4) is a measure of the convective instability of the lower troposphere (Raymond et al. 2015). It is reminiscent of the deep CIN discussed by Raymond et al. It can also be thought as the undiluted component of B
_{
L
}. The second term represents the reduction in buoyancy that a rising plume experiences as it entrains and is diluted by freetropospheric air. In defining Eq. (4), we are assuming that
Figure 3 shows the distribution of rainfall as a function of B _{ L }. Since B _{ L } contains an integrand spanning both positive and negative buoyancy values, it is impacted by fluctuations in both CAPE and CIN. It is clear that rainfall increases rapidly as B _{ L } increases from negative values closer to zero. Similar rapid increases for precipitation conditioned on combinations of environmental thermodynamic variables—for example, column saturation fraction, column water vapor or buoyancy—have been modeled as an exponential curve (Bretherton et al. 2004; Rushley et al. 2018) or as a ramp function (AAN; Kuo et al. 2018). While we do not explicitly fit a function to the precipitation curve in Fig. 3, we note a critical value of B _{ L } ~ −0.02 m s^{−2} that would mark the beginning of a linear precipitation regime. The factors that govern this critical buoyancy value are not yet clear. However, any missing physics from the B _{ L } formulation when treated as stochastic fluctuations in the value of critical B _{ L } can generate curvature suggestive of an exponential precipitation increase (AAN; Stechmann and Neelin 2011).
Normalized distribution of P and B _{ L } for ERA5 output over the 15°N–15°S latitude belt. Values of P are binned every mm day^{−1} while B _{ L } is binned every 0.0015 m s^{−2}. The circles show the mean value of the distribution.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
Normalized distribution of P and B _{ L } for ERA5 output over the 15°N–15°S latitude belt. Values of P are binned every mm day^{−1} while B _{ L } is binned every 0.0015 m s^{−2}. The circles show the mean value of the distribution.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
Normalized distribution of P and B _{ L } for ERA5 output over the 15°N–15°S latitude belt. Values of P are binned every mm day^{−1} while B _{ L } is binned every 0.0015 m s^{−2}. The circles show the mean value of the distribution.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
While a significant scatter exists in the data points in Fig. 3, the points are more clustered in the P–B _{ L } scatterplot than what it is when it is compared to other measures of rainfall such as column water vapor or saturation fraction (not shown) (AAN; Wolding et al. 2020). Further discussion of the joint PDF of precipitation with buoyancy may be found in Kuo et al. (2018), where a cruder estimator of buoyancy (based on column water vapor relative to a temperaturedependent critical value) is used. The robustness of the P–B _{ L } in Fig. 3 and the fact that B _{ L } accounts for fields that are known to play a key role in the evolution of precipitation such as convective instability and moisture justifies its use as a framework to understand precipitation evolution.
Normalized distribution of γ _{ T } [sensitivity of buoyancy to T _{ L } from Eq. (10)] over the 15°N–15°S latitude belt for ERA5. γ _{ T } is binned at intervals of 0.025. The median value of γ = 4.2 is depicted as a vertical dashed line; 95% of points for ERA5 lie within 10% of the median value, denoted as the vertical dotted lines.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
Normalized distribution of γ _{ T } [sensitivity of buoyancy to T _{ L } from Eq. (10)] over the 15°N–15°S latitude belt for ERA5. γ _{ T } is binned at intervals of 0.025. The median value of γ = 4.2 is depicted as a vertical dashed line; 95% of points for ERA5 lie within 10% of the median value, denoted as the vertical dotted lines.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
Normalized distribution of γ _{ T } [sensitivity of buoyancy to T _{ L } from Eq. (10)] over the 15°N–15°S latitude belt for ERA5. γ _{ T } is binned at intervals of 0.025. The median value of γ = 4.2 is depicted as a vertical dashed line; 95% of points for ERA5 lie within 10% of the median value, denoted as the vertical dotted lines.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
It is noteworthy that Eq. (9) exhibits interesting parallels to the way the temperature and moisture equations are combined to create evolution equations for cloud work function or entraining CAPE in certain convective closures (e.g., Arakawa and Schubert 1974; Moorthi and Suarez 1992; Zhang and McFarlane 1995; Neelin et al. 2008).
4. Largescale moist thermodynamics under WTG balance
a. Temperature evolution in the LFT
Scaling of the equations in section 4 based on estimates of their observed variations. For the scaling, the horizontal winds are assumed to have a scale of 10 m s^{−1}, the horizontal scale is assumed to be 10^{6} m, and the variations of c _{ p } T, Lq, and h in the DBL and LFT are rounded from the variances seen in Fig. B1. The first row is in units of Pa s^{−1} while subsequent rows are in units of J kg^{−1} s^{−1}.
(top) Time series from the DYNAMO northern array of LFTaveraged temperature tendency (blue, left y axis) and LFTaveraged heating Q _{1} (red, right y axis). (bottom) As in the top panel, but showing LFTaveraged ω _{ Q } calculated using Eq. (21) and ω _{ a } calculated using Eq. (23). In both panels, a 1day running mean is used to remove diurnal variability from the time series. Note that the scale for ω _{ a } is an order of magnitude smaller than the scale of ω _{ Q }.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
(top) Time series from the DYNAMO northern array of LFTaveraged temperature tendency (blue, left y axis) and LFTaveraged heating Q _{1} (red, right y axis). (bottom) As in the top panel, but showing LFTaveraged ω _{ Q } calculated using Eq. (21) and ω _{ a } calculated using Eq. (23). In both panels, a 1day running mean is used to remove diurnal variability from the time series. Note that the scale for ω _{ a } is an order of magnitude smaller than the scale of ω _{ Q }.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
(top) Time series from the DYNAMO northern array of LFTaveraged temperature tendency (blue, left y axis) and LFTaveraged heating Q _{1} (red, right y axis). (bottom) As in the top panel, but showing LFTaveraged ω _{ Q } calculated using Eq. (21) and ω _{ a } calculated using Eq. (23). In both panels, a 1day running mean is used to remove diurnal variability from the time series. Note that the scale for ω _{ a } is an order of magnitude smaller than the scale of ω _{ Q }.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
b. Moisture evolution in the LFT
It is worth noting that Eq. (29) can be thought of as a secondorder balance equation since the leadingorder balance in precipitating regions is between largescale vertical moisture advection and Q _{2} (see Fig. 1). It is similar to the moisture budgets described by Chikira (2014), Wolding and Maloney (2015) and Wolding et al. (2016), except that diabatic heating is expressed in terms of vertical velocities, and the contribution from ω _{ a } is included in the equation.
From inspection of Eq. (29) we can see that ω _{ a } and ω _{ r } moisten the atmosphere through vertical moisture advection. Convective processes, however, moisten the atmosphere through vertical MSE advection. While ascent driven by adiabatic processes and radiation usually moistens the troposphere, ascent driven by convection only moistens the atmosphere beneath the minimum in MSE, where the vertical MSE gradient is usually positive (beneath the 650 hPa layer in Fig. B1 in appendix B). Above this level, vertical MSE advection by convection usually dries the troposphere.
c. ME in the free troposphere
Equation (32) has several advantages over the traditional MSE budget. First, the contribution of radiative heating to the MSE tendency is included in a single term. In traditional MSE budgets, there is an implicit contribution of radiative heating to vertical MSE advection since ω = ω _{ c } + ω _{ r } + ω _{ a }. Additionally, it is clear that the contribution of ω _{ a } to the evolution of ME is negligible (see Table 3).
It is worth pointing out that if WTG balance is applied strictly, ∂_{ t } T and ω _{ a } are dropped from Eq. (32), and both Eqs. (32) and (29) become the WTGbased moisture equation used by Chikira (2014), Wolding et al. (2016) and others.
d. Moisture and thermodynamic equations averaged over the LFT
e. ME averaged over the DBL
5. The plume buoyancy equation

A: horizontal advection of moisture and temperature in the LFT,

B: horizontal advection of ME in the BL,

C: moistening and cooling by adiabatic vertical motions in the LFT,

D: vertical moisture advection driven by radiative heating in the LFT,

E: vertical MSE advection by convection in the LFT,

F: vertical advection of MSE at the top of the DBL by adiabatic motions,

G: as in F, but for radiatively driven vertical motions,

H: as in F, but for convection at the top of the DBL,

I: sources of ME in the DBL.
a. Scaling and interpretation
Prior to analyzing the evolution of Eq. (46), it is worth interpreting the processes that can lead to the evolution of B _{ L } in tropical motion systems and how this evolution can lead to convective coupling. A schematic describing some of these processes is shown in Fig. 6. Convection is modulated differently depending on the layer in which the process occurs even though their impact on B _{ L } is qualitatively similar. For example, horizontal moisture advection in the LFT (Fig. 6a) reduces the saturation deficit in Eq. (4). Such a reduction implies that ascending plumes dilute less as they ascend through the troposphere. In the case of horizontal moisture advection in the DBL (Fig. 6b), moistening causes ascending updrafts to be more humid and hence have a larger θ _{ e }. As a result, a plume is able to liberate more latent heat from condensation, implying a larger undilute B _{ L }. Changes in surface heat fluxes have the same impact (Fig. 6c). Cloudradiative heating, which impacts both the DBL and the LFT may increase undilute B _{ L } and reduce the saturation deficit (Fig. 6d).
Schematic describing processes that increase B _{ L } in synopticscale tropical motion systems. The way each process modulates the individual terms that contribute to B _{ L } [Eq. (4)] are shown at the top of each panel. (a)–(d),(f) The blue background shading denotes the mean distribution of moisture. (e) The red shading denotes the distribution of temperature. Note that (a), (b), and (e) are latitude–height cross sections that represent meridional moisture and temperature gradients found in the tropics. The other panels are longitude–height cross sections.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
Schematic describing processes that increase B _{ L } in synopticscale tropical motion systems. The way each process modulates the individual terms that contribute to B _{ L } [Eq. (4)] are shown at the top of each panel. (a)–(d),(f) The blue background shading denotes the mean distribution of moisture. (e) The red shading denotes the distribution of temperature. Note that (a), (b), and (e) are latitude–height cross sections that represent meridional moisture and temperature gradients found in the tropics. The other panels are longitude–height cross sections.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
Schematic describing processes that increase B _{ L } in synopticscale tropical motion systems. The way each process modulates the individual terms that contribute to B _{ L } [Eq. (4)] are shown at the top of each panel. (a)–(d),(f) The blue background shading denotes the mean distribution of moisture. (e) The red shading denotes the distribution of temperature. Note that (a), (b), and (e) are latitude–height cross sections that represent meridional moisture and temperature gradients found in the tropics. The other panels are longitude–height cross sections.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
In synopticscale tropical motion systems, departures from WTG balance can arise from largescale adiabatic lifting. When such processes occur Eq. (46) may be preferable to ME or MSE budgets. As discussed in section 4c and shown in Table 3, adiabatic motions contribute little to the evolution of ME but are of leadingorder importance for the evolution of B _{ L } because γ _{ q }∂_{ p } Lq + γ _{ T } S _{ p } ≫ ∂_{ p } m (note that ∂_{ p } m = ∂_{ p } Lq − S _{ p }). One example of adiabatic motions is steady isentropic lifting (Fig. 6e). In this case we have a cancellation between horizontal temperature advection and adiabatic DSE advection −c _{ p }{v ⋅ ∇T}_{ L } = −ω _{ aL } S _{ p }, so that the B _{ L } tendency is purely due to vertical moisture advection from isentropic lifting. As a result, isentropic lifting destabilizes the atmosphere through moistening of the troposphere from adiabatic lifting. Scaling of the horizontal temperature advection term indicates that a temperature gradient on the order of 1 K (1000 km)^{−1} is sufficient for isentropic lifting to be of leadingorder importance. Climatological gradients of this magnitude are observed over the South Asian and African monsoons (Kiladis et al. 2006; Adames and Ming 2018; Russell et al. 2020).
Another example of adiabatic motions arises from gravity waves (Fig. 6f). In this case, adiabatic lifting cools the LFT, increasing undilute B _{ L }. It also moistens the lower troposphere as moist DBL air is advected upward [term C in Eq. (46)]. By comparing the magnitude of the two contributions (Table 1 and Fig. B1), it can be shown that the cooling contribution is ~3.5 times larger, so that it contributes to ~78% of the increase in B _{ L }. While the vertical advection of DSE and Lq are of comparable magnitude, B _{ L } is more sensitive to changes in T _{ L } than q _{ L } (γ _{ T } > γ _{ q }), so that adiabatic cooling has a larger impact on B _{ L } than adiabatic moistening. Adiabatic cooling not only increases undilute B _{ L }, but also reduces the saturation deficit by reducing the saturation specific humidity of the LFT.
b. The plume buoyancy equation in DYNAMO soundings
Moisture and MSE budgets constructed from tropical field campaign data have yielded valuable insights related to the coupling between convection and the largescale environment (Sobel et al. 2014; Inoue and Back 2015; Hannah et al. 2016). However, B _{ L } is more closely related to the occurrence of deep convection, and a B _{ L } budget is therefore expected to offer insights not available from traditional moisture or MSE budgets. Observations or other datasets that provide quantities such as collocated surface fluxes, diabatic heating rates, Q _{1}, and Q _{2}, as well as state variables, can be leveraged to evaluate B _{ L } budgets via Eq. (46), and one example is shown below as a preliminary analysis of the B _{ L } budget using DYNAMO sounding data.
Figure 7 shows a time series of B _{ L }, TRMM 3B42 rainfall and the different terms in Eq. (46) obtained from a lag regression onto a B _{ L } time series calculated from DYNAMO sounding data (section 2c). A clear peak in B _{ L } is shown at lag 0, and TRMM 3B42 precipitation increases at this point. The maximum rainfall lags the B _{ L } maximum by 6 h. The time lag between B _{ L } and precipitation can likely be largely attributed to the widespread occurrence of stratiform precipitation that occurs concurrently with convection while B _{ L } is reduced. Indeed, Zuluaga and Houze (2013) showed that about 1/3 of the precipitation echo seen by radar at the time of peak precipitation during DYNAMO (at 2day time scales) was broad stratiform. TRMM 3B42 may overestimate rainfall in such scenarios where optically thick stratiform cloud contains light precipitation (Xu and Rutledge 2014). The terms that lead to the evolution of B _{ L } [Eq. (46)] are shown in the middle and bottom panels of Fig. 7. The maximum in B _{ L } is preceded by a positive ∂_{ t } B _{ L }, which attains a maximum amplitude at lag −3 h. The 3–12 h time scale for increase in B _{ L } is consistent with the time scale of convective buildup on short time scales during DYNAMO (third column of Table 3 in Powell and Houze 2013).
Lag regression based on a B _{ L } time series obtained from sounding data from DYNAMO northern array. (top) The B _{ L } (blue–green) and TRMM 3b42 rainfall rate (mm day^{−1}) for NSA. (middle) the LFT terms in Eq. (46). (bottom) The DBL terms. In the middle and bottom panels ∂_{ t } B _{ L } is shown as a dark orange dashed line, in units of m s^{−2} day^{−1}. Anomalies are scaled to a 0.025 m s^{−2} perturbation in B _{ L } at lag hour 0. Note that the axis for ∂_{ t } B _{ L } (right y axis) has been compressed to facilitate comparison between the terms in Eq. (46).
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
Lag regression based on a B _{ L } time series obtained from sounding data from DYNAMO northern array. (top) The B _{ L } (blue–green) and TRMM 3b42 rainfall rate (mm day^{−1}) for NSA. (middle) the LFT terms in Eq. (46). (bottom) The DBL terms. In the middle and bottom panels ∂_{ t } B _{ L } is shown as a dark orange dashed line, in units of m s^{−2} day^{−1}. Anomalies are scaled to a 0.025 m s^{−2} perturbation in B _{ L } at lag hour 0. Note that the axis for ∂_{ t } B _{ L } (right y axis) has been compressed to facilitate comparison between the terms in Eq. (46).
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
Lag regression based on a B _{ L } time series obtained from sounding data from DYNAMO northern array. (top) The B _{ L } (blue–green) and TRMM 3b42 rainfall rate (mm day^{−1}) for NSA. (middle) the LFT terms in Eq. (46). (bottom) The DBL terms. In the middle and bottom panels ∂_{ t } B _{ L } is shown as a dark orange dashed line, in units of m s^{−2} day^{−1}. Anomalies are scaled to a 0.025 m s^{−2} perturbation in B _{ L } at lag hour 0. Note that the axis for ∂_{ t } B _{ L } (right y axis) has been compressed to facilitate comparison between the terms in Eq. (46).
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
When considering the processes that lead to the evolution of B _{ L } (middle and bottom panels of Fig. 7), we find that largest contributor was term C. This result is consistent with Fig. 5, which shows that the time series of ω _{ a } and ω _{ Q } are often in quadrature. It indicates that convection was often preceded by adiabatic moistening and cooling of the LFT. The large amplitude of term C and the fast time scale in which destabilization occurred is suggestive of inertiogravity waves (Mapes 2000; Raymond and Fuchs 2007) similarly to features depicted in Fig. 6f and documented by Zuluaga and Houze (2013) and Yu et al. (2018) in DYNAMO radar data.
Following the maximum in B _{ L }, a negative ∂_{ t } B _{ L } is seen, with a maximum amplitude occurring at lag +3 h. At this time, both terms C and E become negative, indicating adiabatic compression and drying of the LFT by downdrafts. Both of these LFT processes contribute to roughly half of the removal of B _{ L } at this time. The other half of the stabilization comes from term H, indicating that downdrafts are stabilizing the DBL, as indicated by Eq. (49). The fact that the total negative B _{ L } tendency is equally split by contributions from the DBL and LFT underscores the importance of processes occurring within the whole lower troposphere in eliminating convective instability. Other terms in Eq. (46) are smaller in amplitude and contribute less to the evolution of B _{ L } at the time scale shown.
It is worth noting that B _{ L } remains positive in the twoday lag regression shown in Fig. 7, suggesting that the maximum in B _{ L } occurs within an envelope of enhanced B _{ L }. We investigate the processes that lead to this enhancement of B _{ L } by including a lag regression from lag day −10 to lag day 10, shown in Fig. 8, and smoothing all the terms by a 2day running mean. The evolution of B _{ L } and TRMM precipitation follow an analogous pattern to that shown in Fig. 7, with TRMM precipitation slightly lagging B _{ L }. The B _{ L } tendency follows a similar evolution to that seen in Fig. 7, but stretched over the 10day period. The tendency is a maximum 1–2 days prior to the maximum in B _{ L }. A minimum in the B _{ L } tendency is observed at lag +1 days. The time scale in which B _{ L } builds up and is eliminated is consistent with the time scale of the MJO events observed during DYNAMO.
As in Fig. 7, but ranging from lag day −10 to lag day 10. All the lines shown are smoothed with a 2day running mean. Note that the limits of the abscissa in the middle and bottom panels are smaller than those in Fig. 7.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
As in Fig. 7, but ranging from lag day −10 to lag day 10. All the lines shown are smoothed with a 2day running mean. Note that the limits of the abscissa in the middle and bottom panels are smaller than those in Fig. 7.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
As in Fig. 7, but ranging from lag day −10 to lag day 10. All the lines shown are smoothed with a 2day running mean. Note that the limits of the abscissa in the middle and bottom panels are smaller than those in Fig. 7.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
When considering the processes that lead to the evolution of B _{ L }, we see that the dominant terms are terms A and C. Both these terms closely follow the B _{ L } tendency. Term C is the largest B _{ L } source term prior to the B _{ L } maximum, followed by term A. In contrast, term A is the largest sink of B _{ L } after B _{ L } is a maximum. Overall, it can be argued that both terms contribute nearly equally to the evolution of B _{ L } at these time scales.
Besides terms A and C, it can be seen that term G plays an important role in the maintenance of the B _{ L } on time scales longer than a day. It is well documented that radiative heating plays an important role in the MJO, but it is interesting that is the DBL radiative heating contribution that dominates at lag day 0. This result is consistent with Fig. 8 of Wolding et al. (2016). A small, but positive contribution from term D is seen at lag days +2–7, likely in association with elevated stratiform convection, in agreement with previous studies (Wang et al. 2016; Ciesielski et al. 2017).
Last, we find that terms H and I are significant, but these DBL contributions along with term G largely cancel one another (bottom panel of Fig. 8) such that the evolution of ∂_{ t } B _{ L } matches the LFT terms A and C fairly closely (middle panel of Fig. 8). That terms H and I nearly cancel one another is qualitatively consistent with the notion that the DBL adjusts to BLQE on a time scale shorter than a day (Raymond 1995; Hansen et al. 2019).
Term E does not contribute significantly to the buildup of B
_{
L
}, but is a significant contributor to the removal of B
_{
L
} from lag day −1 to lag day +7. When terms E, H, and I are considered together as
As in the top panel of Fig. 8, except the black line shows
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
As in the top panel of Fig. 8, except the black line shows
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
As in the top panel of Fig. 8, except the black line shows
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
6. Conclusions
In this study, we build upon a growing literature that indicates that tropical deep convection is sensitive to fluctuations in convective instability (CAPE/CIN) and that moisture plays an important role in its occurrence and organization (Fuchs and Raymond 2002; Frierson et al. 2004; Kuo et al. 2017; among others). The framework is based on the precipitation–buoyancy relationship developed by AN18 and expanded upon by AAN. In this framework, precipitation is proportional to the buoyancy of a plume that rises above the DBL (carrying its θ _{ e }) and entrains environmental air through a deep inflow profile as it ascends in the LFT.
The precipitation–buoyancy relation is used to consolidate the “relaxed” WTG approximation (Sobel et al. 2001; Raymond and Zeng 2005), which allows for temperature fluctuations and prognostic moisture into a single “plume buoyancy equation” [Eq. (46)]. This equation shows that fluctuations in moisture and temperature contribute roughly equally to the evolution of the LFT averaged plume buoyancy B _{ L } [Eq. (47)].
We examined the processes that lead to the evolution of B _{ L } using sounding data from DYNAMO sounding arrays (Ciesielski et al. 2014). On the basis of linear regression analysis, we find that, at the time scale of a day, the buildup of B _{ L } is associated with cooling and moistening of the LFT from departures from WTG, likely in association with adiabatic lifting from gravity waves. The peak in the B _{ L } tendency occurs 3 h prior to the B _{ L } maximum. TRMM 3B42 rainfall peaks up to 6 h after the maximum in B _{ L }. While this temporal shift could be due to discrepancies in the datasets, a delay in the peak and subsequent decay of precipitation suggests that short time scales other than those accounted for by B _{ L } may play a role in rainfall, such as the time for deep convective plume growth, production of stratiform rain, or dynamical effects. After the peak in B _{ L }, warming and drying by departures from WTG, likely from adiabatic compression, and drying of the LFT and DBL by convection eliminate B _{ L }.
When considering variations in B _{ L } at the intraseasonal time scale, we find that horizontal moisture advection and cooling and moistening from adiabatic lifting in the LFT are the dominant terms contributing to ∂_{ t } B _{ L }. The former process has been well documented for the MJO events that occurred during DYNAMO (Sobel et al. 2014; Yokoi and Sobel 2015). The latter is somewhat surprising, and suggests that departures from WTG may play an important role in the MJO. However, it is unclear whether adiabatic lifting is important only for the MJO events that were documented during DYNAMO since these events exhibited a time scale closer to 30 days, a shorter time scale than what is considered typical for the MJO (40–50 days). A scale analysis study by Adames et al. (2019) indicates that the amplitude of adiabatic motions (i.e., the temperature tendency) scale with the square of the wave’s phase speed [see their Eq. (14)]. Thus, adiabatic motions should be much weaker in slowerpropagating MJO events. However, studies by Powell (2016, 2017) and Haertel et al. (2015) have suggested that circumnavigating Kelvin waves, and the cooling of the troposphere associated with them, play an important role in the initiation of at least some MJO events. Considering that the DYNAMO array was located in the region of MJO initiation, then this result may not be surprising, as it may reflect the signature of such a wave before it becomes more coupled to diabatic heating. Nonetheless, this is a result that warrants further study.
There are a few caveats to this study. The centered difference scheme applied to the estimation of the static stability S _{ p } can result in numerical errors in the calculation of ω _{ Q }, which in turn results in uncertainty in the largescale adiabatic vertical motion term ω _{ a }. Studies that use socalled spectral WTG method of calculating diabatically induced vertical motions indicate that traditional means tend to underestimate ω _{ Q } (Herman and Raymond 2014; Wang et al. 2016). Thus, ω _{ Q } may be underestimated and ω _{ a } overestimated in this study. Comparing the difference between the terms calculated via Eq. (46) using the methods shown herein versus using spectral WTG may be useful.
It is important to note that our equations do not have an explicit representation of different cloud types, opting for a simplified representation of convection based on a single deep plume. Even though we do not make this distinction, the equations for B _{ L } bear some similarity to the congestus and stratiform heating equations described by Khouider and Majda (2006, 2008) [their Eqs. (2.7) and (2.8)]. Additionally, our Eq. (6) resembles their equation for precipitation [their Eqs. (2.6) and (2.9)]. The framework presented here does not only validate the foundations of the multicloud models used by Khouider and Majda (2006, 2008) and Khouider and Majda (2016) (among others) but also can be used to tune their model and understand the processes in which these models represent convectively coupled waves.
In addition, we do not include uppertropospheric processes in this framework for simplicity. Schiro et al. (2018) showed that including the upper troposphere’s contribution to the precipitation–buoyancy relation yields similar results to just analyzing the lower troposphere. This result is valid when entrainment in this layer is small, so the effects of uppertropospheric humidity are small. The similarity of these results occurs because temperature tends to be strongly correlated between the upper and lower free troposphere. Integrating the plume buoyancy through the entire troposphere changes its magnitude by including virtual temperature difference between plume environment in the upper layer, but most of the variability is not independent of that included in the lowertropospheric variables. The adjustment processes discussed here in terms of lowerfreetropospheric variables should thus extend to the full troposphere. It is possible that at the next level of refinement, independent fluctuations of temperature, humidity, vertical motions and freezing processes including aerosol effects (Rosenfeld et al. 2008) in the upper troposphere can play a significant role in precipitation, as has been suggested for the observed relationship between the MJO and the QBO (Son et al. 2017; Martin et al. 2019). Developing empirical and theoretical formulations to include these could be a fruitful direction for future work.
Finally, although the B _{ L } equation from which we derive its tendency can be utilized to predict the destabilization of the troposphere that leads to convective onset, the framework presented only considers the thermodynamic structure of the atmosphere on precipitation. However, dynamic factors, such as lowlevel convergence or vertical wind shear, also impact vertical velocities found in convective updrafts, and thus probably contribute to the variance in rain rates that is observed at high B _{ L }. Thus, while our framework is useful for describing how various processes that alter the thermodynamic structure of the atmosphere impact rainfall, a more complete expression for tropical precipitation evolution would include both the B _{ L } tendency and other terms related to largescale dynamic processes.
While the “plume buoyancy equation” may seem complex, especially when compared to budgets such as columnintegrated MSE, it nonetheless provides a lucid picture of convective destabilization. Furthermore, it shows that departures from quasi equilibrium or WTG can play a key role in the destabilization of the troposphere, which is not as clear when analyzing MSE budgets (see Table 3). Related theoretical work in AAN elaborates on the time scales of adjustment of B toward QE or WTG. Thus, using the plume buoyancy equation is complementary to using standard moisture or MSE budgets when considering convectively coupled tropical motion systems, providing insight into the fluctuations that yield strong convective events, while the MSE budget can be helpful for understanding the envelope of these (AAN; Neelin and Yu 1994; Adames et al. 2019; Ahmed et al. 2021). The framework presented here may also be used for intercomparison of model representation of tropical deep convection.
Acknowledgments
ÁFA was supported by the National Science Foundation (NSF) Grant AGS1841559. FA And JDN were supported in part by NSF AGS1540518 and National Atmospheric and Oceanic Administration Grant NA18OAR4310280. ÁFA would like to thank Brandon Wolding, Kuniaki Inoue, Hannah Zanowski, and Daehyun Kim for discussions that helped organize the contents of this manuscript. We also would like to thank George Kiladis for providing us most of the sounding data used in this study.
APPENDIX A
Derivation of the Plume Buoyancy
APPENDIX B
Soundings
The top row in Fig. B1 shows the vertical profiles of variability in T, q, and h. It shows that variability in temperatures in the tropics are on the order of 1 K, while variability in q is on the order of 1 g kg^{−1}. The vertical gradients of DSE, Lq and MSE are shown in the bottom row. The mean gradients are very similar across the three sounding sites analyzed. The DSE and Lq vertical gradients are on the order of 0.5 J kg^{−1} Pa^{−1} over the LFT, while the layeraveraged MSE gradient is on the order of 0.1 J kg^{−1} Pa^{−1}. We can estimate M _{ B } from Fig. B1f to be on the order of 0.3 J kg^{−1}. These estimates are used in the scaling of the equations shown in Table 3 and Eq. (47).
(top) Standard deviation of dailymean (a) internal energy c _{ p } T, (b) latent energy Lq, and (c) moist enthalpy (h) from sounding data from Majuro (7.1°N, 171.4°E; blue), Manaus (3.2°S, 59.9°W; gray) and Gan (0.6°S, 73.1°E; red). (bottom) As in the top row, but showing the climatologicalmean vertical profiles of (d) DSE (s), (e) Lq, and (f) MSE. The mean vertical gradient of the profiles, in units of J kg^{−1} Pa^{−1}, in (d)–(f) is shown as a dot–dashed line.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
(top) Standard deviation of dailymean (a) internal energy c _{ p } T, (b) latent energy Lq, and (c) moist enthalpy (h) from sounding data from Majuro (7.1°N, 171.4°E; blue), Manaus (3.2°S, 59.9°W; gray) and Gan (0.6°S, 73.1°E; red). (bottom) As in the top row, but showing the climatologicalmean vertical profiles of (d) DSE (s), (e) Lq, and (f) MSE. The mean vertical gradient of the profiles, in units of J kg^{−1} Pa^{−1}, in (d)–(f) is shown as a dot–dashed line.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
(top) Standard deviation of dailymean (a) internal energy c _{ p } T, (b) latent energy Lq, and (c) moist enthalpy (h) from sounding data from Majuro (7.1°N, 171.4°E; blue), Manaus (3.2°S, 59.9°W; gray) and Gan (0.6°S, 73.1°E; red). (bottom) As in the top row, but showing the climatologicalmean vertical profiles of (d) DSE (s), (e) Lq, and (f) MSE. The mean vertical gradient of the profiles, in units of J kg^{−1} Pa^{−1}, in (d)–(f) is shown as a dot–dashed line.
Citation: Journal of the Atmospheric Sciences 78, 2; 10.1175/JASD200074.1
APPENDIX C
Interpretation of the Vertical Transport of MSE by Convection
REFERENCES
Adames, Á. F. , 2017: Precipitation budget of the Madden–Julian oscillation. J. Atmos. Sci., 74, 1799–1817, https://doi.org/10.1175/JASD160242.1.
Adames, Á. F. , and Y. Ming , 2018: Moisture and moist static energy budgets of South Asian monsoon low pressure systems in GFDL AM4.0. J. Atmos. Sci., 75, 2107–2123, https://doi.org/10.1175/JASD170309.1.
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, 3863–3881, https://doi.org/10.1175/JASD190121.1.
Ahmed, F. , and J. D. Neelin , 2018: Reverse engineering the tropical precipitation–buoyancy relationship. J. Atmos. Sci., 75, 1587–1608, https://doi.org/10.1175/JASD170333.1.
Ahmed, F. , Á. F. Adames , and J. D. Neelin , 2020: Deep convective adjustment of temperature and moisture. J. Atmos. Sci., 77, 2163–2186, https://doi.org/10.1175/JASD190227.1.
Ahmed, F. , J. D. Neelin , and Á. F. Adames , 2021: Quasiequilibrium and weak temperature gradient balances in an equatorial betaplane model. J. Atmos. Sci., 78, 209–227, https://doi.org/10.1175/JASD200184.1.
Arakawa, A. , and W. H. Schubert , 1974: Interaction of a cumulus cloud ensemble with the largescale environment, part I. J. Atmos. Sci., 31, 674–701, https://doi.org/10.1175/15200469(1974)031<0674:IOACCE>2.0.CO;2.
Betts, A. K. , 1975: Parametric interpretation of tradewind cumulus budget studies. J. Atmos. Sci., 32, 1934–1945, https://doi.org/10.1175/15200469(1975)032<1934:PIOTWC>2.0.CO;2.
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/15200469(1989)046<0740:GWCSAD>2.0.CO;2.
Bretherton, C. S. , and A. H. Sobel , 2002: A simple model of a convectively coupled Walker circulation using the weak temperature gradient approximation. J. Climate, 15, 2907–2920, https://doi.org/10.1175/15200442(2002)015<2907:ASMOAC>2.0.CO;2.
Bretherton, C. S. , and A. H. Sobel , 2003: The Gill model and the weak temperature gradient approximation. J. Atmos. Sci., 60, 451–460, https://doi.org/10.1175/15200469(2003)060<0451:TGMATW>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/15200442(2004)017<1517:RBWVPA>2.0.CO;2.
Charney, J. G. , 1963: A note on largescale motions in the tropics. J. Atmos. Sci., 20, 607–609, https://doi.org/10.1175/15200469(1963)020<0607:ANOLSM>2.0.CO;2.
Chikira, M. , 2014: Eastwardpropagating intraseasonal oscillation represented by Chikira–Sugiyama cumulus parameterization. Part II: Understanding moisture variation under weak temperature gradient balance. J. Atmos. Sci., 71, 615–639, https://doi.org/10.1175/JASD13038.1.
Ciesielski, P. E. , and Coauthors, 2014: Qualitycontrolled upperair sounding dataset for DYNAMO/CINDY/AMIE: Development and corrections. J. Atmos. Oceanic Technol., 31, 741–764, https://doi.org/10.1175/JTECHD1300165.1.
Ciesielski, P. E. , R. H. Johnson , X. Jiang , Y. Zhang , and S. Xie , 2017: Relationships between radiation, clouds, and convection during DYNAMO. J. Geophys. Res. Atmos., 122, 2529–2548, https://doi.org/10.1002/2016JD025965.
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/JCLID1100384.1.
de Rooy, W. C. , and Coauthors, 2013: Entrainment and detrainment in cumulus convection: An overview. Quart. J. Roy. Meteor. Soc., 139, 1–19, https://doi.org/10.1002/qj.1959.
de Szoeke, S. P. , 2018: Variations of the moist static energy budget of the tropical Indian Ocean atmospheric boundary layer. J. Atmos. Sci., 75, 1545–1551, https://doi.org/10.1175/JASD170345.1.
Donner, L. J. , and V. T. Phillips , 2003: Boundary layer control on convective available potential energy: Implications for cumulus parameterization. J. Geophys. Res., 108, 4701, https://doi.org/10.1029/2003JD003773.
Durre, I. , R. S. Vose , and D. B. Wuertz , 2006: Overview of the Integrated Global Radiosonde Archive. J. Climate, 19, 53–68, https://doi.org/10.1175/JCLI3594.1.
Emanuel, K. , 1993: The effect of convective response time on WISHE modes. J. Atmos. Sci., 50, 1763–1776, https://doi.org/10.1175/15200469(1993)050<1763:TEOCRT>2.0.CO;2.
Emanuel, K. , 1995: The behavior of a simple hurricane model using a convective scheme based on subcloudlayer entropy equilibrium. J. Atmos. Sci., 52, 3960–3968, https://doi.org/10.1175/15200469(1995)052<3960:TBOASH>2.0.CO;2.
Emanuel, K. , 2019: Inferences from simple models of slow, convectively coupled processes. J. Atmos. Sci., 76, 195–208, https://doi.org/10.1175/JASD180090.1.
Emanuel, K. , J. D. Neelin , and C. S. Bretherton , 1994: On largescale circulations in convecting atmospheres. Quart. J. Roy. Meteor. Soc., 120, 1111–1143, https://doi.org/10.1002/qj.49712051902.
Frierson, D. M. , A. J. Majda , and O. M. Pauluis , 2004: Large scale dynamics of precipitation fronts in the tropical atmosphere: A novel relaxation limit. Commun. Math. Sci., 2, 591–626, https://doi.org/10.4310/CMS.2004.v2.n4.a3.
Fuchs, Ž. , and D. J. Raymond , 2002: Largescale modes of a nonrotating atmosphere with water vapor and cloud–radiation feedbacks. J. Atmos. Sci., 59, 1669–1679, https://doi.org/10.1175/15200469(2002)059<1669:LSMOAN>2.0.CO;2.
Gao, S. , L. Ran , and X. Li , 2006: Impacts of ice microphysics on rainfall and thermodynamic processes in the tropical deep convective regime: A 2D cloudresolving modeling study. Mon. Wea. Rev., 134, 3015–3024, https://doi.org/10.1175/MWR3220.1.
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 MaddenJulian oscillation. Quart. J. Roy. Meteor. Soc., 141, 1586–1602, http://doi.org/10.1002/qj.2461.
Hannah, W. M. , 2017: Entrainment versus dilution in tropical deep convection. J. Atmos. Sci., 74, 3725–3747, https://doi.org/10.1175/JASD160169.1.
Hannah, W. M. , B. E.<