## 1. Introduction

The Madden–Julian oscillation (MJO) is an atmospheric enigma [Hartmann and Hendon 2007; see Zhang (2005) for a review]. Initiated over the Indian Ocean exhibiting growth of convective activity, an enhanced oceanwide horizontal convergence in the lower troposphere, and an upper-troposphere divergence, this planetary-scale system propagates eastward with a speed of 3–6 m s^{−1} (Madden and Julian 1972). Accelerating its propagation speed as it crosses the date line with convective activity dying out, the MJO completes a full cycle traveling around the globe in about 30–60 days.

Its propagating structure consists of a well-defined pair of solitary vortices, as shown in Fig. 1 (cf. Yanai et al. 2000). Its longitudinal scale is about 3000 km, although the composite studies tend to suggest larger scales [e.g., see Figs. 1–5 of Hendon and Salby (1994), Fig. 2 of Kiladis et al. (2005), Fig. 6 of Adames et al. (2014), and Fig. 1 of Žagar and Franzke (2015)]. As a pair snapshot in Fig. 1 suggests, the phase relationship between convection (OLR field) and the vortex pair tends to fluctuate with time, which contribute to the smoothing effect of the vorticity field, when convection is used as a reference variable. Even in the latter, a core part of the vortex pair tends to represent a comparable longitudinal scale as the snapshots within a factor of 3 (e.g., Fig. 2 of Kiladis et al. 2005). Even a closer correspondence in scale is identified by a pair of the potential-vorticity anomalies in Fig. 9 of Zhang and Ling (2012). Thus, the present paper pursues the explanation of the MJO in terms of its vortex pair structure.

The major challenge of the MJO theories is to explain the slow eastward propagation of this structure. The Kelvin wave, identified as the best candidate under the dry linear limit, propagates too fast (50 m s^{−1}). Thus, a particular additional mechanism must be invoked. Since the tropical atmosphere displays an abundance of moist cumulus convective clouds, it is naturally identified as a prime candidate for providing such a slowdown mechanism. Because of this, vast investments have been made over more than three decades attempting to understand the MJO in terms of an interaction between the large-scale atmospheric flow and moist convection (cf. Wang 2005). The present study, alternatively, suggests that nonlinearity can equally be invoked for obtaining qualitatively an equivalent effect.

Theoretical efforts to explain the MJO in terms of an interaction between the large-scale atmospheric flow and moist convection are numerous to provide an exhaustive list herein, but they include Hayashi (1970), Lindzen (1974), Emanuel (1987), Yano and Emanuel (1991), Fuchs and Raymond (2007), Raymond and Fuchs (2007, 2009), Majda and Stechmann (2009), Yang and Ingersoll (2011, 2013), and Adames and Kim (2016). These theories usually treat the MJO as a linear wave. Unfortunately, these theories have difficulties in explaining both the scale and the propagation speed of MJO. A closer look at some successes often reveals extensive parameter tuning.

Earlier studies (Hayashi 1970; Lindzen 1974) simply fail to identify the MJO as a separate mode from a faster eastward-propagating Kelvin mode. The WISHE instability only predicts the phase speeds faster than 10 m s^{−1} even with its most sophisticated version (Fig. 2 of Yano and Emanuel 1991). The moisture mode (Raymond and Fuchs 2007) is a promising possibility with substantially slower phase velocities. One of the most sophisticated versions predicts the phase velocity of about 6 m s^{−1} (Fig. 2 of Fuchs and Raymond 2007) with the longitudinal wavenumber 2. However, it fails to explain the preferred scale of the MJO (i.e., wavenumber 2) with the growth rate approximately constant independent of the wavenumber. The skeleton model (Majda and Stechmann 2009; Thual and Majda 2016) is apparently successful in explaining the phase speed of 3–6 m s^{−1} at the longitudinal wavenumber 2 (Fig. 2 of Majda and Stechmann 2009). However, a close look at their formulation reveals that this is due to a careful tuning of their free parameter *Q*, which controls the gross moist stability (GMS: effective moist stability of the atmosphere): making the GMS sufficiently small but positive. This choice is rather limited, because observations suggest that GMS could often be negative (López Carrillo and Raymond 2005). Though clearly the convectively coupled wave theories continue to progress (e.g., Yang and Ingersoll 2011, 2013; Adames and Kim 2016), it is fair to say that those theories hardly suggest any tendency for converging.

All these theoretical studies are based on parameterized convection. From this perspective, the difficulties in simulating the MJO in global models (cf. Jiang et al. 2015) may also be attributed to the fact that moist convection is parameterized. Thus, if the horizontal resolution is high enough with mesh sizes less than approximately 50 km, convective parameterization could be turned off, and the MJO would more easily be simulated. A successful simulation of MJO in a high-resolution global model by Japanese Earth Simulator (Miura et al. 2007) supports this expectation.

However, similar efforts with comparable resolutions have turned out to be less successful (Gustafson and Weare 2004; Monier et al. 2010; Hagos et al. 2011; Ray et al. 2011; Holloway et al. 2013). For example, Fig. 1 of Holloway et al. (2013) shows that none of the six simulations reproduces eastward propagation of convective precipitation that is seen in TRMM satellite observation. In contrast, dry large-scale circulations associated with MJO are well simulated, as shown by their Figs. 6–8 for the zonal winds. In spite of their decaying tendencies, eastward propagation closely follows that of the observation in the phase space consisting of the two leading empirical principal components as shown by their Fig. 4. The same point is more succinctly presented by Fig. 2 of Monier et al. (2010), who compare the wavenumber–frequency spectra for the zonal wind (at 850-hPa height) and the top of the atmosphere outgoing longwave (infrared) radiation intensity (OLR), as a measure of convective variability, between the observation and a simulation. The model spectrum for the zonal wind agrees well with the observation, whereas the model OLR spectrum is about 3 times weaker.

