Discontinuous stochastic forcing in Greenland ice core data

TL;DR

This study uses stochastic differential equations to analyze Greenland ice core proxies, revealing monostable temperature dynamics and bistable atmospheric circulation, offering insights into abrupt climate transitions during the last glacial period.

Contribution

It introduces a data-driven stochastic modeling approach to distinguish the dynamical properties of temperature and circulation proxies in Greenland ice cores.

Findings

Greenland temperature proxy exhibits monostable drift.

Atmospheric circulation proxy shows bistable behavior.

Higher-order stochastic terms may be needed for temperature modeling.

Abstract

Paleoclimate proxy records from Greenland ice cores, archiving e.g. $\delta^{18}$O as a proxy for surface temperature, show that sudden climatic shifts called Dansgaard-Oeschger events (DO) occurred repeatedly during the last glacial interval. They comprised substantial warming of the Arctic region from cold to milder conditions. Concomitant abrupt changes in the dust concentrations of the same ice cores suggest that sudden reorganisations of the hemispheric-scale atmospheric circulation have accompanied the warming events. Genuine bistability of the North Atlantic climate system is commonly hypothesised to explain the existence of stadial (cold) and interstadial (milder) periods in Greenland. However, the physical mechanisms that drove abrupt transitions from the stadial to the interstadial state, and more gradual yet still abrupt reverse transitions, remain debated. Here, we conduct a one-dimensional data-driven analysis of the Greenland temperature and atmospheric circulation proxies under the purview of stochastic processes. We take the Kramers-Moyal equation to estimate each proxy's drift and diffusion terms within a Markovian model framework. We then assess noise contributions beyond Gaussian white noise. The resulting stochastic differential equation (SDE) models feature a monostable drift for the Greenland temperature proxy and a bistable one for the atmospheric circulation proxy. Indicators of discontinuity in stochastic processes suggest to include higher-order terms of the Kramers-Moyal equation when modelling the Greenland temperature proxy's evolution. This constitutes a qualitative difference in the characteristics of the two time series, which should be further investigated from the standpoint of climate dynamics.