skip to main content

Title: Extended-range statistical ENSO prediction through operator-theoretic techniques for nonlinear dynamics

Forecasting the El Niño-Southern Oscillation (ENSO) has been a subject of vigorous research due to the important role of the phenomenon in climate dynamics and its worldwide socioeconomic impacts. Over the past decades, numerous models for ENSO prediction have been developed, among which statistical models approximating ENSO evolution by linear dynamics have received significant attention owing to their simplicity and comparable forecast skill to first-principles models at short lead times. Yet, due to highly nonlinear and chaotic dynamics (particularly during ENSO initiation), such models have limited skill for longer-term forecasts beyond half a year. To resolve this limitation, here we employ a new nonparametric statistical approach based on analog forecasting, called kernel analog forecasting (KAF), which avoids assumptions on the underlying dynamics through the use of nonlinear kernel methods for machine learning and dimension reduction of high-dimensional datasets. Through a rigorous connection with Koopman operator theory for dynamical systems, KAF yields statistically optimal predictions of future ENSO states as conditional expectations, given noisy and potentially incomplete data at forecast initialization. Here, using industrial-era Indo-Pacific sea surface temperature (SST) as training data, the method is shown to successfully predict the Niño 3.4 index in a 1998–2017 verification period out to more » a 10-month lead, which corresponds to an increase of 3–8 months (depending on the decade) over a benchmark linear inverse model (LIM), while significantly improving upon the ENSO predictability “spring barrier”. In particular, KAF successfully predicts the historic 2015/16 El Niño at initialization times as early as June 2015, which is comparable to the skill of current dynamical models. An analysis of a 1300-yr control integration of a comprehensive climate model (CCSM4) further demonstrates that the enhanced predictability afforded by KAF holds over potentially much longer leads, extending to 24 months versus 18 months in the benchmark LIM. Probabilistic forecasts for the occurrence of El Niño/La Niña events are also performed and assessed via information-theoretic metrics, showing an improvement of skill over LIM approaches, thus opening an avenue for environmental risk assessment relevant in a variety of contexts.