Thus, in spite of their difficulties simulating convective fields correctly, models somehow reproduce MJO-associated large-scale circulations rather well. It further suggests a possibility of considering the large-scale dynamical component of MJO in isolation without coupling it with convection, because it appears that these large-scale MJO circulations can be reproduced without properly simulating convection.

In fact, it may be argued that the tropical atmosphere is not as convective as commonly believed. Both theoretical (Charney 1963; Yano and Bonazzola 2009) and observational (Yano et al. 2009; Adames et al. 2014; Žagar and Franzke 2015) analyses suggest that the tropical large-scale circulation is rather nondivergent to a leading-order approximation; that is, the tropical large-scale horizontal flows are dominated by vorticity. The point is probably best made with the example of the tropical cyclone: its maintenance fundamentally hinges on moist convection. However, as far as the leading-order dynamics is concerned, it is well understood in terms of pure vortex dynamics (cf. Yanai 1964). Its purely dynamical aspects, such as the trajectory, can well be predicted without moist convection (e.g., Kasahara and Platzman 1963). The present paper will argue that the same is also true for the MJO (cf. Adames et al. 2014; Žagar and Franzke 2015): its maintenance may still fundamentally hinge on moist convection. However, as far as the leading-order dynamics is concerned, the main feature to be explained is the eastward-propagating solitary vortex pair as shown in Fig. 1, or composite fields shown in references cited above, without moist convection.

Our argument above, to which we also follow in the present study, does not intend to dismiss another facet that some of the studies do produce the MJO well by adopting improved or better-tuned convection parameterizations. Especially, the improvement of the MJO simulations in GCMs by introducing superparameterization is dramatic (e.g., Benedict and Randall 2009; Pritchard and Bretherton 2014; Arnold and Randall 2015). However, its success has hardly led to any theoretical consensus nor does it negate the premise of dry dynamics playing the leading role. Thus, there is clearly scope for pursuing a completely different theory for the MJO. The successes of the superparameterization models and of the Earth Simulator hardly explain why so many other high-resolution explicit-convection models have difficulties for simulating the convective component of MJO in spite of relative success in simulating its “dry” component. The present paper presents an explanation for this general tendency.

An attempt to explain the MJO solely in terms of the dry dynamics itself is not an original idea, as this was already pursued by, for example, Straus and Lindzen (2000), Moncrieff (2004), Lin et al. (2007), and Wedi and Smolarkiewicz (2010). Among the aforementioned studies, Wedi and Smolarkiewicz (2010) is closest to the present study in its approach. Moncrieff (2004) proposes a somewhat similar idea, but from a point of view of considering the MJO as an organized convective archetype. Despite invoking the notion of “convection” here, Moncrieff’s paper does not consider convective heating effect explicitly.

In the next section, we will argue from a theoretical basis that the MJO must be considered a nonlinear Rossby wave, if it is to be understood in terms of the dry dynamics. The “modon” theory (Tribbia 1984a; Verkley 1984), identified as a candidate for the theory, is formulated in section 3. Numerical experiments of a global shallow-water system are presented in section 4 to verify the structural stability and the longevity of the modon solutions. The paper concludes with further discussion in section 5.

## 2. Theoretical considerations

How can we, then, explain the MJO in terms the dry dynamics? Under the linear limit, the large-scale tropical atmospheric dynamics consists of a set of the equatorial waves (Matsuno 1966; Yano and Bonazzola 2009). Thus, if the MJO is to be explained in terms of the linear dynamics, it must be explained in terms of the slowest-propagating equatorial planetary-scale waves: the Kelvin or the Rossby wave.

At first sight, the Kelvin wave appears to be more attractive, because it propagates eastward, as does the MJO. However, regrettably, the Kelvin wave is nondispersive, and its propagation speed is fixed by its vertical structure to about 50 m s^{−1}. That is far too fast for explaining the MJO.

Thus, we need to turn to the Rossby wave. The Rossby wave presents more flexibility owing to its dispersiveness, and its propagation speed can be adjusted by changing the horizontal wavenumber.

Equatorial Rossby waves also take the form of a vortex pair trapped along the equator as suggested in Fig. 1, which is fairly nondivergent (Delayen and Yano 2009). An observational study (Žagar and Franzke 2015) more directly shows by a mode decomposition that the MJO is dominated by Rossby modes. Here, note that explaining the equatorially trapped structure of the MJO is much less challenging than its propagation speed, because the former is already consistent with the dry linear equatorial Rossby wave structure.

*β*the latitudinal gradient of the Coriolis parameter,

*K*the horizontal total wavenumber, and

*g*≈ 10 m s

^{−2}is the acceleration of the gravity,

*H*= 10 km is the mean depth (equivalent depth) of the atmosphere, and

*f*

_{0}≈ 10

^{−4}s

^{−1}the Coriolis parameter at a reference latitude.

*K*

^{2}must have a sufficiently large negative value in order for a Rossby wave to propagate eastward. More quantitatively,

*i*is the imaginary unit. Thus, the required imaginary wavelength is compared to the Rossby deformation radius (10

^{3}km). Note that this imaginary scale is much shorter than the typical horizontal scale expected for the MJO, that is, on the order of 10

^{4}km.

This rather odd conclusion can be understood in a straightforward manner by basic knowledge of the wave theories: the nonlinear wave must be identified in the wavenumber–frequency phase space outside the range of the linear-wave dispersion (cf. Grimshaw and Iooss 2003). Thus, the above result simply suggests that the Rossby waves must be nonlinear in order to propagate eastward.

In more physical terms, the above argument suggests that the linear Rossby wave must be evanescent in order to propagate eastward. An immediate application of the general theoretical statement just above is that henceforth, an eastward-propagating Rossby wave must be nonlinear, if it is to be nonevanescence at least in some part. In the following, we are going to apply this conclusion to the wavenumbers both in latitudinal and longitudinal directions; thus, the above rather abstract conclusion from wave theories can be better understood in physical terms. Here, keep in mind that although the following analysis is strictly limited to the linear dispersion relation in Eq. (1), the conclusion always leads to a nonlinear wave based on the general theoretical principle just presented.

