1. Introduction
Most early theories that aimed to explain the existence and properties of the MJO treated it as some form of atmospheric equatorial Kelvin wave, modified by interaction of moist convection with the larger-scale flow. In recent years, it has become clear that convectively coupled Kelvin waves do exist, but that the MJO occupies a distinctly different region of the spectrum (Wheeler and Kiladis 1999). It has also been realized that moisture plays a greater role than it did in the early theories. MJO simulation in particular is improved when convection is made more sensitive to environmental moisture (e.g., Tokioka et al. 1988; Wang and Schlesinger 1999; Maloney and Hartmann 2001; Benedict and Randall 2009; Hannah and Maloney 2011; Kim et al. 2012). Studies with limited-domain cloud-resolving models have shown that, at least under mean conditions of large-scale ascent and relatively frequent deep convection, moisture contains most of the memory in the atmosphere that can regulate convection on time scales longer than a few days (Tulich and Mapes 2010; Kuang 2010).
The idea has emerged that the MJO is a moisture mode. We gave our own definition of this term in a previous study (Sobel and Maloney 2012, hereafter SM12), in which we also proposed a particular very simple model of a moisture mode with the aim of capturing the essential features of the MJO. That model did not appear to be particularly successful. The only linear unstable modes were westward propagating. A nonlinear mode also discussed in SM12 differed substantially in structure from the MJO, having a shocklike jump in moisture at the leading edge of an active phase, whereas the actual MJO has a gradual buildup. Here, we consider extensions to this model to incorporate additional processes that have been hypothesized to be important to the MJO, and which one might expect to cause eastward propagation. These processes include advection of a mean zonal moisture gradient by perturbation zonal winds, modulation of synoptic-scale transient eddy drying by the MJO-scale zonal wind, and frictional convergence.
2. Model framework
a. Basic equations
















b. New processes






Another potentially relevant linear process is zonal MSE advection in the presence of a background eastward MSE gradient, which would contribute a term of the form










The combined effect of perturbations in latent heat flux and the other processes described above—modulation of synoptic-eddy drying, zonal advection, and frictional convergence—is expressed by the quantity Cu − D in (10). If that quantity is positive (as in SM12 where D = 0), low-level westerlies are associated with an MSE source and easterlies with an MSE sink. If it is negative, the converse is true.
Based on both recent studies with reanalysis data (e.g., Kiranmayi and Maloney 2011) and general circulation models (Maloney 2009), it appears that the sum of terms including horizontal advection and vertical advection contributes about 10 W m−2 of MSE increase per 1 m s−1 of zonal-wind anomaly in MJO easterly regions when averaged about the equator. This total exceeds the magnitude of the surface flux anomalies by about a third, so that Cu − D < 0. The reanalysis vertical and horizontal advection contributions are comparable, although vertical advection is somewhat larger, to an extent that depends on the reanalysis product (Kiranmayi and Maloney 2011). In the western Pacific warm pool, the synoptic-eddy advection mechanism dominates the advection by the anomalous wind across the mean zonal humidity gradient, whereas in the Indian Ocean the converse is true. The reanalysis MSE budgets contain relatively large analysis increments that are of the same order as the horizontal and vertical advection terms (Kiranmayi and Maloney 2011), limiting the precision with which these values can be constrained. An implication of our analysis is that tighter observational constraints on these various MSE sources would be valuable.










As the phase relationship between wind and water vapor is central to our discussion, in Fig. 1 we present several aspects of this relationship graphically; some aspects of this were mentioned in SM12. (See Table 1 of SM12 for all parameter values other than those stated here.) In the left panel the quantity G(x | 0) is plotted (with A = 1, for this panel only); this is the wind response to a δ-function heating at the origin. In the middle panel we show the wind response to a sinusoidal heating of wavenumber 3 [P = cos(3πx/Lm) on our domain −Lm < x < Lm, with Lm = 2 × 104 km]; wavenumber 3 is chosen simply as an example. In the right panel we show the amplitude and phase of the function Γ(k); positive phase means that wind lags precipitation, as it does here for all values of k. As k → ∞ (and also as k → 0), the amplitude |Γ(k)| → 0 as the phase goes to π/2. The wavenumber of maximum amplitude and minimum phase lag depends on L, and will be smaller than the value shown in Fig. 1 if L is chosen greater than 1500 km. At the minimum phase lag, the correlation of W and MSE sources correlated with u is greatest, yielding the strongest growth due to these terms if Cu − D > 0 (SM12), or strongest damping by these terms otherwise.

