Why Climate Is Lopsided: A New Lens for Data ✛ ✎

Original: "Beyond Gaussian Assumptions: A Nonlinear Generalization of Linear Inverse Modeling"
· Justin Lien, Hiroyasu Ando
arXiv:2503.24234 · 2025-03-31 · arXiv non-exclusive · independent analysis · 1 min
Standard models see climate as symmetric, but a new lens catches the imbalance between warm El Niño and cold La Niña.

Flat glass doesn't distort the image, but it also doesn't show a tilt. For a long time, climatologists looked at El Niño through such glass and saw symmetric oscillations: up and down, warm and cold. In reality, the warm phase of El Niño is often stronger than the cold La Niña, but the glass didn't notice.

Scientists made a lens with a curve: they added quadratic terms to the equation. It's like giving the glass curvature so it focuses not just the average but also the skew. The shape of the lens is chosen based on how data points relate to each other at different times. If before they only looked at pairwise connections, now they account for connections among three and four time points at once.

On Pacific Ocean data spanning 42 years, the new lens revealed a skew that was previously invisible. The gap between model and reality shrank by nearly half.

Now climatologists can see asymmetry elsewhere: in the atmosphere, the ocean, even in chaotic systems. The lens has become sharper, and the world in it is closer to the real one.

The minimum dimension for chaos in an autonomous differentiable system is three, so the authors first test the method on a two-dimensional system (where chaos is absent) and then on the three-dimensional Lorenz system.

\frac{d}{dt}x_i = \sum_{j,k} B_{ijk} x_j x_k + \sum_j A_{ij} x_j + C_i + \sum_j (\sqrt{2Q})_{ij} \cdot \text{noise}_j
The rate of change of each variable depends on quadratic interactions, linear damping, a constant shift, and random forcing.
0 = \text{Sym}\left( B \times_2 M(0) + A \times_1 K(0) + C \times_0 E^T + Q \right)
In the stationary state, deterministic damping and stochastic forcing balance each other, which allows the noise matrix to be determined.
Q_{C-nLIM} = \arg\min_Q \left( \|R\|_F + \|K(0) - K_{obs}(0)\|_F + \|M(0) - M_{obs}(0)\|_F \right)
The noise matrix is chosen so that synthetic correlations and the balance residual are as close as possible to observed ones.
Concepts
linear inverse model nonlinear laplacian spectral analysis normalizing flows probabilistic model Stochastic modeling gravitational lensing modeling nonlinear dynamics gaussian process Contrastive Learning linear regression quantum-like models Large Language Models Gaussian distribution weighted mean Weakly Supervised Learning matched filter simulated data nonlinearity physics informed neural network lattice simulation