We first apply this argument in the latitudinal direction, because it is often argued that the results of these aforementioned high-resolution simulations (Monier et al. 2010; Hagos et al. 2011; Holloway et al. 2013) are rather trivial, simply being forced from the lateral latitudinal walls of the model channel domain typically placed at 30°–45° latitude. Once the lateral forcing is turned off, the planetary-scale MJO circulations disappear (Ray et al. 2009; Vitart and Jung 2010; Ray and Li 2013). Gustafson and Weare (2004), Ray and Zhang (2010), Ray and Li (2013), and Zhao et al. (2013) suggest that the eastward-propagating Rossby wave train arriving from higher latitudes induces the MJO.

However, the above argument suggests that any eastward-moving structure with a longitudinal wavenumber comparable to that of the MJO would only have an evanescent response away from the forcing wall with a decaying scale of 10^{3} km. So long as a linear theory is concerned, no response would be found at the equator by forcing at 30°–45° latitude. Thus, although the latitudinal sidewall forcing may help to maintain a simulated MJO-like structure, the sole possibility for obtaining a MJO-like feature centered at the equator as its response is by a nonlinearity (cf. Wedi and Smolarkiewicz 2010).

An analogous argument can also be applied in the longitudinal direction. In this case, we prescribe a typical latitudinal scale associated with the MJO (as a real wavenumber) and estimate a required longitudinal wavenumber for obtaining the MJO propagation speed by a linear Rossby wave. This estimate suggests that a linear Rossby wave consistent with the MJO must be evanescent with a decaying scale of about 10^{3} km in the longitudinal direction. Here, it is slightly odd to apply the concept of the evanescent wave to a periodic domain like in the longitudinal direction, because this concept is valid only locally, and only in the decaying direction. In the opposite direction, the evanescent wave actually grows exponentially, and at a certain point it is no longer linear; thus, nonlinearity must be invoked. So, if the MJO is to have nonevanescent structure somewhere along its propagation, it must again be a nonlinear wave.

These two applications of the result from the Rossby wave dispersion in Eq. (1) suggest that (i) the MJO must be a nonlinear wave and (ii) the MJO contains exponentially decaying tails longitudinally. As Yano and Flierl (1994) suggest, a nonlinear wave with exponential tails typically consists of a solitary structure: the main high-amplitude part, which is nonlinear, is isolated in longitude, and the wave consists of exponentially decaying linear tails away from the high-amplitude center. Thus, the most consistent manner for interpreting the MJO by dry dynamics is as a nonlinear solitary Rossby wave.

As for the inherent nonlinearity of the MJO, even a simpler but rather convincing argument can be provided: the observed order of the magnitude (0.2) of horizontal divergence relative to that of vorticity associated with MJO (Yano et al. 2009) is much larger than the values expected for linear Rossby waves of the similar scale, which are 0.03 and 0.09 for the zonal wavenumbers 1 and 2, respectively (Delayen and Yano 2009). The nonlinear nature of the MJO is also emphasized in some of the recent modeling and theoretical studies (e.g., Hagos et al. 2014; Pritchard and Bretherton 2014; Adames and Kim 2016). A scale analysis in appendix A based on the observational data analysis by Zhang and Ling (2012) also suggests strong nonlinearity of MJO.

Here, strong nonlinearity is essential: though the Rossby soliton theory (Boyd 1980) under weak nonlinearity is an attractive alternative explanation,

## 3. Modon formulation

*ψ*is the streamfunction;

*a*is the planetary radius;

*λ*the longitude; and

*F*. For facilitating an analytic progress, the most convenient choice would be to assume a linear functional, but with different forms for the inner and the exterior parts: the inner part represents an internal wave locally characterized by a real wavenumber

*q*.

For further simplification, they divide the inner and the outer regions by a circular boundary with the angular radius

We may set

