Search

Tipping to an undesired state in the climate when a control parameter slowly approaches a critical value is a growing concern with increasing greenhouse gas concentrations. Predictions rely on detecting early warning signals (EWSs) in observations of the system. The primary EWSs are increase in variance, (loss of resilience), and increased autocorrelation (critical slow down). These measures are statistical in nature, which implies that the reliability and statistical significance of the detection depends on the sample size in observations and the magnitude of the change away from the base value prior to the approach to the tipping point. Thus, the possibility of providing useful early warning depends on the relative magnitude of several interdependent time scales in the system. These are (a) the time before the critical value is reached, (b) the (inverse) rate of approach to the tipping point, (c) the size of the time window required to detect a significant change in the EWS and finally, (d) the escape time for noise-induced transition (prior to the tipping). Conditions for early warning of tipping of the Atlantic meridional overturning circulation (AMOC) are marginally fulfilled for the existing past $\sim$150 years of proxy observations where indicators of tipping have recently been reported. Here we provide statistical significance and data driven estimators for the time of tipping. We estimate a collapse of the AMOC to occur around the year 2057 under the assumption of a "business as usual" scenario of future emissions., Comment: 18 pages, 7 figures

Tipping points (TP) in climate sub-systems 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 sub-systems, 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 prior to 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.

The increase in glacial cycle length from approximately $41$ to on average $100$ thousand years around $1$ million years ago, called the Middle Pleistocene Transition (MPT), lacks a conclusive explanation. We describe a dynamical mechanism which we call Ramping with Frequency Locking (RFL), that explains the transition by an interaction between the internal period of a self-sustained oscillator and forcing that contains periodic components. This mechanism naturally explains the abrupt increase in cycle length from approximately $40$ to $80$ thousand years observed in proxy data, unlike some previously proposed mechanisms for the MPT. A rapid increase in durations can be produced by a rapid change in an external parameter, but this assumes rather than explains the abruptness. In contrast, models relying on frequency locking can produce a rapid change in durations assuming only a slow change in an external parameter. We propose a scheme for detecting RFL in complex, computationally expensive models, and motivate the search for climate variables that can gradually increase the internal period of the glacial cycles., Comment: 14 pages, 12 figures

Three-dimensional (3D) turbulence is characterized by a dual forward cascade of both kinetic energy and helicity, a second inviscid flow invariant, from the integral scale of motion to the viscous dissipative scale. In helical flows, however, such as strongly rotating flows with broken mirror symmetry, an inverse energy cascade can be observed analogous to that of two-dimensional turbulence (2D) where a second positive-definite flow invariant, enstrophy, unlike helicity in 3D, effectively blocks the forward cascade of energy. In the spectral-helical decomposition of the Navier-Stokes equation it has previously been show that a subset of three-wave (triad) interactions conserve helicity in 3D in a fashion similar to enstrophy in 2D, thus leading to a 2D-like inverse energy cascade in 3D. In this work, we show both theoretically and numerically that an additional subset of interactions exist conserving a new pseudo-invariant in addition to energy and helicity, which contributes either to a forward or inverse energy cascade depending on the specific triad interaction geometry., Comment: 4 pages, 4 figures

Fully developed homogeneous isotropic turbulence in 2D is fundamentally different from 3D. In 2D, the simultaneous conservation of both energy and enstrophy in the inertial ranges of scales leads to a forward cascade of enstrophy and a reverse cascade of energy. In 3D, helicity, the integral of the scalar product of velocity and vorticity, is also an inviscid flow invariant along with kinetic energy. Unlike enstrophy, however, helicity does not block the cascade of energy to small scales. Energy and helicity are not only globally conserved but also conserved within each non-linear triadic interaction between three plane waves in the spectral form of the Navier--Stokes equation (NSE). By decomposing each plane wave into two helical modes of opposite helicities each triadic interaction is split into a set of eight triadic interactions between helical modes (Waleffe 1992). Biferale et al. (2012) recently found that a subset of these interactions which render both signs of helicity separately conserved (i.e. enstrophy-like) leads to an inverse cascade of (part of) the energy. Motivated by this finding we introduce a new shell model obtained from the NSE expressed in the helical basis (Waleffe 1992). By analysing and integrating the new model we attempt to explain why the dual forward cascade of energy and helicity dominates in 3D turbulence., Comment: 12 pages, 6 figures

Since the beginning of satellite observations, the Arctic sea ice extent has shown a downward trend. The decline has been weaker in the March maximum than in the September minimum and masked by inter-annual fluctuations. One of the less understood aspects of the sea ice response is the persistence times for fluctuations, which could indicate the dominant physical processes behind the sea ice decline. To determine the fluctuation persistence times, however, it is necessary to first filter out the dominant effect of the seasonal cycle. In the current study, we thus develop a statistical model, which accurately decomposes the ice area changes into: (1) a variable seasonal cycle component with a constant shape and (2) a residual (short term) fluctuation. We find the persistence time of fluctuations to be only about three weeks, independently from season, which is substantially shorter than previously reported. Such short time scale points to the dominance of atmospheric forcing. The shape of the seasonal cycle is surprisingly constant for the whole observational record despite the rapid decline. This is in agreement with the suggestion that the asymmetry of the seasonal cycle is an effect of Arctic land-sea geography, which has not changed with climate change. The analysis suggest a jump in the annual sea ice area amplitude occurring in 2007, from which it has not yet recovered, possibly revealing a permanent amplitude shift. In physical sense, this could imply a shift towards the younger, thinner and more susceptible ice cover commencing after the immense 2007 multi-year ice loss., Comment: 8 pages, 4 figures