(left) The projection function G(x, 0)—that is, the wind response to a δ-function heating at x′ = 0. (middle) The wind response to a sinusoidal heating of wavenumber 3, together with that heating itself, as functions of x. (right) The amplitude and phase of the wind response function for sinusoidal heating Γ(k), as a function of wavenumber k. The scale for amplitude is on the left axis, and the scale for phase (°) is on the right axis.
Citation: Journal of the Atmospheric Sciences 70, 1; 10.1175/JAS-D-12-0189.1

(left) The projection function G(x, 0)—that is, the wind response to a δ-function heating at x′ = 0. (middle) The wind response to a sinusoidal heating of wavenumber 3, together with that heating itself, as functions of x. (right) The amplitude and phase of the wind response function for sinusoidal heating Γ(k), as a function of wavenumber k. The scale for amplitude is on the left axis, and the scale for phase (°) is on the right axis.
Citation: Journal of the Atmospheric Sciences 70, 1; 10.1175/JAS-D-12-0189.1
(left) The projection function G(x, 0)—that is, the wind response to a δ-function heating at x′ = 0. (middle) The wind response to a sinusoidal heating of wavenumber 3, together with that heating itself, as functions of x. (right) The amplitude and phase of the wind response function for sinusoidal heating Γ(k), as a function of wavenumber k. The scale for amplitude is on the left axis, and the scale for phase (°) is on the right axis.
Citation: Journal of the Atmospheric Sciences 70, 1; 10.1175/JAS-D-12-0189.1



Now if the sign of this second term in (17) reverses—if the processes represented by D depend on the zonal wind more strongly than the latent heat flux does, so that the net effect of easterlies is moistening and westerlies drying—both these effects reverse. The net effect of all these processes is then to cause weak eastward propagation (relative to the mean flow) but also to induce damping rather than growth. Observations, discussed above, suggest (with considerable uncertainty) that this may be the relevant regime for the Indian Ocean, and perhaps also the far western Pacific. In this regime the wind–evaporation feedback alone is still destabilizing, but not enough so to overcome the other wind-related feedbacks whose signs are opposite to it. Growth of disturbances can still occur, but only if Meff < 0.
The left panel of Fig. 2 shows the growth rate and phase speed as functions of k, for r = 0.15, so that

Phase speed and growth rate for zonal moisture diffusivity (left) kw = 0 and (right) kw = 104 m2 s−1.
Citation: Journal of the Atmospheric Sciences 70, 1; 10.1175/JAS-D-12-0189.1

