Showing 1–50 of 66 results for author: Dijkstra, A

Most Likely Noise-Induced Overturning Circulation Collapse in a 2D Boussinesq Fluid Model

Authors: Jelle Soons, Tobias Grafke, Henk A. Dijkstra

Abstract: There is a reasonable possibility that the present-day Atlantic Meridional Overturning Circulation is in a bi-stable regime and hence it is relevant to compute probabilities and pathways of noise-induced transitions between the stable equilibrium states. Here, the most probable transition pathway of a noise-induced collapse of the northern overturning circulation in a spatially-continuous two-dime… ▽ More There is a reasonable possibility that the present-day Atlantic Meridional Overturning Circulation is in a bi-stable regime and hence it is relevant to compute probabilities and pathways of noise-induced transitions between the stable equilibrium states. Here, the most probable transition pathway of a noise-induced collapse of the northern overturning circulation in a spatially-continuous two-dimensional model with surface temperature and stochastic salinity forcings is directly computed using Large Deviation Theory (LDT). This pathway reveals the fluid dynamical mechanisms of such a collapse. Paradoxically it starts off with a strengthening of the northern overturning circulation before a short but strong salinity pulse induces a second overturning cell. The increased atmospheric energy input of this two-cell configuration cannot be mixed away quickly enough, leading to the collapse of the northern overturning cell and finally resulting in a southern overturning circulation. Additionally, the approach allows us to compare the probability of this collapse under different parameters in the deterministic part of the salinity surface forcing, which quantifies the increase in collapse probability as the bifurcation point of the system is approached. △ Less

Submitted 26 August, 2024; originally announced August 2024.

Machine Learning for the Physics of Climate

Authors: Annalisa Bracco, Julien Brajard, Henk A. Dijkstra, Pedram Hassanzadeh, Christian Lessig, Claire Monteleoni

Abstract: An exponential growth in computing power, which has brought more sophisticated and higher resolution simulations of the climate system, and an exponential increase in observations since the first weather satellite was put in orbit, are revolutionizing climate science. Big data and associated algorithms, coalesced under the field of Machine Learning (ML), offer the opportunity to study the physics… ▽ More An exponential growth in computing power, which has brought more sophisticated and higher resolution simulations of the climate system, and an exponential increase in observations since the first weather satellite was put in orbit, are revolutionizing climate science. Big data and associated algorithms, coalesced under the field of Machine Learning (ML), offer the opportunity to study the physics of the climate system in ways, and with an amount of detail, infeasible few years ago. The inference provided by ML has allowed to ask causal questions and improve prediction skills beyond classical barriers. Furthermore, when paired with modeling experiments or robust research in model parameterizations, ML is accelerating computations, increasing accuracy and allowing for generating very large ensembles at a fraction of the cost. In light of the urgency imposed by climate change and the rapidly growing role of ML, we review its broader accomplishments in climate physics. Decades long standing problems in observational data reconstruction, representation of sub-grid scale phenomena and climate (and weather) prediction are being tackled with new and justified optimism. Ultimately, this review aims at providing a perspective on the benefits and major challenges of exploiting ML in studying complex systems. △ Less

Submitted 18 August, 2024; originally announced August 2024.

Comments: 25 pages, 4 figures

Quantification of salt transports due to exchange flow and tidal flow in estuaries

Authors: Bouke Biemond, Huib E. de Swart, Henk A. Dijkstra

Abstract: To understand mechanisms of salt intrusion in estuaries, we develop a semi-analytical model, that explicitly accounts for salt transport by both exchange flow and tidal flow. This model, after calibration, successfully hindcasts hydrodynamics and salinity dynamics in three estuaries that have strongly different characteristics. We find, from analyzing the model results for these three estuaries, t… ▽ More To understand mechanisms of salt intrusion in estuaries, we develop a semi-analytical model, that explicitly accounts for salt transport by both exchange flow and tidal flow. This model, after calibration, successfully hindcasts hydrodynamics and salinity dynamics in three estuaries that have strongly different characteristics. We find, from analyzing the model results for these three estuaries, that salt transport processes by exchange flow and tidal flow interact through the subtidal stratification. Transport by exchange flow creates stratification, thereby generating a phase shift of tidal salinity with respect to the tidal flow, which is important for the magnitude of the tidal salt transport. Conversely, the strength of tidal currents determines the vertical mixing that breaks down stratification. A new analytical formulation is presented for the component of the salt transport driven by the depth-averaged tidal flow. This salt transport is larger than the component associated with the vertical shear of the tidal current. Finally, a method that yields analytical equations that quantify the importance of different contributions to the salt transport using only primary information is developed using approximate solutions for the subtidal stratification. This method performs well for the estuaries considered. △ Less

Submitted 9 August, 2024; originally announced August 2024.

Substantial Risk of 21st Century AMOC Tipping even under Moderate Climate Change

Authors: René M. van Westen, Elian Y. P. Vanderborght, Michael Kliphuis, Henk A. Dijkstra

Abstract: The Atlantic Meridional Overturning Circulation (AMOC) is a key component of the climate system and considered to be a tipping element. There is still a large uncertainty on the critical global warming level at which the AMOC will start to collapse. Here we analyse targeted climate model simulations, together with observations, reanalysis products and a suite of state-of-the-art climate model resu… ▽ More The Atlantic Meridional Overturning Circulation (AMOC) is a key component of the climate system and considered to be a tipping element. There is still a large uncertainty on the critical global warming level at which the AMOC will start to collapse. Here we analyse targeted climate model simulations, together with observations, reanalysis products and a suite of state-of-the-art climate model results to reassess this critical global warming level. We find a critical threshold of +3C global mean surface temperature increase compared to pre-industrial with a lower bound of +2.2C (10%-Cl). Such global mean surface temperature anomalies are expected to be reached after 2050. This means that the AMOC is more likely than not (> 50%) to tip within the 21st century under a middle-of-the-road climate change scenario and very likely (> 90%) to tip under a high emissions scenario. The AMOC collapse induced cooling is shown to be offset by the regional warming over Northwestern Europe during the 21st century, but will still induce severe impacts on society. △ Less

Submitted 29 July, 2024; originally announced July 2024.

Resilience of the Atlantic Meridional Overturning Circulation

Authors: Valérian Jacques-Dumas, Henk A. Dijkstra, Christian Kuehn

Abstract: We address the issue of resilience of the Atlantic Meridional Overturning Circulation (AMOC) given the many indications that this dynamical system is in a multi-stable regime. A novel approach to resilience based on rare event techniques is presented which leads to a measure capturing `resistance to change` and `ability to return' aspects in a probabilistic way. The application of this measure to… ▽ More We address the issue of resilience of the Atlantic Meridional Overturning Circulation (AMOC) given the many indications that this dynamical system is in a multi-stable regime. A novel approach to resilience based on rare event techniques is presented which leads to a measure capturing `resistance to change` and `ability to return' aspects in a probabilistic way. The application of this measure to a conceptual model demonstrates its suitability for assessing AMOC resilience but also shows its potential use in many other non-autonomous dynamical systems. This framework is then extended to compute the probability that the AMOC undergoes a transition conditioned on an external forcing. Such conditional probability can be estimated by exploiting the information available when computing the resilience of this system. This allows us to provide a probabilistic view on safe operating spaces by defining a conditional safe operating space as a subset of the parameter space of the (possibly transient) imposed forcing. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 35 pages, 8 figures, 1 table

Probability Estimates of a 21st Century AMOC Collapse

Authors: Emma J. V. Smolders, René M. van Westen, Henk A. Dijkstra

Abstract: There is increasing concern that the Atlantic Meridional Overturning Circulation (AMOC) may collapse this century with a disrupting societal impact on large parts of the world. Preliminary estimates of the probability of such an AMOC collapse have so far been based on conceptual models and statistical analyses of proxy data. Here, we provide observationally based estimates of such probabilities fr… ▽ More There is increasing concern that the Atlantic Meridional Overturning Circulation (AMOC) may collapse this century with a disrupting societal impact on large parts of the world. Preliminary estimates of the probability of such an AMOC collapse have so far been based on conceptual models and statistical analyses of proxy data. Here, we provide observationally based estimates of such probabilities from reanalysis data. We first identify optimal observation regions of an AMOC collapse from a recent global climate model simulation. Salinity data near the southern boundary of the Atlantic turn out to be optimal to provide estimates of the time of the AMOC collapse in this model. Based on the reanalysis products, we next determine probability density functions of the AMOC collapse time. The collapse time is estimated between 2037-2064 (10-90% CI) with a mean of 2050 and the probability of an AMOC collapse before the year 2050 is estimated to be 59 +/- 17%. △ Less

Submitted 17 June, 2024; originally announced June 2024.

A detailed analysis of the origin of deep-decoupling oscillations

Authors: John Bailie, Henk A. Dijkstra, Bernd Krauskopf

Abstract: The variability of the strength of the Atlantic Meridional Overturning Circulation is influenced substantially by the formation of deep water in the North Atlantic. In many ocean models, so-called deep-decoupling oscillations have been found, whose timescale depends on the characteristics of convective vertical mixing processes. Their precise origin and sensitivity to the representation of mixing… ▽ More The variability of the strength of the Atlantic Meridional Overturning Circulation is influenced substantially by the formation of deep water in the North Atlantic. In many ocean models, so-called deep-decoupling oscillations have been found, whose timescale depends on the characteristics of convective vertical mixing processes. Their precise origin and sensitivity to the representation of mixing have remained unclear so far. To study this problem, we revisit a conceptual Welander model for the evolution of temperature and salinity in two vertically stacked boxes for surface and deep water, which interact through diffusion and/or convective adjustment. The model is known to exhibit several types of deep-decoupling oscillations, with phases of weak diffusive mixing interspersed with strong convective mixing, when the switching between them is assumed to be instantaneous. We present a comprehensive study of oscillations in Welander's model with non-instantaneous switching between mixing phases, as described by a smooth switching function. A dynamical systems approach allows us to distinguish four types of oscillations, in terms of their phases of diffusive versus convective mixing, and to identify the regions in the relevant parameter plane where they exist. The characteristic deep-decoupling oscillations still exist for non-instantaneous switching, but they require switching that is considerably faster than needed for sustaining oscillatory behaviour. Furthermore, we demonstrate how a gradual freshwater influx can lead to transitions between different vertical mixing oscillations. Notably, the convective mixing phase becomes shorter and even disappears, resulting in long periods of much reduced deep water formation. The results are relevant for the interpretation of ocean-climate variability in models and (proxy) observations. △ Less

