The conductor of an orchestra gradually loses control of the tempo. First, the second violins fall slightly behind, then the cellos speed up, and soon every musician is playing in their own rhythm, even though a minute ago they all followed the same melody. This is roughly what can happen in the El Niño climate system: two independent frequencies—the natural rhythm of the feedback and the imposed seasonal cycle—can coexist for years, then suddenly fall out of sync.

Inside this kitchen, the GZT model operates: a simple delay equation where the rate of change of the thermocline depth depends on its past value. When parameters enter the 2:7 resonance tongue (two feedback cycles over seven years), tori emerge in phase space: surfaces along which the trajectory moves like a thread on a bagel, winding with two different periods.

Theorists long predicted that as parameters change, two such tori can approach and disappear, leaving behind a region of complex dynamics: a Chenciner bubble. But drawing a diagram on paper is one thing; seeing it unfold in an equation with an infinite-dimensional phase space is another. The authors used the software package DDE-Biftool to continue periodic solutions and compute their Floquet multipliers, numbers that indicate how stable each cycle is.

Inside the bubble, for the 2:7 resonance, they found a structure that agrees well with the theoretical picture of Baesens and Mackay: curves of neutral saddles, homoclinic and heteroclinic transitions, and torus folding curves. But there was also a difference: the torus bifurcation turned out to be subcritical, meaning a saddle torus was born rather than a stable one. This implies that when crossing the critical line, the system does not smoothly transition to a new regime but instead snaps into a sharp switch, as if the orchestra did not gradually go out of tune but suddenly broke into cacophony.

Previously, it was unclear how tori break down in delay differential equations and whether the bifurcation structure matches the theoretical scheme for maps. Now, for the first time, the full picture has been reproduced numerically for a delay equation, and it has been shown that the torus bifurcation is subcritical. This was made possible by continuing periodic solutions and computing unstable manifolds of saddle orbits in two-dimensional projections associated with Poincaré sections. This approach allowed tracking all key bifurcations inside the 2:7 resonance tongue.

The results may help in interpreting climate time series and in searching for precursors of abrupt transitions in systems with multiple frequencies, although there is no direct practical application yet. The conductor who loses control is not just a pretty metaphor but also a warning: if we learn to recognize the moment when tori are about to fold, we might just catch the first false notes in the climate orchestra.