The El Niño-Southern Oscillation (ENSO) phenomenon is the main source of the predictability skill in many regions of the world for seasonal and interannual timescales. Longer lead predictability experiments of Niño3.4 Index using simple statistical linear models have shown an important skill loss at longer lead times when the targeted season is summer or autumn. We develop different versions of the model substituting some its variables with others that contain tropical or extratropical information, produce a number of hindcasts with these models using two different predictions schemes and cross validate them. We have identified different sets of tropical or extratropical predictors, which can provide useful values of potential skill. We try to find out the sources of the predictability by comparing the sea surface temperature (SST) and heat content (HC) anomalous fields produced by the successful predictors for the 1980–2012 period. We observe that where tropical predictors are used the prediction reproduces only the equatorial characteristics of the warming (cooling). However, where extratropical predictors are included, the predictions are able to simulate the absorbed warming in the South Pacific Convergence Zone (SPCZ).
The El Niño-Southern Oscillation (ENSO) is the dominant interannual
climatic signal and the main source of the predictability skill in many
regions of the world for seasonal and interannual timescales. Improving the
ENSO understanding and forecast skill is still one of the main goals of
international seasonal forecast programs (Kirtman et al., 2013). The ENSO
impacts in many regions of the world opens an opportunity window for forecast
in those regions. The monthly Climate Diagnostics Bulletin (CDB;
The phenomenon El Niño depends on factors that include interactions between the tropics and extratropics that affect its behavior (Philander, 1999). In a previous work (Tasambay-Salazar et al., 2015a) we have determined the feedbacks between ENSO and teleconnection indexes from other regions of the world. In the present work, the Niño3.4 Index characterizes the ENSO variability and the variability of the other regions is represented by some other indexes. Some of these, as the Southern Oscillation Index, or the Pacific Western Equatorial Heat Content Index, representing the equatorial or the tropical Pacific state, have been used in many studies for ENSO forecasts studies (Barnston and Livezey, 1987; Penland and Magorian, 1993; Zebiak and Cane, 1987). Other, that represent the variability outside of the tropical Pacific basin, were suggested or used as predictors in the most recent literature. Among the ENSO precursors external to the tropical Pacific disclosed by recent studies we count some signals in the tropical Indian ocean (Izumo et al., 2010; Luo et al., 2010) or the tropical Atlantic (Keenlyside et al., 2013; Rodríguez-Fonseca et al., 2009). Other novel research like Stepanov (2009); Terray (2011); Vimont et al. (2009) or Ballester et al. (2011) point to the potential of some extratropical variability the North Pacific, the Antartic current at its Indian, Atlantic or Pacific sectors, for ENSO predictability.
In two recent papers, Tasambay-Salazar et al. (2015a, b) have introduced some
of these novel precursors in a series of simple seasonal stochastic model in
order to study the characteristics of the predictability of the equatorial
Pacific, represented by the Niño3.4, the Niño4 and the Niño1
The models employed in Tasambay-Salazar et al. (2015a, b) were built on the basis of some feedback relationships, identified from the lead-lag correlation coefficient values between the corresponding Niño Index and each of the other indexes. In Tasambay-Salazar et al. (2015a) we compared the feedbacks relationship identified in the previous (1950–1979) with those obtained for the post satellite period. In the case of the Niño3.4 Index, the only variable identified with a predictive potential for both periods was the equatorial Pacific anomalous heat content. The performance of the models was assessed by two parameters, correlation skill and root mean square error.
In the present work, however, we present some other measures of the adequacy of the models, like its performance at hindcasting extreme and normal conditions, or at simulating some basic physical variables, like the sea surface temperature anomalies or the heat content anomalies in the tropical Pacific. Details on the datasets used, the procedures followed for the variables identification and the models are given in Sect. 2. The results of the analysis are presented in Sect. 3, and we finish with some concluding remarks in Sect. 4.
The ENSO state that represents the Niño3.4 Index is obtained as an
average of the monthly SST anomalies in the region (5
In addition to the Niño3.4 Index, the models that incorporate
extratropical predictors use indexes obtained from an analysis of covariance
of the seasonal anomalies of some global oceanic and regional atmospheric
fields. The basic oceanic field is the global SST field from ERA-Interim
reanalysis (Dee et al., 2011) at a spatial grid of 0.5
Acronyms and specifications of the regions used to build the extratropical predictors.
Let all variables of the model
In order to further assess the performance of the models, we compare the patterns of the tropical SST and HC fields simulated with the most successful models with the observed ones. These patterns were obtained by correlation the corresponding hindcasted of observed indexes onto those fields. Additionally, we separate the skill obtained when hindcasting El Niño (or La Niña) situations from the skill of the hindcasts issued for ENSO-neutral conditions. In this way, we are able to identify those configurations that have very good skill with the ENSO anomalous phases from those that correspond to normal situations.
Niño3.4, SOI and WWV indexes in summer (Fig. 1a) and Niño3.4, SOI and WWV indexes in autumn (Fig. 1b).
Seasonal cross-correlation of the Niño3.4 Index
hindcast with the NSW model in
Performance models for normal and extreme conditions. The horizontal axis corresponds to normal skill and the vertical axis to the extreme skill. The green straight dashed line indicates the reference for equal values of the hindcast skill corresponding to anomalous ENSO conditions and those for the ENSO-neutral state. The red dashed straight lines represent an arbitrary threshold for the useful forecast as proposed by Hollingsworth et al. (1980). The acronyms correspond to a model configuration in the different regions, as detailed in Table 1.
The Tropical or T-model used here as a benchmark was chosen as the most successful of several tropical models, with different configurations. We use basically three indexes, the Niño3.4 Index and the SOI Index and the WWV Index (NSW model). These Indexes are represented in Fig. 1a for the summer model and in Fig. 1b for the autumn model. Some other models are built by replacing one of the two basic Indexes other than the Niño3.4 index, as for instance, the NSH model (include a index from heat content field), the NST model (include a index from sea surface temperature field), the NSM model (include a index from middle tropospheric field). The PS scores (seasonal cross-correlation coefficients) of these predictions for Niño3.4 Index at each season and for different lead times are presented in Fig. 2a (summer) and in Fig. 2b (autumn), their seasonal cross-RMSE hindcasted (calculated by RMSE of the seasonal cross-correlation hindcasted previously) are depicted in Fig. 2c and d, respectively. The deterioration of the predictability at lead times greater than two seasons in the case of summer and greater than three seasons in the case of autumn is evident there. We can observe that using the more sophisticated FSM prediction scheme produces useful PS values at leads up to two seasons for summer and at leads up to three seasons in the case of autumn. The measure of the forecast skill given by the root mean square error (RMSE) (represented in Fig. 2c and d) is better for FSM prediction scheme, assuming the persistence of the Niño3.4 values, which are represented by a continuous gray line. The inclusion of extra equatorial information (either surface or subsurface) somewhat improves the PS values for the predictions of the summer Niño3.4 Index at longer lead times. We omit the results in the case of the winter and spring predictions, whereas previously mentioned the PS values remain above 0.7 at lead times up to one year.
Some of the models that use extratropical predictors are able to somehow improve the PS values. Their performance is summarized and depicted for the SST and HC summer fields in the Taylor diagrams depicted in Fig. 3a and b, respectively. We can appreciate that the skill correlation coefficients of the models whose variables were obtained from SAIA or RB regions are above the useful threshold. The relevant traits of the patterns of those PCs included in the configurations that scored higher PS values are represented in the Figs. 4 and 5. In the case of the Niño3.4 Index summer predictions, the 14th PC of the SST global field and the 11th PC of the MTT North Pacific field seem to be the only winter variables that are able to produce useful PS values. The corresponding Empirical Orthogonal Functions (EOF) are represented in Fig. 4a and b. In the oceanic pattern, there are cooling anomalies in the Eastern and Central NP and warming anomalies above. There is also a dipole in the Eastern Pacific, with an incipient cooling in the Niño3 region. In the atmospheric MTT pattern, the horseshoe-cooling pattern over land is similar to the upper wind field associated with the seasonal footprinting by Anderson (2004). Some traits of the oceanic pattern combined with the reduced amount of variance explained support the view that these variables should be included in the T-models as noise, as suggested by Penland and Sardeshmukh (1995).
The summer predictors that score the highest PS values are the 15th PC
of the global SST field and the 3rd and 5th PCs of the GSE MTT field. In the
case of the summer MTT variables, the 3rd, and 5th PCs explain 12 and
7 %, respectively, of the total variance. The SST 15th PC and the GSE MTT
3rd PC are also included in most of the configurations found useful
for the autumn prediction. The corresponding SST EOF, in Fig. 5a shows an
anomalous cooling in the South Pacific Convergence Zone. In Fig. 5b, the 3rd
GSE MTT EOF displays wavelike activity in the Antarctic Circumpolar Current (ACC)
region, with three maxima located near the Shackelton Ridge, the Ross
region and over the Patagonia sector and topped by anomalies of opposite
sign at the lower latitudes. In the 5th GSE MTT EOF, depicted in Fig. 5c,
the maxima, which have opposite signs, are displaced approximately
40
The snapshots in the Figs. 6 and 7 represent the forecasted SST (Fig. 6) and HC (Fig. 7) fields. They were constructed by regressing the forecasted time coefficient into the respective fields. The Figs. 6a and 7a represent the SST and HC fields reconstructed from the observations, respectively. We can appreciate in Fig. 6a the characteristic SST anomaly in the Western Pacific, traveling to the east that signals the first stages of an El Niño event. The HC pattern represented in Fig. 7a corresponds also to an El Niño initial stage. The Figs. 6b and 7b presents the SST and HC fields reconstructed with equatorial predictors, and correspond to an already developed El Niño warming. Notice how in the warm anomalies in the HC field are confined to the equatorial region. The patterns in the Figs. 6c and 7c were obtained using NP predictors, and show a clear El Niño warm episode, that is, however less developed that the one obtained with the equatorial model. In the case of the Figs. 6d and 7d, the patterns were obtained using predictors from the South Atlantic and Indian Antarctic sectors and more similar to the ones obtained with the observation than the preceding ones. The Figs. 6e and 7e were produced with predictors obtained through a statistical analysis of the Ross-Bellingshausen region (RB). The phase of the El Niño warming represented by the last patterns corresponds well to the one obtained with the observations. Furthermore, we have separated the skill values of the hindcast that target anomalous ENSO conditions from the ones issued for ENSO-neutral states.
In Fig. 8, we can observe the performance models for normal and extreme conditions. The horizontal axis corresponds to normal skill and the vertical axis to the extreme skill. The green straight dashed line indicates the reference for equal values of the hindcast skill corresponding to anomalous ENSO conditions and those for the ENSO-neutral state. The red dashed straight lines represent an arbitrary threshold for the useful forecast as proposed by Hollingsworth et al. (1980), the acronyms correspond to a model configuration as detailed in Table 1. We can observe that the best hindcast skill value for the extreme and the normal conditions correspond to predictors from SAI region. A possible explanation of this behaviour could be due to the water mass formation that occurs in this region provides important links to the global ocean circulation, allowing variability in the Southern Ocean to affect global climate conditions. Additionally, the ACC system acts as the only direct link connecting the world's major ocean basins, allowing for signals to be transmitted directly among these basins according to Holland et al. (2005) or Verdy et al. (2006). Another hand, the hindcast skill prediction obtained from GSE region achieves similar values for extreme and normal conditions. This applies also to models that use predictors obtained from the NP region, but with no useful skill for normal conditions.
We have identified that best value hindcast skill of extreme conditions was obtained by the model uses the 7th PC in the SAI sector of the MTT field, while the best value hindcast skill of normal conditions is provided by a model that uses the 3rd PC in the SAI sector of the MTT field. In the configurations that use predictors from the SPA and NP sectors, the model uses the 2nd and the 1st PCs of the MTT field. It is important to highlight that all the retained configurations included among their variables the 15th PC of the SST field.
In this work, we investigate the factors that determine the longer lead ENSO predictability. ENSO state is represented by the seasonal value of the Niño3.4 Index. ENSO predictability is measured by the hindcast skill of a number of simple empirical models. The variables used for the model's configuration were identified by a feedback analysis in a previous study. The skill of the hindcasts obtained by models built from equatorial and tropical Pacific variables serve as benchmarks against which all the others are compared.
In the previous study, we found that the low PS values that characterize the ENSO longer lead predictions are mainly a consequence of the low PS values of the predictions that target the summer and autumn Niño3.4 conditions when using only tropical predictors. To improve these PS values, we build a number of models that include the extratropical processes and regions identified as longer lead ENSO precursors by previous studies. Among these precursors are the NP region, the ACC region, the RB dipole or the Southern Atlantic and Indian Oceans. Although some of these precursors were used for predictions in a previous work (Dominiak and Terray, 2005), no study has provided a detailed comparison of their predictive capabilities. A second novelty of the present work is the inclusion of variables derived from the MTT field in the model.
The results obtained there with the T-model indicate that introducing the variability of the tropical ocean or atmosphere into the model yields better the PS values at longer leads for the summer and autumn Niño3.4. Index. They also show that further improvements of the longer lead PS can be achieved by introducing extratropical information in the model variables. The analysis of many sensitivity experiments revealed that the higher PS scores correspond to configurations that include variables representing global SST and regional MTT anomalies. The results also present a dependence on the seasonal signature of the initial conditions, which is in agreement with the “seasonal footprinting” concept.
The results of previous works also shown that there were ensembles of configurations yielding similar values of PS. The present study investigates other measures than the PS that can help us to select one or other configuration as the most appropriate. One of these is the quality of the reconstruction of the key SST and HC tropical anomalies obtained with those models. Among the configuration analyzed the ones that include MTT variables obtained from the RB or from the SAIA sectors are the ones that provide a better representation of those fields conditions in the Tropical Pacific compared with the observed ones. We notice that when using tropical predictors the prediction reproduces only the equatorial characteristics of the warming (cooling) events while when considering extratropical predictors the predictions are able to simulate the warming in the South Pacific Convergence Zone (SPCZ) that appear in the observations.
The other approach followed here consists in separating the PS values obtained for hindcasts of anomalous (positive or negative) ENSO conditions from those where ENSO is in a neutral state. In this context, we found models that achieve very good performances when targeting the anomalous ENSO conditions and rather loosely ones when hindcasting neutral conditions, and other do the opposite. Moreover, there are other models that have a similarly good performance in both situations. To the first class belong models that include variables that are ahead of the predictand when the maximum value of the correlation coefficient is reached, like the 7th PC of the MTT. To the second class belong those models that have good skill at hindcasting normal conditions although the reasons remain yet unexplained. In the third category (the all-terrain) are the models that include variables that have important correlation for both, negative and positive lags with respect summer Niño3.4 Index.
The ERA-Interim SST data used in this paper are available from the European
Centre for Medium-Range Weather Forecasts
(
We thank the supporting of SENESCYT – Ecuador. Some of the figures were produced using the NCAR Command Language (NCL) software package (version 6.1.2.). Edited by: J. D. Pabón-Caicedo Reviewed by: J. J. Velez and one anonymous referee