Role of convection–circulation coupling in the propagation mechanism of the Madden–Julian Oscillation over the Maritime Continent in a climate model

This study investigates the role of convection–circulation coupling on the simulated eastward propagation of the Madden–Julian Oscillation (MJO) over the Maritime Continent (MC). Experiments are conducted with the European Centre Hamburg Model Version 5 (ECHAM5) coupled with the one-column ocean model—Snow-Ice-Thermocline (SIT) and two different cumulus schemes, Nordeng-Tiedtke (E5SIT-Nord) and Tiedtke (E5SIT-Tied). During the early phase of MJO composites, the E5SIT-Nord simulation reveals stronger intraseasonal anomalies in the apparent heat source (Q1) over the convective center, however, the E5SIT-Tied produces a stronger background Q1, suggesting that deep convection prevails over the MC but does not couple with the MJO circulation. Similarly, in the E5SIT-Tied simulation, in-column moisture is kept mostly by local deep convection over the MC, which is in contrast to the well-correlated relationship between moisture anomaly and MJO circulation in E5SIT-Nord. A case study based on an observational MJO reveals similar biases concerning of convection–circulation coupling emerges within a few days of simulations. The E5SIT-Tied simulation produces weaker heating at the convective center of the MJO than the E5SIT-Nord a few days after model initiation, resulting weaker subsidence to the east and less favorable for propagation. The present findings highlight the instantaneous responses of cumulus parameterization schemes to MJO-related environmental changes can further affect intraseasonal variability through altering convection–circulation coupling over the MC. Physical schemes of moist convection are essential to realistically represent this coupling and thereby improve the simulation of the eastward propagation of the MJO.