Phase speed and growth rate for zonal moisture diffusivity (left) kw = 0 and (right) kw = 104 m2 s−1.
Citation: Journal of the Atmospheric Sciences 70, 1; 10.1175/JAS-D-12-0189.1
Phase speed and growth rate for zonal moisture diffusivity (left) kw = 0 and (right) kw = 104 m2 s−1.
Citation: Journal of the Atmospheric Sciences 70, 1; 10.1175/JAS-D-12-0189.1
The maximum growth rate as k → ∞ for kw = 0 is not an “ultraviolet catastrophe” in that the growth rate does not blow up, but asymptotes to a constant value
As stated above, instability occurs only if
3. Conclusions
We have presented a brief analysis of a simple linear moisture-mode model. This is a modification of that in Sobel and Maloney (2012), including a simple representation of processes that have been proposed as contributing to eastward propagation of the MJO: modulation of synoptic-eddy drying by the MJO-scale zonal winds, advection of a mean zonal moisture gradient, and frictional convergence. All of these are crudely parameterized as sources of moist static energy that are proportional to minus the low-level zonal wind perturbation. Our primary conclusions are as follows:
In order for eastward propagation of moisture modes to occur in this model, the other processes causing moistening in low-level easterlies must be stronger than the drying associated with suppression of surface fluxes in perturbation easterlies.
If the requirement for eastward propagation is met, this also implies that the net effect of all processes directly related to the low-level zonal wind on growth is negative: that is, damping. It is still the case that the surface wind–evaporation feedback by itself is destabilizing, but that destabilization is exceeded by the damping associated with the other zonal wind–related processes.
Eastward-propagating modes can only be unstable because of cloud–radiative feedback or gross moist instability. The effective gross moist stability, which captures the net effect of both, must be negative.
Absent any scale-selective damping, the growth rates of eastward-propagating unstable modes are maximized at the largest and smallest spatial scales in the system. There is a minimum in between at synoptic scales (precisely those where the maximum occurred in SM12).
A modest horizontal diffusion, or other scale-selective damping, causes the growth rate to decrease monotonically with wavenumber so that the largest zonal scales are selected.
Acknowledgments
This research was sponsored by National Science Foundation Grant AGS-1062206 (AHS) AGS-1062161 (EDM) and AGS-1025584 (EDM), National Aeronautics and Space Administration Grant NNX09AK34G (AHS), National Oceanic and Atmospheric Administration Grants NA08OAR4320912 A6R (AHS) and NA08OAR4320893 #7 and #14 (EDM), and Office of Naval Research Grant N00014-12-1-0911 (AHS).
REFERENCES
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.
Benedict, J. J., and D. A. Randall, 2009: Structure of the Madden–Julian oscillation in the superparameterized CAM. J. Atmos. Sci., 66, 3277–3296.
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.
Haertel, P. T., G. N. Kiladis, A. Denno, and T. M. Rickenback, 2008: Vertical-mode decompositions of 2-day waves and the Madden–Julian oscillation. J. Atmos. Sci., 65, 813–833.
Hannah, W. M., and E. D. Maloney, 2011: The role of moisture–convection feedbacks in simulating the Madden–Julian oscillation. J. Climate, 24, 2754–2770.
Hartmann, D. L., and E. D. Maloney, 2001: The Madden–Julian oscillation, barotropic dynamics, and North Pacific tropical cyclone formation. Part I: Observations. J. Atmos. Sci., 58, 2545–2558.
Hendon, H. H., and M. L. Salby, 1994: The life cycle of the Madden–Julian oscillation. J. Atmos. Sci., 51, 2225–2237.
Hsu, P.-C., and T. Li, 2012: Role of the boundary layer moisture asymmetry in causing the eastward propagation of the Madden–Julian oscillation. J. Climate, 25, 4914–4931.
Kim, D., A. H. Sobel, A. D. D. 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.
Kiranmayi, L., and E. D. Maloney, 2011: Intraseasonal moist static energy budget in reanalysis data. J. Geophys. Res.,116, D21117, doi:10.1029/2011JD016 031.
Kuang, Z., 2010: Linear response functions of a cumulus ensemble to temperature and moisture perturbations and implications for the dynamics of convectively coupled waves. J. Atmos. Sci., 67, 941–962.
Lin, J., and B. E. Mapes, 2004: Radiation budget of the tropical intraseasonal oscillation. J. Atmos. Sci., 61, 2050–2062.
Maloney, E. D., 2009: The moist static energy budget of a composite tropical intraseasonal oscillation in a climate model. J. Climate, 22, 711–729.
Maloney, E. D., and D. L. Hartmann, 2001: The Madden–Julian oscillation, barotropic dynamics, and North Pacific tropical cyclone formation. Part I: Observations. J. Atmos. Sci., 58, 2545–2558.
Maloney, E. D., and M. J. Dickinson, 2003: The intraseasonal oscillation and the energetics of summertime tropical western North Pacific synoptic-scale disturbances. J. Atmos. Sci., 60, 2153–2168.
Moskowitz, B. M., and C. S. Bretherton, 2000: An analysis of frictional feedback on a moist equatorial Kelvin mode. J. Atmos. Sci., 57, 1797–1816.
Raymond, D. J., and Z. Fuchs, 2009: Moisture modes and the Madden–Julian oscillation. J. Climate, 22, 3031–3046.
Sobel, A. H., and E. D. Maloney, 2000: Effect of ENSO and the MJO on western North Pacific tropical cyclones. Geophys. Res. Lett., 27, 1739–1742.
Sobel, A. H., and E. D. Maloney, 2012: An idealized semi-empirical framework for modeling the Madden–Julian oscillation. J. Atmos. Sci., 69, 1691–1705.
Su, H., and J. D. Neelin, 2002: Teleconnection mechanisms for tropical Pacific descent anomalies during El Niño. J. Atmos. Sci., 59, 2694–2712.
Tokioka, T., K. Yamazaki, A. Kitoh, and T. Ose, 1988: The equatorial 30–60 day oscillation and the Arakawa–Schubert penetrative cumulus parameterization. J. Meteor. Soc. Japan, 66, 883–901.
Tulich, S., and B. E. Mapes, 2010: Transient environmental sensitivities of explicitly simulated tropical convection. J. Atmos. Sci., 67, 923–940.
Wang, B., 1988: Dynamics of tropical low-frequency waves: An analysis of the moist Kelvin wave. J. Atmos. Sci., 45, 2051–2065.
Wang, B., and H. Rui, 1990: Dynamics of the coupled moist Kelvin–Rossby wave on an equatorial β-plane. J. Atmos. Sci., 47, 397–61.
Wang, W., and M. E. Schlesinger, 1999: The dependence of convective parameterization of the tropical intraseasonal oscillation simulated by the UIUC 11-layer atmospheric GCM. J. Climate, 12, 1423–1457.
Wheeler, M., and G. N. Kiladis, 1999: Convectively coupled equatorial waves: Analysis of clouds and temperature in the wavenumber-frequency domain. J. Atmos. Sci., 56, 374–399.