Note that there are four parameters to be determined in the problem: two constants (

## 4. Numerical experiments

The solitary Rossby–vortex pair modon (Tribbia 1984a; Verkley 1984) is introduced in the last section as a theory for the MJO. The present section tests this theory numerically.

Numerical experiments were performed using a global shallow-water model in order to demonstrate the relevance of this analytic theory. Recall that the primitive equation system reduces to a shallow-water system under a separation of the variables when only the horizontal structure of the flows is in concern. Under this reduction, as detailed in appendix C, a vertical structure of the flow in concern is characterized by a mean depth of the shallow-water system (i.e., equivalent depth).

In choosing the value for the mean depth, note the equivalent barotropic nature of the total flows associated with the MJO, as most carefully elucidated by Nishi (1989), but also as seen in Fig. 2 of Yano et al. (2009), Figs. 6–8 of Holloway et al. (2013), Fig. 8a of Zhang and Ling (2012), and Fig. 2 of Žagar and Franzke (2015). All these figures show that the signals associated with the MJO has a single sign in the vertical direction throughout the troposphere with an opposite sign above the tropopause level, a typical structure of the equivalent barotropic mode. Note the importance of considering the total-flow field owing to a strongly nonlinear nature of the MJO expected here.

Figure 2 of Kasahara and Puri (1981) and Fig. 3 of Fulton and Schubert (1985) show the vertical structure of the modes under the given mean depth for the shallow-water system. As seen therein, the mode presents an equivalent barotropic structure for the mean depth of 1–10 km. In the following, we take the mean depth of 10 km as the standard case. The same sets of experiments are also performed with the mean depths of 1 and 3 km, and the results are discussed later in the section. Note that all the other parameters are set to those of Earth’s with the reference Coriolis parameter defined by *H*, or equivalently the deformation radius

For the initialization, an analytic modon solution constructed under quasigeostrophic theory is translated into a shallow-water model by using a nonlinear normal-mode balance condition (Baer and Tribbia 1977; Kasahara 1982; Tribbia 1984b). More specifically, first, the two horizontal components of the velocity are directly obtained from the streamfunction

Here, the nonlinear normal-mode balance condition not only suppresses the initial tendency for the radiation of gravity waves but also, at least partially, alleviates the wedge instability (Neven 2001) arising from the discontinuity in the gradient of the potential vorticity in the analytic modon solution.

The shallow-water model code used for the runs is also based on the Hough-function expansion (Kasahara 1977; Salby et al. 1990). All the linear computations are performed in Hough-function space. Nonlinearities are evaluated in gridpoint space and converted back to Hough-function space (i.e., transformed spectrum method) and added as an extra forcing term for the evaluations of the temporal tendencies. Linear-wave tendencies are treated analytically using exponentials in the time integration, which allows a longer time step than otherwise. Time stepping is based on a leapfrog technique with an averaging procedure every 100 time steps. Otherwise, the code contains no space–time filtering, dissipation (friction), or forcing either physically or numerically. Thus, the model investigates the pure “free” dynamics.

Computations are performed with R40 truncation in Hough-function space, being truncated at the longitudinal wavenumber

As reviewed in the last section, the modon solution contains the two free parameters: the phase speed

The results are summarized in the phase space of ^{−1}. Under this range, the solution remains stable for more than 80 days. These values compare somewhat favorably with the observed MJO propagation speeds of 3–6 m s^{−1}, especially considering the overall simplicity of the theory. Even the slower-propagating solutions last more than 20 days so long as the phase speeds remain larger than 4 m s^{−1}. This duration is still long enough for crossing the Indian Ocean. Note that obtained preferred phase speeds compare well with those obtained by the existing convectively coupled wave theories reviewed in the introduction.

The most stable size parameter of

Time evolution with the initial condition given in Fig. 2, as an example of steadily eastward propagation of the modon solution, is shown in Fig. 4.

In general, a typical numerical time evolution of a modon solution may be divided into three stages: 1) steady eastward propagation of the solution as predicted by the analytic solution; 2) the modon propagation strolls, and sometimes, begins to propagate backward (westward); and 3) breakdown of the modon solution coherency, either by wave radiation (Flierl and Haines 1994) or local nonlinearity (Neven 2001). Durations until the end of the first stage are shown in Fig. 3. Durations until the end of the last two stages are shown in Fig. 5. Here, the end of each stage is judged by eye with increments of 2.5 days. Note that the transition from stage 1 to stage 3 is relatively short compared to the overall durations. The overall agreement between these three distributions suggests that the results shown in Fig. 3 are robust.

The stability characteristics of the modon solutions obtained by Fig. 3 may be qualitatively understood in the following manner. According to Eq. (1), the modon propagates slower with the larger

Finally, to see a sensitivity of the results to the choice of the mean depth, the same series of runs is repeated for the cases with mean depths of 1 and 3 km. The obtained duration diagram in the same format as in Fig. 3 but with a mean depth of 1 km is shown in Fig. 6. Arguably, the case with the 1-km depth predicts the observed phase speed better with the preferred phase speeds reducing to 2–4 m s^{−1}, but also with a result of a reduced stability and longevity. Also, the preferred stable scale reduces to *c*_{p} = 4 m s^{−1} is shown in Fig. 7: the modon solution relatively stably and steadily propagates for 20 days over 60° longitude, then begins to propagate backward (westward), and after a few days suddenly breaks down. Both the duration and the distance are just enough for crossing the Indian Ocean. The results with the mean depth of 3 km shown in Fig. 8 are essentially intermediate between those with 1 and 10 km.

## 5. Further discussion

The present study suggests that moist convection may not be even an essential part of the dynamics of the MJO. It is shown that a simple nonlinear free Rossby wave solution can generally explain the basic characteristics of the MJO, especially the slow propagation speed, but also the observed dipolar-vortex structure, in terms of the stability of the theoretical solutions. The nonlinear nature of the MJO is already emphasized by, for example, Hagos et al. (2014), Pritchard and Bretherton (2014), and Adames and Kim (2016) in terms of moisture advection. The present study argues, in turn, that the strong nonlinearity is sufficient for explaining the essential feature of the slow eastward propagation of the MJO as well as its preferred scale. Both agreements, within a factor of 3 with the observations, are rather promising with the simplest theory adopted herein.

As discussed in section 2, the success of simulating the dry component of MJO by high-resolution models under channel configurations is likely due to the observed forcing provided at the latitudinal walls. Combination of observational evidences and modeling (Gustafson and Weare 2004; Ray and Zhang 2010; Ray and Li 2013; Zhao et al. 2013) suggest that the eastward-propagating Rossby wave train arriving from higher latitudes induces the MJO. Based on this picture, an earlier nonlinear experiment by Wedi and Smolarkiewicz (2010) explicitly adds the lateral forcing. On the other hand, the present result suggests that the overall feature of the MJO may be understood as a free nonlinear wave without any forcing from the higher latitudes to a leading-order approximation.

Moist convection is clearly another element that is totally neglected in the present study. Here, the present study does not suggest that moist convection is not important for the MJO. We rather simply argue that very basic overall features of MJO can be understood even without moist convection as a leading-order approximation to the full problem. As a whole, it is difficult to overemphasize the importance of such a much-simplified “prototype” model for the MJO: it would provide a solid basis for further theoretical studies as well as for better interpreting full global model results. It is furthermore a true test for actually improving the prediction of the MJO by global models. The basic purpose of the present study is to provide such a prototype.

In this respect, our proposal of the MJO as a “dry” nonlinear free solitary vortex pair may be compared to the well-established dry baroclinic instability theory for explaining midlatitude cyclogenesis. Both phenomena are clearly associated with extensive convective activities and cloud microphysics. Thus, both explanations of these phenomena using dry dynamics are counterintuitive. Indeed, there was a rather long time lapse for baroclinic instability theory to be fully accepted in the meteorological community. For example, the chapter on the synoptic-scale dynamics in Ludlam (1980) has no reference to baroclinic instability.

