Refine
Year of publication
Document Type
- Article (67)
- Postprint (4)
- Doctoral Thesis (1)
- Habilitation Thesis (1)
- Part of Periodical (1)
- Review (1)
Is part of the Bibliography
- yes (75)
Keywords
- Complex networks (4)
- Event synchronization (4)
- Recurrence plot (4)
- Holocene (3)
- precipitation (3)
- synchronization (3)
- Extreme rainfall (2)
- Indian monsoon (2)
- Indian summer monsoon (2)
- classification (2)
Institute
- Institut für Geowissenschaften (36)
- Institut für Physik und Astronomie (32)
- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (4)
- Institut für Umweltwissenschaften und Geographie (3)
- Mathematisch-Naturwissenschaftliche Fakultät (2)
- Potsdam Institute for Climate Impact Research (PIK) e. V. (2)
- Dezernat 2: Studienangelegenheiten (1)
- Institut für Biochemie und Biologie (1)
The rescaling of geological data series to a geological reference time series is of major interest in many investigations. For example, geophysical borehole data should be correlated to a given data series whose time scale is known in order to achieve an age-depth function or the sedimentation rate for the borehole data. Usually this synchronization is performed visually and by hand. Instead of using this wiggle matching by eye, we present the application of cross recurrence plots for such tasks. Using this method, the synchronization and rescaling of geological data to a given time scale is much easier and faster than by hand.
The method of recurrence plots is extended to the cross recurrence plots (CRP), which among others enables the study of synchronization or time differences in two time series. This is emphasized in a distorted main diagonal in the cross recurrence plot, the line of synchronization (LOS). A non-parametrical fit of this LOS can be used to rescale the time axis of the two data series (whereby one of it is e.g. compressed or stretched) so that they are synchronized. An application of this method to geophysical sediment core data illustrates its suitability for real data. The rock magnetic data of two different sediment cores from the Makarov Basin can be adjusted to each other by using this method, so that they are comparable.
We propose a novel approach based on the fluctuation of similarity to identify regimes of distinct dynamical complexity in short time series. A statistical test is developed to estimate the significance of the identified transitions. Our method is verified by uncovering bifurcation structures in several paradigmatic models, providing more complex transitions compared with traditional Lyapunov exponents. In a real-world situation, we apply this method to identify millennial-scale dynamical transitions in Plio-Pleistocene proxy records of the South Asian summer monsoon system. We infer that many of these transitions are induced by the external forcing of the solar insolation and are also affected by internal forcing on Monsoonal dynamics, i.e., the glaciation cycles of the Northern Hemisphere and the onset of the Walker circulation.
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.
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.
High resolution reconstructions of the India Summer Monsoon (ISM) are essential to identify regionally different patterns of climate change and refine predictive models. We find opposing trends of hydrological proxies between northern (Sahiya cave stalagmite) and central India (Lonar Lake) between 100 and 1300 CE with the strongest anti-correlation between 810 and 1300 CE. The apparently contradictory data raise the question if these are related to widely different regional precipitation patterns or reflect human influence in/around the Lonar Lake. By comparing multiproxy data with historical records, we demonstrate that only the organic proxies in the Lonar Lake show evidence of anthropogenic impact. However, evaporite data (mineralogy and delta O-18) are indicative of precipitation/evaporation (P/E) into the Lonar Lake. Back-trajectories of air-mass circulation over northern and central India show that the relative contribution of the Bay of Bengal (BoB) branch of the ISM is crucial for determining the delta O-18 of carbonate proxies only in north India, whereas central India is affected significantly by the Arabian Sea (AS) branch of the ISM. We conclude that the delta O-18 of evaporative carbonates in the Lonar Lake reflects P/E and, in the interval under consideration, is not influenced by source water changes. The opposing trend between central and northern India can be explained by (i) persistent multidecadal droughts over central India between 810 and 1300 CE that provided an effective mechanism for strengthening sub-tropical westerly winds resulting in enhancement of wintertime (non-monsoonal) rainfall over northern parts of the Indian subcontinent, and/or (ii) increased moisture influx to northern India from the depleted BoB source waters.
Encounters with neighbours
(2003)
In this work, different aspects and applications of the recurrence plot analysis are presented. First, a comprehensive overview of recurrence plots and their quantification possibilities is given. New measures of complexity are defined by using geometrical structures of recurrence plots. These measures are capable to find chaos-chaos transitions in processes. Furthermore, a bivariate extension to cross recurrence plots is studied. Cross recurrence plots exhibit characteristic structures which can be used for the study of differences between two processes or for the alignment and search for matching sequences of two data series. The selected applications of the introduced techniques to various kind of data demonstrate their ability. Analysis of recurrence plots can be adopted to the specific problem and thus opens a wide field of potential applications. Regarding the quantification of recurrence plots, chaos-chaos transitions can be found in heart rate variability data before the onset of life threatening cardiac arrhythmias. This may be of importance for the therapy of such cardiac arrhythmias. The quantification of recurrence plots allows to study transitions in brain during cognitive experiments on the base of single trials. Traditionally, for the finding of these transitions the averaging of a collection of single trials is needed. Using cross recurrence plots, the existence of an El Niño/Southern Oscillation-like oscillation is traced in northwestern Argentina 34,000 yrs. ago. In further applications to geological data, cross recurrence plots are used for time scale alignment of different borehole data and for dating a geological profile with a reference data set. Additional examples from molecular biology and speech recognition emphasize the suitability of cross recurrence plots.
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 Nino-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.
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.
Precipitation patterns and extremes are significantly influenced by various climatic factors and large-scale atmospheric circulation patterns. This study uses wavelet coherence analysis to detect significant interannual and interdecadal oscillations in monthly precipitation extremes across India and their teleconnections to three prominent climate indices, namely, Nino 3.4, Pacific Decadal Oscillation, and Indian Ocean Dipole (IOD). Further, partial wavelet coherence analysis is used to estimate the standalone relationship between the climate indices and precipitation after removing the effect of interdependency. The wavelet analysis of monthly precipitation extremes at 30 different locations across India reveals that (a) interannual (2-8 years) and interdecadal (8-32 years) oscillations are statistically significant, and (b) the oscillations vary in both time and space. The results from the partial wavelet coherence analysis reveal that Nino 3.4 and IOD are the significant drivers of Indian precipitation at interannual and interdecadal scales. Intriguingly, the study also confirms that the strength of influence of large-scale atmospheric circulation patterns on Indian precipitation extremes varies with spatial physiography of the region.
Reconciling the paths of extreme rainfall with those of typhoons remains difficult despite advanced forecasting techniques. We use complex networks defined by a nonlinear synchronization measure termed event synchronization to track extreme rainfall over the Japanese islands. Directed networks objectively record patterns of heavy rain brought by frontal storms and typhoons but mask out contributions of local convective storms. We propose a radial rank method to show that paths of extreme rainfall in the typhoon season (August-November, ASON) follow the overall southwest-northeast motion of typhoons and mean rainfall gradient of Japan. The associated eye-of-the-typhoon tracks deviate notably and may thus distort estimates of heavy typhoon rainfall. We mainly found that the lower spread of rainfall tracks in ASON may enable better hindcasting than for westerly-fed frontal storms in June and July.
In recent years, complex network analysis facilitated the identification of universal and unexpected patterns in complex climate systems. However, the analysis and representation of a multiscale complex relationship that exists in the global climate system are limited. A logical first step in addressing this issue is to construct multiple networks over different timescales. Therefore, we propose to apply the wavelet multiscale correlation (WMC) similarity measure, which is a combination of two state-of-the-art methods, viz. wavelet and Pearson’s correlation, for investigating multiscale processes through complex networks. Firstly we decompose the data over different timescales using the wavelet approach and subsequently construct a corresponding network by Pearson’s correlation. The proposed approach is illustrated and tested on two synthetics and one real-world example. The first synthetic case study shows the efficacy of the proposed approach to unravel scale-specific connections, which are often undiscovered at a single scale. The second synthetic case study illustrates that by dividing and constructing a separate network for each time window we can detect significant changes in the signal structure. The real-world example investigates the behavior of the global sea surface temperature (SST) network at different timescales. Intriguingly, we notice that spatial dependent structure in SST evolves temporally. Overall, the proposed measure has an immense potential to provide essential insights on understanding and extending complex multivariate process studies at multiple scales.
A better understanding of precipitation dynamics in the Indian subcontinent is required since India's society depends heavily on reliable monsoon forecasts. We introduce a non-linear, multiscale approach, based on wavelets and event synchronization, for unravelling teleconnection influences on precipitation. We consider those climate patterns with the highest relevance for Indian precipitation. Our results suggest significant influences which are not well captured by only the wavelet coherence analysis, the state-of-the-art method in understanding linkages at multiple timescales. We find substantial variation across India and across timescales. In particular, El Niño–Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD) mainly influence precipitation in the south-east at interannual and decadal scales, respectively, whereas the North Atlantic Oscillation (NAO) has a strong connection to precipitation, particularly in the northern regions. The effect of the Pacific Decadal Oscillation (PDO) stretches across the whole country, whereas the Atlantic Multidecadal Oscillation (AMO) influences precipitation particularly in the central arid and semi-arid regions. The proposed method provides a powerful approach for capturing the dynamics of precipitation and, hence, helps improve precipitation forecasting.
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.
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.
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.
We use the extension of the method of recurrence plots to cross recurrence plots (CRP) which enables a nonlinear analysis of bivariate data. To quantify CRPs, we develop further three measures of complexity mainly basing on diagonal structures in CRPs. The CRP analysis of prototypical model systems with nonlinear interactions demonstrates that this technique enables to find these nonlinear interrelations from bivariate time series, whereas linear correlation tests do not. Applying the CRP analysis to climatological data, we find a complex relationship between rainfall and El Nino data.
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.
Ventricular tachycardia or fibrillation (VT) as fatal cardiac arrhythmias are the main factors triggering sudden cardiac death. The objective of this recurrence quantification analysis approach is to find early signs of sustained VT in patients with an implanted cardioverter-defibrillator (ICD). These devices are able to safeguard patients by returning their hearts to a normal rhythm via strong defibrillatory shocks; additionally, they are able to store at least 1000 beat-to-beat intervals immediately before the onset of a life-threatening arrhythmia. We study the
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.
Three-dimensional quantification of structures in trabecular bone using measures of complexity
(2009)
The study of pathological changes of bone is an important task in diagnostic procedures of patients with metabolic bone diseases such as osteoporosis as well as in monitoring the health state of astronauts during long-term space flights. The recent availability of high-resolution three-dimensional (3D) imaging of bone challenges the development of data analysis techniques able to assess changes of the 3D microarchitecture of trabecular bone. We introduce an approach based on spatial geometrical properties and define structural measures of complexity for 3D image analysis. These measures evaluate different aspects of organization and complexity of 3D structures, such as complexity of its surface or shape variability. We apply these measures to 3D data acquired by high-resolution microcomputed tomography (mu CT) from human proximal tibiae and lumbar vertebrae at different stages of osteoporotic bone loss. The outcome is compared to the results of conventional static histomorphometry and exhibits clear relationships between the analyzed geometrical features of trabecular bone and loss of bone density, but also indicate that the measures reveal additional information about the structural composition of bone, which were not revealed by the static histomorphometry. Finally, we have studied the dependency of the developed measures of complexity on the spatial resolution of the mu CT data sets.
Aims. Sunspot distribution in the northern and southern solar hemispheres exibit striking synchronous behaviour on the scale of a Schwabe cycle. However, sometimes the bilateral symmetry of the Butterfly diagram relative to the solar equatorial plane breaks down. The investigation of this phenomenon is important to explaining the almost-periodic behaviour of solar cycles. Methods. We use cross-recurrence plots for the study of the time-varying phase asymmetry of the northern and southern hemisphere and compare our results with the latitudinal distribution of the sunspots. Results. We observe a long-term persistence of phase leading in one of the hemispheres, which lasts almost 4 solar cycles and probably corresponds to the Gleissberg cycle. Long-term variations in the hemispheric-leading do not demonstrate clear periodicity but are strongly anti-correlated with the long-term variations in the magnetic equator.
In the recent article "Stochastic analysis of recurrence plots with applications to the detection of deterministic signals" (Physica D 237 (2008) 619-629), Rohde et al. stated that the performance of RQA in order to detect deterministic signals would be below traditional and well-known detectors. However, we have concerns about such a general statement. Based on our own studies we cannot confirm their conclusions. Our findings suggest that the measures of complexity provided by RQA are useful detectors outperforming well-known traditional detectors, in particular for the detection of signals of complex systems, with phase differences or signals modified due to the measurement process.
In the recent past, recurrence quantification analysis (RQA) has gained an increasing interest in various research areas. The complexity measures the RQA provides have been useful in describing and analysing a broad range of data. It is known to be rather robust to noise and nonstationarities. Yet, one key question in empirical research concerns the confidence bounds of measured data. In the present Letter we suggest a method for estimating the confidence bounds of recurrence-based complexity measures. We study the applicability of the suggested method with model and real- life data.
Portal alumni
(2005)
Liebe Leserin, lieber Leser, erforschen, was die Welt im Innersten zusammenhält- das ist für viele Studierende ein Traum. Doch welche Opfer muss man bringen, um ihn zu verwirklichen? Welche Bemfsperspektive hat der Bemf Forscher heute noch? Auch viele Absolventen der Universität Potsdam müssen sich diese Fragen beantworten. Zu welchen Antworten einige dabei gekommen sind und welche Probleme sie zu bewältigen haben, vom Spaß am Forschen und von Zukunftsängsten berichten sie in der Rubrik "Forscherkarrieren". Gelder für die Forschung fließen in Deutschland zu spärlich, verglichen mit anderen führenden Industrienationen. So sind die Bedingungen für Forscher hierzulande nicht die besten. Manchen jungen Wissenschaftler zieht es- mitunter notgedrungen- ins Ausland. Wie Deutschland dadurch seine ZukunftsHihigkeit riskiert, thematisiert der Präsident der Fraunhofer-Gesellschaft, Prof. Dr. Hans-Jörg Bullinger, in der Rubrik "wissenstransfer". Auch die Universität ist kein Garant für eine gesicherte Zukunft in der Forschung. Wer sechs Jahre nach der Promotion den Sprung zur Professur nicht geschafft hat, geht einer ungewissen Zukunft als Privatdozent entgegen. Seit einigen Jahren gibt es neben der Habilitation noch einen zweiten Weg zur Professur- die Juniorprofessur. Auch an der Universität Potsdam gibt es seit 2002 Juniorprofessoren, von denen die ersten jetzt evaluiert wurden. Näheres dazu finden Sie ebenfalls in der Rubrik "wissenstransfer". Wer noch nach einer Finanzierungsmöglichkeit für seine Promotion sucht, findet Tipps in der Rubrik "wegweiser". Die Redaktion wünscht Ihnen viel Vergnügen beim Lesen von Portal alumni und freut sich auf zahlreiche Leserbriefe.
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.
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.
Analysis of spatial and temporal extreme monsoonal rainfall over South Asia using complex networks
(2012)
We present a detailed analysis of summer monsoon rainfall over the Indian peninsular using nonlinear spatial correlations. This analysis is carried out employing the tools of complex networks and a measure of nonlinear correlation for point processes such as rainfall, called event synchronization. This study provides valuable insights into the spatial organization, scales, and structure of the 90th and 94th percentile rainfall events during the Indian summer monsoon (June-September). We furthermore analyse the influence of different critical synoptic atmospheric systems and the impact of the steep Himalayan topography on rainfall patterns. The presented method not only helps us in visualising the structure of the extreme-event rainfall fields, but also identifies the water vapor pathways and decadal-scale moisture sinks over the region. Furthermore a simple scheme based on complex networks is presented to decipher the spatial intricacies and temporal evolution of monsoonal rainfall patterns over the last 6 decades.
We investigate a network of influences connected to global mean temperature. Considering various climatic factors known to influence global mean temperature, we evaluate not only the impacts of these factors on temperature but also the directed dependencies among the factors themselves. Based on an existing recurrence-based connectivity measure, we propose a new and more general measure that quantifies the level of dependence between two time series based on joint recurrences at a chosen time delay. The measures estimated in the analysis are tested for statistical significance using twin surrogates. We find, in accordance with earlier studies, the major drivers for global mean temperature to be greenhouse gases, ENSO, volcanic activity, and solar irradiance. We further uncover a feedback between temperature and ENSO. Our results demonstrate the need to involve multiple, delayed interactions within the drivers of temperature in order to develop a more thorough picture of global temperature variations.
Recurrence-plot-based time series analysis is widely used to study changes and transitions in the dynamics of a system or temporal deviations from its overall dynamical regime. However, most studies do not discuss the significance of the detected variations in the recurrence quantification measures. In this letter we propose a novel method to add a confidence measure to the recurrence quantification analysis. We show how this approach can be used to study significant changes in dynamical systems due to a change in control parameters, chaos-order as well as chaos-chaos transitions. Finally we study and discuss climate transitions by analysing a marine proxy record for past sea surface temperature. This paper is dedicated to the 25th anniversary of the introduction of recurrence plots.
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.
The EEG is one of the most commonly used tools in brain research. Though of high relevance in research, the data obtained is very noisy and nonstationary. In the present article we investigate the applicability of a nonlinear data analysis method, the recurrence quantification analysis (RQA), to Such data. The method solely rests on the natural property of recurrence which is a phenomenon inherent to complex systems, such as the brain. We show that this method is indeed suitable for the analysis of EEG data and that it might improve contemporary EEG analysis.
Purpose: We present a new morphometric measure of trabecular bone microarchitecture, called mean node strength (NdStr), which is part of a newly developed approach called long range nodestrut analysis. Our general aim is to describe and quantify the apparent "latticelike" microarchitecture of the trabecular bone network.
Methods: Similar in some ways to the topological node-strut analysis introduced by Garrahan et al. [J. Microsc. 142, 341-349 (1986)], our method is distinguished by an emphasis on long-range trabecular connectivity. Thus, while the topological classification of a pixel (after skeletonization) as a node, strut, or terminus, can be determined from the 3 x 3 neighborhood of that pixel, our method, which does not involve skeletonization, takes into account a much larger neighborhood. In addition, rather than giving a discrete classification of each pixel as a node, strut, or terminus, our method produces a continuous variable, node strength. The node strength is averaged over a region of interest to produce the mean node strength of the region.
Results: We have applied our long range node-strut analysis to a set of 26 high-resolution peripheral quantitative computed tomography (pQCT) axial images of human proximal tibiae acquired 17 mm below the tibial plateau. We found that NdStr has a strong positive correlation with trabecular volumetric bone mineral density (BMD). After an exponential transformation, we obtain a Pearson's correlation coefficient of r - 0.97. Qualitative comparison of images with similar BMD but with very different NdStr values suggests that the latter measure has successfully quantified the prevalence of the "latticelike" microarchitecture apparent in the image. Moreover, we found a strong correlation (r - 0.62) between NdStr and the conventional node-terminus ratio (Nd/Tm) of Garrahan et al. The Nd/Tm ratios were computed using traditional histomorphometry performed on bone biopsies obtained at the same location as the pQCT scans.
Conclusions: The newly introduced morphometric measure allows a quantitative assessment of the long-range connectivity of trabecular bone. One advantage of this method is that it is based on pQCT images that can be obtained noninvasively from patients, i.e., without having to obtain a bone biopsy from the patient.
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.
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.
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.
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.
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.
Identifying causal relations from observational data sets has posed great challenges in data-driven causality inference studies. One of the successful approaches to detect direct coupling in the information theory framework is transfer entropy. However, the core of entropy-based tools lies on the probability estimation of the underlying variables. Herewe propose a data-driven approach for causality inference that incorporates recurrence plot features into the framework of information theory. We define it as the recurrence measure of conditional dependence (RMCD), and we present some applications. The RMCD quantifies the causal dependence between two processes based on joint recurrence patterns between the past of the possible driver and present of the potentially driven, excepting the contribution of the contemporaneous past of the driven variable. Finally, it can unveil the time scale of the influence of the sea-surface temperature of the Pacific Ocean on the precipitation in the Amazonia during recent major droughts.
We present new measures of complexity and their application to event-related potential data. The new measures are based on structures of recurrence plots and makes the identification of chaos-chaos transitions possible. The application of these measures to data from single-trials of the Oddball experiment can identify laminar states therein. This offers a new way of analyzing event-related activity on a single-trial basis
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.
Recurrence plots exhibit line structures which represent typical behaviour of the investigated system. The local slope of these line structures is connected with a specific transformation of the time scales of different segments of the phase-space trajectory. This provides us a better understanding of the structures occurring in recurrence plots. The relationship between the time-scales and line structures are of practical importance in cross recurrence plots. Using this relationship within cross recurrence plots, the time-scales of differently sampled or time- transformed measurements can be adjusted. An application to geophysical measurements illustrates the capability of this method for the adjustment of time-scales in different measurements. (C) 2005 Elsevier B.V. All rights reserved
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.
Higher variability in rainfall and river discharge could be of major importance in landslide generation in the north-western Argentine Andes. Annual layered (varved) deposits of a landslide dammed lake in the Santa Maria Basin (26°S, 66°W) with an age of 30,000 14C years provide an archive of precipitation variability during this time. The comparison of these data with present-day rainfall observations tests the hypothesis that increased rainfall variability played a major role in landslide generation. A potential cause of such variability is the El Niño/ Southern Oscillation (ENSO). The causal link between ENSO and local rainfall is quantified by using a new method of nonlinear data analysis, the quantitative analysis of cross recurrence plots (CRP). This method seeks similarities in the dynamics of two different processes, such as an ocean-atmosphere oscillation and local rainfall. Our analysis reveals significant similarities in the statistics of both modern and palaeo-precipitation data. The similarities in the data suggest that an ENSO-like influence on local rainfall was present at around 30,000 14C years ago. Increased rainfall, which was inferred from a lake balance modeling in a previous study, together with ENSO-like cyclicities could help to explain the clustering of landslides at around 30,000 14C years ago.
The presence of partially folded intermediates along the folding funnel of proteins has been suggested to be a signature of potentially aggregating systems. Many studies have concluded that metastable, highly flexible intermediates are the basic elements of the aggregation process. In a previous paper, we demonstrated how the choice between aggregation and folding behavior was influenced by hydrophobicity distribution patterning along the sequence, as quantified by recurrence quantification analysis (RQA) of the Myiazawa-Jernigan coded primary structures. In the present paper, we tried to unify the "partially folded intermediate" and "hydrophobicity/charge" models of protein aggregation verifying the ability of an empirical relation, developed for rationalizing the effect of different mutations on aggregation propensity of acyl-phosphatase and based on the combination of hydrophobicity RQA and charge descriptors, to discriminate in a statistically significant way two different protein populations: (a) proteins that fold by a process passing by partially folded intermediates and (b) proteins that do not present partially folded intermediates
A statistical model describing the propensity for protein aggregation is presented. Only amino-acid hydrophobicity values and calculated net charge are used for the model. The combined effects of hydrophobic patterns as computed by the signal analysis technique, recurrence quantification, plus calculated net charge were included in a function emphasizing the effect of singular hydrophobic patches which were found to be statistically significant for predicting aggregation propensity as quantified by fluorescence studies obtained from the literature. These results suggest preliminary evidence for a mesoscopic principle for protein folding/aggregation. (C) 2004 Elsevier B.V. All rights reserved
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.
Recurrence-plot-based measures of complexity and its application to heart-rate-variability data
(2002)
The knowledge of transitions between regular, laminar or chaotic behavior is essential to understand the underlying mechanisms behind complex systems. While several linear approaches are often insufficient to describe such processes, there are several nonlinear methods which however require rather long time observations. To overcome these difficulties, we propose measures of complexity based on vertical structures in recurrence plots and apply them to the logistic map as well as to heart rate variability data. For the logistic map these measures enable us not only to detect transitions between chaotic and periodic states, but also to identify laminar states, i.e. chaos-chaos transitions. The traditional recurrence quantification analysis fails to detect the latter transitions. Applying our new measures to the heart rate variability data, we are able to detect and quantify the laminar phases before a life-threatening cardiac arrhythmia occurs thereby facilitating a prediction of such an event. Our findings could be of importance for the therapy of malignant cardiac arrhythmias.