Global reorganization of atmospheric circulation during Dansgaard-Oeschger cycles.

Ice core records from Greenland provide evidence for multiple abrupt cold-warm-cold events recurring at millennial time scales during the last glacial interval. Although climate variations resembling Dansgaard-Oeschger (DO) oscillations have been identified in climate archives across the globe, our understanding of the climate and ecosystem impacts of the Greenland warming events in lower latitudes remains incomplete. Here, we investigate the influence of DO-cold-to-warm transitions on the global atmospheric circulation pattern. We comprehensively analyze δ18O changes during DO transitions in a globally distributed dataset of speleothems and set those in context with simulations of a comprehensive high-resolution climate model featuring internal millennial-scale variations of similar magnitude. Across the globe, speleothem δ18O signals and model results indicate consistent large-scale changes in precipitation amount, moisture source, or seasonality of precipitation associated with the DO transitions, in agreement with northward shifts of the Hadley circulation. Furthermore, we identify a decreasing trend in the amplitude of DO transitions with increasing distances from the North Atlantic region. This provides quantitative observational evidence for previous suggestions of the North Atlantic region being the focal point for these archetypes of past abrupt climate changes.

Climate change-induced population pressure drives high rates of lethal violence in the Prehispanic central Andes.

Understanding the influence of climate change and population pressure on human conflict remains a critically important topic in the social sciences. Long-term records that evaluate these dynamics across multiple centuries and outside the range of modern climatic variation are especially capable of elucidating the relative effect of­and the interaction between­climate and demography. This is crucial given that climate change may structure population growth and carrying capacity, while both climate and population influence per capita resource availability. This study couples paleoclimatic and demographic data with osteological evaluations of lethal trauma from 149 directly accelerator mass spectrometry 14C-dated individuals from the Nasca highland region of Peru. Multiple local and supraregional precipitation proxies are combined with a summed probability distribution of 149 14C dates to estimate population dynamics during a 700-y study window. Counter to previous findings, our analysis reveals a precipitous increase in violent deaths associated with a period of productive and stable climate, but volatile population dynamics. We conclude that favorable local climate conditions fostered population growth that put pressure on the marginal and highly circumscribed resource base, resulting in violent resource competition that manifested in over 450 y of internecine warfare. These findings help support a general theory of intergroup violence, indicating that relative resource scarcity­whether driven by reduced resource abundance or increased competition­can lead to violence in subsistence societies when the outcome is lower per capita resource availability.

Modeling a neurological disorder as the result of an operator acting on the brain: A first sketch based on network channel modeling.

The brain is a complex network, and diseases can alter its structures and connections between regions. Therefore, we can try to formalize the action of diseases by using operators acting on the brain network. Here, we propose a conceptual model of the brain, seen as a multilayer network, whose intra-lobe interactions are formalized as the diagonal blocks of an adjacency matrix. We propose a general and abstract definition of disease as an operator altering the weights of the connections between neural agglomerates, that is, the elements of the brain matrix. As models, we consider examples from three neurological disorders: epilepsy, Alzheimer-Perusini's disease, and schizophrenia. The alteration of neural connections can be seen as alterations of communication pathways, and thus, they can be described with a new channel model.

Resilience of transportation infrastructure networks to road failures.

Anthropogenic climate change drives extreme weather events, leading to significant consequences for both society and the environment. This includes damage to road infrastructure, causing disruptions in transportation, obstructing access to emergency services, and hindering humanitarian organizations after natural disasters. In this study, we develop a novel method for analyzing the impacts of natural hazards on transportation networks rooted in the gravity model of travel, offering a fresh perspective to assess the repercussions of natural hazards on transportation network stability. Applying this approach to the Ahr valley flood of 2021, we discovered that the destruction of bridges and roads caused major bottlenecks, affecting areas considerably distant from the flood's epicenter. Furthermore, the flood-induced damage to the infrastructure also increased the response time of emergency vehicles, severely impeding the accessibility of emergency services. Our findings highlight the need for targeted road repair and reinforcement, with a focus on maintaining traffic flow for emergency responses. This research provides a new perspective that can aid in prioritizing transportation network resilience measures to reduce the economic and social costs of future extreme weather events.