The present paper suggests that the MJO must first of all be considered from the perspective of dry planetary-scale dynamics, rather than that of convective forcing to planetary-scale motions. This suggestion is based on the recent results by high-resolution global channel-domain simulations as well as data analysis. Rather simple theoretical considerations herein further suggest the importance of the inherently nonlinear nature of the MJO dynamics. The scale analysis presented in appendix A further supports this view in a more concrete manner. A more general tropical scale analysis (Yano and Bonazzola 2009) suggests that the same could be true for the whole tropical planetary-scale circulations.

Such a new perspective may have an immediate impact on global model initialization strategy over the tropics: the current basic strategy is an initialization based on a linear equatorial wave decomposition (cf. Žagar et al. 2005). The proposed MJO–modon theory suggests a nonlinear balance initialization (Baer and Tribbia 1977; Kasahara 1982; Tribbia 1984b) as the key for a successful MJO forecast. From a perspective of the present theory, a key for the successful MJO simulation would rather reside on preparing an initial condition that is consistent with the nonlinear balance expected from the modon theory, rather than improving (or tuning) a convection parameterization.

It is commonly considered that a planetary-scale moist convective coherency is the core of the MJO in the same sense as an intensive rainband looks like a core of midlatitude synoptic-scale weather. The present study, in turn, suggests that the core of the MJO is rather a planetary-scale solitary Rossby–vortex pair in the same sense as the core of midlatitude synoptic-scale weather is baroclinic-instability waves. As a case for the latter, it clearly becomes important to analyze the potential-vorticity budget in MJO studies both based on observations and modeling. Specifically, a key for demonstrating the validity of the modon theory would be to identify the nonlinear balance [cf. Eq. (2)] in the potential-vorticity budget as a manifestation of modon dynamics. A preliminary but positive result is already found by the scale analysis presented in appendix A. The analysis in the Hough-function space (cf. Žagar and Franzke 2015) would likely to be helpful for further elucidating the nonlinear balance in a more lucid manner.

For a lucid theoretical demonstration, various higher-order effects are neglected in the present study. An obvious next step is to move to higher-order analyses in asymptotic expansion in the MJO problem. Such an analysis would be likely to identify a necessary constraint on moist convection in order to ensure a stability of the leading-order modon solution, analogous to the solvability condition in weakly nonlinear wave dynamics. Latitudinal lateral forcing may also be included as a higher-order effect. In this manner, the present theory provides a solid basis for constructing a full theoretical picture of the MJO in stepwise manner, which can furthermore be used for interpreting and improving the forecasts by numerical modeling.

Further data analyses of the MJO are also clearly required from the point of view of the potential-vorticity dynamics. Especially, the discrepancy between the vortex-pair size in snapshots, as well as predicted by the present theory and those obtained by previous composite studies, is still to be resolved. We speculate that the discrepancy is due to 1) the smoothing effect of composite by taking convection as a reference variable 2) the fact that a relatively heavy spatial filter is applied in some composite studies (e.g., Hendon and Salby 1994). However, this speculation must be verified by performing an MJO composite by taking the vorticity field as a reference variable.

Moreover, the characteristic equivalent depth (mean depth in the shallow-water system) for the “dynamical” MJO must still be identified observationally. This uncertainty leads to a major uncertainty in interpreting the results of the present study. With the standard parameter with a mean depth of 10 km, the identified preferred propagation speeds, 8–18 m s^{−1}, are too fast compared to the observation (3–6 m s^{−1}). However, simply by reducing the mean depth to 1 km, the preferred propagation speeds can easily be reduced to a range (2–4 m s^{−1}) that is closer to the observation. At the same time, the scale analysis in appendix A suggests that the propagation speed of MJO sensitively depends on diabatic heating associated with moist convection. These issues need to be investigated further.

## Acknowledgments

Continuous discussions with Akira Kasahara during this work are much appreciated. Extensive comments about the experiments by Wim T. M. Verkley and comments on the manuscript by Yanluan Lin, Annick Terpstra, and three anonymous reviewers are also much appreciated. Wenhao Dong has kindly provided Fig. 1. Piotr Smolarkiewicz provided the first author an alternative shallow-water code, which was used for the preliminary study.

## APPENDIX A

### MJO Scale Analysis

The scale analysis of the MJO is presented in this appendix in support of the main text claiming that the MJO is a strongly nonlinear free wave to a leading-order approximation. Note that with only the orders of magnitude as a concern, the rounded numbers are used for the estimates throughout this appendix.

*p*as the vertical coordinate (Holton 2004, chapter 6). Here,

**v**is the horizontal velocity;

*Q*is the diabatic heating.

We refer to the analysis by Zhang and Ling (2012, hereafter ZL12) as an observational reference. As a minor complication, ZL12 use Ertel’s full definition of the potential vorticity, which has a different dimension than the quasigeostrophic potential vorticity. To convert the values from their analysis (potential vorticity units; 1 PVU = 10^{−6} K kg^{−1} m^{2} s^{−1}) to the quasigeostrophic value (s^{−1}), we divide the former by *g* ~ 10 m s^{−2}, the typical increase of the potential temperature over the troposphere ^{−3} s^{−1}.

A typical potential-vorticity value found by ZL12 associated with the MJO, which may also be estimated by relative vorticity *ζ* to a good approximation according to ZL12, is *ζ* ~ 10 PVU ~ 10^{−5} s^{−1}. The fact that potential vorticity can be approximated well by relative vorticity suggests that the contribution of the term with the deformation radius is less important, as is the case for any flow that is close to barotropic.

Their major finding is that the local potential-vorticity tendency, ^{5} s. A characteristic time scale associated with this variability is

*υ*is the meridional velocity. Assuming a typical wind velocity to be

*L*. Rearranging this order-balance condition, we obtain an estimate of the horizontal velocity as