Submitted 2 June, 2024; originally announced June 2024.

Warning Signs for Boundary Noise and their Application to an Ocean Boussinesq Model

Authors: Paolo Bernuzzi, Henk A. Dijkstra, Christian Kuehn

Abstract: In this paper, we construct and discuss early-warning signs of the approach of a parameter to a deterministic bifurcation on a stochastic partial differential equation (SPDE) model with Gaussian white-noise on the boundary. We specifically focus on the qualitative behaviour of the time-asymptotic autocovariance and autocorrelation of the solutions of the linearized system. We also discuss the reli… ▽ More In this paper, we construct and discuss early-warning signs of the approach of a parameter to a deterministic bifurcation on a stochastic partial differential equation (SPDE) model with Gaussian white-noise on the boundary. We specifically focus on the qualitative behaviour of the time-asymptotic autocovariance and autocorrelation of the solutions of the linearized system. We also discuss the reliability of the tools from an analytic perspective and through various examples. Among those, the application of the early-warning signs for an ocean Boussinesq model is explored through numerical simulations. The analytic results obtained expand on previous work and present valuable early-warning signs for various applications. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 24 pages, 40 figures

MSC Class: 60H15 (Primary); 35R60; 37H20 (Secondary)

The Role of Sea-ice Processes on the Probability of AMOC Transitions

Authors: René M. van Westen, Valérian Jacques-Dumas, Amber A. Boot, Henk A. Dijkstra

Abstract: Recent simulations performed with the Community Earth System Model (CESM) have suggested a crucial role of sea-ice processes in AMOC hysteresis behaviour under varying surface freshwater forcing. Here, we further investigate this issue using additional CESM simulations and a novel conceptual ocean-sea-ice box model. The CESM simulations show that the presence of sea ice gives rise to the existence… ▽ More Recent simulations performed with the Community Earth System Model (CESM) have suggested a crucial role of sea-ice processes in AMOC hysteresis behaviour under varying surface freshwater forcing. Here, we further investigate this issue using additional CESM simulations and a novel conceptual ocean-sea-ice box model. The CESM simulations show that the presence of sea ice gives rise to the existence of statistical equilibrium states with a weak AMOC strength. This is confirmed in the conceptual model, which captures the same AMOC hysteresis behaviour as in the CESM simulation and where steady states are computed versus forcing parameters. In the conceptual model, transition probabilities between the different equilibrium states are determined using rare event techniques. The transition probabilities from a strong AMOC state to a weak AMOC state increase when considering sea-ice processes and indicate that sea ice promotes these transitions. On the other hand, sea ice strongly reduces the probabilities of the reverse transition from a weak AMOC state to a strong AMOC state and this implies that sea ice also limits AMOC recovery. The results here indicate that sea-ice processes play a dominant role in AMOC hysteresis width and influence transition probabilities between the different equilibrium states. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Estimation of AMOC transition probabilities using a machine learning based rare-event algorithm

Authors: Valérian Jacques-Dumas, René M. van Westen, Henk A. Dijkstra

Abstract: The Atlantic Meridional Overturning Circulation (AMOC) is an important component of the global climate, known to be a tipping element, as it could collapse under global warming. The main objective of this study is to compute the probability that the AMOC collapses within a specified time window, using a rare-event algorithm called Trajectory-Adaptive Multilevel Splitting (TAMS). However, the effic… ▽ More The Atlantic Meridional Overturning Circulation (AMOC) is an important component of the global climate, known to be a tipping element, as it could collapse under global warming. The main objective of this study is to compute the probability that the AMOC collapses within a specified time window, using a rare-event algorithm called Trajectory-Adaptive Multilevel Splitting (TAMS). However, the efficiency and accuracy of TAMS depend on the choice of the score function. Although the definition of the optimal score function, called ``committor function" is known, it is impossible in general to compute it a priori. Here, we combine TAMS with a Next-Generation Reservoir Computing technique that estimates the committor function from the data generated by the rare-event algorithm. We test this technique in a stochastic box model of the AMOC for which two types of transition exist, the so-called F(ast)-transitions and S(low)-transitions. Results for the F-transtions compare favorably with those in the literature where a physically-informed score function was used. We show that coupling a rare-event algorithm with machine learning allows for a correct estimation of transition probabilities, transition times, and even transition paths for a wide range of model parameters. We then extend these results to the more difficult problem of S-transitions in the same model. In both cases of F-transitions and S-transitions, we also show how the Next-Generation Reservoir Computing technique can be interpreted to retrieve an analytical estimate of the committor function. △ Less

Submitted 2 July, 2024; v1 submitted 19 January, 2024; originally announced January 2024.

Comments: 15 pages, 9 figures

Optimal Transition Paths for AMOC Collapse and Recovery in a Stochastic Box Model

Authors: Jelle Soons, Tobias Grafke, Henk A. Dijkstra

Abstract: There is strong evidence that the present-day Atlantic Meridional Overturning Circulation (AMOC) is in a bi-stable regime and hence it is important to determine probabilities and pathways for noise-induced transitions between its equilibrium states. Here, using Large Deviation Theory (LDT), the most probable transition pathways for the noise-induced collapse and recovery of the AMOC are computed i… ▽ More There is strong evidence that the present-day Atlantic Meridional Overturning Circulation (AMOC) is in a bi-stable regime and hence it is important to determine probabilities and pathways for noise-induced transitions between its equilibrium states. Here, using Large Deviation Theory (LDT), the most probable transition pathways for the noise-induced collapse and recovery of the AMOC are computed in a stochastic box model of the World Ocean. This allows us to determine the physical mechanisms of noise-induced AMOC transitions. We show that the most likely path of an AMOC collapse starts paradoxically with a strengthening of the AMOC followed by an immediate drop within a couple of years due to a short but relatively strong freshwater pulse. The recovery on the other hand is a slow process, where the North Atlantic needs to be gradually salinified over a course of 20 years. The proposed method provides several benefits, including an estimate of probability ratios of collapse between various freshwater noise scenarios, showing that the AMOC is most vulnerable to freshwater forcing into the Atlantic thermocline region. Moreover, a comparison with a quasi-equilibrium approach reveals the contrasts in behavior of a bifurcation-induced and a noise-induced collapse of the AMOC. △ Less

Submitted 26 August, 2024; v1 submitted 21 November, 2023; originally announced November 2023.

Asymmetry of AMOC Hysteresis in a State-of-the-Art Global Climate Model

Authors: René M. van Westen, Henk A. Dijkstra

Abstract: We study hysteresis properties of the Atlantic Meridional Overturning Circulation (AMOC) under a slowly-varying North Atlantic (20$^{\circ}$N -- 50$^{\circ}$N) freshwater flux forcing in state-of-the-art Global Climate Model (GCM), the Community Earth System Model. Results are presented of a full hysteresis simulation ($4,400$ model years) and show that there is a hysteresis width of about $0.4$ S… ▽ More We study hysteresis properties of the Atlantic Meridional Overturning Circulation (AMOC) under a slowly-varying North Atlantic (20$^{\circ}$N -- 50$^{\circ}$N) freshwater flux forcing in state-of-the-art Global Climate Model (GCM), the Community Earth System Model. Results are presented of a full hysteresis simulation ($4,400$ model years) and show that there is a hysteresis width of about $0.4$ Sv. This demonstrates that an AMOC collapse and recovery do not only occur in conceptual and idealised climate models, but also in a state-of-the-art GCM. The AMOC recovery is about a factor six faster than the AMOC collapse and this asymmetry is due to the major effect of the North Atlantic sea-ice distribution on the AMOC recovery. The results have implications for projections of possible future AMOC behaviour and for explaining relatively rapid climate transitions in the geological past. △ Less

Submitted 27 August, 2023; originally announced August 2023.

The effect of model freshwater flux biases on the multi-stable regime of the AMOC

Authors: Henk A. Dijkstra, Rene M. van Westen

Abstract: It is known that global climate models (GCMs) have substantial biases in the surface freshwater flux which forces the ocean component of these models. Using numerical bifurcation analyses on a global ocean model, we study here the effect of a specific freshwater flux bias on the multiple equilibrium regime of the Atlantic Meridional Overturning Circulation (AMOC). We find that a (positive) freshwa… ▽ More It is known that global climate models (GCMs) have substantial biases in the surface freshwater flux which forces the ocean component of these models. Using numerical bifurcation analyses on a global ocean model, we study here the effect of a specific freshwater flux bias on the multiple equilibrium regime of the Atlantic Meridional Overturning Circulation (AMOC). We find that a (positive) freshwater flux bias over the Indian Ocean shifts the multiple equilibrium regime to larger values of North Atlantic freshwater input but hardly affects the associated hysteresis width. The magnitude of this shift depends on the way the anomalous North Atlantic freshwater flux is compensated. We explain the changes in bifurcation diagrams using the freshwater balance over the Atlantic basin. The results suggest that state-of-the-art GCMs may have an AMOC multiple equilibrium regime, but that it is located in a parameter regime that is considered unrealistic and hence is not explored. △ Less

Submitted 22 August, 2023; originally announced August 2023.

New Physics-Based Early Warning Signal shows AMOC is on Tipping Course

