How a Trail in the Sand Revealed El Niño's Future ✎

Original: "Discovery of a Phenomenological Dynamical Model for Predicting the El Niño-Southern Oscillation"
· Shivsai Ajit Dixit, B N Goswami
arXiv:1709.02606 · 2017-09-08 · arXiv non-exclusive · independent analysis · 1 min
Scientists found a simple rule in 135 years of weather records that predicts El Niño better than complex computer models.

Imagine finding a trail of footprints on the beach. From them, you can tell whether the animal walked calmly or ran, and even predict where it will head next. Scientists did something similar with El Niño, the warm Pacific Ocean current that affects weather across the globe. Instead of building a massive ocean-atmosphere model, they took 135 years of observational records and tried to read the rule of its motion from them.

First, they filtered out rapid fluctuations, leaving only slow changes. Then they represented each moment as a point in a space where the coordinates are the record values on different days. It turned out that all points lie almost on a plane, not scattered randomly. This means there is a simple rule among the values, which can be written as a single formula. The formula is nonlinear: it accounts for the fact that warm and cold phases of El Niño behave differently.

The rule was trained on just 40 years of data and tested on the remaining 94 years. It predicts El Niño 6 months ahead with an accuracy of 0.96, while the best modern models achieve about 0.86.

Previously, El Niño forecasts deteriorated sharply every spring, a phenomenon known as the spring barrier. The new rule greatly weakens this barrier. This means we can now prepare in advance for droughts, floods, and other consequences of El Niño. And it all runs on an ordinary computer in seconds, not a supercomputer.

The model is trained on just 486 points (40.5 years of data) but forecasts 94.5 years ahead. Moreover, the computational cost of forecasting is practically zero, unlike the supercomputer calculations of physical models.

y_{i+2\tau} = A + B(y_{i+\tau} - y_{i-\tau}) + C(y_{i+\tau} - y_{i-\tau}) y_i
The series value two time intervals τ ahead is expressed in terms of previous values
\frac{d^3\hat{y}}{d\hat{t}^3} + \frac{P_2}{P_1}\frac{d^2\hat{y}}{d\hat{t}^2} + \frac{P_3}{P_1}\frac{d\hat{y}}{d\hat{t}} + \frac{P_4}{P_1}\hat{y} = 0
Equivalent ODE describing the dynamics of the slow manifold
y_{i+\tau} = A + B y_i + C y_{i-\tau}
Linear analog used as a baseline
Concepts
spring predictability barrier Jin-Timmermann model El Niño-Southern Oscillation ENSO Multivariate ENSO Index El Niño index recharge oscillator model equatorial wave propagation ERSST La Niña Southern Oscillation natural variability relative ONI TAO/TRITON moorings seasonal cycle stochastic dynamical system forecast NOAA CPC NOAA PMEL / TAO Westerly wind burst