Model adaptive phase space reconstruction.

Phase space reconstruction (PSR) methods allow for the analysis of low-dimensional data with methods from dynamical systems theory, but their application to prediction models, such as those from machine learning (ML), is limited. Therefore, we here present a model adaptive phase space reconstruction (MAPSR) method that unifies the process of PSR with the modeling of the dynamical system. MAPSR is a differentiable PSR based on time-delay embedding and enables ML methods for modeling. The quality of the reconstruction is evaluated by the prediction loss. The discrete-time signal is converted into a continuous-time signal to achieve a loss function, which is differentiable with respect to the embedding delays. The delay vector, which stores all potential embedding delays, is updated along with the trainable parameters of the model to minimize prediction loss. Thus, MAPSR does not rely on any threshold or statistical criterion for determining the dimension and the set of delay values for the embedding process. We apply the MAPSR method to uni- and multivariate time series stemming from chaotic dynamical systems and a turbulent combustor. We find that for the Lorenz system, the model trained with the MAPSR method is able to predict chaotic time series for nearly seven to eight Lyapunov time scales, which is found to be much better compared to other PSR methods [AMI-FNN (average mutual information-false nearest neighbor) and PECUZAL (Pecora-Uzal) methods]. For the univariate time series from the turbulent combustor, the long-term cumulative prediction error of the MAPSR method for the regime of chaos stays between other methods, and for the regime of intermittency, MAPSR outperforms other PSR methods.

Power spectral estimate for discrete data.

The identification of cycles in periodic signals is a ubiquitous problem in time series analysis. Many real-world datasets only record a signal as a series of discrete events or symbols. In some cases, only a sequence of (non-equidistant) times can be assessed. Many of these signals are furthermore corrupted by noise and offer a limited number of samples, e.g., cardiac signals, astronomical light curves, stock market data, or extreme weather events. We propose a novel method that provides a power spectral estimate for discrete data. The edit distance is a distance measure that allows us to quantify similarities between non-equidistant event sequences of unequal lengths. However, its potential to quantify the frequency content of discrete signals has so far remained unexplored. We define a measure of serial dependence based on the edit distance, which can be transformed into a power spectral estimate (EDSPEC), analogous to the Wiener-Khinchin theorem for continuous signals. The proposed method is applied to a variety of discrete paradigmatic signals representing random, correlated, chaotic, and periodic occurrences of events. It is effective at detecting periodic cycles even in the presence of noise and for short event series. Finally, we apply the EDSPEC method to a novel catalog of European atmospheric rivers (ARs). ARs are narrow filaments of extensive water vapor transport in the lower troposphere and can cause hazardous extreme precipitation events. Using the EDSPEC method, we conduct the first spectral analysis of European ARs, uncovering seasonal and multi-annual cycles along different spatial domains. The proposed method opens new research avenues in studying of periodic discrete signals in complex real-world systems.

Interpolation and sampling effects on recurrence quantification measures.

The recurrence plot and the recurrence quantification analysis (RQA) are well-established methods for the analysis of data from complex systems. They provide important insights into the nature of the dynamics, periodicity, regime changes, and many more. These methods are used in different fields of research, such as finance, engineering, life, and earth science. To use them, the data have usually to be uniformly sampled, posing difficulties in investigations that provide non-uniformly sampled data, as typical in medical data (e.g., heart-beat based measurements), paleoclimate archives (such as sediment cores or stalagmites), or astrophysics (supernova or pulsar observations). One frequently used solution is interpolation to generate uniform time series. However, this preprocessing step can introduce bias to the RQA measures, particularly those that rely on the diagonal or vertical line structure in the recurrence plot. Using prototypical model systems, we systematically analyze differences in the RQA measure average diagonal line length for data with different sampling and interpolation. For real data, we show that the course of this measure strongly depends on the choice of the sampling rate for interpolation. Furthermore, we suggest a correction scheme, which is capable of correcting the bias introduced by the prepossessing step if the interpolation ratio is an integer.

GPU-based, interactive exploration of large spatiotemporal climate networks.

