Refine
Year of publication
Document Type
- Article (32)
- Postprint (2)
- Habilitation Thesis (1)
- Review (1)
Language
- English (36)
Is part of the Bibliography
- yes (36)
Keywords
- Recurrence plot (4)
- Complex networks (3)
- Event synchronization (3)
- Extreme rainfall (2)
- Holocene (2)
- classification (2)
- climate networks (2)
- events (2)
- identifying influential nodes (2)
- precipitation (2)
- rainfall (2)
- streamflow (2)
- synchronization (2)
- Amino acid (1)
- Andes (1)
- Anisotropy (1)
- Baiu (1)
- Biomedical engineering (1)
- Brain injuries (1)
- Causality (1)
- Climate reconstruction (1)
- Cold air surges (1)
- Complex network (1)
- Computational models (1)
- Computational neuroscience (1)
- Convective storms (1)
- Data acquisition (1)
- Data processing (1)
- Dominant link directions (1)
- Eastern Africa (1)
- Electrical and electronic engineering (1)
- Embedding (1)
- Extreme events (1)
- Hydroclimate (1)
- Indian Summer Monsoon (1)
- Indian monsoon (1)
- Isochrones (1)
- Lake sediment (1)
- Mediterranean Sea (1)
- Mesoscale systems (1)
- Neural circuits (1)
- Optimization (1)
- Paleoclimate dynamics (1)
- Paläoklima (1)
- Pleistocene (1)
- Precipitation (1)
- Precipitation events (1)
- Prediction (1)
- Radiocarbon (1)
- Rainfall (1)
- Rainfall network (1)
- Recurrence analysis (1)
- Recurrence network (1)
- South American monsoon system (1)
- Stable carbon isotope (1)
- Stalagmite (1)
- State space reconstruction (1)
- Subtropical cyclones (1)
- Synchronization (1)
- Time series analysis (1)
- Time-series analysis (1)
- Trace elements (1)
- Tropical storms (1)
- Tropics (1)
- Visual system (1)
- Wiederkehr (1)
- Z-P approach (1)
- artefact avoidance (1)
- carbon isotopes (1)
- climate change (1)
- complex systems (1)
- data archeology (1)
- deglaciation (1)
- komplexe Systeme (1)
- lethal violence (1)
- oxygen isotopes (1)
- palaeoclimate (1)
- phase space time delay embedding reconstruction (1)
- pitfalls (1)
- population pressure (1)
- recurrence (1)
- recurrence network (1)
- recurrence plot (1)
- recurrence quantification analysis (1)
- räumliche Wiederkehr (1)
- sea level (1)
- spatial recurrence (1)
- stalagmite (1)
- tide gauge (1)
- time series analysis (1)
- warfare (1)
Institute
- Institut für Geowissenschaften (36) (remove)
This paper employs a complex network approach to determine the topology and evolution of the network of extreme precipitation that governs the organization of extreme rainfall before, during, and after the Indian Summer Monsoon (ISM) season. We construct networks of extreme rainfall events during the ISM (June-September), post-monsoon (October-December), and pre-monsoon (March-May) periods from satellite-derived (Tropical Rainfall Measurement Mission, TRMM) and rain-gauge interpolated (Asian Precipitation Highly Resolved Observational Data Integration Towards the Evaluation of Water Resources, APHRODITE) data sets. The structure of the networks is determined by the level of synchronization of extreme rainfall events between different grid cells throughout the Indian subcontinent. Through the analysis of various complex-network metrics, we describe typical repetitive patterns in North Pakistan (NP), the Eastern Ghats (EG), and the Tibetan Plateau (TP). These patterns appear during the pre-monsoon season, evolve during the ISM, and disappear during the post-monsoon season. These are important meteorological features that need further attention and that may be useful in ISM timing and strength prediction.
To calibrate delta O-18 time-series from speleothems in the eastern Indian summer monsoon (ISM) region of India, and to understand the moisture regime over the northern Bay of Bengal (BoB) we analyze the delta O-18 and delta D of rainwater, collected in 2007 and 2008 near Cherrapunji, India. delta D values range from + 18.5 parts per thousand to 144.4 parts per thousand, while delta O-18 varies between +0.8 parts per thousand and 18.8 parts per thousand. The Local Meteoric Water Line (LMWL) is found to be indistinguishable from the Global Meteoric Water Line (GMWL). Late ISM (September-October) rainfall exhibits lowest delta O-18 and delta D values, with little relationship to the local precipitation amount. There is a trend to lighter isotope values over the course of the ISM, but it does not correlate with the patterns of temperature and rainfall amount delta O-18 and delta D time-series have to be interpreted with caution in terms of the 'amount effect' in this subtropical region. We find that the temporal trend in delta O-18 reflects increasing transport distance during the ISM, isotopic changes in the northern BoB surface waters during late ISM, and vapor re-equilibration with rain droplets. Using an isotope box model for surface ocean waters, we quantify the potential influence of river runoff on the isotopic composition of the seasonal freshwater plume in the northern BoB. Temporal variations in this source can contribute up to 25% of the observed changes in stable isotopes of precipitation in NE India. To delineate other moisture sources, we use backward trajectory computations and find a strong correlation between source region and isotopic composition. Palaeoclimatic stable isotope time-series from northeast Indian speleothems likely reflect changes in moisture source and transport pathway, as well as the isotopic composition of the BoB surface water, all of which in turn reflect ISM strength. Stalagmite records from the region can therefore be interpreted as integrated measures of the ISM strength.
The South American Andes are frequently exposed to intense rainfall events with varying moisture sources and precipitation-forming processes. In this study, we assess the spatiotemporal characteristics and geographical origins of rainfall over the South American continent. Using high-spatiotemporal resolution satellite data (TRMM 3B42 V7), we define four different types of rainfall events based on their (1) high magnitude, (2) long temporal extent, (3) large spatial extent, and (4) high magnitude, long temporal and large spatial extent combined. In a first step, we analyze the spatiotemporal characteristics of these events over the entire South American continent and integrate their impact for the main Andean hydrologic catchments. Our results indicate that events of type 1 make the overall highest contributions to total seasonal rainfall (up to 50%). However, each consecutive episode of the infrequent events of type 4 still accounts for up to 20% of total seasonal rainfall in the subtropical Argentinean plains. In a second step, we employ complex network theory to unravel possibly non-linear and long-ranged climatic linkages for these four event types on the high-elevation Altiplano-Puna Plateau as well as in the main river catchments along the foothills of the Andes. Our results suggest that one to two particularly large squall lines per season, originating from northern Brazil, indirectly trigger large, long-lasting thunderstorms on the Altiplano Plateau. In general, we observe that extreme rainfall in the catchments north of approximately 20 degrees S typically originates from the Amazon Basin, while extreme rainfall at the eastern Andean foothills south of 20 degrees S and the Puna Plateau originates from southeastern South America.
The appropriate selection of recurrence thresholds is a key problem in applications of recurrence quantification analysis and related methods across disciplines. Here, we discuss the distribution of pairwise distances between state vectors in the studied system’s state space reconstructed by means of time-delay embedding as the key characteristic that should guide the corresponding choice for obtaining an adequate resolution of a recurrence plot. Specifically, we present an empirical description of the distance distribution, focusing on characteristic changes of its shape with increasing embedding dimension. Our results suggest that selecting the recurrence threshold according to a fixed percentile of this distribution reduces the dependence of recurrence characteristics on the embedding dimension in comparison with other commonly used threshold selection methods. Numerical investigations on some paradigmatic model systems with time-dependent parameters support these empirical findings.
Recurrence plots (RPs) provide an intuitive tool for visualizing the (potentially multi-dimensional) trajectory of a dynamical system in state space. In case only univariate observations of the system’s overall state are available, time-delay embedding has become a standard procedure for qualitatively reconstructing the dynamics in state space. The selection of a threshold distance 𝜀
, which distinguishes close from distant pairs of (reconstructed) state vectors, is known to have a substantial impact on the recurrence plot and its quantitative characteristics, but its corresponding interplay with the embedding dimension has not yet been explicitly addressed. Here, we point out that the results of recurrence quantification analysis (RQA) and related methods are qualitatively robust under changes of the (sufficiently high) embedding dimension only if the full distribution of pairwise distances between state vectors is considered for selecting 𝜀, which is achieved by consideration of a fixed recurrence rate.
The habilitation deals with the numerical analysis of the recurrence properties of geological and climatic processes. The recurrence of states of dynamical processes can be analysed with recurrence plots and various recurrence quantification options. In the present work, the meaning of the structures and information contained in recurrence plots are examined and described. New developments have led to extensions that can be used to describe the recurring patterns in both space and time. Other important developments include recurrence plot-based approaches to identify abrupt changes in the system's dynamics, to detect and investigate external influences on the dynamics of a system, the couplings between different systems, as well as a combination of recurrence plots with the methodology of complex networks. Typical problems in geoscientific data analysis, such as irregular sampling and uncertainties, are tackled by specific modifications and additions. The development of a significance test allows the statistical evaluation of quantitative recurrence analysis, especially for the identification of dynamical transitions. Finally, an overview of typical pitfalls that can occur when applying recurrence-based methods is given and guidelines on how to avoid such pitfalls are discussed. In addition to the methodological aspects, the application potential especially for geoscientific research questions is discussed, such as the identification and analysis of transitions in past climates, the study of the influence of external factors to ecological or climatic systems, or the analysis of landuse dynamics based on remote sensing data.
Quantifying the roles of single stations within homogeneous regions using complex network analysis
(2018)
Regionalization and pooling stations to form homogeneous regions or communities are essential for reliable parameter transfer, prediction in ungauged basins, and estimation of missing information. Over the years, several clustering methods have been proposed for regional analysis. Most of these methods are able to quantify the study region in terms of homogeneity but fail to provide microscopic information about the interaction between communities, as well as about each station within the communities. We propose a complex network-based approach to extract this valuable information and demonstrate the potential of our approach using a rainfall network constructed from the Indian gridded daily precipitation data. The communities were identified using the network-theoretical community detection algorithm for maximizing the modularity. Further, the grid points (nodes) were classified into universal roles according to their pattern of within- and between-community connections. The method thus yields zoomed-in details of individual rainfall grids within each community.
Based on high-spatiotemporal-resolution data, the authors perform a climatological study of strong rainfall events propagating from southeastern South America to the eastern slopes of the central Andes during the monsoon season. These events account for up to 70% of total seasonal rainfall in these areas. They are of societal relevance because of associated natural hazards in the form of floods and landslides, and they form an intriguing climatic phenomenon, because they propagate against the direction of the low-level moisture flow from the tropics. The responsible synoptic mechanism is analyzed using suitable composites of the relevant atmospheric variables with high temporal resolution. The results suggest that the low-level inflow from the tropics, while important for maintaining sufficient moisture in the area of rainfall, does not initiate the formation of rainfall clusters. Instead, alternating low and high pressure anomalies in midlatitudes, which are associated with an eastward-moving Rossby wave train, in combination with the northwestern Argentinean low, create favorable pressure and wind conditions for frontogenesis and subsequent precipitation events propagating from southeastern South America toward the Bolivian Andes.
A novel idea for an optimal time delay state space reconstruction from uni- and multivariate time series is presented. The entire embedding process is considered as a game, in which each move corresponds to an embedding cycle and is subject to an evaluation through an objective function. This way the embedding procedure can be modeled as a tree, in which each leaf holds a specific value of the objective function. By using a Monte Carlo ansatz, the proposed algorithm populates the tree with many leafs by computing different possible embedding paths and the final embedding is chosen as that particular path, which ends at the leaf with the lowest achieved value of the objective function. The method aims to prevent getting stuck in a local minimum of the objective function and can be used in a modular way, enabling practitioners to choose a statistic for possible delays in each embedding cycle as well as a suitable objective function themselves. The proposed method guarantees the optimization of the chosen objective function over the parameter space of the delay embedding as long as the tree is sampled sufficiently. As a proof of concept, we demonstrate the superiority of the proposed method over the classical time delay embedding methods using a variety of application examples. We compare recurrence plot-based statistics inferred from reconstructions of a Lorenz-96 system and highlight an improved forecast accuracy for map-like model data as well as for palaeoclimate isotope time series. Finally, we utilize state space reconstruction for the detection of causality and its strength between observables of a gas turbine type thermoacoustic combustor.
Hydrometric networks play a vital role in providing information for decision-making in water resource management. They should be set up optimally to provide as much information as possible that is as accurate as possible and, at the same time, be cost-effective. Although the design of hydrometric networks is a well-identified problem in hydrometeorology and has received considerable attention, there is still scope for further advancement. In this study, we use complex network analysis, defined as a collection of nodes interconnected by links, to propose a new measure that identifies critical nodes of station networks. The approach can support the design and redesign of hydrometric station networks. The science of complex networks is a relatively young field and has gained significant momentum over the last few years in different areas such as brain networks, social networks, technological networks, or climate networks. The identification of influential nodes in complex networks is an important field of research. We propose a new node-ranking measure – the weighted degree–betweenness (WDB) measure – to evaluate the importance of nodes in a network. It is compared to previously proposed measures used on synthetic sample networks and then applied to a real-world rain gauge network comprising 1229 stations across Germany to demonstrate its applicability. The proposed measure is evaluated using the decline rate of the network efficiency and the kriging error. The results suggest that WDB effectively quantifies the importance of rain gauges, although the benefits of the method need to be investigated in more detail.
Hydrometric networks play a vital role in providing information for decision-making in water resource management. They should be set up optimally to provide as much information as possible that is as accurate as possible and, at the same time, be cost-effective. Although the design of hydrometric networks is a well-identified problem in hydrometeorology and has received considerable attention, there is still scope for further advancement. In this study, we use complex network analysis, defined as a collection of nodes interconnected by links, to propose a new measure that identifies critical nodes of station networks. The approach can support the design and redesign of hydrometric station networks. The science of complex networks is a relatively young field and has gained significant momentum over the last few years in different areas such as brain networks, social networks, technological networks, or climate networks. The identification of influential nodes in complex networks is an important field of research. We propose a new node-ranking measure – the weighted degree–betweenness (WDB) measure – to evaluate the importance of nodes in a network. It is compared to previously proposed measures used on synthetic sample networks and then applied to a real-world rain gauge network comprising 1229 stations across Germany to demonstrate its applicability. The proposed measure is evaluated using the decline rate of the network efficiency and the kriging error. The results suggest that WDB effectively quantifies the importance of rain gauges, although the benefits of the method need to be investigated in more detail.
Non-linear time series analysis of precipitation events using regional climate networks for Germany
(2016)
Synchronous occurrences of heavy rainfall events and the study of their relation in time and space are of large socio-economical relevance, for instance for the agricultural and insurance sectors, but also for the general well-being of the population. In this study, the spatial synchronization structure is analyzed as a regional climate network constructed from precipitation event series. The similarity between event series is determined by the number of synchronous occurrences. We propose a novel standardization of this number that results in synchronization scores which are not biased by the number of events in the respective time series. Additionally, we introduce a new version of the network measure directionality that measures the spatial directionality of weighted links by also taking account of the effects of the spatial embedding of the network. This measure provides an estimate of heavy precipitation isochrones by pointing out directions along which rainfall events synchronize. We propose a climatological interpretation of this measure in terms of propagating fronts or event traces and confirm it for Germany by comparing our results to known atmospheric circulation patterns.
Sea surface temperature (SST) patterns can – as surface climate forcing – affect weather and climate at large distances. One example is El Niño-Southern Oscillation (ENSO) that causes climate anomalies around the globe via teleconnections. Although several studies identified and characterized these teleconnections, our understanding of climate processes remains incomplete, since interactions and feedbacks are typically exhibited at unique or multiple temporal and spatial scales. This study characterizes the interactions between the cells of a global SST data set at different temporal and spatial scales using climate networks. These networks are constructed using wavelet multi-scale correlation that investigate the correlation between the SST time series at a range of scales allowing instantaneously deeper insights into the correlation patterns compared to traditional methods like empirical orthogonal functions or classical correlation analysis. This allows us to identify and visualise regions of – at a certain timescale – similarly evolving SSTs and distinguish them from those with long-range teleconnections to other ocean regions. Our findings re-confirm accepted knowledge about known highly linked SST patterns like ENSO and the Pacific Decadal Oscillation, but also suggest new insights into the characteristics and origins of long-range teleconnections like the connection between ENSO and Indian Ocean Dipole.
Sea surface temperature (SST) patterns can – as surface climate forcing – affect weather and climate at large distances. One example is El Niño-Southern Oscillation (ENSO) that causes climate anomalies around the globe via teleconnections. Although several studies identified and characterized these teleconnections, our understanding of climate processes remains incomplete, since interactions and feedbacks are typically exhibited at unique or multiple temporal and spatial scales. This study characterizes the interactions between the cells of a global SST data set at different temporal and spatial scales using climate networks. These networks are constructed using wavelet multi-scale correlation that investigate the correlation between the SST time series at a range of scales allowing instantaneously deeper insights into the correlation patterns compared to traditional methods like empirical orthogonal functions or classical correlation analysis. This allows us to identify and visualise regions of – at a certain timescale – similarly evolving SSTs and distinguish them from those with long-range teleconnections to other ocean regions. Our findings re-confirm accepted knowledge about known highly linked SST patterns like ENSO and the Pacific Decadal Oscillation, but also suggest new insights into the characteristics and origins of long-range teleconnections like the connection between ENSO and Indian Ocean Dipole.
Multiple landslide clusters record quaternary climate changes in the northwestern Argentine andes
(2003)
The chronology of multiple landslide deposits and related lake sediments in the semi-arid eastern Argentine Cordillera suggests that major mass movements cluster in two time periods during the Quaternary, i.e. between 40 and 25 and after 5 14C kyr BP. These clusters may correspond to the Minchin (maximum at around 28-27 14C kyr BP) and Titicaca wet periods (after 3.9 14C kyr BP). The more humid conditions apparently caused enhanced landsliding in this environment. In contrast, no landslide-related damming and associated lake sediments occurred during the Coipasa (11.5- 10 14C yr BP) and Tauca wet periods (14.5-11 14C yr BP). The two clusters at 40-25 and after 5 14C kyr BP may correspond to periods where the El Niño-Southern Oscillation (ENSO) and Tropical Atlantic Sea Surface Temperature Variability (TAV) were active. This, however, was not the case during the Coipasa and Tauca wet periods. Lake-balance modelling of a landslide-dammed lake suggests a 10-15% increase in precipitation and a 3-4 ° C decrease in temperature at ~30 14C kyr BP as compared to the present. In addition, time-series analysis reveals a strong ENSO and TAV during that time. The landslide clusters in northwestern Argentina are therefore best explained by periods of more humid and more variable climates.
The temporal dynamics of climate processes are spread across different timescales and, as such, the study of these processes at only one selected timescale might not reveal the complete mechanisms and interactions within and between the (sub-) processes. To capture the non-linear interactions between climatic events, the method of event synchronization has found increasing attention recently. The main drawback with the present estimation of event synchronization is its restriction to analysing the time series at one reference timescale only. The study of event synchronization at multiple scales would be of great interest to comprehend the dynamics of the investigated climate processes. In this paper, the wavelet-based multi-scale event synchronization (MSES) method is proposed by combining the wavelet transform and event synchronization. Wavelets are used extensively to comprehend multi-scale processes and the dynamics of processes across various timescales. The proposed method allows the study of spatio-temporal patterns across different timescales. The method is tested on synthetic and real-world time series in order to check its replicability and applicability. The results indicate that MSES is able to capture relationships that exist between processes at different timescales.
We present the results of biogeochemical and mineralogical analyses on a sediment core that covers the Holocene sedimentation history of the climatically sensitive, closed, saline, and alkaline Lonar Lake in the core monsoon zone in central India. We compare our results of C/N ratios, stable carbon and nitrogen isotopes, grain-size, as well as amino acid derived degradation proxies with climatically sensitive proxies of other records from South Asia and the North Atlantic region. The comparison reveals some more or less contemporaneous climate shifts. At Lonar Lake, a general long term climate transition from wet conditions during the early Holocene to drier conditions during the late Holocene, delineating the insolation curve, can be reconstructed. In addition to the previously identified periods of prolonged drought during 4.6-3.9 and 2.0-0.6 cal ka that have been attributed to temperature changes in the Indo Pacific Warm Pool, several additional phases of shorter term climate alteration superimposed upon the general climate trend can be identified. These correlate with cold phases in the North Atlantic region. The most pronounced climate deteriorations indicated by our data occurred during 62-5.2,4.6-3.9, and 2.0-0.6 cal ka BP. The strong dry phase between 4.6 and 3.9 cal ka BP at Lonar Lake corroborates the hypothesis that severe climate deterioration contributed to the decline of the Indus Civilisation about 3.9 ka BP. (C) 2014 Elsevier B.V. All rights reserved.
As an effort to reduce parameter uncertainties in constructing recurrence plots, and in particular to avoid potential artefacts, this paper presents a technique to derive artefact-safe region of parameter sets. This technique exploits both deterministic (incl. chaos) and stochastic signal characteristics of recurrence quantification (i.e. diagonal structures). It is useful when the evaluated signal is known to be deterministic. This study focuses on the recurrence plot generated from the reconstructed phase space in order to represent many real application scenarios when not all variables to describe a system are available (data scarcity). The technique involves random shuffling of the original signal to destroy its original deterministic characteristics. Its purpose is to evaluate whether the determinism values of the original and the shuffled signal remain closely together, and therefore suggesting that the recurrence plot might comprise artefacts. The use of such determinism-sensitive region shall be accompanied by standard embedding optimization approaches, e.g. using indices like false nearest neighbor and mutual information, to result in a more reliable recurrence plot parameterization.
The analysis of palaeoclimate time series is usually affected by severe methodological problems, resulting primarily from non-equidistant sampling and uncertain age models. As an alternative to existing methods of time series analysis, in this paper we argue that the statistical properties of recurrence networks - a recently developed approach - are promising candidates for characterising the system's nonlinear dynamics and quantifying structural changes in its reconstructed phase space as time evolves. In a first order approximation, the results of recurrence network analysis are invariant to changes in the age model and are not directly affected by non-equidistant sampling of the data. Specifically, we investigate the behaviour of recurrence network measures for both paradigmatic model systems with non-stationary parameters and four marine records of long-term palaeoclimate variations. We show that the obtained results are qualitatively robust under changes of the relevant parameters of our method, including detrending, size of the running window used for analysis, and embedding delay. We demonstrate that recurrence network analysis is able to detect relevant regime shifts in synthetic data as well as in problematic geoscientific time series. This suggests its application as a general exploratory tool of time series analysis complementing existing methods.
Hydrological and climatological controls on radiocarbon concentrations in a tropical stalagmite
(2016)
Precisely-dated stalagmites are increasingly important archives for the reconstruction of terrestrial paleoclimate at very high temporal resolution. In-depth understanding of local conditions at the cave site and of the processes driving stalagmite deposition is of paramount importance for interpreting proxy signals incorporated in stalagmite carbonate. Here we present a sub-decadally resolved dead carbon fraction (DCF) record for a stalagmite from Yok Balum Cave (southern Belize). The record is coupled to parallel stable carbon isotope (delta C-13) and U/Ca measurements, as well as radiocarbon (C-14) measurements from soils overlying the cave system. Using a karst carbon cycle model we disentangle the importance of soil and karst processes on stalagmite DCF incorporation, revealing a dominant host rock dissolution control on total DCF. Covariation between DCF, delta C-13, and U/Ca indicates that karst processes are a common driver of all three parameters, suggesting possible use of delta C-13 and trace element ratios to independently quantify DCF variability. A statistically significant multi-decadal lag of variable length exists between DCF and reconstructed solar activity, suggesting that solar activity influenced regional precipitation in Mesoamerica over the past 1500 years, but that the relationship was non-static. Although the precise nature of the observed lag is unclear, solar-induced changes in North Atlantic oceanic and atmospheric dynamics may play a role. (C) 2016 Elsevier Ltd. All rights reserved.
Recurrence plots and recurrence quantification analysis have become popular in the last two decades. Recurrence based methods have on the one hand a deep foundation in the theory of dynamical systems and are on the other hand powerful tools for the investigation of a variety of problems. The increasing interest encompasses the growing risk of misuse and uncritical application of these methods. Therefore, we point out potential problems and pitfalls related to different aspects of the application of recurrence plots and recurrence quantification analysis.
In this paper, we present the new frequency spectrum recurrence analysis technique by means of electro-encephalon signals (EES) analyses. The technique is suitable for time series analysis with noise and disturbances. EES were collected, and alpha waves of the occipital region were analysed by comparing the signals from participants in two states, eyes open and eyes closed. Firstly, EES were characterized and analysed by means of techniques already known to compare with the results of the innovative technique that we present here. We verified that, standard recurrence quantification analysis by means of EES time series cannot statistically distinguish the two states. However, the new frequency spectrum recurrence quantification exhibit quantitatively whether the participants have their eyes open or closed. In sequence, new quantifiers are created for analysing the recurrence concentration on frequency bands. These analyses show that EES with similar frequency spectrum have different recurrence levels revealing different behaviours of the nervous system. The technique can be used to deepen the study on depression, stress, concentration level and other neurological issues and also can be used in any complex system.
Recurrence-plot-based recurrence networks are an approach used to analyze time series using a complex networks theory. In both approaches - recurrence plots and recurrence networks -, a threshold to identify recurrent states is required. The selection of the threshold is important in order to avoid bias of the recurrence network results. In this paper, we propose a novel method to choose a recurrence threshold adaptively. We show a comparison between the constant threshold and adaptive threshold cases to study period-chaos and even period-period transitions in the dynamics of a prototypical model system. This novel method is then used to identify climate transitions from a lake sediment record.
Extreme Rainfall of the South American Monsoon System: A Dataset Comparison Using Complex Networks
(2015)
In this study, the authors compare six different rainfall datasets for South America with a focus on their representation of extreme rainfall during the monsoon season (December February): the gauge-calibrated TRMM 3B42 V7 satellite product; the near-real-time TRMM 3B42 V7 RT, the GPCP 1 degrees daily (1DD) V1.2 satellite gauge combination product, the Interim ECMWF Re-Analysis (ERA-Interim) product; output of a high-spatial-resolution run of the ECHAM6 global circulation model; and output of the regional climate model Eta. For the latter three, this study can be understood as a model evaluation. In addition to statistical values of local rainfall distributions, the authors focus on the spatial characteristics of extreme rainfall covariability. Since traditional approaches based on principal component analysis are not applicable in the context of extreme events, they apply and further develop methods based on complex network theory. This way, the authors uncover substantial differences in extreme rainfall patterns between the different datasets: (i) The three model-derived datasets yield very different results than the satellite gauge combinations regarding the main climatological propagation pathways of extreme events as well as the main convergence zones of the monsoon system. (ii) Large discrepancies are found for the development of mesoscale convective systems in southeastern South America. (iii) Both TRMM datasets and ECHAM6 indicate a linkage of extreme rainfall events between the central Amazon basin and the eastern slopes of the central Andes, but this pattern is not reproduced by the remaining datasets. The authors' study suggests that none of the three model-derived datasets adequately captures extreme rainfall patterns in South America.
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.
Sedimentary proxy records constitute a significant portion of the recorded evidence that allows us to investigate paleoclimatic conditions and variability. However, uncertainties in the dating of proxy archives limit our ability to fix the timing of past events and interpret proxy record intercomparisons. While there are various age-modeling approaches to improve the estimation of the age-depth relations of archives, relatively little focus has been placed on the propagation of the age (and radiocarbon calibration) uncertainties into the final proxy record.
We present a generic Bayesian framework to estimate proxy records along with their associated uncertainty, starting with the radiometric age-depth and proxy-depth measurements, and a radiometric calibration curve if required. We provide analytical expressions for the posterior proxy probability distributions at any given calendar age, from which the expected proxy values and their uncertainty can be estimated. We illustrate our method using two synthetic data sets and then use it to construct the proxy records for groundwater inflow and surface erosion from Lonar lake in central India.
Our analysis reveals interrelations between the uncertainty of the proxy record over time and the variance of proxies along the depth of the archive. For the Lonar lake proxies, we show that, rather than the age uncertainties, it is the proxy variance combined with calibration uncertainty that accounts for most of the final uncertainty. We represent the proxy records as probability distributions on a precise, error-free timescale that makes further time series analyses and intercomparisons of proxies relatively simple and clear. Our approach provides a coherent understanding of age uncertainties within sedimentary proxy records that involve radiometric dating. It can be potentially used within existing age modeling structures to bring forth a reliable and consistent framework for proxy record estimation.
Climatic changes are of major importance in landslide generation in the Argentine Andes. Increased humidity as a potential influential factor was inferred from the temporal clustering of landslide deposits during a period of significantly wetter climate, 30,000 years ago. A change in seasonality was tested by comparing past (inferred from annual-layered lake deposits, 30,000 years old) and modern (present-day observations) precipitation changes. Quantitative analysis of cross recurrence plots were developed to compare the influence of the El Nino/Southern Oscillation (ENSO) on present and past rainfall variations. This analysis has shown the stronger influence of NE trades in the location of landslide deposits in the intra-andean basin and valleys, what caused a higher contrast between summer and winter rainfall and an increasing of precipitation in La Nina years. This is believed to reduce thresholds for landslide generation in the arid to semiarid intra-andean basins and valleys.
Reliable age models are fundamental for any palaeoclimate reconstruction. Available interpolation procedures between age control points are often inadequately reported, and very few translate age uncertainties to proxy uncertainties. Most available modeling algorithms do not allow incorporation of layer counted intervals to improve the confidence limits of the age model in question.
We present a framework that allows detection and interactive handling of age reversals and hiatuses, depth-age modeling, and proxy-record reconstruction. Monte Carlo simulation and a translation procedure are used to assign a precise time scale to climate proxies and to translate dating uncertainties to uncertainties in the proxy values. The presented framework allows integration of incremental relative dating information to improve the final age model. The free software package COPRA1.0 facilitates easy interactive usage.
Northeastern (NE) India experiences extraordinarily pronounced seasonal climate, governed by the Indian summer monsoon (ISM). The vulnerability of this region to floods and droughts calls for detailed and highly resolved paleoclimate reconstructions to assess the recurrence rate and driving factors of ISM changes. We use stable oxygen and carbon isotope ratios (delta O-18 and delta C-13) from stalagmite MAW-6 from Mawmluh Cave to infer climate and environmental conditions in NE India over the last deglaciation (16-6ka). We interpret stalagmite delta O-18 as reflecting ISM strength, whereas delta C-13 appears to be driven by local hydroclimate conditions. Pronounced shifts in ISM strength over the deglaciation are apparent from the delta O-18 record, similarly to other records from monsoonal Asia. The ISM is weaker during the late glacial (LG) period and the Younger Dryas, and stronger during the BOlling-Allerod and Holocene. Local conditions inferred from the delta C-13 record appear to have changed less substantially over time, possibly related to the masking effect of changing precipitation seasonality. Time series analysis of the delta O-18 record reveals more chaotic conditions during the late glacial and higher predictability during the Holocene, likely related to the strengthening of the seasonal recurrence of the ISM with the onset of the Holocene.
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 C-14-dated individuals from the Nasca highland region of Peru. Multiple local and supraregional precipitation proxies are combined with a summed probability distribution of 149 C-14 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.
The Chew Bahir Drilling Project (CBDP) aims to test possible linkages between climate and evolution in Africa through the analysis of sediment cores that have recorded environmental changes in the Chew Bahir basin. In this statistical project we consider the Chew Bahir palaeolake to be a dynamical system consisting of interactions between its different components, such as the waterbody, the sediment beneath lake, and the organisms living within and around the lake. Recurrence is a common feature of such dynamical systems, with recurring patterns in the state of the system reflecting typical influences. Identifying and defining these influences contributes significantly to our understanding of the dynamics of the system. Different recurring changes in precipitation, evaporation, and wind speed in the Chew Bahir basin could result in similar (but not identical) conditions in the lake (e.g., depth and area of the lake, alkalinity and salinity of the lake water, species assemblages in the water body, and diagenesis in the sediments). Recurrence plots (RPs) are graphic displays of such recurring states within a system. Measures of complexity were subsequently introduced to complement the visual inspection of recurrence plots, and provide quantitative descriptions for use in recurrence quantification analysis (RQA). We present and discuss herein results from an RQA on the environmental record from six short (< 17 m) sediment cores collected during the CBDP, spanning the last 45 kyrs. The different types of variability and transitions in these records were classified to improve our understanding of the response of the biosphere to climate change, and especially the response of humans in the area.
In this paper a change-point detection method is proposed by extending the singular spectrum transformation (SST) developed as one of the capabilities of singular spectrum analysis (SSA). The method uncovers change points related with trends and periodicities. The potential of the proposed method is demonstrated by analysing simple model time series including linear functions and sine functions as well as real world data (precipitation data in Kenya). A statistical test of the results is proposed based on a Monte Carlo simulation with surrogate methods. As a result, the successful estimation of change points as inherent properties in the representative time series of both trend and harmonics is shown. With regards to the application, we find change points in the precipitation data of Kenyan towns (Nakuru, Naivasha, Narok, and Kisumu) which coincide with the variability of the Indian Ocean Dipole (IOD) suggesting its impact of extreme climate in East Africa.
Identifying abrupt transitions is a key question in various disciplines. Existing transition detection methods, however, do not rigorously account for time series uncertainties, often neglecting them altogether or assuming them to be independent and qualitatively similar. Here, we introduce a novel approach suited to handle uncertainties by representing the time series as a time-ordered sequence of probability density functions. We show how to detect abrupt transitions in such a sequence using the community structure of networks representing probabilities of recurrence. Using our approach, we detect transitions in global stock indices related to well-known periods of politico-economic volatility. We further uncover transitions in the El Niño-Southern Oscillation which coincide with periods of phase locking with the Pacific Decadal Oscillation. Finally, we provide for the first time an ‘uncertainty-aware’ framework which validates the hypothesis that ice-rafting events in the North Atlantic during the Holocene were synchronous with a weakened Asian summer monsoon.
Quantitative estimates of sea-level rise in the Mediterranean Basin become increasingly accurate thanks to detailed satellite monitoring. However, such measuring campaigns cover several years to decades, while longer-term sea-level records are rare for the Mediterranean. We used a data archeological approach to reanalyze monthly mean sea-level data of the Antalya-I (1935–1977) tide gauge to fill this gap. We checked the accuracy and reliability of these data before merging them with the more recent records of the Antalya-II (1985–2009) tide gauge, accounting for an eight-year hiatus. We obtain a composite time series of monthly and annual mean sea levels spanning some 75 years, providing the longest record for the eastern Mediterranean Basin, and thus an essential tool for studying the region's recent sea-level trends. We estimate a relative mean sea-level rise of 2.2 ± 0.5 mm/year between 1935 and 2008, with an annual variability (expressed here as the standard deviation of the residuals, σresiduals = 41.4 mm) above that at the closest tide gauges (e.g., Thessaloniki, Greece, σresiduals = 29.0 mm). Relative sea-level rise accelerated to 6.0 ± 1.5 mm/year at Antalya-II; we attribute roughly half of this rate (~3.6 mm/year) to tectonic crustal motion and anthropogenic land subsidence. Our study highlights the value of data archeology for recovering and integrating historic tide gauge data for long-term sea-level and climate studies.
Frequent and intense rainfall events demand innovative techniques to better predict the extreme rainfall dynamics. This task requires essentially the assessment of the basic types of atmospheric processes that trigger extreme rainfall, and then to examine the differences between those processes, which may help to identify key patterns to improve predictive algorithms. We employ tools from network theory to compare the spatial features of extreme rainfall over the Japanese archipelago and surrounding areas caused by two atmospheric processes: the Baiu front, which occurs mainly in June and July (JJ), and the tropical storms from August to November (ASON). We infer from complex networks of satellite-derived rainfall data, which are based on the nonlinear correlation measure of event synchronization. We compare the spatial scales involved in both systems and identify different regions which receive rainfall due to the large spatial scale of the Baiu and tropical storm systems. We observed that the spatial scales involved in the Baiu driven rainfall extremes, including the synoptic processes behind the frontal development, are larger than tropical storms, which even have long tracks during extratropical transitions. We further delineate regions of coherent rainfall during the two seasons based on network communities, identifying the horizontal (east-west) rainfall bands during JJ over the Japanese archipelago, while during ASON these bands align with the island arc of Japan.
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.