Authors: René M. van Westen, Michael Kliphuis, Henk A. Dijkstra

Abstract: One of the most prominent climate tipping elements is the Atlantic Meridional Overturning Circulation (AMOC), which can potentially collapse due to the input of fresh water in the North Atlantic. Although AMOC collapses have been induced in state-of-the-art Global Climate Models (GCMs) by strong freshwater forcing, an AMOC tipping event has so far not been found. Here, we show results of the first… ▽ More One of the most prominent climate tipping elements is the Atlantic Meridional Overturning Circulation (AMOC), which can potentially collapse due to the input of fresh water in the North Atlantic. Although AMOC collapses have been induced in state-of-the-art Global Climate Models (GCMs) by strong freshwater forcing, an AMOC tipping event has so far not been found. Here, we show results of the first AMOC tipping event in such a GCM, the Community Earth System Model, including the large climate impacts of the associated AMOC collapse. Using these results, we develop a new physics-based and observable early warning signal of AMOC tipping: the minimum of the AMOC induced freshwater transport at the southern boundary of the Atlantic. Available observational and reanalysis data indicate that the present-day AMOC is on route to tipping. The new early warning signal is a useful alternative to classical statistical ones which, when applied to our simulated AMOC tipping event, turn out to be sensitive to the analysed time interval before tipping. △ Less

Submitted 3 August, 2023; originally announced August 2023.

Data-driven methods to estimate the committor function in conceptual ocean models

Authors: Valérian Jacques-Dumas, René M. van Westen, Freddy Bouchet, Henk A. Dijkstra

