main concepts scientists sections
SimplePopularAdvanced

Climate as an Orchestra: How to Hear the Response to Sharp Sounds ✛

Original: "Kolmogorov Modes and Linear Response of Jump-Diffusion Models"
· Mickaël D. Chekroun, Niccolò Zagli, Valerio Lucarini
arXiv:2411.14769 · 2024-11-22 · CC BY 4.0 · 1 min
Chaotic Dynamics Statistical Mech Math Physics Atmospheric and Oceanic Physics
Scientists have found a way to predict how the climate will respond to sudden events, such as a sharp gust of wind.

Climate is like an orchestra. It has its own natural rhythms, like the parts of violins and winds. Sometimes sharp sounds break into the music: a sudden gust of wind or a heavy downpour. Previously, scientists could only predict how the system would respond to gradual changes, while sharp sounds remained a mystery.

Now, formulas have emerged that help understand the response to any sharp jolt. The idea is simple: decompose the response into the same rhythms that are playing in the calm state. Then the reaction to a disturbance becomes a familiar melody, not a cacophony.

In the El Niño model, which describes temperature fluctuations in the Pacific Ocean, the theory accurately predicted the response to a change in the ocean-atmosphere coupling. It's as if the orchestra kept playing without a wrong note, even when the conductor slightly changed the tempo.

Testing on another climate model showed that the calculations match reality, although this requires more computation. But now we can assess how the climate will respond to extreme events and better prepare for them.

Whole graph →

Ulam's method, used to estimate modes, was originally proposed by Stanislaw Ulam in the 1960s to study ergodic properties of dynamical systems, and now it helps analyze climate models with jumps.

\partial_t p(x,t) = -\nabla \cdot (F(x) p(x,t)) + \frac{1}{2} \sum_{i,j=1}^d \partial_{ij} (a_{ij}(x) p(x,t)) + \int_{\mathbb{R}^d \setminus \{0\}} [p(x-y,t) - p(x,t) + y \cdot \nabla p(x,t) \mathbf{1}_{\{\|y\|<1\}}] \nu(dy)
The change in probability density is the sum of drift, diffusion, and a nonlocal jump contribution described by an integral over the Lévy measure ν.
\delta^{(1)}[\Psi](t) = \epsilon \int_{-\infty}^{t} G_{\Psi,G}(t-s) g(s) ds + \epsilon \int_{-\infty}^{t} ds \int_{\mathbb{R}^d \setminus \{0\}} G_{\Psi,\nu}(t-s,y) \delta\nu_s(dy)
The response is the sum of two contributions: a local one from the drift change and a nonlocal one from the change in the jump law, each expressed through its own Green's function.
Concepts
fluctuation-dissipation theorem coupled general circulation model linear inverse model Kolmogorov–Arnold–Moser theorem Lotka-Volterra equation linear regression game theory Jin-Timmermann model