Entropy of Climate Whiplash ✛

Original: "Statistics of stochastic entropy for recorded transitions between ENSO states"
· Sílvio M. Duarte Queirós
arXiv:2507.10516 · 2025-07-14 · CC BY 4.0 · 3 min
Transitions between El Niño and La Niña turn out to be statistically ordinary events, not extreme entropy spikes.

Context

El Niño and La Niña shift rainfall and winds across half the planet, yet their alternation is still described mostly in phenomenological terms. Nearby in the archive are studies on predicting El Niño a year and a half ahead, on randomness in its evolution, and on reviving old weather records. All of them treat climate either as a deterministic system or as time-series statistics. This paper approaches from a different angle: it treats the Southern Oscillation Index as the trajectory of a stochastic system driven by a slow protocol, and asks how unusual the phase transitions themselves are in terms of entropy production.

Methods

The author takes daily SOI values from 1991 to 2023, totaling 11,744 records. First, a multitaper spectral estimate extracts a background protocol with seven significant periods, including the annual cycle and the Madden–Julian oscillation. The residual is modeled with a stochastic differential equation with multiplicative noise, whose coefficients are estimated from empirical Kramers–Moyal moments. To compute transition probabilities, the author uses a transformation to additive noise and an optimal-path approximation via the Euler–Lagrange equation. Entropy change is computed in two ways: along the full realized trajectory and using only the initial and final states.

Results

Of the 11 recorded transitions, only three deviate by more than two standard deviations. The transition from the strong La Niña of 1999–2000 to the moderate El Niño of 2002–2003 reaches four sigma when the full trajectory is considered and three sigma when only the endpoints are used. This is the only candidate for an extreme event, and it barely crosses the threshold. All transitions satisfy the integral fluctuation theorem, which is the probabilistic version of the second law of thermodynamics. The mean rate of entropy change does not correlate with phase strength, which separates the index's informational entropy from the ocean's thermal entropy.

Implications

The work sits at the periphery: its nearest neighbors are busy predicting El Niño, studying randomness in evolution, or chasing quantum tails. Here, for the first time, the formalism of fluctuation theorems, usually applied to small systems, is brought to a climate index. The difference from neighbors is that ENSO transitions are treated not as predictable events but as statistical realizations of a stochastic process. This shifts the question from 'when will it happen' to 'how likely is such a path'.

Future development

The method could be extended to multivariate ENSO indicators: sea surface temperature, outgoing longwave radiation, precipitation. Then entropy analysis would no longer depend on a single index and would reveal which system components truly carry information about the phase transition. A move to self-organizing maps of climate patterns, with entropy defined over sequences of states, is also possible.

Impact

The result matters for climate diagnostics and for the statistical physics of nonequilibrium systems. It provides a benchmark for what ordinary transitions look like in a large-scale climate system.

Next steps

Test fluctuation relations of the Speck–Seifert and Hatano–Sasa type for individual ENSO phases. Extend the analysis to multivariate data and compare entropy measures across different indicators.

Key open problems

In climatology, the question remains open how adequately the SOI reflects the state of the entire ocean–atmosphere system. This work shows that the index's informational entropy is weakly linked to the phenomenon's thermal entropy, suggesting the indicator may miss important details of transitions.

The knowledge machine's view

This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.

What the work does

The author models daily SOI values as a stochastic process with multiplicative noise, extracting the slow protocol via multitaper spectral estimation and estimating the coefficients of the stochastic differential equation from Kramers–Moyal moments. The entropy change along ENSO transition trajectories is computed in two ways, and the integral fluctuation theorem is tested. Of the 11 transitions, only the 1999–2002 transition approaches the extremality threshold.

Where its strength is

Bringing the formalism of fluctuation theorems to a climate index opens a fresh perspective: ENSO transitions are treated as statistical realizations rather than predictable events. Separating information entropy from thermal entropy allows discussing the unusualness of paths without invoking ocean physics.

Where it could go next

Worth trying

  • Test whether the distribution of entropy production changes when individual spectral components of the protocol (e.g., the annual cycle) are excluded.
  • Apply the method to other climate indices (e.g., PDO or AMO) to see how universal the found extremality threshold is.
  • Compare entropy measures for transitions of different durations to reveal a possible dependence on time scale.

El Niño literally translates from Spanish as 'the boy,' and La Niña as 'the girl.' The names arose because the warm current off Peru often appeared around Christmas, associating it with the Christ child.

\Delta S(\vec{s}) = -\ln\frac{p(s_{t_f}; \{s_{t_i},\ldots,s_{t_{f-1}}\})}{p(s_{t_i}; \{s_{t_{i+1}},\ldots,s_{t_f}\})}
Logarithm of the ratio of probabilities of forward and backward SOI trajectories
p_t^*(s) = \frac{1}{Z_t^*}\exp\left[\frac{a\beta+2b\alpha^2}{\alpha^2\sqrt{4\alpha^2\gamma^2-\beta^2}}\arctan\left(\frac{\beta+2\alpha^2(s-\Phi(t))}{\sqrt{4\alpha^2\gamma^2-\beta^2}}\right)\right]\left(\gamma^2+\beta(s-\Phi(t))+\alpha^2(s-\Phi(t))^2\right)^{-(1+\frac{a}{2\alpha^2})}
Probability distribution of SOI values given the background protocol
ds = D_1(s,t)\,dt + \sqrt{2D_2(s,t)}\,dW_t
SOI dynamics with drift and multiplicative noise

Key numbers

  • transitions analyzed: 11
  • daily SOI records: 11744
  • data period: 1991–2023
  • maximum entropy deviation: 4 sigma
  • significant periods in protocol: 7
Scientists
Ludwig BoltzmannRudolf ClausiusJacob BekensteinStephen HawkingJames Clerk MaxwellEdward Witten
Concepts
ENSO El Niño index El Niño-Southern Oscillation Multivariate ENSO Index ERSST spring predictability barrier Southern Oscillation La Niña Jin-Timmermann model recharge oscillator model NOAA CPC interannual variability TAO/TRITON moorings natural variability Warm water volume relative ONI stochastic dynamical system seasonal cycle mixed mode oscillations Westerly wind burst