A Ceiling of Randomness ✛ ✎

Original: "Bound on forecasting skill for models of North Atlantic tropical cyclone counts"
arXiv:2410.05503 · 2024-10-07 · arXiv non-exclusive · independent analysis · 3 min
Models of Atlantic cyclone counts have already hit the lower bound of error set by the data's own nature.

Each Atlantic hurricane season is like a game of dice, where the number of spots depends on three major climate forces: El Niño, ocean temperature in the 'hurricane nursery,' and the North Atlantic Oscillation. These three factors explain roughly half of the year-to-year variability in cyclone counts, and for a long time it seemed that adding more data or cleverer algorithms could sharpen the forecasts. But a new study shows that models have already hit a ceiling imposed by the very nature of randomness.

The researchers took cyclone count data from 1880 to 2019, carefully adjusting for the improvements in observation over the century and a half. The baseline model was a Poisson regression with three indices: the temperature of the main development region in August–October, the Niño 3.4 index in December–February, and the North Atlantic Oscillation in December–March. They then tested whether they could squeeze out more accuracy by adding quadratic and cross terms, or by using an elastic net on global sea surface temperature maps with a 2-degree resolution. All comparisons used five-fold cross-validation, with adjacent years kept in the same fold to prevent information leakage from autocorrelation.

The theoretical lower bound on mean absolute error for Poisson models was 2.51 cyclones. The baseline model achieved an error of 2.46 ± 0.10, and the elastic net 2.51 ± 0.14, meaning both reach the limit within uncertainty.

Validation on real data confirmed agreement with a Poisson distribution: a chi-square test gave p = 0.84, and residual autocorrelation was only 11% (z = 1.3). Adding nonlinear interactions did not improve the result, except for a weak effect from the product of the Niño 3.4 and North Atlantic Oscillation indices, which reduced the error by 0.03 cyclones but did not change the overall picture. In other words, no matter how much you complicate the model, you cannot predict hurricane counts more accurately within this approach.

What sets this work apart from neighboring studies is that instead of hunting for more powerful predictors, it proves the existence of an insurmountable statistical barrier. Neighboring work on El Niño and neural-network weather forecasts aims to improve accuracy, whereas here it is shown that for this class of models, that is impossible. This shifts the focus: further efforts should go not into complicating regressions, but into searching for hidden dependencies between years or new sources of information that break the assumption of independent Poisson counts.

The average number of cyclones per season in the Atlantic is about ten, but because of the Poisson nature of the process, the year-to-year spread is roughly three cyclones. This means that even a perfect model cannot predict the exact number better than with an error of about 2.5 cyclones on average.

A promising direction is the search for subtle correlations between adjacent seasons that standard tests do not catch. If such dependencies exist, they might be driven by slow processes in the ocean or atmosphere, such as multi-year temperature oscillations. Another path is to abandon the Poisson assumption in favor of distributions with heavier tails or models with memory. However, so far all checks point to agreement with a Poisson process, so the likelihood of a breakthrough is unclear.

The ceiling of randomness that models have hit does not spell the end of forecasting. It shows where the predictable ends and pure chance begins. For insurance companies and coastal communities, this is important knowledge: do not expect a more powerful computer or a deeper neural network to give a precise answer to the question 'how many hurricanes will there be next season?' Instead, it is wiser to prepare for the spread that is inherent in the very nature of tropical cyclones.

The average number of cyclones per season in the Atlantic is about ten, but due to the Poisson nature of the process, the year-to-year spread is roughly three cyclones. This means even a perfect model cannot predict the exact number better than with an error of about 2.5 cyclones on average.

In Kim Stanley Robinson's novel 'Forty Signs of Rain', scientists try to improve extreme weather forecasts while grappling with atmospheric chaos. Here the situation is similar: even with complete knowledge of climate indices, an irreducible randomness remains that no computation can overcome.

\mathrm{err}(z, \lambda) = \sum_{y \ge 0} |y - z| \, P(y|\lambda)
Average deviation of the predicted number of events from the observed one, averaged over the Poisson distribution.
b(z) = \min_{\lambda} \mathrm{err}(z, \lambda)
Minimum possible error for unknown but fixed intensity.
B = \frac{1}{n} \sum_{t=1}^{n} b(y_t)
Averaging over all observations gives the theoretical minimum mean absolute error for any Poisson model.
Scientists
Max PlanckClaude ShannonJohn von NeumannLudwig BoltzmannJosiah Willard GibbsAndrey Markov
Concepts
spring predictability barrier NOAA CPC El Niño index ERSST ENSO Oceanic Niño index Jin-Timmermann model TAO/TRITON moorings relative ONI tropical cyclone Southern Oscillation La Niña North Equatorial Countercurrent recharge oscillator model Pacific decadal oscillation natural variability seasonal cycle NOAA PMEL / TAO forecast skill Warm water volume