This paper introduces the Graphics Processing Unit (GPU)-based tool Geo-Temporal eXplorer (GTX), integrating a set of highly interactive techniques for visual analytics of large geo-referenced complex networks from the climate research domain. The visual exploration of these networks faces a multitude of challenges related to the geo-reference and the size of these networks with up to several million edges and the manifold types of such networks. In this paper, solutions for the interactive visual analysis for several distinct types of large complex networks will be discussed, in particular, time-dependent, multi-scale, and multi-layered ensemble networks. Custom-tailored for climate researchers, the GTX tool supports heterogeneous tasks based on interactive, GPU-based solutions for on-the-fly large network data processing, analysis, and visualization. These solutions are illustrated for two use cases: multi-scale climatic process and climate infection risk networks. This tool helps one to reduce the complexity of the highly interrelated climate information and unveils hidden and temporal links in the climate system, not available using standard and linear tools (such as empirical orthogonal function analysis).

Study of interaction and complete merging of binary cyclones using complex networks.

Cyclones are among the most hazardous extreme weather events on Earth. In certain scenarios, two co-rotating cyclones in close proximity to one another can drift closer and completely merge into a single cyclonic system. Identifying the dynamic transitions during such an interaction period of binary cyclones and predicting the complete merger (CM) event are challenging for weather forecasters. In this work, we suggest an innovative approach to understand the evolving vortical interactions between the cyclones during two such CM events (Noru-Kulap and Seroja-Odette) using time-evolving induced velocity-based unweighted directed networks. We find that network-based indicators, namely, in-degree and out-degree, quantify the changes in the interaction between the two cyclones and are excellent candidates to classify the interaction stages before a CM. The network indicators also help to identify the dominant cyclone during the period of interaction and quantify the variation of the strength of the dominating and merged cyclones. Finally, we show that the network measures also provide an early indication of the CM event well before its occurrence.

Stickiness and recurrence plots: An entropy-based approach.

The stickiness effect is a fundamental feature of quasi-integrable Hamiltonian systems. We propose the use of an entropy-based measure of the recurrence plots (RPs), namely, the entropy of the distribution of the recurrence times (estimated from the RP), to characterize the dynamics of a typical quasi-integrable Hamiltonian system with coexisting regular and chaotic regions. We show that the recurrence time entropy (RTE) is positively correlated to the largest Lyapunov exponent, with a high correlation coefficient. We obtain a multi-modal distribution of the finite-time RTE and find that each mode corresponds to the motion around islands of different hierarchical levels.

Revealing recurrent regimes of mid-latitude atmospheric variability using novel machine learning method.

The low-frequency variability of the extratropical atmosphere involves hemispheric-scale recurring, often persistent, states known as teleconnection patterns or regimes, which can have a profound impact on predictability on intra-seasonal and longer timescales. However, reliable data-driven identification and dynamical representation of such states are still challenging problems in modeling the dynamics of the atmosphere. We present a new method, which allows us both to detect recurring regimes of atmospheric variability and to obtain dynamical variables serving as an embedding for these regimes. The method combines two approaches from nonlinear data analysis: partitioning a network of recurrent states with studying its properties by the recurrence quantification analysis and the kernel principal component analysis. We apply the method to study teleconnection patterns in a quasi-geostrophical model of atmospheric circulation over the extratropical hemisphere as well as to reanalysis data of geopotential height anomalies in the mid-latitudes of the Northern Hemisphere atmosphere in the winter seasons from 1981 to the present. It is shown that the detected regimes as well as the obtained set of dynamical variables explain large-scale weather patterns, which are associated, in particular, with severe winters over Eurasia and North America. The method presented opens prospects for improving empirical modeling and long-term forecasting of large-scale atmospheric circulation regimes.

A complex network approach to study the extreme precipitation patterns in a river basin.

