An empirical model of decadal ENSO variability

This paper assesses potential predictability of decadal variations in the El Ni?o/Southern Oscillation (ENSO) characteristics by constructing and performing simulations using an empirical nonlinear stochastic model of an ENSO index. The model employs decomposition of global sea-surface temperature (SST) anomalies into the modes that maximize the ratio of interdecadal-to-subdecadal SST variance to define low-frequency predictors called the canonical variates (CVs). When the whole available SST time series is so processed, the leading canonical variate (CV-1) is found to be well correlated with the area-averaged SST time series which exhibits a non-uniform warming trend, while the next two (CV-2 and CV-3) describe secular variability arguably associated with a combination of Atlantic Multidecadal Oscillation (AMO) and Pacific Decadal Oscillation (PDO) signals. The corresponding ENSO model that uses either all three (CVs 1–3) or only AMO/PDO-related (CVs 2 and 3) predictors captures well the observed autocorrelation function, probability density function, seasonal dependence of ENSO, and, most importantly, the observed interdecadal modulation of ENSO variance. The latter modulation, and its dependence on CVs, is shown to be inconsistent with the null hypothesis of random decadal ENSO variations simulated by multivariate linear inverse models. Cross-validated hindcasts of ENSO variance suggest a potential useful skill at decadal lead times. These findings thus argue that decadal modulations of ENSO variability may be predictable subject to our ability to forecast AMO/PDO-type climate modes; the latter forecasts may need to be based on simulations of dynamical models, rather than on a purely statistical scheme as in the present paper.