Points of No Return in El Niño: How Delay Breeds Climate Chaos ✛

Original: "Effective Reduced Models from Delay Differential Equations: Bifurcations, Tipping Solution Paths, and ENSO variability"
· Mickaël D. Chekroun, Honghu Liu
arXiv:2402.03318 · 2024-01-02 · CC BY 4.0 · 2 min
A simplified delayed model shows that random kicks and a slow drift of the delay turn regular oscillations into the complex patterns of El Niño.

The Pacific Ocean remembers. Water warmed off the coast of Peru drifts slowly westward, and six months later its heat returns to South America. Scientists call this six-month journey the delay. It turns a simple climate model into a maze where the route depends not only on your current position but also on where you were half a year ago.

The Suarez–Schopf model describes ocean surface temperature with a single delay equation. Despite its simplicity, it reproduces the irregularity of El Niño.

When the delay changes, the maze reshapes itself. At a delay of about 1.74 years, the system loses its stable equilibrium and begins to oscillate. But the oscillations do not become chaotic right away. First, an unstable periodic orbit is born: a closed trajectory that the system repels. If a random kick pushes the state aside, it may wander for a long time before returning.

The authors reduced the complex equation to a two-dimensional system and found three key points: the onset of oscillations, the appearance of a loop connecting an unstable state to itself, and the disappearance of stable cycles. These points, like invisible landmarks, define the boundaries of possible routes.

At a delay of 1.59 years, a homoclinic loop emerges: a trajectory that leaves a saddle and returns to it. Near this loop, the system spends a lot of time, giving rise to decadal variability.

Then noise was added to the model. Random disturbances, like gusts of wind, constantly knock the system off course. Meanwhile, a slow drift of the delay, mimicking changes in wave travel time, makes the maze slowly change shape. As a result, episodes resembling real El Niño and La Niña appear: warm phases give way to cold ones, and the rhythm is irregular.

Spectral analysis revealed two peaks: around 5.5 years for interannual variability and around 21 years for decadal variability. The decadal peak is linked to unstable periodic orbits near the homoclinic loop. Wandering through the maze, the system occasionally approaches these orbits and lingers near them for several cycles before being repelled again.

The maze of invisible orbits explains why El Niño is so hard to predict. Even knowing all the bifurcations, we cannot say exactly where the system will turn next. But understanding which routes are possible helps us recognize when we are approaching a point of no return. This could form the basis for early warnings of strong events, and for more accurate climate models that account not only for the current state but also for the ocean's memory.

The Suárez–Schopf model was proposed in 1988 and has since served as a testing ground for studying El Niño. Despite its simplicity, it reproduces key features of the phenomenon, including irregularity and decadal modulations.

The idea that small random perturbations can switch a system between qualitatively different regimes echoes the concept of 'bifurcation points' in science fiction, for instance, in Isaac Asimov's novel 'The End of Eternity', where small changes in the past lead to unpredictable consequences in the future.

\frac{dT}{dt} = T(t) - \alpha T(t-\tau) - T^3(t)
The rate of change of the anomaly equals linear amplification, delayed negative feedback, and nonlinear saturation.
d\theta(t) = [a\theta(t) - \alpha\theta(t-\tau(t)) - b\theta^2(t) - \theta^3(t)]dt + \frac{\sigma}{1+|\theta^2(t)|}dW_t
The deterministic part mirrors the Suárez–Schopf model, while the noise term is stronger for small deviations.
Concepts
Jin-Timmermann model ENSO El Niño-Southern Oscillation spring predictability barrier TAO/TRITON moorings Multivariate ENSO Index El Niño index recharge oscillator model equatorial wave propagation delay differential equation mixed mode oscillations Pacific decadal oscillation natural variability hopf bifurcation NOAA PMEL / TAO delayed oscillator theory seasonal cycle saddle node bifurcation decadal variability differential equations