The quantification of spatial propagation of extreme precipitation events is vital in water resources planning and disaster mitigation. However, quantifying these extreme events has always been challenging as many traditional methods are insufficient to capture the nonlinear interrelationships between extreme event time series. Therefore, it is crucial to develop suitable methods for analyzing the dynamics of extreme events over a river basin with a diverse climate and complicated topography. Over the last decade, complex network analysis emerged as a powerful tool to study the intricate spatiotemporal relationship between many variables in a compact way. In this study, we employ two nonlinear concepts of event synchronization and edit distance to investigate the extreme precipitation pattern in the Ganga river basin. We use the network degree to understand the spatial synchronization pattern of extreme rainfall and identify essential sites in the river basin with respect to potential prediction skills. The study also attempts to quantify the influence of precipitation seasonality and topography on extreme events. The findings of the study reveal that (1) the network degree is decreased in the southwest to northwest direction, (2) the timing of 50th percentile precipitation within a year influences the spatial distribution of degree, (3) the timing is inversely related to elevation, and (4) the lower elevation greatly influences connectivity of the sites. The study highlights that edit distance could be a promising alternative to analyze event-like data by incorporating event time and amplitude and constructing complex networks of climate extremes.

Spike Spectra for Recurrences.

In recurrence analysis, the τ-recurrence rate encodes the periods of the cycles of the underlying high-dimensional time series. It, thus, plays a similar role to the autocorrelation for scalar time-series in encoding temporal correlations. However, its Fourier decomposition does not have a clean interpretation. Thus, there is no satisfactory analogue to the power spectrum in recurrence analysis. We introduce a novel method to decompose the τ-recurrence rate using an over-complete basis of Dirac combs together with sparsity regularization. We show that this decomposition, the inter-spike spectrum, naturally provides an analogue to the power spectrum for recurrence analysis in the sense that it reveals the dominant periodicities of the underlying time series. We show that the inter-spike spectrum correctly identifies patterns and transitions in the underlying system in a wide variety of examples and is robust to measurement noise.

Recurrence-Based Synchronization Analysis of Weakly Coupled Bursting Neurons under External ELF Fields.

We investigate the response characteristics of a two-dimensional neuron model exposed to an externally applied extremely low frequency (ELF) sinusoidal electric field and the synchronization of neurons weakly coupled with gap junction. We find, by numerical simulations, that neurons can exhibit different spiking patterns, which are well observed in the structure of the recurrence plot (RP). We further study the synchronization between weakly coupled neurons in chaotic regimes under the influence of a weak ELF electric field. In general, detecting the phases of chaotic spiky signals is not easy by using standard methods. Recurrence analysis provides a reliable tool for defining phases even for noncoherent regimes or spiky signals. Recurrence-based synchronization analysis reveals that, even in the range of weak coupling, phase synchronization of the coupled neurons occurs and, by adding an ELF electric field, this synchronization increases depending on the amplitude of the externally applied ELF electric field. We further suggest a novel measure for RP-based phase synchronization analysis, which better takes into account the probabilities of recurrences.

Joint Trends in Flood Magnitudes and Spatial Extents Across Europe.

The magnitudes of river floods in Europe have been observed to change, but their alignment with changes in the spatial coverage or extent of individual floods has not been clear. We analyze flood magnitudes and extents for 3,872 hydrometric stations across Europe over the past five decades and classify each flood based on antecedent weather conditions. We find positive correlations between flood magnitudes and extents for 95% of the stations. In central Europe and the British Isles, the association of increasing trends in magnitudes and extents is due to a magnitude-extent correlation of precipitation and soil moisture along with a shift in the flood generating processes. The alignment of trends in flood magnitudes and extents highlights the increasing importance of transnational flood risk management.

Wavelet entropy-based evaluation of intrinsic predictability of time series.

Intrinsic predictability is imperative to quantify inherent information contained in a time series and assists in evaluating the performance of different forecasting methods to get the best possible prediction. Model forecasting performance is the measure of the probability of success. Nevertheless, model performance or the model does not provide understanding for improvement in prediction. Intuitively, intrinsic predictability delivers the highest level of predictability for a time series and informative in unfolding whether the system is unpredictable or the chosen model is a poor choice. We introduce a novel measure, the Wavelet Entropy Energy Measure (WEEM), based on wavelet transformation and information entropy for quantification of intrinsic predictability of time series. To investigate the efficiency and reliability of the proposed measure, model forecast performance was evaluated via a wavelet networks approach. The proposed measure uses the wavelet energy distribution of a time series at different scales and compares it with the wavelet energy distribution of white noise to quantify a time series as deterministic or random. We test the WEEM using a wide variety of time series ranging from deterministic, non-stationary, and ones contaminated with white noise with different noise-signal ratios. Furthermore, a relationship is developed between the WEEM and Nash-Sutcliffe Efficiency, one of the widely known measures of forecast performance. The reliability of WEEM is demonstrated by exploring the relationship to logistic map and real-world data.