The variation in the Arctic sea ice is dominated by the seasonal cycle with little inter-annual correlation. Though the mean sea ice area has decreased steadily in the period of satellite observations, a dramatic transition in the dynamics was initiated with the record low September ice area in 2007. The change is much more pronounced in the amplitude of the seasonal cycle than in the annual mean ice area. The shape of the seasonal cycle is surprisingly constant for the whole observational record despite the general decline. A simple explanation, independent of the increased greenhouse warming, for the shape of the seasonal cycle is offered. Thus the dramatic climate change in arctic ice area is seen in the amplitude of the cycle and to a lesser extend the annual mean and the summer ice extend. The reason why the climate change is most pronounced in the amplitude is related to the rapid reduction in perennial ice and thus a thinning of the ice. The analysis shows that a tipping point for the arctic ice area was crossed in 2007., Comment: 9 pages, 4 figures

The dual cascade of energy and enstrophy in 2D turbulence cannot easily be understood in terms of an analog to the Richardson-Kolmogorov scenario describing the energy cascade in 3D turbulence. The coherent up- and downscale fluxes points to non-locality of interactions in spectral space, and thus the specific spatial structure of the flow could be important. Shell models, which lack spacial structure and have only local interactions in spectral space, indeed fail in reproducing the correct scaling for the inverse cascade of energy. In order to exclude the possibility that non-locality of interactions in spectral space is crucial for the dual cascade, we introduce a stochastic spectral model of the cascades which is local in spectral space and which shows the correct scaling for both the direct enstrophy - and the inverse energy cascade., Comment: 4 pages, 3 figures

Here we present a simple stochastic threshold model consisting of a deterministic slowly decaying term and a fast stochastic noise term. The process shows a pseudo-resonance, in the sense that for small and large intensities of the noise the signal is irregular and the distribution of threshold crossings is broad, while for a tuned intermediate value of noise intensity the signal becomes quasi-periodic and the distribution of threshold crossings is narrow. The mechanism captured by the model might be relevant for explaining apparent quasi-periodicity of observed climatic variations where no internal or external periodicities can be identified., Comment: 8 pages, 4 figures, to appear in Stochastics and Dynamics

The glacial cycles are attributed to the climatic response of the orbital changes in the irradiance to the Earth. These changes in the forcing are to small to explain the observed climate variations as simple linear responses. Non-linear amplifications are necessary to account for the glacial cycles. Here an empirical model of the non-linear response is presented. From the model it is possible to assess the role of stochastic noise in comparison to the deterministic orbital forcing of the ice ages. The model is based on the bifurcation structure derived from the climate history. It indicates the dynamical origin of the Mid-Pleistocene transition (MPT) from the '41 kyr world' to the '100 kyr world'. The dominant forcing in the latter is still the 41 kyr obliquity cycle, but the bifurcation structure of the climate system is changed. The model indicates that transitions between glacial and interglacial climate are assisted by internal stochastic noise in the period prior to the last five glacial cycles, while the last five cycles are deterministic responses to the orbital forcing., Comment: 30 pages, 6 figures, submitted to Paleoceanography

By comparing the high-resolution isotopic records from the GRIP and NGRIP icecores, we approximately separate the climate signal from local noise to obtain an objective criterion for defining Dansgaard-Oeschger events. Our analysis identifies several additional short lasting events, increasing the total number of DO events to 27 in the period 12-90 kyr BP. The quasi-regular occurrence of the DO events could indicate a stochastic or coherent resonance mechanism governing their origin. From the distribution of waiting times we obtain a statistical upper bound on the strength of a possible periodic forcing. This finding indicates that the climate shifts are purely noise driven with no underlying periodicity., Comment: 9 figures

The mean surface temperature on Earth and other planets with atmospheres is determined by the radiative balance between the non-reflected incoming solar radiation and the outgoing long-wave black-body radiation from the atmosphere. The surface temperature is higher than the black-body temperature due to the greenhouse warming. Balancing the ice-albedo cooling and the greenhouse warming gives rise to two stable climate states. A cold climate state with a completelyice-covered planet, called Snowball Earth, and a warm state similar to our present climate where greenhouse warming prevents the total glacition. The warm state has dominated Earth in most of its geological history despite a 30 % fainter young Sun. The warming could have been controlled by a greenhouse thermostat operating by temperature control of the weathering process depleting the atmosphere from $CO_2$. This temperature control has permitted life to evolve as early as the end of the heavy bombartment 4 billion years ago., Comment: 4 figures, Proceedings, NORDITA conf. Astrobiology 2004