AI Discovers Hidden Equations Governing North Atlantic Climate Swings

For decades, scientists have struggled to untangle the slow, mysterious rhythms of the North Atlantic, where ocean temperatures, atmospheric pressure patterns and rainfall rise and fall in ways that shape hurricanes, droughts and European winters. Now a team of researchers from the University of Oxford, the UK Centre for Ecology and Hydrology and the Met Office Hadley Centre has turned to artificial intelligence not as a black-box forecaster, but as an equation hunter. By applying a machine-learning technique called equation discovery to more than seventy years of climate data, they have distilled the sprawling complexity of the North Atlantic climate system into just three coupled ordinary differential equations. The findings, published in Earth System Dynamics, suggest that rainfall plays a far more active role in driving decadal climate variability than previously appreciated, and that these simple equations can skillfully predict ocean temperature, pressure and precipitation patterns years to a decade ahead.

The two dominant players in North Atlantic climate are the Atlantic Multidecadal Variability, or AMV, a basin-wide fluctuation in sea surface temperatures, and the North Atlantic Oscillation, or NAO, the pressure seesaw between the Azores high and the Icelandic low. Both are major sources of predictability on seasonal to decadal timescales, influencing Atlantic hurricane frequency, shifts in the Intertropical Convergence Zone, drought in the Sahel, and summer climate across Europe and North America. Yet the physical processes and feedback mechanisms linking the oceanic AMV to the atmospheric NAO have remained fiercely debated, with proposed drivers ranging from internal ocean circulation to anthropogenic aerosols, greenhouse gases and stochastic atmospheric forcing. What has been missing, the researchers argue, is a compact mathematical description that captures how these variables actually interact over decades.

To find one, the team employed the Sparse Identification of Nonlinear Dynamics algorithm, known as SINDy, a sparse regression method that learns governing equations directly from time-series data. The researchers trained the algorithm on monthly indices of the AMV, the NAO and North Atlantic precipitation anomalies derived from the ERA5 reanalysis dataset spanning January 1950 to December 2022. A twelve-month running average was applied to filter out high-frequency noise, and the candidate library of functions included constant, linear and quadratic polynomial terms, yielding thirty possible terms across three equations. A LASSO optimiser with sparsity thresholding forced the regression to retain only the most important terms, while a bootstrapping ensemble of one thousand resampled datasets ensured the identified coefficients were robust. Validation and testing used independent segments of the ERA20C reanalysis from 1900 to 1949, data the models had never seen during training.

The procedure generated 344 unique candidate models of varying complexity, and the team faced a subtle selection problem. Simpler models with thirteen or fewer terms either collapsed onto stable fixed points, converged on rigid limit cycles, or blew up entirely, none of which resembles the restless variability of the real North Atlantic. Models with fourteen to twenty terms were stable but produced unnaturally sparse power spectra with delta-like peaks, unlike the broad-band spectra seen in observations. Only models with twenty-one or more terms exhibited chaotic attractors and physically realistic variability. The researchers ultimately selected a twenty-six-term model sitting on the Pareto front, the optimal trade-off between accuracy and complexity, whose fifty-thousand-year simulated trajectories reproduce the irregular, Lorenz-like character of the observed climate system.

Remarkably, the decadal timescales embedded in these equations were never explicitly taught to the model. Because SINDy fits instantaneous tendencies, the emergence of a pronounced quasi-periodic twenty-year oscillation across all three variables is an entirely emergent property. The model’s power spectra show strong peaks near twenty years and thirty-five years, matching significant signals in the observed NAO and precipitation records, along with weaker peaks near 13.5 and 15.5 years. The team interprets the twenty-year mode as a damped oscillatory ocean mode, set by internal ocean dynamics but excited by decadal NAO forcing, consistent with earlier inverse modelling studies that identified oscillatory eigenmodes of North Atlantic sea surface temperature at periods of roughly twenty and thirty-seven years.