Recurrence analysis of slow-fast systems.

Many complex systems exhibit periodic oscillations comprising slow-fast timescales. In such slow-fast systems, the slow and fast timescales compete to determine the dynamics. In this study, we perform a recurrence analysis on simulated signals from paradigmatic model systems as well as signals obtained from experiments, each of which exhibit slow-fast oscillations. We find that slow-fast systems exhibit characteristic patterns along the diagonal lines in the corresponding recurrence plot (RP). We discern that the hairpin trajectories in the phase space lead to the formation of line segments perpendicular to the diagonal line in the RP for a periodic signal. Next, we compute the recurrence networks (RNs) of these slow-fast systems and uncover that they contain additional features such as clustering and protrusions on top of the closed-ring structure. We show that slow-fast systems and single timescale systems can be distinguished by computing the distance between consecutive state points on the phase space trajectory and the degree of the nodes in the RNs. Such a recurrence analysis substantially strengthens our understanding of slow-fast systems, which do not have any accepted functional forms.

Disentangling the multi-scale effects of sea-surface temperatures on global precipitation: A coupled networks approach.

The oceans and atmosphere interact via a multiplicity of feedback mechanisms, shaping to a large extent the global climate and its variability. To deepen our knowledge of the global climate system, characterizing and investigating this interdependence is an important task of contemporary research. However, our present understanding of the underlying large-scale processes is greatly limited due to the manifold interactions between essential climatic variables at different temporal scales. To address this problem, we here propose to extend the application of complex network techniques to capture the interdependence between global fields of sea-surface temperature (SST) and precipitation (P) at multiple temporal scales. For this purpose, we combine time-scale decomposition by means of a discrete wavelet transform with the concept of coupled climate network analysis. Our results demonstrate the potential of the proposed approach to unravel the scale-specific interdependences between atmosphere and ocean and, thus, shed light on the emerging multiscale processes inherent to the climate system, which traditionally remain undiscovered when investigating the system only at the native resolution of existing climate data sets. Moreover, we show how the relevant spatial interdependence structures between SST and P evolve across time-scales. Most notably, the strongest mutual correlations between SST and P at annual scale (8-16 months) concentrate mainly over the Pacific Ocean, while the corresponding spatial patterns progressively disappear when moving toward longer time-scales.

Extended recurrence plot and quantification for noisy continuous dynamical systems.

One main challenge in constructing a reliable recurrence plot (RP) and, hence, its quantification [recurrence quantification analysis (RQA)] of a continuous dynamical system is the induced noise that is commonly found in observation time series. This induced noise is known to cause disrupted and deviated diagonal lines despite the known deterministic features and, hence, biases the diagonal line based RQA measures and can lead to misleading conclusions. Although discontinuous lines can be further connected by increasing the recurrence threshold, such an approach triggers thick lines in the plot. However, thick lines also influence the RQA measures by artificially increasing the number of diagonals and the length of vertical lines [e.g., Determinism ( DET ) and Laminarity ( LAM ) become artificially higher]. To take on this challenge, an extended RQA approach for accounting disrupted and deviated diagonal lines is proposed. The approach uses the concept of a sliding diagonal window with minimal window size that tolerates the mentioned deviated lines and also considers a specified minimal lag between points as connected. This is meant to derive a similar determinism indicator for noisy signal where conventional RQA fails to capture. Additionally, an extended local minima approach to construct RP is also proposed to further reduce artificial block structures and vertical lines that potentially increase the associated RQA like LAM. The methodology and applicability of the extended local minima approach and DET equivalent measure are presented and discussed, respectively.