How to Find an El Niño Precursor Faster ✎

Original: "The Sampling Method for Optimal Precursors of ENSO Events"
· Bin Shi, Junjie Ma
arXiv:2308.13830 · 2023-08-26 · CC BY-NC-SA 4.0 · independent analysis · 1 min
Random sampling instead of complex calculations finds the same warm patch in the Pacific Ocean five times faster.

A blindfolded sapper searches for a mine by tossing pebbles at random. He can't see where the charge is buried. But if you toss many pebbles and listen for the one that lands with a duller thud, you can locate the danger. That's how researchers hunt for the precursor of El Niño, the warm current that alters weather across the entire planet.

El Niño begins with a warm patch in the Pacific Ocean. If the patch grows quickly, drought hits Australia and rains hit Peru a few months later. The goal is to find that initial patch.

Previously, this required building an inverse model. It's like a translator who back-translates a speech into the original language to check every word. That's difficult and time-consuming. The new method is simpler: the computer repeatedly tweaks the initial conditions and watches how fast the patch grows. Each trial is like a tossed pebble. From hundreds of trials, you can infer where the 'mine' is without ever seeing it. The calculation takes 3 seconds instead of 15, five times faster, and the result is nearly identical.

The sampling method is based on Stokes' theorem, which relates the integral over a volume to the integral over its boundary. Here, it allows replacing the gradient with an average of function values on a sphere: an elegant mathematical trick.

In Isaac Asimov's story 'The Last Question', a supercomputer searches for a way to reverse entropy by trying possibilities. The sampling method also tries random directions, but instead of an exhaustive search, it uses statistics to quickly find the optimal perturbation.

\nabla \hat{J}(u_0) = \frac{d}{\epsilon} \mathbb{E}_{v \in S^{d-1}} [J(u_0 + \epsilon v) v]
The gradient of the smoothed function equals the average over the sphere of the product of the function value and the direction, multiplied by the dimension and the inverse radius.
\nabla \hat{J}(u_0) \approx \frac{d}{n \epsilon} \sum_{i=1}^n [J(u_0 + \epsilon v_i) v_i]
The gradient is estimated as the arithmetic mean over n random directions, which avoids computing the exact gradient.
Concepts
zebiak cane model ENSO El Niño index Jin-Timmermann model spring predictability barrier Multivariate ENSO Index El Niño-Southern Oscillation Southern Oscillation NOAA CPC TAO/TRITON moorings La Niña ERSST NOAA PMEL / TAO relative ONI equatorial wave propagation forecast North Equatorial Countercurrent Warm water volume Monte Carlo methods ensemble forecast