Propagating a Storm: How the TEAMS Algorithm Grows Rare El Niños ✛

Original: "Using a rare event sampling technique to quantify extreme El Niño event statistics"
arXiv:2512.23038 · 2025-12-28 · CC BY 4.0 · 3 min
The TEAMS algorithm reproduces the statistics of extreme El Niños in the Zebiak-Cane model at one-fifth the computational cost of direct simulation.

In 1997, the Pacific Ocean warmed so much that corals bleached from Indonesia to the Galápagos, and floods and droughts swept across five continents. Such events, which climate scientists call extreme El Niños, occur every few decades, and to understand how often they will arrive in the future, we need thousands of years of observations. Direct climate simulation over such timescales requires computational power comparable to forecasting weather a hundred thousand years ahead. The authors of a new study propose a workaround: instead of waiting for a rare event, they propagate its trajectories, much as a gardener propagates a rare plant from cuttings.

The TEAMS algorithm (Transition-Enhanced Adaptive Multilevel Splitting) works like a plant breeder. It starts with an ensemble of 16 short simulations, each 25 years long, and at each step retains only the quarter of trajectories with the highest temperature peaks in the eastern tropical Pacific. From these “parents,” it generates “offspring” by adding a tiny random perturbation to the sea surface temperature (amplitude 0.1°C) 7.5 years before the threshold crossing. The process repeats 50 times, and the total cost of all runs is 100,000 model years, five times less than the benchmark direct simulation of 500,000 years.

The key trick is probability weights. Each offspring inherits the parent's weight divided by the number of offspring, so the final statistics remain unbiased even though the sample is artificially enriched with extreme events.

Comparison with the benchmark showed that TEAMS reproduces return time estimates for extreme El Niños up to 60,000 years. The distribution of block maxima, weighted by branching probabilities, matched the benchmark well, especially in the tail where the strongest events reside. Interestingly, the temperature threshold grew rapidly in the first generations but then leveled off at about 4.2°C, suggesting a possible internal limit to extremity in the model. The success rate of generating offspring that exceed the threshold was about 77%, and trajectory analysis showed that a small perturbation can first trigger a La Niña and then amplify the subsequent El Niño through non-normal growth of perturbations.

The TEAMS cost of 100,000 model years is equivalent to simulating climate from 1997 to about the year 102,000 CE, but achieved by branching just 50 short runs.

Before this work, it was unclear whether the statistics of extreme El Niños could be obtained without prohibitively long simulations. Now it is shown that the algorithm reproduces the distribution and return times in the Zebiak-Cane model at a cost five times lower than the benchmark. The trick is simple and elegant: branching trajectories with small perturbations produces offspring that more often reach extreme thresholds, and probability weights correct for sampling bias. The method can be transferred to expensive global climate models to assess the risks of rare El Niños today.

Propagating a storm from a single trajectory, which began with cuttings in 1997, now yields a whole family of extreme events ready for statistical analysis. Instead of waiting for a rare guest to knock on the door, we have learned to invite it hundreds of times and record every visit.

TEAMS' cost of 100,000 model years is equivalent to simulating climate from 1997 to about the year 102,000 CE, but achieved by branching only 50 short runs.

The idea of branching trajectories to explore rare outcomes echoes the multiverse in Ted Chiang's 'Story of Your Life', where different branches of the future are viewed simultaneously.

\tau(\ell) = -\frac{T}{\log[1-p_T(\ell)]}
The return time for level ℓ equals the block length T divided by the negative logarithm of the probability of not exceeding the threshold in a block.
\mathrm{CDF}(\theta) = \frac{1}{M}\sum_{m=1}^{M}\frac{1}{N}\sum_{j=1}^{J_m} w_{m,j} \mathbb{I}\{\theta_k \le \theta\}
The cumulative distribution function is computed as the average over M runs of weighted indicators of threshold exceedance.
Scientists
Jascha Sohl-Dickstein
Concepts
Multivariate ENSO Index El Niño zebiak cane model Jin-Timmermann model ENSO El Niño-Southern Oscillation spring predictability barrier ERSST TAO/TRITON moorings ensemble forecast NOAA PMEL / TAO relative ONI rare events event selection Warm water volume rare event search Oceanic Niño index IRI