Climate variability in the tropical Indo-Pacific sector has undergone dramatic changes under global ocean warming. Extreme Indian Ocean dipole (IOD) events occurred repeatedly in recent decades with an unprecedented series of three consecutive episodes during 2006–08, causing vast climate and socioeconomic effects worldwide and weakening the historic El Niño–Indian monsoon relationship. Major attention has been paid to the El Niño influence on the Indian Ocean, but how the IOD influences El Niño and its predictability remained an important issue to be understood. On the basis of various forecast experiments activating and suppressing air–sea coupling in the individual tropical ocean basins using a state-of-the-art coupled ocean–atmosphere model with demonstrated predictive capability, the present study shows that the extreme IOD plays a key role in driving the 1994 pseudo–El Niño, in contrast with traditional El Niño theory. The pseudo–El Niño is more frequently observed in recent decades, coincident with a weakened atmospheric Walker circulation in response to anthropogenic forcing. The study’s results suggest that extreme IOD may significantly enhance El Niño and its onset forecast, which has being a long-standing challenge, and El Niño in turn enhances IOD and its long-range predictability. The intrinsic El Niño–IOD interaction found here provides hope for enhanced prediction skill of both of these climate modes, and it sheds new light on the tropical climate variations and their changes under the influence of global warming.
El Niño and extreme Indian Ocean dipole (IOD) are two dominant drivers for year-to-year climate variability on earth. Predicting those climate modes is of great value because of their large environmental and societal effects, both globally and regionally. El Niño is now generally predictable at a lead time of several seasons (e.g., Palmer et al. 2004; Luo et al. 2005a; Saha et al. 2006; Jin et al. 2008), and it may be predicted even up to two years in advance for specific events (Chen et al. 2004; Luo et al. 2008b). However, predicting the exact onset of El Niño remains a long-standing challenge, for either the strongest 1997/98 episode or the recent weak-to-moderate 2006/07 event (e.g., Barnston et al. 1999; Landsea and Knaff 2000; McPhaden 2006; see also the real-time multimodel ENSO forecasts at http://iri.columbia.edu/climate/ENSO/currentinfo/QuickLook.html). This is consistent with the observational finding (e.g., Kessler 2002); unlike its La Niña counterpart, the El Niño onset does not follow the theoretical self-sustained oscillatory behavior (e.g., Jin 1997; Neelin et al. 1998). External forcing is often necessary to initiate El Niño onsets.
Considerable observational, model, and theoretical efforts were made to understand the influences of intraseasonal disturbances, especially westerly wind bursts in the equatorial western Pacific, on the El Niño onset and growth (e.g., McPhaden 1999; Fedorov et al. 2003; Gebbie et al. 2007). The importance of stochastic forcing was recognized particularly after monitoring the evolution of the very strong 1997/98 El Niño, in which westerly wind bursts apparently contributed to its rapid development [e.g., McPhaden 1999; see also the Pacific Triangle Tropical Atmosphere–Ocean/Trans-Ocean Buoy Network (TAO–TRITON) buoy observations online at http://www.pmel.noaa.gov/tao/jsdisplay/]. The unpredictable intraseasonal wind forcing, however, was believed to impose a fundamental limit to El Niño predictability (e.g., Flügel et al. 2004); a pessimistic result for El Niño forecasters. Recent model studies suggested that the Indian Ocean climate variability might have an important influence on the Pacific El Niño–Southern Oscillation (ENSO) variance and frequency (e.g., Yu et al. 2002; Wu and Kirtman 2004; Behera et al. 2006), though their results appear to be inconsistent and model dependent. Under certain circumstances, surface wind anomalies in the western Pacific related to the Indian Ocean basin-wide cooling or warming, as a delayed response to La Niña or El Niño, may in turn affect the El Niño onset or decay during specific seasons (e.g., Kug et al. 2005, 2006). The potential contribution of just this type of Indian Ocean signal to ENSO predictability, however, appears to be limited because the source for the predictability of the former is essentially determined by ENSO itself.
In addition to the basin-wide signal forced by ENSO predominantly through a fast atmospheric Walker circulation adjustment (e.g., Klein et al. 1999), internal air–sea coupling in the tropical Indian Ocean generates another important climate mode, called IOD (e.g., Saji et al. 1999; Webster et al. 1999). Extreme IOD events characterized by strong cooling in the eastern Indian Ocean (EIO) can induce large convective diabatic heating anomalies and therefore affect the equatorial atmospheric Walker circulation and ENSO evolution. El Niño and IOD repeatedly co-occurred in recent observations since the mid-1970s (e.g., Annamalai et al. 2005), indicating the interactive nature of the two major climate modes. We note that, unlike the ENSO-induced basin-wide signal in the Indian Ocean, IOD has its own dynamics and independent source for its predictability (e.g., Behera et al. 2006; Luo et al. 2007, 2008a). Thus, exploring the IOD influence on El Niño would provide additional space for our better understanding and prediction of El Niño.
In this study, we examine potential contributions of the Indo-Pacific interbasin coupling to the predictions of both El Niño and extreme IOD, in particular their onsets. This is not just because both climate modes have large climate and socioeconomic effects worldwide and hence enormous societal benefits from their predictions. Improved prediction of either mode could enhance the predictive skill of the other in terms of their interacting feedbacks. By conducting various forecast experiments with a fully coupled global ocean–atmosphere model, we demonstrate that extreme IOD has important influence on El Niño and its predictability, and vice versa. Thus the interbasin coupling via the atmospheric Walker circulation is crucial to both El Niño and IOD evolutions and predictions. The coupled model and seasonal prediction experiments are described in section 2. Retrospective forecast and “perfect model” prediction results for both El Niño and extreme IOD are presented in sections 3 and 4. Summary and discussions are given in section 5.
2. Model, retrospective forecasts, and perfect model predictions
a. The model
The numerical model used in the present study is the Scale Interaction Experiment–Frontier Research Center for Global Change (SINTEX-F) fully coupled global ocean–atmosphere general circulation model (Luo et al. 2003, 2005b; Masson et al. 2005), which was developed at Frontier Research Center for Global Change under the European Union–Japan collaboration (Gualdi et al. 2003; Guilyardi et al. 2003). The atmospheric component (ECHAM4.6) has the resolution of 1.1° × 1.1° (T106) with 19 vertical levels (Roeckner et al. 1996). The oceanic component [Ocean Parallelise 8 (OPA8)] has a relatively coarse resolution with a 2° Mercator horizontal mesh and 31 layers in vertical (Madec et al. 1998). Its meridional resolution was increased to 0.5° near the equator to properly capture the equatorial wave dynamics. Those two components are coupled every two hours without any flux corrections using a standardized interface (Valcke et al. 2000). The SINTEX-F model has been applied to various climate studies, and it has shown good performance in simulating and predicting both ENSO and IOD (e.g., Yamagata et al. 2004; Luo et al. 2005a,b; Behera et al. 2005, 2006; Luo et al. 2007, 2008a,b; Jin et al. 2008). Here, we refer to IOD as its positive phase (strong cooling in the eastern Indian Ocean and weak warming in the west). We note that a large asymmetry exists between the surface warming and cooling intensity in the eastern Indian Ocean associated with the negative and positive IOD events (Hong et al. 2008a,b). Negative IOD events do not appear to evolve into strong air–sea coupled processes in the Indian Ocean, and therefore their peak magnitudes are weak with lower predictability in general (Luo et al. 2007).
b. Retrospective forecast experiments
Nine-member ensemble retrospective forecasts for 12 target months from the first day of each month during 1982–2006 were performed based on a semi-multimodel ensemble approach; both model coupling physics and initial conditions were perturbed separately in three different ways [see Luo et al. (2005a, 2007) for detailed descriptions]. A simple but effective way was adopted to produce realistic and well-balanced ocean–atmosphere initial conditions by assimilating only satellite-observed sea surface temperatures (SSTs) into the coupled model (Luo et al. 2005a). The well-balanced initial conditions are important to achieving skillful long-lead ENSO prediction out to two years ahead (Chen et al. 2004; Luo et al. 2008b). Interannual oceanic thermocline variations in the equatorial Pacific, the source of memory on which useful seasonal predictions of ENSO could be based, were reproduced well in general using the coupled SST-nudging initialization approach (see Luo et al. 2005a, 2007). However, significant errors exist for specific seasons (cf. left and right panels in Fig. 1), particularly when intraseasonal Kelvin waves passed. This suggests the importance of assimilating oceanic subsurface information to improve the model initial conditions. The same model initial conditions for the retrospective forecasts were used for the two sensitivity forecast experiments described below.
The east–west Walker circulation in the equatorial atmosphere, particularly the easterly (westerly) trade winds in the Pacific (Indian Ocean) driven by the strong convection over the warm waters in the Maritime Continent–western Pacific, could be regulated by and play an important role in the redistribution of the warm waters in the Indo-Pacific region, in association with ENSO and IOD evolution. The influence of the interbasin coupling on the predictions of El Niño and extreme IOD was examined by suppressing air–sea coupling in the tropical Indian Ocean (hereafter referred to as dIO) and Pacific Ocean (dPO), respectively. Monthly-mean climatological SSTs, based on the satellite observations during 1983–2006 (Reynolds et al. 2002), were prescribed in the individual ocean basins between 25°S and 25°N.1 By doing this over the basin without air–sea coupling, the atmosphere there will respond to the climatological observed SST, rather than predicted SST, during the 12-month forecast period. That is, oceanic feedbacks to the atmosphere there are suppressed, and thereby the interacting feedbacks between the Indian Ocean and Pacific SST variations via the atmospheric bridge (particularly the Walker circulation) are suppressed.2
c. Perfect model prediction experiments
We first integrated the free coupled model for 520 years, starting from Levitus annual mean climatology without motion (Luo et al. 2005b). The first 20 years had a rapid spinup in global mean SST and thus were discarded. Within the remaining 500-yr simulation, 25 co-occurrences of El Niño and extreme IOD events occurred. Consistent with the observational definitions, El Niño in the model was defined when the boreal winter (December–February) mean SST anomaly in the Niño-3.4 region (5°S–5°N, 120°–170°W) exceeds 0.5°C. An average SST anomaly lower than −1°C in the EIO (10°S–0°, 90°–110°E) during September–November was used as the criterion for the extreme IOD definition.
To assess the interbasin coupling effects on the El Niño and extreme IOD predictions, the same air–sea decoupling approach was adopted except that the 500-yr climatological model SST was prescribed in the decoupled basins. The same initial conditions produced by the free coupled model simulation were used in the two perfect model experiments for predicting each of the 25 El Niño and extreme IOD events. The predictions were conducted in the same manner as was done in the retrospective forecast experiments for 12 target months, starting from the first day of each month during July of the year before the onsets of El Niño and IOD to December of the year when the El Niño and IOD co-occurred.
3. Retrospective forecast results
a. Influence of extreme IOD on El Niño
Over last two decades, one pseudo–El Niño in 1994 and two El Niño in 1997 and 2006 appeared concurrently with extreme IOD. This combination caused severe climate effects over the globe, including the exacerbated drought in Southeast Asia–Australia and floods in East Africa in 1994 and 2006. Figure 2 shows our model retrospective forecasts of the three El Niño with and without the Indian Ocean influence (see section 2b for the methods). Apparently, predictions of these El Niño onsets appear to completely fail in the latter case (cf. the left and right panels in Fig. 2). The pseudo-El Niño in 1994/95 is characterized by major surface warming in the central Pacific rather than in the east during canonical El Niño years. Because of this, its associated global climate effects were found to be distinct from those of typical El Niño (e.g., Weng et al. 2007). The warming in the central Pacific was also found to be more effective in forcing drought over India and Australia (Kumar et al. 2006; Wang and Hendon 2007). Interestingly, along the equatorial Pacific thermocline, no significant warm signal preceded its onset in October 1994, except for a weak and shallow warming confined in the central area (see Fig. 1, top panels). This suggests the necessity of external forcing for its growth. The extreme IOD event that started early in 1994 appears to have suppressed the development of a La Niña signal in the eastern Pacific and played a key role in driving the evolution of this pseudo-El Niño (the green and orange lines in Figs. 2a and 2b; see also Fig. 4a).3
The SINTEX-F model shows a limited success in predicting the 1997/98 El Niño onset (Fig. 2c, green lines), despite the substantial influence of the unpredictable intraseasonal westerly wind bursts from late 1996 to early 1997 (e.g., McPhaden 1999; see also the Pacific TAO–TRITON buoy observations online at http://www.pmel.noaa.gov/tao/jsdisplay/) and cold biases in model initial subsurface conditions in the equatorial Pacific (Fig. 1, middle panels). Without the Indian Ocean influence, its onset predictions initiated from 1 October 1996 to 1 March 1997 basically fail, and the predicted peak magnitudes are below 0.5°C (see the green lines in Fig. 2d). Model forecasts from April to May 1997 and later months, however, are not much affected by the Indian Ocean. This is probably because the surface and subsurface warming in the eastern Pacific had already become sufficiently strong in April–May 1997 (Figs. 1d and 1j) to initiate unstable air–sea coupled processes in the Pacific for the El Niño growth.
Existing models with a variety of complexities failed to predict the weak-to-moderate 2006/07 El Niño event (McPhaden 2006; see also the real-time multimodel ENSO forecasts online at http://iri.columbia.edu/climate/ENSO/currentinfo/QuickLook.html), possibly owing to the influences of westerly wind bursts in the western Pacific and intraseasonal oceanic Kelvin waves during the course of its evolution (e.g., Luo 2007; McPhaden 2008). Despite this difficulty, our model correctly predicts the phase transition from La Niña to El Niño in 2006 boreal summer but underestimates its peak magnitude in November–December 2006 [the green and orange lines in Fig. 2e; see also Luo et al. (2007) for our real-time forecasts], presumably as a result of the cold biases in the model initial subsurface conditions (Fig. 1, bottom panels). Again, suppressing the Indian Ocean effect leads to a long-lasting La Niña throughout 2006 and a complete failure in predicting the El Niño onset (Fig. 2f). This is related to the influence of cold SST anomalies (though small) in the eastern Indian Ocean from the boreal spring to the fall of 2006 in the model predictions (see the green lines in Fig. 6e; possible reasons are discussed later).
The underlying mechanisms for enhanced El Niño prediction by resolving the Indian Ocean signal can be understood in the following way: Early in the year of the El Niño onset, colder than normal SSTs throughout the tropical Indian Ocean sometimes exist as a result of atmospheric response to the precedent La Niña in the Pacific (e.g., Klein et al. 1999). This induces anomalous westerlies in the western Pacific, which force an eastward-propagating equatorial downwelling (warm) Kelvin wave along the thermocline in the Pacific (Figs. 3a and 3b; see also Kug et al. 2005). This tends to terminate the La Niña condition and initiate surface warming in the eastern Pacific, which sometimes may trigger unstable air–sea interactions or Bjerknes feedback (Bjerknes 1969) responsible for the El Niño development. The El Niño growth is further enhanced by a subsequent warm Kelvin wave in boreal summer, which is driven by the anomalous westerlies related to the surface cooling and downdraft air over the eastern pole of the growing IOD (Figs. 3a and 3b).
Starting from the middle of the El Niño onset year, well-developed cooling in the eastern Indian Ocean enhances westerly winds in the western central Pacific (e.g., Saji and Yamagata 2003; Annamalai et al. 2005), thus deepening the central eastern Pacific thermocline through an eastward-propagating warm Kelvin wave (Figs. 3c and 3d). This eventually accelerates the growth of El Niño or pseudo–El Niño, particularly in 1994 and 2006 (Fig. 4; see also the orange lines in Fig. 2) when the local Bjerknes feedback in the Pacific is not strong. In the absence of subsurface warming in the east, the major surface warming in the central Pacific associated with the 1994/95 pseudo–El Niño is related to the eastward displacement of warm pool via anomalous surface current advection (Fig. 4a; see also Picaut et al. 1997) driven by the IOD-induced anomalous westerlies in the western Pacific. The anomalous westerlies also tend to weaken the equatorial upwelling and reduce the entrainment of subsurface cold water into the surface layer in the equatorial western and central Pacific (not shown).
In the latter part of the El Niño onset year, contributions from the extreme IOD events to the El Niño growth become generally small (Fig. 2, blue lines). This is not in agreement with previous studies (e.g., Saji and Yamagata 2003; Annamalai et al. 2005), in which they speculated that the peak IOD signal in boreal fall might have a large influence on El Niño growth. Contributions of the external forcing to El Niño growth in the latter part of the year appear to be essentially constrained by the intrinsic phase locking of ENSO to the annual cycle (e.g., Rasmusson and Carpenter 1982). We note that the IOD influence on the El Niño growth in 1997 is much limited after the El Niño onset in boreal summer (see the orange lines in Figs. 2c, 2d and 4b). This is probably because the local Bjerknes feedback in the Pacific after the El Niño onset might be already strong enough to support the rapid growth of the 1997/98 El Niño, and that the model forecast of the 1997 IOD is less successful (see the orange lines in Fig. 6c, discussed later). The El Niño decay during the following boreal spring, however, appears to be slightly delayed without the Indian Ocean influence (Fig. 2, blue lines). This is particularly true for the 1998 case, attributed to the simultaneous strong basin-wide surface warming in the Indian Ocean (e.g., Kug et al. 2006).
To further examine the influence of the Indian Ocean on El Niño prediction, we performed additional forecast experiments of other El Niño—1982/83, 1987/88, 1991/92, and 2002/03—with and without the Indian Ocean influence (Fig. 5). In the absence of extreme IOD influence, the prediction differences of the four El Niño between the two experiments appear to be much smaller compared to those of the three El Niño in 1994/95, 1997/98, and 2006/07. The results support the present finding that extreme IOD may have an important influence on El Niño development. Without the Indian Ocean influence, a systematic underestimation of the intensity of the 1982/83 El Niño can also be seen (Figs. 5a and 5b, green and orange lines). This is probably related to the effect of an overpredicted IOD in 1982 (see Luo et al. 2007). We note that model prediction of the El Niño onset in 1986 appears to be enhanced by the Indian Ocean basin-wide cooling in late 1985 (not shown); the latter is closely related to the long-lasting La Niña condition in the Pacific during late 1983 to early 1986. This is also consistent with the present and existing studies (e.g., Kug et al. 2005, 2006).
b. Influence of ENSO on extreme IOD
Low-level westerly winds blow in the equatorial Indian Ocean throughout most of the year. The winds drive surface eastward currents along the equator, particularly the semiannual Wyrtki jets (Wyrtki 1973), pushing the surface water eastward and thereby forming a deeper thermocline and warm water pool in the east; this is a unique structure compared to the Pacific and Atlantic Oceans (CLIVAR–Goos Indian Ocean Panel et al. 2006). Because of the deep thermocline there, air–sea feedbacks in the Indian Ocean are generally weak compared to ENSO-related air–sea coupling in the Pacific. However, strong Bjerknes feedback may occur during boreal summer when the southeasterly monsoon winds prevail along the west coast of Sumatra and raise the coastal thermocline to be shallow enough—a prerequisite condition—for unstable growth of extreme cooling in the eastern Indian Ocean.
Figure 6 shows our model forecasts of the past three extreme IOD events with and without the Pacific influence (see section 2b for the methods). At long lead times (see the green lines in left panels of Fig. 6), the model shows a limited success in forecasting the extreme IOD, especially in 1994 and 2006.4 When the Pacific influence is suppressed, surprisingly the model’s long-lead predictions appear to be improved (see the green lines in the right panels of Fig. 6). This is particularly the case for the peak intensities in 1994 and 2006 boreal fall. The model experimental results, however, do not contradict the common notion that ENSO has large influences on the Indian Ocean (e.g., Klein et al. 1999). In fact, the underestimated peak intensity of the extreme IOD in the model retrospective forecasts is primarily caused by a systematic La Niña–like model climate drift in the Pacific (Figs. 7a and 8a; see also Luo et al. 2005a, 2008b). This model cold bias, via the atmospheric Walker circulation, induces excessive westerly winds in the Indian Ocean, which hamper the shoaling of thermocline in the east and thereby damp the IOD growth (Figs. 7a and 7b). This suggests the importance of reducing the model climate drift in the Pacific during forecast. Present results also indicate that a part of long-lead predictability of extreme IOD resides in the Indian Ocean itself, in support of our previous studies (Luo et al. 2007, 2008a). At mid-lead times, the Pacific El Niño influence on the extreme IOD growth appears to be rather limited (orange lines in Figs. 6 and 7c), indicating a dominant control by the local Bjerknes feedback over the Indian Ocean region.
In summary, our model forecast experiments show that the extreme IOD has important influences on El Niño evolution and its onset predictability; and in specific years it may induce a different kind of El Niño (i.e., pseudo–El Niño) with distinctive climate effects. Besides, the Indo-Pacific interbasin coupling appears to have substantial effects on the onset and evolution of both El Niño and extreme IOD via the atmospheric Walker circulation. To further confirm this, we have conducted two perfect model prediction experiments using a 500-yr coupled model simulation (see section 2c for the methods). A total of 25 cases in which extreme IOD co-occurred with El Niño are adopted (Fig. 9). The evolution and spatial patterns of the two climate modes are reproduced realistically in the model.
4. Perfect model prediction results
Applicability of the results based on the retrospective forecast experiments (see the previous section) might be limited because of the few co-occurrences of El Niño and extreme IOD in the short observational records. In addition, model biases and errors in initial conditions reduce the predictive skill of El Niño and IOD; and this may lead to some unrealistic estimations for the effects of the interbasin coupling. In the perfect model prediction experiments, errors in both model physics and initial conditions were assumed to be removed. It is worth noting that, although intrinsic model biases might still affect the estimation, some confidence arises from the model’s good performance in simulating and predicting the tropical climate (e.g., Yamagata et al. 2004; Luo et al. 2005a,b; Behera et al. 2005, 2006; Luo et al. 2007, 2008a,b; Jin et al. 2008).
Apparently, the evolution of both El Niño and extreme IOD cannot be correctly reproduced if the interbasin coupling is suppressed, particularly during their onset from April to July (gray and green lines in Figs. 10a and 10b). Thus, contributions from the internal air–sea coupled processes in the individual ocean basins are not enough to support their rapid onset and growth. It is worth noting that the internal precursor for the limited long-range predictability of the extreme IOD in the model (gray and green lines in Fig. 10b) comes from the strong anomalous subsurface cooling in the southwestern Indian Ocean (Fig. 11). The cold subsurface anomaly in the southwestern Indian Ocean in boreal winter, which is a response to the cyclonic-like wind forcing in the east about one season ahead (Fig. 11a; see also Xie et al. 2002; Rao et al. 2002), will propagate westward to the western boundary and then reflect as eastward-propagating equatorial upwelling Kelvin waves in the following seasons (Fig. 11b; see also Rao et al. 2002; Behera et al. 2006), providing a long-lead precursor for the extreme IOD evolution. This is consistent with recent observations and model forecast experiments (Horii et al. 2008; Luo et al. 2007, 2008a).5
The underlying physical mechanisms for the interbasin coupling in the perfect model experiments are essentially consistent with those found in the retrospective forecast experiments. Weak basin-wide surface cooling in boreal winter and strong cooling in the eastern Indian Ocean in the following seasons associated with IOD development induce anomalous westerly winds in the western Pacific, which tend to drive two pulses of equatorial warm Kelvin waves in the Pacific during boreal winter and summer, enhancing El Niño onset and growth significantly (Figs. 10c and 10d). SST differences in the equatorial Pacific enhanced by the Indian Ocean signal reach as much as 0.6°–0.8°C in late year (Fig. 10c). In turn, the growing El Niño in the Pacific induces anomalous easterly winds in the Indian Ocean during boreal spring via the atmospheric Walker circulation, raising the thermocline west of Sumatra and thereby reinforcing the onset and growth of IOD significantly (Figs. 10e and 10f). The resemblance of the appearance of these anomaly patterns in the equatorial Indo-Pacific region (cf. the middle and bottom panels in Fig. 10) indicates an intrinsic interaction between El Niño and extreme IOD via the Walker circulation. For the predictions initiated after their onsets, the interbasin coupling also has significant but makes smaller contributions to their ultimate magnitudes (orange lines in Figs. 10a, 10b and 12), in accordance with the retrospective forecast experiment results (Figs. 3 and 7).
5. Summary and discussion
Both El Niño and extreme IOD have large environmental and societal effects over the globe. Therefore, improving the forecasts of El Niño and extreme IOD, particularly their onsets, has enormous societal benefits. The important influence of ENSO on the Indian Ocean climate (particularly the basin-wide mode) has been well recognized over the past decades. However, little has been known for the IOD influence on El Niño and its predictability. Our results based on the retrospective forecast and “perfect model” prediction experiments suggest that extreme IOD has significant contributions to El Niño onset and its long-lead predictability, and hence it may have large indirect climate effects worldwide in addition to its directly related ones. It was found that the IOD also acts to weaken the El Niño–induced drought over India (e.g., Kumar et al. 1999; Ashok et al. 2004). Our results indicate that El Niño and extreme IOD intimately interact with each other via the atmospheric east–west Walker circulation. The indo-Pacific interbasin coupling is crucial to the evolution of both El Niño and extreme IOD and their predictions at long lead times. After their onsets, however, contributions of the interbasin coupling to their subsequent growth become limited, owing to the dominant role of the local Bjerknes feedback in the individual ocean basins.
In contrast to the pessimism that stresses the role of stochastic forcing in the ENSO predictability (e.g., Flügel et al. 2004), the results presented here provide some optimism: Improving the seasonal forecast of the predictable Indian Ocean climate variability may eventually lead to more skillful El Niño forecasts, and vice versa. The net gain in skill comes from the interbasin coupling and independent sources for the predictability of the two climate modes. Thus, current efforts in establishing a long-term monitoring system in the Indian Ocean (CLIVAR–GOOS Indian Ocean Panel et al. 2006) will contribute to better understanding and prediction of not only the Indian Ocean climate but also ENSO in the Pacific. It is surprising that the El Niño–like signal can be fully generated by extreme IOD as in 1994, in contrast with classical ENSO theory [see Neelin et al. (1998) for a review]. Better understanding of how El Niño and IOD might evolve and influence each other under global warming may have important implications for the future projection of the climate on earth. Noticing the more frequent occurrences of extreme IOD and pseudo–El Niño in recent decades (e.g., Trenberth and Stepaniak 2001; Ashok et al. 2007), probably in association with the weakened Walker circulation in response to anthropogenic forcing (Vecchi et al. 2006) and decadal El Niño–like variability in the Pacific (e.g., Luo and Yamagata 2001), it is conceivable that the intensified IOD activity (Ihara et al. 2008; Abram et al. 2008) would play a more important role in El Niño evolution in the present global warming trend. This may have implications for our future projection of ENSO under the influence of global warming.
Satellite SST observations were provided by NOAA/OAR/ESRL (available online at http://www.cdc.noaa.gov/cdc/data.noaa.oisst.v2.html). All model experiments were carried out on the JAMSTEC Earth Simulator. We thank S.-P. Xie for the helpful discussions and two anonymous reviewers for their valuable comments, which helped to improve the manuscript.
* Additional affiliation: Department of Earth and Planetary Science, University of Tokyo, Tokyo, Japan.
Corresponding author address: Jing-Jia Luo, Research Institute for Global Change, JAMSTEC, 3173-25 Showa-machi, Kanazawa-ku, Yokohama, Kanagawa 236-0001, Japan. Email: firstname.lastname@example.org
A “sponge layer” of 10° width in longitude and 5° in latitude, where weights of the free air–sea coupling increase linearly toward the outside ocean was applied to avoid unrealistic instability near the boundary. Climate variations in the South China Sea, influenced by and may in turn influence ENSO and IOD, were not suppressed in the present study [see also Behera et al. (2006) for the free coupled model sensitivity experiments].
We note that, using a simple two-predictor regression model based on the equatorial Pacific warm water volume and the Madden–Julian oscillation forcing in the western Pacific, McPhaden et al. (2006) achieved good skill in predicting most of the ENSO events during 1981–2005. This linear model, however, completely failed to predict the pseudo–El Niño in 1994/95 (see their Fig. 3). This appears to support our finding here.
The prediction of the 1997 IOD was less successful owing to systematic model biases and bad initial conditions (Luo et al. 2005a, 2007). Besides, the warm intraseasonal event in 1997 boreal spring caused warm subsurface initial conditions, leading to false alarms of a negative IOD for the predictions at midlead times (orange lines in Figs. 6c and 6d).
It is interesting to note that one recent study based on a similar perfect model experiment using the Geophysical Fluid Dynamics Laboratory (GFDL) coupled model (Song et al. 2008) found that predictable IOD in that model is largely preconditioned by the tropical Pacific subsurface signal, while the internal subsurface memory in the Indian Ocean plays a little role. Further model intercomparison study might be helpful to understand the discrepancy.