Introduction
The Madden-Julian Oscillation (MJO) is the dominant pattern of atmospheric intraseasonal variability in the tropics (Madden and Julian 1972;Zhang 2005;Jiang et al. 2020). Numerous observational and modeling studies have explored the fundamental physics regarding the eastward propagation and development of the MJO, including the interaction with the frictional boundary layer (Wang and Rui 1990;Salby and Hendon 1994;Maloney and Hartmann 1998;Hsu et al. 2004;Kang et al. 2013), convective intraseasonal air-sea interactions (Flatau et al. 1997;Waliser et al. 1999), ocean surface flux (Maloney and Sobel 2004;Maloney 2009;Kiranmayi and Maloney 2011;Andersen and Kuang 2012), moisture transport by cumulus (Benedict and Randall 2007;Jiang et al. 2016;Jiang 2017), and air-sea interaction (DeMott et al. 2015(DeMott et al. , 2019. Despite the considerable progress in state-of-the-art general circulation models (GCMs) in recent decades, realistic simulations of the MJO remain difficult.
One challenge to improving models is the MJO is organized by strong coupling between MJO circulations, convection, and ocean thermodynamics. Studies have addressed numerous crucial physical processes for reproducing the MJO in GCMs, including the mean state (Kim et al. 2011;Kim 2017;Jiang et al. 2018;Klingaman and Demott 2020), convection simulated by cumulus parameterization schemes (CPS; Liu et al. 2005;Deng and Wu 2010;Zhou et al. 2012), convectively coupled tropical waves (Kiladis et al. 2005;Janiga et al. 2018), ocean-atmosphere interaction (DeMott et al. 2015(DeMott et al. , 2019Tseng et al. 2015), the diurnal cycle of sea surface temperature (SST; Bernie et al. 2005;Klingaman et al. 2011), and cloud radiative feedback (Jiang et al. 2011;DelGenio and Chen 2015;Ciesielski et al. 2017). These findings jointly reflect the complexity of the physical processes involved in understanding and simulating the MJO.
To untangle this multifaceted problem, all coupling mechanisms must be carefully examined. One key aspect of the MJO behavior manifested of such coupling is MJO's eastward propagation from the Indian Ocean to the Maritime Continent (MC) and onward to the western Pacific (Zhang 2005). The eastward propagation of the MJO over the MC has been challenging to simulate in GCMs, and remains so, even with improvements found in the Coupled Model Intercomparison Project Phase 6 models (CMIP6; Ahn et al. 2020a, b;Fiedler et al. 2020). After the MJO appears over the Indian Ocean, deep convection in the oscillation induces the dynamic processes associated with the equatorial wavelike perturbations to the east over the MC, further enhancing both low-level moisture convergence and moistening (Zhang 2005;Jiang et al. 2020). Observational studies have suggested that the shallow convection with bottom-heavy heating profiles that occurs before the propagation of the MJO convective center helps precondition the development of deep convection (Kemball-Cook and Wang 2001;Kikuchi and Takayabu 2004). Such shallow convection that is associated with low-level convergence and ascending motion moistens the lower troposphere and enhances coupling with intraseasonal perturbations (Benedict and Randall 2011;Lappen and Schumacher 2014). In a model experiment with prescribed heating profiles, Lappen and Schumacher (2014) also demonstrated the essential role of the low-level heating in the development and maintenance of the MJO.
Multiple studies have also reported that including air-sea coupling aids MJO simulation in climate models (Tseng et al. 2015;Jiang et al. 2015Jiang et al. , 2019DeMott et al. 2018DeMott et al. , 2019. Crueger et al. (2013) explored the role of CPS, ocean coupling, and model resolution by evaluating the simulation of the MJO with the European Centre Hamburg Model (ECHAM) Version 6 under various configurations. They have found that CPS plays an essential role in simulating MJO signatures, along with ocean coupling. The researchers attributed the discrepancy to differences in the sensitivity of the two CPS to tropospheric moisture (Crueger et al. 2013). However, they did not discuss the role of convection-circulation coupling in their analysis. By jointly using the ECHAM5 and their self-developed one-column ocean model, the Snow-Ice-Thermocline (SIT) model, Tseng et al. (2015) successfully simulated the eastward propagation of the MJO, observing that the low-level moisture convergence prior to the active MJO phase over the MC played a pivotal role in preconditioning the propagation of deep convection and circulation. They found that the upper oceanic warm layer to the east of organized deep convection gets warmer with more incident solar radiation and reduced evaporative cooling over the subsidence region. The resulting warm SST ahead of the MJO was suggested to destabilize the boundary layer and enhance frictional convergence, underscoring the importance of shallow moistening for the further development of the MJO (Tseng et al. 2015). Comparing between models, Jiang et al. (2015) reported that the ECHAM5-SIT model was one of the best GCMs in MJO simulations. In another study, that model was used to determine the role of air-sea coupling in MJO simulations (DeMott et al. 2019). Although both studies have evaluated the performance of the model in simulating the MJO, the contribution of convection has not been addressed. Thus, the key mechanism underlying convection-circulation coupling, which may improve MJO simulations, remains unclear.
To examine the role of convection in the mechanism of MJO propagation over the MC, we use the model setup described by Tseng et al. (2015) to investigate the convection-circulation coupling in the ECHAM5 model experiments with two CPSs, namely the Nordeng-Tiedtke scheme (E5SIT-Nord) and the Tiedtke scheme (E5SIT-Tied). By using the model setup with a high resolution one column ocean model, we can investigate the role of convection-circulation coupling for simulating the eastward propagation of MJO in GCMs without interferences of possible biases from air-sea coupling. Our analysis will mainly focus on the role of oceanic convection over the MC regions for MJO propagation. Besides the advantage provided by ECHAM5-SIT, this analysis focus is based on the fact that oceanic convection over MC plays an important role in sustaining MJO-scale circulation, while the land convection is found to compete with MJO-related convection in both observations and models (Hagos et al. 2016;Ahn et al. 2020a). How the land convective activity may affect the MJO propagation will be a great topic for future study.
In this study, we adopted the budget analysis of apparent heat sources (Q 1 ) and apparent moisture sinks (Q 2 ), as well as the moist static energy (MSE) budget, to the 25-year climate simulations of two model setup, with the central aim of identifying the contribution of the moisture buildup over the MC in the propagation of the MJO. To further clarify the effects of the representation of CPS on MJO propagation, we conduct a case study by using the initial conditions from the ERA-Interim reanalysis and examining the MJO evolution simulated by the two CPSs. By observing how the model biases evolve in case study provides more information about how biases in climate runs may link to the fast physics from CPS related to convection-circulation coupling. The remainder of this paper is organized as follows. The data, model, and methodology are described in Sect. 2. Sections 3 and 4 present the analytical results and the results from the case study, respectively. The discussion and conclusion are in Sects. 5 and 6, respectively.

Data, models, and methodology
In the present study, the ECHAM5 (Roeckner et al. 2003(Roeckner et al. , 2006 is used in combination with the SIT (Tu and Tsuang 2005;Tsuang et al. 2009). As its name implies, the ECHAM5 is the fifth version of the ECHAM, a GCM developed at the Max Planck Institute for Meteorology. The horizontal resolution used is T63 (~ 1.8°), with 31 vertical layers and a model top at 10 hPa (approximately 30 km). Although not the most updated version of ECHAM, previous studies have found that, when coupling with SIT, ECHAM5 can realistically simulate many critical features and processes of MJO (Tseng et al. 2015;DeMott et al. 2019). ECHAM5-SIT also had been used to participate in the model intercomparison study hosted by the MJO task Force, under the auspices of Year of Tropical Convection and later the Working Group on Numerical Experimentation, and was identified as one of the better skill models for MJO simulations in many studies (Jiang et al. 2015;Wang et al. 2017;Zhu and Li 2020;Heath et al. 2021).
SIT simulates the SSTs and variability in upper ocean temperature, including that attributable to the cool-skin and diurnal warm-layer effects in the upper ocean. The turbulent kinetic energy of a water column is also simulated (Gaspar et al. 1990). In our experiments, SIT has 42 vertical layers, 12 of which are in the upper 10 m. In one study (Tu and Tsuang 2005), the SIT, with a 1 m resolution in the upper 10 m and a layer at 0.05 mm for reproducing the cool-skin effects at the ocean surface, realistically simulated the near-surface warm layer (Tu and Tsuang 2005). It is not conventional to combine such a high-vertical-resolution ocean model of turbulent kinetic energy to an atmospheric GCM (AGCM). To account for horizontal processes, the ocean model is weakly nudged (with a 30-day time scale) at depths exceeding 10 m and strong nudged (with a 1-day time scale) below 100 m to the observed climatological ocean temperature. The SIT and ECHAM exchange SST and surface fluxes at every time step (12 min) in the tropics (30° S-30° N). Elsewhere, climatological SST drives the AGCM.
Two CPS used in the present study are the Nordeng-Tiedtke scheme, the default scheme used in the ECHAM5 (Nordeng 1994), and the Tiedtke scheme (Tiedtke 1989). The Nordeng-Tiedtke scheme, an improved version of the Tiedtke mass flux convection scheme, considers organized entrainment and detrainment in buoyancy-related penetrative convection (Nordeng 1994). The differences between the two schemes lie mainly in the representation of deep convection regimes, particularly closures and entrainment rates (Möbis and Stevens 2012). For closures, the Tiedtke scheme uses the moisture convergence within the subcloud layer to make approximation for convective intensity (Tiedtke 1989). By contrast, the Nordeng-Tiedtke scheme uses the quasiequilibrium closure for the relaxation of the vertically integrated buoyancy (Nordeng 1994). For the entrainment rate, the Nordeng-Tiedtke scheme uses additional lateral mixing; specifically, the buoyancy-driven entrainment for deep convection is included in the organized entrainment module (Nordeng 1994).
In the present study, the simulation results are compared to data from the Global Precipitation Climatology Project (GPCP: Adler et al. 2003), data on outgoing longwave radiation and daily Optimum Interpolation SST from the National Oceanic and Atmosphere Administration (Banzon et al. 2014), and atmospheric states from the ERA-Interim reanalysis (Dee et al. 2011). We used 22-year ERA-Interim reanalysis from 1989 to 2010 and 14-year GPCP dataset from 1997 to 2010. The simulation diagnostics package developed by the Climate Variability and Predictability MJO Working Group (CLIVAR Madden-Julian Oscillation Working Group 2009) and a 20-to-100-day filter are used to isolate and determine intraseasonal variability. The MJO phases are classified on the basis of the Real-time Multivariate MJO (RMM) index, which uses the principal components (PCs) of the two leading pairs of empirical orthogonal function modes of intraseasonal outgoing longwave radiation and zonal winds at 850 and 200 hPa between 15° S and 15° N to represent MJO features (Wheeler and Hendon 2004). Because of the strong eastward propagation tendency of the MJO, the analysis is centered on the boreal cool season (November-April). The total simulation length of each experiment is 25 years, with data in the last 24 years subjected to analysis.
To understand the simulated convection, we calculate the Q 1 and Q 2 associated with the MJO (Yanai et al. 1973), which are presented in the following equations: where v represents horizontal winds, represents vertical pressure velocity, c and e represent condensation and evaporation, respectively; q is the water vapor mixing ratio; s is the dry static energy (s = c p T + gz, where c p is specific heat capacity, T is temperature, g is gravitational constant, and z is height); and L v is the latent heat of vaporization. As for − (⋅) , it denotes the grid scale quantity relative to the subgrid components. For the convection-dominated regions, convective activities could be inferred by comparing the Q 1 and Q 2 vertical profiles to reflect the vertical transport of the MSE by convection. For example, deep convection regimes often have top-heavy Q 1 vertical profiles with peaks at approximately 300 hPa and Q 2 profiles with peaks at 700 hPa, indicative of upward convective transport from the lower troposphere (Yanai et al. 1973). By contrast, Q 1 vertical profiles in shallow convection regimes tend to be bottomheavy, with peaks at lower atmospheric pressures (Nitta and Esbensen 1974). The vertically integrated MSE budget is defined as follows: where h is the MSE (h = c p T + gz + L v q); u and v are the zonal and meridional velocities, respectively; ω is the vertical pressure velocity; LW and SW are the longwave and shortwave radiation fluxes, respectively; and LHF and SHF are the latent and sensible surface heat fluxes, respectively. The mass-weighted vertical integration from the surface to 200 hPa is denoted by 〈⋅〉, and intraseasonal anomalies are represented by 〈⋅〉'. All fields are isolated using a 20-to-100-day band-pass Lanczos filter (Duchon 1979).

ECHAM-SIT simulations of the MJO
The simulations of the MJO in two experiments involving the E5SIT-Nord and E5SIT-Tied are examined in this section. Figure 1a-c presents the wavenumber-frequency spectra of simulated 850 hPa zonal winds (shading) and precipitation (contours). The E5SIT-Nord realistically simulates the 30-to-80-day eastward-propagating signals of planetary wavenumber 1 in circulation and wavenumber 1-3 in precipitation ( Fig. 1b), which is consistent with the finding of Tseng et al. (2015). By contrast, the spectra of the E5SIT-Tied simulation exhibit weaker signals in both circulation and precipitation over wavenumber 1 (Fig. 1c). The remaining parts of the Fig. 1d-f illustrate the eastward propagation of intraseasonal fluctuations with the lag-longitude diagram of precipitation and 10 m zonal winds regressed to the rainfall over the tropical Indian Ocean (10° S-5° N, 75° E-100° E). The E5SIT-Nord simulation reproduces the distinctly observed eastward propagation, albeit slightly more slowly, of precipitation and surface winds (Fig. 1e), whereas the E5SIT-Tied simulation presents weak intraseasonal fluctuations in precipitation which is decoupled from zonal winds at the MC, especially for the region of 120°-150° E (Fig. 1f). Therefore, our analysis focuses on the convection-circulation coupling associated with MJO over the MC region. The mean states of zonal winds and precipitation rate are shown in Fig. S1, and maps of the ratio of the intraseasonal rainfall variability to the total rainfall variability are presented in Fig. S2. In line with the results of Tseng et al. (2015), the two model simulations have relatively similar mean states (Fig. S1). Both models simulate a higher ratio of intraseasonal rainfall variability to total variability than was observed, with the E5SIT-Nord model simulating the highest variability among the three datasets (Fig. S2b). By contrast, the variability simulated by the E5SIT-Tied model is weaker than what was observed (Fig. S2c). While the contribution of the mean distribution of lowlevel moisture to MJO propagation has been emphasized in numerous studies (e.g., Jiang 2017), it appears not to be a major factor in our simulations, considering, as shown in Fig. S3a-c, the dryer mean state of the E5SIT-Nord simulation (relative to the E5SIT-Tied simulation), especially below 700 hPa (Fig. S3d). An examination of model pattern correlations with ERA-Interim (from 1000 to 200 hPa) also reveals that both models effectively simulate distributions of horizonal moisture; the coefficient of correlation in the both runs exceeds 0.95 (Fig. S3e). Other studies have also shown that the vertical moisture distribution is not sufficient conditions for simulating MJO, and suggest the net vertical moistening associated with convective processes which coupled with large-scale environments are a key element for MJO simulation (Jiang et al. 2015;Peters et al. 2017). While the mean moisture pattern appears not to be the key differences between the two simulations, the convective processes that result in more accurate simulation of MJO propagation by the E5SIT-Nord are worthy of further exploration.
The autocorrelation analysis of daily rainfall of two runs and GPCP in Fig. S4 shows the coherence behaviors of convection in the tropics. In observations of GPCP, tropical rainfall signal is composed of strong alternating spells of wet and dry period (i.e. intermittency), which is shown by the reduction of autocorrelation function with lag days. Examining rainfall coherence can help understand the tendency of raining processes in climate models (Klingaman et al. 2017;Peters et al. 2017). Compared with GPCP and E5SIT-Tied, E5SIT-Tied shows weaker intermittency with tropical daily rainfall, especially for timescale longer than 5 days. This result suggests convection in E5SIT-Tied tends to be less organized. Such characteristics may hamper development of MJO organization in E5SIT-Tied run.
To examine the convection response associated with MJO circulations, the following analysis in Figs. 2, 3, 4 and Fig.  S5 will focus on the region of 10° S-EQ and 120°-150° E, spanning the area from the Banda Sea to the Maluku Islands and Java, where the strongest convective activities on intraseasonal scale in the MC region (Fig. S2). It is also the region where the MJO convective center often makes a southward detour when passing through the MC (Tseng et al. 2017). The evolution of convection activity of this region is closely coupled with MJO circulations during the MJO passage.
The evolution of essential atmospheric variables for the initiation and development of convection over the region of 10° S-EQ and 120°-150° E are plotted with respect to the MJO life cycle in Fig. 2. All variables underwent 20-to-100-day band-pass filtering and are plotted against the simulated MJO phases, per the procedures of Wheeler and Hendon (2004). Figure 2a-c shows the moisture convergence associated with the simulated MJO phases. The low-level convergence over the MC is the strongest below 850 hPa in phases 1 and 2, and it continues shifting upward into the free troposphere starting from phase 3 (Fig. 2a). Both simulations capture the structure of low-level moisture convergence in the first two phases, but at smaller magnitudes. As shown in Fig. 2b, c, the simulated moisture convergence shifts to the free troposphere relatively quickly between phases 3 and 4 and becomes detached from the boundary layer. Moisture convergence is stronger and extends to higher levels in the E5SIT-Nord simulation than in the E5SIT-Tied simulation. The associated Q 1 and Q 2 over the MC region are also observable in Fig. 2d-f and g-i, respectively. E5SIT-Nord has a sharp contrast between convective heating and cooling during the MJO life cycle that is greater in amplitude than that in the E5SIT-Tied simulation and in the observation data. Consistent with the wave spectrum signals in Fig. 1, the magnitudes of Q 1 and Q 2 in the E5SIT-Tied simulation are substantially weaker than those in the E5SIT-Nord simulation. Notably, the observed Q 2 exhibits a distinct tilt, and both models simulate a rapid transition from the near-surface features to the deep convective structure in the middle troposphere. Figure 2j-l presents the evolution of the MSE profiles over the MC associated with MJO phases. The vertical development of the MSE, specifically through the upward transport of moisture from the near-surface to the low-level troposphere, results in conditions conducive to deep convection (Fig. 2j). The E5SIT-Nord simulates the gradual upward development of MSE with time, whereas the E5SIT-Tied (Fig. 2k, l) simulates weak preconditioning with an unstable near-surface layer, followed by the sudden onset of high MSE and deep convection in mature phases. As noted by Tseng et al. (2015), this moistening process, characterized by the tilting of MSE, was only captured when the ECHAM5 was used in concert with the SIT.
As tropical convection has a multiscale nature, convective activities on other time scales can affect MJO-related convection development. Figure 3 presents the unfiltered Q 1 and Q 2 in eight MJO phases, corresponding to the intraseasonal ones shown in Fig. 2d-i. In the reanalysis, the Q 1 and Q 2 profiles (Fig. 3a, d) reveal the evolution of the convective activities in the MJO, from the shallow convection in the earlier phases of phase 1-2 to the deep convection in the mature phases of phase 3-5. Both models simulate heating or moistening over all phases that are stronger in magnitude than their observed counterparts. In the E5SIT-Nord experiment, the simulated Q 1 exhibits heating during the mature phases that is comparable to its observed counterpart; however, in phase 4, heating is greatly overestimated by as much as 4 K/day (Fig. 3b). By contrast, in the E5SIT-Tied experiment, Q 1 heating persists from phase 3 to phase 8 (Fig. 3c). Notably, the mean Q 1 in the E5SIT-Tied simulation is larger than that in the E5SIT-Nord simulation when averaged over the entire MJO life cycle, the opposite of the Q 1 projection on the intraseasonal time scale (Fig. 2). Unlike the prominent Q 1 associated with MJO convective center in the observation data and the E5SIT-Nord simulation, the Q 1 heating in the E5SIT-Tied simulation persists from phases 1-8 of the MJO, demonstrating the maintenance of deep convection over the MC. This result and its relationship with lower simulation skill of MJO is explored in later sections. Fig. 2 Vertical profiles of intraseasonal anomalies (i.e., with 20-to-100-day band-pass filtering) in a-c moisture convergence (10 −6 g kg −1 s −1 ); d-f the apparent heat source (Q 1 ; K day −1 ); g-i the apparent moisture sink (Q 2 ; K day −1 ); and j-l moist static energy (10 3 J kg −1 s −1 ), with respect to Madden-Julian Oscillation phases averaged over 10° S-EQ, 120°-150° E. The variables are based on observational data from the ERA-Interim reanalysis, the Nordeng-Tiedtke scheme (E5SIT-Nord) simulation, and the Tiedtke scheme (E5SIT-Tied) simulation We focus on the convection activities over the MC during phases 1 and 2 of the MJO, when the deep convection center is over the Indian Ocean. Figure 4 shows the Q 1 /Q 2 profiles over 10° S-EQ and 120°-150° E during these phases in the ERA-Interim analysis, in the E5SIT-Nord simulation, and in the E5SIT-Tied simulation. In phase 1, the peak of Q 1 in the reanalysis exhibits relatively homogeneous heating at a maximum rate of 1 K/day, suggesting a relative clear atmosphere with radiative heating and shallow convective heating dominate Q 1 (red line in Fig. 4a). The Q 2 profile peaks at 700 hPa, indicating drying within the planetary boundary layer (PBL) and moistening in the lower troposphere. Such profiles suggest that during this period, moisture is transported from the PBL to the lower troposphere over the MC (blue line in Fig. 4a) through bottom-heavy convection. In the E5SIT-Nord, during phase 1, the Q 1 value peaks at approximately 700 hPa, with a bottom-heavy profile. At the same time, the Q 1 and Q 2 peaks correspond to each other, implying that convection is shallow at this point, mainly moistening the low-level troposphere (Fig. 4b). By contrast, in the E5SIT-Tied simulation, the Q 1 and Q 2 profiles peak near 500 and 700 hPa, respectively, suggesting deep, intense convection over the MC (Fig. 4c). As shown in Fig. 4d-f, deep convective heating at approximately 500 hPa is higher in phase 2 in both the ERA-Interim analysis and the two simulations. As is the case in phase 1, heating is deeper and more intense in the E5SIT-Tied simulation than in the E5SIT-Nord simulation or the ERA-Interim reanalysis. Figure S5, which presents an analysis of the MSE budget over the region of 10° S-EQ, 120°-150° E during phases 1 and 2, elucidates the differences between the two simulations and the ERA-Interim reanalysis. Among the budget terms in E5SIT-Nord, those concerning horizontal and vertical advection terms contribute positively to the tendency term, which is partially cancelled out by longwave radiation. By contrast, the E5SIT-Tied simulates less vertical advection and very weak negative horizontal advection. The cancellation by longwave radiation results in an almost negligible tendency. The principal differences between the E5SIT-Nord and E5SIT-Tied simulations appears to be in the horizontal advection term and vertical advection terms, which further result in the poorer ability of E5SIT-Tied to simulate the propagation of the MJO. In phase 2, differences in horizontal advection that become the dominant term in the MSE budget increase as convection moves closer to the MC (Fig. S5b). Although the budget terms computed from the reanalysis have large MSE residuals, as in previous studies (Kiranmayi and Maloney 2011;Jiang 2017), vertical and horizontal advection likely to be crucial contributors to MSE tendency; this is highly consistent with that calculated from the E5SIT-Nord simulation. Figure 5 shows the MJO structures in phases 1 and 2, with overturning circulation (vectors), Q 1 (shading), and moisture convergence (green contours). The horizontal maps represent the precipitation and sea level pressure, respectively (bottom panel in Fig. 5). In phases 1, the MJO convective center is located over the central Indian Ocean (around 80°-90° E), with an upward motion driven by deep convective heating. To the east, subsidence is induced above most of the MC (120°-150° E) with lowlevel moisture convergence developed below 850 hPa (Fig. 5a). During this phase, E5SIT-Nord simulates strong subsidence and low-level moisture convergence over the MC and low-level easterlies over the eastern Indian Ocean (Fig. 5b). By contrast, E5SIT-Tied simulates very weak subsidence and easterlies over the MC and the eastern Indian Ocean, respectively, indicating weak convective responses to MJO-scale circulations (Fig. 5c). In phase 2, MJO convection propagates eastward with more intense precipitation covering Indian Ocean and Sumatra island. Both subsidence and Kelvin-wave structure are enhanced in response to stronger convective activities (Fig. 5d). Similar to the phase 1, E5SIT-Tied exhibits much weaker deep convective activities, weaker overturning structure, and little eastward propagation compared with observations and E5SIT-Nord (Fig. 5d-f).
To further visualize the coupling of moisture buildup and the MJO over the MC, we examine the large-scale circulation patterns and thermodynamics during periods when the moisture buildup over the MC is the most substantial. We first define intraseasonal moisture buildup by characterizing the duration of moisture residence in the tropospheric column before removal through rain ([∫ p T p B qdp ]/[P]) over the entire MC (10° S-10° N, 120°-150° E), where [·] is the intraseasonal filtering, q is moisture, P is the rain rate, p T is the pressure at the cloud top, and p B is the pressure at the cloud bottom. The normalized convective time scale over the southern MC are calculated with respect to the lag-regressed fields to the series 1 of RMM (RMM1) index with ERA-Interim and two model runs in Figure S6. The RMM1 is defined as the PC of the leading EOF mode representing the coupled structure of convection and circulations when the convection center of MJO occurs at the longitude of MC (Wheeler and Hendon 2004). In ERA-Interim, convective time scale over the MC increases gradually from phase 1 (i.e., − 20 days), reaching the maximum in phase 2 (i.e., − 12 days; green line in Fig.  S6), indicating moisture buildup peaks around 10 days before the deep convection activity on day 0 and under dominant shallow convection (from Q 1 /Q 2 in Fig. 2). The longitude-height cross-sections (averaged over 10°S-EQ) of the MJO-scaled wind circulation (vector, u: ms −1 , omega: 10 −2 Pa s −1 ), Q 1 (shaded area; K day −1 ), and the horizontal moisture convergence (green contour; the contour lines are shown only in 0 and 3; 10 −6 g kg −1 s −1 ) in the a, d ERA-Interim rea-nalysis; b, e Nordeng-Tiedtke scheme (E5SIT-Nord) simulation, and c, f Tiedtke scheme (E5SIT-Tied) simulation. The contour interval of the moisture convergence is 8 × 10 −6 g kg −1 s −1 . The solid line is positive. Precipitation (shaded area; mm day −1 ) and sea level pressure (contour; hPa). The contour interval is 30 hPa. The dashed line is negative Fig. 6 Regressed structure of circulation patterns, the apparent moisture sink (Q 2 ), and moist static energy (MSE) on intraseasonal convective time scale (shown in Fig.S6) in the a ERA-Interim reanalysis, b Nordeng-Tiedtke scheme (E5SIT-Nord) simulation, and c Tiedtke scheme (E5SIT-Tied) simulation. The MSE (J kg −1 s −1 ) is shown in the shaded area, Q 2 (K day −1 ) is shown in black contours, and the circulation winds (ms −1 ) are shown as vectors where the value less than 0.1 is not shown. The left panel shows the horizontal maps of 850 hPa fields, and the right panel shows the height-longitude crosssections averaged over 10° S-EQ. The black contour intervals in the left and right panels are 0.2 and 0.1, respectively and E5SIT-Tied simulation. A distinct feature of MJO with convective center at the Indian Ocean and dynamical responses are shown in both the reanalysis and the E5SIT-Nord simulation. Notably, a realistic structure similar to Kelvin waves with two off-equatorial MSE maxima over the eastern Indian Ocean under strong easterlies at the equatorial MC is present only in the E5SIT-Nord simulation. On the other hand, the missing of the deep convection signal over the east Indian Ocean and west MC are shown in the E5SIT-Tied run (Fig. 6c). Also in reanalysis and E5SIT-Nord, the zonal-height cross-section exhibits a tilted structure of well-collocated MSE and Q 2 , which depicts a MJO structure in the developing phase (at approximately phase 2) (Fig. 6a, b). The tilted structure involves deep convection near 90° E, relatively shallow convection at 120° E, and low-level moistening by convective processes between 120° E and 150° E. The Q 2 tilt in the E5SIT-Nord is more rapid than that in the ERA-Interim reanalysis (as indicated by black contours in the right panels of Fig. 6a, b). Such westward-tilting vertical structures have been identified in numerous observational studies (Kiladis et al. 2005;Kim et al. 2009;Tseng et al. 2015;Jiang et al. 2015). By contrast, the E5SIT-Tied simulation reveals local deep convection and a MSE maximum over the MC (120°-150° E), implying the MSE anomaly over MC couples only with local convection activity, not with the large-scale circulation patterns that are further linked to deep convection in the Indian Ocean. The meridional cross-section over 120°-150° E also indicates the anomalies in the deep Q 2 and the MSE over the MC in the E5SIT-Tied simulation (Fig. S7). By contrast, both the E5SIT-Nord simulation and the ERA-Interim reanalysis exhibit a shallower Q 2 and MSE structure, with the maximum at approximately 800 hPa. As in the ERA-Interim reanalysis and the E5SIT-Nord simulation, this is a precondition for the development of deep convection over the MC and the eastward propagation of the MJO when convection is still over the Indian Ocean.
The characteristics of convective moistening (i.e., Q 2 ) to rainfall categories and SST are presented as follows to investigate the impacts of convective process on environments.
We especially focused on oceanic convection in our analysis to take the advantage of ECHAM5-SIT model setup in better representation of the atmosphere-ocean coupling. Also, many studies have found the land convection over the MC competes moisture with MJO-scale convection, having a blocking effect rather than supportive effect for MJO propagation (Hagos et al. 2016;Ahn et al. 2020a). Examining interactions between land convective activity and the MJO propagation is beyond the scope of this study and will be a great topic for future study. Figure 7 shows the Q 2 profiles stratified by the daily rain rate averaged over the tropical oceans (10° S-5° N, 80°-160° E) in the two model runs and their differences. For low rain rates (1-10 mm/day), the E5SIT-Nord simulates a more bottom-heavy Q 2 profile than the E5SIT-Tied does, transporting moistening from PBL into the low troposphere. For high rain rates, E5SIT-Nord simulates stronger drying (more latent heating) in the middle troposphere associated with strong convective activity, compared with E5SIT-Tied. The two characteristics (i.e., stronger low-level moistening and stronger deep drying) of the E5SIT-Nord simulation help create more favorable environmental conditions for the propagation of the MJO. Figure 8 shows the percentage of rainfall occurrences stratified with respect to SST and rainfall intensity. The highest rainfall is noted between 28 and 30 °C (Fig. 8a, b) in both the observation and the E5SIT-Nord simulation. By contrast, the E5SIT-Tied simulates more frequent rainfall when the SST is higher than 30 °C, suggesting the E5SIT-Tied tends to be overresponsive to high SST, producing an excessive amount of local convective rainfall (Fig. 8c).

Case study on the MJO involving two CPSs
The difference in the representation of convection-circulation coupling in the two simulations may not be directly ascribable to the representation of the convective processes in the CPS. To examine the relationship between the instantaneous convective response to the MJO circulation and convection-circulation biases found in climate runs, we conduct a simulation of an MJO case, using initial conditions from Fig. 7 Vertical profiles of Q 2 (K day −1 ) composited from categories of rainfall intensity (mm day −1 ) over the tropical oceanic regions covering10° S-5° N, 80°-160° E simulated by the a Nordeng-Tiedtke scheme (E5SIT-Nord) and b Tiedtke scheme (E5SIT-Tied); c shows the differences between the two simulations. Only the ocean grid is calculated the ERA-Interim reanalysis. Having the setup of the two schemes, the simulation is initiated at the start point of October 31, 2011, when a well-developed MJO is observed over the Indian Ocean. The two experiments are then integrated for 3 days and analyzed for their differences in days 1 and 2. This setup is designed to compare the response of moist convection processes and ensure that the large-scale circulation patterns remain close to the initial conditions. The diagnostics are useful for understanding the biases caused by the model representation of convective processes from diurnal to interannual time scales (Ma et al. 2021). Figure 9 shows the differences in circulation patterns and sea level pressure between the two simulations, in the form of a height-longitude cross-section and horizontal maps, on the first and second days after model initiation. On day 1, the convection center of the MJO is over the tropical Indian Ocean (near 90° E), with weak subsidence in the east over the MC and the western Pacific. After 1-day integration, the E5SIT-Nord simulation produces (1) stronger upward motion associated with convective center at 90° E and (2) a stronger basin-wide subsidence coupled with positive anomaly of sea level pressure over the MC and the western Pacific from 120° E to 180° E (Fig. 9a). Only spotty upward motion can be found over islands of New Guinea and Celebes as increased solar radiation under subsidence creates favored environments for local land convection (Fig. 9a). In accompany with circulations, stronger convective drying represented by Q2 in E5SIT-Nord at the convective center compared with E5SIT-Tied (shown as black contours in Fig. 9a). To the east of deep convective center of MJO, the convection over the MC in the E5SIT-Nord simulation in response to the subsidence is much shallower than that in the E5SIT-Tied simulation, especially over the Borneo, near 120° E (black contours in Fig. 9a). Notably, in the E5SIT-Nord simulation, more moisture is available in the lower troposphere and less moisture is available in the PBL below 850 hPa over the MC and the western Pacific (represented by MSE; shading in Fig. 9a); this is indicative of stronger shallow convection that helps create a favorable environment for deep convection. In the meantime, the E5SIT-Nord presents a horizontal structure similar to Kelvin waves which appears in the zonal dipole of sea level pressure, forming easterlies over the equatorial MC and the eastern equatorial Indian Ocean (bottom panel in Fig. 9a). This implies that the deep convection heating in the E5SIT-Nord simulation produces stronger dynamic responses than E5SIT-Tied simulation.
On day 2, when the convective center moves toward the western MC, different responses attributable to the two schemes become more substantial. Relative to the E5SIT-Tied simulation, the stronger upward motion at the convective center of the E5SIT-Nord simulation leads to stronger subsidence and enhanced tropical easterlies extending to the western Pacific (arrows in the bottom panel of Fig. 9b). Notably, a more prominent structure of Kelvin waves appears in response to both winds and sea level pressure (near 130° E, 10° N; bottom panel of Fig. 9b). In addition to the enhanced easterlies, the meridional winds at the north side of the equator also become stronger, which in turn enhances low-level convergence at the equator. Accompanying the vertical transport of shallow convection, more moisture accumulates over the low levels between 800 and 500 hPa in the E5SIT-Nord simulation, providing conditions more conductive for MJO propagation, as mentioned in other studies (Kiranmayi and Maloney 2011;Adames 2017;Maloney et al. 2019). Consistent with the findings on climatological MJO composites in Fig. 7, these features also suggest that E5SIT-Nord can simulate stronger convective heating in the There are 25 bins for both rainfall and SST. Only the ocean grid is calculated deep convection regime and greater moistening in the shallow convection regime than the E5SIT-Tied can.

Discussion
The comparison of the two model simulations against observational data revealed the importance of proper coupling between convection and large-scale circulation patterns in climate models to the realistic simulation of the MJO. In the E5SIT-Tied simulation, the persistent but fleeting period of deep convection dries out the moisture over the MC even during MJO development, resulting in conditions that are less favorable for MJO propagation. By contrast, E5SIT-Nord better simulates the moisture buildup over the MC during MJO development, facilitating the propagation of the MJO over deep convective regions. Figure 10 presents the differences between the model simulations with regard to the convection-circulation coupling during MJO development.
Based on our results, we suggest two aspects of convection-circulation coupling that are vital to improving simulations of MJO propagation: (1) the ability to generate deep convective heating and induce subsidence to the east of convective center and (2) the ability to induce moisture buildup in the subsidence region and create an environment that is conducive to the later development of deep convection. During the early phases of the MJO, when the convective center is located in the Indian Ocean, deep heating induces Kelvin-wave-like perturbations (e.g., easterlies and lowpressure anomalies) that extend eastward with a subsidence branch to the MC and the western Pacific (Milliff and Madden 1996;Hendon and Salby 1996;Cravatte et al. 2003;Kiladis et al. 2005). If the deep convective heating is only weakly simulated, the induced subsidence becomes too weak to suppress the deep convection over the MC region where the SST is high. The preconditioning for the propagation of the MJO to the MC also requires the buildup of low-level moisture, which is often caused by the moistening effect of shallow convection (from the transport of moisture from the PBL). Such moisture buildup is often controlled by the representation of moist convection in climate models, especially in the determination of the type of convection and the associated heating/moistening. When the convection in the model is overly sensitive to high SST, convection can be frequently triggered, depleting moisture and precluding the realistic simulation of the moistening process over the MC when MJO is approaching. Our analysis also suggests that climate models must be capable of producing the contrast between light and intense rainfall events-that is, to simulate deep heating with intense rainfall and shallow moistening with light rainfall. E5SIT-Tied is more likely to present weaker and more frequent convection than its observed counterpart and also lighter than that simulated by E5SIT-Nord. Such a tendency results in a more stable atmosphere, which makes it more challenging to model the convection-circulation coupling with MJO in E5SIT-Tied.
Numerous studies have indicated that the design of convection schemes is essential to accurately simulating the MJO, including the sensitivity of CPS to the tropospheric moisture by entrainment or detrainment (Hannah and Maloney 2014), momentum transport (Wu et al. 2007;Miyakawa et al. 2012), microphysical processes, intensity closure (Zhang and Mu 2005;Peters et al. 2017), triggering (Peters et al. 2017), the representation of stratiform convection (Fu and Wang 2009), and shallow convection (Zhang and Song 2009). The discrepancies between simulations observed in the present study could be attributable to various differences in the design of the Nordeng-Tiedtke and Tiedtke schemes, such as the sensitivity to tropospheric humidity caused by entrainment design (Möbis and Stevens 2012). We take a diagnostic approach, interpreting these differences in a context of convection-circulation coupling. To ensure the successful simulation of MJO propagation, the model representations of convection must be consistent. The similarity in the biases concerning convection-circulation coupling as observed in both the case study and the climate simulations suggests that biases in fast convective response can further induce intraseasonal biases. To improve the simulation of intraseasonal variability by better representing convection-circulation coupling, it is necessary to carefully design responses of CPS to environmental characteristics, such as wind, humidity, and SST. The process-based diagnostic used in the present study is useful for understanding model performance with regard to convection-circulation coupling and can also be used to evaluate designs of CPS for MJO simulation.
While our study focuses on how the convection responds to MJO-scale circulations, air-sea coupling is a key process to enhance MJO simulations in the ECAHM5 in previous studies. Tseng et al. (2015) has found by coupling SIT to ECHAM5, i.e. E5SIT-Nord setup in our study, the model can better resolve variation of upper ocean temperature structure, especially in the warm layer from the surface to 10 m. When MJO moves to Indian Ocean, its convection center induces subsidence to its east leading edge, causing clear sky condition and weakening surface winds there. With SIT, surface warm layer is resolved and warmed due to more incident solar radiation associated with the enhanced subsidence and weaker evaporative cooling due to weakened winds. The warmer SST ahead of the MJO further destabilizes the boundary layer and enhances frictional convergence to help organize and strengthen deep convection, which is hypothesized to lead to a more realistic eastward-propagating MJO. More vigorous atmosphere-ocean interaction can thus strengthen intraseasonal variations in both SST and atmospheric circulation. Our preliminary study has shown that with AGCM setup the ECHAM5 models hardly capture the eastward propagating intraseasonal signals, which hampers further study of convection response to MJO-scale circulations (Fig. S8).
Although we focus on the CPS, MJO-related convection-circulation coupling actually depends on the combined effects of moist convection designs, including shallow convection, the PBL, and the parameterization of air-sea interaction. For example, with regard to application in models derived from the Community Atmosphere Model 2, Liu et al. (2005) concluded that the Tiedtke scheme outperforms the Nordeng-Tiedtke scheme because of its ability to trigger convection with low-level convergence. Considering that the ECHAM5 used here had different PBL schemes and shallow convection schemes from Community Atmosphere Model 2, the moisture buildup processes that are involved are likely to differ substantially. The simulations in the present study considered air-sea coupling, which, according to DeMott et al. (2019), are crucial contributors to the low-level moistening observed in successful MJO simulations. Although it is beyond the scope of this study, the short period of convection in the E5SIT-Tied simulation likely allows more shortwave radiation to warm ocean regions, enhancing the SST, and further destabilizes atmosphere, causing positive feedback between convection and SST. This feedback strengthens the convective responses to SST and hinders the coupling of convection with large-scale circulation patterns such as those in the MJO. Therefore, the holistic consideration of moist convection and air-sea coupling is necessary to improve the representation of the MJO in climate models. This premise warrants further investigation.

Conclusion
We conducted experiments using the ECHAM model with one-column high resolution ocean model to investigate the key processes for MJO propagation over the MC in climate simulations. With models with better representation of air-sea coupling, we demonstrated that convection-circulation coupling is crucial to the modeling of MJO propagation over the MC. E5SIT-Tied simulates less moisture buildup over the MC when the deep convection in phases 1 and 2 is necessary for the propagation of the MJO over the MC. The E5SIT-Tied simulation, which does not model MJO propagation, produces local MSE anomaly which is greatly correlated with local deep convection over the MC, in contrast to the shallow convection simulated by E5SIT-Nord. Such rapid development of deep convection over the MC creates an environment that is less conducive to MJO propagation. Differences in low-level moistening in relation to light rain events between the two simulations indicated that E5SIT-Tied tends to remove more moisture through convective processes in the presence of subsidence than does E5SIT-Nord. The hindcast case study of an observed MJO indicated that the environments that are less favorable for the eastward propagation of the MJO in the E5SIT-Tied simulation may be formed through weaker deep convective heating at the convective center of the MJO, which in turn leads to weaker subsidence and moisture buildup over the MC and the western Pacific. The present findings suggest that even with air-sea coupling, MJO propagation over the MC is highly dependent on the capacity of CPS to produce MJO-related convection-circulation coupling, including (1) inducing strong deep convective heating associated with convective center over the Indian Ocean and (2) forming a favorable environment characterized by low-level moistening under the effects of a subsidence branch to the east of convective center. If CPSs are capable of this, convection-circulation coupling is reinforced, and the maintenance of convective processes of the MJO over the MC is facilitated. Our analysis indicates that the instantaneous response of CPS to MJO-related environmental changes can further affect intraseasonal variability in climate models; thus, the response of CPS must be carefully designed. Our analysis may be used as a process-based diagnostic for evaluating convection-circulation coupling in the MJO. Our findings can help model developers resolve the challenge of the "MJO barrier" over the MC, helping them better forecast MJO in their climate models.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.