A central Indian Ocean (CIO) mode is found to play a critical role in driving the heavy precipitation during the Indian summer monsoon (ISM). It is typically denoted with a combination of intraseasonal sea surface temperature (SST) anomalies and intraseasonal wind anomalies over the central Indian Ocean, and it preserves the mechanistic links among various dynamic and thermodynamic fields. Like a T junction, it controls the propagation direction of the intraseasonal variabilities (ISVs) originating in the western Indian Ocean. During the ISM, the CIO mode creates an environment favorable for the northward-propagating mesoscale variabilities. These results unveil the relation between the subseasonal monsoonal precipitation and the CIO mode in the ocean–atmosphere system in the Indian Ocean. The identification of the CIO mode deepens our understanding of the coupled monsoon system and brightens the prospects for better simulation and prediction of monsoonal precipitation in the affected countries.
Precipitation during the Indian summer monsoon (ISM) is the lifeline for billions of people living on the rim of the Indian Ocean. The monsoonal precipitation is dominated by intraseasonal variabilities (ISVs), which are commonly known as the monsoon intraseasonal oscillation (MISO; Goswami 2005; Shukla 2014). For example, the variance of the intraseasonal precipitation accounts for about 60% of the total variance of precipitation over the Bay of Bengal (BoB). A comprehensive description of MISO is given by Annamalai and Slingo (2001), Waliser (2006), Annamalai and Sperber (2016), and Jones (2016). The MISO can be represented with a multivariate index that was established by Lee et al. (2013) and Suhas et al. (2013). The MISO is believed to be closely related to the northward-propagating ISVs. Yasunari (1980) was one of the pioneers to demonstrate that MISO was likely to originate from the ISVs over the equatorial Indian Ocean. Wang and Rui (1990) categorized the ISVs into three groups and found that the northward-propagating ISVs prevailed during boreal summer. Several mechanisms have been proposed for MISO. Webster (1983) hypothesized that the increased surface heat flux ahead of the convection over the continent destabilized the atmosphere, leading to the northward propagation. However, in reality, the ISVs mainly propagate northward over the ocean, not over land. Wang and Xie (1997) argued that the apparent northward propagation of ISVs was actually a westward-tilting convection front, which emanated from the tropical Pacific as Rossby waves. This hypothesis was supported by the evidence that many northward-propagating ISVs coincided with the eastward-propagating ones. However, the independent northward-propagating ISVs, as categorized in Wang and Rui (1990), cannot be explained. Furthermore, the convective momentum transport (CMT; Kang et al. 2010) and the interaction between the perturbations and the background easterly shear (Jiang et al. 2004) were both found to be important for MISO. Using an idealized land–atmosphere coupled model, Bellon (2011) showed that the small heat capacity of the land greatly damped the ISVs over the Indian subcontinent, which could explain the small intraseasonal variance over land during ISM. Recently, the moisture effects during tropical convection has received much attention. For example, Kemball-Cook and Wang (2001) showed that a relation exists between the moist Kelvin–Rossby coupled waves and MISO; Abhik et al. (2013) argued that the atmospheric instability is attributable to the moisture convergence due to the low-level advection and subgrid eddies. Ocean–atmosphere interaction is another important mechanism since the warm sea surface temperature (SST) anomalies can cause moisture convergence in the atmospheric boundary layer, destabilizing the atmosphere to the north of convection (Kemball-Cook and Wang 2001; Vecchi and Harrison 2002; Fu and Wang 2004; Krishnan et al. 2011). In addition, the ocean is able to play an active role during MISO (Yano and McBride 1998; Zhou and Murtugudde 2014). In spite of the differences in details between the various mechanisms, one consensus is that the MISO stems from the tropical ISVs, which originally propagate eastward from the western equatorial Indian Ocean. It is still unclear what mechanism causes such tropical ISVs to take a “left turn” and propagate northward during ISM. This question is addressed in this study.
In addition, the relationship between the monsoonal precipitation and the intrinsic low-frequency modes in the coupled ocean–atmosphere system, mainly El Niño–Southern Oscillation (ENSO; e.g., Kumar et al. 1999; Krishnamurthy and Kirtman 2003; Gadgil et al. 2004; Kumar et al. 2006; Gill et al. 2015) and the Indian Ocean dipole/zonal mode (IODZM; e.g., Allan et al. 2001; Ashok et al. 2001; Abram et al. 2008; Sabeerali et al. 2012; Sabeerali et al. 2014), has been widely studied on the subseasonal-to-interannual and longer time scales. Although the variabilities on the intraseasonal time scale are pronounced during ISM, no mode is yet identified for the ISVs. In this study, a relation between the MISO and an intrinsic mode, the central Indian Ocean (CIO) mode, in the atmosphere–ocean system is established at intraseasonal time scales. Diagnostic results show that the CIO mode captures the mechanistic links among various dynamic and thermodynamic fields and plays an important role in influencing the propagation direction of the ISVs during the ISM.
The atmospheric variables such as wind, precipitation, and relative humidity are obtained from the daily National Centers for Environmental Prediction (NCEP)–National Center for Atmospheric Research (NCAR) reanalysis (Kalnay et al. 1996), from 1982 to 2014. The atmospheric variables are also obtained from daily ERA-Interim data, the global atmospheric reanalysis produced by the European Centre for Medium-Range Weather Forecasts (ECMWF; Uppala et al. 2005) for the verification of the robustness of the results. The robustness of precipitation from a reanalysis product may be a concern. However, by comparing with the Tropical Rainfall Measuring Mission (TRMM) rainfall data, it can be shown that the NCEP–NCAR reanalysis has a reasonable rendition of both the total and the intraseasonal monsoonal rainfall. SST data are obtained from the NOAA ¼° daily Optimum Interpolation SST (Reynolds et al. 2007), from 1982 to 2014 and the outgoing longwave radiation (OLR) data are obtained from the NOAA satellite data (Liebmann and Smith 1996), from 1982 to 2013. All intraseasonal variabilities are obtained with a 20–100-day bandpass Butterworth filter.
3. Phenomenon of the central Indian Ocean mode
The mechanism of MISO is examined by contrasting the composites of two groups categorized with respect to the intraseasonal precipitation anomalies, which are the averages within 80°–95°E and 15°–24°N during the traditional ISM period from June to September (Fig. 1). One group consists of all intraseasonal precipitation anomalies that are larger than their mean (the mean of ISVs is almost zero; the solid black line in Fig. 1) plus their standard deviation (STD; the top dashed line in Fig. 1), which is referred to as the large-rainfall group hereafter, and the other group consists of all intraseasonal precipitation anomalies smaller than the mean minus the STD (the bottom dashed line in Fig. 1), which is referred to as the small-rainfall group hereafter. The differences in the intraseasonal precipitation during ISM between the two groups are shown in Fig. 2a. The significant differences mainly occur over the northern BoB, where the monsoonal precipitation is pronounced. The differences in the intraseasonal SST anomalies and intraseasonal wind anomalies in the lower troposphere at 850 hPa between the two groups are shown in Fig. 3a. During the days with large intraseasonal precipitation, there are warm SST anomalies over the central Indian Ocean along the equator. The negative SST anomalies occur off the Somali coast and cross the Arabian Sea and BoB. For the days with small intraseasonal precipitation, the ISV patterns are largely opposite to what are shown in Fig. 3a (i.e., negative SST anomalies in the central Indian Ocean and positive SST anomalies crossing the Arabian Sea and BoB). Corresponding to these SST anomalies, there is an anticyclonic gyre in the lower troposphere, as shown with black vectors in Fig. 3a. This anticyclonic gyre is similar to the mean atmospheric general circulation during ISM (Webster et al. 1998; Goswami and Mohan 2001). However, there are fundamental differences in the physical mechanism. First, the two gyres are separated in time scales; the former is in the intraseasonal band, while the latter is in the seasonal band that provides a background state for the former. Second, the structures and dynamics of the air–sea interaction are different. Seasonal mean easterlies occur in the Southern Hemisphere along the isotherms on the edge of the warm pool over the eastern Indian Ocean, connecting to the Findlater jet (Schott and McCreary 2001) and contributing to the energy transfer from the Southern Hemisphere to the Northern Hemisphere during the ISM (Schott and McCreary 2001; Schneider et al. 2014). In contrast, as shown in Figs. 3b,c, the intraseasonal anticyclone is mainly attributable to the conservation of the potential vorticity (PV; Hoskins et al. 1985; Zhang and Ling 2012). Zhang and Ling (2012) showed that PV can represent the major properties of the ISVs. In the isobaric (pressure) coordinate, PV is defined as P = −g(fk + ∇ × V) ⋅ ∇θ, where g is the gravitational acceleration; f is the Coriolis parameter; V = ui + υj + ωk, where u, υ, and ω are the eastward, northward, and vertical pressure velocities, respectively; i, j, and k are the unit vectors pointing to the east, north, and vertical, respectively; θ is the air potential temperature; and ∇ = i∂/∂x + j∂/∂y + k∂/∂p. Usually, since the vertical component of P dominates, only the vertical component of PV is used in this study, which is calculated as −g(f + ζ)∂θ/∂p, where ζ = ∂υ/∂x + ∂u/∂y is the relative vorticity of the air. PV is generally conserved unless it is modified by diabatic heating or friction. Given the PV conservation, the distribution of PV on an isobaric surface can diagnostically deduce most dynamical fields, such as the horizontal wind anomalies (Hoskins et al. 1985). In the lower troposphere (Fig. 3b), the easterly wind anomalies enhanced by the underlying warm SST anomalies along the equator turn northward around 70°E, roughly following the contours of PV. In the northern BoB, there is a positive PV center, which is consistent with the cyclone there. In the upper troposphere at 300 hPa (Fig. 3c), the anticyclone (cyclone) in the southern (northern) BoB is coherently constrained by the conservation of PV since the diabatic heating and friction, which are the major source and sink of PV, are also small in the upper troposphere.
The above structures associated with MISO represent an intrinsic mode in the Indian Ocean, which incorporates both the atmospheric and oceanic properties. The empirical orthogonal function (EOF) analysis is conducted in the Indian Ocean to isolate this intrinsic mode. A combined field of the daily intraseasonal SST data and daily intraseasonal zonal winds at 850 hPa from 1982 to 2014 is employed for the EOF analysis. The domain is from 40° to 120°E and from 20°N to 20°S, which is what is commonly used (e.g., Murtugudde and Busalacchi 1999; Saji et al. 1999; Murtugudde et al. 2000). Each field is normalized by subtracting its regional mean and then divided by its regional variance within the above domain before doing the EOF analysis, so that the SST and winds can be compared and combined. It has been well known that the first EOF mode of the SST anomalies in the Indian Ocean is the basin mode (Yang et al. 2007, 2009) and the second EOF mode is the IODZM (Saji et al. 1999; Li et al. 2003; Saji and Yamagata 2003; Annamalai and Murtugudde 2004), which has been widely studied since its discovery. The EOF analysis was also conducted using the ISVs in OLR, zonal winds in the lower troposphere, and the multivariate fields (Goswami and Mohan 2001; Wheeler and Hendon 2004; Lee et al. 2013; Suhas et al. 2013; Shukla 2014). As shown in Fig. 3d, the atmospheric and oceanic features of the MISO are captured by the first combined EOF mode of intraseasonal SST anomalies and intraseasonal zonal winds at 850 hPa, which explains 13.6% of the total variance. The positive SST node of the combined EOF mode captures warm intraseasonal SST anomalies along the central Indian Ocean. Correspondingly, the negative wind node denoting the equatorial easterlies and the positive wind node denoting the westerlies from the Arabian Sea to the South China Sea (SCS) indicate an anticyclonic gyre between the tropics and subtropics, which is consistent with the low-level winds in Fig. 3a (black vectors). Hereafter, we refer to the first EOF mode shown in Fig. 3d as the CIO mode, since the key feature and a major node of this mode is the consistent structure between the intraseasonal dynamic and thermodynamic fields over the central Indian Ocean. The positive (negative) CIO mode is defined when there are pronounced warm (cool) SST anomalies and an anticyclone (cyclone) over the central Indian Ocean and the first principal component (PC1) is significantly positive (negative). The differences in composite intraseasonal precipitation between the positive CIO mode and the negative CIO mode are shown in Fig. 2b, which are similar to the differences between the large-rainfall and the small-rainfall groups shown in Fig. 2a. The northward propagation associated with the CIO mode is shown with the differences in the intraseasonal precipitation between the positive and the negative CIO mode phases (Fig. 4). Day 0 is the day when PC1 reaches its local maximum (minimum) for the positive (negative) CIO mode, and negative days in Fig. 4 indicate the days before day 0. In addition, the composites of intraseasonal SST anomalies and intraseasonal wind anomalies for the positive (negative) CIO mode are also similar to the ones for the large (small) rainfall group (not shown). Overall, the properties associated with the CIO mode well capture the main features of the MISO.
4. Dynamics of the central Indian Ocean mode
The downdraft over the central Indian Ocean associated with deep convection over the western Pacific warm pool is the trigger for the intrinsic CIO mode. The clear sky over the central Indian Ocean due to the downdraft increases the incident solar radiation at the sea surface, leading to warm SST anomalies (i.e., the positive CIO mode). During the ISM, the warm pool (where the SSTs are larger than 28°C) extends from the equator to 20°N over the western Pacific Ocean (Fig. 5a). Correspondingly, the mean ITCZ centers around 10°N (Fig. 5b). However, the synoptic convection does not follow the mean warm SST exactly, since deep convection is bounded by the large-scale dynamical subsidence (Takayabu et al. 2010). The differences in midtropospheric intraseasonal specific humidity and the differences in intraseasonal OLR anomalies, both of which are good indicators of tropical convection between the large-rainfall group and small-rainfall group, straddle the equator (Fig. 5). The differences in the intraseasonal zonal and vertical momentum between the two groups are shown in Fig. 6. Because of the air density ρ differences in the upper and lower troposphere, the momentum (ρu, ρω) rather than (u, ω), where u and ω are the zonal and vertical velocities, is used in Fig. 6. Deep convection over the western Pacific is stronger for the large-rainfall group. Correspondingly, the downdraft occurs from central to eastern Indian Ocean between 70° and 120°E. The variability in vertical motion can be clearly detected from the vertical profile of specific humidity, which is critical for atmospheric instability and tropical convection. As a result, for the positive CIO mode, specific humidity increases in the midtroposphere over the tropical western Pacific because of the local enhanced deep convection, while the central Indian Ocean becomes dry because of the downdraft associated with this deep convection. Note that the zonal convection cell shown in Fig. 6 is separated from the well-known Walker cell on both temporal and spatial scales. The convection cell associated with the subseasonal CIO mode is significantly smaller in its zonal extent than the Walker cell, which is closely related to ENSO, IODZM, and the Indo-Pacific tripole (Chen and Cane 2008; Chen 2011; Lian et al. 2014) that dominate the interannual variability over the whole globe. The multiscale nature of the tropical convection system has been explored before (Hendon and Liebmann 1994; Biello and Majda 2005; Tulich and Mapes 2008; Kikuchi and Wang 2010), but here we capture for the first time the intrinsic mode that steers the poleward propagation originating in the Indian Ocean on intraseasonal time scales. The updraft branch of Walker circulation occurs over the Indo-Pacific warm pool ranging from the eastern Indian Ocean to the western Pacific Ocean (roughly from 90° to 160°E, albeit a little weaker over the Maritime Continent) and the downdraft on the Indian Ocean side mainly occurs over the western Indian Ocean to the west of 60°E, which is largely consistent with the mean SST distribution shown in Fig. 5a. In contrast, the convection cell associated with the CIO mode (Fig. 6) is narrower and is mainly consistent with the synoptic convection from the western Pacific to the central Indian Ocean. In addition, the convection cell with the CIO mode is very much local and does not have a branch over the Pacific Ocean as the Walker cell does.
Over the central tropical Indian Ocean, the mean SSTs during ISM are actually slightly cooler than the adjoining areas (Fig. 5a). However, the intraseasonal warm SST anomalies are pronounced during the positive CIO mode and create a favorable environment for MISO. The different relative SST anomalies between the central Indian Ocean and the adjacent regions indicate a separation between the subseasonal scale and the seasonal scale, which is probably due to the ocean–atmosphere interaction dynamics associated with the CIO mode. In reality, the origin of ISVs usually resides in the western Indian Ocean (Seo and Kim 2003; Zhao et al. 2013). There is a year-around shallow convection and a moderate increase of humidity in the lower troposphere, which is critical for the initialization of the Madden–Julian oscillation (MJO), which is the major component of eastward-propagating ISVs in the tropics (Zhang 2005). Convection (denoted with the negative OLR anomalies) has a very similar statistical distribution as shown in Fig. 7, which indicates that convection over the western Indian Ocean can always occur following a very similar statistical distribution, although the mechanism for the generation of convection over the western Indian Ocean is still debated (Kemball-Cook and Weare 2001; Matthews 2008; Straub 2013). In boreal winter, the ISVs (mainly MJO) mostly propagate eastward across the central and eastern Indian Ocean. But during ISM, most ISVs veer northward, which is attributed to the atmospheric circulation modified by the CIO mode. Although the downdraft leads to anticyclones in both the upper and the lower troposphere, the wind anomalies actually enhance the easterly vertical wind shear. Usually, the vertical wind shear is calculated using the wind in the upper atmosphere minus the wind in the lower atmosphere. The easterly wind shear is favorable for the northward propagation of ISV via the CMT (Kang et al. 2010). The downdraft over the central Indian Ocean brings easterly momentum from the upper troposphere to the lower troposphere. Hence, an eastward Coriolis force is required to balance the westward CMT forcing. As a result, a secondary meridional convection cell is induced, with a northward (southward) branch in the lower (upper) troposphere (see Fig. 3 in Kang et al. 2010). Similarly, using a two-layer idealized model, Jiang et al. (2004) showed that the barotropic vorticity and the baroclinic divergence are coupled by the easterly vertical wind shear. The baroclinic divergence, which coincides with the convection center, can induce a positive barotropic vorticity to the north, which consequently enhances the convergence in the planetary boundary layer (PBL) and increases the moisture transport from PBL to the free troposphere to the north of convection, inducing the northward propagation of the ISVs. According to both mechanisms, the mesoscale atmospheric circulation modified by the warm SST anomalies shifts the ISVs originating in the western Indian Ocean to the north, leading to MISO and the supply of momentum and moisture for the onset and heavy precipitation during ISM (Fu et al. 2003; Zhou and Murtugudde 2014). To show the consistency between the warm SST anomalies associated with the positive CIO mode and the precipitation over BoB, we define a convection event over the western Indian Ocean when the consecutive days with intraseasonal OLR anomalies lower than −25 W m−2 last for 10 days with fewer than 3 break days (the day with an intraseasonal OLR anomaly larger than −25 W m−2). The scatterplot of the intraseasonal SST in the central Indian Ocean (averaged between 5°N and 5°S and between 80° and 90°E) and precipitation anomalies over the northern BoB (averaged between 15° and 24°N and between 85° and 95°E) after the convection events are shown in Fig. 8a. All anomalies in Fig. 8a are averaged between 5 and 10 days after the central day of the convection events, which is a reasonable period for the ISVs propagating from the western Indian Ocean to the central Indian Ocean. The correlation coefficient is 0.43 (statistically significant at a 97% confidence level), which indicates that warm intraseasonal SST anomalies of the positive CIO mode are necessary to steer the ISVs originating from the western Indian Ocean to the north and lead to heavy monsoonal rainfall over BoB. The linear relation is robust and is not sensitive to the time period used for the averaging. Note that the effective sample size (ESS) is reduced by the bandpass filtering. When performing the significance test on the correlation coefficient, the ESS is adjusted following Bretherton et al. (1999), so that the influence of high autocorrelation of each bandpass-filtered time series is removed. The CIO mode represents the consistent structure between the dynamic and thermodynamic fields in the ocean and the atmosphere. Thus, it has a better linear correlation (the correlation coefficient is 0.69 at a 99% confidence level) with the monsoonal precipitation as shown in Fig. 8b. In addition, the correlation coefficients between the principal component of CIO mode and the intraseasonal monsoonal precipitation can reach a maximum of 0.8 over the northern BoB (Fig. 8b), which is much larger than the monsoon correlations with the common indices for ENSO and IODZM (Fig. 9). ENSO is denoted with the trans-Niño index (TNI; Trenberth 1997; Trenberth and Stepaniak 2001; Fig. 9c) and Niño-3.4 index (Fig. 9d). IODZM is denoted with the dipole model index (DMI; Saji et al. 1999; Fig. 9b). Usually, the indices for ENSO and IODZM are made with monthly data. But, in order to retain the ISVs in all indices, daily data are used for the calculations and no smoothing in time is applied in this study. Recently, considerable efforts have been made to improve the prediction of the ISM rainfall (e.g., Webster and Hoyos 2004; Abhilash et al. 2013). Although some successes have been achieved, the predictability is still very limited, partly as a result of the weak correlations between the ISM monsoon rainfall and the traditional weather and climate indices (Webster et al. 1998; Wang et al. 2009, 2015). Thus, the large correlations of the CIO mode imply a potential advance in predictive understanding of MISO as well.
The transit between the positive and negative CIO mode occurs as a result of the ocean–atmosphere interactions over the eastern Indian Ocean. The composite Hovmöller diagram of intraseasonal SST anomalies for the positive CIO mode is shown in Fig. 10. The intraseasonal SST anomalies are averaged between 5°N and 5°S. Day 0 is the day with heavy precipitation (i.e., the day when the intraseasonal rainfall is larger than the mean plus the STD). There are 644 such days from 1982 to 2014. Warm SST anomalies occur about 5 days before the heavy precipitation over BoB (Xi et al. 2015). The warm SST anomalies continue to increase since the easterly winds, the equatorial arm of the anticyclone (black vectors in Fig. 3a), continue to transport warm water from the east to the west. In addition, the westerly mean winds during the ISM are weakened by the easterly anomalies, which reduce the latent heat flux and contribute to the warm SST anomalies. Note that Fig. 10 shows the average in the deep tropics between 5°S and 5°N. Along the equator, usually the eastward-propagating Kelvin waves dominate (Valsala 2008; Zhou and Murtugudde 2010). However, the westward propagation is obvious in Fig. 10, which is probably driven by the easterly wind anomalies. Around day 10, negative SST anomalies begin to occur off Sumatra, which is likely to be generated by the coastal upwelling. Then, the cold water is brought to the west also by the easterlies to finally erase the warm anomalies over the central Indian Ocean. Such negative feedback is likely to control the life cycle of the CIO mode, which depends on the background state. If there is a relatively thick layer of warm water in some years off Sumatra, it will take relatively longer for the dissipation of warm SSTs and the CIO phase transition to occur. On the contrary, if the warm water layer off Sumatra is thin, it will be eroded faster. This indicates an irregularity of the CIO mode and the low-frequency variability (such as the interannual variability) of the CIO mode. Moreover, the strong ocean mixing of surface and subsurface waters over the central Indian Ocean may also play a role during the CIO mode, as shown with observations by Pujiana et al. (2015). Most present-day oceanic reanalysis products are not good enough to capture the ISVs in the ocean interior (Zhang et al. 2016). Thus, a solid conclusion about the relation between the oceanic processes and the CIO mode requires more observations and detailed analysis in the future. In addition, the eastward deepening of the thermocline in the Indian Ocean, which is quite unlike the equatorial Pacific and Atlantic Oceans, is well known (Murtugudde and Busalacchi 1999). The interannual variability of the thermocline is of interest in the context of IODZM, but interactions of these lower frequencies with the intraseasonal time scale can occur through the CIO mode. This would offer a new avenue to investigate the interannual and lower-frequency variabilities of MISO and the monsoon itself, likely including the pathways for IODZM and ENSO influences on the monsoon.
5. Conclusions and discussion
In this study, we find an intrinsic mode in the Indian Ocean, the CIO mode, which is closely related to the heavy precipitation during the Indian summer monsoon season. The warm intraseasonal SST anomalies and the associated anticyclone are the key features of the positive CIO mode, which modifies the atmospheric circulation over the Indo-Pacific region and creates a favorable environment for MISO. The CIO mode is like a T junction that controls whether the ISVs can propagate northward into BoB. The entire sequence of processes is summarized in Fig. 11. The CIO mode has a high simultaneous correlation with MISO during ISM, which is expected to contribute to improve the predictive understanding of MISO. It is necessary to make comprehensive evaluations of the CIO mode contribution to the MISO prediction in the future. In addition, the oceanic responses and feedbacks to the CIO mode require further study, probably using an ocean–atmosphere coupled model.
The CIO mode is a piece of the puzzle in the regional characteristics of the global ITCZ variability. On annual and longer time scales (including the geological epochs), the ITCZ usually locates in the warm hemisphere, maintaining the atmospheric energy balance (Schneider et al. 2014). On seasonal time scale, the energetic constraint is still operational. However, the jump of the ITCZ from the Southern Hemisphere to the Northern Hemisphere during ISM is an exception (Fig. 12a), and several mechanisms have been proposed to address this sudden migration (Webster et al. 1998; Chao and Chen 2001; Bordoni and Schneider 2008; Boos and Kuang 2010; Wu et al. 2012). On the subseasonal time scale, the seasonal mean position of ITCZ does not change much and becomes the background state (black line in Fig. 12b). The CIO mode (blue line in Fig. 12b) leads to the northward-propagating MISO (color shading in Fig. 12b), which can overshoot the seasonal mean ITCZ, transferring moisture and momentum to far-northern locations. In the multiscale but “seamless” climate system, the CIO mode is certainly influenced by the background states, such as the ITCZ, convection over the Pacific warm pool, and ENSO, and it may also have feedbacks to the background state. In a multiscale system, the coexistence of the various spatiotemporal scales cannot always be neatly deconvolved, but the existence of the CIO mode points to its role in this complex multiscale system. Its scale selection and its relation to the other members of this important MJO–monsoon–ENSO system require further studies. Nonetheless, the phenology and mechanism of CIO mode we report here are advancements in the process understanding of the monsoon intraseasonal variability and should eventually lead to predictive understanding as well with further analysis of the CIO mode and its interactions with and feedbacks to the large-scale and lower-frequency modes of coupled ocean–atmosphere interactions.
This work is supported by grants from the National Natural Science Foundation of China (41376034, 41321004, and 41690121), the National Basic Research Program (2013CB430302), and the IPOVAR Project (GASI-IPOVAI-01-02 and GASI-IPOVAI-02). RM gratefully acknowledges the CYGNSS grant from NASA and the National Monsoon Mission funds for partial support.