Abstract: In recent years, several climate subsystems have been identified that may undergo a relatively rapid transition compared to the changes in their forcing. Such transitions are rare events in general, and simulating long-enough trajectories in order to gather sufficient data to determine transition statistics would be too expensive. Conversely, rare events algorithms like TAMS (trajectory-adaptive m… ▽ More In recent years, several climate subsystems have been identified that may undergo a relatively rapid transition compared to the changes in their forcing. Such transitions are rare events in general, and simulating long-enough trajectories in order to gather sufficient data to determine transition statistics would be too expensive. Conversely, rare events algorithms like TAMS (trajectory-adaptive multilevel sampling) encourage the transition while keeping track of the model statistics. However, this algorithm relies on a score function whose choice is crucial to ensure its efficiency. The optimal score function, called the committor function, is in practice very difficult to compute. In this paper, we compare different data-based methods (analog Markov chains, neural networks, reservoir computing, dynamical Galerkin approximation) to estimate the committor from trajectory data. We apply these methods on two models of the Atlantic Ocean circulation featuring very different dynamical behavior. We compare these methods in terms of two measures, evaluating how close the estimate is from the true committor and in terms of the computational time. We find that all methods are able to extract information from the data in order to provide a good estimate of the committor. Analog Markov Chains provide a very reliable estimate of the true committor in simple models but prove not so robust when applied to systems with a more complex phase space. Neural network methods clearly stand out by their relatively low testing time, and their training time scales more favorably with the complexity of the model than the other methods. In particular, feedforward neural networks consistently achieve the best performance when trained with enough data, making this method promising for committor estimation in sophisticated climate models. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: 22 pages, 5 figures

Journal ref: Nonlin. Processes Geophys. 30 (2023) 195-216

Bifurcation analysis of a North Atlantic Ocean box model with two deep-water formation sites

Authors: Alannah Neff, Andrew Keane, Henk A. Dijkstra, Bernd Krauskopf

Abstract: The tipping of the Atlantic Meridional Overturning Circulation (AMOC) to a 'shutdown' state due to changes in the freshwater forcing of the ocean is of particular interest and concern due to its widespread ramifications, including a dramatic climatic shift for much of Europe. A clear understanding of how such a shutdown would unfold requires analyses of models from across the complexity spectrum.… ▽ More The tipping of the Atlantic Meridional Overturning Circulation (AMOC) to a 'shutdown' state due to changes in the freshwater forcing of the ocean is of particular interest and concern due to its widespread ramifications, including a dramatic climatic shift for much of Europe. A clear understanding of how such a shutdown would unfold requires analyses of models from across the complexity spectrum. For example, detailed simulations of sophisticated Earth System Models have identified scenarios in which deep-water formation first ceases in the Labrador Sea before ceasing in the Nordic Seas, en route to a complete circulation shutdown. Here, we study a simple ocean box model with two polar boxes designed to represent deep-water formation at these two distinct sites. A bifurcation analysis reveals how, depending on the differences of freshwater and thermal forcing between the two polar boxes, transitions to 'partial shutdown' states are possible. Our results shed light on the nature of the tipping of AMOC and clarify dynamical features observed in more sophisticated models. △ Less

Submitted 8 June, 2023; v1 submitted 19 May, 2023; originally announced May 2023.

Comments: 13 pages, 13 figures

Multistability and Intermediate Tipping of the Atlantic Ocean Circulation

Authors: Johannes Lohmann, Henk A. Dijkstra, Markus Jochum, Valerio Lucarini, Peter D. Ditlevsen

Abstract: 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 low… ▽ More 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. △ Less

Submitted 15 November, 2023; v1 submitted 12 April, 2023; originally announced April 2023.

Physics captured by data-based methods in El Niño prediction

Authors: G. Lancia, I. J. Goede, C. Spitoni, H. A. Dijkstra

Abstract: On average once every four years, the Tropical Pacific warms considerably during events called El Niño, leading to weather disruptions over many regions on Earth. Recent machine-learning approaches to El Niño prediction, in particular Convolutional Neural Networks (CNNs), have shown a surprisingly high skill at relatively long lead times. In an attempt to understand this high skill, we here use da… ▽ More On average once every four years, the Tropical Pacific warms considerably during events called El Niño, leading to weather disruptions over many regions on Earth. Recent machine-learning approaches to El Niño prediction, in particular Convolutional Neural Networks (CNNs), have shown a surprisingly high skill at relatively long lead times. In an attempt to understand this high skill, we here use data from distorted physics simulations with an intermediate complexity El Niño model to determine what aspects of El Niño physics are represented in a specific CNN-based classification method. We find that the CNN can adequately correct for distortions in the ocean adjustment processes, but that the machine-learning method has far more trouble to deal with distortions in upwelling feedback strength. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 17 pages, 12 figures

A simple mechanism for stable oscillations in the large-scale ocean circulation

Authors: Andrew Keane, Alexandre Pohl, Henk A. Dijkstra, Andy Ridgwell

Abstract: The global ocean circulation plays a pivotal role in the regulation of the Earth's climate. The specific pattern and strength of circulation also determines how carbon and nutrients are cycled and via the resulting distribution of dissolved oxygen, where habitats suitable for marine animals occur. However, evidence from both geological data and models suggests that state transitions in circulation… ▽ More The global ocean circulation plays a pivotal role in the regulation of the Earth's climate. The specific pattern and strength of circulation also determines how carbon and nutrients are cycled and via the resulting distribution of dissolved oxygen, where habitats suitable for marine animals occur. However, evidence from both geological data and models suggests that state transitions in circulation patterns have occurred in the past. Understanding the controls on marine environmental conditions and biological productivity then requires a full appreciation of the nature and drivers of such transitions. Here we present an analysis of stable millennial oscillations that occur in the ocean circulation of an Earth System Model of intermediate complexity in the presence of a circumpolar current. To demonstrate that a circumpolar current can act as a driver of stable oscillations, we adapt a simple ocean box model to include a delayed feedback to represent a circumpolar current and investigate the stable millennial oscillatory solutions that arise in the box model by bifurcation analysis. Our results provide new insights into the nature of oscillations that could have occurred under certain continental configurations in the geological past, and also highlight the potential influence of changing circumpolar current speed on the stability of the ocean's meridional overturning circulation. △ Less

Submitted 18 November, 2022; v1 submitted 19 January, 2022; originally announced January 2022.

Comments: 13 pages, 8 figures

Partial tipping in a spatially heterogeneous world

Authors: Robbin Bastiaansen, Henk A. Dijkstra, Anna S. von der Heydt

Abstract: Many climate subsystems are thought to be susceptible to tipping - and some might be close to a tipping point. The general belief and intuition, based on simple conceptual models of tipping elements, is that tipping leads to reorganization of the full (sub)system. Here, we explore tipping in conceptual, but spatially extended and spatially heterogenous models. These are extensions of conceptual mo… ▽ More Many climate subsystems are thought to be susceptible to tipping - and some might be close to a tipping point. The general belief and intuition, based on simple conceptual models of tipping elements, is that tipping leads to reorganization of the full (sub)system. Here, we explore tipping in conceptual, but spatially extended and spatially heterogenous models. These are extensions of conceptual models taken from all sorts of climate system components on multiple spatial scales. By analysis of the bifurcation structure of such systems, special stable equilibrium states are revealed: coexistence states with part of the spatial domain in one state, and part in another, with a spatial interface between these regions. These coexistence states critically depend on the size and the spatial heterogeneity of the (sub)system. In particular, in these systems a tipping point might lead to a partial tipping of the full (sub)system, in which only part of the spatial domain undergoes reorganization, limiting the impact of these events on the system's functioning. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Comments: 18 pages, 7 figures

Projections of the Transient State-Dependency of Climate Feedbacks

Authors: Robbin Bastiaansen, Henk A. Dijkstra, Anna S. von der Heydt

Abstract: When the climate system is forced, e.g. by emission of greenhouse gases, it responds on multiple time scales. As temperatures rise, feedback processes might intensify or weaken. Current methods to analyze feedback strength, however, do not take such state dependency into account; they only consider changes in (global mean) temperature and assume all feedbacks are linearly related to that. This mak… ▽ More When the climate system is forced, e.g. by emission of greenhouse gases, it responds on multiple time scales. As temperatures rise, feedback processes might intensify or weaken. Current methods to analyze feedback strength, however, do not take such state dependency into account; they only consider changes in (global mean) temperature and assume all feedbacks are linearly related to that. This makes (transient) changes in feedback strengths almost intangible and generally leads to underestimation of future warming. Here, we present a multivariate (and spatially explicit) framework that facilitates dissection of climate feedbacks over time scales. Using this framework, information on the composition of projected (transient) future climates and feedback strengths can be obtained. Moreover, it can be used to make projections for many emission scenarios through linear response theory. The new framework is illustrated using the Community Earth System Model version 2 (CESM2). △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: main text: 11 pages, 4 figures, 1 table Supporting Information: 14 pages, 17 figures, 1 table, 8 movies

Eddy memory as an explanation of intra-seasonal periodic behavior in baroclinic eddies

Authors: Woosok Moon, Georgy E. Manucharyan, Henk A. Dijkstra

Abstract: The baroclinic annular mode (BAM) is a leading-order mode of the eddy-kinetic energy in the Southern Hemisphere exhibiting. oscillatory behavior at intra-seasonal time scales. The oscillation mechanism has been linked to transient eddy-mean flow interactions that remain poorly understood. Here we demonstrate that the finite memory effect in eddy-heat flux dependence on the large-scale flow can exp… ▽ More The baroclinic annular mode (BAM) is a leading-order mode of the eddy-kinetic energy in the Southern Hemisphere exhibiting. oscillatory behavior at intra-seasonal time scales. The oscillation mechanism has been linked to transient eddy-mean flow interactions that remain poorly understood. Here we demonstrate that the finite memory effect in eddy-heat flux dependence on the large-scale flow can explain the origin of the BAM's oscillatory behavior. We represent the eddy memory effect by a delayed integral kernel that leads to a generalized Langevin equation for the planetary-scale heat equation. Using a mathematical framework for the interactions between planetary and synoptic-scale motions, we derive a reduced dynamical model of the BAM - a stochastically-forced oscillator with a period proportional to the geometric mean between the eddy-memory time scale and the diffusive eddy equilibration timescale. Our model provides a formal justification for the previously proposed phenomenological model of the BAM and could be used to explicitly diagnose the memory kernel and improve our understanding of transient eddy-mean flow interactions in the atmosphere. △ Less

Submitted 8 February, 2021; originally announced February 2021.

Comments: 12 pages, 1 figure

Quantification and interpretation of the climate variability record

Authors: Anna S. von der Heydt, Peter Ashwin, Charles D. Camp, Michel Crucifix, Henk A. Dijkstra, Peter Ditlevsen, Timothy M. Lenton

Abstract: The spectral view of variability is a compelling and adaptable tool for understanding variability of the climate. In Mitchell (1976) seminal paper, it was used to express, on one graph with log scales, a very wide range of climate variations from millions of years to days. The spectral approach is particularly useful for suggesting causal links between forcing variability and climate response vari… ▽ More The spectral view of variability is a compelling and adaptable tool for understanding variability of the climate. In Mitchell (1976) seminal paper, it was used to express, on one graph with log scales, a very wide range of climate variations from millions of years to days. The spectral approach is particularly useful for suggesting causal links between forcing variability and climate response variability. However, a substantial degree of variability is intrinsic and the Earth system may respond to external forcing in a complex manner. There has been an enormous amount of work on understanding climate variability over the last decades. Hence in this paper, we address the question: Can we (after 40 years) update the Mitchell (1976) diagram and provide it with a better interpretation? By reviewing both the extended observations available for such a diagram and new methodological developments in the study of the interaction between internal and forced variability over a wide range of timescales, we give a positive answer to this question. In addition, we review alternative approaches to the spectral decomposition and pose some challenges for a more detailed quantification of climate variability. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Journal ref: Global and Planetary Change (2021), Vol. 197, article number 103399

Application of Adaptive Multilevel Splitting to High-Dimensional Dynamical Systems

Authors: S. Baars, D. Castellana, F. W. Wubs, H. A. Dijkstra

Abstract: Stochastic nonlinear dynamical systems can undergo rapid transitions relative to the change in their forcing, for example due to the occurrence of multiple equilibrium solutions for a specific interval of parameters. In this paper, we modify one of the methods developed to compute probabilities of such transitions, Trajectory-Adaptive Multilevel Sampling (TAMS), to be able to apply it to high-dime… ▽ More Stochastic nonlinear dynamical systems can undergo rapid transitions relative to the change in their forcing, for example due to the occurrence of multiple equilibrium solutions for a specific interval of parameters. In this paper, we modify one of the methods developed to compute probabilities of such transitions, Trajectory-Adaptive Multilevel Sampling (TAMS), to be able to apply it to high-dimensional systems. The key innovation is a projected time-stepping approach, which leads to a strong reduction in computational costs, in particular memory usage. The performance of this new implementation of TAMS is studied through an example of the collapse of the Atlantic Ocean Circulation. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Continuation of Probability Density Functions using a Generalized Lyapunov Approach

Authors: S. Baars, J. P. Viebahn, T. E. Mulder, C. Kuehn, F. W. Wubs, H. A. Dijkstra

Abstract: Techniques from numerical bifurcation theory are very useful to study transitions between steady fluid flow patterns and the instabilities involved. Here, we provide computational methodology to use parameter continuation in determining probability density functions of systems of stochastic partial differential equations near fixed points, under a small noise approximation. Key innovation is the e… ▽ More Techniques from numerical bifurcation theory are very useful to study transitions between steady fluid flow patterns and the instabilities involved. Here, we provide computational methodology to use parameter continuation in determining probability density functions of systems of stochastic partial differential equations near fixed points, under a small noise approximation. Key innovation is the efficient solution of a generalized Lyapunov equation using an iterative method involving low-rank approximations. We apply and illustrate the capabilities of the method using a problem in physical oceanography, i.e. the occurrence of multiple steady states of the Atlantic Ocean circulation. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Journal ref: Journal of Computational Physics 336 (2017) 627-643

Modelling size distributions of marine plastics under the influence of continuous cascading fragmentation

Authors: Mikael L. A. Kaandorp, Henk A. Dijkstra, Erik van Sebille

Abstract: Field studies have shown that plastic fragments make up the majority of plastic pollution in the oceans in terms of abundance. How quickly environmental plastics fragment is not well understood, however. Here, we study this fragmentation process by considering a model which captures continuous fragmentation of particles over time in a cascading fashion. With this cascading fragmentation model we s… ▽ More Field studies have shown that plastic fragments make up the majority of plastic pollution in the oceans in terms of abundance. How quickly environmental plastics fragment is not well understood, however. Here, we study this fragmentation process by considering a model which captures continuous fragmentation of particles over time in a cascading fashion. With this cascading fragmentation model we simulate particle size distributions (PSDs), specifying the abundance or mass of particles for different size classes. The fragmentation model is coupled to an environmental box model, simulating the distributions of plastic particles in the ocean, coastal waters, and on the beach. We compare the modelled PSDs to available observations, and use the results to illustrate the effect of size-selective processes such as vertical mixing in the water column and resuspension of particles from the beach into coastal waters. The model quantifies the role of fragmentation on the marine plastic mass budget: while fragmentation is a major source of (secondary) plastic particles in terms of abundance, it seems to have a minor effect on the total mass of particles larger than 0.1 mm. Future comparison to observed PSD data allow us to understand size-selective plastic transport in the environment, and potentially inform us on plastic longevity. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Multivariate Estimations of Equilibrium Climate Sensitivity from Short Transient Warming Simulations

Authors: Robbin Bastiaansen, Henk A. Dijkstra, Anna S. von der Heydt

Abstract: One of the most used metrics to gauge the effects of climate change is the equilibrium climate sensitivity, defined as the long-term (equilibrium) temperature increase resulting from instantaneous doubling of atmospheric CO$_2$. Since global climate models cannot be fully equilibrated in practice, extrapolation techniques are used to estimate the equilibrium state from transient warming simulation… ▽ More One of the most used metrics to gauge the effects of climate change is the equilibrium climate sensitivity, defined as the long-term (equilibrium) temperature increase resulting from instantaneous doubling of atmospheric CO$_2$. Since global climate models cannot be fully equilibrated in practice, extrapolation techniques are used to estimate the equilibrium state from transient warming simulations. Because of the abundance of climate feedbacks - spanning a wide range of temporal scales - it is hard to extract long-term behaviour from short-time series; predominantly used techniques are only capable of detecting the single most dominant eigenmode, thus hampering their ability to give accurate long-term estimates. Here, we present an extension to those methods by incorporating data from multiple observables in a multi-component linear regression model. This way, not only the dominant but also the next-dominant eigenmodes of the climate system are captured, leading to better long-term estimates from short, non-equilibrated time series. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Main Text (10 pages, 4 figures) plus Supporting Information (36 pages, 18 figures, 1 table)

Journal ref: Geophysical Research Letters 48 (2021) e2020GL091090

A Delay Equation Model for the Atlantic Multidecadal Oscillation

Authors: Swinda K. J. Falkena, Courtney Quinn, Jan Sieber, Henk A. Dijkstra

Abstract: A new technique to derive delay models from systems of partial differential equations, based on the Mori-Zwanzig formalism, is used to derive a delay difference equation model for the Atlantic Multidecadal Oscillation. The Mori-Zwanzig formalism gives a rewriting of the original system of equations which contains a memory term. This memory term can be related to a delay term in a resulting delay e… ▽ More A new technique to derive delay models from systems of partial differential equations, based on the Mori-Zwanzig formalism, is used to derive a delay difference equation model for the Atlantic Multidecadal Oscillation. The Mori-Zwanzig formalism gives a rewriting of the original system of equations which contains a memory term. This memory term can be related to a delay term in a resulting delay equation. Here the technique is applied to an idealized, but spatially extended, model of the Atlantic Multidecadal Oscillation. The resulting delay difference model is of a different type than the delay differential model which has been used to describe the El Niño- Southern Oscillation. In addition to this model, which can also be obtained by integration along characteristics, error terms for a smoothing approximation of the model have been derived from the Mori-Zwanzig formalism. Our new method of deriving delay models from spatially extended models has a large potential to use delay models to study a range of climate variability phenomena. △ Less

Submitted 6 January, 2021; v1 submitted 17 August, 2020; originally announced August 2020.

Comments: Submitted to Proceedings of the Royal Society A - revised version

Moored observations of turbulent mixing events in deep Lake Garda (I)

Authors: Hans van Haren, Sebastiano Piccolroaz, Marina Amadori, Marco Toffolon, Henk A. Dijkstra

Abstract: Deep water circulation and mixing processes in deep lakes are largely unknown, although they are responsible for the transport of matter, nutrients and pollutants. Such a lack of knowledge cannot be reliably provided by numerical hydrodynamic modelling studies because detailed observations are typically not available to validate them. To overcome some of these deficiencies, a dedicated yearlong mo… ▽ More Deep water circulation and mixing processes in deep lakes are largely unknown, although they are responsible for the transport of matter, nutrients and pollutants. Such a lack of knowledge cannot be reliably provided by numerical hydrodynamic modelling studies because detailed observations are typically not available to validate them. To overcome some of these deficiencies, a dedicated yearlong mooring comprising 100 high-resolution temperature sensors and a single current meter were located in the deeper half of the 344 m deepest point of the subalpine Lake Garda (Italy). The observations show peaks and calms of turbulent exchange, besides ubiquitous internal wave activity. In late winter, northerly winds activate episodic deep convective overturning, the dense water being subsequently advected along the lake-floor. Besides deep convection, such winds also set-up seiches and inertial waves that are associated with about 100 times larger turbulence dissipation rates than that by semidiurnal internal wave breaking observed in summer. In the lower 60 m above the lake-floor however, the average turbulence dissipation rate is approximately constant in value year-around, being about 10 times larger than open-ocean values, except during deep convection episodes. △ Less

Submitted 4 August, 2020; originally announced August 2020.

Comments: 42 pages, 10 figures

Mid-infrared emission and absorption in strained and relaxed direct bandgap GeSn semiconductors

Authors: Simone Assali, Alain Dijkstra, Anis Attiaoui, Étienne Bouthillier, Jos E. M. Haverkort, Oussama Moutanabbir

Abstract: By independently engineering strain and composition, this work demonstrates and investigates direct band gap emission in the mid-infrared range from GeSn layers grown on silicon. We extend the room-temperature emission wavelength above ~4.0 μm upon post-growth strain relaxation in layers with uniform Sn content of 17 at.%. The fundamental mechanisms governing the optical emission are discussed bas… ▽ More By independently engineering strain and composition, this work demonstrates and investigates direct band gap emission in the mid-infrared range from GeSn layers grown on silicon. We extend the room-temperature emission wavelength above ~4.0 μm upon post-growth strain relaxation in layers with uniform Sn content of 17 at.%. The fundamental mechanisms governing the optical emission are discussed based on temperature-dependent photoluminescence, absorption measurements, and theoretical simulations. Regardless of strain and composition, these analyses confirm that single-peak emission is always observed in the probed temperature range of 4-300 K, ruling out defect- and impurity-related emission. Moreover, carrier losses into thermally-activated non-radiative recombination channels are found to be greatly minimized as a result of strain relaxation. Absorption measurements validate the direct band gap absorption in strained and relaxed samples at energies closely matching photoluminescence data. These results highlight the strong potential of GeSn semiconductors as versatile building blocks for scalable, compact, and silicon-compatible mid-infrared photonics and quantum opto-electronics. △ Less

Submitted 22 December, 2020; v1 submitted 28 April, 2020; originally announced April 2020.

Journal ref: Phys. Rev. Applied 15, 024031 (2021)

Resolution dependency of sinking Lagrangian particles in ocean general circulation models

Authors: Peter D. Nooteboom, Philippe Delandmeter, Erik van Sebille, Peter K. Bijl, Henk A. Dijkstra, Anna S. von der Heydt

Abstract: Any type of non-buoyant material in the ocean is transported horizontally by currents during its sinking journey. This lateral transport can be far from negligible for small sinking velocities. To estimate its magnitude and direction, the material is often modelled as a set of Lagrangian particles advected by current velocities that are obtained from Ocean General Circulation Models (OGCMs). State… ▽ More Any type of non-buoyant material in the ocean is transported horizontally by currents during its sinking journey. This lateral transport can be far from negligible for small sinking velocities. To estimate its magnitude and direction, the material is often modelled as a set of Lagrangian particles advected by current velocities that are obtained from Ocean General Circulation Models (OGCMs). State-of-the-art OGCMs are strongly eddying, similar to the real ocean, providing results with a spatial resolution on the order of 10 km on a daily frequency. While the importance of eddies in OGCMs is well-appreciated in the physical oceanographic community, other marine research communities may not. To demonstrate how much the absence of mesoscale features in low-resolution models influences the Lagrangian particle transport, we simulate the transport of sinking Lagrangian particles using low- and high-resolution global OGCMs, and assess the lateral transport differences resulting from the difference in spatial and temporal model resolution. We find major differences between the transport in the non-eddying OGCM and in the eddying OGCM. Addition of stochastic noise to the particle trajectories in the non-eddying OGCM parameterises the effect of eddies well in some cases. The effect of a coarser temporal resolution (5-daily) is smaller compared to a coarser spatial resolution (0.1$^{\circ}$ versus 1$^{\circ}$ horizontally). We recommend to use sinking Lagrangian particles, representing e.g. marine snow, microplankton or sinking plastic, only with velocity fields from eddying OGCMs, requiring high-resolution models in e.g. paleoceanographic studies. To increase the accessibility of our particle trace simulations, we launch planktondrift.science.uu.nl, an online tool to reconstruct the surface origin of sedimentary particles in a specific location. △ Less

Submitted 14 September, 2020; v1 submitted 15 April, 2020; originally announced April 2020.

Convection under internal waves in a mountain lake

Authors: Hans van Haren, Henk A. Dijkstra

Abstract: Turbulent mixing processes in deep alpine Lake Garda (I) have not extensively been observed. Knowledge about drivers of turbulent fluxes are important for insights in the transport of matter, nutrients and pollutants, in the lake and in natural water bodies in general. In this paper, the occurrence of internal wave induced turbulent convection is addressed as opposed to the more common shear-induc… ▽ More Turbulent mixing processes in deep alpine Lake Garda (I) have not extensively been observed. Knowledge about drivers of turbulent fluxes are important for insights in the transport of matter, nutrients and pollutants, in the lake and in natural water bodies in general. In this paper, the occurrence of internal wave induced turbulent convection is addressed as opposed to the more common shear-induced turbulence in a density stratified environment. Observations are analyzed from a dedicated yearlong mooring holding 100 high-resolution temperature sensors at 1.5 m intervals under a single current meter in the deeper half of the 344 m deep lake-center. Episodically, the weakly density stratified waters in the lower 50 m above the lake floor show spectral slope and coherence evidence of short-term (15 to 30 minutes) convective motions under internal waves that are supported by the stronger stratified waters above. The near-homogeneous conditions are not attributable to frictional Ekman dynamics, but to large-scale internal wave crests. △ Less

Submitted 22 April, 2021; v1 submitted 7 February, 2020; originally announced February 2020.

Comments: 19 pages, 4 figures

Journal ref: environmental fluid dynamics, 2021, 21, 305-316

Ruelle-Pollicott Resonances of Stochastic Systems in Reduced State Space. Part I: Theory

Authors: Mickaël Chekroun, Alexis Tantet, Henk Anton Dijkstra, J. David Neelin

Abstract: A theory of Ruelle-Pollicott (RP) resonances for stochastic differential systems is presented. These resonances are defined as the eigenvalues of the generator (Kolmogorov operator) of a given stochastic system. By relying on the theory of Markov semigroups, decomposition formulas of correlation functions and power spectral densities (PSDs) in terms of RP resonances are then derived. These formula… ▽ More A theory of Ruelle-Pollicott (RP) resonances for stochastic differential systems is presented. These resonances are defined as the eigenvalues of the generator (Kolmogorov operator) of a given stochastic system. By relying on the theory of Markov semigroups, decomposition formulas of correlation functions and power spectral densities (PSDs) in terms of RP resonances are then derived. These formulas describe, for a broad class of stochastic differential equations (SDEs), how the RP resonances characterize the decay of correlations as well as the signal's oscillatory components manifested by peaks in the PSD.It is then shown that a notion reduced RP resonances can be rigorously defined, as soon as the dynamics is partially observed within a reduced state space V . These reduced resonances are obtained from the spectral elements of reduced Markov operators acting on functions of the state space V , and can be estimated from series. They inform us about the spectral elements of some coarse-grained version of the SDE generator. When the time-lag at which the transitions are collected from partial observations in V , is either sufficiently small or large, it is shown that the reduced RP resonances approximate the (weak) RP resonances of the generator of the conditional expectation in V , i.e. the optimal reduced system in V obtained by averaging out the contribution of the unobserved variables. The approach is illustrated on a stochastic slow-fast system for which it is shown that the reduced RP resonances allow for a good reconstruction of the correlation functions and PSDs, even when the time-scale separation is weak.The companions articles, Part II and Part III, deal with further practical aspects of the theory presented in this contribution. One important byproduct consists of the diagnosis usefulness of stochastic dynamics that RP resonances provide. This is illustrated in the case of a stochastic Hopf bifurcation in Part II. There, it is shown that such a bifurcation has a clear manifestation in terms of a geometric organization of the RP resonances along discrete parabolas in the left half plane. Such geometric features formed by (reduced) RP resonances are extractable from time series and allow thus for providing an unambiguous "signature" of nonlinear oscillations embedded within a stochastic background. By relying then on the theory of reduced RP resonances presented in this contribution, Part III addresses the question of detection and characterization of such oscillations in a high-dimensional stochastic system, namely the Cane-Zebiak model of El Ni{ñ}o-Southern Oscillation subject to noise modeling fast atmospheric fluctuations. △ Less

Submitted 20 March, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Journal ref: Journal of Statistical Physics, Springer Verlag

Direct Bandgap Emission from Hexagonal Ge and SiGe Alloys

Authors: E. M. T. Fadaly, A. Dijkstra, J. R. Suckert, D. Ziss, M. A. J. v. Tilburg, C. Mao, Y. Ren, V. T. v. Lange, S. Kölling, M. A. Verheijen, D. Busse, C. Rödl, J. Furthmüller, F. Bechstedt, J. Stangl, J. J. Finley, S. Botti, J. E. M. Haverkort, E. P. A. M. Bakkers

Abstract: Silicon crystallized in the usual cubic (diamond) lattice structure has dominated the electronics industry for more than half a century. However, cubic silicon (Si), germanium (Ge) and SiGe-alloys are all indirect bandgap semiconductors that cannot emit light efficiently. Accordingly, achieving efficient light emission from group-IV materials has been a holy grail in silicon technology for decades… ▽ More Silicon crystallized in the usual cubic (diamond) lattice structure has dominated the electronics industry for more than half a century. However, cubic silicon (Si), germanium (Ge) and SiGe-alloys are all indirect bandgap semiconductors that cannot emit light efficiently. Accordingly, achieving efficient light emission from group-IV materials has been a holy grail in silicon technology for decades and, despite tremendous efforts, it has remained elusive. Here, we demonstrate efficient light emission from direct bandgap hexagonal Ge and SiGe alloys. We measure a subnanosecond, temperature-insensitive radiative recombination lifetime and observe a similar emission yield to direct bandgap III-V semiconductors. Moreover, we demonstrate how by controlling the composition of the hexagonal SiGe alloy, the emission wavelength can be continuously tuned in a broad range, while preserving a direct bandgap. Our experimental findings are shown to be in excellent quantitative agreement with the ab initio theory. Hexagonal SiGe embodies an ideal material system to fully unite electronic and optoelectronic functionalities on a single chip, opening the way towards novel device concepts and information processing technologies. △ Less

Submitted 2 November, 2019; originally announced November 2019.

Comments: 25 pages,5 main figures, 7 supplementary figures

Quantum dissipative systems beyond the standard harmonic model: features of linear absorption and dynamics

Authors: Luke D. Smith, Arend G. Dijkstra

Abstract: Current simulations of ultraviolet-visible absorption lineshapes, and dynamics of condensed phase systems, largely adopt a harmonic description to model vibrations. Often, this involves a model of displaced harmonic oscillators that have the same curvature. Although convenient, for many realistic molecular systems this approximation no longer suffices. We elucidate non-standard harmonic, and anhar… ▽ More Current simulations of ultraviolet-visible absorption lineshapes, and dynamics of condensed phase systems, largely adopt a harmonic description to model vibrations. Often, this involves a model of displaced harmonic oscillators that have the same curvature. Although convenient, for many realistic molecular systems this approximation no longer suffices. We elucidate non-standard harmonic, and anharmonic effects, on linear absorption and dynamics using a stochastic Schrödinger equation approach to account for the environment. Firstly, a harmonic oscillator model with ground and excited potentials that differ in curvature is utilised. Using this model, it is shown that curvature difference gives rise to an additional sub-structure in the vibronic progression of absorption spectra. This effect is explained, and subsequently quantified, via a derived expression for the Franck-Condon coefficients. Subsequently, anharmonic features in dissipative systems are studied, using a Morse potential, and parameters that correspond to the diatomic molecule $H_{2}$ for differing displacements and environment interaction. Lastly using a model potential, the population dynamics and absorption spectra for the stiff-stilbene photoswitch is presented and features are explained by a combination of curvature difference and anharmonicity in the form of potential energy barriers on the excited potential. △ Less

Submitted 27 August, 2019; originally announced August 2019.

Comments: 16 pages, 9 figures

Kinetic Control of Morphology and Composition in Ge/GeSn Core/Shell Nanowires

Authors: Simone Assali, Roberto Bergamaschini, Emilio Scalise, Marcel A. Verheijen, Marco Albani, Alain Dijkstra, Ang Li, Sebastian Koelling, Erik P. A. M. Bakkers, Francesco Montalenti, Leo Miglio

Abstract: The growth of Sn-rich group-IV semiconductors at the nanoscale provides new paths for understanding the fundamental properties of metastable GeSn alloys. Here, we demonstrate the effect of the growth conditions on the morphology and composition of Ge/GeSn core/shell nanowires by correlating the experimental observations with a theoretical interpretation based on a multi-scale approach. We show tha… ▽ More The growth of Sn-rich group-IV semiconductors at the nanoscale provides new paths for understanding the fundamental properties of metastable GeSn alloys. Here, we demonstrate the effect of the growth conditions on the morphology and composition of Ge/GeSn core/shell nanowires by correlating the experimental observations with a theoretical interpretation based on a multi-scale approach. We show that the cross-sectional morphology of Ge/GeSn core/shell nanowires changes from hexagonal to dodecagonal upon increasing the supply of the Sn precursor. This transformation strongly influences the Sn distribution as a higher Sn content is measured under the {112} growth front. Ab-initio DFT calculations provide an atomic-scale explanation by showing that Sn incorporation is favored at the {112} surfaces, where the Ge bonds are tensile-strained. A phase-field continuum model was developed to reproduce the morphological transformation and the Sn distribution within the wire, shedding light on the complex growth mechanism and unveiling the relation between segregation and faceting. The tunability of the photoluminescence emission with the change in composition and morphology of the GeSn shell highlights the potential of the core/shell nanowire system for opto-electronic devices operating at mid-infrared wavelengths. △ Less

Submitted 16 January, 2020; v1 submitted 27 June, 2019; originally announced June 2019.

Efficient long-distance energy transport in molecular systems through adiabatic passage

Authors: Arend G. Dijkstra, Almut Beige

Abstract: The efficiencies of light-harvesting complexes in biological systems can be much higher than the current efficiencies of artificial solar cells. In this paper, we therefore propose and analyse an energy transport mechanism which employs adiabatic passages between the states of an artificially designed antenna molecular system to significantly enhance the conversion of incoming light into internal… ▽ More The efficiencies of light-harvesting complexes in biological systems can be much higher than the current efficiencies of artificial solar cells. In this paper, we therefore propose and analyse an energy transport mechanism which employs adiabatic passages between the states of an artificially designed antenna molecular system to significantly enhance the conversion of incoming light into internal energy. It is shown that the proposed transport mechanism is relatively robust against spontaneous emission and de-phasing, while also being able to take advantage of collective effects. Our aim is to provide new insight into the energy transport in molecular complexes and to improve the design of solar cells. △ Less

Submitted 9 March, 2020; v1 submitted 14 June, 2019; originally announced June 2019.

Comments: 8 pages, 5 figures, typos in Section III.B corrected

Journal ref: J. Chem. Phys. 151, 034114 (2019)

Modeling and analyzing a photo-driven molecular motor system: Ratchet dynamics and non-linear optical spectra

Authors: Tatsushi Ikeda, Arend G. Dijkstra, Yoshitaka Tanimura

Abstract: A light-driven molecular motor system is investigated using a multi-state Brownian ratchet model described by a single effective coordinate with multiple electronic states in a dissipative environment. The rotational motion of the motor system is investigated on the basis of wavepacket dynamics. A current determined from the interplay between a fast photochemical isomerization (photoisomerization)… ▽ More A light-driven molecular motor system is investigated using a multi-state Brownian ratchet model described by a single effective coordinate with multiple electronic states in a dissipative environment. The rotational motion of the motor system is investigated on the basis of wavepacket dynamics. A current determined from the interplay between a fast photochemical isomerization (photoisomerization) process triggered by pulses and a slow thermal isomerization (thermalization) process arising from an overdamped environment is numerically evaluated. For this purpose, we employ the multi-state low-temperature quantum Smoluchowski equations that allow us to simulate the fast quantum electronic dynamics in the overdamped environment. We analyze the motor efficiency by numerically integrating the equations of motion for a rotator system driven by repeatedly impulsive excitations. When the timescales of the pulse repetition, photoisomerization, and thermalization processes are separated, the average rotational speed of the motor is determined by the timescale of thermalization. In this regime, the average rotational current can be described by a simple equation derived from a rate equation for the thermalization process. When laser pulses are applied repeatedly and the timescales of the photoisomerization and pulse repetition are close, the details of the photoisomerization process become important to analyze the entire rotational process. We examine the possibility of observing the photoisomerization and the thermalization processes associated with stationary rotating dynamics of the motor system by spectroscopic means, e.g. pump-probe, transient absorption, and two-dimensional electronic spectroscopy techniques. △ Less

Submitted 21 March, 2019; v1 submitted 19 February, 2019; originally announced February 2019.

Comments: 18 pages, 13 figures

Journal ref: J. Chem. Phys. 150, 114103 (2019)

Derivation of Delay Equation Climate Models Using the Mori-Zwanzig Formalism

Authors: Swinda K. J. Falkena, Courtney Quinn, Jan Sieber, Jason Frank, Henk A. Dijkstra

Abstract: Models incorporating delay have been frequently used to understand climate variability phenomena, but often the delay is introduced through an ad-hoc physical reasoning, such as the propagation time of waves. In this paper, the Mori-Zwanzig formalism is introduced as a way to systematically derive delay models from systems of partial differential equations and hence provides a better justification… ▽ More Models incorporating delay have been frequently used to understand climate variability phenomena, but often the delay is introduced through an ad-hoc physical reasoning, such as the propagation time of waves. In this paper, the Mori-Zwanzig formalism is introduced as a way to systematically derive delay models from systems of partial differential equations and hence provides a better justification for using these delay-type models. The Mori-Zwanzig technique gives a formal rewriting of the system using a projection onto a set of resolved variables, where the rewritten system contains a memory term. The computation of this memory term requires solving the orthogonal dynamics equation, which represents the unresolved dynamics. For nonlinear systems, it is often not possible to obtain an analytical solution to the orthogonal dynamics and an approximate solution needs to be found. Here, we demonstrate the Mori-Zwanzig technique for a two-strip model of the El Nino Southern Oscillation (ENSO) and explore methods to solve the orthogonal dynamics. The resulting nonlinear delay model contains an additional term compared to previously proposed ad-hoc conceptual models. This new term leads to a larger ENSO period, which is closer to that seen in observations. △ Less

Submitted 17 May, 2019; v1 submitted 8 February, 2019; originally announced February 2019.

Comments: Submitted to Proceedings of the Royal Society A, 25 pages, 10 figures

Controlling a Quantum System via its Boundary Conditions

Authors: Christian Duffin, Arend G. Dijkstra

Abstract: We numerically study a particle in a box with moving walls. In the case where the walls are oscillating sinusoidally with small amplitude, we show that states up to the fourth state can be populated with more than 80 percent population, while higher-lying states can also be selectively excited. This work introduces a way of controlling quantum systems which does not rely on (dipole) selection rule… ▽ More We numerically study a particle in a box with moving walls. In the case where the walls are oscillating sinusoidally with small amplitude, we show that states up to the fourth state can be populated with more than 80 percent population, while higher-lying states can also be selectively excited. This work introduces a way of controlling quantum systems which does not rely on (dipole) selection rules. △ Less

Submitted 1 October, 2018; originally announced October 2018.

Numerical Bifurcation Analysis of Marine Ice Sheet Models

Authors: Thomas E. Mulder, Henk A. Dijkstra, Fred W. Wubs

Abstract: The climate variability associated with the Pleistocene Ice Ages is one of the most fascinating puzzles in the Earth Sciences still awaiting a satisfactory explanation. In particular, the explanation of the dominant 100 kyr period of the glacial cycles over the last million years is a long-standing problem. Based on bifurcation analyses of low-order models, many theories have been suggested to exp… ▽ More The climate variability associated with the Pleistocene Ice Ages is one of the most fascinating puzzles in the Earth Sciences still awaiting a satisfactory explanation. In particular, the explanation of the dominant 100 kyr period of the glacial cycles over the last million years is a long-standing problem. Based on bifurcation analyses of low-order models, many theories have been suggested to explain these cycles and their frequency. The new aspect in this contribution is that, for the first time, numerical bifurcation analysis is applied to a two-dimensional marine ice sheet model with a dynamic grounding line. In this model, we find Hopf bifurcations with an oscillation period of about 100 kyr which may be relevant to glacial cycles. △ Less

Submitted 3 May, 2018; originally announced May 2018.

Comments: in Computational Modelling of Bifurcations and Instabilities in Fluid Dynamics, ed. Alexander Gelfgat (Springer, 2018)

Using Network Theory and Machine Learning to predict El Niño

Authors: Peter D. Nooteboom, Qing Yi Feng, Cristóbal López, Emilio Hernández-García, Henk A. Dijkstra

Abstract: The skill of current predictions of the warm phase of the El Niño Southern Oscillation (ENSO) reduces significantly beyond a lag of six months. In this paper, we aim to increase this prediction skill at lags up to one year. The new method to do so combines a classical Autoregressive Integrated Moving Average technique with a modern machine learning approach (through an Artificial Neural Network).… ▽ More The skill of current predictions of the warm phase of the El Niño Southern Oscillation (ENSO) reduces significantly beyond a lag of six months. In this paper, we aim to increase this prediction skill at lags up to one year. The new method to do so combines a classical Autoregressive Integrated Moving Average technique with a modern machine learning approach (through an Artificial Neural Network). The attributes in such a neural network are derived from topological properties of Climate Networks and are tested on both a Zebiak-Cane-type model and observations. For predictions up to six months ahead, the results of the hybrid model give a better skill than the CFSv2 ensemble prediction by the National Centers for Environmental Prediction (NCEP). Moreover, results for a twelve month lead time prediction have a similar skill as the shorter lead time predictions. △ Less

Submitted 10 August, 2018; v1 submitted 27 March, 2018; originally announced March 2018.

Comments: Earth Syst. Dynam., 2018

Journal ref: Earth System Dynamics 9, 969-983 (2018)

Resonances in a Chaotic Attractor Crisis of the Lorenz Flow

Authors: Alexis Tantet, Valerio Lucarini, Henk A. Dijkstra

Abstract: Local bifurcations of stationary points and limit cycles have successfully been characterized in terms of the critical exponents of these solutions. Lyapunov exponents and their associated covariant Lyapunov vectors have been proposed as tools for supporting the understanding of critical transitions in chaotic dynamical systems. However, it is in general not clear how the statistical properties of… ▽ More Local bifurcations of stationary points and limit cycles have successfully been characterized in terms of the critical exponents of these solutions. Lyapunov exponents and their associated covariant Lyapunov vectors have been proposed as tools for supporting the understanding of critical transitions in chaotic dynamical systems. However, it is in general not clear how the statistical properties of dynamical systems change across a boundary crisis during which a chaotic attractor collides with a saddle. This behavior is investigated here for a boundary crisis in the Lorenz flow, for which neither the Lyapunov exponents nor the covariant Lyapunov vectors provide a criterion for the crisis. Instead, the convergence of the time evolution of probability densities to the invariant measure, governed by the semigroup of transfer operators, is expected to slow down at the approach of the crisis. Such convergence is described by the eigenvalues of the generator of this semigroup, which can be divided into two families, referred to as the stable and unstable Ruelle-Pollicott resonances, respectively. The former describes the convergence of densities to the attractor (or escape from a repeller) and is estimated from many short time series sampling the phase space. The latter is responsible for the decay of correlations, or mixing, and can be estimated from a long times series, invoking ergodicity. It is found numerically for the Lorenz flow that the stable resonances do approach the imaginary axis during the crisis, as is indicative of the loss of global stability of the attractor. On the other hand, the unstable resonances, and a fortiori the decay of correlations, do not ag the proximity of the crisis, thus questioning the usual design of early warning indicators of boundary crises of chaotic attractors and the applicability of response theory close to such crises. △ Less

Submitted 23 May, 2017; originally announced May 2017.

Ruelle-Pollicott Resonances of Stochastic Systems in Reduced State Space. Part II: Stochastic Hopf Bifurcation

Authors: Alexis Tantet, Mickaël D. Chekroun, Henk A. Dijkstra, J. David Neelin

Abstract: The spectrum of the generator (Kolmogorov operator) of a diffusion process, referred to as the Ruelle-Pollicott (RP) spectrum, provides a detailed characterization of correlation functions and power spectra of stochastic systems via decomposition formulas in terms of RP resonances. Stochastic analysis techniques relying on the theory of Markov semigroups for the study of the RP spectrum and a rigo… ▽ More The spectrum of the generator (Kolmogorov operator) of a diffusion process, referred to as the Ruelle-Pollicott (RP) spectrum, provides a detailed characterization of correlation functions and power spectra of stochastic systems via decomposition formulas in terms of RP resonances. Stochastic analysis techniques relying on the theory of Markov semigroups for the study of the RP spectrum and a rigorous reduction method is presented in Part I. This framework is here applied to study a stochastic Hopf bifurcation in view of characterizing the statistical properties of nonlinear oscillators perturbed by noise, depending on their stability. In light of the Hörmander theorem, it is first shown that the geometry of the unperturbed limit cycle, in particular its isochrons, is essential to understand the effect of noise and the phenomenon of phase diffusion. In addition, it is shown that the spectrum has a spectral gap, even at the bifurcation point, and that correlations decay exponentially fast. Explicit small-noise expansions of the RP eigenvalues and eigenfunctions are then obtained, away from the bifurcation point, based on the knowledge of the linearized deterministic dynamics and the characteristics of the noise. These formulas allow one to understand how the interaction of the noise with the deterministic dynamics affect the decay of correlations. Numerical results complement the study of the RP spectrum at the bifurcation, revealing useful scaling laws. The analysis of the Markov semigroup for stochastic bifurcations is thus promising in providing a complementary approach to the more geometric random dynamical system approach. This approach is not limited to low-dimensional systems and the reduction method presented in part I is applied to a stochastic model relevant to climate dynamics in part III. △ Less

Submitted 9 March, 2020; v1 submitted 22 May, 2017; originally announced May 2017.

A Lazy Language Needs a Lazy Type System: Introducing Polymorphic Contexts

Authors: S. Doaitse Swierstra, Marcos Viera, Atze Dijkstra

Abstract: Most type systems that support polymorphic functions are based on a version of System-F. We argue that this limits useful programming paradigms for languages with lazy evaluation. We motivate an extension of System-F alleviating this limitation. First, using a sequence of examples, we show that for lazily evaluated languages current type systems may force one to write a program in an unnatural w… ▽ More Most type systems that support polymorphic functions are based on a version of System-F. We argue that this limits useful programming paradigms for languages with lazy evaluation. We motivate an extension of System-F alleviating this limitation. First, using a sequence of examples, we show that for lazily evaluated languages current type systems may force one to write a program in an unnatural way; we in particular argue that in such languages the relationship between polymorphic and existential types can be made more systematic by allowing to pass back (part of) an existential result of a function call as an argument to the the function call that produced that value. After presenting our extension to System-F we show how we can implement the strict-state thread monad $\mathrm{ST}$ by using a returned existential type in specialising the polymorphic function which returns that type. Currently this monad is built-in into the runtime system of GHC and as such has become part of the language. Our proposed language extension, i.e. the introduction of polymorphic contexts, reverses the relationship between the context of a function call and the called function with respect to where it is decided with which type to instantiate a type variable. △ Less

Submitted 14 December, 2016; originally announced December 2016.

Report number: UU-CS-2016-012 ACM Class: D.3.3; D.1.1

How two-dimensional brick layer J-aggregates differ from linear ones: excitonic properties and line broadening mechanisms

Authors: Arend G. Dijkstra, Hong-Guang Duan, Jasper Knoester, Keith A. Nelson, Jianshu Cao

Abstract: We study the excitonic coupling and homogeneous spectral line width of brick layer J-aggregate films. We begin by analysing the structural information revealed by the two-exciton states probed in two-dimensional spectra. Our first main result is that the relation between the excitonic couplings and the spectral shift in a two-dimensional structure is different (larger shift for the same nearest ne… ▽ More We study the excitonic coupling and homogeneous spectral line width of brick layer J-aggregate films. We begin by analysing the structural information revealed by the two-exciton states probed in two-dimensional spectra. Our first main result is that the relation between the excitonic couplings and the spectral shift in a two-dimensional structure is different (larger shift for the same nearest neighbour coupling) from that in a one-dimensional structure, which leads to an estimation of dipolar coupling in two-dimensional lattices. We next investigate the mechanisms of homogeneous broadening - population relaxation and pure dephasing - and evaluate their relative importance in linear and two-dimensional aggregates. Our second main result is that pure dephasing dominates the line width in two-dimensional systems up to a crossover temperature, which explains the linear temperature dependence of the homogeneous line width. This is directly related to the decreased density of states at the band edge when compared with linear aggregates, thus reducing the contribution of population relaxation to dephasing. Pump-probe experiments are suggested to directly measure the lifetime of the bright state and can therefore support the proposed model. △ Less

Submitted 16 March, 2016; originally announced March 2016.

Comments: accepted for publication in J. Chem. Phys

Energy transfer in structured and unstructured environments: master equations beyond the Born-Markov approximations

Authors: Jake Iles-Smith, Arend G. Dijkstra, Neill Lambert, Ahsan Nazir

Abstract: We explore excitonic energy transfer dynamics in a molecular dimer system coupled to both structured and unstructured oscillator environments. By extending the reaction coordinate master equation technique developed in [J. Iles-Smith, N. Lambert, and A. Nazir, Phys. Rev. A 90, 032114 (2014)], we go beyond the commonly used Born-Markov approximations to incorporate system-environment correlations a… ▽ More We explore excitonic energy transfer dynamics in a molecular dimer system coupled to both structured and unstructured oscillator environments. By extending the reaction coordinate master equation technique developed in [J. Iles-Smith, N. Lambert, and A. Nazir, Phys. Rev. A 90, 032114 (2014)], we go beyond the commonly used Born-Markov approximations to incorporate system-environment correlations and the resultant non-Markovian dynamical effects. We obtain energy transfer dynamics for both underdamped and overdamped oscillator environments that are in perfect agreement with the numerical hierarchical equations of motion over a wide range of parameters. Furthermore, we show that the Zusman equations, which may be obtained in a semiclassical limit of the reaction coordinate model, are often incapable of describing the correct dynamical behaviour. This demonstrates the necessity of properly accounting for quantum correlations generated between the system and its environment when the Born-Markov approximations no longer hold. Finally, we apply the reaction coordinate formalism to the case of a structured environment comprising of both underdamped (i.e. sharply peaked) and overdamped (broad) components simultaneously. We find that though an enhancement of the dimer energy transfer rate can be obtained when compared to an unstructured environment, its magnitude is rather sensitive to both the dimer-peak resonance conditions and the relative strengths of the underdamped and overdamped contributions. △ Less

Submitted 4 April, 2016; v1 submitted 16 November, 2015; originally announced November 2015.

Comments: 14 pages, 11 figures. v2 - close to published version

Journal ref: J. Chem. Phys. 144, 044110 (2016)

Effects of Drake Passage on a strongly eddying global ocean

Authors: Jan P. Viebahn, Anna S. von der Heydt, Dewi Le Bars, Henk A. Dijkstra

Abstract: The climate impact of ocean gateway openings during the Eocene-Oligocene transition is still under debate. Previous model studies employed grid resolutions at which the impact of mesoscale eddies has to be parameterized. We present results of a state-of-the-art eddy-resolving global ocean model with a closed Drake Passage, and compare with results of the same model at non-eddying resolution. An an… ▽ More The climate impact of ocean gateway openings during the Eocene-Oligocene transition is still under debate. Previous model studies employed grid resolutions at which the impact of mesoscale eddies has to be parameterized. We present results of a state-of-the-art eddy-resolving global ocean model with a closed Drake Passage, and compare with results of the same model at non-eddying resolution. An analysis of the pathways of heat by decomposing the meridional heat transport into eddy, horizontal, and overturning circulation components indicates that the model behavior on the large scale is qualitatively similar at both resolutions. Closing Drake Passage induces (i) sea surface warming around Antarctica due to changes in the horizontal circulation of the Southern Ocean, (ii) the collapse of the overturning circulation related to North Atlantic Deep Water formation leading to surface cooling in the North Atlantic, (iii) significant equatorward eddy heat transport near Antarctica. However, quantitative details significantly depend on the chosen resolution. The warming around Antarctica is substantially larger for the non-eddying configuration (5.5{degree sign}C) than for the eddying configuration (2.5{degree sign}C). This is a consequence of the subpolar mean flow which partitions differently into gyres and circumpolar current at different resolutions. We conclude that for a deciphering of the different mechanisms active in Eocene-Oligocene climate change detailed analyses of the pathways of heat in the different climate subsystems are crucial in order to clearly identify the physical processes actually at work. △ Less

Submitted 14 October, 2015; originally announced October 2015.

Crisis of the Chaotic Attractor of a Climate Model: A Transfer Operator Approach

Authors: Alexis Tantet, Valerio Lucarini, Frank Lunkeit, Henk A. Dijkstra

Abstract: The destruction of a chaotic attractor leading to rough changes in the dynamics of a dynamical system is studied. Local bifurcations are characterised by a single or a pair of characteristic exponents crossing the imaginary axis. The approach of such bifurcations in the presence of noise can be inferred from the slowing down of the correlation decay. On the other hand, little is known about global… ▽ More The destruction of a chaotic attractor leading to rough changes in the dynamics of a dynamical system is studied. Local bifurcations are characterised by a single or a pair of characteristic exponents crossing the imaginary axis. The approach of such bifurcations in the presence of noise can be inferred from the slowing down of the correlation decay. On the other hand, little is known about global bifurcations involving high-dimensional attractors with positive Lyapunov exponents. The global stability of chaotic attractors may be characterised by the spectral properties of the Koopman or the transfer operators governing the evolution of statistical ensembles. It has recently been shown that a boundary crisis in the Lorenz flow coincides with the approach to the unit circle of the eigenvalues of these operators associated with motions about the attractor, the stable resonances. A second type of resonances, the unstable resonances, is responsible for the decay of correlations and mixing on the attractor. In the deterministic case, those cannot be expected to be affected by general boundary crises. Here, however, we give an example of chaotic system in which slowing down of the decay of correlations of some observables does occur at the approach of a boundary crisis. The system considered is a high-dimensional, chaotic climate model of physical relevance. Moreover, coarse-grained approximations of the transfer operators on a reduced space, constructed from a long time series of the system, give evidence that this behaviour is due to the approach of unstable resonances to the unit circle. That the unstable resonances are affected by the crisis can be physically understood from the fact that the process responsible for the instability, the ice-albedo feedback, is also active on the attractor. Implications regarding response theory and the design of early-warning signals are discussed. △ Less

Submitted 13 December, 2017; v1 submitted 8 July, 2015; originally announced July 2015.

Unified functional network and nonlinear time series analysis for complex systems science: The pyunicorn package

Authors: Jonathan F. Donges, Jobst Heitzig, Boyan Beronov, Marc Wiedermann, Jakob Runge, Qing Yi Feng, Liubov Tupikina, Veronika Stolbova, Reik V. Donner, Norbert Marwan, Henk A. Dijkstra, Jürgen Kurths

Abstract: We introduce the \texttt{pyunicorn} (Pythonic unified complex network and recurrence analysis toolbox) open source software package for applying and combining modern methods of data analysis and modeling from complex network theory and nonlinear time series analysis. \texttt{pyunicorn} is a fully object-oriented and easily parallelizable package written in the language Python. It allows for the co… ▽ More We introduce the \texttt{pyunicorn} (Pythonic unified complex network and recurrence analysis toolbox) open source software package for applying and combining modern methods of data analysis and modeling from complex network theory and nonlinear time series analysis. \texttt{pyunicorn} is a fully object-oriented and easily parallelizable package written in the language Python. It allows for the construction of functional networks such as climate networks in climatology or functional brain networks in neuroscience representing the structure of statistical interrelationships in large data sets of time series and, subsequently, investigating this structure using advanced methods of complex network theory such as measures and models for spatial networks, networks of interacting networks, node-weighted statistics or network surrogates. Additionally, \texttt{pyunicorn} provides insights into the nonlinear dynamics of complex systems as recorded in uni- and multivariate time series from a non-traditional perspective by means of recurrence quantification analysis (RQA), recurrence networks, visibility graphs and construction of surrogate time series. The range of possible applications of the library is outlined, drawing on several examples mainly from the field of climatology. △ Less

Submitted 6 April, 2016; v1 submitted 2 July, 2015; originally announced July 2015.

Comments: 28 pages, 17 figures

Journal ref: Chaos 25, 113101 (2015)