Estimation of AMOC transition probabilities using a machine learning based rare-event algorithm

TL;DR

This paper introduces a novel approach combining a rare-event algorithm with machine learning to accurately estimate the probability of Atlantic Meridional Overturning Circulation (AMOC) collapse, improving understanding of climate tipping points.

Contribution

It develops a method integrating Trajectory-Adaptive Multilevel Splitting with Reservoir Computing to estimate transition probabilities and paths in climate models, overcoming limitations of traditional score functions.

Findings

Accurately estimates AMOC transition probabilities and paths.

Demonstrates effectiveness for both fast and slow transition types.

Provides a way to approximate the committor function analytically.

Abstract

The Atlantic Meridional Overturning Circulation (AMOC) is an important component of the global climate, known to be a tipping element, as it could collapse under global warming. The main objective of this study is to compute the probability that the AMOC collapses within a specified time window, using a rare-event algorithm called Trajectory-Adaptive Multilevel Splitting (TAMS). However, the efficiency and accuracy of TAMS depend on the choice of the score function. Although the definition of the optimal score function, called ``committor function" is known, it is impossible in general to compute it a priori. Here, we combine TAMS with a Next-Generation Reservoir Computing technique that estimates the committor function from the data generated by the rare-event algorithm. We test this technique in a stochastic box model of the AMOC for which two types of transition exist, the so-called F(ast)-transitions and S(low)-transitions. Results for the F-transtions compare favorably with those in the literature where a physically-informed score function was used. We show that coupling a rare-event algorithm with machine learning allows for a correct estimation of transition probabilities, transition times, and even transition paths for a wide range of model parameters. We then extend these results to the more difficult problem of S-transitions in the same model. In both cases of F-transitions and S-transitions, we also show how the Next-Generation Reservoir Computing technique can be interpreted to retrieve an analytical estimate of the committor function.