The most striking discovery lies in the structure of the equations themselves. Cross-terms multiplying precipitation with either the AMV or the NAO, written as AP and NP products, appear among the largest and most robust coefficients across the entire ensemble of 344 candidate models, persisting even in the sparsest formulations. Because precipitation serves as a proxy for both latent heat release in the atmosphere and freshwater fluxes into the ocean, these terms point to diabatic processes acting as a dynamical bridge between ocean and atmosphere. In the AMV equation, the negative NP term suggests that heavy rainfall over the Labrador Sea region, combined with a positive NAO, freshens the surface ocean, increases stratification, weakens deep convection and thereby dampens the overturning circulation and the AMV itself.

Feedbacks onto the atmosphere appear equally important. In the NAO equation, the positive NP term indicates that enhanced convective precipitation near the southern node of the NAO reinforces a negative NAO state through diabatic heating that lowers surface pressure, echoing recent causal studies of oceanic feedbacks onto the winter NAO. Conversely, the model suggests that the positive feedback of North Atlantic sea surface temperature onto the NAO, in which the familiar SST tripole strengthens the meridional temperature gradient and the zonal winds, is strongest precisely when the precipitation pattern favours a positive NAO state. In other words, the response of the atmosphere to the ocean appears to be conditioned on the sign of rainfall anomalies, a state-dependent coupling that conventional linear analyses would miss entirely.

The predictive credentials of the twenty-six-term model are substantial. Initialised every six months from observed data and integrated forward with a fourth-order Runge-Kutta scheme, the model achieves anomaly correlation coefficients of 0.86 for the AMV and 0.54 for the NAO at one-year lead times on training data, and 0.93 and 0.52 respectively on the unseen early-twentieth-century test data. Most impressively, the model shows significant skill in predicting North Atlantic precipitation anomalies out to a decade ahead, with correlations reaching 0.7 at one-year leads on test data. The authors attribute this precipitation skill to the fact that the decadal precipitation index varies in phase with the damped twenty-year ocean oscillation, so that slow, predictable ocean responses to NAO forcing propagate directly into rainfall patterns over the ocean and adjacent continents.

The researchers are careful to note the limits of their framework. The equations encode dynamical couplings but do not by themselves distinguish cause from effect, and the mechanistic pathways they propose are supported by consistency with existing theory rather than uniquely proven. The characteristic fifty-to-seventy-year signal of the AMV is not well captured, likely because the training record is too short, and external influences such as volcanic eruptions and lunar cycles were not considered. The authors also acknowledge that stochastic processes present in the real system are only implicitly represented through the nonlinear terms, and they propose comparing their models against Linear Inverse Models and coupled climate simulations as future work.

Even so, the study represents a compelling demonstration that machine learning can do more than forecast: it can illuminate physics. As the journal’s editorial statement observes, the paper exemplifies how equation discovery transforms AI from a black box into a hypothesis-generating scientific instrument. If the identified feedbacks hold up under scrutiny, the implications extend beyond academic understanding. Correctly initialised versions of these low-order models could sharpen decadal forecasts of rainfall over the North Atlantic and neighbouring regions such as Europe and North America, offering societies a rare glimpse of climate conditions years in advance. In an era when decadal prediction remains one of climate science’s hardest challenges, three short equations learned by an algorithm may prove to be an unexpectedly powerful key.

Subject of Research: Data-driven discovery of coupled ocean-atmosphere equations governing decadal North Atlantic climate variability

Article Title: New insights into decadal climate variability in the North Atlantic revealed by data-driven dynamical models

Article References: Nicoll, A. J., Christensen, H. M., Huntingford, C., & Smith, D. (2026). New insights into decadal climate variability in the North Atlantic revealed by data-driven dynamical models. Earth System Dynamics, 17(4), 1061-1079. https://doi.org/10.5194/esd-17-1061-2026

Image Credits: AI Generated

DOI: 10.5194/esd-17-1061-2026

Keywords: North Atlantic, Atlantic Multidecadal Variability, North Atlantic Oscillation, equation discovery, SINDy, machine learning, decadal prediction, precipitation feedbacks, climate variability, ocean-atmosphere coupling, ERA5 reanalysis, latent heat