The climate sensitivity of the Madden–Julian oscillation (MJO) is measured across a broad range of temperatures (1°–35°C) using a convection-permitting global climate model with homogenous sea surface temperatures. An MJO-like signal is found to be resilient in all simulations. These results are used to investigate two ideas related to the modern “moisture mode” view of MJO dynamics. The first hypothesis is that the MJO has dynamics analogous to a form of radiative convective self-aggregation in which longwave energy maintenance mechanisms shut down for SST ≪ 25°C. Inconsistent with this hypothesis, the explicitly simulated MJO survives cooling and retains leading moist static energy (MSE) budget terms associated with longwave destabilization even at SST < 10°C. Thus, if the MJO is a form of longwave-assisted self-aggregation, it is not one that is temperature critical, as is observed in some cases of radiative–convective equilibrium (RCE) self-aggregation. The second hypothesis is that the MJO is propagated by horizontal advection of column MSE. Inconsistent with this view, the simulated MJO survives reversal of meridional moisture gradients in the basic state and a striking role for horizontal MSE advection in its propagation energy budget cannot be detected. Rather, its eastward motion is balanced by vertical MSE advection reminiscent of gravity or Kelvin wave dynamics. These findings could suggest a tight relation between the MJO and classic equatorial waves, which would tend to challenge moisture mode views of MJO dynamics that assume horizontal moisture advection as the MJO’s propagator. The simulation suite provides new opportunities for testing predictions from MJO theory across a broad climate regime.
The Madden–Julian oscillation (MJO) is a slow eastward-moving mode that is the most energetic intraseasonal cycle in the tropical atmosphere. Although critical to regional weather and climate, its physics remain in debate (Zhang 2005; Raymond and Fuchs 2009; Majda and Stechmann 2009; Yang and Ingersoll 2013, 2014). Until recently explicit MJO simulation has been difficult but new advances in modern multiscale atmospheric simulation (Grabowski and Smolarkiewicz 1999; Randall et al. 2003) have enabled realistic MJO simulation with minimal approximations (Khairoutdinov et al. 2005, 2008; Benedict and Randall 2009). Combined with new computational acceleration techniques (Pritchard et al. 2014; Jones et al. 2015), this technology now allows for ideas about MJO physics to be efficiently stress tested explicitly in extreme settings.
One idea worth exploring is the “moisture mode” view that the MJO is fundamentally energized and propagated by thermodynamic feedbacks acting to regulate the tropical moisture field (Raymond and Fuchs 2009; Sobel and Maloney 2012)—in turn a good proxy for the column moist static energy (MSE) under the weak temperature gradient approximation to tropical dynamics (Sobel et al. 2001).
Two specific feedbacks involving longwave radiation and horizontal advection are commonly emphasized as critical to MJO maintenance and propagation from this view. The first idea is that radiation from high clouds helps energize the MJO. The second idea is that the eastward motion of the MJO is supported through the advection of MSE by its horizontal circulation anomalies. Column MSE budget analysis of the MJO’s horizontal structure has supported both of these views, as applied to reanalyses, satellite data, and MJO-permitting global models.
We explore two key hypotheses to this view. The first hypothesis H1 is that the existence of the MJO depends on a radiative–convective instability—a positive feedback between atmospheric thermal radiation and water vapor anomalies. This is motivated by the fact that longwave heating is a leading-order diabatic balance of the MJO variance column energy budget in MJO-permitting models (Maloney 2009; Andersen and Kuang 2012, hereafter AK12; Kim et al. 2015) and in reanalysis (Kiranmayi and Maloney 2011) and satellite (Ma and Kuang 2011) data. The importance of a longwave feedback mechanism to MJO dynamics is especially motivated by the fact that in MJO-permitting global models, mechanism denial experiments that spatially homogenize longwave heating anomalies can destroy the simulated MJO (Landu and Maloney 2011; Andersen 2012).
The first hypothesis is relevant to an analogy A1 that has been made between the MJO and convective self-aggregation (spontaneous conglomeration of convective updrafts into large, organized clusters) in the idealized limit of radiative–convective equilibrium (RCE). Here RCE refers to a state in which convective heating balances radiative cooling, neglecting large-scale lateral energy transport into the domain of interest. The analogy between RCE self-aggregation and the MJO is that cloud resolving simulations under RCE can produce MJO-like aggregation and correlations between fields such as MSE, precipitation, surface flux, and longwave radiation anomalies (Bretherton et al. 2005).
It is logical to wonder under H1 (if A1 is valid) whether the MJO might cease to exist below a critical temperature threshold. This is based on the fact that the self-aggregation longwave feedback mechanism has the capacity to shut down under a critical temperature threshold somewhere below 20°C (Khairoutdinov and Emanuel 2010; Emanuel et al. 2014; Wing and Emanuel 2014). The existence and magnitude of a critical temperature threshold for RCE self-aggregation and its relevance to MJO dynamics are in debate (Emanuel et al. 2014; Wing and Cronin 2016; Abbot 2014), motivating the need for extreme climate sensitivity testing of MJO dynamics.
Our second hypothesis H2 is that MJO eddies help it move east by mixing moisture horizontally across sharp meridional gradients. The second hypothesis is motivated by consistent diagnostic phase relationships seen in reanalyses (Benedict and Randall 2007; Kiranmayi and Maloney 2011) and global climate models (Maloney 2009; Maloney et al. 2010; AK12). It is especially supported by MJO phase speed sensitivities found in controlled climate model experiments that have interfered with the horizontal moisture advection process either through the basic state (Maloney et al. 2010) or directly (Pritchard and Bretherton 2014). The idea of H2 can be encapsulated diffusively in a semiempirical analytic MJO model (Sobel and Maloney 2012). One way to think about this advective propagation mechanism is that the observed MJO is dried on the trailing (western) flank of its moist core by the meridional advection of relatively dry free-tropospheric subtropical air by strong MJO Rossby gyres, and that this western drying is consistent with the MJO’s eastward motion if moisture is fundamental to its dynamics.
Whether this process is actually fundamental to the eastward travel of the MJO remains in debate. In their analysis of the reanalysis-derived moisture budget, Hsu and Li (2012) argued for equatorial-wave-driven vertical moisture advection in the PBL as a more prominent controlling factor for eastward propagation, in contrast with the interpretation of Kiranmayi and Maloney (2011), which emphasized free-tropospheric moisture variations. Yang and Ingersoll (2013) argue that the MJO’s eastward propagation can be explained by the zonal asymmetry between eastward and westward inertia–gravity waves. Some GCM mechanism denial studies even suggest a resilience of MJO signals to denying the free-tropospheric moisture advection process through basic-state interference and nudging (Hsu et al. 2014).
As a new test of both H1 and H2, and the RCE analogy (i.e., A1), we use a global climate model with a convincing intrinsic MJO signal and apply it in an idealized setting that has not been previously explored to our knowledge. The first idealization is to vary temperature across an extreme range—from 35° down to 1°C. The expectation from H1, if RCE temperature criticality is relevant, is that the MJO should cease to exist below a critical temperature below 20°C. The second idealization is to force the model unusually with globally uniform sea surface temperatures (SSTs) to achieve an exotic basic state that can be exploited to test H2. For our model we will show this produces a planet with a background near-tropical meridional MSE gradient that has the opposite sign to what is observed in real-world settings, with an MSE minimum (rather than a maximum) on the equator. The expectation in this case from H2 is that an associated reversal of the MSE mixing by MJO-related Rossby gyres should interfere with the ability of the eastward-propagating MJO signal to exist at a realistic phase speed. As we will demonstrate, the results on this front are surprising.
We apply the Superparameterized Community Atmosphere Model, version 3.0 (SPCAM3.0), a proven multiscale global atmospheric model that avoids approximating moist turbulence by explicitly resolving it in >10 000 nested cloud system resolving subdomains. SPCAM3.0 is known to produce an unusually realistic MJO signal while making minimal assumptions about deep convection, in both real-geography and idealized settings (Benedict and Randall 2009, 2011; AK12; Arnold et al. 2013; Pritchard et al. 2014; Pritchard and Bretherton 2014). As such it is a useful tool for MJO hypothesis stress testing.
For computational efficiency we use an eightfold accelerated version of SPCAM3.0, using fourfold acceleration from cloud resolving model (CRM) extent reduction (Pritchard et al. 2014), and an additional twofold acceleration using a new time acceleration algorithm (Jones et al. 2015), neither of which significantly impacts SPCAM3.0’s intrinsic MJO. This enables a wide ensemble of experiments.
Eight simulations are run with sea surface temperatures spanning an extreme range including many values far below the critical temperature threshold for RCE self-aggregation (1°, 5°, 10°, 15°, 20°, 25°, 30°, and 35°C). The simulations are idealized by removing continents, turning off the sun (perpetual night) and—through homogeneous SSTs—by muting meridional temperature gradients, which deemphasizes extratropical modes of variability and initiates an unusual basic-state climate useful for hypothesis testing. The simulations are analyzed for a 6-yr period following a 2-yr spinup, except for the two warmest cases, which are analyzed for a 5-yr period following a 3-month spinup.
The basic state in the experiments is summarized in Figs. 1–4. Water vapor and rainfall increase with in-line expectations for extreme climate change (Fig. 1a). An exotic global mass circulation results from using idealized globally homogeneous SSTs in aquaplanet SPCAM3.0 across the range of cold climates examined. For SST < 25°C, each simulated climate exhibits a double-ITCZ-like precipitation structure (Fig. 1b) associated with two circulation cells consisting of ascending, moist air adjacent to the equator near 15°S and 15°N (Fig. 2). For these coldest climates, associated with the double-ITCZ-like circulation anchor is a subsidence structure directly over the equator (Fig. 2), such that relatively dry and MSE-poor tropical air is located on the equator.
For the purpose of testing H2, the key property of this idealized circulation is that for SST < 25°C it happens to reverse the sign of the meridional gradient of basic-state MSE near the equator relative to what is observed in the real climate (Fig. 1b). Following H2, this should in turn reverse the ability of MJO eddies to sustain MJO eastward propagation through their advection of basic-state MSE in our coldest simulated climates. That is, for an MJO to exist in these conditions its trailing gyres would have to advect in relatively moist air on the MJO’s western flank, which would become inconsistent with eastward propagation if horizontal advection of moisture is fundamental to it.
As an aside, we note that exotic basic-state general circulations similar to this are not unusual and have been seen by others in GCMs when using globally uniform SSTs (Sumi 1992; Kirtman and Schneider 2000; Chao 2000; Chao and Chen 2004). They are possible because despite the fact that removing sharp extratropical gradients of background MSE drastically damps the overall mass circulation, throttles baroclinic instability, and reduces eddy activity, the existence of rotation alone is nonetheless sufficient to set up a weak basic-state circulation (Chao and Chen 2004). Furthermore, internal variability is such that surface homogenization of SSTs does not ensure perfectly flat isotherms throughout the atmosphere, though it comes close (Fig. 3).
In the coldest climate simulations, the basic-state zonal-mean wind on the equator is very weak and mostly easterly throughout the troposphere but becomes stronger with increasing vertical shear as SST is increased from 1° to 35°C (Fig. 4). Strong vertical shear including super-rotating zonal wind in the upper troposphere is only observed for the warmest cases.
a. Existence of an MJO signal
MJO signal analysis is performed for daily mean anomalies of key state variables (zonal winds, precipitation, and outgoing longwave radiation) averaged within a narrow band directly on the equator, from 5°S to 5°N. We use standard methods (Hayashi 1971) to discriminate eastward- versus westward-propagating tropical signals in observations following the implementation by Wheeler and Kiladis (1999, hereafter WK99), beginning from daily mean anomalies. The only deviation from the procedure of WK99 is that we use a 250-day windowing interval, with overlap of 150 days, in place of the standard nonoverlapping 96-day window. This adjustment accounts for the fact that raw signal analysis reveals the eastward-moving anomaly propagation to be slower than the real-world MJO signal in some of our coldest simulations (Fig. 5).
Figure 6 shows the resulting wavenumber–frequency analysis for the 850-hPa zonal wind U850 anomalies. Consistent with the existence of an MJO signal, in each simulated climate a distinct mode of eastward-moving (positive zonal wavenumber) equatorial wind variance is found maximizing at time scales longer than 30 days. In some climates (5°, 20°, 25°, 30°, and 35°C) it is the dominant mode of tropical variability. In other climates (1°, 10°, and 15°C) it is the second most dominant mode, rivaled by a westward-moving intraseasonal mode of equal or greater magnitude. Analysis of raw Hovmöller diagrams of 850 hPa (Fig. 5) confirms that in the cold climates these distinct variance maxima are not the result of a standing wave, as can occur using the wavenumber–frequency decomposition methodology (Zhang and Hendon 1997). Rather, unfiltered fields tend to confirm that the model produces a superposition of eastward- and westward-moving modes.
Is the distinct slow eastward-moving mode MJO-like in the cold-climate simulations? This demands multivariate analysis. Figure 7a illustrates conceptually the expected archetypal multivariate phase progression signature of the observed MJO in the time–longitude plane: an eastward-propagating negative (easterly) 850-hPa zonal wind anomaly is followed by three closely linked thermodynamic variations (positive vapor anomaly, positive precipitation anomaly, and negative outgoing longwave radiation anomaly), which are in turned followed by a positive (westerly) 850-hPa anomaly. Figure 7a demonstrates how this archetypal MJO phase relationship can be concisely summarized visually using a smooth color gradation pattern defined across the multiple MJO-modulated state variables involved.
We exploit this new visualization strategy to test whether an MJO-like multivariate phase relationship (Fig. 7a) exists in our simulations in association with the eastward-moving variance mode identified in Fig. 6. Figure 7b visualizes plus/minus one standard deviation of spectrally filtered fields in the hypothesized MJO band (defined initially as zonal wavenumber k = 1–10 and time scale τ = 20–250 days in order to encompass both the warm- and cold-climate intraseasonal spectral features in Fig. 6). Consistent with the existence of an MJO-like signal, most of the long-lived zonal wind anomalies that constitute the variance node in Fig. 6 are linked to structures that resemble the MJO’s archetypal phase progression of multiple dynamic and thermodynamic state variables. The MJO-associated thermodynamic and dynamic anomalies propagate eastward at the same speed.
Figures 8 and 9 investigate the three-dimensional dynamic flow structure of the slow eastward-moving signal. Time series of velocity potential and streamfunction on two pressure levels are lag regressed in horizontal space against a local MJO precipitation index time series (defined as the anomaly of latitudinally averaged precipitation near the equator, filtered spectrally to an MJO band). This process is repeated at many reference longitudes, which leverages zonal symmetry for statistical clarity, following the methods of AK12. To avoid sensitivity to assumed parameters in this method, we vary the latitudinal extent Φ and the spectral bandpass filter settings (maximum zonal wavelength kmax, minimum and maximum time scales τ1 and τ2) used to define the precipitation index variable over a 16-member permutation ensemble (Φ = 5°, 10°; kmax = 10, 15; τ1 = 20, 40 days; τ2 = 200, 300 days) and analyze statistically significant components of the ensemble mean.
The results confirm that the slow eastward-moving mode resembles an MJO with reduced zonal extent. In all climates, the familiar MJO 3D flow structure is clear, consisting of a central velocity potential vertical dipole anomaly in phase with precipitation (Fig. 8) and straddled horizontally by vertically antiphase quadrupole streamfunction anomalies (Fig. 9). The essence of this MJO-like 3D flow pattern survives cooling to temperatures far below 20°–25°C. In the cold climates the zonal scale of the mode is reduced. That the eastward-moving mode is not a Doppler-shifted Rossby wave, which could also have quadrupole dynamic flow anomalies, is confirmed by offline analysis of the mass-weighted barotropic flow, which is insufficient to sustain this (not shown), as also indicated by the mean zonal flow in Fig. 4. Rather, the mode appears to propagate to the east relative to the mean flow as occurs for the MJO.
In summary, because of the above evidence, we conclude that an MJO-like signal exists in each of our uniform-SST aquaplanet simulations, even when SST < 20°C. Figure 10a shows the amplification of the MJO with SST, which scales to first order with mean rainfall across the climates. That is, when normalized by climatological rainfall (Fig. 10b) the relative amplitude of the MJO is remarkably climate invariant.
b. Column moist static energy analysis
The existence of a cold MJO coexisting with reversed basic-state meridional moisture gradients in our simulations for SST < 25°C is inconsistent with expectations from H2, and, thus, could be viewed as challenging moisture mode views of MJO dynamics that emphasize horizontal advection as the key propagation mechanism.
We, therefore, analyze the column MSE budget of the simulated MJO—a diagnostic that has been a cornerstone of the moisture mode interpretation of the MJO (Raymond et al. 2009; Maloney 2009). This is by definition a purely diagnostic framework that highlights correlations between the horizontal structure of the MJO’s MSE anomaly versus terms contributing to its tendency. Such relationships can be viewed as suggestive of the MJO’s fundamental dynamics, though in and of themselves cannot offer causative proof of the processes that control MJO amplitude and phase speed.
To isolate the column MSE budget perturbations associated with the slow eastward-moving rainfall mode, we follow the established methods of AK12. That is, a filtered, latitudinally averaged precipitation index time series is used as a reference against which unfiltered terms in the column MSE budget are lag regressed spatially. By leveraging zonal symmetry, many reference longitudes can be used to enhance statistical significance of the resulting MJO-regression slopes. The regression slopes are in turn scaled by the zonal-mean standard deviation of the MJO index to yield “MJO regressed” horizontal maps of each term in the column MSE budget in physical energy units (see Figs. S1–S4 in the supplemental material). Finally, summary statistics are calculated for each of these terms by projecting their distinct MJO-regression map onto the map of MJO-regressed MSE, and onto the map of MJO-regressed MSE tendency. Thus, a single compact measure helps discriminate and contrast MSE budget terms contributing to MJO “maintenance” (in phase with MSE) versus “propagation” (in phase with MSE tendency). In practice, the method can be sensitive to several parameter choices such as Φ and the spectral filter settings used to define the MJO index variable (kmax, τ1, and τ2), and the latitudinal extent over which projections are accumulated Φp. Therefore, we repeat this calculation 48 times, varying each of these parameters over a 48-member permutation ensemble (Φ = 5°, 10°; kmax = 10, 15; τ1 = 20, 40 days; τ2 = 200, 300 days; Φp = 5°, 10°, 15°). The resulting MJO MSE budget summary statistics are shown in Fig. 11, and the underlying full MJO-regression maps of column MSE and key tendency contributions are available in Figs. S1–S4 in the supplemental material. We note that the horizontal pattern of precipitation and outgoing longwave radiation anomalies associated with the signal that we have isolated closely resemble the column MSE anomaly in Fig. S1 in the supplemental material.
1) Propagation dynamics
Horizontal MSE advection cannot explain the eastward propagation of the MJO-like mode simulated by SPCAM3.0 under constant SST boundary conditions (Fig. 11a). That is, projecting MJO-regressed maps of horizontal MSE advection does not yield leading positive projection values relative to other terms in the MSE budget, as is observed for the MJO in reanalysis and as occurs in real-geography SPCAM3.0 simulations with warm SSTs. Rather, the propagation projection of horizontal MSE advection is difficult to distinguish from zero within the parameter uncertainty of the regression–projection method. Consistently, the horizontal structure of the MJO-regression map of horizontal column MSE advection does not imply a clear propagating role (see Fig. S4 in the supplemental material). This tends to argue against H2 that horizontal moisture advection plays an important role in controlling the eastward propagation of the MJO. Rather, the propagation of the MJO in our simulations is balanced primarily through vertical MSE advection. This may suggest that gravity acts as a restoring force and leads to the eastward propagation. The possibility of an overarching connection between the cold MJO and these wave classes is explored in further detail in a companion paper.
It is reasonable to consider an alternate interpretation of Fig. 11a—that the modes we have simulated are not in fact MJO-like. This follows if one assumes H2 is true—that the propagation of the MJO is already known to be due to horizontal moisture advection based on the range of motivating evidence reviewed in the introduction (e.g., Maloney 2009; Maloney et al. 2010; Kiranmayi and Maloney 2011; Sobel and Maloney 2012; AK12; Pritchard and Bretherton 2014). Under this assumption, lack of strong budget projections implicating horizontal advection, or the predominance of vertical advective terms in our coldest climates, could be viewed as inconsistent with the conclusion that the modes in our simulations are MJO-like. To some extent this is a philosophical question about the extent to which MJO propagation dynamics are well understood, and what it means diagnostically to call a slow eastward-moving intraseasonal mode MJO-like.
Our own view is that the question of why the MJO propagates to the east at its observed rate is not yet fully settled, and that the evidence in Figs. 7–9 is sufficient to claim the modes we have simulated are MJO-like. From this view, if one accepts that the causes of correlations in the MJO’s MSE budget are difficult to fully prove, and if one accepts the possibility that another process may play a controlling role in regulating its eastward propagation, then the existence of an MJO-like signal without a corresponding advective propagation term may be an important, if previously unencountered, counterexample to the advective propagation hypothesis.
2) Maintenance dynamics
Longwave heating anomalies provide the largest and most robustly detectable contribution to MJO MSE maintenance across the entire range of climates (Fig. 11b). For warm climates, this is consistent with ideas from modern moisture mode theory that radiative effects of high clouds and water vapor play an important role in locally energizing the MJO. But for cold climates, this provides a negative answer to our initial question motivated by H1 about whether the MJO might shutdown due to temperature criticalities of longwave feedbacks of RCE self-aggregation, assuming they are an analogy to the MJO. That is, if SPCAM3.0’s MJO is a form of self-aggregation, it is not a temperature-critical one. This is consistent with ideas in other studies that aggregation involving significant longwave effects can occur on small spatial scales even at very cold temperatures under RCE in limited-domain cloud resolving simulations (Abbot 2014) or via shortwave feedbacks in tropical channel simulations (Wing and Cronin 2016). This is also relevant to questions raised by Emanuel et al. (2014) about the generality of a minimum temperature threshold for aggregation. Our results support the idea that longwave radiative heating anomalies can associate with signals of convective self-aggregation in very cold climates, and demonstrate that this can occur on global scales in the form of cold-MJO-like signals.
Our results suggest that physics other than the efficiency of vertical MSE advection may play a role in orchestrating the amplification of the MJO with SST warming, an idea that is further explored in a companion paper. The climate sensitivity of the MJO in the vicinity of current temperatures has been linked to systematic changes in the efficiency of vertical MSE advection in its maintenance dynamics (Arnold et al. 2013). Consistent with an extension of this idea across a broader SST range, the efficiency of this damping term is strongly reduced as SST is warmed from 1° to 30°C (Fig. 11b). However, one difference in our findings is that we do not detect a transition from a dissipative to an amplifying role (sign flip) for the vertical advection term in the 25°–35°C SST range, indicating that this aspect of the MJO response may be basic-state dependent. We also find an interesting counterexample to this mechanism in that the efficiency of the vertical advection term is not detectably sensitive to increasing SST from 30° to 35°C in our simulations, despite MJO amplification in this climate regime (Fig. 10).
We have performed idealized experiments using a novel global climate model (SPCAM3.0) that produces a realistic MJO signal while making minimal assumptions about how moist convection interacts with large-scale dynamics. We have used this model to test two hypotheses (H1 and H2) motivated by a modern view of the “moisture mode” theory of the MJO and an analogy (A1) that the MJO is an analog to temperature-critical RCE self-aggregation. Inconsistent with A1, we have found that the MJO does not shut down even for very cold climates down to 1°C. Evidence of a cold-MJO signal exists from a variety of diagnostic vantage points. Apparently the intrinsic MJO signal in SPCAM3.0 is a surprisingly resilient mode that can survive extreme cooling of the climate. Inconsistent with H2, the MJO in SPCAM3.0 can also survive a reversal of near-equatorial meridional MSE gradients that occurs in our coldest simulations due to exotic basic-state responses to homogeneous SSTs.
The existence of a cold MJO in our simulations might raise some troubling concerns about the validity of a moisture mode interpretation, including ideas that the MJO is propagated by horizontal moisture advection or critically dependent on SST > 25°C. The fact that the cold MJO can occur even in a weakly sheared easterly basic state is also hard to reconcile with ideas related to vertical shearing of its longwave heating anomaly as an MJO propagation mechanism. Resilience of familiar phase relationships among terms in the MJO-column MSE budget in our simulations (leading longwave maintenance and vertical MSE advective dissipation) could indicate some broader cause requires such correlations to exist.
Convective self-aggregation arguments could nonetheless be invoked to explain the MJO signals in our simulations via dynamics that are unfamiliar insofar as they would rely on positive longwave feedback even at cold SST and could not rely on horizontal MSE advection as the primary propagation mechanism. The fact that longwave heating anomalies are in phase with MSE anomalies is consistent with H1 and could explain the maintenance of the MJO through aggregation at low effective gross moist stability. A horizontal phase offset in vertical motion and balanced equatorial beta-plane dynamics could then require moist anomalies to move east or west independent of horizontal advection (C. Bretherton 2015, personal communication). Detailed analysis of the initial development of moist anomalies might help elucidate if such dynamics help cause a cold MJO in a horizontally uniform SST setting.
An alternate hypothesis is that dynamics more fundamentally tied to high-frequency equatorial waves and scale interactions provide an overarching control on the MJO across our simulations, consistent with the gravity wave–like signatures in its composite MSE balance in Fig. 11a. This is also consistent with the fact that the MJO can be simulated without interactive radiation at constant SST in global models with explicit convection (Grabowski 2001) and without meridional moisture gradients in conventionally parameterized aquaplanet experiments (Hsu et al. 2014). In a companion paper, we test predictions from the recent theory of Yang and Ingersoll (2013, 2014) and find they exhibit skill in quantitatively predicting both the zonal enlargement and speedup of the MJO seen in our simulations, consistent with their idea that the MJO may be a large-scale envelope of high-frequency small-scale tropical waves.
Meanwhile, this simulation suite provides a useful test bed for modern theories of the MJO since it quantifies MJO amplification, enlargement, and speedup across the 1°–35°C climate regime, in a relatively simple basic state. It would be worth testing predictions from more MJO theories against these simulation results, and testing their resilience across multiple MJO-permitting global climate models.
Funding for this work was provided by the U.S. Department of Energy under DE-SC0012152 and DE-SC0012548 and by the National Science Foundation (NSF) under AGS-1419518. Da Yang is supported by a fellowship from the Miller Institute for Basic Research in Science at UC Berkeley. Computational resources were provided by the Extreme Science and Engineering Discovery Environment, which is supported by NSF Grant OCI-1053575, under Allocation TG-ATM120034.
Supplemental information related to this paper is available at the Journals Online website: http://dx.doi.org/10.1175/JCLI-D-15-0790.s1.
Publisher’s Note: This article was revised on 28 June 2016 to fix a typographical error in the title.