In this manner, we demonstrate by scale analysis that the nonlinear self-advection balance is the key in the MJO dynamics, as we also assume is the main ingredient of the modon dynamics in the main text. An important corollary from the scale analysis here is that to a leading-order approximation, the MJO may be considered stationary (i.e.,

## APPENDIX B

### Full Modon Solution

*H*(cf. Abramowitz and Stegun 1965, chapter 15) by

## APPENDIX C

### Derivation of the Shallow-Water System from the Primitive Equation System

The primitive equation system under the linearization reduces to the shallow-water system by applying a separation of the variables to the former (Kasahara and Puri 1981; Fulton and Schubert 1985).

*p*as the vertical coordinate [e.g., Eq. (1) of Yano et al. 2009]. Here,

*ξ*is the geopotential,

**k**is the normal vector in the vertical direction,

*R*the gas constant,

*T*is the temperature also defined as a deviation from a reference state,

*ω*is the vertical velocity (the pressure velocity), and

*H*and a reference pressure

*W*and

Note that Eqs. (C3) and (C10) constitute the linear shallow-water system with mean depth *H*. Conceptually, the nonlinear advection terms may be further added and consider it as an analog model for full atmospheric circulation under a prescribed vertical structure. Strictly speaking, such a shallow-water analog model is justified only if the vertical structure of the flow is reasonably homogeneous. The equivalent barotropic structure assumed in the main text is consistent with this requirement. A special case for the above is when a flow is well confined to a shallow layer. For example, the ocean currents are often strongly confined to a surface layer; thus, the nonlinear shallow-water model can extensively be used for the basic investigations of the ocean dynamics (cf. Charney and Flierl 1981). For similar reasons, the nonlinear shallow-water analog modeling has been applied to the studies of the tropical-atmospheric dynamics (e.g., Yano et al. 1995).

The vertical structure under the mean depth *H* is defined by Eq. (C9b) separately. Vertical structure solutions

## REFERENCES

Abramowitz, M., and I. A. Stegun, 1965:

*Handbook of Mathematical Functions*. Dover, 1046 pp.Adames, Á. F., and D. Kim, 2016: The MJO as a dispersive, convectively coupled moisture wave: Theory and observations.

,*J. Atmos. Sci.***73**, 913–941, doi:10.1175/JAS-D-15-0170.1.Adames, Á. F., J. Patoux, and R. C. Foster, 2014: The contribution of extratropical waves to the MJO wind field.

,*J. Atmos. Sci.***71**, 155–176, doi:10.1175/JAS-D-13-084.1.Arnold, N. P., and D. A. Randall, 2015: Global-scale convective aggregation: Implications for the Madden-Julian Oscillation.

,*J. Adv. Model. Earth Syst.***7**, 1499–1518, doi:10.1002/2015MS000498.Baer, F., and J. J. Tribbia, 1977: On complete filtering of gravity modes through nonlinear initialization.

,*Mon. Wea. Rev.***105**, 1536–1539, doi:10.1175/1520-0493(1977)105<1536:OCFOGM>2.0.CO;2.Benedict, J. J., and D. A. Randall, 2009: Structure of the Madden–Julian oscillation in the superparameterized CAM.

,*J. Atmos. Sci.***66**, 3277–3296, doi:10.1175/2009JAS3030.1.Boyd, J. P., 1980: Equatorial solitary waves. Part I: Rossby solitons.

,*J. Phys. Oceanogr.***10**, 1699–1717, doi:10.1175/1520-0485(1980)010<1699:ESWPIR>2.0.CO;2.Charney, J. G., 1963: A note on large-scale motions in the tropics.

,*J. Atmos. Sci.***20**, 607–609, doi:10.1175/1520-0469(1963)020<0607:ANOLSM>2.0.CO;2.Charney, J. G., and F. R. Flierl, 1981: Oceanic analogues in large-scale atmospheric motions.

*Evolution of Physical Oceanography*, B. Warren and C. Wunsch, Eds., MIT Press, 502–546.Delayen, K., and J.-I. Yano, 2009: Is asymptotic non-divergence of the large-scale tropical atmosphere consistent with equatorial wave theories?

,*Tellus***61A**, 491–497, doi:10.1111/j.1600-0870.2009.00404.x.Emanuel, K. A., 1987: An air–sea interaction model of intraseasonal oscillations in the tropics.

,*J. Atmos. Sci.***44**, 2324–2340, doi:10.1175/1520-0469(1987)044<2324:AASIMO>2.0.CO;2.Flierl, G. R., and K. Haines, 1994: The decay of modons due to Rossby wave radiation.

,*Phys. Fluids***6**, 3487–3497, doi:10.1063/1.868405.Fuchs, Ž., and D. J. Raymond, 2007: A simple, vertically resolved model of tropical disturbances with a humidity closure.

,*Tellus***59A**, 344–354, doi:10.1111/j.1600-0870.2007.00230.x.Fulton, S. R., and W. H. Schubert, 1985: Vertical normal mode transforms: Theory and application.

,*Mon. Wea. Rev.***113**, 647–658, doi:10.1175/1520-0493(1985)113<0647:VNMTTA>2.0.CO;2.Grimshaw, R., and G. Iooss, 2003: Solitary waves of a coupled Korteweg–de Vries system.

,*Math. Comput. Simul.***62**, 31–40, doi:10.1016/S0378-4754(02)00189-1.Gustafson, W. I., Jr., and B. C. Weare, 2004: MM5 modeling of the Madden–Julian oscillation in the Indian West Pacific Oceans: Model description and control run results.

,*J. Climate***17**, 1320–1337, doi:10.1175/1520-0442(2004)017<1320:MMOTMO>2.0.CO;2.Hagos, S., R. L. Leung, and J. Dudhia, 2011: Thermodynamics of the Madden–Julian oscillation in a regional model with constrained moisture.

,*J. Atmos. Sci.***68**, 1974–1989, doi:10.1175/2011JAS3592.1.Hagos, S., Z. Feng, K. Landu, and C. N. Long, 2014: Advection, moistening, and shallow-to-deep convection transitions during the initiation and propagation of Madden-Julian Oscillation.

,*J. Adv. Model. Earth Syst.***6**, 938–949, doi:10.1002/2014MS000335.Hartmann, D. L., and H. H. Hendon, 2007: Resolving an atmospheric enigma.

,*Science***318**, 1731, doi:10.1126/science.1152502.Hayashi, Y., 1970: A theory of large-scale equatorial waves generated by condensation heat and accelerating the zonal wind.

,*J. Meteor. Soc. Japan***48**, 140–160, doi:10.2151/jmsj1965.48.2_140.Hendon, H. H., and M. L. Salby, 1994: The life cycle of the Madden–Julian oscillation.

,*J. Atmos. Sci.***51**, 2225–2237, doi:10.1175/1520-0469(1994)051<2225:TLCOTM>2.0.CO;2.Holloway, C. E., S. J. Woolnough, and G. M. S. Lister, 2013: The effects of explicit versus parameterized convection on the MJO in a large-domain high-resolution tropical case study. Part I: Characterization of large-scale organization and propagation.

,*J. Atmos. Sci.***70**, 1342–1369, doi:10.1175/JAS-D-12-0227.1.Holton, J. R., 2004:

*An Introduction to Dynamic Meteorology*. 4th ed. International Geophysical Series, Vol. 88, Academic Press, 507 pp.Jiang, X., and Coauthors, 2015: Vertical structure and physical processes of the Madden–Julian oscillation: Exploring key model physics in climate simulations.

,*J. Geophys. Res. Atmos.***120**, 4718–4748, doi:10.1002/2015JD023196.Kasahara, A., 1977: Numerical integration of the global barotropic primitive equations with Hough harmonic expansions.

,*J. Atmos. Sci.***34**, 687–701, doi:10.1175/1520-0469(1977)034<0687:NIOTGB>2.0.CO;2.Kasahara, A., 1982: Nonlinear normal mode initialization and the bounded derivative method.

,*Rev. Geophys.***20**, 385–397, doi:10.1029/RG020i003p00385.Kasahara, A., and G. W. Platzman, 1963: Interaction of a hurricane with the steering flow and its effect upon the hurricane trajectory.

,*Tellus***15**, 321–335, doi:10.3402/tellusa.v15i4.8863.Kasahara, A., and K. Puri, 1981: Spectral representation of three-dimensional global data by expansion in normal mode functions.

,*Mon. Wea. Rev.***109**, 37–51, doi:10.1175/1520-0493(1981)109<0037:SROTDG>2.0.CO;2.Kiladis, G. N., K. H. Straub, and P. T. Haertel, 2005: Zonal and vertical structure of the Madden–Julian oscillation.

,*J. Atmos. Sci.***62**, 2790–2809, doi:10.1175/JAS3520.1.Lin, H., G. Brunet, and J. Derome, 2007: Instraseasonal variability in a dry atmospheric model.

,*J. Atmos. Sci.***64**, 2422–2441, doi:10.1175/JAS3955.1.Lindzen, R. S., 1974: Wave-CISK in the tropics.

,*J. Atmos. Sci.***31**, 156–179, doi:10.1175/1520-0469(1974)031<0156:WCITT>2.0.CO;2.López Carrillo, C., and D. J. Raymond, 2005: Moisture tendency equation in a tropical atmosphere.

,*J. Atmos. Sci.***62**, 1601–1613, doi:10.1175/JAS3424.1.Ludlam, F. H., 1980:

*Clouds and Storms*. Penn State University Press, 488 pp.Madden, R. A., and P. R. Julian, 1972: Description of global-scale circulation cells in the tropics with a 40–50 day period.

,*J. Atmos. Sci.***29**, 1464–1469, doi:10.1175/1520-0469(1972)029<1109:DOGSCC>2.0.CO;2.Majda, A. J., and S. N. Stechmann, 2009: The skeleton of tropical intraseasonal oscillations.

,*Proc. Natl. Acad. Sci. USA***106**, 8417–8422, doi:10.1073/pnas.0903367106.Matsuno, T., 1966: Quasi-geostrophic motions in the equatorial area.

,*J. Meteor. Soc. Japan***44**, 25–43, doi:10.2151/jmsj1965.44.1_25.Miura, H., M. Satoh, T. Nasuno, A. T. Noda, and K. Oouchi, 2007: A Madden–Julian oscillation event realistically simulated by a global cloud-resolving model.

,*Science***318**, 1763–1765, doi:10.1126/science.1148443.Moncrieff, M. W., 2004: Analytic representation of the large-scale organization of tropical convection.

,*J. Atmos. Sci.***61**, 1521–1538, doi:10.1175/1520-0469(2004)061<1521:AROTLO>2.0.CO;2.Monier, E., B. C. Weare, and W. I. Gustafson Jr., 2010: The Madden–Julian oscillation wind–convection coupling and the role of moisture processes in the MM5 model.

,*Climate Dyn.***35**, 435–447, doi:10.1007/s00382-009-0626-4.Neven, E. C., 2001: Linear stability of modons on a sphere.

,*J. Atmos. Sci.***58**, 2280–2305, doi:10.1175/1520-0469(2001)058<2280:LSOMOA>2.0.CO;2.Nishi, N., 1989: Observational study on the 30-60 day variations in the geopotential and temperature fields in the equatorial region.

,*J. Meteor. Soc. Japan***67**, 187–203, doi:10.2151/jmsj1965.67.2_187.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, doi:10.1175/JAS-D-13-0119.1.Ray, P., and C. Zhang, 2010: A case study of the mechanics of extratropical influence on the initiation of the Madden–Julian oscillation.

,*J. Atmos. Sci.***67**, 515–528, doi:10.1175/2009JAS3059.1.Ray, P., and T. Li, 2013: Relative roles of circumnavigating waves and extratropics on the MJO and its relationship with the mean state.

,*J. Atmos. Sci.***70**, 876–893, doi:10.1175/JAS-D-12-0153.1.Ray, P., C. Zhang, J. Dudhia, and S. S. Chen, 2009: A numerical case study on the initiation of the Madden–Julian oscillation.

,*J. Atmos. Sci.***66**, 310–331, doi:10.1175/2008JAS2701.1.Ray, P., C. Zhang, M. W. Moncrieff, J. Dudhia, J. M. Caron, L. R. Leung, and C. Bruyere, 2011: Role of the atmospheric mean state on the initiation of the Madden–Julian oscillation in a tropical channel model.

,*Climate Dyn.***36**, 161–184, doi:10.1007/s00382-010-0859-2.Raymond, D. J., and Ž. Fuchs, 2007: Convectively coupled gravity and moisture modes in a simple atmospheric model.

,*Tellus***59A**, 627–640, doi: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, doi:10.1175/2008JCLI2739.1.Salby, M. L., R. R. Garcia, D. O’Sullivan, and J. J. Tribbia, 1990: Global transport calculations with an equivalent barotropic system.

,*J. Atmos. Sci.***47**, 188–214, doi:10.1175/1520-0469(1990)047<0188:GTCWAE>2.0.CO;2.Straus, D. M., and R. S. Lindzen, 2000: Planetary-scale baroclinic instability and the MJO.

,*J. Atmos. Sci.***57**, 3609–3626, doi:10.1175/1520-0469(2000)057<3609:PSBIAT>2.0.CO;2.Thual, S., and A. J. Majda, 2016: A skeleton model for the MJO with refined vertical structure.

,*Climate Dyn.***46**, 2773–2786, doi:10.1007/s00382-015-2731-x.Tribbia, J. J., 1984a: Modons in spherical geometry.

,*Geophys. Astrophys. Fluid Dyn.***30**, 131–168, doi:10.1080/03091928408210080.Tribbia, J. J., 1984b: A simple scheme for high-order nonlinear normal mode initialization.

,*Mon. Wea. Rev.***112**, 278–284, doi:10.1175/1520-0493(1984)112<0278:ASSFHO>2.0.CO;2.Verkley, W. T. M., 1984: The construction of barotropic modons on a sphere.

,*J. Atmos. Sci.***41**, 2492–2504, doi:10.1175/1520-0469(1984)041<2492:TCOBMO>2.0.CO;2.Vitart, F., and T. Jung, 2010: Impact of the northern hemisphere extratropics on the skill in predicting the Madden Julian Oscillation.

,*Geophys. Res. Lett.***37**, M23805, doi:10.1029/2010GL045465.Wang, B., 2005: Theory.

*Intraseasonal Variability in the Atmosphere–Ocean System*, W. K. Lau and D. E. Waliser, Eds., Springer, 307–360.Wedi, N. P., and P. K. Smolarkiewicz, 2010: A nonlinear perspective on the dynamics of the MJO: Idealized large-eddy simulations.

,*J. Atmos. Sci.***67**, 1202–1217, doi:10.1175/2009JAS3160.1.Yanai, M., 1964: Formation of tropical cyclones.

,*Rev. Geophys.***2**, 367–414, doi:10.1029/RG002i002p00367.Yanai, M., B. Chen, and W.-W. Tung, 2000: The Madden–Julian oscillation observed during the TOGA COARE IOP: Global view.

,*J. Atmos. Sci.***57**, 2374–2396, doi:10.1175/1520-0469(2000)057<2374:TMJOOD>2.0.CO;2.Yang, D., and A. P. Ingersoll, 2011: Testing the hypothesis that the MJO is a mixed Rossby–gravity wave packet.

,*J. Atmos. Sci.***68**, 226–239, doi:10.1175/2010JAS3563.1.Yang, D., and A. P. Ingersoll, 2013: Triggered convection, gravity waves, and the MJO: A shallow-water model.

,*J. Atmos. Sci.***70**, 2576–2486, doi:10.1175/JAS-D-12-0255.1.Yano, J.-I., and K. A. Emanuel, 1991: An improved model of the equatorial troposphere and its coupling with the stratosphere.

,*J. Atmos. Sci.***48**, 377–389, doi:10.1175/1520-0469(1991)048<0377:AIMOTE>2.0.CO;2.Yano, J.-I., and G. R. Flierl, 1994: Jupiter’s Great Red Spot: Compacting conditions, stabilities.

,*Ann. Geophys.***12**, 1–18, doi:10.1007/s00585-994-0001-z.Yano, J.-I., and M. Bonazzola, 2009: Scale analysis for the large-scale tropical atmospheric dynamics.

,*J. Atmos. Sci.***66**, 159–172, doi:10.1175/2008JAS2687.1.Yano, J.-I., J. C. McWilliams, M. W. Moncrieff, and K. A. Emanuel, 1995: Hierarchical tropical cloud systems in an analog shallow-water model.

,*J. Atmos. Sci.***52**, 1723–1742, doi:10.1175/1520-0469(1995)052<1723:HTCSIA>2.0.CO;2.Yano, J.-I., S. Mulet, and M. Bonazzola, 2009: Tropical large-scale circulations: Asymptotically non-divergent?

,*Tellus***61A**, 417–427, doi:10.1111/j.1600-0870.2009.00397.x.Žagar, N., and C. L. E. Franzke, 2015: Systematic decomposition of the Madden-Julian Oscillation into balanced and inertio-gravity composition.

,*Geophys. Res. Lett.***42**, 6829–6835, doi:10.1002/2015GL065130.Žagar, N., E. Andersson, and M. Fisher, 2005: Balanced tropical data assimilation based on a study of equatorial waves in ECMWF short-range forecast errors.

,*Quart. J. Roy. Meteor. Soc.***131**, 987–1011, doi:10.1256/qj.04.54.Zhang, C., 2005: Madden–Julian oscillation.

,*Rev. Geophys.***43**, RG2003, doi:10.1029/2004RG000158.Zhang, C., and J. Ling, 2012: Potential vorticity of the Madden–Julian oscillation.

,*J. Atmos. Sci.***69**, 65–78, doi:10.1175/JAS-D-11-081.1.Zhao, C., T. Li, and T. Zhou, 2013: Precursor signals and processes associated with MJO initiation over the tropical Indian Ocean.

,*J. Climate***26**, 291–307, doi:10.1175/JCLI-D-12-00113.1.