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.

Early warning signals of the termination of the African Humid Period(s).

The transition from a humid green Sahara to today's hyperarid conditions in northern Africa ~5.5 thousand years ago shows the dramatic environmental change to which human societies were exposed and had to adapt to. In this work, we show that in the 620,000-year environmental record from the Chew Bahir basin in the southern Ethiopian Rift, with its decadal resolution, this one thousand year long transition is particularly well documented, along with 20-80 year long droughts, recurring every ~160 years, as possible early warnings. Together with events of extreme wetness at the end of the transition, these droughts form a pronounced climate "flickering", which can be simulated in climate models and is also present in earlier climate transitions in the Chew Bahir environmental record, indicating that transitions with flickering are characteristic of this region.

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.

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.

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.

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.

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.

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.

Decline in seasonal predictability potentially destabilized Classic Maya societies.

Classic Maya populations living in peri-urban states were highly dependent on seasonally distributed rainfall for reliable surplus crop yields. Despite intense study of the potential impact of decadal to centennial-scale climatic changes on the demise of Classic Maya sociopolitical institutions (750-950 CE), its direct importance remains debated. We provide a detailed analysis of a precisely dated speleothem record from Yok Balum cave, Belize, that reflects local hydroclimatic changes at seasonal scale over the past 1600 years. We find that the initial disintegration of Maya sociopolitical institutions and population decline occurred in the context of a pronounced decrease in the predictability of seasonal rainfall and severe drought between 700 and 800 CE. The failure of Classic Maya societies to successfully adapt to volatile seasonal rainfall dynamics likely contributed to gradual but widespread processes of sociopolitical disintegration. We propose that the complex abandonment of Classic Maya population centres was not solely driven by protracted drought but also aggravated by year-to-year decreases in rainfall predictability, potentially caused by a regional reduction in coherent Intertropical Convergence Zone-driven rainfall.

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.

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.

Drought-Induced Civil Conflict Among the Ancient Maya.

The influence of climate change on civil conflict and societal instability in the premodern world is a subject of much debate, in part because of the limited temporal or disciplinary scope of case studies. We present a transdisciplinary case study that combines archeological, historical, and paleoclimate datasets to explore the dynamic, shifting relationships among climate change, civil conflict, and political collapse at Mayapan, the largest Postclassic Maya capital of the Yucatán Peninsula in the thirteenth and fourteenth centuries CE. Multiple data sources indicate that civil conflict increased significantly and generalized linear modeling correlates strife in the city with drought conditions between 1400 and 1450 cal. CE. We argue that prolonged drought escalated rival factional tensions, but subsequent adaptations reveal regional-scale resiliency, ensuring that Maya political and economic structures endured until European contact in the early sixteenth century CE.

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.

Sampling rate-corrected analysis of irregularly sampled time series.

The analysis of irregularly sampled time series remains a challenging task requiring methods that account for continuous and abrupt changes of sampling resolution without introducing additional biases. The edit distance is an effective metric to quantitatively compare time series segments of unequal length by computing the cost of transforming one segment into the other. We show that transformation costs generally exhibit a nontrivial relationship with local sampling rate. If the sampling resolution undergoes strong variations, this effect impedes unbiased comparison between different time episodes. We study the impact of this effect on recurrence quantification analysis, a framework that is well suited for identifying regime shifts in nonlinear time series. A constrained randomization approach is put forward to correct for the biased recurrence quantification measures. This strategy involves the generation of a type of time series and time axis surrogates which we call sampling-rate-constrained (SRC) surrogates. We demonstrate the effectiveness of the proposed approach with a synthetic example and an irregularly sampled speleothem proxy record from Niue island in the central tropical Pacific. Application of the proposed correction scheme identifies a spurious transition that is solely imposed by an abrupt shift in sampling rate and uncovers periods of reduced seasonal rainfall predictability associated with enhanced El Niño-Southern Oscillation and tropical cyclone activity.

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.

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.

Publisher's Note: Entropy of weighted recurrence plots [Phys. Rev. E 90, 042919 (2014)].

This corrects the article DOI: 10.1103/PhysRevE.90.042919.