Code for the coupled Stommel–L84 model

Anticipating critical transitions in the Earth system is of great societal relevance, yet there may be intrinsic limitations to their predictability. For instance, from the theory of dynamical systems possessing multiple chaotic attractors, it is known that the asymptotic state depends sensitively on the initial condition in the proximity of a fractal basin boundary. Here, we approach the problem of final-state sensitivity of the Atlantic Meridional Overturning Circulation (AMOC) using a conceptual climate model, composed of a slow bistable ocean coupled to a fast chaotic atmosphere. First, we explore the occurrence of long chaotic transients in the monostable regime, which can mask a loss of stability near bifurcations. In the bistable regime, we explicitly construct the chaotic saddle using the edge tracking technique. Quantifying the final-state sensitivity through the maximum Lyapunov exponent and the lifetime of the saddle, we find that the system exhibits a fractal basin boundary with almost full phase space dimension, implying vanishing predictability of the second kind near the basin boundary. Our results demonstrate the usefulness of studying non-attracting chaotic sets in the context of predicting climatic tipping points, and provide guidance for the interpretation of higher-dimensional models such as general circulation models.