Why El Niño Loves December ✛ ✎

Original: "Resonance structure of a periodically forced delay differential equation model for the El Niño--Southern Oscillation"
arXiv:2607.00386 · 2026-07-01 · CC BY-NC-ND 4.0 · independent analysis · 1 min
The ocean remembers its past, and the seasons nudge it into a certain rhythm.

Every winter, the water off the coast of Peru turns warmer or cooler than usual. The warm state is called El Niño, the cold one La Niña. Both phenomena influence weather across the globe: they bring droughts in some regions and floods in others. And both most often reach their greatest strength in December.

The ocean behaves like a swing with a delay: if the water warms up today, the effects won't show up right away but several months later. That's why scientists describe this phenomenon with an equation that has memory. In such a model, the rate of temperature change depends not only on what's happening now but also on what happened before. The seasonal cycle, that is, the changing of the seasons, acts as a rhythmic external push.

The season sets the rhythm but doesn't control the motion. A swing can oscillate at its own pace, or it can sync with an external rhythm if the frequencies match.

The mathematicians added seasonal forcing to the equation and mapped how the system behaves under different conditions. It turned out that the regions of rhythm synchronization are very narrow. That's why real El Niño doesn't repeat strictly every year, although it does tend to peak in December. The scientists also found that the model doesn't capture all the chaos of the real climate, so it needs to be made more complex for accurate forecasts.

The rotation number on an invariant torus can be rational (then the dynamics are periodic) or irrational (then the trajectory densely fills the torus). In this work, computing the rotation number required about 768,000 calculations, taking roughly 8 hours on an Apple M1 computer.

u'(t) = u(t) - u^3(t) - \alpha u(t-\tau)
The rate of change of the sea surface temperature anomaly depends on the current value, cubic saturation, and delayed negative feedback.
u'(t) = u(t) - u^3(t) - \alpha u(t-\tau) + c \cos(\omega t)
A periodic forcing with amplitude c and angular frequency ω, modeling the seasonal cycle, is added to the autonomous model.
\Phi(p:q) = p:(q+2p)
When the delay is increased by 1, a resonance tongue with ratio p:q maps to a tongue with ratio p:(q+2p), corresponding to a double Dehn twist on the torus.
Scientists
Ludwig BoltzmannJohn von NeumannJosiah Willard GibbsHenri Poincaré
Concepts
El Niño-Southern Oscillation recharge oscillator model ENSO Jin-Timmermann model Multivariate ENSO Index Southern Oscillation index spring predictability barrier TAO/TRITON moorings seasonal forcing La Niña Bjerknes feedback delayed oscillator theory tropospheric biennial oscillation ERSST