Source author record

Henk A. Dijkstra

Henk A. Dijkstra appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2021arXiv

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

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.

preprint2021arXiv

Quantification and interpretation of the climate variability record

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.

preprint2020arXiv

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

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.

preprint2020arXiv

Multivariate Estimations of Equilibrium Climate Sensitivity from Short Transient Warming Simulations

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.

preprint2020arXiv

Resolution dependency of sinking Lagrangian particles in ocean general circulation models

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.

preprint2020arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

An early warning indicator for atmospheric blocking events using transfer operators

The existence of persistent midlatitude atmospheric flow regimes with time-scales larger than 5-10 days and indications of preferred transitions between them motivates to develop early warning indicators for such regime transitions. In this paper, we use a hemispheric barotropic model together with estimates of transfer operators on a reduced phase space to develop an early warning indicator of the zonal to blocked flow transition in this model. It is shown that, the spectrum of the transfer operators can be used to study the slow dynamics of the flow as well as the non-Markovian character of the reduction. The slowest motions are thereby found to have time scales of three to six weeks and to be associated with meta-stable regimes (and their transitions) which can be detected as almost-invariant sets of the transfer operator. From the energy budget of the model, we are able to explain the meta-stability of the regimes and the existence of preferred transition paths. Even though the model is highly simplified, the skill of the early warning indicator is promising, suggesting that the transfer operator approach can be used in parallel to an operational deterministic model for stochastic prediction or to assess forecast uncertainty.

preprint2015arXiv

Effects of Drake Passage on a strongly eddying global ocean

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.

preprint2014arXiv

Deep ocean early warning signals of an Atlantic MOC collapse

The Atlantic Meridional Overturning Circulation (MOC) is a crucial part of the climate system because of its associated northward heat transport. The present-day MOC is sensitive to freshwater anomalies and may collapse to a state with a strongly reduced northward heat transport. A future collapse of the Atlantic MOC has been identified as one of the most dangerous tipping points in the climate system. It is therefore crucial to develop early warning indicators for such a potential collapse based on relatively short time series. So far, attempts to use indicators based on critical slowdown have been marginally successful. Based on complex climate network reconstruction, we here present a promising new indicator for the MOC collapse that efficiently monitors spatial changes in deep ocean circulation. Through our analysis of the performance of this indicator we formulate optimal locations of measurement of the MOC to provide early warning signals of a collapse. Our results imply that an increase in spatial resolution of the Atlantic MOC observations (i.e., at more sections) can improve early detection, because the spatial coherence in the deep ocean arising near the transition is better captured.

preprint2014arXiv

On the state dependency of fast feedback processes in (palaeo) climate sensitivity

Palaeo data have been frequently used to determine the equilibrium (Charney) climate sensitivity $S^a$, and - if slow feedback processes (e.g. land ice-albedo) are adequately taken into account - they indicate a similar range as estimates based on instrumental data and climate model results. Most studies implicitly assume the (fast) feedback processes to be independent of the background climate state, e.g., equally strong during warm and cold periods. Here we assess the dependency of the fast feedback processes on the background climate state using data of the last 800 kyr and a conceptual climate model for interpretation. Applying a new method to account for background state dependency, we find $S^a=0.61\pm0.06$ K(Wm$^{-2}$)$^{-1}$ using the latest LGM temperature reconstruction and significantly lower climate sensitivity during glacial climates. Due to uncertainties in reconstructing the LGM temperature anomaly, $S^a$ is estimated in the range $S^a=0.55-0.95$ K(Wm$^{-2}$)$^{-1}$.

preprint2012arXiv

A global hybrid coupled model based on Atmosphere-SST feedbacks

A global hybrid coupled model is developed, with the aim of studying the effects of ocean-atmosphere feedbacks on the stability of the Atlantic meridional overturning circulation. The model includes a global ocean general circulation model and a statistical atmosphere model. The statistical atmosphere model is based on linear regressions of data from a fully coupled climate model on sea surface temperature both locally and hemispherically averaged, being the footprint of Atlantic meridional overturning variability. It provides dynamic boundary conditions to the ocean model for heat, freshwater and wind-stress. A basic but consistent representation of ocean-atmosphere feedbacks is captured in the hybrid coupled model and it is more than ten times faster than the fully coupled climate model. The hybrid coupled model reaches a steady state with a climate close to the one of the fully coupled climate model, and the two models also have a similar response (collapse) of the Atlantic meridional overturning circulation to a freshwater hosing applied in the northern North Atlantic.

preprint2012arXiv

High-Performance Distributed Multi-Model / Multi-Kernel Simulations: A Case-Study in Jungle Computing

High-performance scientific applications require more and more compute power. The concurrent use of multiple distributed compute resources is vital for making scientific progress. The resulting distributed system, a so-called Jungle Computing System, is both highly heterogeneous and hierarchical, potentially consisting of grids, clouds, stand-alone machines, clusters, desktop grids, mobile devices, and supercomputers, possibly with accelerators such as GPUs. One striking example of applications that can benefit greatly of Jungle Computing Systems are Multi-Model / Multi-Kernel simulations. In these simulations, multiple models, possibly implemented using different techniques and programming models, are coupled into a single simulation of a physical system. Examples include the domain of computational astrophysics and climate modeling. In this paper we investigate the use of Jungle Computing Systems for such Multi-Model / Multi-Kernel simulations. We make use of the software developed in the Ibis project, which addresses many of the problems faced when running applications on Jungle Computing Systems. We create a prototype Jungle-aware version of AMUSE, an astrophysical simulation framework. We show preliminary experiments with the resulting system, using clusters, grids, stand-alone machines, and GPUs.

preprint2012arXiv

Meridional overturning circulation: stability and ocean feedbacks in a box model

A box model of the inter-hemispheric Atlantic meridional overturning circulation is developed, including a variable pycnocline depth for the tropical and subtropical regions. The circulation is forced by winds over a periodic channel in the south and by freshwater forcing at the surface. The model is aimed at investigating the ocean feedbacks related to perturbations in freshwater forcing from the atmosphere, and to changes in freshwater transport in the ocean. These feedbacks are closely connected with the stability properties of the meridional overturning circulation, in particular in response to freshwater perturbations.

preprint2012arXiv

Sensitivity of the Atlantic Meridional Overturning Circulation to South Atlantic freshwater anomalies

The sensitivity of the Atlantic Meridional Overturning Circulation (AMOC) to changes in basin integrated net evaporation is highly dependent on the zonal salinity contrast at the southern border of the Atlantic. Biases in the freshwater budget strongly affect the stability of the AMOC in numerical models. The impact of these biases is investigated, by adding local anomaly patterns in the South Atlantic to the freshwater fluxes at the surface. These anomalies impact the freshwater and salt transport by the different components of the ocean circulation, in particular the basin--scale salt--advection feedback, completely changing the response of the AMOC to arbitrary perturbations. It is found that an appropriate dipole anomaly pattern at the southern border of the Atlantic Ocean can collapse the AMOC entirely even without a further hosing. The results suggest a new view on the stability of the AMOC, controlled by processes in the South Atlantic.