« less
; ;
Award ID(s):
1854383 1842538
Publication Date:
Journal Name:
Scientific Reports
Nature Publishing Group
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Skillfully predicting the North Atlantic Oscillation (NAO), and the closely related northern annular mode (NAM), on ‘subseasonal’ (weeks to less than a season) timescales is a high priority for operational forecasting centers, because of the NAO’s association with high-impact weather events, particularly during winter. Unfortunately, the relatively fast, weather-related processes dominating total NAO variability are unpredictable beyond about two weeks. On longer timescales, the tropical troposphere and the stratosphere provide some predictability, but they contribute relatively little to total NAO variance. Moreover, subseasonal forecasts are only sporadically skillful, suggesting the practical need to identify the fewer potentially predictable events at the time of forecast. Here we construct an observationally based linear inverse model (LIM) that predicts when, and diagnoses why, subseasonal NAO forecasts will be most skillful. We use the LIM to identify those dynamical modes that, despite capturing only a fraction of overall NAO variability, are largely responsible for extended-range NAO skill. Predictable NAO events stem from the linear superposition of these modes, which represent joint tropical sea-surface temperature-lower stratosphere variability plus a single mode capturing downward propagation from the upper stratosphere. Our method has broad applicability because both the LIM and the state-of-the-art European Centre for Medium-Rangemore »Weather Forecasts Integrated Forecast System (IFS) have higher (and comparable) skill for the same set of predicted high skill forecast events, suggesting that the low-dimensional predictable subspace identified by the LIM is relevant to real-world subseasonal NAO predictions.

    « less
  2. Abstract This study assesses the predictive skill of eight North American Multimodel Ensemble (NMME) models in predicting the Indian Ocean dipole (IOD). We find that the forecasted ensemble-mean IOD–El Niño–Southern Oscillation (ENSO) relationship deteriorates away from the observed relationship with increasing lead time, which might be one reason that limits the IOD predictive skill in coupled models. We are able to improve the IOD predictive skill using a recently developed stochastic dynamical model (SDM) forced by forecasted ENSO conditions. The results are consistent with the previous result that operational IOD predictability beyond persistence at lead times beyond one season is mostly controlled by ENSO predictability and the signal-to-noise ratio of the Indo-Pacific climate system. The multimodel ensemble (MME) investigated here is found to be of superior skill compared to each individual model at most lead times. Importantly, the skill of the SDM IOD predictions forced with forecasted ENSO conditions were either similar or better than those of the MME IOD forecasts. Moreover, the SDM forced with observed ENSO conditions exhibits significantly higher IOD prediction skill than the MME at longer lead times, suggesting the large potential skill increase that could be achieved by improving operational ENSO forecasts. We find thatmore »both cold and warm biases of the predicted Niño-3.4 index may cause false alarms of negative and positive IOD events, respectively, in NMME models. Many false alarms for IOD forecasts at lead times longer than one season in the original forecasts disappear or are significantly reduced in the SDM forced by forecasted ENSO conditions.« less
  3. null (Ed.)
    Abstract Using a high-resolution atmospheric general circulation model simulation of unprecedented ensemble size, we examine potential predictability of monthly anomalies under El Niño Southern Oscillation (ENSO) forcing and back-ground internal variability. This study reveals the pronounced month-to-month evolution of both the ENSO forcing signal and internal variability. Internal variance in upper-level geopotential height decreases (∼ 10%) over the North Pacific during El Niño as the westerly jet extends eastward, allowing forced signals to account for a greater fraction of the total variability, and leading to increased potential predictability. We identify February and March of El Niño years as the most predictable months using a signal-to-noise analysis. In contrast, December, a month typically included in teleconnection studies, shows little-to-no potential predictability. We show that the seasonal evolution of SST forcing and variability leads to significant signal-to-noise relationships that can be directly linked to both upper-level and surface variable predictability for a given month. The stark changes in forced response, internal variability, and thus signal-to-noise across an ENSO season indicate that subseasonal fields should be used to diagnose potential predictability over North America associated with ENSO teleconnections. Using surface air temperature and precipitation as examples, this study provides motivation to pursue ‘windowsmore »of forecast opportunity’, in which statistical skill can be developed, tested, and leveraged to determine times and regions in which this skill may be elevated.« less
  4. Abstract

    The mean-state bias and the associated forecast errors of the El Niño–Southern Oscillation (ENSO) are investigated in a suite of 2-yr-lead retrospective forecasts conducted with the Community Earth System Model, version 1, for 1954–2015. The equatorial Pacific cold tongue in the forecasts is too strong and extends excessively westward due to a combination of the model’s inherent climatological bias, initialization imbalance, and errors in initial ocean data. The forecasts show a stronger cold tongue bias in the first year than that inherent to the model due to the imbalance between initial subsurface oceanic states and model dynamics. The cold tongue bias affects not only the pattern and amplitude but also the duration of ENSO in the forecasts by altering ocean–atmosphere feedbacks. The predicted sea surface temperature anomalies related to ENSO extend to the far western equatorial Pacific during boreal summer when the cold tongue bias is strong, and the predicted ENSO anomalies are too weak in the central-eastern equatorial Pacific. The forecast errors of pattern and amplitude subsequently lead to errors in ENSO phase transition by affecting the amplitude of the negative thermocline feedback in the equatorial Pacific and tropical interbasin adjustments during the mature phase of ENSO. Thesemore »ENSO forecast errors further degrade the predictions of wintertime atmospheric teleconnections, land surface air temperature, and rainfall anomalies over the Northern Hemisphere. These mean-state and ENSO forecast biases are more pronounced in forecasts initialized in boreal spring–summer than other seasons due to the seasonal intensification of the Bjerknes feedback.

    « less
  5. Abstract Current climate models have relatively high skills in predicting El Niño–Southern Oscillation (ENSO) phase (i.e., El Niño, neutral, and La Niña), once leaping over the spring predictability barrier. However, it is still a big challenge to realistically forecast the ENSO amplitude, for instance, whether a predicted event will be strong, moderate, or weak. Here we demonstrate that the accumulated westerly wind events (WWEs)/easterly wind surges (EWSs) and oceanic recharged/discharged states are both of importance in accurate ENSO amplitude forecasts. El Niño and La Niña events exhibit asymmetric temporal and spatial features in the atmospheric and oceanic preconditions. El Niño amplitude at the peak season is closely associated with the accumulated WWEs over the eastern equatorial Pacific from the previous December to May and the recharged state in the western equatorial Pacific during February. In contrast, the amplitude of La Niña events is sensitive to the accumulated EWSs over the equatorial western Pacific from the previous November to April and the discharged state extending from the equatorial western to central Pacific during February. Considering these asymmetric atmospheric and oceanic preconditions of El Niño and La Niña cases, a statistical model is established to accurately forecast the ENSO amplitude at itsmore »mature phase during 1982–2018, which is validated to be robust based on a 1-yr cross-validation and independent sample tests. The feasibility and the limitation of the established statistical model are also discussed by examining its practical utility.« less