Climate tipping elements are large-scale subsystems of the Earth that may transgress critical thresholds (tipping points) under ongoing global warming, with substantial impacts on the biosphere and human societies. Frequently studied examples of such tipping elements include the Greenland Ice Sheet, the Atlantic Meridional Overturning Circulation (AMOC), permafrost, monsoon systems, and the Amazon rainforest. While recent scientific efforts have improved our knowledge about individual tipping elements, the interactions between them are less well understood. Also, the potential of individual tipping events to induce additional tipping elsewhere or stabilize other tipping elements is largely unknown. Here, we map out the current state of the literature on the interactions between climate tipping elements and review the influences between them. To do so, we gathered evidence from model simulations, observations, and conceptual understanding, as well as examples of paleoclimate reconstructions where multi-component or spatially propagating transitions were potentially at play. While uncertainties are large, we find indications that many of the interactions between tipping elements are destabilizing. Therefore, we conclude that tipping elements should not only be studied in isolation, but also more emphasis has to be put on potential interactions. This means that tipping cascades cannot be ruled out on centennial to millennial timescales at global warming levels between 1.5 and 2.0 ^∘C or on shorter timescales if global warming surpassed 2.0 ^∘C. At these higher levels of global warming, tipping cascades may then include fast tipping elements such as the AMOC or the Amazon rainforest. To address crucial knowledge gaps in tipping element interactions, we propose four strategies combining observation-based approaches, Earth system modeling expertise, computational advances, and expert knowledge.

Significance: Ongoing greenhouse gas emissions put elements of the Earth system at risk for crossing critical thresholds (tipping points), leading to abrupt irreversible climate change. Measures for reducing emissions should keep Earth in the safe operating space away from tipping points. Here we show that increasing rates of change of ice melt can induce a collapse of the Atlantic Meridional Overturning Circulation in a global ocean model, while no critical threshold in ice melt is crossed and slower increases to the same level of ice melt do not induce tipping. Moreover, the chaotic dynamics of the climate make such a collapse hard to predict. This shows that the safe operating space of the Earth system might be smaller than previously thought.

Large volcanic eruptions occurring in the last glacial period can be detected by their accompanying sulfuric acid deposition in continuous ice cores. Here we employ continuous sulfate and sulfur records from three Greenland and three Antarctic ice cores to estimate the emission strength, the frequency and the climatic forcing of large volcanic eruptions that occurred during the second half of the last glacial period and the early Holocene, 60–9 kyr before 2000 CE (b2k). Over most of the investigated interval the ice cores are synchronized, making it possible to distinguish large eruptions with a global sulfate distribution from eruptions detectable in one hemisphere only. Due to limited data resolution and large variability in the sulfate background signal, particularly in the Greenland glacial climate, we only list Greenland sulfate depositions larger than 20 kg km⁻² and Antarctic sulfate depositions larger than 10 kg km⁻². With those restrictions, we identify 1113 volcanic eruptions in Greenland and 737 eruptions in Antarctica within the 51 kyr period – for which the sulfate deposition of 85 eruptions is found at both poles (bipolar eruptions). Based on the ratio of Greenland and Antarctic sulfate deposition, we estimate the latitudinal band of the bipolar eruptions and assess their approximate climatic forcing based on established methods. A total of 25 of the identified bipolar eruptions are larger than any volcanic eruption occurring in the last 2500 years, and 69 eruptions are estimated to have larger sulfur emission strengths than the Tambora, Indonesia, eruption (1815 CE). Throughout the investigated period, the frequency of volcanic eruptions is rather constant and comparable to that of recent times. During the deglacial period (16–9 ka b2k), however, there is a notable increase in the frequency of volcanic events recorded in Greenland and an obvious increase in the fraction of very large eruptions. For Antarctica, the deglacial period cannot be distinguished from other periods. This confirms the suggestion that the isostatic unloading of the Northern Hemisphere (NH) ice sheets may be related to the enhanced NH volcanic activity. Our ice-core-based volcanic sulfate records provide the atmospheric sulfate burden and estimates of climate forcing for further research on climate impact and understanding the mechanism of the Earth system.

Tipping points (TP) in climate subsystems are usually thought to occur at a well-defined, critical forcing parameter threshold, via destabilization of the system state by a single, dominant positive feedback. However, coupling to other subsystems, additional feedbacks, and spatial heterogeneity may promote further small-amplitude, abrupt reorganizations of geophysical flows at forcing levels lower than the critical threshold. Using a primitive-equation ocean model, we simulate a collapse of the Atlantic Meridional Overturning Circulation (AMOC) due to increasing glacial melt. Considerably before the collapse, various abrupt, qualitative changes in AMOC variability occur. These intermediate tipping points (ITP) are transitions between multiple stable circulation states. Using 2.75 million years of model simulations, we uncover a very rugged stability landscape featuring parameter regions of up to nine coexisting stable states. The path to an AMOC collapse via a sequence of ITPs depends on the rate of change of the meltwater input. This challenges our ability to predict and define safe limits for TPs.

We propose a conceptual model comprising a cascade of tipping points as a mechanism for past abrupt climate changes. In the model, changes in a control parameter, which could for instance be related to changes in the atmospheric circulation, induce sequential tipping of sea ice cover and the ocean's meridional overturning circulation. The ocean component, represented by the well-known Stommel box model, is shown to display so-called rate-induced tipping. Here, an abrupt resurgence of the overturning circulation is induced before a bifurcation point is reached due to the fast rate of change of the sea ice. Because of the multi-scale nature of the climate system, this type of tipping cascade may also be a risk concerning future global warming. The relatively short timescales involved make it challenging to detect these tipping points from observations. However, with our conceptual model we find that there can be a significant delay in the tipping because the system is attracted by the stable manifold of a saddle during the rate-induced transition before escaping towards the undesired state. This opens up the possibility for an early warning of the impending abrupt transition via detection of the changing linear stability in the vicinity of the saddle. To do so, we propose estimating the Jacobian from the noisy time series. This is shown to be a useful generic precursor to detect rate-induced tipping.

The strongest mode of centennial to millennial climate variability in the paleoclimatic record is represented by Dansgaard–Oeschger (DO) cycles. Despite decades of research, their dynamics and physical mechanisms remain poorly understood. Valuable insights can be obtained by studying high-resolution Greenland ice core proxies, such as the NGRIP δ¹⁸O record. However, conventional statistical analysis is complicated by the high noise level, the cause of which is partly due to glaciological effects unrelated to climate and which is furthermore changing over time. We remove the high-frequency noise and extract the most robust features of the DO cycles, such as rapid warming and interstadial cooling rates, by fitting a consistent piecewise linear model to Greenland ice core records. With statistical hypothesis tests we aim to obtain an empirical understanding of what controls the amplitudes and durations of the DO cycles. To this end, we investigate distributions and correlations between different features, as well as modulations in time by external climate factors, such as CO₂ and insolation. Our analysis suggests different mechanisms underlying warming and cooling transitions due to contrasting distributions and external influences of the stadial and interstadial durations, as well as the fact that the interstadial durations can be predicted to some degree by linear cooling rates already shortly after interstadial onset.