The team focused on causal inference methods, a class of statistical techniques that go beyond simple correlation to reconstruct the directed web of influences among variables in a dynamic system. Correlation alone cannot distinguish whether sea ice loss drives ocean circulation changes, whether circulation changes drive sea ice loss, or whether both are merely responding to a common background factor such as rising global temperatures. Causal methods attempt to resolve this ambiguity by exploiting the temporal structure of the data, testing whether knowledge of one variable’s past improves predictions of another’s present in ways that cannot be explained by other candidate causes. Their promise for tipping point research is considerable, because they could, in principle, extract interaction evidence directly from observational time series without waiting for fully coupled Earth system models to represent every relevant process dynamically.

Yet applying these methods to climate tipping elements is fraught with difficulty, and the researchers identified four central challenges. First, the modern observational record is short relative to the timescales of tipping elements, so the number of available samples is often below one thousand. Second, interactions between elements may be weak or highly delayed, since effects must propagate through atmospheric or oceanic transport. Third, researchers may need to analyze large and dense networks of variables, for example when regional tipping patterns are of interest. Fourth, global warming acts as a confounder, influencing all tipping elements simultaneously and potentially introducing nonlinear, noisy trends that can masquerade as causal links between them. Any method that fails under these conditions could produce misleading conclusions about the stability of the climate system.

To evaluate how well existing techniques cope, the authors generated synthetic data from networks of cubic stochastic differential equations, a mathematical form that reproduces the hallmark behavior of tipping elements: hysteresis between two stable states and an abrupt transition once a forcing threshold is crossed. They then fed these time series into three widely used multivariate causal inference methods and scored each method’s ability to recover the true underlying network of interactions. The benchmark metric was the Matthews Correlation Coefficient, chosen because it rewards both correct detections and correct rejections symmetrically, avoiding biases that would favor methods in sparse or dense networks. Each experimental configuration was repeated one hundred times to quantify the variability of the results.

The three methods tested represent distinct philosophies of causality. The Liang–Kleeman Information Flow (LKIF) takes an information-theoretic approach, measuring how the entropy of one variable would change if another were removed, and fits a low-complexity linear stochastic model to the data. The Peter–Clark Momentary Conditional Independence algorithm (PCMCI) iteratively prunes a fully connected candidate network using conditional independence tests, and offers unusual flexibility: users can impose known connections, prohibit implausible ones, and mask out sections of the data, for instance to restrict analysis to particular seasons. Granger Causality for State Space Models (GCSS), a method more familiar in neuroscience than in climate science, fits a latent state space model in which hidden variables can implicitly encode time-shifted information, giving it a natural capacity for handling delayed interactions at the cost of higher data demands.

The benchmark results revealed clear niches for each method. With limited data, in the range of a few hundred samples, LKIF outperformed the alternatives, converging early but plateauing at imperfect accuracy because its strict linear model assumptions cannot fully capture the nonlinear dynamics. GCSS, by contrast, achieved nearly perfect detection when given large sample counts or strong interactions, but its performance degraded sharply in larger networks. Time delays proved to be a critical discriminator: LKIF’s accuracy collapsed even at delays of a single sampling step, because its underlying model cannot represent delayed feedback loops, while GCSS handled delays of up to five samples robustly and PCMCI declined only gradually. The authors’ first recommendation follows directly from this finding: the sampling rate of observations should match the expected timescale of the interaction delays, and should not be orders of magnitude finer than the internal dynamics of the systems involved.

The confounder experiments carried perhaps the most consequential message for applied work. When a warming-like forcing was applied to all variables but excluded from the causal analysis, false positives rose significantly for LKIF and true positives dropped for GCSS. Including the forcing variable as an explicit node in the analysis largely eliminated these problems, as long as the systems had not yet entered an actual tipping process. Once tipping events occurred in the data, however, all three methods deteriorated markedly, with LKIF producing false positive rates above twenty-five percent, a level that in sparse physical systems could yield more spurious links than genuine ones. The authors therefore advise against trusting any of these methods when an active tipping process is present in the data, and stress that including global temperature as a confounder is crucial whenever a destabilizing forcing acts on the system.

Armed with these guidelines, the team turned to a real-world application that has long divided the literature: the interaction between Arctic summer sea ice and the AMOC. Model studies suggest that a weakening AMOC, which transports less heat northward, should stabilize Arctic sea ice. The reverse direction is contested, because melting sea ice injects freshwater into the North Atlantic, which tends to destabilize the AMOC by reducing buoyancy in the convection regions, while the increased area of exposed ocean surface can lose more heat to the atmosphere, which stabilizes it. Using reanalysis data, an established sea surface temperature fingerprint of the AMOC, and Arctic temperature records as a potential confounder, the researchers applied PCMCI and LKIF with careful preprocessing: detrending, deseasonalizing, spatial filtering, and seasonal masking that restricted the analysis to the March-to-September period when Arctic sea ice is most dynamic.

The results were striking. PCMCI detected a bidirectional stabilizing interaction: a weaker AMOC would increase Arctic summer sea ice concentration, and a loss of sea ice would strengthen the AMOC in the short term, with the strongest effect arriving after a delay of just one month. The strength estimates implied that every ten percentage points of sea ice concentration loss would strengthen the AMOC by roughly 0.61 Sverdrups one month later, while the AMOC’s influence on sea ice, though real, was very weak, amounting to about 0.1 percentage points of sea ice concentration per Sverdrup of circulation change. PCMCI also detected a weaker, delayed destabilizing link at five months, consistent with the slower freshwater mechanism proposed in the literature, although this link did not survive all robustness tests with alternative sea surface temperature datasets. LKIF, hampered by its inability to handle delays and the reduced sample count imposed by masking, detected only the link from the AMOC to sea ice. The authors judged PCMCI the more reliable tool for this application, and the overall picture, a bidirectional stabilizing interaction on monthly timescales, agrees with the physical mechanisms identified by domain experts and model experiments.

The study’s implications extend well beyond the Arctic. It provides the first systematic assessment of how causal inference methods behave on nonlinear data resembling tipping element dynamics, and its recommendations, match sampling rates to interaction delays, always include confounders such as global temperature, avoid analyzing data containing active tipping events, and choose the method whose assumptions fit the problem, offer a practical roadmap for the field. The authors note that the detected short-term effects likely underestimate the full magnitude of the interactions, since slower components of the physical coupling fall outside the observational window, and they point to Earth system model experiments extending beyond 2100 as a promising target for future causal analysis. As the Intergovernmental Panel on Climate Change prepares a dedicated tipping points chapter for its next assessment, work of this kind supplies a much-needed foundation of methodological rigor for a question on which the stability of the climate system may ultimately depend.

Subject of Research: Quantitative comparison of causal inference methods for detecting interactions between climate tipping elements, applied to Arctic sea ice and the Atlantic Meridional Overturning Circulation

Article Title: Quantitative comparison of causal inference methods for climate tipping points

Article References: Lohmann, N., Strahl, D., Högner, A., Huiskamp, W., Boehm, M., & Wunderling, N. (2026). Quantitative comparison of causal inference methods for climate tipping points. Nonlinear Processes in Geophysics, 33(2), 313-334. https://doi.org/10.5194/npg-33-313-2026

Image Credits: AI Generated

DOI: 10.5194/npg-33-313-2026

Keywords: causal inference, climate tipping points, AMOC, Arctic sea ice, PCMCI, Liang-Kleeman information flow, Granger causality, confounders, time series analysis, Earth system, nonlinear dynamics, reanalysis data