N. Marwan:
Recurrence Quantification Analysis
for Understanding Complex Systems,
CTCS Seminar series,
IIT Madras (India)/ Online,
Sep 25, 2023,
Talk invited.
» Abstract
Recurrence is a ubiquitous and fundamental feature in many real world processes. It is present at many scales in time and space, such as in celestial mechanics, alternating sediment layers, thermoacoustic oscillation, cardiac variability, and numerous other contexts. The study of recurrence properties (such as frequency analysis) can provide deeper insights into the underlying dynamical processes in general. A rather novel approach for the study of recurrences is the recurrence plot and its quantification, rooted in the theory of dynamical systems. In this lecture, the basic concepts and the major extensions applicable to various research questions are introduced and demonstrated. Discussed examples include the temporal change of recurrence properties for identification of regime shifts, for classification, and bivariate extensions for synchronization/coupling analysis.
N. Marwan:
Palaeoclimate research in RD4,
PIK cross-RD Paleoclimate and Long-term Climate Evolution Seminar,
Potsdam (Germany),
Sep 18, 2023,
Talk.
J. Wassmer, B. Merz, N. Marwan:
Resilience of transportation networks to road failures,
Dynamics Days Europe 2023,
Naples (Italy),
Sep 6, 2023,
Talk.
» Abstract
Damages to road infrastructure can cause disruptions in transportation and obstruct access to emergency services. With anthropogenic climate change increasing the probability of extreme weather events, such as floods or storms, the necessity for a resilient road infrastructure gets even more important.
In our research, we identify roads that play a crucial role in maintaining the stability of the transportation network. To this end we develop a framework that is built on a traffic-based centrality measure that can be interpreted as individual vehicles traversing the network. We then apply methods from the field of energy system analysis to derive the significance of each road segment in the network. The benefit of this framework is that it exclusively depends on openly accessible data sources such as OpenStreetMap, hence making it straightforward to apply and extend to different geographic locations of varying scales. As a case study, we analyse the impacts of the Ahr valley flood in Germany in 2021.
Our findings indicate that the road damages caused by the flood event led to an increase in the severity and frequency of congestion, result- ing in a deterioration of the accessibility to emergency services for a substantial portion of the population. We are further able to identify roads that are especially important for the resilience of the network. A broader application of our methodology can help decision makers to reduce costs by prioritising mitigation and reconstruction measures on important road sections.
M. Mannone, N. Marwan, A. Chella, V. Seidita, A. Giacometti, P. Fazio:
Entangled Gondolas. Design of Multi-Layer Networks of Quantum-Driven Robotic Swarms,
XVII International Workshop on Artificial Life and Evolutionary Computation (WIVACE 2023),
Venice (Italy),
Sep 6, 2023,
Talk.
» Abstract
Swarms of robots can be thought of as networks, using the tools from telecommunications and network theory. A recent study designed sets of aquatic swarms of robots to clean the canals of Venice, interacting with computers on gondolas. The interaction between gondolas is one level higher in the hierarchy of communication. In other studies, pairwise communications between the robots in robotic swarms have been modeled via quantum computing. Here, we first apply quantum computing to the telecommunication-based model of an aquatic robotic swarm. Then, we use multilayer networks to model interactions within the overall system. Finally, we apply quantum entanglement to formalize the interaction and synchronization between "heads" of the swarms, that is, between gondolas.
Our study can foster new strategies for search-and-rescue robotic-swarm missions, strengthening the connection between different areas of research in physics and engineering.
A. Syta, J. Czarnigowski, P. Jaklinski, N. Marwan:
Recurrence quantificators in misfire detection in a small aircraft engine,
10th International Symposium on Recurrence Plots,
Tsukuba (Japan),
Aug 28, 2023,
Talk.
» Abstract
Misfires in internal combustion engines are a frequent issue where one or more cylinders fail to ignite properly, leading to decreased engine performance, increased fuel consumption, and potential damage. A piston failure can be perceived as a disturbance in the repeatability of the engine operation, leading to changes in non-linearity. To study the dynamic behavior of the system over time, Recurrence Quantification Analysis (RQA) is employed in nonlinear time series analysis. The technique involves creating a recurrence plot (RP) from the time series data, which illustrates the temporal evolution of the system's states. RQA is capable of detecting changes in the system's behavior, such as transitions from regular to chaotic or from stable to unstable states. In the context of detecting piston failure, RQA can be used to identify patterns in the engine's vibration signals that are indicative of such failure. By placing sensors at different locations in the engine, vibrations can be recorded corresponding to separate engine states, including all cylinders working correctly and one of the cylinders being switched off, at various engine speeds. The RQA indices can then be used as non-linear features to classify the engine condition utilizing a linear model. Determining the RQA statistics on component signals with frequencies centered around the dominant ones increases the dimension of features and leads to higher accuracy in damage detection and identification of a cylinder with a misfire.
M. R. Sales, M. Mugnaine, J. D. Szezech Jr., R. L. Viana, I. L. Caldas, N. Marwan, J. Kurths:
Characterization of stickiness in quasi-integrable Hamiltonian systems by an entropy-based measure of the recurrence plots,
10th International Symposium on Recurrence Plots,
Tsukuba (Japan),
Aug 28, 2023,
Talk.
» Abstract
The stickiness effect is a fundamental feature of quasi-integrable Hamiltonian systems, characterized by the long time spend by a chaotic orbit when near enough a periodic island. We propose the use of an entropy-based measure of the recurrence plots (RPs), namely, the entropy of the distribution of the recurrence times (estimated from the RP), to characterize the dynamics of a typical quasi-integrable Hamiltonian system with coexisting regular and chaotic regions, the Chirikov-Taylor standard map. We show that the recurrence time entropy (RTE) is positively correlated to the largest Lyapunov exponent with a high correlation coefficient. We obtain a multi-modal distribution of the finite-time RTE and find that each mode corresponds to the motion around islands of different hierarchical levels.
B. G. Straiotto, N. Marwan, P. J. Seeley:
Exploring synchronisation in lower limb coordination in a rhythmic body movement: A quantitative analysis,
10th International Symposium on Recurrence Plots,
Tsukuba (Japan),
Aug 28, 2023,
Talk.
» Abstract
Studies of human movement often concern movement quality and that quality may be represented by the everyday term coordination. Research reports alternatively use terms such as correlation and synchronisation. We have developed a previous study of a martial arts movement pattern through study of synchronisation within and between the lower limbs. We explored synchronisation in taekwondo players who utilise repetitive backwards-forwards movements to mount attacks on their opponents and operate speedy retreats, movements that are developed in both training and competition. Eighteen players (nine elite and nine non-elite) performed backwards-forwards movements in a simulated training environment whilst their actions were recorded in detail via motion capture using multiple cameras. Recurrence analysis involved re-representing the time-dependent signals in multidimensional space and then characterising the revisits of a movement trajectory to different sub-regions of that space. The joint probability of recurrence index ($p_j$) was then calculated for centres of mass of limb segments (foot, shank, thigh) in relation to orthogonal movement coordinates (medio-lateral, anterior-posterior, vertical directions). Application of surrogation to the recurrence data indicated that derived $p_j$ values for elite and non-elite groups were deterministic in origin and not the result of data noise ($p < 0.01$). Interlimb pairwise segment relations yielded $p_j$ values in the range 0.23 to 0.29; intralimb relations in the range 0.24 to 0.40. Nonparametric statistical analysis combining Mann-Whitney and Kruskal-Wallis tests along with Bonferroni corrections directly indicated statistically significant differences between elite and non-elite groups for interlimb $p_j$ values and analogous differences for some comparisons for intralimb segment use ($0.05 > p > 0.01$). The potential of recurrence analysis for studies of limb segment synchronisation is revealed by this study. The method may be fruitfully extended in application not only to athletic movement but to transitions in coordination for general members of the public caused by ageing and pathology.
R. Pánis, K. Adámek, N. Marwan:
Averaged recurrence quantification analysis – Method omitting the recurrence threshold choice,
10th International Symposium on Recurrence Plots,
Tsukuba (Japan),
Aug 28, 2023,
Talk.
» Abstract
Recurrence quantification analysis (RQA) is a well established method of nonlinear data analysis. In this work we present a new strategy for an almost parameter-free RQA. The approach finally omits the choice of the threshold parameter by calculating the RQA measures for a range of thresholds (in fact recurrence rates). Specifically, we test the ability of the RQA measure determinism, to sort data with respect to their signal to noise ratios. We consider a periodic signal, simple chaotic logistic equation, and Lorenz system in thetested data set with different and even very small signal to noise ratios of lengths 102, 103, 104, and 105. To make the calculations possible a new effective algorithm was developed for streamlining of the numerical operations on Graphics Processing Unit (GPU).
N. Marwan, T. Braun, K. H. Kraemer, A. Banerjee, D. Eroglu:
Recurrence plots for analysing extreme events data,
10th International Symposium on Recurrence Plots,
Tsukuba (Japan),
Aug 28, 2023,
Talk.
» Abstract
The analysis of time series of extreme events is a challenging task. Many research questions, such as synchronisation analysis or power spectrum estimation, are challenging for linear tools. We demonstrate some recent extensions of the recurrence plot approach for various applications in the field of extreme events data. We demonstrate their potential for synchronisation analysis between signals of extreme events and signals with continuous and slower variations, for estimation of power spectra of spiky signals, and for analysing data with irregular sampling.
N. Marwan, T. Braun:
Power spectrum estimation for extreme events data,
10th International Congress on Industrial and Applied Mathematics (ICIAM),
Tokyo (Japan),
Aug 23, 2023,
» Talk (PDF, 5.27M)
.
» Abstract
The estimation of power spectral density (PSD) of time series is an important task in many quantitative scientific disciplines. However, the estimation of PSD from discrete data, such as extreme event series is challenging. We present a novel approach for the estimation of a PSD of discrete data. Combining the edit distance metric with the Wiener-Khinchin theorem provides a simple yet powerful PSD analysis for discrete time series (e.g., extreme events). This method works directly with the event time series without interpolation or transformation to continuous data. We demonstrate the method's potential on some prototypical examples and on event sequences of atmospheric rivers (AR), narrow filaments of extensive water vapor transport in the lower troposphere. Considering the spatial-temporal event series of ARs over Europe, we investigate the presence of a seasonal cycle as well as periodicities in the multi-annual range for specific regions, likely related to the North-Atlantic Oscillation (NAO).
M. L. Fischer, V. Foerster, F. Schaebitz, N. Marwan, S. Kaboth-Bahr, W. Schwanghart, M. H. Trauth:
A pan-African spatiotemporal framework of the past one million years using advanced multi-record time-series analysis,
XXI INQUA Conference,
Rome (Italy),
Jul 19, 2023,
Poster.
» Abstract
For several decades, eastern Africa was considered the origin ofH. sapiens, documented by the oldest fossil finds from Omo Kibish (233±22 ka BP) and Herto (160–154 ka BP), from where the species was thought to have spread across the rest of the continent and beyond. Recent finds of human fossils and related stone tools in several parts of Africa between roughly 315–75 ka ago, i.e. Jebel Irhoud, southern Africa, Arabia, and eastern Africa initiated a lively discussion of a multiregional model of the origin and development of H. sapiens. The chronology and diversity of human fossils and archaeological remains, associated with a pan-African cultural patchwork, are underpinned by the availability of suitable and connected environments offering enough resources for our species to survive and reproduce. As new paleoanthropological research expands into poorly understood regions of Africa, the key to understanding emerging patterns of mobility and dispersal within and out of Africa is strongly linked to accurately reconstructed climate and environmental conditions in time and space. Here, we aim to create a spatiotemporal paleoclimatic framework for testing current hypotheses about a multiregional origin of our species. To do so, we collect and review suitable climate archives that cover the time since the Mid-Pleistocene Transition, the Mid-Bruhnes Event, and the late Pleistocene. Prerequisites for site selection are: (1) an age model without major gaps, (2) proxy data with a sufficient temporal resolution, precision, and accuracy, (3) a good understanding of the mechanisms that are represented by the proxy data, and (4) together they offer good geographical coverage of Africa's most important climates. We compare records using correlation analysis, such as windowed Spearman correlation, principal component analysis, fast Fourier transformation, and wavelet-based cross-spectral analysis. Furthermore, we analyze long-term trends, shifts, and transition types, that may have provided a catalyst for evolutionary changes, cultural innovation, and expansion/ migration, using e.g. breakfit regression, running Mann-Whitney and Ansari-Bradley test, and recurrence-based transition tests, such as recurrence quantification and recurrence networks. Here, we show the first results of our experiments.
J. Klose, D. Scholz, M. Weber, H. Vonhof, B. Plessen, S. Breitenbach, N. Marwan:
Timing and progression of Dansgaard-Oeschger events in Central Europe based on three precisely dated speleothems from Bleßberg Cave, Germany,
XXI INQUA Conference,
Rome (Italy),
Jul 19, 2023,
Poster.
» Abstract
Speleothems can be dated with unprecedented precision using U-series disequilibrium methods and provide numerous climate proxies, such as stable oxygen (δ18O) and carbon isotopes (δ13C) or trace elements, resulting in long, sometimes continuous climate proxy records. Therefore, speleothems have great potential for reconstruction of past climate variability during Marine Isotope Stage (MIS) 3 and precise determination of the timing and duration of Dansgaard-Oeschger (D/O) events. While first discovered in Greenland ice cores, various speleothem records around the globe provided clear evidence for the supra-regional character of the D/O events. However, MIS 3 speleothem records from Central Europe are very limited. Here we present three spleothem (flowstone) MIS 3 records from Bleßberg Cave, Germany.
All flowstones show episodic growth with distinctive, partially very thin (<2 mm) growth phases, interrupted by visible hiatuses consisting of detrital material. Precise and accurate 230Th/U dating of the individual growth phases is challenging due to potential detrital contamination from these layers. Combining different sampling and analytical techniques, we were able to date even the thinnest growth layers with very high precision, i.e., 2σ-age uncertainties of at most a few hundred years.
The timing of the growth phases aligns with several D/O events, which have not been recorded in other Central European speleothems yet. The δ18O and δ13C records of all three flowstones are highly correlated which suggests a dominant process influencing both isotope systems. Comparison with the Sr and Mg records provides evidence for a strong influence of Prior Calcite Precipitation (PCP) in the aquifer above and inside the cave on the stable isotope and trace element signals. In addition, all proxy records are interpreted as evidence for past changes in precipitation and vegetation density and document a clear trend from more humid climate during early MIS 3 (ca. 57 – 50 ka) to less humid conditions during mid and late MIS 3 (ca. 45 – 30 ka).
Our multi-proxy approach thus allows us not only to precisely determine the timing, duration, and progression of several D/O events, but also to deepen our general understanding of climate variability during MIS 3 in Central Europe.
V. Skiba, C. Spötl, M. Trüssel, A. Schröder-Ritzrau, B. Plessen, N. Frank, R. Eichstädter, R. Tjallingii, N. Marwan, X. Zhang, J. Fohlmeister:
High-elevation speleothems suggest close coupling between North Atlantic millennial-scale variability and Alpine glacier dynamics during Marine Isotope Stage 8,
XXI INQUA Conference,
Rome (Italy),
Jul 18, 2023,
Poster.
» Abstract
Processes triggering abrupt climate transitions during glacial periods are still not fully understood. Most research has focused on the last glacial cycle, limiting our test bed for studying the occurrence and absence of millennial-scale variability and, thus, our understanding of these large-scale reorganisations of the climate system under different background conditions.
Here, we present new stalagmite oxygen and carbon data from high-elevation caves in central Switzerland covering the period from 300 to 200 ka. We demonstrate that millennial-scale variability recorded by these speleothems is representative of Northern Hemisphere interstadial-stadial variability. We use isotope-enabled fully-coupled ocean-atmosphere model simulations to show that the δ18O value of meteoric precipitation was higher by 1 ‰ during interstadials compared to stadials. This agrees with interstadial-stadial amplitudes of the last glacial cycle recorded by stalagmites from other caves in the Alps and is likely the result of North Atlantic seawater δ18O millennial-scale variability.
We find that the effect of prior carbonate precipitation (PCP) is superimposed on the meteoric δ18O signal, amplifying the isotope signal captured by Alpine speleothems on interstadial-stadial timescales. We propose that PCP variability provides a new proxy for milllennial-scale dynamics of warm-based paleoglaciers above these caves.
B. Keenan, J. Collins, B. Aichner, F. Schenk, S. Engels, C. Lane, W. Hoek, I. Neugebauer, T. Grunwald, F. Ott, M. Slowinski, S. Wulf, B. Goswami, N. Marwan, A. Brauer, D. Sachse:
Atmospheric blocking as a stabiliser during abrupt climate change in eastern Europe during the Last Deglaciation,
XXI INQUA Conference,
Rome (Italy),
Jul 14, 2023,
Talk.
» Abstract
Abrupt climate change has occurred frequently in Earth history, most notably during the termination of major glaciations in the Quaternary. Changes occurred over decadal timescales and destabilised or transformed landscapes and ecosystems. We use the Younger Dryas as a natural experiment to better understand the regional propagation of abrupt change. We applied hydrogen isotope analyses (δD) of plant wax lipid biomarkers as indicators of hydrological change and moisture origin from five lake records across Central Europe, namely lakes Meerfelder Maar, Steisslingen, Hämelsee, and Rehwiese palaeolake in Germany, and Lake Czechowskie (Trzechowskie palaeolake) in Poland. Using recurrence analysis, a method to detect and classify time-series and characterise dynamical regime shifts (tipping points), we identify the transition from warm to cold states, or from the relatively warm Allerød to the cold Younger Dryas, and return to the Holocene warm state. Further, isotopic gradients from the Allerød, Younger Dryas and early Holocene are compared with modern gradients, revealing spatial differences in the timing of the onset of the Younger Dryas event, in the magnitude and variability of change as well as the structure of the Younger Dryas event. We show that the pattern of δD responses during the YD were not consistent geographically, and that variation at different locations suggests greater climatic stability and smaller degree of change in the east. These spatiotemporal patterns are compared with modelling results and infer that atmospheric blocking over the Fennoscandian Ice Sheet (FIS) was likely the main driver of spatiotemporal patterns, creating a sustained high-pressure system over Fennoscandia deflecting the flow of westerly winds. This blocking effect only happens during summer as there are strong westerlies from September to April or May. Despite the abrupt climate change at the end of the Deglaciation, this resulted in a more stable climate during the YD event in eastern Europe, while western Europe was affected by major climate fluctuations during the second half of the YD because of weaker summer atmospheric blocking at the end of the YD meant more Atlantic inflow reaching MFM in the west but not in the east.
S. M. Vallejo-Bernal1, L. Luna, F. Wolf, N. Marwan, N. Boers, J. Kurths:
Atmospheric rivers are the drivers of precipitation-triggered landslides in western North America,
28th IUGG General Assembly,
Berlin (Germany),
Jul 14, 2023,
Talk.
» Abstract
Characterized by their specific geometry, atmospheric rivers (ARs) are narrow, long, and transient channels of intensive water vapor transport in the lower troposphere. They play an essential role in the water supply for precipitation in the mid-latitudes but can also trigger natural hazards such as floods and landslides by facilitating heavy precipitation events. In this study, we link the occurrence of landslides in western North America (NA) during the past decades to the precipitation triggered by land-falling ARs hitting the western coastline of the region. For this, we use a landslide inventory, rainfall estimates with a daily temporal resolution, and a catalog of land-falling ARs characterized in terms of strength and persistence based on the AR scale by Ralph et al., 2019. We employ two attribution models to relate rainfall to ARs and then landslides to AR-induced rainfall. Our results show that ARs precede between 60% and 100% of the landslides reported along the western coast of North America. Intense and persistent ARs are the most common precursors. As a further analysis, we study the synchronization pattern of landslides and ARs to determine if their association is unique and significant. In the coastal regions, the precedence relation of ARs leading to landslides is statistically significant. Further inland, landslides are less likely, but those that do occur are significantly correlated with very intense and persistent ARs. Understanding and revealing the impacts of ARs on landslides in western North America will lead to better forecasts and risk assessments of these natural hazards.
N. Marwan, K. H. Kraemer:
Recurrences of movement,
Workshop on Biomarkers arising from nonlinear analysis of movement variability,
Zurich (Switzerland),
July 5, 2023,
Lecture and workshop.
» Abstract
Hands-on workshop on recurrence analysis.
N. Marwan:
Recurrence analysis for complex systems,
Workshop on Biomarkers arising from nonlinear analysis of movement variability,
Zurich (Switzerland),
July 5, 2023,
Talk invited.
» Abstract
Recurrence is a ubiquitous and fundamental feature in many real-world processes. The study of recurrence properties (such as frequency analysis) can provide deeper insights into the dynamical processes in general. A rather novel approach for the study of recurrences is the framework based on recurrence plots and their quantification, rooted in the theory of dynamical systems. This talk will introduce the basic concepts and major extensions applicable to various research questions. Discussed examples include the temporal change of recurrence properties for identification of regime shifts; spatio-temporal recurrences for classification tasks; and bivariate extensions for synchronization/coupling analysis. Methodological and numerical challenges and pitfalls will also be discussed.
J. Zhang, J. Klose, M. Sierralta, S. Tsukamoto, D. Scholz, N. Marwan, S. Breitenbach:
Isothermal thermoluminescence (ITL) dating of a speleothem from Bleßberg Cave,
17th International Luminescence and Electron Spin Resonance Dating conference (LED2023),
Copenhagen (Denmark),
Jun 29, 2023,
Talk.
» Abstract
Their amenability to radiometric dating (U-series) makes speleothems (secondary cave carbonate deposits) a key archive of past climatic and environmental changes. However, incorporation of non-radiogenic thorium can hamper U-series dating, and samples older than ca. 600 ka remain out-of-reach for U-Th dating. Our aim is to develop isothermal thermoluminescence (ITL) dating as alternative approach for otherwise ‘undatable’ samples.
The calcite thermoluminescence (TL) signal (280 °C peak) saturates at much higher doses (saturation dose up to 5000 Gy) compared to quartz and feldspar, which shows great potential to extend the dating limit. However, spurious TL signal occurred at the high temperature range hindered its application. The conventional multiple-aliquot additive-dose (MAAD) protocol used for TL dating applies extrapolation for equivalent dose (De) estimation, which also has large error. Isothermal TL (ITL) dating with the single-aliquot regenerative-dose (SAR) protocol might be a promising way as it reduces the influence of the spurious TL signal, and it applies interpolation to obtain the De. However, this protocol has not been tested on samples with independent age control.
This study tests the ITL SAR dating protocol on a speleothem sample from Bleßberg cave, which has been accurately dated with 230Th/U (ca. 320–425 ka). ITL measurement at 235 °C for 200 °C can remove the 280 °C TL peak completely without TL contribution from higher temperature range. ITL De shows a plateau when the ITL temperature varies between 230 °C and 240 °C. Peak shifting and isothermal annealing tests indicate the 280 °C TL peak has a lifetime of tens of millions years at 10 °C, which is stable enough for the age range of this speleothem sample. The accurate alpha efficiency (α-value) and the U, Th distribution within the sample are measured to estimate the dose rate. The dose rate variation with time due to U-series disequilibrium is corrected for. The ITL ages are compared with the 230Th/U ages to evaluate the performance of the ITL dating protocol.
S. M. Vallejo-Bernal, T. Braun, N. Marwan, J. Kurths:
AR-tracks: A new comprehensive global catalog of atmospheric rivers,
EGU General Assembly,
Vienna (Austria),
Apr 25, 2023,
DOI:10.5194/egusphere-egu23-9251,
Talk.
» Abstract
Atmospheric rivers (ARs) are filaments of extensive water vapor transport in the lower troposphere. They play a crucial role in the global water cycle and are a main source of fresh water for the mid-latitudes. However, very intense and persistent ARs are important triggers of heavy rainfall events and have been associated with natural and economical damage. Further motivated by their high impacts, in the last decade occurrences of ARs have been intensively studied, detection algorithms have been developed, and multiple AR catalogs have been produced. As a common approach, the detection of ARs is based on localizing anomalous atmospheric transport of moisture, usually by setting an absolute threshold on vertically integrated vapor transport (IVT) and/or vertically integrated water vapor (IWV) fields. Behind this methodology, there is the implicit assumption of stationary atmospheric moisture levels, which is not necessarily true for long periods under the context of a warming atmosphere. Also, these thresholds have proven to vary regionally which results in often excluded low-level ARs.
Here, we introduce AR-tracks, a global, high-resolution catalog of atmospheric rivers that we have developed based on the Image-Processing-based Atmospheric River Tracking (IPART) algorithm, using IVT estimates of the ERA5 reanalysis data set. As opposed to conventional detection methods, IPART calculates anomalies of the IVT field at the synoptical spatiotemporal scale of ARs and is, therefore, free from magnitude thresholds and stationarity assumptions. The resulting catalog displays a list of AR events, with a spatial resolution of 0.75° x 0.75° and a temporal resolution of 6 hours, covering the period between 1979 and 2019. For each AR, we provide common parameters such as the time and location of the landfall, the respective IVT value, the area, the width, and the length of the AR. Moreover, we also track the contour and the axis of each AR, the position of the centroid, and the proportion of the AR that is located over ocean and land, and over the different continents.
To show the potential of this new catalog, we study the spatiotemporal variability of European ARs between 1979-2019, analyzing the robustness of our results for distinct parameter choices in the definition of AR-tracks. We also use a novel power spectral measure to identify periodic cycles in the occurrence of European ARs, revealing spatially heterogeneous seasonal and multi-annual periodicities. Finally, we discuss the role of land-falling ARs as a trigger of heavy precipitation events in the regional domain.
With the extensive data we provide in this new catalog, we aim at contributing to the further understanding of the role of ARs in global climate dynamics, as long-lived ARs having cross-continent tracks can be reliably traced through their tropical/subtropical origins to high-latitude landfall, and novel topics such as inland penetration of ARs can be studied.
M. H. Trauth, A. Asrat, M. L. Fischer, P. O. Hopcroft, V. Foerster, S. Kaboth-Bahr, H. F. Lamb, N. Marwan, M. A. Maslin, F. Schäbitz, P. J. Valdes:
Early Warning Signals for the Termination of the African Humid Period(s),
EGU General Assembly,
Vienna (Austria),
Apr 24, 2023,
DOI:10.5194/egusphere-egu23-5277,
Talk.
» Abstract
The study of the mid-Holocene climate tipping point in tropical and subtropical Africa is the subject of current research, not only because there is a comparatively simple but nonlinear relationship between the change in cause (orbital forcing) and the accelerated response of the monsoon system, but also because the African monsoon is an example of a potentially positive evolution of living conditions for humans: modeling results suggest that the Sahel is expanding northward in the wake of human-induced recent global warming, with green belts spreading northward. New literature distinguishes tipping elements such as the African monsoon according to the nature of the cause and the response of the climate system. Research here focuses primarily on tipping points of the type, which is characterized by a critical slowing down and a decreasing recovery from perturbations. The African monsoon, on the other hand, is an example of the tipping point of the type, which is characterized by flickering before the transition. The two types also differ in the nature of their early warning signals (EWS). These EWS are increasingly becoming the focus of research, as they are particularly important for predicting possible tipping of climate in the future of our planet. For the African monsoon system, flickering between two stable states near the transition has been predicted by modeling, but has not yet been demonstrated on paleoclimate time series.
The paleoenvironmental record from the Chew Bahir Basin in the southern Ethiopian Rift, which documents the climate history of eastern Africa of the past 620 ka with decadal resolution in some parts provides the possibility to examine the termination of the African Humid Period (AHP, 15–5 kyr BP) with regard to the possible occurrence of EWS. Thanks to six well-dated short sediment cores (<17 m, <47 kyr BP) and two long cores ( 290 m, <620 ka BP) we can not only study the last climate transition at 5.5 kyr BP in detail, but also similar transitions including possible EWS long before the first occurrence of Homo sapiens at 318 ka BP on the African continent. The analysis of the Chew Bahir record reveals a rapid ( 880 yr) change of climate at 5.5 kyr BP in response to a relatively modest change in orbital forcing that appears to be typical of climate tipping points. If this is the case then 14 dry events at the end of the AHP and 7 wet events after the transition, each of them 20–80 yrs long and recurring every 160±40 yrs, could indeed indicate a pronounced flickering between wet and dry conditions at the end of the AHP, providing significant EWS of an imminent tipping point. Compared to the low-frequency cyclicity of climate variability before and after the termination of the AHP, the flickering occurs on time scales equivalent to a few human generations and it is very likely (albeit speculative) that people were conscious of these changes and adapted their lifestyles to the rapid changes in water and food availability.
J. Wassmer, B. Merz, N. Marwan:
Resilience of emergency infrastructure networks after flooding events,
EGU General Assembly,
Vienna (Austria),
Apr 26, 2023,
DOI:10.5194/egusphere-egu23-1383,
Talk.
» Abstract
Extreme weather events can drastically influence the dynamics and stability of networked infrastructure systems like transportation networks or power grids. Climate change is increasing the frequency of such events, making their impact on human society and ecosystems increasingly relevant. Prominent examples include damage of critical infrastructure caused by heavy rainfalls and landslides. The devastating floods that struck Germany’s Ahr valley in 2021 are yet another reminder of the threat posed by such extreme events. Due to washed-out roads and further severe infrastructure damages, critical bottlenecks effectively cut off a substantial share of the population from assistance, hampering or even impeding their rescue.
In this study, we investigate the impact of flood events on transportation networks where stability is particularly important in order to ensure the accessibility of emergency services. Local changes in the underlying network dynamics can affect the whole road network and, in the worst case, cause a total collapse of the system through cascading failures. Because of the severe consequences of cascading events, we aim to recognise such spreading processes at an early stage and, in a further step, be able to prevent them. To this end, we set up a gravity model of travel to simulate the changes of the traffic load after flooding events to identify vulnerabilities in the system. We further analyse how the accessibility of emergency services is affected and if the population can be effectively reached in time.
N. Marwan, T. Braun:
Power Spectrum Estimation for (Extreme) Events Data,
EGU General Assembly,
Vienna (Austria),
Apr 25, 2023,
DOI:10.5194/egusphere-egu23-7730,
» Talk (PDF, 6.39M)
.
» Abstract
The estimation of power spectral density (PSD) of time series is an important task in many quantitative scientific disciplines. However, the estimation of PSD from discrete data, such as extreme event series is challenging. We present a novel approach for the estimation of a PSD of discrete data. Combining the edit distance metric with the Wiener-Khinchin theorem provides a simple yet powerful PSD analysis for discrete time series (e.g., extreme events). This method works directly with the event time series without interpolation. We demonstrate the method's potential on some prototypical examples and on event sequences of atmospheric rivers (AR), narrow filaments of extensive water vapor transport in the lower troposphere. Considering the spatial-temporal event series of ARs over Europe, we investigate the presence of a seasonal cycle as well as periodicities in the multi-annual range for specific regions, likely related to the North-Atlantic Oscillation (NAO).
N. Marwan, T. Braun:
Power Spectrum Estimation for (Extreme) Events Data,
PIK cross-RD Climate and weather extremes seminar,
Potsdam (Germany),
Apr 20, 2023,
Talk.
» Abstract
The estimation of power spectral density (PSD) of time series is an important task in many quantitative scientific disciplines. However, the estimation of PSD from discrete data, such as extreme event series is challenging. We present a novel approach for the estimation of a PSD of discrete data. Combining the edit distance metric with the Wiener-Khinchin theorem provides a simple yet powerful PSD analysis for discrete time series (e.g., extreme events). This method works directly with the event time series without interpolation or transformation to continuous data. We demonstrate the method's potential on some prototypical examples and on event sequences of atmospheric rivers (AR), narrow filaments of extensive water vapor transport in the lower troposphere. Considering the spatial-temporal event series of ARs over Europe, we investigate the presence of a seasonal cycle as well as periodicities in the multi-annual range for specific regions, likely related to the North-Atlantic Oscillation (NAO).
S. Gupta, Z. Su, A. Banerjee, N. Boers, N. Marwan, L. Magnusson, C. Lopez, E. Hernandez-Garcia, F. Pappenberger, J. Kurths:
Spatial synchronization patterns of extreme rainfall and convection in the Asian Summer Monsoon region,
Conference on Nonlinear Data Analysis and Modeling: Advances, Applications, Perspectives,
PIK Potsdam (Germany),
Mar 17, 2023,
Poster.
» Abstract
A deeper knowledge about the spatially coherent patterns of extreme rainfall events in the South and East Asian regions is of utmost importance for substantially improving the forecasts of extreme rainfall as their agro-based economies predominantly rely on the monsoon. In our work, we use a combination of a nonlinear synchronization measure and complex networks to investigate the spatial characteristics of extreme rainfall synchronicity in the Asian Summer Monsoon (ASM) region and gain a comprehensive understanding of the intricate relationship between its Indian and East Asian counterparts. We identify two modes of synchronization between the Indian Summer Monsoon (ISM) and the East Asian Summer Monsoon (EASM) – a southern mode between the Arabian Sea and south-eastern China in June which relates the onset of monsoon in the two locations, and a northern mode between the core ISM zone and northern China which occurs in July. Thereafter, we determine the specific times of high extreme rainfall synchronization, and identify the distinctively different large-scale atmospheric circulation, convection and moisture transport patterns associated with each mode. Furthermore, we discover that the intraseasonal variability of the ISM-EASM interconnection may be influenced by the different modes of the tropical intraseasonal oscillation (ISO). Our findings show that certain phases of the Madden-Julian oscillation and the boreal summer ISO favour the synchronization of extreme rainfall events in the June-July-August season between ISM and EASM. The impact of El Nino-Southern Oscillation on the convective sources of the two monsoon subsystems, and thus their interannual variability is investigated.
M. Gelbrecht, K. H. Krämer, N. Marwan:
TreeEmbedding: Optimal state space reconstruction via Monte Carlo decision tree search,
Conference on Nonlinear Data Analysis and Modeling: Advances, Applications, Perspectives,
PIK Potsdam (Germany),
Mar 17, 2023,
Poster.
» Abstract
TreeEmbedding is a novel method for an optimal time delay state space reconstruction from uni- and multivariate time series. The embedding process is considered as a decision tree, in which each leaf corresponds to an embedding cycle and is subject to an evaluation through an 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 that minimises the objective function. The Monte Carlo approach aims to prevent getting stuck in a local minimum of the objective function and can be used in a modular way: Practitioners can choose suitable statistics for delay-preselection and the 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. To showcase the method, we demonstrate its improvements over the classical time delay embedding methods on various 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. The method is ready to use in the form of an accompanying Julia package TreeEmbedding.jl.
S. M. Vallejo-Bernal, T. Braun, N. Marwan, J. Kurths:
Synchronized heavy rainfall events in Europe: the role of atmospheric rivers,
Conference on Nonlinear Data Analysis and Modeling: Advances, Applications, Perspectives,
PIK Potsdam (Germany),
Mar 17, 2023,
Poster.
» Abstract
Atmospheric rivers (ARs) are channels of enhanced water vapor transport in the lower troposphere. They play a crucial role in the fresh water supply of Europe, contributing to up to 30% of the rainfall budget in some regions along the western coast. However, very intense and persistent ARs are important triggers of heavy rainfall events and have been associated with natural and economical damage. Here, we investigate the large-scale spatio-temporal synchronization patterns between heavy rainfall events and landfalling ARs over Europe, during the period from 1979 to 2019. For that, we employ ARtracks, a novel global catalog of ARs, and select the AR events whose footprint intercept Europe. Then, we use an AR-intensity scale to rank the ARs in terms of strength and persistence. Based on ERA5 daily precipitation estimates, we obtain binary time series indicating the absence or presence of heavy rainfall by thresholding the daily precipitation intensity at the 95th percentile. Subsequently, we utilize event synchronization incorporating varying delays to reveal the temporal evolution of spatial patterns of heavy rainfall events in the aftermath of land-falling ARs. Finally, using composites of integrated water vapor transport, geopotential height, upper-level meridional wind, and rainfall, we attribute the formation of the synchronization patterns to well-known atmospheric circulation configurations, depending on the intensity level of the ARs. Our results reveal the role of ARs in the distribution of heavy rainfall events over Europe and advance the understanding of inland heavy precipitation by revealing the characteristic circulation patterns and the main climatic drivers associated to the synchronization patterns.
T. Braun, S. M. Vallejo-Bernal, D. Traxl, N. Marwan, J. Kurths:
A spatio-temporal analysis of global atmospheric rivers,
Conference on Nonlinear Data Analysis and Modeling: Advances, Applications, Perspectives,
PIK Potsdam (Germany),
Mar 17, 2023,
Poster.
» Abstract
Atmospheric rivers (ARs) are narrow, transient corridors of extensive water vapor transport in the lower troposphere. The role ARs play in the global water cycle can be regarded as a double-edged sword: while low-intensity ARs provide vital supply of freshwater and are rarely associated with heavy precipitation events (HPEs), high-level ARs can cause detrimental impacts when they land-fall. Detection of ARs is based on localizing anomalous atmospheric transport of moisture. Many approaches define a threshold to identify local anomalies in integrated vapor transport (IVT) in order to obtain catalogues of ARs, effectively assuming stationary atmospheric moisture levels and often excluding low-level ARs.
Here, we employ an AR-detection framework (`ARtracks') based on global ERA5 reanalysis data that utilizes image processing techniques (using the IPART algorithm). Our approach allows us to study the spatio-temporal variability of globally distributed AR tracks and potential changes due to increasing atmospheric moisture levels on a warming planet. We implement a scale that characterizes ARs based on their strength and persistence, distinguishing between ARs with potentially beneficial and detrimental impacts. A recent study has demonstrated the scope of this categorized AR catalogue for the analysis of synchronization of ARs and HPEs in North America. We analyse the robustness of our results for distinct parameter choices in the definition of AR tracks. A novel power spectral measure for the analysis of event-like time series enables us to identify significant cycles in AR occurrence. Finally, we discuss the role of land-falling ARs as a trigger of HPEs on a global scale.
T. Haselhoff, T. Braun, N. Marwan, S. Moebus:
Complex networks for the urban acoustic environment,
Conference on Nonlinear Data Analysis and Modeling: Advances, Applications, Perspectives,
PIK Potsdam (Germany),
Mar 17, 2023,
Poster.
» Abstract
The urban acoustic environment (AE) plays an underestimated role in the daily life of residents inhabiting metropolitan regions. The urban AE contains valuable information on complex sub-systems of urban areas, such as traffic, infrastructure and biodiversity. Associations between noise exposure and the mental or physical health of urban residents are an important subject of ongoing research. Despite the extensive information that is recorded by modern acoustic sensors, few approaches are designed to capture the rich complexity embedded in the time-frequency domain of the urban AE. The decreasing costs of acoustic sensors and rapid growth of storage space and computational power have led to an increase of acoustical data to be processed. Quantitative methods need to account for this complexity, while effectively reducing the high dimensionality of terabytes of audio data.
We take this as an opportunity to introduce complex networks to the field of urban acoustics. We use one of the world's most extensive longitudinal audio datasets from the SALVE study to systematically characterize the urban AE. SALVE is an ongoing study since 2019, in which 3-min acoustic recordings are made twice per hour at 23 locations in Bochum, Germany. The recorded acoustic samples exhibit a clear diel cycle and reveal site-dependent communities of interlinked frequencies. We demonstrate the utility of frequency-correlation matrices (FCMs) to effectively capture these communities. Based on these results, we construct (functional) networks of day time-specific 3-min audio recordings from 05.2019 to 03.2020 (n = 319,385 = 665 days). We show that the average shortest path length of an acoustic frequency network informs on site- and time-specific distinctiveness of frequency dynamics in the urban AE. To validate our findings, we use the land use mix around each site as a proxy for the AE, as the acoustic environment is heavily impacted by the built environment. The proposed method enables us to clearly identify 4-5 clusters of distinct urban AEs based on hourly variations in the distinctiveness of frequency dynamics. Our results indicate that complex networks represent a promising approach to analyse large-scale audio data and help to understand the time-frequency domain of the urban acoustic environment.
J. Wassmer, B. Merz, N. Marwan:
Resilience of emergency infrastructure networks after flooding events,
Conference on Nonlinear Data Analysis and Modeling: Advances, Applications, Perspectives,
PIK Potsdam (Germany),
Mar 17, 2023,
Poster.
» Abstract
Extreme weather events can drastically influence the dynamics and stability of networked infrastructure systems like transportation networks or power grids. Climate change is increasing the frequency of such events, making their impact on human society and ecosystems increasingly relevant. Prominent examples include damage of critical infrastructure caused by heavy rainfalls and landslides. The devastating floods that struck Germany’s Ahr valley in 2021 are yet another reminder of the threat posed by such extreme events. Due to washed-out roads and further severe infrastructure damages, critical bottlenecks effectively cut off a substantial share of the population from assistance, hampering or even impeding their rescue.
In this study, we investigate the impact of flood events on transportation networks where stability is particularly important in order to ensure the accessibility of emergency services. Local changes in the underlying network dynamics can affect the whole road network and, in the worst case, cause a total collapse of the system through cascading failures. Because of the severe consequences of cascading events, we aim to recognise such spreading processes at an early stage and, in a further step, be able to prevent them. To this end, we set up a gravity model of travel to simulate the changes of the traffic load after flooding events to identify vulnerabilities in the system. We further analyse how the accessibility of emergency services is affected and if the population can be effectively reached in time.
N. Antary, N. Marwan:
Interpolation effects an RQA measures,
Conference on Nonlinear Data Analysis and Modeling: Advances, Applications, Perspectives,
PIK Potsdam (Germany),
Mar 17, 2023,
Poster.
» Abstract
The recurrence plot and recurrence quantification analysis (RQA) are well established methods for the analysis of data from complex systems. They provide import insides about the nature of the dynamics, periodicity, regime changes, and many more. This method is used in different fields of research like finance, engineering, life and earth science. In order to use this method the data has usually to be uniformly sampled. This poses a difficulty for data, which is taken from palaeoclimate archives like sediment cores or stalagmite. One frequently used solution is interpolation to generate uniform time series. However, this prepossessing changes the RQA measures like DET, LAM, or the average line length. Using auto-regression processes, we systematically analyse how these measures increase when interpolating the data. For other systems which show a smoother behavior there is only an effect if the interpolation takes place on a time scale close to the characteristic timescale of the system, like the period lengths. For the Roessler system, the RQA measures decrease when approaching this timescale and show a very irregular behavior below. For real data, we show that the course of the DET measure strongly depends on the choice of interpolation.
M. R. Sales, M. Mugnaine, J. D. Szezech Jr., R. L. Viana, I. L. Caldas, N. Marwan, J. Kurths:
Characterizing stickiness using recurrence time entropy,
Conference on Nonlinear Data Analysis and Modeling: Advances, Applications, Perspectives,
PIK Potsdam (Germany),
Mar 17, 2023,
Poster.
» Abstract
The stickiness effect is a fundamental characteristic of quasi-integrable Hamiltonian systems. We propose the use of an entropy-based measure of recurrence plots (RP), namely, the entropy of the distribution of the recurrence times (estimated from the RP), to characterize the dynamics of a typical quasi-integrable Hamiltonian system with coexisting regular and chaotic regions. We show that the recurrence time entropy (RTE) is positively correlated to the largest Lyapunov exponent, with a high correlation coefficient. We obtain a multi-modal distribution of the finite-time RTE and show that each mode corresponds to the motion around islands of different hierarchical levels.
C. Özdes, D. Eroglu, N. Marwan, T. Braun:
Multi-stable synchronization patterns and switching dynamics of paleoclimate networks,
Conference on Nonlinear Data Analysis and Modeling: Advances, Applications, Perspectives,
PIK Potsdam (Germany),
Mar 17, 2023,
Talk.
» Abstract
To improve our understanding of climate dynamics, we first need to deeply understand the climate’s past if we hope to mitigate and adapt to oncoming critical climate change. Understanding the past climate dynamics depends on the interpretation of paleo proxies. Blending dynamical systems theory, recurrence theorem, multi-stability, and synchronization with complex networks theory and machine learning techniques have become instrumental for a more profound understanding of climate dynamics in the last few decades. However, these techniques are not directly applicable to paleoclimate research since the proxy data is subject to different distortions. The paleoclimate proxy measurements carry uncertainty in nominal and temporal dimensions, and also the choice of proxy and varying effects of local and global interactions matter.
Paleoclimate proxies typically represent the climate dynamics of large spatial regions and long periods. Furthermore, the proxies contain many switching transitions between droughts and wet seasons, showing that paleoclimate dynamics have multi-stability. To mimic paleoclimate dynamics, we introduce a multi-layer network model of coupled chaotic maps where multiple chimera configurations of synchronized subsystems co-exist as stable states. This multi-stable system goes through a series of critical transitions into another stable state through noise induction. We collect only the mean field of the state variables from each layer to imitate the spatial sparsity of paleoclimate measurements. Using this limited information, we developed a methodology to reconstruct paleoclimate networks and identify the critical switching of dynamical patterns.
Our paleoclimate network approach pivots around the recurrent property of climate system states. After suitable transformations, recurrence quantification analyses (RQA) of proxy series are shown to be robust indicators of the dynamical properties of represented dynamics in the form of time series. We construct a functional network from these series with nodes representing proxy sources using the time evolution of individual series. This allows us to classify the system state with respect to the visible relational dynamics between nodes. We also extended our studies to real paleoclimate datasets around Northern Africa and found the dominant dynamical patterns associated with known periods.
N. Marwan, T. Braun, K. H. Kraemer, A. Banerjee, D. Eroglu:
Recurrence plots for analysing extreme events data,
Conference on Nonlinear Data Analysis and Modeling: Advances, Applications, Perspectives,
PIK Potsdam (Germany),
Mar 17, 2023,
Talk.
» Abstract
The analysis of time series of extreme events is a challenging task. Many research questions, such as synchronisation analysis or power spectrum estimation, are challenging for linear tools. We demonstrate some recent extensions of the recurrence plot approach for various applications in the field of extreme events data. We demonstrate their potential for synchronisation analysis between signals of extreme events and signals with continuous and slower variations, for estimation of power spectra of spiky signals, and for analysing data with irregular sampling.
S. De, S. Gupta, V. R. Unni, R. Ravindran, P. Kasthuri, N. Marwan, J. Kurths, R. I. Sujith:
Implications of a Complex Network-Based Approach to the Analysis of Cyclone Merger,
Conference on Nonlinear Systems & Dynamics,
IISER Pune (India),
Dec 17, 2022,
Talk.
» Abstract
When cyclones are formed in close proximity, they can interact. Such an interaction is termed as “Fujiwhara effect” [1]. Due to this effect, the mutual distance between the cyclones decreases, which triggers a variety of interactions such as elastic interaction, partial merger, and complete merger. However, the interaction between the cyclones leading to a completer merger is a rare event in nature. The complete merger between cyclones can result in a more intense, long-lived cyclone. However, understanding the dynamics of cyclonic interactions is presently challenging for weather forecasters, making the prediction of a cyclone merger more difficult. The main reason attributed to the prediction inaccuracy is that, to date, cyclone forecasting models have not completely incorporated the Fujiwhara effect due to a lack of knowledge [2,3]. As a consequence, inaccurate cyclone merger predictions may result in substantial economic losses and fatalities. Hence, we require a method that can enable us to obtain profound insights into the dynamics of cyclone interaction that leads to a complete cyclone merger.
S. Kulkarni, U. Öztürk, N. Marwan, J. Kurths, B. Merz, A. Agarwal:
Complex networks in hydrologic sciences,
AGU 2022 Fall Meeting,
New Orleans (USA),
December 12–16, 2022,
Talk.
» Abstract
Complex network science is a stemming interdisciplinary field of research spreading to diverse branches of sciences such as physics, engineering, social science, and Earth Sciences. The complex systems can be represented as graphs with individual components called nodes and the links representing the interaction between nodes. Regardless of their physical nature, complex networks of different systems exhibit common structural properties that distinguish them from purely random graphs. The application of complex networks in hydrology and water resources management is in its infancy but overgrowing. However, this innovative approach has already led to important insights in hydrology. In this review, we first provide a comprehensive overview of the multiple aspects of complex networks and their measures; further, we summarize applications of complex networks in water science to offer a current picture of state-of-the-art; and lastly, we highlight arising open problems and new directions. Our work, with the help of examples, advocates that complex network science can be a generic theory to understand different hydrologic systems.
A. Manapat, J. L. Oster, F. Lechleitner, H. Cheng, J. F. Adkins, S. M. Bernasconi, W. D. Sharp, N. Marwan, B. Plessen, S. F. M. Breitenbach:
Stalagmite record of Indian Summer Monsoon variability during Marine Isotope Stage 3,
AGU 2022 Fall Meeting,
New Orleans (USA),
December 12–16, 2022,
Poster.
» Abstract
Variation in the strength of the Indian Summer Monsoon (ISM) affects the food security and livelihood of around one-third of the world’s population. However, many of the causes of variation in ISM strength remain poorly understood. Variation in the stable oxygen and carbon isotope ratios (d18O and d13C) in speleothems have been used as proxies to understand pre-instrumental variation in ISM intensity- with δ18O variation linked to changes in precipitation dynamics and δ13C reflecting changes in the local rainfall amounts.
We present a stable isotope record from MAW-3, a stalagmite from Mawmluh cave, northeast India. This site receives between 70-80% of its annual rainfall between June and September. Precipitation, cave drip monitoring, and studies of modern stalagmites from this site indicate that stalagmite d18O reflects changes in ISM strength linked to large-scale atmospheric dynamics. We discuss the interval of growth for stalagmite MAW-3 that occurred between 44 ka BP and 28 ka BP, corresponding to Marine Isotope Stage 3 (MIS3). The resultant δ18O record displays significant millennial-scale oscillations of between 1.5-3.0 ‰, possibly corresponding to Dansgaard-Oeschger (D-O) events 5-8. Continuous Wavelet Transform (CWT) indicates a 1350-year cycle in δ18O. Stalagmite δ13C also decreases by 2.0-4.0 ‰ during D-O interstadials, likely due to reduced prior carbonate precipitation (PCP) during periods of higher rainfall. The MAW3 δ18O shifts to more negative values during the D-O warmings (interstadials) noted in the North Greenland Ice Core Project (NGRIP) oxygen isotope record. However, the δ18O variation in MAW-3 is much more gradual than that of NGRIP, indicating differences in how these D-O events affect climate variation outside of the North Atlantic region. Variations in δ18O in MAW3 are also in phase with δ18O variations in the Hulu Cave record over this interval indicating synchronous variability between the ISM and East Asian Monsoon over D-O events 5-8. The MAW3 record further strengthens evidence for climate teleconnections between the North Atlantic and the Indian and East Asian monsoon systems. Reconstructing such variation is important in understanding how the monsoons may change in a warming planet.
N. Marwan:
Methods for Extreme Events Time Series,
PIK cross-RD Climate and weather extremes seminar,
Potsdam (Germany),
Dec 8, 2022,
Talk.
N. Marwan:
Investigating palaeoclimate conditions with nonlinear time series analysis,
MARUM research seminar,
MARUM University of Bremen (Germany),
Dec 5, 2022,
Talk invited.
» Abstract
The study of palaeoclimate data is related with specific challenges, such as nonstationarities, nonlinear feedbacks, irregular sampling, or different kinds of uncertainties. Recurrence analysis and complex networks are concepts based on nonlinear dynamics and complex systems science that can help to identify regime transitions and couplings in palaeoclimate data. I demonstrate their potential for studying variations and couplings for selected palaeoclimate research questions.
N. Marwan:
Ways to Quantitative Recurrence Plot Analysis,
Complexity Lab Seminar,
Rochester Institute of Technology (US),
Nov 15, 2022,
Lecture.
» Abstract
Recurrence is a ubiquitous and fundamental feature in many real world processes. Recurrence plots are versatile tools for studying such phenomena. The lecture introduces the basic concepts and major extensions of quantifying recurrence plots. Discussed examples include the temporal change of recurrence properties for identification of regime shifts, methodological and numerical challenges, as well as potential pitfalls.
V. Skiba, M. Trüssel, B. Plessen, C. Spötl, R. Eichstädter, A. Schröder-Ritzrau, T. Braun, T. Mitsui, N. Frank, N. Boers, N. Marwan, J. Fohlmeister, R. Tjallingii, X. Zhang:
On the forcing of glacial abrupt climate transitions of the last 300,000 years,
DEUQUA-Tagung 2022,
Potsdam (Germany),
Sep 26, 2022,
Talk.
» Abstract
Abrupt stadial-interstadial transitions, are a prominent feature of the last glacial as recorded in Greenland ice core records (Dansgaard-Oeschger events). Event abruptness and presence of statistical early warning signals before these transitions indicate that they involve crossing of a tipping point of the climate system. However, only little information is available for periods before the last glacial period as Greenland ice cores and many other high-resolution records do not extent beyond the last glacial cycle. Given the lack of understanding of the triggering mechanisms responsible for glacial abrupt climate transitions with palaeoclimate data from the last glacial, it is essential to investigate this phenomenon during earlier glacial periods.
Here, we present a new highly resolved, precisely U-Th-dated speleothem oxygen isotope record from the Northern European Alps for the penultimate glacial (MIS7-MIS8), a region which has been shown to record similar climate variability as Greenland ice core records. Together with previously obtained speleothem data from this cave site for MIS5-MIS7 and Greenland ice core data (NGRIP, MIS1-4), we investigate background climate conditions which favour occurrence of abrupt climate transitions using regression analysis. Besides intermediate background conditions (sea level, CO2 and CH4) and low precession, we find either relatively low or high obliquity to favour glacial abrupt climate transitions, perhaps depending on the initial mode of the Atlantic Meridional Overturning Circulation before these occurrences.
T. Braun, N. Marwan:
A recurrence flow based approach to state space reconstruction,
Dynamics Days Europe 2022,
University of Aberdeen (UK),
Aug 23, 2022,
Talk.
» Abstract
In the study of nonlinear observational time series, reconstructing the system’s state space via time-delay embedding represents the basis for many widely-used analyses. Recurrence plots indicate the appropriateness of the underlying embedding parameters by the presence of well-formed diagonal lines that represent the predictability of the system's evolution. However, an approach that systematically exploits this information for optimal state space reconstruction is so far missing. In this talk, we propose a recurrence based framework for state space reconstruction. The concept is based on a novel recurrence quantification measure that captures how well a fictive fluid can permeate an RP diagonally, the recurrence flow. The recurrence flow can be regarded as a nonlinear dependence measure that quantifies the relationship between multiple time series based on the predictability of their joint evolution. We demonstrate the effectiveness of the proposed method in detecting nonlinear multiscale relations and informing on the choice of optimal embedding parameters for complex real-world time series.
N. Marwan, K. H. Kraemer:
Recent Exciting Developments in Recurrence Plot Analysis,
Dynamics Days Europe 2022,
University of Aberdeen (UK),
Aug 23, 2022,
Talk.
» Abstract
The last decade has witnessed a number of important and exciting developments that had been achieved for improving recurrence plot based data analysis and to widen its application potential. I will give a brief overview about important and innovative developments, such as computational improvements, alternative recurrence definitions (event-like, multiscale, heterogeneous, and spatio-temporal recurrences) and ideas for parameter selection, theoretical considerations of recurrence quantification measures, new recurrence quantifiers (e.g., for transition detection and causality detection), and correction schemes. Moreover, new perspectives have recently been opened by combining recurrence plots with machine learning.
N. Marwan:
Reconstructing Complex Networks from Data,
(Virtual) Workshop on the Application of Complex Networks to Fluid Mechanics,
IIT Madras (India),
Aug 15, 2022,
Talk.
» Abstract
Complex networks provide an interesting tool to investigate spatio-temporal data. The first step is to reconstruct a (functional) network from data. I will show different reconstruction approaches depending on the research question and the nature of the data. The procedure is illustrated with applications on climate data.
S. Gupta, Z. Su, N. Boers, J. Kurths, N. Marwan, F. Pappenberger:
Spatial Synchronization Patterns of Extreme Rainfall Events in the Asian Summer Monsoon Region,
AOGS2022 Virtual 19th Annual Meeting,
(online meeting),
Aug 1, 2022,
Talk.
» Abstract
A deeper knowledge about the spatially coherent patterns of extreme rainfall events in the South and East Asian regions is of utmost importance for substantially improving the forecasts of extreme rainfall as their agro-based economies predominantly rely on the monsoon. In our work, we use a combination of a nonlinear synchronization measure and complex networks to investigate the spatial characteristics of extreme rainfall synchronicity in the Asian Summer Monsoon (ASM) region and gain a comprehensive understanding of the intricate relationship between its Indian and East Asian counterparts. We identify two modes of synchronization between the Indian Summer Monsoon (ISM) and the East Asian Summer Monsoon (EASM) – a southern mode between the Arabian Sea and south-eastern China in June which relates the onset of monsoon in the two locations, and a northern mode between the core ISM zone and northern China which occurs in July. Thereafter, we determine the specific times of high extreme rainfall synchronization, and identify the distinctively different large-scale atmospheric circulation, convection and moisture transport patterns associated with each mode. Furthermore, we discover that the intraseasonal variability of the ISM-EASM interconnection may be influenced by the different modes of the tropical intraseasonal oscillation (ISO). Our findings show that certain phases of the Madden-Julian oscillation and the boreal summer ISO favour the synchronization of extreme rainfall events in the June-July-August season between ISM and EASM. This work is funded by the CAFE project which has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 813844.
J. Klose, M. Weber, H. Vonhof, B. Plessen, S. F. M. Breitenbach, N. Marwan, D. Scholz:
Timing of Dansgaard-Oeschger events in Central Europe based on three precisely dated speleothems from Bleßberg Cave, Germany,
Climate Change, The Karst Record IX (KR9),
Innsbruck (Austria),
July 18, 2022,
Talk.
» Abstract
The last glacial period and especially Marine Isotope stage 3 (MIS 3, ca. 57 - 27 ka) was characterized by various climate oscillations (i.e., rapid increases in temperature, followed by a gradual cooling, the Dansgaard-Oeschger (D/O) events), which were first discovered in Greenland ice cores. Although their causes are still not fully understood, clear evidence for their supra-regional character was found in various climate records around the globe. However, European speleothem samples, which grew during MIS 3, are limited and mainly restricted to alpine regions, where glacier meltwater enabled speleothem growth, and to south/south-western parts of Europe characterised by a generally warmer climate. This led to the opinion that it was too cold and/or too dry in central Europe for speleothem growth. Here we present three speleothem (flowstone) records from Bleßberg Cave, Germany, which grew during MIS 3.
All flowstones show episodical growth patterns with distinctive, thin growth phases. Potential contamination deriving form detrital material deposited during hiatuses between individual growth phases, open-system behaviour around the hiatuses and the limited thickness of the growth layers are the biggest challenges during sampling for 230Th/U dating. By combination of different sampling techniques (i.e., laser ablation and micro-milling) in addition to the common approach of handheld drilling and due to the relatively high 238U concentration of the samples (approx. 0.4 - 1 μg/g), we were able to date even the thinnest growth layers (< 2 mm) of the Bleßberg flowstones with a very high precision (i.e., with 2σ-age uncertainties of a few hundred years or even lower).
The timing of the growth phases of the Bleßberg flowstones correlates with several D/O events recorded in the Greenland ice cores. This proves that at least some phases of MIS 3 had favourable climate conditions for speleothem growth in Central Europe. In addition, the analysis of the stable oxygen and carbon isotopes (δ18O and δ13C) for all three flowstones revealed several D/O events, which have not been recorded in any other speleothem from central Europe so far. This will enhance our understanding of climate variability during MIS 3 and specific D/O events in central Europe.
N. Marwan:
Investigating palaeoclimate conditions with nonlinear time series analysis,
Colloquium on Complex and Biological Systems, University of Potsdam,
Potsdam (Germany),
Jun 24, 2022,
Talk.
» Abstract
The study of palaeoclimate data is related with specific challenges, such as irregular sampling or different kinds of uncertainties. Recurrence analysis and complex networks are concepts based on nonlinear dynamics and complex systems science that can help to identify regime transitions and couplings. I demonstrate their potential for studying variations and couplings in the palaeoclimate Monsoon system.
N. Marwan:
Das Sägistal - Alpines Karstgebiet im Berner Oberland,
60. Jahrestagung des VdHK,
Truckenthal (Germany),
Jun 17, 2022,
Talk.
N. Marwan:
Bleßberghöhle – Schatzkammer für die Wissenschaft,
VdHK-Symposium: Wissenschaft unter Tage – Höhlenforschung im Dialog,
Truckenthal (Germany),
Jun 15, 2022,
» Talk (PDF, 8.96M)
.
N. Marwan, J. F. Donges, R. V. Donner, D. Eroglu:
Integrative multivariate study of past African climate variability,
EGU General Assembly,
Vienna (Austria),
May 24, 2022,
DOI:10.5194/egusphere-egu22-6559,
» Talk (PDF, 7.01M)
.
» Abstract
Based on a set of various marine palaeoclimate proxy records, we investigate African climate variations during the past 5 million years. We use a collection of modern approaches from non-linear time series analysis to identify and characterise dynamical regime shifts as changes in signal predictability, regularity, complexity, and higher-order stochastic properties such as multi-stability. We observe notable nonlinear transitions and important climate events in the African palaeoclimate, which can be attributed to phases of intensified Walker circulation, marine isotope stage M2, the onset of northern hemisphere glaciation, and the mid-Pleistocene transition, and relate them to variations of the Earth's orbital parameters.
A. Giesche, D. A. Hodell, C. A. Petrie, G. H. Haug, J. F. Adkins, B. Plessen, N. Marwan, H. J. Bradbury, A. Hartland, A. D. French, S. F. M. Breitenbach:
Northwest Indian stalagmite shows evidence for recurring summer and winter droughts after 4.2 ka BP,
EGU General Assembly,
Vienna (Austria),
May 24, 2022,
DOI:10.5194/egusphere-egu22-396,
Talk.
» Abstract
We reconstructed changes in summer and winter precipitation using a well-dated (±18 years 2σ error) speleothem spanning 4.2-3.1 ka BP from Dharamjali Cave in the central Himalaya. The record was sampled at a sub-annual resolution for a suite of trace elements, as well oxygen and carbon stable isotopes. Calcium isotopes at decadal resolution provide additional hydroclimatic evidence. This DHAR-1 stalagmite records a 230-year period of increased drought frequency in both the summer and winter seasons after 4.2 ka BP, with aridity events centered on 4.19, 4.11 and 4.02 ka BP each lasting between 25 and 90 years. The data after 3.97 ka BP support a recovery in summer monsoon rainfall, peaking around 3.7 ka BP. The significance of this new record includes remarkable coherence between the moisture proxies over 4.2-3.97 ka BP in a well-dated record, which provides confidence in the duration of droughts and timing of monsoon recovery. It also places seasonal climate variability on a timescale relevant to human decision-making, which is particularly significant for this region nearby the Indus River Basin. The Indus Civilization reached its urban apex by 4.2 ka BP, and archaeologists have documented a shift in settlement locations, population, health, and agricultural strategies thereafter for a period of several centuries. This stalagmite record provides valuable insights into seasonal precipitation availability during a critical climatic and cultural transition phase.
J. Wassmer, N. Marwan, B. Merz:
Impact of extreme events on topological robustness of infrastructure networks,
EGU General Assembly,
Vienna (Austria),
May 26, 2022,
DOI:10.5194/egusphere-egu22-2705,
Talk.
» Abstract
Climate change is increasing the frequency of extreme weather events such as floods, making their impact on human society and ecosystems increasingly relevant. Extreme weather events can drastically influence the dynamics and stability of networked infrastructure systems like transportation networks or power grids. Local changes in the dynamics can affect the whole network and, in the worst case, cause a total collapse of the system through cascading failures. Hence, methods are needed to understand and prevent such collapses.
In this project, we analyse the influences of flooding events on transportation networks using the Ahr valley flood of July 2021 in Germany as a case study. To this end, we set up a gravity model for road networks to compute the traffic loads. We use satellite data provided by EU Copernicus programme to access information about the state of the road network right after the flooding event. By removing flooded roads from the model, we can estimate the effect on the traffic load and identify secondarily affected roads. This approach enables us to identify and optimise critical links to ensure that affected areas are not isolated after extreme weather events and can receive disaster assistance from surrounding areas.
F. Wolf, S. M. Vallejo-Bernal, N. Boers, N. Marwan, D. Traxl, J. Kurths:
Spatio-temporal synchronization of heavy rainfall events triggered by atmospheric rivers in North America,
EGU General Assembly,
Vienna (Austria),
May 24, 2022,
DOI:10.5194/egusphere-egu22-2993,
» Talk (PDF, 758.19K)
.
» Abstract
Atmospheric rivers (ARs) are filaments of extensive water vapor transport in the lower troposphere. They are important triggers of heavy rainfall events, contributing to more than 50% of the rainfall sums in some regions along the western coast of North America. ARs play a crucial role in the distribution of water, but can also cause natural and economical damage by facilitating heavy rainfall. Here, we investigate the large-scale spatio-temporal synchronization patterns of heavy rainfall triggered by ARs over the western coast and the continental regions of North America.
For our work, we employ daily ERA5 rainfall estimates at a spatial resolution of 0.25°x0.25° latitude and longitude which we threshold at the 95th percentile to obtain binary time series indicating the absence or presence of heavy rainfall. Subsequently, we separate periods with ARs and periods without ARs and investigate the differing spatial synchronization pattern of heavy rainfall. To establish that our results are not dependent on the chosen AR catalog, this is conducted in two different ways: first based on a recently published catalog by Gershunov et al. (2017) , and second based on a catalog constructed using the IPART algorithm (Xu et al, 2020). For both approaches, we subsequently utilize event synchronization and a complex network framework to reveal distinct spatial patterns of heavy rainfall events for periods with and without active ARs. Using composites of upper-level meridional wind, we attribute the formation of the rainfall synchronization patterns to well-known atmospheric circulation configurations, whose intensity scales with the strength of the ARs. Furthermore, we demonstrate that enhanced AR activity is going in hand with a suppressed seasonal shift of the characteristic meridional wind pattern. To verify and illustrate how small changes of the high-level meridional wind affect the distribution of heavy rainfall, we, additionally, perform a case study focusing on the boreal winter.
Our results indicate the strong sensitivity of the intensity, location, frequency, and pattern of synchronized heavy rainfall events related to ARs to small changes in the large-scale circulation.
S. M. Vallejo-Bernal, F. Wolf, L. Luna, N. Boers, N. Marwan, J. Kurths:
Relationship between atmospheric rivers and landslides in western North America,
EGU General Assembly,
Vienna (Austria),
May 25, 2022,
DOI:10.5194/egusphere-egu22-8096,
» Talk (PDF, 1.79M)
.
» Abstract
In this study, we investigate the relationship between land-falling atmospheric rivers (ARs) and landslides in western North America. ARs are channels of enhanced water vapor flux in the atmosphere and play an essential role in the water supply for precipitation in the midlatitudes. However, they can also trigger natural hazards such as floods and landslides. Our objective is to determine if the occurrence of landslides in western North America can be attributed to ARs hitting the western coastline and causing rainfall at the locations of the landslides and to characterize the strength and persistence of the ARs that lead to landslides. To that aim, we use landslide records with daily temporal resolution along with daily rainfall estimates from the ERA5 reanalysis, for the period between 1996 and 2018. We propose and run two attribution models to relate landslides to rainfall and rainfall to ARs and subsequently verify statistically if there is a unique and significant association between the landslides and the ARs. Our results show that the majority of the landslides reported along the western coast of North America are preceded by an AR. In the coastal regions, ARs and landslides are significantly correlated. Further inland, landslides are less likely, but those that do occur are significantly correlated with very intense ARs. Understanding and revealing the impacts of ARs on landslides in western North America will lead to better forecasts and risk assessments of these natural hazards.
S. Gupta, Z. Su, N. Boers, J. Kurths, N. Marwan, F. Pappenberger:
Interrelation between the Indian and East Asian Summer Monsoon: A complex network-based approach,
EGU General Assembly,
Vienna (Austria),
May 23, 2022,
DOI:10.5194/egusphere-egu22-8626,
» Talk (PDF, 1.96M)
.
» Abstract
The Indian Summer Monsoon (ISM) and the East Asian Summer monsoon (EASM) are two integral components of the Asian Summer Monsoon system, largely influencing the agro-based economy of the densely populated southern and eastern parts of Asia. In our study, we use a complex network based approach to investigate the spatial coherence of extreme precipitation in the Asian Summer Monsoon region and gain a deep insight into the complex nature of the interaction between the ISM and the EASM. We identify two dominant modes of ISM-EASM interaction – (a) a southern mode connecting onset of the ISM over the Arabian Sea and southern India in June to the onset of Meiyu over south-eastern China, i.e., lower and middle reaches of the Yangtze river valley, and (b) a northern mode relating the occurrence and intensity of rainfall over the northern and central parts of India to that in northern China during July. Through determination of specific times of high synchronization of extreme precipitation, we distinctly identify the particular large-scale atmospheric circulation and moisture transport patterns associated with each mode. Thereafter, we investigate the role of the different components of the tropical intraseasonal oscillations, such as the Madden-Julian Oscillation and the boreal summer intraseasonal oscillation, in the intraseasonal variability of the relationship between the ISM and the EASM.
This work is funded by the CAFE project which has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 813844.
T. Braun, K. H. Kraemer, N. Marwan:
A Recurrence Flow based Approach to Attractor Reconstruction,
EGU General Assembly,
Vienna (Austria),
May 25, 2022,
DOI:10.5194/egusphere-egu22-9626,
Talk.
» Abstract
In the study of nonlinear observational time series, reconstructing the system’s state space represents the basis for many widely-used analyses. From the perspective of dynamical system’s theory, Taken’s theorem states that under benign conditions, the reconstructed state space preserves the most fundamental properties of the real, unknown system’s attractor. Through many applications, time delay embedding (TDE) has established itself as the most popular approach for state space reconstruction. However, standard TDE cannot account for multiscale properties of the system and many of the more sophisticated approaches either require heuristic choice for a high number of parameters, fail when the signals are corrupted by noise or obstruct analysis due to their very high complexity.
We present a novel semi-automated, recurrence based method for the problem of attractor reconstruction. The proposed method is based on recurrence plots (RPs), a computationally simple yet effective 2D-representation of a univariate time series. In a recent study, the quantification of RPs has been extended by transferring the well-known box-counting algorithm to recurrence analysis. We build on this novel formalism by introducing another box-counting measure that was originally put forward by B. Mandelbrot, namely succolarity. Succolarity quantifies how well a fluid can permeate a binary texture. We employ this measure by flooding a RP with a (fictional) fluid along its diagonals and computing succolarity as a measure of diagonal flow through the RP. Since a non-optimal choice of embedding parameters impedes the formation of diagonal lines in the RP and generally results in spurious patterns that block the fluid, the attractor reconstruction problem can be formulated as a maximization of diagonal recurrence flow.
The proposed state space reconstruction algorithm allows for non-uniform embedding delays to account for multiscale dynamics. It is conceptually and computationally simple and (nearly) parameter-free. Even in presence of moderate to high noise intensity, reliable results are obtained. We compare the method’s performance to existing techniques and showcase its effectiveness in applications to paradigmatic examples and nonlinear geoscientific time series.
V. Skiba, M. Trüssel, B. Plessen, C. Spötl, R. Eichstädter, A. Schröder-Ritzrau, T. Braun, T. Mitsui, N. Frank, N. Boers, N. Marwan, J. Fohlmeister:
Abrupt climate events recorded in speleothems from the ante penultimate glacial,
EGU General Assembly,
Vienna (Austria),
May 23, 2022,
DOI:10.5194/egusphere-egu22-11671,
» Talk (PDF, 1.04M)
.
» Abstract
Millennial-scale climate variability, especially abrupt stadial-interstadial transitions, are a prominent feature of the last glacial as recorded in Greenland ice core records (Dansgaard-Oeschger events). Event abruptness and presence of statistical early warning signals before these transitions indicate that they involve repeated crossing of a tipping point of the climate system. However, only little information is available for periods before the last glacial period as Greenland ice cores and many other high-resolution records do not extent beyond the last glacial cycle. Given the lack of understanding of the triggering mechanism responsible for glacial millennial-scale variability with palaeoclimate data from the last glacial, it is essential to investigate this phenomenon during earlier glacial periods.
Here, we present a new highly resolved, precisely U-Th-dated speleothem oxygen isotope record from the Northern European Alps, a region which has been previously shown to resemble the glacial millennial-scale climate variability obtained from Greenland ice core records very well. Our new data covers the time interval from the ante-penultimate glacial to the penultimate glacial (MIS8-MIS6) with a high degree of replication. For both glacial periods, we find phases of pronounced millennial-scale variability but also several, 10 ka long phases with the climate system being exclusively in stadial conditions. We compare our data with conceptual model results and investigate the occurrence and absence of abrupt climate transitions of the last 300,000 a.
T. Braun, C. N. Fernandez, D. Eroglu, A. Hartland, S. F. M. Breitenbach, N. Marwan:
Sampling rate-corrected time series analysis of irregularly sampled palaeoclimate proxy records,
6th PAGES Open Science Meeting,
(online meeting),
May 18, 2022,
Talk.
» Abstract
Irregular sampling remains a challenge in the analysis of time series from different palaeoclimate archives. Aside from rendering most standard time series analysis methods inapplicable, changes in the sampling rate of a record entail significant biases that can not be corrected by basic pre-processing procedures such as linear interpolation. Yet, sampling rates frequently show non-stationary characteristics, e.g. for speleothems, as they are coupled to environmental parameters via their growth rate.
Consequently, methods that account for continuous and abrupt changes of sampling resolution without introducing additional biases are required. In several applications, the edit–distance has proven to be an effective metric to quantitatively compare time series segments of unequal length by computing the cost of transforming one segment into the other. We demonstrate that transformation costs generally exhibit a non-trivial relationship with local sampling rate. If the sampling rate undergoes significant variations, this dependence rules out an unbiased comparison between different time episodes. We study the impact of this effect on recurrence quantification analysis, a framework that is well-suited for identifying regime shifts in nonlinear time series. A constrained randomization procedure is proposed as a bias correction for recurrence quantification analysis.
We demonstrate the effectiveness of the proposed approach in the analysis of an irregularly sampled speleothem proxy record with seasonal laminae from Niue island in the central tropical Pacific. Application of the proposed correction scheme identifies a spurious transition that is solely imposed by an abrupt shift in sampling rate and uncovers periods of reduced seasonal rainfall predictability associated with enhanced ENSO and tropical cyclone activity.
M. H. Trauth, A. Asrat, A. S. Cohen, W. Duesing, V. Foerster, S. Kaboth-Bahr, K. H. Kraemer, H. F. Lamb, N. Marwan, M. A. Maslin, F. Schäbitz:
Tipping points in the 620 kyr proxy record from Chew Bahir, S Ethiopia,
6th PAGES Open Science Meeting,
(online meeting),
May 18, 2022,
Talk.
» Abstract
We have used a change point analysis (CPA) and a recurrence plot/recurrence quantification analysis (RP/RQA) on a 300 m / 620 kyr lake-sediment record from the Chew Bahir basin in the southern Ethiopian Rift to determine the amplitude and duration of past climate transitions. In this record, there are numerous transitions from wet to dry conditions, as well as from dry to wet, which show the typical characteristics of a tipping point, where the change is always faster than the forcing and the actual transition is preceded by possible precursor events. One of the most interesting transition examined with the CPA and RP/RQA was the termination of the African Humid period (15–5 kyr BP). The rapid ( 880 yr) change of climate in response to a relatively modest change in orbital forcing appears to be typical of tipping points in complex systems such as the Chew Bahir basin. If this is the case then 14 dry events at the end of the AHP at 5.5 kyr BP, each of them 20–80 yrs long and recurring every 160±40 yrs as documented in the Chew Bahir cores could represent precursors of an imminent tipping point which, if properly interpreted, would allow predictions to be made of future climate change in the Chew Bahir basin. Compared to the low-frequency cyclicity of climate variability before and after the termination of the AHP, this type of cyclicity occurs on time scales equivalent to a few human generations. In other words, it is very likely (albeit speculative) that people were conscious of these changes and adapted their lifestyles to the consequent changes in water and food availability. A deeper analysis of our data is however required to understand whether the wet-dry climate transition in the area was due to a saddle-node bifurcation in the structural stability of the climate, or whether it was induced by a stochastic fluctuation.
N. Marwan:
Palaeoclimate variability in Central America in the last 2 millenia,
Thematic Einstein Semester "The Mathematics of Complex Social Systems: Past, Present, and Future",
Berlin (Germany),
April 25, 2022,
Talk.
» Abstract
Description: (i) Stable isotope (d18O) climate record derived from stalagmite YOK-I from Yok Balum Cave in Belize, representing regional palaeoclimate variation between 40 BC and 2006 AD. (ii) Reconstruction of tropical Atlantic sea surface temperatures (SSTs) spanning the last 2000 years using seasonally representative foraminifera from the Cariaco Basin.
Background: The stalagmite based stable isotope climate record from Belize represents variability in tropical rainfall over the last 2000 years. This variability is related to a displacement of the Intertropical Convergence Zone (ITCZ) which seems to be controlled by the North Atlantic Oscillation (NAO) or changes in the tropical Atlantic sea surface temperatures (SSTs)
Original Purpose: The stable isotope record of stalagmite YOK-I is a reference record of regional past rainfall variability, used to investigate, e.g., droughts and their impact on politics, war, and population fluctuations of the Mayans. Related studies: 10.1126/science.1226299, 10.1038/srep45809, 10.1002/2013GL058458.
Questions: Is there a relationship between the SST variability of the tropical Atlantic and the rainfall variability in Belize (including leads and lags)? Does this relationship change over time? Consider uncertainties in the dating procedure and include them in the analysis.
N. Marwan:
Nonlinear Data Analysis Concepts,
Graduate School NatRiskChange, University of Potsdam,
Potsdam (Germany),
March 28-29, 2022,
Lecture.
» Abstract
The lecture introduces the basic concepts of nonlinear dynamics and chaos and how they can be applied for the study of complex systems, spatiotemporal data, and nonlinear interrelationships in geosciences. The specific topics contain
- Basic terminology, dynamical systems, and simple prototypical models
- Dimensions, fractals
- Concept of symbolic dynamics
- Concept of phase space, phase space reconstruction, Lyapunov exponent and correlation sum
- Concept of recurrence in phase space, recurrence plots, recurrence quantification analysis
- Detection of regime transitions, statistical tests
- Concept of synchronization, coupling analysis
- Spatial and spatio-temporal data analysis using recurrence features
- Complex networks, network models, measures, network representations
- Functional networks, reconstruction of networks, climate networks
- Complex networks based time series analysis
R. Krishnan, M. Singh, T. P. Sabin, B. Goswami, A. D. Choudhury, P. Swapna, R. Vellore, A. G. Prajeesh, N. Sandeep, C. Venkataraman, R. V. Donner, N. Marwan, J. Kurths:
Implications of volcanic aerosols for seasonal forecasting of the Indian monsoon in a changing climate,
Seventh WMO International Workshop on Monsoons (IWM-7),
New Delhi (India),
March 22-26, 2022,
Talk invited.
» Abstract
There is unequivocal evidence that human-induced climate change, in particular greenhouse gas (GHG) emissions, has been the main driver of the observed intensification of heavy precipitation over the land regions across the globe, and has also contributed to increases in agricultural droughts in some regions, which are further projected to enhance with additional warming during the 21st century (IPCC AR6 WG1, 2021). In addition to GHG forcing, anthropogenic aerosol emissions from the Northern Hemisphere (NH) are recognized to have influenced monsoon precipitation changes over the West African, South Asian and East Asian monsoon regions, since the second half of the 20th century (IPCC AR6 WG1, 2021). In particular, the expected enhancement of the South Asian monsoon precipitation by GHG forcing since 1950s has been offset by precipitation reduction caused by the NH anthropogenic aerosols (IPCC AR6 WG1, 2021).
Near-term climate projections for the period 2021-2040 indicate that the South Asian monsoon will be dominated by the effects of internal variability, but will increase in the long-term (IPCC AR6 WG1, 2021). In this context, it must be highlighted that uncertainties due to unpredictable natural forcings such as large volcanic eruptions can lower the degree of confidence in projecting near-term monsoonal changes. This talk is aimed to provide some insights into the role of large volcanic eruptions on the tropical atmosphere-ocean coupled system and the Indian monsoon, with implications for monsoon seasonal forecasting.
N. Marwan:
Transparent and efficient data storage,
Graduate School NatRiskChange University of Potsdam,
Potsdam (Germany),
March 14-15, 2022,
Lecture.
» Abstract
The lecture provides an overview of the need for sustainable storage of scientific data, various concepts of data storage and archiving, their planning and practical implementation, both at the personal, institutional, and public levels in publicly accessible data archives. Specific topics discussed include reproducibility and transparency, important data formats, data integrity, standards, encryption, backup, coding conventions, documentation/meta-data, and version control.
M. Gadhawe, R. Guntu, A. Banerjee, N. Marwan, A. Agarwal:
A complex network approach to study the extreme precipitation patterns in a river basin,
AGU 2021 Fall Meeting,
New Orleans (USA),
December 13–17, 2021,
DOI:10.1002/essoar.10509273.1,
Poster.
» Abstract
The spatiotemporal patterns of precipitation are critical for understanding the underlying mechanism of many hydrological and climate phenomena. Over the last decade, applications of the complex network theory as a data-driven technique has contributed significantly to study the intricate relationship between many variable in a compact way. In our work, we conduct a study to compare an extreme precipitation pattern in Ganga River Basin, by constructing the networks using two nonlinear methods - event synchronization (ES) and edit distance (ED). Event synchronization has been frequently used to measure the synchronicity between the climate extremes like extreme precipitation by calculating the number of synchronized events between two events like time series. Edit distance measures the similarity/dissimilarity between the events by reducing the number of operations required to convert one segment to another, that consider the events’ occurrence and amplitude. Here, we compare the extreme precipitation patterns obtained from both network construction methods based on different network’s characteristics. We used degree to understand network topology and identify important nodes in the networks. We also attempted to quantify the impact of precipitation seasonality and topography on extreme events. The study outcomes suggested that the degree is decreased in the southwest to the northwest direction and the timing of peak precipitation influences it. We also found an inverse relationship between elevation and timing of peak precipitation exists and the lower elevation greatly influences the connectivity of the stations. The study highlights that Edit distance better captures the network’s topology without getting affected by artificial boundaries.
T. Braun, S. Breitenbach, E. Ray, J. U. L. Baldini, L. M. Baldini, F. Lechleitner, Y. Asmerom, K. M. Prufer, N. Marwan:
Two millennia of seasonal rainfall predictability in the neotropics with repercussions for agricultural societies,
AGU 2021 Fall Meeting,
New Orleans (USA),
December 13–21, 2021,
Talk.
» Abstract
The reconstruction and analysis of palaeoseasonality from speleothem records remains a notoriously challenging task. Although the seasonal cycle is obscured by noise, dating uncertainties and irregular sampling, its extraction can identify regime transitions and enhance the understanding of long-term climate variability. Shifts in seasonal predictability of hydroclimatic conditions have immediate and serious repercussions for agricultural societies. We present a highly resolved speleothem record (ca. 0.22 years temporal resolution with episodes twice as high) of palaeoseasonality from Yok Balum cave in Belize covering the Common Era (400-2006 CE) and demonstrate how seasonal-scale hydrological variability can be extracted from δ13C and δ18O isotope records. We employ a Monte-Carlo based framework in which dating uncertainties are transferred into magnitude uncertainty and propagated. Regional historical proxy data enable us to relate climate variability to agricultural disasters throughout the Little Ice Age and population size variability during the Terminal Classic Maya collapse.
Spectral analysis reveals the seasonal cycle as well as nonstationary ENSO- and multi-decadal-scale variability. The degree to which farmers can reliably predict crop yield is both affected by long-term hydroclimate conditions and short-term variations in the subannual distribution of rainfall. A recurrence analysis reveals transitions in seasonal-scale predictability and links them to the mean hydroclimate. Predictability of seasonal rainfall variations was rendered progressively less predictable during the classic Maya collapse. These results are discussed in the context of their implications for rainfall dependent agricultural societies.
N. Marwan:
Nonlinear Time Series Analysis in Geosciences,
Kolloquium of Institute of Geoscience, University of Potsdam,
Potsdam (Germany),
December 13, 2021,
Lecture.
N. Marwan, H. Kraemer:
Hands-on workshop on complex systems,
GFZ Potsdam,
Potsdam (Germany),
November 4–5, 2021,
Lecture and workshop.
M. Kemter, N. Marwan, G. Villarini, B. Merz:
United States Flood Trends and Their Drivers,
Second International Conference on Natural Hazards and Risks in a Changing World,
Potsdam (Germany),
October 5-6, 2021,
Poster.
» Abstract
Climate change has already altered the magnitude and frequency of river floods around the world and if these trends persist in the future, our flood risk management will have to adapt to the changing conditions. However, to predict future changes, we must first understand how floods have changed in the past. Here we present a study of more than 4000 river gauges across the USA with time series of annual maximum streamflow between 1960–2010. We use a novel clustering approach to find 12 hydro-meteorologically and spatially distinct clusters of catchments with similar flood behavior. Based on re-analysis data we calculate more than 30 hydro-climatological and land-use variables to use them as predictors for 12 separate Random Forest models, one for each cluster. With these we search for the drivers of past trends in flood magnitudes, differences between common and rare floods as well as in the synchrony of floods in different catchments. We use Accumulated Local Effect plots to understand how each of the predictors affected the different trends. We find that in many regions changes in precipitation (e.g. annual precipitation sum, flood generating precipitation) translated to similar changes in flood magnitudes. Furthermore, we show that static land use conditions were of unexpected importance for flood trends, especially in the form of canopy cover, reservoirs, and impervious surfaces. We show that forests and reservoirs have mitigated some of the effects of a changing climate on floods, while urbanization has amplified them. We find varying importance and occasionally opposing effects of the same predictors in different clusters, showing that flood trends are highly dependent on the hydro-climatological circumstances of the catchments. Our results highlight the importance of a holistic approach to the analysis of flood trends and their drivers, considering a wide and consistent range of variables while taking into account the regional variability in flood generation.
A. Banerjee, B. Goswami, N. Marwan, B. Merz, B., J. Kurths:
Complex network approach to study the impacts of ENSO over the United States,
Second International Conference on Natural Hazards and Risks in a Changing World,
Potsdam (Germany),
October 5-6, 2021,
Poster.
» Abstract
Complex network analysis is a powerful tool that encodes the intricate relationship between the many components of a complex system. Functional climate network analysis is particularly designed to study the complex interactions between the different components of the Earth’s climate system. As the weather changes, the dynamical interaction between the grid points or the nodes should also change. We have used an evolving climate network approach to study this evolution of interaction between the nodes over time. In our work, we study temperature and precipitation patterns in the United States using this framework. El Niño–Southern Oscillation (ENSO) is one of the most important sources of annual global climate variability, associated with characteristic patterns of rainfall and temperature, includinng extreme events such as floods and droughts. The United States is one of the most susceptible regions to severe weather outbreaks due to the ENSO. We use Pearson correlation and edit-distance methods as similarity measures to construct temperature and extreme precipitation networks respectively.
We study the course of evolution of the link patterns using network measures such as robust links, link density, and transitivity during the different phases of the ENSO. Through our analysis, we are able to distinguish between the different phases of the ENSO, and hence, identify the different large-scale atmospheric circulation patterns associated with them.
S. M. Vallejo-Bernal, F. Wolf, N. Boers, N. Marwan, J. Kurths:
Synchronicity of heavy rainfall induced by atmospheric rivers over North America,
Second International Conference on Natural Hazards and Risks in a Changing World,
Potsdam (Germany),
October 5-6, 2021,
Poster.
» Abstract
Atmospheric rivers (ARs) are dynamical features of the low atmosphere, responsible for much of the moisture transport in the midlatitudes, that can produce copious amounts of precipitation as long as an external uplifting mechanism is available. In particular, land-falling ARs are strongly linked to heavy precipitation over the orographically complex western coast of the United States. Recently, complex network approaches have proven to effectively extract spatiotemporal variability patterns from climate data and have contributed to significant advances in the understanding and prediction of extreme weather events. In this study, we investigate the synchrony and interdependency of heavy rainfall occurrences related to ARs along the west coast of North America and the underlying physical mechanisms. We use the SIO R1 catalog of land-falling ARs and the daily rainfall estimates of the ERA5 reanalysis project during the period 1980–2018 to construct complex climate networks by applying the nonlinear Event Synchronization measure to the rainfall events above the 95th percentile. Our results reveal substantially different spatiotemporal rainfall patterns depending on the presence (or absence) of ARs. On the one hand, we find that ARs are linked with highly interconnected regions of synchronous heavy rainfall along the coastline and the adjacent Pacific Ocean. On the other hand, weaker but significant connections are observed over the continental North America in the absence of land-falling ARs. Also, the underlying atmospheric conditions differ visibly and exhibit a robust decadal pattern that is, however, highly variable for seasonal means. Resolving the typical synchrony structures of heavy rainfall related to the land-falling ARs, should lead to improved understanding of hydroclimate variability, likely leading to improved seasonal predictability of extreme precipitation.
N. Marwan:
Recent exciting developments in recurrence plot analysis,
9th International Symposium on Recurrence Plots,
Lublin (Poland),
September 22–24, 2021,
Talk.
» Abstract
I reported on recent achievements in recurrence plot research ten years ago. It is time to look back at the last ten years of exciting developments that had been achieved for improving recurrence plot analysis and to widen its application potential. I will give a brief overview about important and innovative developments, such as conceptual recurrence plots, ideas for parameter selection, recurrence grammars, event-like, multiscale, and heterogeneous recurrences, and correction schemes. New perspectives have recently been opened by combining recurrence plots with machine learning.
A. Das, A. P. Nandan, N. Marwan, A. Koseska:
Identifying transient metastable states from live-cell imaging data,
9th International Symposium on Recurrence Plots,
Lublin (Poland),
September 22–24, 2021,
Talk.
» Abstract
The current models of cellular information processing are formulated with the assumption that the relevant dynamics occur near or at the steady-state. One of the advantages when focusing on the asymptotic behavior at or near a steady-state is that it simplifies the analysis of the system. Recently, however, we proposed theoretical and practical reasons that steady-state analysis misses essential behavior of living systems, especially when describing how single cells process non-stationary signals. Using fluorescence measurements of real-time protein activities in single cells, we demonstrate that protein networks utilize transient dynamics away from steady-state to maintain a memory of previously encountered signals and thereby navigate in complex environments. To identify such out-of-equilibrium transients from experimental data in general, we develop a recurrence-based method that can be applied to a broad class of systems, without prior knowledge of the underlying dynamics of the system.
I. Pavithran, V. R. Unni, R. I. Sujith, J. Kurths, N. Marwan:
Recurrence condensation during critical transitions in complex systems,
9th International Symposium on Recurrence Plots,
Lublin (Poland),
September 22–24, 2021,
Talk.
» Abstract
Many dynamical systems exhibit critical transitions, which can result in catastrophic changes to the state of the system. Spontaneous emergence of an ordered dynamics from disorder is common among several systems. Even inherently turbulent systems can undergo such transition to ordered periodic dynamics. A recent discovery shows that the emergence of self-sustained periodic oscillations from an initially disordered state in various systems is accompanied by the phenomenon of spectral condensation. It is the sharpening of the peak in the amplitude spectrum accompanied by the growth of amplitude of the dominant oscillatory mode.
In the present work, we identify an analogous phenomenon known as recurrence condensation in recurrence plots (RPs) of the system variables during the transition to self-sustained periodic oscillations. We construct RPs from the time series of state variables by fixing the recurrence rate to be constant. During the transition from a chaotic state to periodic oscillations via intermittency, the RP changes from a disordered arrangement of short, broken diagonal lines to patches of ordered short diagonal lines and then to a pattern with long diagonal lines. Here, the scattered RP changes to diagonal lines, and their spacing corresponds to the dominant time scale in the system.
We aim to quantify recurrence condensation using recurrence quantification measures, an appropriate method for analyzing nonlinear, high-dimensional complex systems. Here, we use these measures to study the emergence of periodicity, which is reflected as long diagonal lines in RPs. The recurrence measures such as determinism, entropy, laminarity, and trapping time exhibit a gradual variation during recurrence condensation. Further, these recurrence measures follow a power-law scaling with the control parameter on approaching the transition. These scaling relations hold only till a critical point during the transition.
The best power-law fit with minimum scatter is observed when considering the measures only up to a particular control parameter. The fitting error increases afterwards, indicating a deviation from the scaling. Utilizing this property, we find a critical value of the control parameter up to which the recurrence measures follow scalings. Further, we construct RPs of synthetic data from a noisy Hopf bifurcation model and uncover that the value of the critical point detected is the same as the Hopf point. Our analysis indicates that the finding of the optimal power law can be used to define critical points in noisy systems with gradual transitions, where the transition point is not easily identifiable.
T. Braun, V. R. Unni, R. I. Sujith, J. Kurths, N. Marwan:
Multiscale recurrence quantification with recurrence lacunarity,
9th International Symposium on Recurrence Plots,
Lublin (Poland),
September 22–24, 2021,
Talk.
» Abstract
We propose Recurrence Lacunarity (RL) as a novel recurrence quantification measure. Lacunarity, as originally introduced by B.Mandelbrot (1983), is usually interpreted as a measure of heterogeneity or translational invariance of some spatial pattern. Compared to traditional recurrence quantifiers, RL offers the advantage that it is sensitive to more general geometric features beyond line structures of a recurrence plot. It is demonstrated how RL can be used to characterize the recurrence between dynamical states at all relevant time scales of a complex system. An application to time series of acoustic pressure fluctuations from a turbulent combustor showcases the great potential of RL in detecting abrupt regime shifts of different origin in nonlinear dynamical systems. Finally, potential future applications and conceptual extensions of RL are discussed.
K. H. Kraemer, F. Hellmann, J. Kurths, N. Marwan:
Recurrence powerspectra,
9th International Symposium on Recurrence Plots,
Lublin (Poland),
September 22–24, 2021,
Talk.
» Abstract
A novel kind of powerspectrum is constructed, the textitspike powerspectrum, which transforms spike-train-like signals into their frequency domain. This method clearly shows the apparent cycles in the data and overcomes the problems when using the obvious idea of Fourier-transforming it. We invent this instructive approach with the idea of transforming the $\tau$-recurrence rate of a recurrence plot (RP), which often has a spiky appearance. The $\tau$-recurrence rate is the density of recurrence points along diagonals of the RP, which are parallel to the main diagonal with a distance of $\tau$. In this context the spike powerspectrum can be interpreted as a nonlinear power spectrum of a potentially high dimensional system which constitutes the RP. The proposed measure is simple to compute, robust to noise and is able to detect bifurcations inducing regular-regular, regular-chaos as well as chaos-chaos transitions.
C. Ozdes, S. Breitenbach, D. Eroglu, N. Marwan, N:
Revealing past climate networks from data,
9th International Symposium on Recurrence Plots,
Lublin (Poland),
September 22–24, 2021,
Talk.
» Abstract
Real-world complex systems such as the climate system of interacting spatially-diverse weather patterns or ecological communities are essential parts of our everyday lives. Such systems are composed of a large number of units represented by the nodes of a network enacting an intricate interaction structure. These structures are the topic of network science, a very active research field significantly improving the quality of living standards by forecasting the dynamical behavior of complex systems. These interaction structures are also very hard to detect in climate networks, where a reconstruction procedure requires dealing with sparse and noisy data coming from multi-scale spatiotemporal dynamics. Paleoclimate data sets are also not regularly sampled, and the measurements are not precise. To model this error process, we fit a trained transformation-cost time series to paleoclimate data using the metric of segment similarity, which maximizes prior normality. This transformation-cost series can act as a regularly sampled proxy to be used in detecting the dynamical properties of data sets using recurrence quantification methods. This procedure reveals interactions that are often not detectable in the original data. We then use sliding windows analyses of these recurrence measures to obtain the temporal correlation networks that define interaction between our data sources. We use this method to reconstruct the evolution of the climate interaction network for Asia to reveal its transient and stable interaction patterns in the past.
B. G. Straiotto, N. Marwan, D. C. James, P. J. Seeley:
A combination of principal component and recurrence analyses discriminates between closely similar movement patterns,
9th International Symposium on Recurrence Plots,
Lublin (Poland),
September 22–24, 2021,
Talk.
» Abstract
Investigation of the movement of body segments and their coordination can provide insight into underpinning motor control strategies. We aimed to determine whether combined application of principal component and recurrence quantification analyses might discriminate between spatial and temporal aspects of apparently similar movement patterns. Backwards-forwards movements of elite (n = 9) and non-elite (n = 9) taekwondo players employed in both defensive and attacking actions were recorded using motion capture techniques, and features of whole-body movement defined at segment level were investigated by principal component analysis. For both groups of players, four movement components explained > 90% of the variability in the data.
The time series derived from scores for each of the principal components were subsequently subjected to recurrence quantification analysis, player by player. For the first component, statistically significant differences between groups were detected for the recurrence measurement determinism (p < 0.05). For the third component, statistically significant differences were detected for the recurrence measurements laminarity and maxline (p < 0.01). Application of surrogation to recurrence data indicated that differences between elite and non-elite groups were deterministic in origin and not the result of data noise ( p < 0.01).
Contrasting operation of taekwondo movements within and between our player groups were revealed qualitatively in the coefficients derived from principal component analysis and quantitatively through combined principal component and recurrence techniques. Variations in TKD player execution of this apparently simple movement pattern may represent more skilled motor control in elite players that is related to the functional importance of backwards-forwards movements in sparring and competition.
N. Marwan:
Co-authorship network of the recurrence plot domain,
9th International Symposium on Recurrence Plots,
Lublin (Poland),
September 22–24, 2021,
» Poster (PDF, 58.87M)
.
N. Marwan, T. Braun:
Module 5 – Complex Systems, Recurrence and Networks in Climate,
Summer school on Trends, rhythms and events in the Earth's climate system,
Online/ Potsdam (Germany),
August 30 – September 3, 2021,
Lecture.
» Abstract
Many processes in the Earth's climate system are complex, behave in an unpredictable way, and are nonlinearly coupled, which is why linear methods of time series analysis fail to characterize them or to uncover their complex interrelationships. This module will teach fundamental properties of complex systems and the basics of dynamical systems theory. Examples from atmospheric science and paleoclimate variability will be used to demonstrate the challenges in investigating such systems. Modern methods from nonlinear dynamics will be introduced, such as recurrence-based methods and complex networks. The students will learn in which applications these methods will provide complimentary information, how to reliably apply these methods, and how to develop corresponding statistical tests. Moreover, concepts to address common challenges in (paleo-)climate data analysis (such as heavy-tailed and event-like data, dating uncertainties, and irregular sampling) will be introduced. A focus will be on detecting subtle regime changes and complex interrelationships, where linear methods fail.
M. Kemter, N. Marwan, G. Villarini, B. Merz:
Drivers of Flood Trends in the United States,
Tag der Hydrologie,
Potsdam (Germany),
August 30 – September 1, 2021,
Poster.
I. Pavithran, V. R. Unni, A. Saha, A. J. Varghese, R. I. Sujith, J. Kurths, N. Marwan:
Predicting the Amplitude of Thermoacoustic Instability Using Universal Scaling Behaviour,
ASME Turbo Expo 2022,
Rotterdam (The Netherlands),
June 9, 2021,
Talk.
» Abstract
The complex interaction between the turbulent flow, combustion and the acoustic field in gas turbine engines often results in thermoacoustic instability that produces ruinously high-amplitude pressure oscillations. These self-sustained periodic oscillations may result in a sudden failure of engine components and associated electronics, and increased thermal and vibrational loads. Estimating the amplitude of the limit cycle oscillations that are expected during thermoacoustic instability helps in devising strategies to mitigate and to limit the possible damages due to thermoacoustic instability. We propose two methodologies to estimate the amplitude using only the pressure measurements acquired during stable operation. First, we use the universal scaling relation of the amplitude of the dominant mode of oscillations with the Hurst exponent to predict the amplitude of the limit cycle oscillations. We also present a methodology to estimate the amplitudes of different modes of oscillations separately using “spectral measures,” which quantify the sharpening of peaks in the amplitude spectrum. The scaling relation enables us to predict the peak amplitude at thermoacoustic instability, given the data during the safe operating condition. The accuracy of prediction is tested for both methods, using the data acquired from a laboratory-scale turbulent combustor. The estimates are in good agreement with the actual amplitudes.
A. Banerjee, B. Goswami, N. Marwan, B. Merz, J. Kurths:
Recurrence based coupling analysis between event-like data and continuous data,
(Virtual) EGU General Assembly,
Vienna (Austria),
April 19–30, 2021,
DOI:10.5194/egusphere-egu21-14831,
Talk.
» Abstract
Extreme events such as earthquakes, tsunamis, heat weaves, droughts, floods, heavy precipitation, or tornados – affect the human communities and cause tremendous loss of property and wealth, but can be related to multiple and complex sources. For example, a flood is a natural event caused by many drivers such as extreme precipitation, soil moisture, or temperature. We are interested in understanding the direct and indirect coupling between flood events with different climatological and hydrological drivers such as soil moisture and temperature.
We use multivariate recurrence plot and recurrence quantification analysis as a powerful framework to study the couplings between the different systems, especially the direction of coupling. The standard delay-embedding method is not a suitable for the recurrence analysis of event-like data. Therefore, we apply the novel edit-distance method to compute recurrence plots of time series of flood events and use the standard recurrence plot method for the continuous varying time series such as soil moisture and temperature. The coupling analysis is performed using the mean conditional probabilities of recurrence derived from the different recurrence plots. We demonstrate this approach on a prototype system and apply it on the hydrological data. Using this approach we are able to indicate the coupling direction and lag between the different coupled systems.
T. Braun, S. Breitenbach, E. Ray, J. U. L. Baldini, L. M. Baldini, F. Lechleitner, Y. Asmerom, K. M. Prufer, N. Marwan:
Two millennia of seasonal rainfall predictability in the neotropics with repercussions for agricultural societies,
(Virtual) EGU General Assembly,
Vienna (Austria),
April 19–30, 2021,
DOI:10.5194/egusphere-egu21-11012,
Talk.
» Abstract
The reconstruction and analysis of palaeoseasonality from speleothem records remains a notoriously challenging task. Although the seasonal cycle is obscured by noise, dating uncertainties and irregular sampling, its extraction can identify regime transitions and enhance the understanding of long-term climate variability. Shifts in seasonal predictability of hydroclimatic conditions have immediate and serious repercussions for agricultural societies.
We present a highly resolved speleothem record (ca. 0.22 years temporal resolution with episodes twice as high) of palaeoseasonality from Yok Balum cave in Belize covering the Common Era (400-2006 CE) and demonstrate how seasonal-scale hydrological variability can be extracted from δ13C and δ18O isotope records. We employ a Monte-Carlo based framework in which dating uncertainties are transferred into magnitude uncertainty and propagated. Regional historical proxy data enable us to relate climate variability to agricultural disasters throughout the Little Ice Age and population size variability during the Terminal Classic Maya collapse.
Spectral analysis reveals the seasonal cycle as well as nonstationary ENSO- and multi-decadal-scale variability. Variations in both the subannual distribution of rainfall and mean average hydroclimate pose limitations on how reliably farmers can predict crop yield. A characterization of year-to-year predictability as well as the complexity of seasonal patterns unconver shifts in the seasonal-scale variability. These are discussed in the context of their implications for rainfall dependent agricultural societies.
Z. Su, S. Gupta, N. Marwan, N. Boers, J. Kurths:
A comparative study of extreme precipitation patterns using complex networks,
(Virtual) EGU General Assembly,
Vienna (Austria),
April 19–30, 2021,
DOI:10.5194/egusphere-egu21-8740,
Talk.
» Abstract
The spatio-temporal patterns of precipitation are of considerable relevance in the context of understanding the underlying mechanism of climate phenomena. The application of the complex network paradigm as a data-driven technique for the investigation of the climate system has contributed significantly to identifying the key regions influencing the climate variability of a target region of interest and, in particular, to improving the predictability of extreme events. In our work, we conduct a comparative study of precipitation patterns by constructing functional climate networks using two nonlinear event similarity measures – event synchronization (ES) and edit-distance (ED). Event synchronization has been widely applied to identify interactions between occurrences of different climate phenomena by counting the number of synchronized events between two event series. Edit-distance measures the similarity between sequences by minimizing the number of operations required to transform one sequence to another. We suggest edit-distance as an alternative approach for network reconstruction that can measure similarity between two event series by incorporating not only event occurrences but also event amplitudes. Here, we compare the global extreme precipitation patterns obtained from both reconstruction methods based on the topological characteristics of the resulting networks. As a case study, we compare selected features of network representations of East Asian heavy precipitation events obtained using both ES and ED. Our results reveal the complex nature of the interaction between the Indian Summer Monsoon (ISM) and the East Asian Summer Monsoon (EASM) systems. Through a systematic comparison, we explore the limitations of both measures and show the robustness of the network structures.
H. Kraemer, G. Datseris, J. Kurths, I. Kiss, J. L. Ocampo-Espindola, N. Marwan:
A unified and automated approach to attractor reconstruction,
(Virtual) EGU General Assembly,
Vienna (Austria),
April 19–30, 2021,
DOI:10.5194/egusphere-egu21-1495,
Talk.
» Abstract
Since acquisition costs for sensors and data collection decrease rapidly especially in the geo-scientific fields, researchers often have to deal with a large amount of multivariable data, which they would need to automatically analyze in an appropriate way. In nonlinear time series analysis, phase space reconstruction often makes the very first step of any sophisticated analysis, but the established methods are either unable to reliably automate the process or they can not handle multivariate time series input. Here we present a fully automated method for the optimal state space reconstruction from univariate and multivariate time series. The proposed methodology generalizes the time delay embedding procedure by unifying two promising ideas in a symbiotic fashion. Using non-uniform delays allows the successful reconstruction of systems inheriting different time scales. In contrast to the established methods, the minimization of an appropriate cost function determines the embedding dimension without using a threshold parameter. Moreover, the method is capable of detecting stochastic time series and, thus, can handle noise contaminated input without adjusting parameters. The superiority of the proposed method is shown on some paradigmatic models and experimental data.
J. Fohlmeister, N. Sekhon, A. Columbu, K. Rehfeld, L. Sime, C. Veige-Pires, N. Marwan, N. Boers:
Global reorganization of atmospheric circulation during Dansgaard-Oschger cycles,
(Virtual) EGU General Assembly,
Vienna (Austria),
April 19–30, 2021,
DOI:10.5194/egusphere-egu21-9433,
Talk.
» Abstract
Ice core records from Greenland provide evidence for multiple abrupt warming events recurring at millennial time scales during the last glacial interval. Although climate transitions strongly resembling these Dansgaard-Oeschger (DO) transitions have been identified in several speleothem records, our understanding of the climate and ecosystem impacts of the Greenland warming events in lower latitudes remains incomplete.
Here, we investigate the influence of DO transitions on the global atmospheric circulation pattern. We comprehensively analyse d18O changes during DO transitions in a globally distributed dataset of speleothems (SISALv2; Comas-Bru et al., 2020). Speleothem d18O signals mostly reflect changes in precipitation amount and moisture source. Thereby this proxy allows us to infer spatially resolved changes in global atmospheric dynamics that are characteristically linked to DO transitions. We confirm the previously proposed shift of the Intertropical Convergence Zone towards more northerly positions. In addition, we find evidence for a similar northward shift of the westerly winds of the Northern Hemisphere. Furthermore, we identify a decreasing trend in the transition amplitudes with increasing distances from the North Atlantic region. This confirms previous suggestions of this region being the core and origin of these past abrupt climate changes.
M. Singh, R. Krishnan, B. Goswami, A. Dey Choudhury, S. Panickal, R. Vellore, P. A. Gopinathan, S. Narayanasetti, C. Venkataraman, R. Donner, N. Marwan, J. Kurths:
Fingerprint of volcanic forcing on the ENSO–Indian monsoon coupling,
(Virtual) EGU General Assembly,
Vienna (Austria),
April 19–30, 2021,
DOI:10.5194/egusphere-egu21-9059,
Talk.
» Abstract
The coupling between the El Niño–Southern Oscillation (ENSO) and Indian Monsoon (IM) plays a significant role in the summer rainfall over the Indian subcontinent. In this study, we provide insights into the IM variability with regard to the degree of ENSO variability and radiative forcing from large volcanic eruptions (LVEs). Volcanic dust and gas injected into the stratosphere during major eruptions influence the ENSO from seasonal to interannual timescales. However, the effects of LVEs on the ENSO-IM coupling remain unclear. The relationship between ENSO and IM systems in the context of LVEs is examined using a panoply of datasets and advanced statistical analysis techniques in this study. We find that there is a significant enhancement of the phase-synchronization between ENSO and IM oscillations due to increase in angular frequency of ENSO in the last millennium. Twin surrogates-based statistical significance testing is also used to affirm this result and similar evidence is found in the combinations of 14 ENSO and 11 IM paleoclimate proxy records in the last millennium. Bayesian probabilities conditioned with and without LVEs show LVEs lead to a strong ENSO-IM phase-coupling, with the probabilities remaining higher till the fourth year from the eruption. A large-ensemble climate model experiment with and without the 1883 Krakatoa eruption is conducted using the IITM-ESM, and also with varied volcanic radiative forcing (VRF) depending on the evolved state of ENSO. The simulations show that LVEs force the ENSO-IM systems into a coupled state, and increase (decrease) in the VRF leads to an enhanced (decreased) probability of the phase synchronisation of ENSO-IM systems with a high chance of El Niño-IM drought in the year following the LVE. Our results promisingly pave a way not only for improving the seasonal monsoon prediction improvements but also for the regional impact assessment from the proposed geo-engineering activities over the South Asian region.
N. Marwan:
Reconstructing Complex Networks from Data,
Seminar Networks Unit, IMT School For Advanced Studies,
Lucca (Italy),
February 4, 2021,
Talk invited.
» Abstract
Complex networks provide an interesting tool to investigate spatio-temporal data. The first step is to reconstruct a (functional) network from data. I will show different reconstruction approaches depending on the research question and the nature of the data. The procedure is illustrated with applications on climate data.
M. Kemter, B. Merz, N. Marwan:
Trends in flood magnitudes, flood extents and their drivers,
AGU Fall Meeting,
online,
December 1–17, 2020,
Talk.
» Abstract
When rivers flood, they often do so simultaneously with surrounding rivers. We call the area across which this happens the flood extent. We study trends in flood extents in Europe and the US, and analyze how they are linked to trends in flood magnitude. For the time period 1960-2010 we analyze the annual maximum floods of 5000 US- and 4000 European hydrometric stations and investigate generating processes of all floods based on climate re-analysis data. We use event synchronization and complex networks to find groups of stations that have similar flood behavior. For both regions, we find a positive correlation between extents and magnitudes for 93% of the stations. While the trends of both variables are well aligned in Europe (increasing in the west, decreasing in the south), the same is not true for the US. The station groups help us to understand these differences. We show that the changing influence of snowmelt on flood generation is crucial for the relationship of magnitude and extent trends.
A. Agarwal, N. Marwan, R. Maheswaran, U. Ozturk, J. Kurths, B. Merz:
Optimal design of hydrometric station networks based on complex network analysis,
AGU Fall Meeting,
online,
December 1–17, 2020,
Talk.
» Abstract
Hydrometric networks play a vital role in providing information for decision-making in water resources management. They should be set up optimally to provide as much and as accurate information 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, yet there is 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 in the last 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), to evaluate the importance of nodes in a network. It is compared to previously proposed measures 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 network efficiency and the kriging error. The results suggest that WDB effectively quantifies the importance of rain gauges, although the benefit of the method needs to be investigated in more detail.
T. Westerhold, N. Marwan, A. Joy Drury, D. Liebrand, C. Agnini, E. Anagnostou, J. Barnet, S. M. Bohaty, D. De Vleeschouwer, F. Florindo, T. Frederichs, D. A. Hodell, A. Holbourn, D. Kroon, V. Lauretano, K. Littler, L. J. Lourens, M. W. Lyle, H. Paelike, U. Roehl, J. Tian, R. Wilkens, Paul. A. Wilson, J. C. Zachos:
Changing state of Earth’s climate for the last 66 million years,
AGU Fall Meeting,
online,
December 1–17, 2020,
Talk.
» Abstract
We combined the best available high-resolution ocean drilling records with newly generated data to produce a continuous, astronomically tuned 66-million-year record of global climate that can serve as a new Cenozoic global reference benthic foraminiferal carbon and oxygen isotope dataset (CENOGRID). The CENOGRID represents the first community effort to systematically assemble a high-fidelity deep-sea isotope record that captures the high and low frequency variations of the climate system on a global scale, providing the framework required for understanding the nature of the major climate states, transitions and events. Within the constraints imposed by temporal resolution in older segments and potential artefacts related to differences in the regional response to orbital forcing, (e.g., obliquity), this singly unique record reveals how long-term gradual shifts in boundary conditions, principally paleogeography, ice-volume, and the mean level of greenhouse gases create state dependent shifts in the sensitivity of the climate system to periodic oscillations in solar forcing. We believe that this phenomenon is related to feedbacks associated with the carbon cycle supporting the critical role of GHG in enhancing climate sensitivity to solar radiative forcing.
N. Marwan:
Measuring complexity of recurrence plots,
GMT Morning Workshop on Nonlinear Dynamics and Statistics,
virtual,
December 3, 2020,
Talk.
N. Marwan:
Investigation of Recurrence Phenomena in the Earth System,
Seminar "Big Data Platform" at DLR Institute of Data Science,
Jena (Germany),
August 6, 2020,
Lecture.
» Abstract
Recurrence is a ubiquitous and fundamental feature in many real world processes. Here, I will focus on recurrence in the system Earth, where it is present at many scales in time and space, such as the rock cycle, activity of an active geyser, celestial mechanics, repeating patterns in a landscape, cycles of glaciation, epochs of geomagnetic polarity, or alternating sediment layers. The study of recurrence properties (such as frequency analysis) can provide deeper insights into the dynamical processes in general. A rather novel approach for the study of recurrences is the so-called recurrence plot and its quantification, rooted in the theory of dynamical systems. In this talk, I will present the basic concept and the major extensions applicable to various research questions. Discussed examples include the temporal change of recurrence properties for identification of regime shifts in climate; spatio-temporal recurrences for classification of landuse dynamics; and bivariate extensions for synchronization/coupling analysis for time scale alignment of palaeoclimate observations. Methodological and numerical Challenges and pitfalls will also be discussed.
N. Marwan:
Investigation of Recurrence Phenomena in the Earth System,
HEIBRIDS Lecture Series,
Berlin (Germany),
July 1, 2020,
Lecture.
» Abstract
Recurrence is a ubiquitous and fundamental feature in many real world processes. Here, I will focus on recurrence in the system Earth, where it is present at many scales in time and space, such as the rock cycle, activity of an active geyser, celestial mechanics, repeating patterns in a landscape, cycles of glaciation, epochs of geomagnetic polarity, or alternating sediment layers. The study of recurrence properties (such as frequency analysis) can provide deeper insights into the dynamical processes in general. A rather novel approach for the study of recurrences is the so-called recurrence plot and its quantification, rooted in the theory of dynamical systems. In this talk, I will present the basic concept and the major extensions applicable to various research questions. Discussed examples include the temporal change of recurrence properties for identification of regime shifts in climate; spatio-temporal recurrences for classification of landuse dynamics; and bivariate extensions for synchronization/coupling analysis for time scale alignment of palaeoclimate observations. Methodological and numerical Challenges and pitfalls will also be discussed.
B. Goswami, A. Hartland, C. Hu, S. Hoepker, B. R. S Fox, N. Marwan, S. F. M. Breitenbach:
Paleo-drip rates from trace metal concentrations in stalagmites: An inverse modeling problem with data uncertainties,
(Virtual) EGU General Assembly,
Vienna (Austria),
May 4-8, 2020,
DOI:10.5194/egusphere-egu2020-19959,
Talk.
» Abstract
The concentration of trace elements such as Ni, Co, and Cu in a stalagmite is determined by (i) the amount of these elements present in so-called organic-metal complexes (OMCs) that trap the ionic forms of such elements in the dripwater, and (ii) the amount that is able to decay from the OMCs into the aqueous phase, from where the elements can adsorb to the growing stalagmite surface (and remain captured within the stalagmite crystal structure). A statistical treatment of the decay of a population of trace element ions from OMCs allow us to model the rates at which the dripwater dropped from the roof of the cave on to the stalagmite's surface. The problem is however made challenging due to: (i) the lack of reliable monitoring data that quantifies the relationship between OMC trace metal ion concentration and stalagmite trace metal ion concentration, and (ii) the presence of chronological uncertainties in our estimates of trace element concentrations at past time points from the depth-based measurements along the stalagmite. We present here a semi-heuristic, semi-theoretical approach that estimates dripwater rates using a theoretical model based on the population-level chemical kinetics of trace element decay from OMCs, and a heuristic choice of calibration data sets based on precipitation and temperature from nearby weather station data. Our approach is applied to trace metal data from the Heshang Cave in southeastern China, and we are able to reconstruct a driprate proxy time series – a first quantitative hydrological proxy record presented along with well-defined estimates of uncertainty.
T. Braun, N. Marwan, V. R. Unni, R. I. Sujith, J. Kurths:
Detection of dynamical regime transitions with lacunarity as a multiscale recurrence quantification measure,
(Virtual) EGU General Assembly,
Vienna (Austria),
May 4-8, 2020,
DOI:10.5194/egusphere-egu2020-3475,
Talk.
» Abstract
We propose Lacunarity as a novel recurrence quantification measure and apply it in the context of dynamical regime transitions. Many complex real-world systems exhibit abrupt regime shifts. We carry out a recurrence plot based analysis for different paradigmatic systems and thermoacoustic combustion time series in order to demonstrate the ability of our method to detect dynamical transitions on variable temporal scales. Lacunarity is usually interpreted as a measure of "gappiness" of an arbitrary spatial pattern. In application to recurrence plots, it quantifies the degree of heterogenity in the temporal recurrent patterns. Our method succeeds to distinguish states of varying dynamical complexity in presence of noise and short time series length. In contrast to traditional recurrence quantifiers, no specification of minimal line lengths is required and features beyond the scope of line structures can be accounted for. Applied to acoustic pressure fluctuation time series, it captures both the rich variability in dynamical complexity and detects shifts of characteristic time scales.
K. Prufer, S. F. M. Breitenbach, J. Baldini, T. Braun, E. Ray, L. Baldini, V. Polyak, F. Lechleitner, N. Marwan, D. Kennett, Y. Asmerom:
A 1,600 year record of paleoseasonality from the neotropics of Central America and its implications for rainfall predictability in agricultural societies,
(Virtual) EGU General Assembly,
Vienna (Austria),
May 4-8, 2020,
DOI:10.5194/egusphere-egu2020-18100,
Talk.
» Abstract
For millions of people living in the humid neotropics seasonally predictable rainfall is crucial for agricultural success and food security. Understanding long-term stability and volatility of seasonal rainfall distributions should be of concern to researchers and policy makers. However, reconstructions of paleorainfall seasonality in the neotropics have been constrained by a lack of precisely dated and sub-annually resolved records. We present a 1,600-year rainfall paleoseasonality reconstruction from speleothem sample Yok G, from Yok Balum Cave located in southern Belize, Central America. Yok G grew continuously from 400 C.E. to 2,006 C.E. and its age is constrained by 52 U-series dates with a mean error of 7 years. The isotope record consists of 7,151 δ18O and δ13C measurements at 0.22-year resolution allowing us to detect the presence and amplitude of annual wet-dry cycles. In Belize rainfall distribution and seasonality controls are currently dominated by the annual migration of the intertropical convergence zone (ITCZ) with marked meridional contrast. The Yok G record suggest distinct changes in seasonality at multi-centennial intervals. The earliest portion of the record (400- 850 C.E.) shows little intra-annual seasonal variation, the period from 850-1400 C.E. has highly variable annual oscillations and periods of low seasonality, while the period from 1,400-2,006 C.E. shows well developed seasonal signals. Element ratios (Mg/Ca, Sr/Ca, and U/Ca) are used to assess Prior Carbonate Precipitation in the epikarst system. We review these changes and the isotopic record from Yok G and discuss tools for interpreting the stability and volatility in seasonal rainfall distributions and possible implications for past and modern agricultural societies.
A. Giesche, S. F. M. Breitenbach, N. Marwan, A. Hartland, B. Plessen, J. F. Adkins, H. H. Haug, A. French, C. A. Petrie, D. A. Hodell:
Rainfall seasonality changes in northern India across the 4.2 ka event,
(Virtual) EGU General Assembly,
Vienna (Austria),
May 4-8, 2020,
DOI:10.5194/egusphere-egu2020-16898,
Talk.
» Abstract
For millions of people living in the humid neotropics seasonally predictable rainfall is crucial for agricultural success and food security. Understanding long-term stability and volatility of seasonal rainfall distributions should be of concern to researchers and policy makers. However, reconstructions of paleorainfall seasonality in the neotropics have been constrained by a lack of precisely dated and sub-annually resolved records. We present a 1,600-year rainfall paleoseasonality reconstruction from speleothem sample Yok G, from Yok Balum Cave located in southern Belize, Central America. Yok G grew continuously from 400 C.E. to 2,006 C.E. and its age is constrained by 52 U-series dates with a mean error of 7 years. The isotope record consists of 7,151 δ18O and δ13C measurements at 0.22-year resolution allowing us to detect the presence and amplitude of annual wet-dry cycles. In Belize rainfall distribution and seasonality controls are currently dominated by the annual migration of the intertropical convergence zone (ITCZ) with marked meridional contrast. The Yok G record suggest distinct changes in seasonality at multi-centennial intervals. The earliest portion of the record (400- 850 C.E.) shows little intra-annual seasonal variation, the period from 850-1400 C.E. has highly variable annual oscillations and periods of low seasonality, while the period from 1,400-2,006 C.E. shows well developed seasonal signals. Element ratios (Mg/Ca, Sr/Ca, and U/Ca) are used to assess Prior Carbonate Precipitation in the epikarst system. We review these changes and the isotopic record from Yok G and discuss tools for interpreting the stability and volatility in seasonal rainfall distributions and possible implications for past and modern agricultural societies.
K. Riechers, N. Boers, J. Fohlmeister, N. Marwan:
Hypothesis testing and uncertainty propagation in paleo climate proxy data evidencing abrupt climate shifts,
(Virtual) EGU General Assembly,
Vienna (Austria),
May 4-8, 2020,
DOI:10.5194/egusphere-egu2020-15148,
Talk.
» Abstract
Reconstruction of ancient climate variability relies on inference from paleoclimate proxy data. However, such data often suffers from large uncertainties in particular concerning the age assigned to measured proxy values, which makes the derivation of clear conclusions challenging. Especially in the study of abrupt climatic shifts, dating uncertainties in the proxy archives merit increased attention, since they frequently happen to be of the same order of magnitude as the dynamics of interest. Yet, analyses of paleoclimate proxy reconstructions tend to focus on mean values and thereby conceal the full range of uncertainty. In addition, the statistical significance of the reported results is sometimes not or at least not accurately tested. Here we discuss both, methods for rigorous propagation of uncertainties and for hypothesis testing with applications to the Dansgaard-Oeschger (DO) events of the last glacial interval and their varying timings in different proxy variables and archives. We scrutinized the mathematical analysis of different paleoclimate records evidencing the DO events and provide results that take into account the full range of uncertainties. We discuss several possibilities of testing the significance of apparent leads and lags between transitions found in proxy data evidencing DO events within and across different ice core archives from Greenland and Antarctica.
D.-D. Rousseau, S. Barbosa, W. Bagniewski, N. Boers, E. Cook, J. Fohlmeister, B. Goswami, N. Marwan, S. O. Rasmussen, L. Sime, A. Svensson:
Data quality in different paleo archives and covering different time scales: A key issue in studying tipping elements,
(Virtual) EGU General Assembly,
Vienna (Austria),
May 4-8, 2020,
DOI:10.5194/egusphere-egu2020-14267,
Talk.
» Abstract
Although the Earth system is described to react relatively abruptly to present anthropogenic forcings, the notion of abruptness remains questionable as it refers to a time scale that is difficult to constrain properly. Recognizing this issue, the tipping elements as listed in Lenton et al. (2008) rely on long-term observations under controlled conditions, which enabled the associated tipping points to be identified. For example, there is evidence nowadays that if the rate of deforestation from forest fires and the climate change does not decrease, the Amazonian forest will reach a tipping point towards savanna (Nobre, 2019), which would impact the regional and global climate systems as well as various other ecosystems, directly or indirectly. However, if the present tipping elements, which are now evidenced, are mostly related to the present climate change and thus directly or indirectly related to anthropogenic forcing, their interpretation must still rely on former cases detected in the past, and especially from studies of abrupt climatic transitions evidenced in paleoclimate proxy records. Moreover, recent studies of past changes have shown that addressing abrupt transitions in the past raises the issue of data quality of individual records, including the precision of the time scale and the quantification of associated uncertainties. Investigating past abrupt transitions and the mechanisms involved requires the best data quality possible. This can be a serious limitation when considering the sparse spatial coverage of high resolution paleo-records where dating is critical and corresponding errors often challenging to control. In theory, this would therefore almost limit our investigations to ice-core records of the last climate cycle, because they offer the best possible time resolution. However, evidence shows that abrupt transitions can also be identified in deeper time with lower resolution records, but still revealing changes or transitions that have impacted the dynamics of the Earth system globally. TiPES Work Package 1 will address these issues and collect paleorecords permitting to evidence the temporal behavior of tipping elements in past climates, including several examples.
K. H. Kraemer, N. Marwan, K. Wiesner, J. Kurths:
Recurrence Plot based entropies and their ability to detect transitions,
(Virtual) EGU General Assembly,
Vienna (Austria),
May 4-8, 2020,
DOI:10.5194/egusphere-egu2020-10861,
Talk.
» Abstract
Many dynamical processes in Earth Sciences are the product of many interacting components and have often limited predictability, not least because they can exhibit regime transitions (e.g. tipping points).To quantify complexity, entropy measures such as the Shannon entropy of the value distribution are widely used. Amongst other more sophisticated ideas, a number of entropy measures based on recurrence plots have been suggested. Because different structures, e.g. diagonal lines, of the recurrence plot are used for the estimation of probabilities, these entropy measures represent different aspects of the analyzed system and, thus, behave differently. In the past, this fact has led to difficulties in interpreting and understanding those measures. We review the definitions, the motivation and interpretation of these entropy measures, compare their differences and discuss some of the pitfalls when using them.
Finally, we illustrate their potential in an application on paleoclimate time series. Using the presented entropy measures, changes and transitions in the climate dynamics in the past can be identified and interpreted.
J. Fohlmeister, N. Bores, N. Marwan, A. Columbu, K. Rehfeld, N. Sekhon, L. Sime, C. Veiga-Pires:
Composite data set of last glacial Dansgaard/Oeschger events obtained from stable oxygen isotopes in speleothems,
(Virtual) EGU General Assembly,
Vienna (Austria),
May 4-8, 2020,
DOI:10.5194/egusphere-egu2020-6894,
Talk.
» Abstract
Millennial scale climate variations called Dansgaard-Oeschger cycles occurred frequently during the last glacial, with their central impact on climate in the North Atlantic region. These events are, for example, well captured by the stable oxygen isotope composition in continental ice from Greenland, but also in records from other regions. Recently, it has been shown that a water isotope enabled general circulation model is able to reproduce those millennial-scale oxygen isotope changes from Greenland (Sime et al., 2019). On a global scale, this isotope-enabled model has not been tested in its performance, as stable oxygen isotope records covering those millennial scale variability were so far missing or not systematically compiled.
In the continental realm, speleothems provide an excellent archive to store the oxygen isotope composition in precipitation during those rapid events. Here, we use a newly established speleothem data base (SISAL, Atsawawaranunt et al., 2018) from which we extracted 126 speleothems, growing in some interval during the last glacial period. We established an automated method for identification of the rapid onsets of interstadials. While the applied method seems to be not sensitive enough to capture all warming events due to the diverse characteristics of speleothem data (temporal resolution, growth stops and dating uncertainties) and low signal-to-noise-ratio, we are confident that our method is not detecting variations in stable oxygen isotopes that do not reflect stadial-interstadial transitions. Finally, all found transitions were stacked for individual speleothem records in order to provide a mean stadial-interstadial transition for various continental locations. This data set could be useful for future comparison of isotope enabled model simulations and corresponding observations, and to test their ability in modelling millennial scale variability.
R. van Dongen, D. Scherler, D. Wendi, E. Deal, C. Meier, N. Marwan, L. Mao:
El Niño-Southern Oscillation (ENSO) controls on mean streamflow and streamflow variability in Central Chile,
(Virtual) EGU General Assembly,
Vienna (Austria),
May 4-8, 2020,
DOI:10.5194/egusphere-egu2020-14011,
Talk.
» Abstract
Understanding hydrological extremes is becoming increasingly important for future adaptation strategies to global warming. Hydrologic extremes affect food security, water resources, natural hazards, and play an important role in the context of erosional processes and landscape evolution. The Pacific region is strongly affected by large-scale climatic anomalies induced by the El Niño-Southern Oscillation (ENSO). How these climatic anomalies translate into hydrological extremes is complex, because both temperature and precipitation deviate from normal conditions and the effect of this simultaneous change on hydrological processes in river catchments (e.g., snowmelt, evapotranspiration) is challenging to understand.
In this study, we investigate the effect of ENSO on mean precipitation, mean temperature, mean stream flow, and streamflow variability in Chile. We have applied extensive quality control on a large hydrological dataset from the Dirección General de Aguas in Chile, resulting in 200 good quality streamflow stations. The dataset envelopes the extent from semi-arid climate in the north ( 28°S) to humid climate in the south ( 42°S). Additionally, the dataset includes low elevation catchments located in the Coastal Cordillera and high elevation catchments in the Andes. We used the monthly Multivariate ENSO Index (MEI) to classify the 5 strongest El Niño and La Niña years, and 5 non-ENSO years after 1975. Changes in mean streamflow and streamflow variability were calculated based on the monitored data from the streamflow stations. For each river catchment, we calculated mean seasonal precipitation using the 0.25°-resolution gridded dataset from the Global Precipitation Climatology Centre (GPCC) and mean seasonal temperature using the 0.5°-resolution global temperature dataset from the Climatology Prediction Centre (CPC).
The precipitation, temperature, and discharge patterns show seasonal variation, varying in strength over the north-south gradient and between low and high elevation catchments. Mean annual precipitation generally increases significantly during El Niño events, and slightly decreases during La Niña events. For both El Niño and La Niña events the mean temperature predominantly changes between 28°S and 35°S and shows increasing temperatures in the Andes and decreasing temperatures in the low elevation Coastal Cordillera. The mean annual streamflow increases during El Niño events, and shows similarities to the pattern of increased mean annual precipitation. However, at the seasonal level, there is a time-lag between precipitation and streamflow, which is regulated by slower snowmelt processes. During La Niña events, the mean annual streamflow increases in the north (28°S-34°S) and decreases in the south (34°S-42°S). Interestingly, the mean annual precipitation and mean annual streamflow patterns behave inversely in the northern Andes. Mean streamflow increases, whereas mean precipitation decreases. This possibly results from enhanced snowmelt because of increased temperatures, but this needs to be further investigated. Finally, the magnitude and frequency of extreme floods predominantly increases in the northern Andean catchments and decreases towards the south for both El Niño and La Niña events. This study shows that large-scale climatic phenomena like ENSO affect catchment hydrology through both anomalies in precipitation and temperature.
W. Duesing, A. Asrat, A. S. Cohen, V. S. Foerster, S. Kaboth-Bahr, K. H. Kraemer, H. F. Lamb, N. Marwan, H. M. Roberts, F. Schaebitz:
Climate beats from Africa: a statistical analysis of the 620 kyr Chew Bahir climate record, eastern Africa,
(Virtual) EGU General Assembly,
Vienna (Austria),
May 4-8, 2020,
DOI:10.5194/egusphere-egu2020-13026,
Talk.
» Abstract
The sediment cores of the Chew Bahir drilling project, part of the Hominin Sites and Paleolakes Drilling Project (HSPDP), from southern Ethiopia, were used to reconstruct climatic changes by analyzing the sediment geochemistry with high-resolution XRF scanning. To interpret the multidimensional XRF dataset we computed a principal component analysis. We used the first principal component (PC1) to detect changes in variability by running a windowed standard deviation analysis and additionally a change point analysis to detect the exact timing of variability changes.
Additionally we used the established Chew Bahir log(K/Zr) aridity proxy, representing clay mineral chemistry-detrital input ratio and compared it to a new Chew Bahir climate indicator, the log(Ca/Ti) proxy, an evaporation signal that is probably inversely related to lake level stands. We find that the log(Ca/Ti) record is also an exceptionally good climate indicator because, compared to the established log(K/Zr) proxy, it reacts with greater amplitude to insolation-controlled signals such as orbital precession. This is confirmed by the log (Ca/Ti) record showing a very clear signal during the African Humid Period, which is however less pronounced in the log(K/Zr) record.
To gain a deeper understanding of the climate cycles and their temporal evolution, we computed a continuous wavelet transformation (CWT) for each of the climate proxies, and studied temporal changes in their cyclicity. Our results indicate that in addition to the precession cycle ( 20 kyr), the Chew Bahir climate record contains earth eccentricity cycles ( 100 kyr), as well as half-precession cycles during high eccentricity. During low eccentricity (450-350 kyr ago), we find reduced variability, three of five changes in standard deviation, damped precession and half precession cycles, and an abrupt transition from dry to wet climate, possibly due to climatic change in high latitudes which may be related to the Mid-Bruhnes event (MBE).
The results confirm that during high eccentricity the tropics are insolation controlled, largely independent of the high latitudes, whereas during low eccentricity the climate of tropical eastern Africa is sensitive to climatic drivers other than precession, possibly originating from high latitudes. Such a period occurring 450 to 350 kyr ago could have led to large regional differences in moisture availability and may have affected early humans by habitat separation, which by isolating populations, resulted in technological diversification. This possible scenario may help to explain the technological transition from Middle Stone Age (MSA) to Acheulean technology that was documented in the Olorgesailie basin during the same time period.
M. H. Trauth, A. Asrat, A. S. Cohen, W. Duesing, V. Foerster, S. Kaboth-Bahr, K. H. Kraemer, H. Lamb, N. Marwan, M. A. Maslin, F. Schaebitz:
Recurrence quantification analysis of the \sim 620 kyr record of climate change from the Chew Bahir basin, southern Ethiopia,
(Virtual) EGU General Assembly,
Vienna (Austria),
May 4-8, 2020,
DOI:10.5194/egusphere-egu2020-4660,
Talk.
» Abstract
The Chew Bahir Drilling Project (CBDP) aims to test possible linkages between climate and evolution in Africa through the analysis of sediment cores that record Quaternary environmental changes in the Chew Bahir basin. In this statistical project we used recurrence plots (PRs) together with a recurrence quantification analysis (RQA) to distinguish two types of variability and transitions in Chew Bahir and compare them with the ODP 967 wetness index from the eastern Mediterranean. The first type of variability are slow variations with cycles of 20 kyr and subharmonics of this cycle. In addition to the these cyclical wet-dry fluctuations in the area, extreme events often occur, i.e. short wet or dry episodes, lasting for several centuries or even millennia, with rapid transitions between wet and dry episodes. The second type of variability is characterized by relatively low variation on orbital time scales, but significant century-to-millennium-scale variations with increasing frequency in the course of an episode of type 2 variability. Within this type of variability there are extremely fast transitions between dry and wet, and vice versa, within a few decades or years, in contrast to those within type 1 which have transitions lasting several hundred years. Type 1 variability probably reflects the influence of precessional forcing in the lower latitudes at times of increased eccentricity, with the tendency towards extreme events, whereas type 2 variability seems to be linked with minimum values of the long (400 kyr) eccentricity cycle, and there does not seem to be a link with atmospheric CO2 levels. The different types of variability and transitions certainly had a completely different influence on the availability of water, food and shelter, and hence eastern Africa's biotic environment, including the habitat of H. sapiens.
M. Kemter, B. Merz, N. Marwan, S. Vorogushyn, G. Blöschl:
Mutual increases in flood extents and magnitudes intensify flood hazard in Central and Western Europe,
(Virtual) EGU General Assembly,
Vienna (Austria),
May 4-8, 2020,
DOI:10.5194/egusphere-egu2020-4731,
Talk.
» Abstract
Climate change has led to changing flood synchrony scales (extents) and flood magnitudes across Europe. We discovered a tight alignment between extents and magnitudes and found the drivers of their joint trends. We analyzed the annual maximum floods of 3872 hydrometric stations across Europe from 1960-2010 and classified all floods in terms of their generating processes based on antecedent weather conditions. There is a positive correlation between flood extents and magnitudes for 95% of the stations. While both parameters increased in Central and Western Europe, they jointly decreased in the East. This widespread magnitude extent correlation is caused by similar correlations for precipitation, soil moisture and snowmelt. We found trends in the relevance of the different flood generation processes, which explain the regional flood trends. The aligned increases of flood extents and magnitudes emphasize the growing importance of transnational flood risk management.
A. Krishnan, R. Manikandan, P. R. Midhun, K. V. Reeja, V. R. Unni, R. I. Sujith, N. Marwan, J. Kurths:
Suppression of oscillatory instability in turbulent reactive flows using network theory,
Complex Dynamical Systems and Applications 2020 (CDSA 2020), Central University of Rajasthan,
Ajmer (India),
February 22, 2020,
Talk.
A. Hartland, B. Goswami, C. Hu, B. Fox, N. Marwan, S. F. M. Breitenbach:
Stalactite flow rates in a central Chinese cave over the last 9000 years,
AGU Fall Meeting,
San Francisco (USA),
December 9–13, 2019,
Talk.
» Abstract
Chinese stalagmites represent unparalleled archives of past East Asian monsoon dynamics. However, the so-called 'amount effect' which links stalagmite oxygen isotope (δ18O) values to monsoon intensity is increasingly questioned. This study introduces a new method to quantitatively reconstruct stalactite discharge (drip rates) to stalagmites, representing a direct hydrological proxy that measures karst aquifer recharge, and by extension effective rainfall. Heshang Cave stalagmite HS4, from central China, grew under a perennial drip point which fluctuates in direct response to annual effective rainfall. The cave system is among the best studied and monitored in the world, providing a strong foundation for calibrating and testing this new proxy. Our analysis shows that drip discharge is highly correlated with stalagmite δ18O, broadly following insolation throughout the Holocene, with centennial-scale variability mainly in phase with δ18O. Flow rates had peak values between 8000-6000 years BP, subsequently declining to minimum flow around 300 years BP. Drip discharge decoupled from δ18O at important intervals in the early, mid and late Holocene. While fluctuations in monsoon rains are clearly coupled to previously identified forcings, we suggest that paleoclimatic drip rates have the potential to redefine our understanding of 'monsoon failure' and to test the drivers embedded in traditional but ambiguous proxies.
N. Marwan:
Datenmanagement in den Geowissenschaften,
Institute of Geoscience, University of Potsdam,
Potsdam (Germany),
November 27, 2019,
Lecture.
» Abstract
Die Vorlesung gibt einen Überblick über die Notwendigkeit nachhaltiger Aufbewahrung wissenschaftlicher Daten, über verschiedene Konzepte zur Datenaufbewahrung, deren Planung und praktischen Umsetzung, sowohl auf persönlicher, institutioneller, als auch auf öffentlicher Ebene in öffentlich zugänglichen Datenarchiven. Konkret werden u. a. Reproduzierbarkeit und Transparenz, wichtige Datenformate, Datenintegrität, Standards, Verschlüsselung, Backup, Coding-Konventionen, Dokumentation/ Meta-Daten und Versionskontrolle besprochen
N. Marwan:
Recurrence Plot Techniques for the Investigation of Recurring Phenomena in the System Earth,
GFZ Earth Surface Seminar Series "Nonlinear Dynamics Workshop",
Potsdam (Germany),
November 11, 2019,
Talk.
N. Marwan:
Recurrence Plot Techniques for the Investigation of Recurring Phenomena in the System Earth,
Habilitation Kolloquium Institute of Geoscience, University of Potsdam,
Potsdam (Germany),
October 16, 2019,
Talk.
» Abstract
Das immer wiederkehrende Auftreten von ähnlichen Zuständen ist eine grundlegende Eigenschaft von Prozessen, die unsere belebte und unbelebte Welt formen und beeinflussen. So gibt es auch zahlreiche Beispiele geologischer und klimatischer Vorgänge sowohl auf kurzen als auch langen Zeit- und Raumskalen, wie dem El Niño-Klimaphänomen, das alle drei bis fünf Jahre auftritt, den Milanković-Zyklen, die in regelmäßigen Abständen zu Eiszeiten führen, die regelmäßige Aktivität von Geysiren, oder das mehr unregelmäßige aber trotzdem wiederkehrende Auftreten von Erdbeben. Die Wiederkehr von Zuständen solch eines dynamischen Prozesses erzeugt ein typisches Wiederkehrmuster, das mit einem Verfahren aus der nichtlinearen Zeitreihenanalyse untersucht werden kann, den sogenannten recurrence plots. In der Habilitation werden fortgeschrittene Aspekte dieser Methodik besprochen. Diese beinhalten
- die Bedeutung der Strukturen und Informationen in recurrence plots,
- die Erweiterung zur Beschreibung räumlich und raumzeitlich wiederkehrender Muster, deren Anwendung zur Identifizierung von
- abrupten Änderungen in der Dynamik und
- äußeren Einflüssen auf die Dynamik eines Systems als auch
- Kopplungen zwischen verschiedenen Systemen,
- eine Kombination mit der Methodik der komplexen Netzwerke,
- Modifikationen zur Behandlung typischer Probleme mit geowissenschaftlichen Daten, wie unregelmäßiges Datensampling und Unsicherheiten in den Daten,
- die Entwicklung eines Signifikanztests und schließlich
- einen Überblick typischer Fehler, die im Zusammenhang mit dieser Methode auftreten können und wie man diese vermeidet.
Neben den methodischen Aspekten werden Anwendungsmöglichkeiten vor allem für geowissenschaftliche Fragestellungen vorgestellt, wie die Analyse von Klimaänderungen, von externen Einflußfaktoren auf ökologische oder klimatische Systeme, oder der Landnutzungsdynamik anhand von Fernerkundungsdaten.
H. Kraemer, N. Marwan:
Border effect corrections for diagonal line based Recurrence Quantification Analysis measures,
16th International Workshop on Complex Systems and Networks,
Berlin (Germany),
September 16, 2019,
Talk.
» Abstract
Recurrence quantification analysis (RQA) is a powerful tool for the identification of characteristic dynamics and of regime changes. This approach is successfully applied in many scientific disciplines. Several measures of complexity are defined on features (such as diagonal and vertical lines) in the recurrence plot (RP) and the corresponding recurrence network (RN). These line structures represent typical dynamical behavior and can be related to certain properties of the dynamical system, e.g., chaotic or periodic dynamics. Therefore, their quantitative study by the RQA measures within sliding windows is a frequently used task for the detection of regime changes. However, as some RQA measures rely on the probability distribution of the lengths of the diagonal lines in an RP, the artificial alteration of these lines due to border effects, insufficient embedding, or a certain sampling setting can have a significant impact on these measures. A few ideas have been suggested to overcome problems. Here we review these ideas, propose novel correction schemes, and systematically compare them. Specifically, we investigate the proper estimation of the diagonal line length entropy for exemplary systems (discrete and continuous). We propose corrections schemes, which yield less biased estimates, especially under noise.
N. Marwan, H. Kraemer, B. Goswami, D. Eroglu, S. Breitenbach, J. Leonhardt:
Quantifiers from Recurrence Plots,
Quantitative paleoenvironments from speleothems workshop, Waikato University,
Cambridge (New Zealand),
September 10, 2019,
Talk.
N. Marwan:
Entropies from Recurrence Plots,
Seminar, Fudan University,
Shanghai (China),
August 26, 2019,
Talk.
N. Marwan, H. Kraemer, K. Wiesner, S. Breitenbach, J. Leonhardt:
Recurrence based entropies,
8th International Symposium on Recurrence Plots,
Zhenjiang (China),
August 21-23, 2019,
Talk.
» Abstract
Dynamical processes in Earth Sciences are generally considered to be of complex nature. The term complexity is frequently used for processes that are either unpredictable (e.g. nonlinear dynamics), consist of many different components, or exhibit regime transitions (e.g. tipping points). To measure complexity, the Shannon entropy is mainly used.
Here we present various entropy measures that have been defined on the base of the recurrence plot. Because of the different features used, these entropy measures represent different aspects of the analysed system and, thus, behave differently. In the past, this fact has lead to difficulties in interpreting and understanding those measures. We summarize the definitions, the motivation and interpretation of these entropy measures, compare their differences and discuss some of the pitfalls when using them.
Finally, we illustrate their potential by applying them on a speleotheme-based palaeoclimate record from Blessberg Cave (Germany). Using entropy measures, the alternating influence of continental versus maritime climate in past central Europe can be identified.
H. Kraemer, N. Marwan:
Border effect corrections for diagonal line based recurrence quantification analysis measures,
8th International Symposium on Recurrence Plots,
Zhenjiang (China),
August 21-23, 2019,
Talk.
» Abstract
Recurrence quantification analysis (RQA) is a powerful tool for the identification of characteristic dynamics and of regime changes. This approach is successfully applied in many scientific disciplines. Several measures of complexity are defined on features (such as diagonal and vertical lines) in the recurrence plot (RP), which represents time points $j$ when a state $\vec{x}_i$ at time $i$ recurs. These line structures represent typical dynamical behavior and can be related to certain properties of the dynamical system, e.g., chaotic or periodic dynamics. Therefore, their quantitative study by the RQA measures within sliding windows is a frequently used task for the detection of regime changes. However, as some RQA measures rely on the probability distribution of the lengths of the diagonal lines in an RP, the artificial alteration of these lines due to border effects, insufficient embedding, or a certain sampling setting can have a significant impact on these measures. A few ideas have been suggested to overcome the mentioned problems. Here we review these ideas, propose novel correction schemes, and systematically compare them.
Specifically, we investigate the proper estimation of the diagonal line length entropy for exemplary systems (discrete and continuous). We propose corrections schemes, which yield less biased estimates, especially under noise.
I. Pavithran, P. Kasthuri, A. Krishnan, S. A. Pawar, R. I. Sujith, R. Gejji, W. E. Anderson, N. Marwan, J. Kurths:
Recurrence networks of spiky signals,
8th International Symposium on Recurrence Plots,
Zhenjiang (China),
August 21-23, 2019,
Talk.
» Abstract
Thermoacoustic instability is a challenging problem faced in gas turbines and rockets. It is a state of self-sustained large amplitude periodic oscillations arising due to the positive feedback between the unsteady heat release rate oscillations and the acoustic field in a confinement. The large amplitude oscillations occurring during thermoacoustic instability are detrimental and their presence can lead to increased heat transfer, violent vibrations causing fatigue failure of the components and even mission failure in rockets. Thermoacoustic system involves interaction between processes at different timescales and length scales leading to complex spatio-temporal dynamics.
In order to study the transitions to such oscillatory instabilities in a complex system, we use a complex networks approach. Complex network is an efficient tool to study systems composed of different interacting entities. Various types of networks have been used in literature such as correlation networks, visibility graphs, recurrence networks, cycle networks, etc. Among these networks, recurrence-based complex networks (RNs) provide information about the topology of the attractor in high dimensional phase space. We can interpret the characteristics of the networks in terms of geometric properties of the phase space.
In the present work, we construct recurrence networks from the time series of acoustic pressure fluctuations obtained during different dynamical regimes from a liquid rocket combustor. We notice that the dynamics of acoustic pressure during thermoacoustic instability is akin to a spiky periodic signal. Such behaviour of the pressure signal is in contrary to the sinusoidal variation observed during thermoacoustic instability in gas turbine combustors. Previous studies have shown that the RN during thermoacoustic instability in gas turbines display a ring-like structure. We unravel a different pattern in RN for a rocket combustor which shows protrusions at different locations on the ring-like structure. These extra patterns are found to be exhibited due to the spiky nature of the time series. We create synthetic time series similar to the experimental data to explain this particular topology. The reconstructed phase space of such largely unexplored signals allows us to get deeper insights about the underlying dynamics of the system. Further, complex networks based on recurrences are an appropriate method for analysing highly nonlinear, high dimensional complex systems.
A. Banerjee, B. Goswami, N. Marwan, B. Merz, J. Kurths:
Recurrence analysis of flood events,
8th International Symposium on Recurrence Plots,
Zhenjiang (China),
August 21-23, 2019,
Poster.
» Abstract
Extreme hydrological events such as floods severely affect the communities living in the corresponding river basins and result in tremendous loss of property and wealth.
The aim of this work is to investigate flood behavior with respect to local effects, e.g. implementation of flood retention basins, and external controls by using recurrence analysis. Flood events occur at irregular time intervals and have a heavy tailed distribution, hence, such data often require data preprocessing and special methods able to analyze data with heavy tailed distribution. In this study, we use the edit distance approach in combination with recurrence plots and recurrence quantification analysis to investigate flood events.
The edit distance approach allows us to use the recurrence-based characteristics to quantify how the dynamics of the flood occurrence has changed over time. We apply our approach to the river discharge data from the river Elbe and study the dynamical interactions of different variables such as precipitation, temperature and catchment wetness.
H. Kraemer, R. Donner, J. Heitzig, N. Marwan:
Recurrence threshold selection for obtaining robust recurrence characteristics in different embedding dimensions,
8th International Symposium on Recurrence Plots,
Zhenjiang (China),
August 21-23, 2019,
Poster.
» Abstract
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.
N. Marwan:
Co-authorship network of the recurrence plot domain,
8th International Symposium on Recurrence Plots,
Zhenjiang (China),
August 21-23, 2019,
Poster.
» Abstract
A presentation of the co-authorship network of publications on recurrence plots, recurrence networks and recurrence quantification analysis.
N. Marwan:
Karstgebiet Sägistal,
14. Nationaler Höhlenforscher-Kongress Sinterlaken19,
Interlaken (Switzerland),
August 9-11, 2019,
Talk.
» Abstract
Das Sägistal ist ein abgelegenes Hochtal der Berner Voralpen mit typischen Karsterscheinungen. Die Erforschung der Höhlen begann in den 1970er Jahren durch die SGH Interlaken und wird seit 1988 durch die Internationale Speläologische Arbeitsgruppe Alpiner Karst (ISAAK) unter Beteiligung zahlreicher Höhlenforschergruppen aus verschiedenen Ländern organisiert. Mittlerweile wurden über 400 Höhlen gefunden mit dem "Oberländer-Chessiloch"-System als größtem Objekt (2346 m Länge, -488 m Tiefe).
D. Wendi, B. Merz, N. Marwan:
Novel quantification method for hydrograph similarity,
SimHydro 2019,
Sophia Antipolis (France),
June 12, 2019,
Talk.
» Abstract
We propose an additional elaborate hydrological signature index to quantify similarity (and dissimilarity) between recurring flood dynamics and between observation and model simulation as implied by their phase space trajectories. These phase space trajectories are reconstructed from their corresponding hydrographs (i.e., event time series) using Taken's time delay embedding method. This reconstructed phase space allows multi-dimensional relationship between observation points (i.e., at different time of the event) to be analyzed. Such approach considers the relationships of set of magnitude points in their unique time sequence that are relevant to the complex temporal cascading processes in flood. In a simpler terms, the new index considers the characteristics shape dynamics of a hydrograph and optionally the antecedent discharge conditions that may implicitly cascade to the subsequent rainfall-runoff event and cause an extreme or unusual hydrograph shape. This new similarity index can be used to comprehensively assess the recurrence of extreme event characteristics, change of flood dynamics, shift of seasonality, and as additional metric or objective function to evaluate and calibrate hydrological and hydraulics models.
N. Marwan:
Recurrence Plots for Time Series Analysis,
Geological Remote Sensing Seminar, University of Potsdam,
Potsdam (Germany),
June 4, 2019,
Lecture.
U. Ozturk, N. Malik, K. Cheung, N. Marwan, J. Kurths:
Tracking tropical and frontal storms driven extreme rainfalls over Japan using complex networks,
SIAM Workshop on Network Science 2019,
Snowbird (USA),
May 22-23, 2019,
Talk.
» Abstract
Predicting extreme rainfall is a challenging but a necessary task due to the concomitant natural hazards, such as flash floods or landslides. Theory of network science offers alternative tools to explore the spatiotemporal properties of extreme rainfall, which might reveal the predictive behavior of those extremes. In this case study, we use complex network metrics in conjunction with a nonlinear correlation measures of event synchronization to study extreme rainfall generated by the tropical and frontal (Baiu) storms over Japan. These two weather systems trigger extreme events in two discrete seasons; the Baiu front dominates the rainfall events from June to July, whereas tropical storms activity peak at August, and are active until November. We found that the spatial scales involved in the Baiu driven rainfall extremes are consistently more extensive than the extremes due to tropical storms. We further delineate an east-west extending horizontal region of coherent rainfall during Baiu season based on network communities, whereas nearly all Japan fall in one single coherent rainfall community during tropical storm season.
N. Marwan, H. Kraemer, K. Wiesner, S. Breitenbach, J. Leonhardt:
Recurrence based entropies,
Fourth International Conference on Recent Advances in Nonlinear Mechanics,
Łódz (Poland),
May 7-10, 2019,
Talk.
N. Marwan, H. Kraemer, K. Wiesner, S. Breitenbach, J. Leonhardt:
Recurrence based entropies,
EGU General Assembly,
Vienna (Austria),
April 7-12, 2019,
» Poster (PDF, 6.42M)
.
» Abstract
Dynamical processes in Earth sciences are often considered to be of complex nature. The term complexity is often used for processes that are either unpredictable (e.g. nonlinear dynamics), consist of many different components, or exhibit regime transitions (e.g. tipping points). To measure complexity, the Shannon entropy is often used. Here we present various entropy measures that have been defined on the base of the recurrence plot. Because of the different features that are used, these entropy measures represent different aspects of the analysed system and, thus, behave differently. In the past, this fact has lead to difficulties in interpreting and understanding those measures. We summarize the definitions, the motivation and interpretation of these entropy measures, compare their differences and discuss some of the pitfalls when using them.
Finally, we illustrate their potential in an application on palaeoclimate time series. Using entropy measures, changes and transitions in the climate dynamics in the past can be identified and interpreted.
W. Düsing, H. Kraemer, A. Asrat, M. Chapot, A. Cohen, A. Deino, V. Foerster, H. Lamb, N. Marwan, C. Lane, M. Maslin, C. Ramsey, H. Roberts, F. Schaebitz, M. Trauth, C. Vidal:
Differentiating local from regional climate signals using the 600 ka Chew Bahir paleoclimate record from South Ethiopia,
EGU General Assembly,
Vienna (Austria),
April 7-12, 2019,
Talk.
» Abstract
Cores from terrestrial archives, such as the lacustrine sediments from the Chew Bahir basin in southern Ethiopia, which cover the last 600 ka, often reflect both local, regional and global climate influences. In our analysis we were able to identify several time windows in which the Chew Bahir climate is in resonance with regional and global climate change.
As a contribution to understanding and differentiating these connections recorded in the Chew Bahir sediments, we have correlated the 2nd principal component of the MSCL based color reflectance values representing wet conditions in the Chew Bahir basin, with the wetness index from ocean core ODP 967 from the eastern Mediterranean Sea. The correlation between these two time series was calculated using the Spearman correlation coefficient in a sliding window. Episodes with high correlation between the two records of wetness could indicate a strong link between both regions, possibly through an increased outflow of the river Nile into the eastern Mediterranean Sea due to higher precipitation values on the Ethiopian plateau.
Our preliminary results show that when correlating the two records, two distinct temporal units can be distinguished. Between 570 ka and 350 ka the correlation is dominated by cycles that correspond with orbital precession whereas the second unit (after 350 ka) reveals a strong influence of atmospheric CO2. This observation suggests that both orbital precession and atmospheric CO2. may cause a synchronization of different regions in the African climate system, possibly depending on boundary conditions which are still to be identified.
As a next step we'll investigate the nonlinear relationships between the two records by focusing on the transition between the two main observed phases. The transition around 350 kyrs however, is not only highly interesting from a climatic perspective, but it is also a noteworthy period for human cultural evolution as a transition from Acheulean to Middle Stone Age (MSA) technologies takes place at this time. So far our results outline that during this climatically and evolutionary relevant episode a relatively stable, long-lasting, pan-African wet phase may have existed, with possible green corridors connecting the habitats of hominins, and ample resources supporting large population sizes.
A. Banerjee, B. Goswami, N. Marwan, B. Merz, J. Kurths:
Recurrence Analysis of Flood Events,
EGU General Assembly,
Vienna (Austria),
April 7-12, 2019,
Poster.
» Abstract
Extreme hydrological events such as floods severely affect the communities living in the corresponding river basins and result in tremendous loss of property and wealth.
The aim of this work is to investigate flood behavior with respect to local effects, e.g. implementation of flood retention basins, and external controls by using recurrence analysis. Flood events occur at irregular time intervals and hence such data often requires data preprocessing. In this study, we use the TACTS approach in combination with recurrence plot and recurrence quantification analysis to investigate flood events, which occur on irregular time scales.
The TACTS approach allows us to construct a recurrence plot from the irregularly spaced flood event series, and the recurrence-based characteristics help to quantify how the dynamics of the flood occurrence have changed over time. We apply our approach to the 150-year river discharge data from the river Elbe and study the dynamic interactions of different variables such as precipitation, temperature and catchment wetness.
A. Agarwal, L. Caesar, N. Marwan, R. Maheswaran, B. Merz, J. Kurths:
Detection of short- and long-range teleconnections in SST patterns on different time scales,
EGU General Assembly,
Vienna (Austria),
April 7-12, 2019,
Poster.
» Abstract
Sea surface temperature (SST) anomaly patterns can - as surface climate forcing -affect the weather at large distances. This is why following an El Niño event major global climate anomalies occur. This paper characterizes the links between the cells of a global SST grid data set at different temporal and spatial scales with the help of climate networks. These networks are constructed using wavelet multi-scale correlation. This way we identify and visualise the SST patterns that develop very similarly over time and distinguish them from those that have long-range teleconnections to other ocean regions. Our findings re-confirm accepted knowledge about known highly linked SST patterns like El Niño Southern Oscillation and the Pacific Decadal Oscillation, but also suggest new insights into the characteristics and origins of long-range teleconnections.
M. Kemter, B. Merz, N. Marwan:
Using multi-layer complex networks to understand interrelationships and changes in extreme flood generation,
EGU General Assembly,
Vienna (Austria),
April 7-12, 2019,
Poster.
» Abstract
The generation of extreme flood events is influenced by a multitude of parameters that interact in complex ways. To understand their temporal and spatial relationships as well as changes in this system we need adequate tools. We therefore use multilayer complex networks and extreme event statistics to discover and interpret relationships between flood influencing parameters (e.g. precipitation, catchment wetness, discharge). Complex networks have formerly been successfully used for climatic and hydrological representations. A multilayer approach enables us to find inter-relationships between the different influences. We use non-linear similarity measures to generate the network connections. By analysis of variations of the network appearance and metrics with time, we can reconstruct temporal changes in the underlying processes.We are investigating several hundred river gauges across Europe over a timeframe of 70 years.
M. Trauth, A. Asrat, C. Bronk Ramsey, M. Chapot, A. Cohen, A. Deino, W. Duesing, V. Foerster, H. Kraemer, H. Lamb, C. Lane, N. Marwan, M. Maslin, H. Roberts, F. Schaebitz, C. Vidal:
Recurring types of variability and transitions in the \~280 m long (\~600 kyr) sediment core from the Chew Bahir basin, southern Ethiopia,
EGU General Assembly,
Vienna (Austria),
April 7-12, 2019,
Poster.
» Abstract
The Chew Bahir Drilling Project (CBDP) aims to test hypothesized linkages between climate and human evolution, dispersal and technological innovation by the acquisition and analysis of long (280 m) sediment cores that have recorded environmental change in the Chew Bahir basin, southern Ethiopia. In this time-series analysis 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, and humans within the lake catchment. 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.
We use methods of linear and nonlinear time series analysis, such as change point detection, semblance analysis and recurrence plots, to identify and classify recurring types of variability and transitions on the time scales of human life spans. For example, we investigate the rapidness of transitions, possible precursor events, and tipping points in our palaeoenvironmental data and discuss their possible impact on the living conditions of humans in the region. First results of the analysis show that we indeed find, as an example, recurring threshold-type transitions, when the Chew Bahir system switched from one stable mode to another, such as from stable wet to dry conditions. Such a rapid change of climate in response to a relatively modest change in forcing appears to be typical of tipping points in complex systems such as the Chew Bahir. If this is the case then the 14 dry events idenfified at the end of the African Humid Period (15-5 kyr BP) could represent precursors of an imminent tipping point that, if properly interpreted, would allow predictions to be made of future climate change in the Chew Bahir basin.
N. Marwan:
Recurrence Plots for Time Series Analysis,
Seminar at MPI Molecular Physiology,
Dortmund (Germany),
March 18, 2019,
Talk.
N. Marwan:
Transparent and efficient data storage,
Graduate School NatRiskChange University of Potsdam,
Potsdam (Germany),
March 13-14, 2019,
Lecture.
» Abstract
The lecture provides an overview of the need for sustainable storage of scientific data, various concepts of data storage and archiving, their planning and practical implementation, both at the personal, institutional, and public levels in publicly accessible data archives. Specific topics discussed include reproducibility and transparency, important data formats, data integrity, standards, encryption, backup, coding conventions, documentation/meta-data, and version control.
W. Duesing, J. Thom, A. Deino, A. Asrat, V. Foerster, H. Kraemer, N. Marwan, H. Lamb, F. Schaebitz, M. H. Trauth:
Human evolution and climate change: What can we learn from the 0-630 kyrs BP paleoclimate record from the Chew Bahir basin in eastern Africa?,
AGU Fall Meeting,
Washington DC (USA),
December 10-14, 2018,
Talk.
» Abstract
Numerous authors have developed hypotheses linking climate, the environment and human evolution, expansion and technological innovation in eastern Africa (e.g., Potts, 2013, Potts et al. 2018; Maslin et al., 2015). The Chew Bahir Drilling Project (CBDP) as part of the Hominin Sites Paleolakes Drilling Project (HSPDP) aims to test some of these hypotheses by providing a long, continuous and high-resolution paleoclimate records of climatic and environmental change through critical intervals of human evolution.
Testing such hypothesis requires a very accurate age model both for the paleoclimate record, but also for the archeological/anthropological evidence. We therefore developed a MATLAB-based Multiband Wavelet Based Age Modeling Technique (mubawa), which generates an orbital tuned age model that comprises uncertainties. Next, we use a piecewise correlation using a set of different sliding windows to compare the Chew Bahir paleo records on different time scales with Indian Ocean SSTs, Terrestrial Dust and Nile-Outflow records. To classify variability and transitions we use recurrence plots/recurrence quantification analysis. Application of this method detects nonlinear features such as tipping points, which are often coherent with changes in variability and strong precursor events. The recurrences of such precursor events often lie within the life span of hominins. For individuals living during that time these drastic climate shift were perceptive and probably provoked new survival strategies that may have preserved in the archeological record.
In a last step we evaluate if these detected climate and environmental shifts can be related to archeological sights, human evolution, technological innovation, migration and dispersal events.
E. Macau, A. M. Ramos, J. Kurths, N. Marwan:
Detecting causal relations from real data experiments by using recurrence,
Dynamics Days Latin America and the Caribbean 2018,
Punta del Este (Uruguay),
November 26-30, 2018,
Talk.
» Abstract
In this work, we present the Recurrence Measure of Conditional Dependence (RMCD), a recent data-driven causality inference method using the framework of recurrence plots. The RMCD incorporates the recurrence behavior into the transfer entropy theory. We apply this methodology to some paradigmatic models and to investigate the possible influence of the Pacific Ocean temperatures on the South West Amazon for the 2010 and 2005 droughts. The results reveal that for the 2005 drought there is not a significant signal of dependence from the Pacific Ocean and that for 2010 there is a signal of dependence of around 200 days. These outcomes are confirmed by the traditional climatological analysis of these episodes available in the literature and show the accuracy of RMCD inferring causal relations in climate systems.
T. Vantuch, I. Zelinka, A. Adamatzky, N. Marwan:
Detecting causal relations from real data experiments by using recurrence,
17th International Conference on Unconventional Computation and Natural Computation,
Fontainebleau (France),
June 26, 2018,
Talk.
» Abstract
Natural systems often exhibit chaotic behavior in their space-time evolution. Systems transiting between chaos and order manifest a potential to compute, as shown with cellular automata and artificial neural networks. We demonstrate that swarms optimisation algorithms also exhibit transitions from chaos, analogous to motion of gas molecules, when particles explore solution space disorderly, to order, when particles follow a leader, similar to molecules propagating along diffusion gradients in liquid solutions of reagents. We analyse these ‘phase-like’ transitions in swarm optimization algorithms using recurrence quantification analysis and Lempel-Ziv complexity estimation. We demonstrate that converging and non-converging iterations of the optimization algorithms are statistically different in a view of applied chaos, complexity and predictability estimating indicators.
W. Duesing, A. Asrat, V. E. Foerster, H. Kraemer, H. F. Lamb, N. Marwan, F. Schaebitz, M. H. Trauth:
Trends, rhythms and transitions during the Late Quaternary in southern Ethiopia,
EGU General Assembly,
Vienna (Austria),
April 8-13, 2018,
Poster.
» Abstract
This project aims at statistically analyzing the long ( 278 m) sediment record of the Chew Bahir basin, as part of the ICDP-funded Hominin Sites and Paleolakes Drilling Project (HSPDP). The aim of the project is (1) to establish a robust age-depth model for the sediment cores, (2) to correlate the Chew Bahir record with other records within and outside HSPDP, (3) to detect trends, rhythms and events in the environmental record of the basin, and (4) identify recurrent, characteristic types of climate transitions in the time series, as compared with the ones of the other HSPDP sites and climate records outside HSPDP. The work presented here will provide first results of age-depth modelling, including cyclostratigraphy, of the long Chew Bahir cores. Second, it gives an overview of the first results from evolutionary spectral analysis to detect changes in the response of the Chew Bahir to orbital forcing during the last 550 kyr. Third, the results of a change point analysis will be presented to define the amplitude and duration of past climate transitions and their possible influence on the development of early modern human cultures.
M. H. Trauth, A. Asrat, W. Duesing, V. Foerster, H. Kraemer, H. Lamb, N. Marwan, M. A. Maslin, F. Schaebitz::
Classifying past climate variation in the Chew Bahir basin, southern Ethiopia, using recurrence quantification analysis,
EGU General Assembly,
Vienna (Austria),
April 8-13, 2018,
Poster.
» Abstract
The Chew Bahir Drilling Project (CBDP) aims to test hypothesized linkages between climate and mammalian (including hominin) evolution in tropical-subtropical eastern Africa by the acquisition and analysis of long ( 280 m) sediment cores that have recorded environmental change in the Chew Bahir basin. In our statistical project, we describe the Chew Bahir paleolake as a dynamical system composed of interacting components, such as the water body, the sediment below the bottom of the (paleo-)lake, and the organisms living in the lake and its surroundings. A common feature of dynamical systems is the property of recurrence, where patterns of recurring states reflect typical system characteristics whose description contribute significantly to understanding its dynamics. In our example it could be a recurrence of changes in the state variables precipitation, evaporation and wind speed, which lead to similar (but not identical) conditions in the lake (e.g., depth and size of the lake, alkalinity and salinity of the lake water, species assemblage in the water body, diagenesis in the sediment). A recurrence plot (RP), first introduced by J.P. Eckmann in 1987, is a graphical display of such recurring states of the system, calculated from the distance (e.g. Euclidean) between all pairs of observations x(t), within a cutoff limit. To complement the visual inspection of recurrence plots, measures of complexity were introduced for their quantitative description to perform the recurrence quantification analysis (RQA). Here we present and discuss preliminary results of a RQA of the 550 kyr long environmental record from the Chew Bahir basin.
A. Agarwal, N. Marwan, M. Rathinasamy, U. Ozturk, B. Merz, J. Kurths:
Complex network-based approach for identification of influential and expandable station across rainfall network,
EGU General Assembly,
Vienna (Austria),
April 8-13, 2018,
Poster.
» Abstract
The complex network has gained significant momentum in last decades and has found application wide areas ranging from biological networks to climate networks. In analysing physical complex networks, identification of key influential nodes is l an important field of research. In this study, we propose a new effective node ranking method based on network measure degree and betweenness values. The proposed method is tested and compared to previously proposed node ranking methods on synthetic sample networks and then applied to a real-world raingauge network of 1229 stations from Germany to check its replicability and applicability. Raingauge networks play a vital role in providing information for making crucial decisions in water resources management and resources estimation. The network of operating raingauges should be set up optimally to provide as much and as accurate information as possible and at the same time cost-effective. The proposed method is evaluated using decline rate efficiency and kriging error. The results of the study show that the proposed method based on complex network theory for ranking the raingauges is robust and can be used for design and redesign the raingauge network. The method is very useful in identifying the highly influential station which needs high attention and expendable stations which either can be relocated, uninstalled or removed without much effect on the overall accuracy of the observations provided by the raingauge network.
D. Wendi, N. Marwan, B. Merz:
The recurrence of unseasonable and rare flood dynamics,
EGU General Assembly,
Vienna (Austria),
April 8-13, 2018,
» Poster (PDF, 2.26M)
.
» Abstract
The question whether a certain flood is rare/ unusual or not is often evaluated from the frequency curve analysis (also called growth curve) of flood discharge peaks and their corresponding recurrence interval (return period). Flood discharge peak and maximum depth are some of the hydrological signatures (i.e. an element of hydrograph) and are popular choice for flood risk assessment due to their close relationship with socio-economic impact of flood and often used for damage modeling. However what if our question now is whether the flood process dynamics is rare (influenced by unusual or more driving mechanism, e.g. ice jam, dam break, clogged drain, etc) and especially if hydrological boundary condition is no longer the same as before.
The confinement of flood peaks reduces the information about the flood dynamics inferred by the shape of hydrograph, especially should one be interested to evaluate a rarity/ extremity of a flood process dynamics since flood peak is just one element of a flood hydrograph. Although other indices derived from hydrological signature (e.g. volume, slopes, base flow index, etc) are useful descriptors of a process dynamics, most of them are still either just a part of hydrograph, or derived as an aggregate (e.g. slopes and volume) and therefore unable to provide bigger picture of the flood dynamics and suffer from statistical uncertainty. Furthermore, with singular descriptor from the mentioned, different flood dynamics (i.e. resulted from different processes/ boundary conditions) could be mistaken as the same and might lead to misinterpretation (e.g. snow melt and rainfall triggered runoff may easily share similar flood peak and volume). Moreover, stationarity in season is often assumed in flood frequency analysis, that different flood processes are classified to follow strict calendar month seasons. Such practice could fail to analyze the occurrence of unusual climatic event such as early or late snow melt and unseasonably heavy rainfall in winter.
In this study, we focus on the utilization of hydrological signature to characterize temporal flood event dynamics with the objective to analyze their recurrences and to be able to evaluate if a process dynamic of a certain flood is rare or perhaps unprecedented. We propose using the analysis of phase space trajectories reconstructed through time delay embedding of a time series to characterize different flood events. To allow the visualization and analysis of high embedding dimension (i.e. above 3), we suggest the use of recurrence plot (RP) and quantification (RQA) as similarity measures between the flood dynamics of one event to another and allow non-stationarity occurrence of their typology.
B. Goswami, S. Breitenbach, F. Lechleitner, J. Baldini, H. Cheng, N. Marwan:
Is this an event? – Detecting abrupt changes in palaeoclimate records,
EGU General Assembly,
Vienna (Austria),
April 8-13, 2018,
» Poster (PDF, 772.16K)
.
» Abstract
Abrupt shifts in a certain climate state is a pertinent question in palaeoclimate studies and is crucial for determining leads or lags between spatially disperse observations. Such events have also a direct bearing on the vulnerability of society to drastic – and possibly difficult to mitigate – changes. Determining whether or not, and when exactly, abrupt changes in climate occurred is made challenging by temporal resolution and uncertainties associated with determining the age of climate proxy measurements.
In this study, we present a robust, ‘uncertainty-aware’ approach to determine periods of abrupt change from palaeoclimate proxy records. Our method is based on a new representation of time series and it utilises the recurrence properties of the proxy record to distinguish time points of abrupt climate change. We first validate our approach with a synthetic example, and thereafter, we apply our approach to speleothem records from China and India. Our results reveal a highly non-trivial spatio-temporal pattern of the detected events in the Asian monsoon domain.
U. Öztürk, N. Malik, N. Marwan, J. Kurths:
Comparison of tropical and frontal storms using complex networks,
EGU General Assembly,
Vienna (Austria),
April 8-13, 2018,
Poster.
» Abstract
Complex network analysis supports exploring spatiotemporal dynamics of significant climate phenomena, such as heavy precipitation. Complex networks are able to capture the spreading and concentration of extreme rainfall by using event-synchronization, such as rainfall propagation patterns of the Indian Summer Monsoon by parameterizing the delay from precipitation time-series. Despite much advancement in monitoring extreme rainfall, capturing spatiotemporal dynamics of the fast-evolving atmospheric events (e.g., tropical storms) is still a challenge. Quantifying spatial scales of extreme rainfall will aid mitigating concomitant flood and landslide hazards.
We use network analysis to compare spatial features of extreme rainfall over Japan using satellite-derived rainfall data (TRMM-3B42V7). We first divide the time series into two subsets: June to July (JJ) and August to November (ASON) to concentrate on the Baiu front season (JJ) and the tropical storms season (ASON). We assess the spatial scales involved in the two distinct mechanisms and define regions of coherent rainfall during the two seasons. We additionally propose using radial statistics to trace the network flux over long distances, which allows us to observe the general pattern of extreme rainfall tracks. Extreme rainfall associated with tropical storms show smaller spatial scales (in the range of 100 km) compared to Baiu linked extremes. We also discovered a consistent deviation of the extreme rainfall from the eye of the tropical storm tracks.
H. Kraemer, R. V. Donner, N. Marwan, M. H. Trauth:
Detecting abrupt transitions during the Late Quaternary in southern Ethiopia using Recurrence Quantification Analyses,
EGU General Assembly,
Vienna (Austria),
April 8-13, 2018,
Poster.
» Abstract
In many data driven fields of research, categorizing abrupt transitions / regime changes is of high interest. The different aspects of temporal recurrence patterns of previous states can help to identify and characterize subtle changes in systems dynamics. Besides the identification of transitions, recurrence methods can provide a better understanding of the process underlying these transitions by statistically describing the dynamical characteristics, e.g. the predictability, determinism and complexity of the dynamical system. For example, the characteristic block structures in the recurrence plot can be used to identify different types of intermittency. In general, changes between different dynamical regimes are visually well expressed in recurrence plots. The introduction of selected recurrence quantifiers (such as recurrence rate, determinism, or laminarity) together with a running window approach has paved the way for a quantitative recurrence analysis of transitions and therefore should be able to provide a classification of different transition types.
In order to achieve such a classification there is necessity for developing a method which is capable to statistically analyze the behavior of recurrence quantifiers at transitions. In this work, we show how to make statements about the significance of estimated values of recurrence quantifiers using a bootstrap approach. We also highlight the specific technical problems related to that task. The presented method also allows gaining information about the duration of a transition. Here we demonstrate potentials of the proposed approach to detect abrupt transitions in (1) prototypical models of transitions as well as in (2) real data of past climate variations in the Chew Bahir basin (South Ethiopia), investigated within the Hominin Sites and Paleolakes Drilling Project (HSPDP).
B. Goswami, N. Boers, A. Rheinwalt, N. Marwan, J. Heitzig, S. Breitenbach, J. Kurths:
Identifying sudden dynamical shifts in time series with uncertainties,
EGU General Assembly,
Vienna (Austria),
April 8-13, 2018,
Poster.
S. Breitenbach, B. Plessen, S. Waltgenbach, R. Tjallingii, J. Leonhardt, K.-P. Jochum, H. Meyer, N. Marwan, D. Scholz:
Tracing past shifts of the boundary between maritime and continental climate over Central Europe,
EGU General Assembly,
Vienna (Austria),
April 8-13, 2018,
Poster.
» Abstract
European climate is characterized by heterogeneous climate conditions, with distinct boundaries between zones that can be classified according to the Köppen classification (Peel et al. 2007), and detected using climate network techniques (Rheinwalt et al. 2016). These boundaries are not stationary, but shift geographically, depending on large scale atmospheric conditions.
Central European climate is strongly influenced by intricately linked North Atlantic Oscillation and Siberian High (SH), which govern precipitation and temperature over Europe. Shifts of these climatic boundaries in response to global warming and circulation changes might lead to more frequent extreme weather patterns like heat waves, with significant repercussions for society (Cohen et al. 2014).
Speleothem-based palaeoclimate reconstructions enable us to understand underlying forcing mechanisms and speed of climatic reorganizations. Here we present a first reconstruction of multi-centennial shifts of the boundary between western European maritime Cfb climate and continental Dfb climate through the last ca. 5,000 years using speleothems from Bleßberg Cave, Thuringia, Central Europe.
Thanks to its location near the Cfb-Dfb climatic boundary, Bleßberg Cave is ideally suited to reconstruct past W-E shifts of this divide longitudinally crossing Central Europe. We compare a decadally resolved stalagmite δ18O record with data from Bunker Cave (Mischel et al. 2017), western Germany, and an NAO reconstruction from Greenland (Olsen et al. 2012).
Over the last 5,000 years, the boundary between Cfb and Dfb climate shifted repeatedly. When the Cfb-Dfb border was east (west) of Bleßberg (Bunker) Cave maritime (continental) climate prevailed at both sites. Discrepancies between investigated proxy records are found when the boundary is located between the two caves. Comparison with the Greenland NAO record shows that a westerly shifted boundary is often associated with a strong SH and a negative NAO. An easterly shift, in contrast, is found to be linked with weak a SH and a positive NAO.
D. Wendi, N. Marwan, B. Merz:
The importance of hydrological signature and its recurring dynamics,
AGU Fall Meeting,
New Orleans (USA),
December 1–17, 2017,
Poster.
» Abstract
Temporal changes in hydrology are known to be challenging to detect and attribute due to multiple drivers that include complex processes that are non-stationary and highly variable. These drivers, such as human-induced climate change, natural climate variability, implementation of flood defense, river training, and land use change, could impact variably on space-time scales and influence or mask each other. Besides, data depicting these drivers are often not available.
One conventional approach of analyzing the change is based on discrete points of magnitude (e.g. the frequency of recurring extreme discharge) and often linearly quantified and hence do not reveal the potential change in the hydrological process. Moreover, discharge series are often subject to measurement errors, such as rating curve error especially in the case of flood peaks where observation are derived through extrapolation.
In this study, the system dynamics inferred from the hydrological signature (i.e. the shape of hydrograph) is being emphasized. One example is to see if certain flood dynamics (instead of flood peak) in the recent years, had also occurred in the past (or rather extraordinary), and if so what is its recurring rate and if there had been a shift in its occurrence in time or seasonality (e.g. earlier snow melt dominant flood). The utilization of hydrological signature here is extended beyond those of classical hydrology such as base flow index, recession and rising limb slope, and time to peak. It is in fact all these characteristics combined i.e. from the start until the end of the hydrograph. Recurrence plot is used as a method to quantify and visualize the recurring hydrological signature through its phase space trajectories, and usually in the order of dimension above 2. Such phase space trajectories are constructed by embedding the time series into a series of variables (i.e. number of dimension) corresponding to the time delay. Since the method is rather novel in hydrological community, the study presents an overview and a guideline to the method with an application example on analyzing the change of hydrological signature and discussion of its benefits and flaws.
N. Marwan:
Recurrence Plots for Data Analysis,
QUEST Workshop on palaeoclimate time series analysis and statistics,
Potsdam (Germany),
November 3, 2017,
Lecture and workshop.
N. Marwan:
Transparent and efficient data storage,
QUEST Workshop on palaeoclimate time series analysis and statistics,
Potsdam (Germany),
November 3, 2017,
Lecture.
N. Marwan:
Höhlen als wissenschaftliche Archive,
20th Anniversary of Speleo Club Berlin,
Kienitz (Germany),
September 23, 2017,
Talk.
N. Marwan:
What can you see? Perception problems of recurrence plots,
7th International Symposium on Recurrence Plots,
São Paulo (Brazil),
August 23-25, 2017,
» Talk (PDF, 3.44M)
.
» Abstract
Recurrence plots exhibit features and patterns which are characteristic for typical dynamics. How difficult does it be to visually recognize the dynamics from the recurrence plot? Does everybody see the same or judges the different patterns with similar importance? In this (interactive) talk we will figure out the subjective nature of visual inspection and discuss the difficulties. It finally underlines the importance of applying objective quantifiers such as recurrence quantification analysis.
D. Eroglu, N. Marwan:
Multiplex Recurrence Networks,
7th International Symposium on Recurrence Plots,
São Paulo (Brazil),
August 23-25, 2017,
Poster.
» Abstract
The complex nature of a variety of phenomena in physical, biological, or earth sciences is driven by a large number of degrees of freedom which are strongly interconnected. Although the evolution of such systems is described by multivariate time series (MTS), so far research mostly focuses on analyzing these components one by one.
Recurrence based analyses are powerful methods to understand the underlying dynamics of a dynamical system and have been used for many successful applications including examples from earth science, economics, or chemical reactions. The backbone of these techniques is creating the phase space of the system. However, increasing the dimension of a system requires increasing the length of the time series in order get significant and reliable results. This requirement is one of the challenges in many disciplines, in particular in palaeoclimate, thus, it is not easy to create a phase space from measured MTS due to the limited number of available obervations (samples). To overcome this problem, we suggest to create recurrence networks from each component of the system and combine them into a multiplex network structure, the it multiplex recurrence network (MRN). We test the MRN by using prototypical mathematical models and demonstrate its use by studying high-dimensional palaeoclimate dynamics derived from pollen data from the Bear Lake (Utah, US). By using the MRN, we can distinguish typical climate transition events, e.g., such between Marine Isotope Stages.
H. Kraemer, N. Marwan, M. H. Trauth:
Classifying abrupt transitions in IPCC climate models and paleoclimate proxy data using recurrence quantification analysis,
7th International Symposium on Recurrence Plots,
São Paulo (Brazil),
August 23-25, 2017,
Talk.
» Abstract
In many data driven disciplines, categorising abrupt transitions / regime changes are of high interest. The different aspects of recurrence can help to identify and characterize subtle changes in systems dynamics. Besides the identification of transitions, recurrence methods can help to provide a better understanding of the underlying process of these transitions by statistically describing the dynamical characteristics, e.g. the predictability, determinism and complexity of the dynamical system. For example, the characteristic block structures in the recurrence plot can be used to identify different types of intermittency. In general, changes between different dynamics are visually well expressed in recurrence plots. The introduction of selected recurrence quantifiers (such as recurrence rate, determinism, or laminarity) together with a running window approach has paved the way for a quantitative recurrence analysis of transitions and therefore allow a classification of different transition types.
In this work first results of such recurrence based classification is shown. We demonstrate it by analysing prototypical models of transitions as well as on real world data related to palaeoclimate. The prototypical models are selected from a catalogue of transition types which have been used and discussed in models presented in the reports of the Intergovernmental Panel on Climate Change (IPCC)[1]. In the palaeoclimate example we consider two Potassium time series of two drilling cores from the Chew Bahir Bassin, which is part of the Hominin Sites and Paleolakes Drilling Project (HSPDP).
D. Wendi, N. Marwan, B. Merz, J. Kurths:
Change in flood hazard dynamics from recurrence perspective,
7th International Symposium on Recurrence Plots,
São Paulo (Brazil),
August 23-25, 2017,
Talk.
» Abstract
Temporal changes in flood hazard systems are known to be difficult to detect and attribute due to multiple drivers that include processes that are non-stationary and highly variable. Often such analysis of change is quantified from single points perspective (i.e. extreme values) that may subject to high errors and uncertainties. In contrast, the hydrological signature derived from the time series could provide a better picture of a process characteristic resulting from the drivers and hence a step closer to understanding the change of process and is less prone to artifacts caused by single point analysis.
This study focuses on the application of recurrence based data analysis techniques (recurrence plot) for understanding and quantifying spatio-temporal changes in flood hazards in Germany through its hydrological signature. The recurrence plot is known as an effective tool to visualize the dynamics of phase space trajectories i.e. constructed from a time series by using an embedding dimension and a time delay, and it is known to be effective in analyzing non-stationary and non-linear time series. The emphasis will be on the identification of characteristic recurrence properties that could associate typical dynamic behavior to certain flood situations.
A. M. T. Ramos, A. Builes-Jaramillo, G. Poveda, B. Goswami, E. E. N. Macau, J. Kurths, N. Marwan:
Causality detection based on recurrence plot,
7th International Symposium on Recurrence Plots,
São Paulo (Brazil),
August 23-25, 2017,
Talk.
» Abstract
We will present the Recurrence Measure of Conditional Dependence (RMCD), a recent data-driven causality inference method using the framework of recurrence plots. The RMCD incorporates the recurrence behavior into the transfer entropy theory. We will discuss how this methodology can reveal the lagged coupling of some paradigmatic models and how it reveals causal relations of climate systems. For instance, RMCD detects the influence of the Pacific Ocean temperatures on the South West Amazon rainfall during the 2010 droughts, as well as its influence absence during 2005.
D. Eroglu, N. Marwan:
Multiplex Recurrence Networks,
EGU General Assembly,
Vienna (Austria),
April 23-28, 2017,
» Poster (PDF, 5.64M)
.
» Abstract
The complex nature of a variety of phenomena in physical, biological, or earth sciences is driven by a large number of degrees of freedom which are strongly interconnected. Although the evolution of such systems is described by multivariate time series (MTS), so far research mostly focuses on analyzing these components one by one.
Recurrence based analyses are powerful methods to understand the underlying dynamics of a dynamical system and have been used for many successful applications including examples from earth science, economics, or chemical reactions. The backbone of these techniques is creating the phase space of the system. However, increasing the dimension of a system requires increasing the length of the time series in order get significant and reliable results. This requirement is one of the challenges in many disciplines, in particular in palaeoclimate, thus, it is not easy to create a phase space from measured MTS due to the limited number of available obervations (samples). To overcome this problem, we suggest to create recurrence networks from each component of the system and combine them into a multiplex network structure, the multiplex recurrence network (MRN). We test the MRN by using prototypical mathematical models and demonstrate its use by studying high-dimensional palaeoclimate dynamics derived from pollen data from the Bear Lake (Utah, US). By using the MRN, we can distinguish typical climate transition events, e.g., such between Marine Isotope Stages.
F. Brenner, N. Marwan, P. Hoffmann:
Modelling fast spreading patterns of airborne infectious diseases using complex networks,
EGU General Assembly,
Vienna (Austria),
April 23-28, 2017,
Talk.
» Abstract
The pandemics of SARS (2002/2003) and H1N1 (2009) have impressively shown the potential of epidemic outbreaks of infectious diseases in a world that is strongly connected. Global air travelling established an easy and fast opportunity for pathogens to migrate globally in only a few days. This made epidemiological prediction harder. By understanding this complex development and its link to climate change we can suggest actions to control a part of global human health affairs.
In this study we combine the following data components to simulate the outbreak of an airborne infectious disease that is directly transmitted from human to human:
- GlobalAirTrafficNetwork(fromopenflights.org) with information on airports, airportlocation, directflight connection, airplane type
- Global population dataset (from SEDAC, NASA)
- Susceptible-Infected-Recovered (SIR) compartmental model to simulate disease spreading in the vicinity of airports. A modified Susceptible-Exposed-Infected-Recovered (SEIR) model to analyze the impact of the incubation period.
- WATCH-Forcing-Data-ERA-Interim(WFDEI) climatedata: temperature, specific humidity, surface air pressure, and water vapor pressure
These elements are implemented into a complex network. Nodes inside the network represent airports. Each single node is equipped with its own SIR/SEIR compartmental model with node specific attributes. Edges between those nodes represent direct flight connections that allow infected individuals to move between linked nodes. Therefore the interaction of the set of unique SIR models creates the model dynamics we will analyze.
To better figure out the influence on climate change on disease spreading patterns, we focus on Influenza-like-Illnesses (ILI). The transmission rate of ILI has a dependency on climate parameters like humidity and temperature. Even small changes of environmental variables can trigger significant differences in the global outbreak behavior. Apart from the direct effect of climate change on the transmission of airborne diseases, there are indirect ramifications that alter spreading patterns. An example is seasonal human mobility behavior which will change with varied climate conditions. The direct and indirect effects of climate change on disease spreading patterns will be discussed in this study.
A. Agarwal, N. Marwan, M. Rathinasamy, U. Oeztuerk, B. Merz, J. Kurths:
Multiscale complex network analysis: An approach to study spatiotemporal rainfall pattern in south Germany,
EGU General Assembly,
Vienna (Austria),
April 23-28, 2017,
Poster.
» Abstract
Understanding of the climate sytems has been of tremendous importance to different branches such as agriculture, flood, drought and water resources management etc. In this regard, complex networks analysis and time series analysis attracted considerable attention, owing to their potential role in understanding the climate system through characteristic properties. One of the basic requirements in studying climate network dynamics is to identify connections in space or time or space-time, depending upon the purpose. Although a wide variety of approaches have been developed and applied to identify and analyse spatio-temporal relationships by climate networks, there is still further need for improvements in particular when considering precipitation time series or interactions on different scales. In this regard, recent developments in the area of network theory, especially complex networks, offer new avenues, both for their generality about systems and for their holistic perspective about spatio-temporal relationships.
The present study has made an attempt to apply the ideas developed in the field of complex networks to examine connections in regional climate networks with particular focus on multiscale spatiotemporal connections. This paper proposes a novel multiscale understanding of regional climate networks using wavelets. The proposed approach is applied to daily precipitation records observed at 543 selected stations from south Germany for a period of 110 years (1901-2010). Further, multiscale community mining is performed on the same study region to shed more light on the underlying processes at different time scales.
Various network measure and tools so far employed provide micro-level (individual station) and macro-level (community structure) information of the network. It is interesting to investigate how the result of this study can be useful for future climate predictions and for evaluating climate models on their implementation regarding heavy precipitation.
U. Ozturk, N. Marwan, J. Kurths:
Identifying typhoon tracks based on event synchronization derived spatially embedded climate networks,
EGU General Assembly,
Vienna (Austria),
April 23-28, 2017,
Poster.
» Abstract
Complex networks are commonly used for investigating spatiotemporal dynamics of complex systems, e.g. extreme rainfall. Especially directed networks are very effective tools in identifying climatic patterns on spatially embedded networks. They can capture the network flux, so as the principal dynamics of spreading significant phenomena. Network measures, such as network divergence, bare the source-receptor relation of the directed networks. However, it is still a challenge how to catch fast evolving atmospheric events, i.e. typhoons.
In this study, we propose a new technique, namely Radial Ranks, to detect the general pattern of typhoons forward direction based on the strength parameter of the event synchronization over Japan. We suggest to subset a circular zone of high correlation around the selected grid based on the strength parameter. Radial sums of the strength parameter along vectors within this zone, radial ranks are measured for potential directions, which allows us to trace the network flux over long distances. We employed also the delay parameter of event synchronization to identify and separate the frontal storms’ and typhoons’ individual behaviors.
U. Ozturk, N. Marwan, O. Korup, J. Jensen:
Completing the record of 20th century sea level rise in the Eastern Mediterranean,
EGU General Assembly,
Vienna (Austria),
April 23-28, 2017,
Poster.
» Abstract
Quantitative studies of sea-level rise in the Mediterranean are becoming more and more accurate thanks to detailed satellite monitoring campaigns. However, these studies cover several years to a couple of decades at best, while longer-term sea-level records for the area are rare. Long-term sea-level measurements are essential in order to derive accurate trends free of conspicuous oscillations in shorter records. We use an approach from data archaeology to meet this shortcoming, and to offer a more complete record of sea-level rise cross-checked among several tide gauges. Specifically, we investigate monthly mean sea-level data of the Antalya-I (1935-1977) tide gauge provided by the Turkish National Mapping Agency. We checked how accurately and reliably these monthly records were digitized, quality-controlled, and tied to a common datum. We then merged these data with the more recent records of the nearby Antalya-II (1985-2010) tide gauge, obtaining a composite time series of monthly and annual mean sea levels spanning approximately 75 years. We thus offer the hitherto longest record in the Eastern Mediterranean Basin as an essential tool for studying the region’s sea-level trends. We estimate a relative mean sea-level rise of 2.46 ± 1.65 mm/yr between 1935 and 2010, with a sub-decadal variability (σresiduals = 49.47 mm) that is higher than at nearby tide gauges (e.g. Thessaloniki, Greece, σresiduals = 28.71 mm). Our study highlights the value of data archaeology for recovering and integrating early tide-gauge data for long-term sea-level and climate studies.
N. Marwan:
Recurrence Plots for the Analysis of Complex Systems,
Chapman Chair Seminar "Complex Systems Science Meets Arctic Science", University of Fairbanks,
Fairbanks (USA),
March 2, 2017,
Lecture.
N. Marwan:
Caves as scientific archives,
Chapman Chair Seminar "Complex Systems Science Meets Arctic Science", University of Fairbanks,
Fairbanks (USA),
March 2, 2017,
Lecture.
A. M. D. T. Ramos, L. A. Builes-Jaramillo, G. Poveda, B. Goswami, E. E. N. Macau, J. Kurths, N. Marwan:
Non-linear interactions between Niño region 3 and the Southern Amazon,
AGU Fall Meeting,
San Francisco (USA),
December 12-16, 2016,
» Talk (PDF, 2.17M)
.
» Abstract
Identifying causal relations from the observational dataset has posed great challenges in data-driven inference study. However, complex system framework offers promising approaches to tackle such problems. Here we propose a new data-driven causality inference method using the framework of recurrence plots. We present the Recurrence Measure of Conditional Dependence (RMCD) and its applications. The RMCD incorporates the recurrence behavior into the transfer entropy theory. Therefore, it quantifies the causal dependence between two processes based on joint recurrence patterns between the past of the potential driver and present on the potential driven, except for any contribution that has already been in the past of the driven. We apply this methodology to some paradigmatic models and to investigate the possible influence of the Pacific Ocean temperatures on the South West Amazon for the 2010 and 2005 droughts. The results reveal that for the 2005 drought there is not a significant signal of dependence from the Pacific Ocean and that for 2010 there is a signal of dependence of around 200 days. These outcomes are confirmed by the traditional climatological analysis of these episodes available in the literature and show the accuracy of RMCD inferring causal relations in climate systems.
A. Agarwal, N. Marwan, R. Maheswaran, U. Ozturk, B. Merz, J. Kurths:
Multiscale event synchronization analysis for unraveling climate processes: A wavelet-based approach,
AGU Fall Meeting,
San Francisco (USA),
December 12-16, 2016,
» Poster (PDF, 2.43M)
.
» Abstract
The temporal dynamics of natural processes (climatic/hydrological) are spread across different time scales and, as such, the study of these processes only at a given scale would not reveal all the underlying governing processes spanning over different scales. As needs be, it is vital to investigate such processes at various time-scales. Wavelets have been used extensively to comprehend the multiscale process and have been appeared to be exceptionally reliable and useful in understanding dynamics of the process across various time scales as these evolve in time. Recently, event synchronization has picked up an interest in capturing the nonlinear interactions between different climatic signals based on synchronization and time delays of events (characterized as local maxima) in the timeseries. At present, the event synchronization analyses the relationship at a single scale. However, visualization of the time evolution of delay and synchronization at multiscale would be of great interest to comprehend the dynamics of the signal.
In this paper, a wavelet based multi-scale event synchronization (MSES) methodology has been proposed and tested. The main advantage of this method is that it provide a quantitative measure of process across different scales. The applicability of the proposed method has explored using various case studies –both real as well as synthetic. The study also shows that proposed methodology works well also for model evaluation.
D. Wendi, N. Marwan, B. Merz:
Identifying changes of complex flood dynamics with recurrence analysis,
AGU Fall Meeting,
San Francisco (USA),
December 12-16, 2016,
Talk invited.
» Abstract
Temporal changes in flood hazard system are known to be difficult to detect and attribute due to multiple drivers that include complex processes that are non-stationary and highly variable. These drivers, such as human-induced climate change, natural climate variability, implementation of flood defense, river training, or land use change, could impact variably on space-time scales and influence or mask each other. Flood time series may show complex behavior that vary at a range of time scales and may cluster in time. Moreover hydrological time series (i.e. discharge) are often subject to measurement errors, such as rating curve error especially in the case of extremes where observation are actually derived through extrapolation.
This study focuses on the application of recurrence based data analysis techniques (recurrence plot) for understanding and quantifying spatio-temporal changes in flood hazard in Germany. The recurrence plot is known as an effective tool to visualize the dynamics of phase space trajectories i.e. constructed from a time series by using an embedding dimension and a time delay, and it is known to be effective in analyzing non-stationary and non-linear time series. Sensitivity of the common measurement errors and noise on recurrence analysis will also be analyzed and evaluated against conventional methods. The emphasis will be on the identification of characteristic recurrence properties that could associate typical dynamic to certain flood events.
F. Brenner, P. Hoffmann, N. Marwan:
Combining a complex network approach and a SEIR compartmental model to link fast spreading of infectious diseases with climate change,
AGU Fall Meeting,
San Francisco (USA),
December 12-16, 2016,
Talk.
» Abstract
Infectious diseases are a major threat to human health. The spreading of airborne diseases has become fast and hard to predict. Global air travelling created a network which allows a pathogen to migrate worldwide in only a few days. Pandemics of SARS (2002/03) and H1N1 (2009) have impressively shown the epidemiological danger in a strongly connected world.
In this study we simulate the outbreak of an airborne infectious disease that is directly transmitted from human to human. We use a regular Susceptible-Infected-Recovered (SIR) model and a modified Susceptible-Exposed-Infected-Recovered (SEIR) compartmental approach with the basis of a complex network built by global air traffic data (from openflights.org). Local Disease propagation is modeled with a global population dataset (from SEDAC and MaxMind) and parameterizations of human behavior regarding mobility, contacts and awareness. As a final component we combine the worldwide outbreak simulation with daily averaged climate data from WATCH-Forcing-Data-ERA-Interim (WFDEI) and Coupled Model Intercomparison Project Phase 5 (CMIP5).
Here we focus on Influenza-like illnesses (ILI), whose transmission rate has a dependency on relative humidity and temperature. Even small changes in relative humidity are sufficient to trigger significant differences in the global outbreak behavior. Apart from the direct effect of climate change on the transmission of airborne diseases, there are indirect ramifications that alter spreading patterns. For example seasonal changing human mobility is influenced by climate settings.
N. Marwan:
Hot Topics in Recurrence Plot Analysis,
Humboldt-Kolleg,
Yaoundé (Cameroon),
November 22-24, 2016,
Lecture.
N. Marwan:
Recurrence Analysis,
Graduate School NatRiskChange, University of Potsdam,
Potsdam (Germany),
November 8, 2016,
Lecture.
N. Marwan:
Analysing spatially extended high-dimensional dynamics by Recurrence plots,
Conference Perspectives in Nonlinear Dynamics (PNLD2016),
Berlin (Germany),
July 24-26, 2016,
Talk.
» Abstract
Recurrence plot based measures of complexity are capable tools for characterizing complex dynamics. We show the potential of selected recurrence plot measures for the investigation of even high-dimensional dynamics. We apply this method on spatially extended chaos, such as derived from the Lorenz96 model and show that the recurrence plot based measures can qualitatively characterize typical dynamical properties such as chaotic or periodic dynamics. Moreover, we demonstrate its power by analyzing satellite image time series of vegetation cover with contrasting dynamics as a spatially extended and potentially high-dimensional example from the real world.
N. Marwan:
Modern approaches for nonlinear time series analysis,
Course Global Change Management Studies, University of Applied Science Eberswalde,
Eberswalde (Germany),
June 7, 2016,
Lecture.
N. Marwan, S. Foerster, J. Kurths:
Analysing spatially extended high-dimensional dynamics by recurrence plots,
EGU General Assembly,
Vienna (Austria),
April 17-22, 2016,
» Poster (PDF, 2.92M)
.
» Abstract
Recurrence plot based measures of complexity are capable tools for characterizing complex dynamics. We show the potential of selected recurrence plot measures for the investigation of even high-dimensional dynamics. We apply this method on spatially extended chaos, such as derived from the Lorenz96 model and show that the recurrence plot based measures can qualitatively characterize typical dynamical properties such as chaotic or periodic dynamics. Moreover, we demonstrate its power by analyzing satellite image time series of vegetation cover with contrasting dynamics as a spatially extended and potentially high-dimensional example from the real world.
N. Marwan, N. Malik, N. Boers, A. Rheinwalt, V. Stolbova, B. Bookhagen, J. Kurths:
Potentials of complex network analysis of regional rainfall co-variability,
EGU General Assembly,
Vienna (Austria),
April 17-22, 2016,
Poster.
» Abstract
Regional variability of extreme precipitation is characterised by distinctive patterns of statistical interrelations. The propagation, structure, and complexity of such patterns can be studied by means of complex network theory. We demonstrate the potential of this analytical framework for identification of teleconnections, propagation patterns, or finding prediction schemes for extreme rainfall.
T. Nocke, S. Buschmann, J. Donges, N. Marwan:
Visualization techniques and tools for large geo-physical networks,
EGU General Assembly,
Vienna (Austria),
April 17-22, 2016,
Talk.
» Abstract
Network analysis is an important approach in studying complex phenomena within geophysical observation and simulation data. This field produces increasing numbers of large geo-referenced networks to be analyzed. Particular focus lies on the network analysis of the complex statistical interrelationship structure within climatological fields. The typical procedure for such network analyzes is the extraction of network measures in combination with static standard visualization methods.
To analyze the visualization challenges within this field, we performed a questionnaire with climate and complex system scientists, and identified a strong requirement for solutions visualizing large and very large geo-referenced networks by providing alternative mappings for static plots and allowing for interactive visualization for networks with 100.000 or even millions of edges. In addition, the questionnaire revealed, that existing interactive visualization methods and tools for geo-referenced network exploration are often either not known to the analyst or their potential is not fully exploited.
Within this presentation, we illustrate how interactive visual analytics methods in combination with geo-visualisation can be tailored for visual large climate network investigation (see as well Nocke et al. 2015). Therefore, we present a survey of requirements of network analysts and the related challenges and, as an overview for the interested practitioner, we review the state-of-the-art in climate network visualization techniques and tools, underpinned with concrete examples from climate network research and innovative solutions (e.g. alternative projections, 3D layered networks) implemented within the network visualization system GTX.
J. Runge, V. Petoukhov, J. Donges, J. Hlinka, N. Jajcay, M. Vejmelka, D. Hartman, N. Marwan, M. Palus, J. Kurths:
Identifying causal gateways and mediators in complex spatio-temporal systems,
EGU General Assembly,
Vienna (Austria),
April 17-22, 2016,
Talk.
» Abstract
Identifying regions important for spreading and mediating perturbations is crucial to assess the susceptibilities of spatio-temporal complex systems such as the Earth’s climate to volcanic eruptions, extreme events or geoengineering. Here a data-driven approach is introduced based on a dimension reduction, causal reconstruction, and novel network measures based on causal effect theory that go beyond standard complex network tools by distinguishing direct from indirect pathways. Applied to a data set of atmospheric dynamics, the method identifies several strongly uplifting regions acting as major gateways of perturbations spreading in the atmosphere. Additionally, the method provides a stricter statistical approach to pathways of atmospheric teleconnections, yielding insights into the Pacific–Indian Ocean interaction relevant for monsoonal dynamics. The novel causal interaction perspective provides a complementary approach to simulations or experiments for understanding the functioning of complex spatio-temporal systems with potential applications in increasing their resilience to shocks or extreme events.
T. Nocke, S. Buschmann, J. Donges, N. Marwan:
New solutions for climate network visualization,
EGU General Assembly,
Vienna (Austria),
April 17-22, 2016,
Talk.
» Abstract
An increasing amount of climate and climate impact research methods deals with geo-referenced networks, including energy, trade, supply-chain, disease dissemination and climatic tele-connection networks. At the same time, the size and complexity of these networks increases, resulting in networks of more than hundred thousand or even millions of edges, which are often temporally evolving, have additional data at nodes and edges, and can consist of multiple layers even in real 3D. This gives challenges to both the static representation and the interactive exploration of these networks, first of all avoiding edge clutter (“edge spagetti”) and allowing interactivity even for unfiltered networks.
Within this presentation, we illustrate potential solutions to these challenges. Therefore, we give a glimpse on a questionnaire performed with climate and complex system scientists with respect to their network visualization requirements, and on a review of available state-of-the-art visualization techniques and tools for this purpose (see as well Nocke et al., 2015). In the main part, we present alternative visualization solutions for several use cases (global, regional, and multi-layered climate networks) including alternative geographic projections, edge bundling, and 3-D network support (based on CGV and GTX tools), and implementation details to reach interactive frame rates.
S. F. M. Breitenbach, B. Plessen, S. Wenz, J. Leonhardt, R. Tjallingii, D. Scholz, K. P. Jochum, N. Marwan:
A multi-proxy reconstruction of Holocene climate change from Blessberg Cave,
EGU General Assembly,
Vienna (Austria),
April 17-22, 2016,
» Poster (PDF, 2.85M)
.
» Abstract
Although Holocene climate dynamics were relatively stable compared to glacial conditions, climatic changes had significant impact on ecosystems and human society on various timescales (Mayewski et al. 2004, Donges et al. 2015, Tan et al. 2015). Precious few high-resolution records on Holocene temperature and precipitation conditions in Central Europe are available (e.g., von Grafenstein et al. 1999, Fohlmeister et al. 2012).
Here we present a speleothem-based reconstruction of past climate dynamics from Blessberg Cave, Thuringia, central Germany. Three calcitic stalagmites were recovered when the cave was discovered during tunneling operations in 2008. Samples BB-1, -2 and -3 were precisely dated by the 230Th/U-method, with errors between 10 and 160 years (2σ). The combined record covers large parts of the Holocene (10 – 0.4 ka BP). δ13C and δ18O were analysed at 100 μm resolution. To gain additional insights in infiltration conditions, Sr/Ca and S/Ca were measured on BB-1 and BB-3 using an Röntgenanalytik Eagle XXL μXRF scanner.
Differences to other central European records (e.g., von Grafenstein et al. 1999, Fohlmeister et al. 2012) suggest complex interaction between multiple factors influencing speleothem δ18O in Blessberg Cave. Furthermore, no clear influence of the North Atlantic Oscillation on our proxies is found. However, a link across the N Atlantic realm is indicated by a centennial-scale correlation between Blessberg δ18O values and minerogenic input into lake SS1220 in Greenland over the last 5 ka (Olsen et al. 2012). In addition, recurrence analysis indicates an imprint of Atlantic Bond events on Blessberg δ18O values (Marwan et al. 2014), corroborating the suggested link with high northern latitudes. Larger runoff into the Greenland lake seems to be associated with lower δ18O, higher δ13C and S/Ca ratios, as well as lower Sr/Ca ratios in Blessberg Cave speleothems. This might be linked to lower local temperature and/or changes in precipitation seasonality. Opposing millennial scale trends with lowering S/Ca ratios and δ13C values but increasing Sr/Ca ratios calls for more than one controlling factor. Most likely, δ13C decreased through the Holocene due to afforestation, which in turn might have increased sulphate retention in the thickening soil cover (Frisia et al. 2005) and limited sulphur flux into the cave. Alternatively, marine sulfur flux could have diminished with winter wind intensities. However, additional data is required to clarify this hypothesis. A positive Sr/Ca trend through the Holocene might result from increasing prior calcite precipitation induced by a negative moisture balance in summer.
D. Wendi, B. Merz, N. Marwan:
Novel flood detection and analysis method using recurrence property,
EGU General Assembly,
Vienna (Austria),
April 17-22, 2016,
Poster.
» Abstract
Temporal changes in flood hazard are known to be difficult to detect and attribute due to multiple drivers that include processes that are non-stationary and highly variable. These drivers, such as human-induced climate change, natural climate variability, implementation of flood defence, river training, or land use change, could impact variably on space-time scales and influence or mask each other. Flood time series may show complex behavior that vary at a range of time scales and may cluster in time.
This study focuses on the application of recurrence based data analysis techniques (recurrence plot) for understanding and quantifying spatio-temporal changes in flood hazard in Germany. The recurrence plot is known as an effective tool to visualize the dynamics of phase space trajectories i.e. constructed from a time series by using an embedding dimension and a time delay, and it is known to be effective in analyzing non-stationary and non-linear time series. The emphasis will be on the identification of characteristic recurrence properties that could associate typical dynamic behavior to certain flood situations.
J. Donges, J. Heitzig, B. Beronov, M. Wiedermann, J. Runge, Q. Y. Feng, L. Tupikina, V. Stolbova, R. Donner, N. Marwan, H. Dijkstra, J. Kurths:
Unified functional network and nonlinear time series analysis for complex systems science: The pyunicorn package,
EGU General Assembly,
Vienna (Austria),
April 17-22, 2016,
Poster.
» Abstract
We introduce the pyunicorn (Pythonic unified complex network and recurrence analysis toolbox) open source software package for applying and combining modern methods of data analysis and modeling from complex network theory and nonlinear time series analysis. pyunicorn is a fully object-oriented and easily parallelizable package written in the language Python. It allows for the construction of functional networks such as climate networks in climatology or functional brain networks in neuroscience representing the structure of statistical interrelationships in large data sets of time series and, subsequently, investigating this structure using advanced methods of complex network theory such as measures and models for spatial networks, networks of interacting networks, node-weighted statistics, or network surrogates. Additionally, pyunicorn provides insights into the nonlinear dynamics of complex systems as recorded in uni- and multivariate time series from a non-traditional perspective by means of recurrence quantification analysis, recurrence networks, visibility graphs, and construction of surrogate time series. The range of possible applications of the library is outlined, drawing on several examples mainly from the field of climatology. pyunicorn is available online at https://github.com/pik-copan/pyunicorn.
D. Eroglu, N. Marwan, I. Ozken, T. Stemler, J. Kurths:
How to overcome data based difficulties in geoscience?,
EGU General Assembly,
Vienna (Austria),
April 17-22, 2016,
Poster.
» Abstract
There are many challenges in the field of time series analysis such as cumulative trends or sampling irregularities. Geophysical time series, particularly paleo-climate ones, have such problems almost in all proxies. The novel TrAnsformation-Cost Time-Series (TACTS) method is a suitable approach to overcome these challenges of cumulative trends and irregular sampling without degenerating the quality of the data set by, e.g., interpolation.
The standard method to regularize time sampling of time series is interpolation, but it collapses the quality of the proxies. Moreover, there are many different approaches to de-trend time series such as Gaussian high-pass filter, the de-trended fluctuation analysis. At the same time, the TACTS is able to de-trend and regularize the time series at the same time with keeping the quality of time series rather high. After applying the TACTS method the resulting cost time series shows regular sampling and can be further analyzed using standard methods.
The TACTS method has been developed and tested by using prototypical mathematical models. We have demonstrated its use by studying paleoclimate dynamics derived from speleothem data from the Secret Cave in Borneo, the KNI-51 cave in North Australia, and Dongge Cave, East China. By using the TACTS, we could distinguish all extreme transition events and found interesting alternating monsoon dynamics between North Australia and East China.
M. Riedl, N. Marwan, J. Kurths:
Extended quantification of the generalized recurrence plot,
EGU General Assembly,
Vienna (Austria),
April 17-22, 2016,
Poster.
» Abstract
The generalized recurrence plot is a modern tool for quantification of complex spatial patterns. Its application spans the analysis of trabecular bone structures, Turing structures, turbulent spatial plankton patterns, and fractals. But, it is also successfully applied to the description of spatio-temporal dynamics and the detection of regime shifts, such as in the complex Ginzburg-Landau-equation. The recurrence plot based determinism is a central measure in this framework quantifying the level of regularities in temporal and spatial structures. We extend this measure for the generalized recurrence plot considering additional operations of symmetry than the simple translation. It is tested not only on two-dimensional regular patterns and noise but also on complex spatial patterns reconstructing the parameter space of the complex Ginzburg-Landau-equation. The extended version of the determinism resulted in values which are consistent to the original recurrence plot approach. Furthermore, the proposed method allows a split of the determinism into parts which based on laminar and non-laminar regions of the two-dimensional pattern of the complex Ginzburg-Landau-equation. A comparison of these parts with a standard method of image classification, the co-occurrence matrix approach, shows differences especially in the description of patterns associated with turbulence. In that case, it seems that the extended version of the determinism allows a distinction of phase turbulence and defect turbulence by means of their spatial patterns. This ability of the proposed method promise new insights in other systems with turbulent dynamics coming from climatology, biology, ecology, and social sciences, for example.
O. Kwiecien, N. Marwan:
What drives the δ18O signal in European carbonates during LGM and Termination II – multi-archive approach revised,
EGU General Assembly,
Vienna (Austria),
April 17-22, 2016,
Poster.
» Abstract
Due to their natural occurrence in a wide spectrum of environmental settings (marine, lacustrine, terrestrial) carbonates constitute one of the most widely used archives of (paleo)environmental information Although the thermodynamic laws controlling isotopic fractionation are universal, each of the different environments leaves its characteristic imprint on the carbonate chemistry. Consequently, individual carbonate-based records are a combination of regional-scale and archive- or site-specific factors. Environmental heterogeneity, even within a confined geographical area, remains a challenge in teasing out the regional or local nature of the recorded climatic signal. Similarly, interpretation of new data often pays tribute to outdated assumptions.
Here we compare and contrast published oxygen isotope records from different archives (lake sediments and stalagmites) from central Europe and the Black Sea region covering the time span between 26 and 8 ka BP. This exercise aims at (1) identification of the regional-scale mechanisms responsible for common trends and (2) better understanding of archive- and site-specific factors accounting for differences in the records. Additionally, this approach allows for testing ‘closed vs outflowing’ scenarios for the glacial Black Sea basin.
It appears that changes in isotopic composition of atmospheric precipitation are a common denominator for the analyzed records. Site-specific factors include moisture source (stalagmites) and volume of the basin (lakes). Both, comparison of available geochemical records and data-based theoretical calculations suggest that since the LGM and until reconnection with Mediterranean at ca. 8 ka BP the Black Sea was an open system.
D. Eroglu, I. Ozken, F. McRobie, T. Stemler, N. Marwan, K.-H. Wyrwoll, J. Kurths:
A solar variability driven monsoon see-saw: switching relationships of the Holocene East Asian-Australian summer monsoons,
EGU General Assembly,
Vienna (Austria),
April 17-22, 2016,
» Poster (PDF, 842.68K)
.
» Abstract
The East Asian-Indonesian-Australian monsoon is the predominant low latitude monsoon system, providing a major global scale heat source. Here we apply newly developed non-linear time series techniques on speleothem climate proxies, from eastern China and northwestern Australia and establish relationships between the two summer monsoon regimes over the last ∼9000 years. We identify significant variations in monsoonal activity, both dry and wet phases, at millennial to multi-centennial time scales and demonstrate for the first time the existence of a see-saw antiphase relationship between the two regional monsoon systems. Our analysis attributes this inter-hemispheric linkage to the solar variability that is effecting both monsoon systems.
N. Marwan, Y. Zou, M. Riedl, N. Wessel, J. Kurths:
Bestimmung der Kopplungsrichtungen im kardiorespiratorischen System anhand wiederkehrbasierter Methoden,
Workshop Biosignalverarbeitung 2016,
Berlin (Germany),
April 7, 2016,
Talk.
» Abstract
The asymmetry of coupling between complex systems can be studied by conditional probabilities of recurrence, which can be estimated by joint recurrence plots. This approach is applied for the first time on experimental data: time series of the human cardiorespiratory system in order to investigate the couplings between heart rate, mean arterial blood pressure and respiration. We find that the respiratory system couples towards the heart rate, and the heart rate towards the mean arterial blood pressure. However, our analysis could not detect a clear coupling direction between the mean arterial blood pressure and respiration.
N. Marwan:
Recurrence Plots for the Analysis of Complex Systems,
Seminar Department of Knowledge Engineering, University of Maastricht,
Maastricht (The Netherlands),
February 17, 2016,
Talk.
N. Boers, A. Rheinwalt, B. Bokkhagen, N. Marwan, J. Kurths, H. M. J. Barbosa, J. Marengo:
Complex network analysis of extreme rainfall in South America: Climatic analysis – model evaluation – prediction,
International Workshop "Analysis of dynamic networks and data driven modelling of the climate",
Potsdam (Germany),
October 12–14, 2015,
Talk.
» Abstract
We construct networks from synchronization of extreme rainfall events over South America for the monsoon season from December to February, using 6 different datasets. This methodology is designed to complement PCA-based techniques in the case of extreme events. A climatological interpretation of various network measures (1) reveals the most important features of the South American Monsoon System (SAMS), (II) allows to compare the representation of the spatiotemporal co-variability of extreme rainfall in different datasets and models, and (III) leads to a statistical forecasting framework of extreme events.
B. Goswami, A. Rheinwalt, N. Boers, N. Marwan, J. Heitzig, S. F. M. Breitenbach, J. Kurths:
Using complex networks to detect abrupt climate change in paleoclimate datasets,
International Workshop "Analysis of dynamic networks and data driven modelling of the climate",
Potsdam (Germany),
October 12–14, 2015,
Talk.
» Abstract
A fundamental obstacle in the analysis of paleoclimate datasets is the uncertainty involved in determining the liming of past climatic events. Principled approaches to such data necessarily represents them as a probability distribution at eacli time point of the past rather than a point estimate. Here, we provide an approach that helps infer dynamical shifts from paleoclimate proxy records, considering them as a sequence of possibly correlated time-ordered marginal probability distributions instead of a point-estimate time series. Using bounds on the difference distriburions we define a probability of recurrence matrix with which we estimate the recurrence network for the data. We propose that the identification of communities in the recurrence network corresponds to the detection of sudden dynamical shifts. We demonstrate our approach with a synthetic paleomonsoon dataset and use it thereafter to infer periods of sudden climatic change in Asian monsoon records that span the past 9000 years.
M. Sips, T. Rawald, C. Witt, N. Marwan:
Towards multi-scale RQA using Visual Analytics,
6th International Symposium on Recurrence Plots,
Grenoble (France),
June 17-19, 2015,
Talk.
» Abstract
Time series often capture a variety of sub processes on different temporal scales. Applying the recurrence quantification analysis (RQA) to time series describes the behavior of all processes simultaneously. In many application scenarios, users want to restrict the RQA only on relevant processes. Therefore, there is increasing interest to extract relevant processes from time series.
Previous work has shown that multi-scale methods, such as wavelet transformation, are capable of isolating relevant processes. Our research goal is to extent the traditional RQA into a multi-scale RQA method. Our current research efforts show that multi-scale RQA involves two major challenges. First, the decomposition of a time series produces a huge set of time series; each time series is associated with a particular temporal scale. Second, users need to locate the temporal scales at which relevant processes operate. To address these two challenges, we currently work on a Visual Analytics approach that combines multi-scale decomposition, the fast computation of RQA for all temporal scales, and the clustering of these RQA results.
In our presentation, we will discuss the three main components of our Visual Analytics approach. First, our approach supports a broad range of multi-scale methods. Second, our fast RQA computation is based on subdividing the RQA computation and distributing the computational work across multiple graphics processing units. Third, the clustering of RQA measures allow users to focus on characteristic RQA results. We will present our latest research efforts and discuss the challenges involved to extent traditional RQA into a multi-scale method.
D. Eroglu, T. K. D. Peron, N. Marwan, F. A. Rodrigues, L. da F. Costa, M. Sebek, I. Z. Kiss, J. Kurths:
Entropy of weighted recurrence plots,
6th International Symposium on Recurrence Plots,
Grenoble (France),
June 17-19, 2015,
Talk.
» Abstract
Phase-space binning and using information in the bins is the standard method to establish complexity of a dynamical system. Methods such as the Shannon entropy, the Hausdorff dimension, the Kolmogorov complexity etc., can be used to quantify transitions between different dynamical regimes. In order to create a recurrence plot, we divide the phase space into equi-distant but overlapping bins. The information in the bins shows us similar features as in the mentioned measures before. Here we suggest a method based on weighted recurrence plots and show that the associated Shannon entropy is positively correlated with the largest Lyapunov exponent. We demonstrate the potential on a prototypical example as well as on experimental data of a chemical experiment.
D. Schultz, S. Spiegel, N. Marwan, S. Albayrak:
Approximation of diagonal line based measures in recurrence quantification analysis,
6th International Symposium on Recurrence Plots,
Grenoble (France),
June 17-19, 2015,
Talk.
» Abstract
Given a trajectory of length N, recurrence quantification analysis (RQA) traditionally operates on the recurrence plot, whose calculation requires quadratic time and space (O(N2)), leading to expensive computations and high memory usage for large N. However, if the similarity threshold ε is zero, we show that the recurrence rate (RR), the determinism (DET) and other diagonal line based RQA-measures can be obtained algorithmically taking O(N log(N)) time and O(N) space. Furthermore, for the case of ε > 0 we propose approximations to the RQA-measures that are computable with same complexity. Simulations with autoregressive systems, the logistic map and a Lorenz attractor suggest that the approximation error is small if the dimension of the trajectory and the minimum diagonal line length are small. When applying the approximate determinism to the problem of detecting dynamical transitions we observe that it performs as well as the exact determinism measure.
N. Marwan, S. Foerster, J. Kurths:
Recurrence plot analysis of spatially extended high-dimensional dynamics,
6th International Symposium on Recurrence Plots,
Grenoble (France),
June 17-19, 2015,
» Talk (PDF, 5.77M)
.
» Abstract
Recurrence plot based measures of complexity are capable tools for characterizing complex dynamics. We show the potential of selected recurrence plot measures for the investigation of spatially extended high-dimensional dynamics by applying them to data from the Lorenz96 model. The recurrence plot based measures are able to qualitatively characterize typical dynamical properties such as chaotic or periodic dynamics. Moreover, we demonstrate its power by analyzing satellite image time series of vegetation cover with contrasting dynamics as a spatially extended and potentially high-dimensional example from the real world.
M. Riedl, J. Kurths, N. Marwan:
Spatial-temporal recurrence analysis based on a global measure of spatial similarity,
6th International Symposium on Recurrence Plots,
Grenoble (France),
June 17-19, 2015,
Talk.
» Abstract
The analysis of spatial-temporal data is still a methodological question of current research, e.g. in earth science. In this work an extension of the recurrence plot is suggested in order to track regime shifts in spatial patterns. For this purpose, the global state, i.e. all data in one time step, is compared to each other by means of a similarity measure that is known from tracking algorithms. This routine consists of a digitalization, blurring and bin wise comparison of the single global states. In its roughest version, the algorithm is equivalent to the kappa-statistic which is widely used to assess similarities of spatial pattern in ecosystems. The use of this similarity measure in the framework of recurrence plots enables a visualization and quantification of the temporal evolution of the spatial extended system. In the end, the extended recurrence plot and its quantification allows a detection of changes in the spatial-temporal dynamic of the observed system.
O. Afsar, N. Marwan, J. Kurths:
Scaling relations from recurrence quantification analysis for the Logistic map at the edge of chaos: Connection with universal Huberman-Rudnick scaling law,
6th International Symposium on Recurrence Plots,
Grenoble (France),
June 17-19, 2015,
Poster.
» Abstract
Huberman-Rudnick universal Lyapunov scaling law is a kind of predictive capability if a system become chaotic through a sequence of period doublings, then one could be predict how it will be as a function of the control parameter [R.C. Hilborn, Chaos and Nonlinear Dynamics (Oxford University Press, NewYork, 1994)]. As chaos threshold is approached within this scaling law, it is fact that the Lyapunov exponent exhibits a power law behaviour depending to distance of the chaos threshold (a-a_c) possessing the a universal critical exponent with ν=ln2/ln δ ∼ 0.45 (δ is the Feigenbaum constant) [B. A. Huberman and J. Rudnick, PRL 45 (1980) 154]. We numerically introduce the new relationships and scaling laws related to the recurrence rate (RR), the average length of the diagonal lines (L), the determinism (DET) and exponantial divergence of phase space trajectory (DIV) as the measures from Recurrence Quantification Analysis (RQA) as we approach the chaos threshold of the logistic map with Huberman-Rudnick scaling procedure. After we determine the critical values (RRc, Lc, DETc and DIVc) of these measures on the chaos threshold of the logistic map, firstly, we verify that a scaling law of type RR-RRc ∝ (a-ac)α is evident with the critical exponent α=0.47 ±0.01. Secondly, we show that the quantity L scales as | L - Lc| ∝ (a-ac)β, where the exponent is β=0.46±0.02. Thirdly, we numerically verify that DET exhibits a scaling law of type | DET-DETc | ∝ (a-ac)γ, where the exponent is γ=0.27±0.01. Finally, we numerically show relation between the Huberman-Rudnick universal Lyapunov scaling law and Divergence scaling which behaves as λ ∝ (a-ac)ν and DIV ∝ (a-ac)κ, where the exponents are ν=0.449±0.001 and κ=0.50±0.03.
B. Goswami, A. Rheinwalt, N. Boers, N. Marwan, J. Heitzig, S. Breitenbach, J. Kurths:
Detecting paleoclimate transitions of the East Asian Summer Monsoon with recurrence networks,
NetSci2015,
Zaragoza (Spain),
June 1-5, 2015,
Talk.
N. Boers, B. Bookhagen, H. Barbosa, N. Marwan, J. Kurths, J. Marengo:
Prediction of the most extreme rainfall events in the South American Andes: A statistical forecast based on complex networks,
EGU General Assembly,
Vienna (Austria),
April 12-17, 2015,
Talk.
» Abstract
During the monsoon season, the subtropical Andes in South America are exposed to spatially extensive extreme rainfall events that frequently lead to flashfloods and landslides with severe socio-economic impacts. Since dynamical weather forecast has substantial problems with predicting the most extreme events (above the 99th percentile), alternative forecast methods are called for. Based on complex network theory, we developed a general mathematical framework for statistical prediction of extreme events in significantly interrelated time series. The key idea of our approach is to make the internal synchronization structure of extreme events mathematically accessible in terms of the topology of a network which is constructed from measuring the synchronization of extreme events at different locations. The application of our method to high-spatiotemporal resolution rainfall data (TRMM 3B42) reveals a migration pattern of large convective systems from southeastern South America towards the Argentinean and Bolivian Andes, against the direction of the northwesterly low-level moisture flow from the Amazon Basin. Once these systems reach the Andes, they lead to spatially extensive extreme events up to elevations above 4000m, leading to substantial risks of associated natural hazards. Based on atmospheric composites, we could identify an intricate interplay of frontal systems approaching from the South, low-level moisture flow from the Amazon Basin to the North, and the Andean orography as responsible climatic mechanism. These insights allow to formulate a simple forecast rule predicting 60% (90% during El Niño conditions) of extreme rainfall events at the eastern slopes of the subtropical Andes. The rule can be computed from readily available rainfall and pressure data and is already being tested by local institutions for disaster preparation.
N. Boers, N. Marwan, H. Barbosa, J. Kurths:
How Amazonian deforestation can alter the South American circulation regime: Insights from a non-linear moisture transport model,
EGU General Assembly,
Vienna (Austria),
April 12-17, 2015,
Talk.
» Abstract
A key driver of South American climate are the low-level trade winds from the tropical Atlantic Ocean towards the continent. After crossing the Amazon Basin, they are blocked by the Andes mountain range, and forced southward to the subtropics. These winds are crucial for the atmospheric moisture supply in most parts of South America. In particular, the hydrology of the two largest river basins of the Continent, namely the Amazon and the La Plata Basins, strongly depend on the moisture inflow provided by the trade winds. In turn, the Amazon rainforest can be assumed to have a strong influence on this low-level moisture circulation over South America by exchanging moisture with the atmosphere through precipitation and evapotranspiration. A pronounced positive feedback in this context is established through precipitation-induced release of latent heat over the Amazon Basin, which significantly enhances the moisture inflow from the tropical Atlantic Ocean toward the continent and can thus be considered to be crucial for the existence of today’s South American climate. Ongoing deforestation and resulting reduction in evapotranspiration rates in particular in the eastern Amazon carry the risk of a strongly nonlinear response in these interactions with the low-level atmosphere. We propose a simple differential transport model describing the cascading moisture transport from the eastern coast of South America across the Amazon Basin to the Andes, taking into account the nonlinearity associated with the release of latent heat. The results of the model suggest that the system is indeed very sensitive to relatively small reductions of the evapotranspiration rates in the eastern Amazon Basin. These reductions increase river runoff, but limit the moisture availability farther west. This leads to a reduction in precipitation rates and thereby diminishes the release of latent heat which, in turn, reduces the overall moisture inflow. We show that, according to our model, there exist critical thresholds on the spatial extents and intensities of deforestation. Beyond these thresholds, the positive feedback between the Amazon rainforest and the low-level circulation would collapse, resulting in substantial reductions in moisture available for precipitation in the western part of the Amazon Basin and further downstream of the low-level flow, including most of subtropical South America.
D. Eroglu, I. Ozken, T. Stemler, N. Marwan, K. -H. Wyrwoll, J. Kurths:
How to analyse irregularly sampled geophysical time series?,
EGU General Assembly,
Vienna (Austria),
April 12-17, 2015,
Talk.
» Abstract
One of the challenges of time series analysis is to detect dynamical changes in the dynamics of the underlying system.There are numerous methods that can be used to detect such regime changes in regular sampled times series. Here we present a new approach, that can be applied, when the time series is irregular sampled. Such data sets occur frequently in real world applications as in paleo climate proxy records.
The basic idea follows Victor and Purpura [1] and considers segments of the time series. For each segment we compute the cost of transforming the segment into the following one. If the time series is from one dynamical regime the cost of transformation should be similar for each segment of the data. Dramatic changes in the cost time series indicate a change in the underlying dynamics. Any kind of analysis can be applicable to the cost time series since it is a regularly sampled time series. While recurrence plots are not the best choice for irregular sampled data with some measurement noise component, we show that a recurrence plot analysis based on the cost time series can successfully identify the changes in the dynamics of the system.
We tested this method using synthetically created time series and will use these results to highlight the performance of our method. Furthermore we present our analysis of a suite of calcite and aragonite stalagmites located in the eastern Kimberley region of tropical Western Australia. This oxygen isotopic data is a proxy for the monsoon activity over the last 8,000 years. In this time series our method picks up several so far undetected changes from wet to dry in the monsoon system and therefore enables us to get a better understanding of the monsoon dynamics in the North-East of Australia over the last couple of thousand years.
N. Marwan, P. Koethur, C. Witt, S. F. M. Breitenbach, M. Sips:
Analysing the degree of replication of palaeoclimate records,
EGU General Assembly,
Vienna (Austria),
April 12-17, 2015,
» Poster (PDF, 2.97M)
.
» Abstract
Palaeoclimate proxy records (such as time series derived from ice cores or stalagmites) from the same or nearby location would be expected to represent similar climate variation. This is called replication of proxy records but is often difficult to achieve, because either the proxies are not reflecting the paleoclimate variation, external factors overprint the climate signal in the proxy record, or chronological uncertainties cause a serious mismatch between the individual records. In order to minimize the later issue and take the chronological uncertainties into account, we combine a Monte Carlo based approach (COPRA) with an ensemble based windowed cross-correlation analysis. This allows the investigation of potential replication of proxy records from a statistical perspective. We demonstrate this approach by comparing two stalagmite δ18O records from Heshang cave and Sanbao cave, both strongly influenced by the East Asian Summer Monsoon and covering the period between 9000 yr BP and 500 yrBP. We find that both proxy records reproduce well, although not perfectly. Main issues are differences between the records caused by unresolved geochemical processes influencing the U-series system and possibly kinetic fractionation in the oxygen isotope system. Overall, the proposed approach can provide a means to extract a correction function which reduces the uncertainties in the dating procedure. This method is a precursory step towards composite reconstructions that are based on multiple, replicating, time series.
N. Marwan, S. Foerster, J. Kurths:
Recurrence plot analysis of spatially extended high-dimensional dynamics,
EGU General Assembly,
Vienna (Austria),
April 12-17, 2015,
» Poster (PDF, 1.96M)
.
» Abstract
Recurrence plot based measures of complexity are capable tools for characterizing complex dynamics. We show the potential of selected recurrence plot measures for the investigation of spatially extended high-dimensional dynamics by applying them to data from the Lorenz96 model. The recurrence plot based measures are able to qualitatively characterize typical dynamical properties such as chaotic or periodic dynamics. Moreover, we demonstrate its power by analyzing satellite image time series of vegetation cover with contrasting dynamics as a spatially extended and potentially high-dimensional example from the real world.
J. Runge, J. Donges, J. Hlinka, N. Jajcay, N. Marwan, M. Palus, J. Kurths:
Quantifying causal pathways of interactions in the complex tropical climate system,
EGU General Assembly,
Vienna (Austria),
April 12-17, 2015,
Talk.
» Abstract
The focus of this work is to better understand the complex interplay between different subprocesses in the climate system, especially how tropical processes such as El Nino-Southern Oscillation (ENSO), the Indian Ocean Dipole, Tropical Atlantic Variability, and the tropical monsoons affect global climate.
Here a novel data-driven method is proposed based on: (1) a dimension reduction of the global surface pressure field yielding components that represent various known subprocesses such as ENSO or the North Atlantic Oscillation, (2) a causal reconstruction algorithm to detect which subprocesses are only indirectly interacting or are only spuriously correlated due to common drivers, and (3) measures to identify causal pathways in the reconstructed interaction network.
Two main results will be presented: (1) an hypothesis of a mechanism by which ENSO influences the Indian Monsoon within the surface pressure field. (2) In an explorative analysis it is shown that the method correctly identifies the major regions of upwelling convergence in the tropical oceans and also regions of strong downwelling. The approach provides a novel causal interaction perspective on complex spatio-temporal systems.
L. Tupikina, N. Molkentin, C. Lopez, E. Hernandez-Garcia, N. Marwan, J. Kurths:
Time-dependent flow-networks,
EGU General Assembly,
Vienna (Austria),
April 12-17, 2015,
Poster.
» Abstract
Complex networks have been successfully applied to various systems such as society, technology, and recently climate. Links in a climate network are defined between two geographical locations if the correlation between the time series of some climate variable is higher than a threshold. Therefore, network links are considered to imply information or heat exchange. However, the relationship between the oceanic and atmospheric flows and the climate network's structure is still unclear. Recently, a theoretical approach verifying the correlation between ocean currents and surface air temperature networks has been introduced, where the Pearson correlation networks were constructed from advection-diffusion dynamics on an underlying flow. Since the continuous approach has its limitations, i.e. high computational complexity and fixed variety of the flows in the underlying system, we introduce a new, method of flow-networks for changing in time velocity fields including external forcing in the system, noise and temperature-decay. Method of the flow-network construction can be divided into several steps: first we obtain the linear recursive equation for the temperature time-series. Then we compute the correlation matrix for time-series averaging the tensor product over all realizations of the noise, which we interpret as a weighted adjacency matrix of the flow-network and analyze using network measures. We apply the method to different types of moving flows with geographical relevance such as meandering flow. Analyzing the flow-networks using network measures we find that our approach can highlight zones of high velocity by degree and transition zones by betweenness, while the combination of these network measures can uncover how the flow propagates within time. Flow-networks can be powerful tool to understand the connection between system's dynamics and network's topology analyzed using network measures in order to shed light on different climatic phenomena.
B. Goswami, A. Rheinwalt, N. Boers, N. Marwan, J. Heitzig, S. F. M. Breitenbach, J. Kurths:
Different stages of the East Asian Summer Monsoon in the Holocene,
EGU General Assembly,
Vienna (Austria),
April 12-17, 2015,
Poster.
» Abstract
Paleoclimate proxy reconstructions (e.g. from lake sediments, ice cores, or stalagmites) have inherent age uncertainties resulting in a time-ordered series of correlated probability distributions rather than a precisely measured time series. Correlated errors make it challenging to analyze and extract valuable climate information from such records.
We show how (a) dynamical recurrences can be estimated in the presence of correlated noise, and (b) the modularity of recurrence networks can be used to infer dynamical transitions. We demonstrate our approach with a simple model and apply it to two isotope records from China to identify dynamical transitions of the East Asian Summer Monsoon (EASM) in the last 9000 years. Our results suggest that the Holocene EASM proceeded in four consecutive stages, becoming abruptly weaker at around 6400, 4400 and 3000 years ago. These transitions are known from earlier studies as abrupt, dry periods of weak monsoon superimposed on a long-term trend of gradual monsoon weakening. However, our results indicate that these events are critical shifts between four basins of stability of the Holocene EASM. We propose that these shifts could be triggered by small changes in the rate of solar insolation, which are then amplified by regional feedbacks.
J. Donges, I. Petrova, A. Löw, N. Marwan, J. Kurths:
How complex climate networks complement eigen techniques for the statistical analysis of climatological data,
EGU General Assembly,
Vienna (Austria),
April 12-17, 2015,
Poster.
» Abstract
Eigen techniques such as empirical orthogonal function (EOF) or coupled pattern (CP) / maximum covariance analysis have been frequently used for detecting patterns in multivariate climatological data sets. Recently, statistical methods originating from the theory of complex networks have been employed for the very same purpose of spatio-temporal analysis. This climate network (CN) analysis is usually based on the same set of similarity matrices as is used in classical EOF or CP analysis, e.g., the correlation matrix of a single climatological field or the cross-correlation matrix between two distinct climatological fields. In this study, formal relationships as well as conceptual differences between both eigen and network approaches are derived and illustrated using global precipitation, evaporation and surface air temperature data sets. These results allow us to pinpoint that CN analysis can complement classical eigen techniques and provides additional information on the higher-order structure of statistical interrelationships in climatological data. Hence, CNs are a valuable supplement to the statistical toolbox of the climatologist, particularly for making sense out of very large data sets such as those generated by satellite observations and climate model intercomparison exercises.
S. F. M. Breitenbach, H. E. Ridley, F. A. Lechleitner, Y. Asmerom, K. Rehfeld, K. M. Prufer, D. J. Kennett, V. V. Aquino, V. Polyak, B. Goswami, N. Marwan, G. H. Haug, and J. U. L. Baldini:
Volcanic forcing of the North Atlantic Oscillation over the last 2,000 years,
EGU General Assembly,
Vienna (Austria),
April 12-17, 2015,
Poster.
» Abstract
The North Atlantic Oscillation (NAO) is a principal mode of atmospheric circulation in the North Atlantic realm (Hurrell et al. 2003) and influences rainfall distribution over Europe, North Africa and North America. Although observational data inform us on multi-annual variability of the NAO, long and detailed paleoclimate datasets are required to understand the mechanisms and full range of its variability and the spatial extent of its influence. Chronologies of available proxy-based NAO reconstructions are often interdependent and cover only the last 1,100 years, while longer records are characterized by low sampling resolution and chronological constraints. This complicates the reconstruction of regional responses to NAO changes.
We present data from a 2,000 year long sub-annual carbon isotope record from speleothem YOK-I from Yok Balum Cave, Belize, Central America. YOK-I has been extensively dated using U-series (Kennett et al. 2012). Monitoring shows that stalagmite δ13C in Yok Balum cave is governed by infiltration changes associated with tropical wet season rainfall. Higher (lower) δ13C values reflect drier (wetter) conditions related to Intertropical Convergence Zone position and trade winds intensity.
Comparison with NAO reconstructions (Proctor et al. 2000, Trouet et al. 2009, Wassenburg et al. 2013) reveals that YOK-I δ13C sensitively records NAO-related rainfall dynamics over Belize. The Median Absolute Deviation (MAD) of δ13C extends NAO reconstructions to the last 2,000 years and indicates that high latitude volcanic aerosols force negative NAO phases.
We infer that volcanic aerosols modify inter-hemispheric temperature contrasts at multi-annual scale, resulting in meridional relocation of the ITCZ and the Bermuda-Azores High, altering NAO and tropical rainfall patterns. Decade-long dry periods in the 11th and the late 18th century relate to major high northern latitude eruptions and exemplify the climatic response to volcanic forcing by reorganization of atmospheric circulation over the North Atlantic.
P. Martin, V. Stolbova, B. Bookhagen, N. Marwan, J. Kurths:
Topology and seasonal evolution of the network of extreme precipitation over the Indian subcontinent and Sri Lanka,
EGU General Assembly,
Vienna (Austria),
April 12-17, 2015,
Poster.
» Abstract
The Indian Summer Monsoon (ISM) is one of the active components of the global climate system, and its behavior and variability is of great interest to climate researchers around the world. Here, we examine the topology and evolution of extreme rainfall across the Indian subcontinent by constructing a complex network of extreme rainfall events in the region for three periods - pre-monsoon (March - May), ISM (June - August), and post-monsoon (October - December). Networks are constructed using a synchronization measure between grid cells for a satellite-derived data set (TRMM) and a rain-gauge interpolated data set (APHRODITE). Through the analysis of various complex network metrics, such as degree, betweenness, and average link length, we describe typical repetitive patterns in North Pakistan, the Eastern Ghats, and the Tibetan Plateau. These patterns appear during the premonsoon season, evolve during the ISM, and disappear during the post-monsoon season. Our findings suggest that these are important meteorological features that deserve further attention and may be useful for the prediction of the strength and timing of the ISM.
K. Guhathakurta, N. Marwan:
Investigating chaos in emerging and developed stock markets using recurrence network analysis,
4th India Finance Conference 2014,
Bangalore (India),
December 17-19, 2014,
Talk.
P. Martin, V. Stolbova, B. Bookhagen, N. Marwan, J. Kurths:
Topology and Seasonal Evolution of the Network of Extreme Precipitation over the Indian Subcontinent and Sri Lanka,
AGU Fall Meeting,
San Francisco (USA),
December 15-19, 2014,
Talk.
» Abstract
The Indian Summer Monsoon (ISM) is one of the active components of the global climate system, and its behavior and variability is of great interest to climate researchers around the world. Here, we examine the topology and evolution of extreme rainfall across the Indian subcontinent by constructing a complex network of extreme rainfall events in the region for three periods - pre-monsoon (March - May), ISM (June - September), and post-monsoon (October - December). Networks are constructed using a synchronization measure between grid cells for a satellite-derived data set (TRMM) and a rain-gauge interpolated data set (APHRODITE). Through the analysis of various complex network metrics, such as degree, betweenness, and average link length, we describe typical repetitive patterns in North Pakistan, the Eastern Ghats, and the Tibetan Plateau. These patterns appear during the pre-monsoon season, evolve during the ISM, and disappear during the post-monsoon season. Our findings suggest that these are important meteorological features that deserve further attention and may be useful for the prediction of the strength and timing of the ISM.
N. Marwan:
Modern Approaches for Nonlinear Time Series Analysis,
Seminar at the Institute for Geology, Mineralogy and Geophysics, Ruhr University Bochum,
Bochum (Germany),
December 3, 2014,
Lecture.
N. Marwan:
Recurrence plots for time series analysis,
Workshop on Network Analysis and Data Driven Modelling of the Climate,
Potsdam (Germany),
October 27-29, 2014,
Talk.
N. Marwan:
Caves as Scientific Archives,
PIK Science & Pretzels,
Potsdam (Germany),
October 8, 2014,
Lecture.
A. Rheinwalt, B. Goswami, N. Boers, J. Heitzig, N. Marwan, R. Krishnan, J. Kurths:
A network of networks approach to investigate the influence of sea surface temperature variability on monsoon systems,
International Symposium Topical Problems of Nonlinear Wave Physics,
Nizhny Novgorod/Saratov (Russia),
July 17-23, 2014,
Talk.
» Abstract
In this study we analyze large-scale inter-dependences between Sea Surface Temperature (SST) and rainfall variability using climate networks. On account of this analysis, we coarse-grain gridded SST and rainfall datasets by merging grid points that are dynamically similar to each other. We consider the SST and rainfall systems as two distinct climate networks and use established cross-network measures to understand their interrelations. As a first step, the spatial distributions of these cross-network measures illustrate regions which are of particular importance in the interaction between SST and rainfall. Secondly, we go into further detail by investigating the cross-network topology explicitly for these regions. Here, strong influences from regions in the SST system in relation to other regions in the rainfall system are detected. These influences structured in a spatially embedded directed network describe important mechanisms behind monsoon systems. For example, behind the Indian Summer Monsoon, which is known to be controlled by SST variability over the adjacent Indian Ocean.
N. Boers, B. Bookhagen, H. M. J. Barbosa, N. Marwan, J. A. Marengo, J. Kurths:
A Complex Network approach to investigate the spatiotemporal co-variability of extreme rainfall,
4th International Workshop on Climate Informatics,
Boulder (USA),
September 25-26, 2014,
» Poster (PDF, 15.60M)
.
» Abstract
The analysis of spatial patterns of covariability of extreme rainfall is challenging because traditional, PCA-based techniques only involve the first two statistical moments of the data distribution, and are thus not able to capture the behavior in the right tails of the distributions. Here, we propose an alternative to these techniques which is based on the combination of a non-linear synchronization measure and complex network theory. This approach allows to derive spatial patterns encoding the co-variability of extreme rainfall at different locations. By introducing suitable network measures, the methodology can be used to perform climatological analysis, dataset and model intercomparisons, as well as prediction of extreme rainfall events. We introduce the methodological framework as well as applications to highspatiotemporal resolution rainfall data (TRMM 3B42) over South America.
R. Donner, J. Feldhoff, J. Donges, N. Marwan, J. Kurths:
Multivariate Extensions of Recurrence Networks Reveal Geometric Signatures of Coupling Between Nonlinear Systems,
NOLTA Conference,
Luzern (Switzerland),
September 14-18, 2014,
Talk.
» Abstract
Recurrence networks have recently proven their great potential for characterizing important properties of dynamical systems. However, in the real-world such systems typically do not evolve completely isolated from each other, but exhibit mutual interactions with their neighborhood. Here, we extend the recent view on isolated systems towards an coupled network approach to interacting systems. Specifically, we illustrate how to modify the concept of recurrence networks for studying dynamical interrelationships between two or more coupled nonlinear dynamical systems exclusively based on their attractors' geometric structures in phase space.
T. Rawald, M. Sips, N. Marwan, D. Dransch:
Fast Recurrence Quantification Analysis on GPUs,
NOLTA Conference,
Luzern (Switzerland),
September 14-18, 2014,
Talk.
» Abstract
We present a novel computing approach for recurrence quantification analysis (RQA). We refer to the concepts of "Divide and Recombine" subdividing a recurrence matrix into multiple sub matrices, computing basic RQA measures for each sub matrix in a massively parallel manner, and recombining the individual results into a valid global RQA result. Providing an implementation running on multiple graphics card processors at the same time, we are able to reduce the runtime for analyzing a temperature profile consisting of more than one million data points from over six hours to roughly five minutes.
N. Marwan, J. Donges, R. Donner, J. Kurths:
Recurrence plots and complex networks for time series analysis,
International Symposium Topical Problems of Nonlinear Wave Physics,
Nizhny Novgorod/Saratov (Russia),
July 17-23, 2014,
Talk invited.
» Abstract
Recurrence plots and derived techniques are powerful and modern time series analysis tools with a wide applicability. Recent developments have linked recurrence plots with complex networks analysis, thus, providing new and complementary measures of complexity for time series analysis. The complex network measures are related to geometrical and topological properties of the phase space representation of the dynamics. Recent applications have demonstrated the potential for data classification (e.g., for medical diagnosis), transition analysis (e.g., for detecting paleoclimatic regime transitions), or coupling analysis (e.g., for identifying coupling directions or indirect couplings).
N. Boers, B Bookhagen, H. Barbosa, N. Marwan, J. Kurths, J. Marengo:
Prediction of extreme floods in the Central Andes by means of Complex Networks,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Talk.
» Abstract
Based on a non-linear synchronisation measure and complex network theory, we present a novel framework for the prediction of extreme events of spatially embedded, interrelated time series. This method is general in the sense that it can be applied to any type of spatially sampled time series with significant interrelations, ranging from climate observables to biological or stock market data.
In this presentation, we apply our method to extreme rainfall in South America and show how this leads to the prediction of more than 60% (90% during El Niño conditions) of extreme rainfall events in the eastern Central Andes of Bolivia and northern Argentina, with only 1% false alarms. From paleoclimatic to decadal time scales, the Central Andes continue to be subject to pronounced changes in climatic conditions. In particular, our and past work shows that frequency as well as magnitudes of extreme rainfall events have increased significantly during past decades, calling for a better understanding of the involved climatic mechanisms. Due to their large spatial extend and occurrence at high elevations, these extreme events often lead to severe floods and landslides with disastrous socioeconomic impacts. They regularly affect tens of thousands of people and produce estimated costs of the order of several hundred million USD.
Alongside with the societal value of predicting natural hazards, our study provides insights into the responsible climatic features and suggests interactions between Rossby waves in polar regions and large scale (sub-)tropical moisture transport as a driver of subseasonal variability of the South American monsoon system. Predictable extreme events result from the propagation of extreme rainfall from the region of Buenos Aires towards the Central Andes given characteristic atmospheric conditions. Our results indicate that the role of frontal systems originating from Rossby waves in polar latitudes is much more dominant for controlling extreme rainfall in subtropical South America than has been assumed so far: These cold fronts cause abundant rainfall in southeastern South America, but they also dictate the direction of low-level flow from the Amazon to the subtropics. The low-level flow provides moisture for extreme rainfall which subsequently propagates from the La Plata Basin towards the eastern slopes of the Central Andes. These events become particularly interesting in view of the important role of tropical to extra-tropical couplings for assessing impacts of global warming to regional climate systems.
J. Donges, R. Donner, N. Marwan, S. Breitenbach, K. Rehfeld, J. Kurths:
A multi-archive method for robustly detecting nonlinear regime shifts in palaeoclimate dynamics under consideration of dating uncertainties,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Talk.
» Abstract
The number of available high-resolution palaeoclimate records, e.g., such as those from speleothems, is currently growing quickly. This wealth of data calls for robust methods of time series analysis that can detect common signals in a set of several records and at the same time cope with dating uncertainties and irregular sampling. In this contribution, a multi-archive method is proposed for identifying nonlinear regime shifts that are common to several palaeoclimate records. In a first step, the COPRA framework (Breitenbach et al., Clim. Past 8, 1765-1779, 2012) is applied to convert each individual irregularly sampled record into an ensemble of regularly sampled time series that are all consistent with the given dating uncertainties. Next, recurrence network analysis (RNA, Donges et al., PNAS 108, 20422-20427, 2011) is used to detect epochs with significant nonlinear deviations from the dominant dynamical regime for all records and ensemble members. Finally, we employ a Monte Carlo approach to identify the relevant time periods during which a significant fraction of all available records shows a regime shift according to RNA. We apply the proposed methodology for revealing continental-scale nonlinear transitions in the Asian monsoon system during the Holocene based on oxygen isotope records from speleothems.
J. Donges, R. Donner, N. Marwan, S. Breitenbach, K. Rehfeld, J. Kurths:
Regime shifts in Holocene Asian monsoon dynamics inferred from speleothems: Potential impacts on cultural change and migratory patterns,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Talk.
» Abstract
The Asian monsoon system has been recognized as an important potential tipping element in Earth’s climate. A global warming-driven change in monsoonal circulation, potentially towards a drier and more irregular regime, would profoundly affect up to 60% of the global human population. Hence, to improve our understanding of this major climate system, it is mandatory to investigate evidence for nonlinear transitions in past monsoonal dynamics and the underlying mechanisms that are contained in the available palaeoclimatic record. For this purpose, speleothems are among the best available high-resolution archives of Asian palaeomonsoonal variability during the Holocene and well beyond.
In this work, we apply recurrence networks, a recently developed technique for nonlinear time series analysis of palaeoclimate data (Donges et al., PNAS 108, 20422-20427, 2011), for detecting episodes with pronounced changes in Asian monsoon dynamics during the last 10 ka in oxygen isotope records from spatially distributed cave deposits covering the different branches of the Asian monsoon system. Our methodology includes multiple archives, explicit consideration of dating uncertainties with the COPRA approach and rigorous significance testing to ensure the robust detection of continental-scale changes in monsoonal dynamics.
We identify several periods characterised by nonlinear changes in Asian monsoon dynamics (e.g., ∼0.5, 2.2-2.8, 3.6-4.1, 5.4-5.7, and 8.0-8.5 ka before present [BP]), the timing of which suggests a connection to extra-tropical Bond events and rapid climate change (RCC) episodes during the Holocene. Interestingly, we furthermore detect an epoch of significantly increased regularity of monsoonal variations around 7.3 ka BP, a timing that is consistent with the typical 1.0-1.5 ka periodicity of Bond events but has been rarely reported in the literature so far. Furthermore, we find that the detected epochs of nonlinear regime shifts in Asian monsoon dynamics partly coincide with known major periods of migration, pronounced cultural changes, and the collapse of ancient human societies from the archaeological record. These findings point to a possible causal mechanism, which indicates that also future changes in monsoonal dynamics could significantly contribute to potentially severe socio-economic impacts of climate change in the Asian monsoon domain.
J. Hlinka, D. Hartman, N. Jajcay, M. Vejmelka, R. Donner, N. Marwan, J. Kurths, M. Palus:
Regional and inter-regional effects in evolving climate network,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Talk.
» Abstract
Real-world systems composed of many interacting subsystems are frequently studied as complex networks. Studied systems are thus represented by graphs composed of nodes standing for the subsystems and edges denoting interactions present among the nodes; the characteristic properties of the graph are subsequently studied and related to the system’s behavior. Potential time-dependency of edges is conveniently captured in so-called evolving networks. There is a growing interest in the application of complex network analysis approach to climate data. Use of evolving networks is a promising technique in this research area due to non-stationarity of the climate dynamics. Recently, it has been shown that an evolving climate network can be used to disentangle different types of El-Nino episodes described in the literature. In particular, an evolving network was constructed as thresholded correlation matrix of a year-long daily surface air temperature data from the NCEP/NCAR reanalysis dataset remapped onto a 10242-point equidistant geodesic grid. The time evolution of several graph characteristics, including density, clustering coefficient or average path length, has been compared with the intervals of El Niño and La Niña episodes. In the current study we identify the sources of the evolving network characteristics by considering a reduceddimensionality description of the climate system. First, we have used low density geodesic grid remapping as well as rotated principal component analysis to define the network nodes. In a more detailed analysis, the uncovered components were used to segment the whole globe surface into 68 regions. The time evolution of temperature correlation structures in local intra-component networks was studied and compared to evolving inter-component connectivity. This detailed analysis showed that the evolution of graph properties of the global network can be mostly attributed to the evolution of the intra-regional connectivity of the ENSO area and adjacent tropical regions and of the inter-regional connectivity between those.
N. Marwan, S. Breitenbach, B. Plessen, D. Scholz, J. Leonhardt:
Recurrence properties as signatures for abrupt climate change,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Talk.
» Abstract
The study of recurrence properties of dynamical systems has been shown to be very successful in characterising typical dynamical behaviour, finding regime transitions, or detecting couplings and synchronisations, even for short, noisy, and nonstationary data, as typical in Earth Sciences. Recurrence plots and their quantifications are powerful techniques for the investigation of recurrence and increasingly attract attention in recent years.
We demonstrate the potential of the newly introduced extension of recurrence plot analysis by complex network measures for the detection of abrupt dynamical changes. This method is applied on a Holocene palaeoclimate data set from Central Europe derived from a stalagmite from Blessberg Cave, Thuringia, Germany. The stalagmite δ18O proxy record covers the middle to late Holocene (6000–400 years BP). Dating uncertainties are considered by an ensemble approach derived from the COPRA framework. Characteristic changes in the recurrence properties reflecting regular dynamics coincide well with the occurrence of the Bond events 1, 2, and 3. During Bond events the Central European climate variability appears more regular. The analysis presented here examplifies the potency of quantitative recurrence methods in detecting climatic events, which otherwise remain hidden in the raw proxy time series.
D. Eroglu, N. Marwan, S. Prasad, J. Kurths:
How to use recurrence networks for geophysical time series analysis?,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Poster.
» Abstract
Recurrence plot based recurrence networks are an approach to analyze time series using complex networks theory. In both approaches, recurrence plot and the recurrence networks, we define a threshold to identify recurrent states. The selection of the threshold is very important for analysing the time series correctly via recurrence networks approach. In this talk we contribute a novel method to choose a threshold adaptively for time series. We show comparison between constant threshold and adaptive threshold cases to study transitions in the dynamics due to a change in the control parameters. This novel methods enables us to identify climate transitions from a lake sediment record.
T. Rawald, M. Sips, N. Marwan, D. Dransch:
Fast computation of recurrences in long time series,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Poster.
» Abstract
The quadratic time complexity of calculating basic RQA measures, doubling the size of the input time series leads to a quadrupling in operations, impairs the fast computation of RQA in many application scenarios. As an example, we analyze the Potsdamer Reihe, an ongoing non-interrupted hourly temperature profile since 1893, consisting of 1,043,112 data points. Using an optimized single-threaded CPU implementation this analysis requires about six hours. Our approach conducts RQA for the Potsdamer Reihe in five minutes.
We automatically split a long time series into smaller chunks (Divide) and distribute the computation of RQA measures across multiple GPU devices. To guarantee valid RQA results, we employ carryover buffers that allow sharing information between pairs of chunks (Recombine). We demonstrate the capabilities of our Divide and Recombine approach to process long time series by comparing the runtime of our implementation to existing RQA tools.
We support a variety of platforms by employing the computing framework OpenCL. Our current implementation supports the computation of standard RQA measures (recurrence rate, determinism, laminarity, ratio, average diagonal line length, trapping time, longest diagonal line, longest vertical line, divergence, entropy, trend) and also calculates recurrence times. To utilize the potential of our approach for a number of applications, we plan to release our implementation under an Open Source software license. It will be available at http://www.gfz-potsdam.de/fast-rqa/.
Since our approach allows to compute RQA measures for a long time series fast, we plan to extend our implementation to support multi-scale RQA.
L. Tupikina, N. Molkenthin, N. Marwan, J. Kurths:
Flow networks for ocean currents,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Poster.
» Abstract
Complex networks have been successfully applied to various systems such as society, technology, and recently climate. Links in a climate network are defined between two geographical locations if the correlation between the time series of some climate variable is higher than a threshold. Therefore, network links are considered to imply heat exchange. However, the relationship between the oceanic and atmospheric flows and the climate network’s structure is still unclear. Recently, a theoretical approach verifying the correlation between ocean currents and surface air temperature networks has been introduced, where the Pearson correlation networks were constructed from advection-diffusion dynamics on an underlying flow. Since the continuous approach has its limitations, i.e. by its high computational complexity, we here introduce a new, discrete construction of flow-networks, which is then applied to static and dynamic velocity fields. Analyzing the flow-networks of prototypical flows we find that our approach can highlight the zones of high velocity by degree and transition zones by betweenness, while the combination of these network measures can uncover how the flow propagates within time. We also apply the method to time series data of the Equatorial Pacific Ocean Current and the Gulf Stream ocean current for the changing velocity fields, which could not been done before, and analyse the properties of the dynamical system. Flow-networks can be powerful tools to theoretically understand the step from system’s dynamics to network’s topology that can be analyzed using network measures and is used for shading light on different climatic phenomena.
A. Rheinwalt, B. Goswami, N. Boers, J. Heitzig, N. Marwan, J. Kurths:
A network of networks approach to investigate the influence of sea surface temperature variability on monsoon systems,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Poster.
B. Goswami, J. Heitzig, N. Marwan, J. Kurths:
Recurrence networks, age uncertainties, and joint distributions: An analysis motivated by problems in paleoclimate studies,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Poster.
» Abstract
Chronological uncertainties are a typical feature of paleoclimate datasets, making their analysis non-trivial. Recent advances in age modeling and paleoclimate proxy reconstruction methodologies have enabled the representation of age uncertainties as uncertainties in the proxy value, resulting in an error-free time axis of measurement. This results in an ensemble of proxy measurements which are effectively described by a marginal probability density of the proxy value at each time instant of consideration.
However, most complex systems methods that characterize time series are not easily applicable to analyze such probability densities, primarily because the covariance between the densities at any two given points is typically unavailable. The knowledge of the covariance (or equivalently the joint probability) of the proxy values, at pairs of time points of interest, is crucial for estimating system characteristics such as the auto-correlation, power spectrum, and recurrence properties.
We present here some ideas on how to analyze such a given sequence of time-ordered probability densities in the absence of the covariance information between the densities at two distinct times. We focus mainly on the recurrence properties of the system and estimate relevant recurrence network measures along with their uncertainties of estimation. We also present a test case where the joint distribution of such probability densities could be made available and how one can analyze the recurrences in such a case.
Y. Zou, R. Donner, N. Marwan, M. Small, J. Kurths:
Long-term changes in the North-South asymmetry of solar activity: a nonlinear dynamics characterization using visibility graphs,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Poster.
» Abstract
Solar activity is characterized by complex dynamics superimposed to an almost periodic, about 11-years cycle. One of its main features is the presence of a marked, time-varying hemispheric asymmetry, the deeper reasons of which have not yet been completely uncovered. Traditionally, this asymmetry has been studied by considering amplitude and phase differences. Here, we use visibility graphs, a novel tool of nonlinear time series analysis, to obtain complementary information on hemispheric asymmetries in dynamical properties. Our analysis provides deep insights into the potentials and limitations of this method, revealing a complex interplay between factors relating to statistical and dynamical properties, i.e. effects due to the probability distribution and the regularity of observed fluctuations. We demonstrate that temporal changes in the hemispheric predominance of the graph properties lag those directly associated with the total hemispheric sunspot areas. Our findings open a new dynamical perspective on studying the North-South sunspot asymmetry, which is to be further explored in future work.
V. Stolbova, B. Bookhagen, N. Marwan, J. Kurths:
Geographic patterns of networks derived from extreme precipitation over the Indian subcontinent,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Poster.
» Abstract
Complex networks (CN) and event synchronization (ES) methods have been applied to study a number of climate phenomena such as Indian Summer Monsoon (ISM), South-American Monsoon, and African Monsoon. These methods proved to be powerful tools to infer interdependencies in climate dynamics between geographical sites, spatial structures, and key regions of the considered climate phenomenon. Here, we use these methods to study the spatial temporal variability of the extreme rainfall over the Indian subcontinent, in order to filter the data by coarse-graining the network, and to identify geographic patterns that are signature features (spatial signatures) of the ISM. We find four main geographic patterns of networks derived from extreme precipitation over the Indian subcontinent using up-to-date satellite-derived, and high temporal and spatial resolution rain-gauge interpolated daily rainfall datasets. In order to prove that our results are also relevant for other climatic variables like pressure and temperature, we use re-analysis data provided by the National Center for Environmental Prediction and Na- tional Center for Atmospheric Research (NCEP/NCAR). We find that two of the patterns revealed from the CN extreme rainfall analysis coincide with those obtained for the pressure and temperature fields, and all four above mentioned patterns can be explained by topography, winds, and monsoon circulation. CN and ES enable to select the most informative regions for the ISM, providing realistic description of the ISM dynamics with fewer data, and also help to infer geographic pattern that are spatial signatures of the ISM. These patterns deserve a special attention for the meteorologists and can be used as markers of the ISM variability.
N. Boers, B. Bookhagen, N. Marwan, J. Kurths, J. Marengo:
Complex networks identify spatial patterns of extreme rainfall of the South American monsoon system,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Poster.
» Abstract
In this study, we investigate the spatial characteristics of extreme rainfall synchronicity of the South American Monsoon System (SAMS) by means of Complex Networks. We first show how this approach leads to the identification of linkages between large-scale atmospheric conditions and natural hazards occurring at the earth’s surface. Thereafter, we exemplify how our methodology can be used to compare different datasets and to test the performance of climate models.
In recent years, complex networks have attracted great attention for analyzing the spatial characteristics of interrelations of various time series. Outstanding examples in this context are functional brain networks as well as so-called climate networks. In most approaches, the basic idea is to represent time series at different locations by network nodes, which will be connected by network links if the corresponding time series behave similar. Information on the spatial characteristics of these similarities can be inferred by network measures quantifying different aspects of the networks’ topology.
By combining several network measures and interpreting them in a climatic context, we investigate climatic linkages and classify the spatial characteristics of extreme rainfall synchronicity. Although our approach is based on only one variable (high spatiotemporal resolution rainfall), it reveals the most important features of the SAMS, such as the main moisture pathways, areas with frequent development of Mesoscale Convective Systems, and the major convergence zones. We will show that these features are only partially reproduced by reanalysis and (regional and global) climate model data.
V. Stolbova, B. Bookhagen, N. Marwan, J. Kurths:
Indian Monsoon: complex network analysis, spatial patterns and the prospects for prediction,
EGU General Assembly,
Vienna (Austria),
April 27-May 2, 2014,
Poster.
» Abstract
The Indian Summer Monsoon (ISM) is a global climate phenomenon that affects half of the world’s population. The prediction of the Indian Summer Monsoon rainfall and its extremes remains an important concern. In our study we aim to determine spatial distribution of patterns of extreme rainfall and their synchronization, because the understanding of the structure of the spatial heterogeneity of extreme rainfall is crucial for Indian agriculture and economy.
We use complex networks to identify dominant spatial patterns that govern the organization of extreme rainfall during the ISM season. We construct networks of extreme rainfall events during the ISM, the pre-monsoon, and the post-monsoon period from satellite-derived (TRMM, Tropical Rainfall Measurement Mission, product 3B42 V7) and rain-gauge interpolated (APHRODITE) datasets. 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 that can be used as indicators of the ISM variability: North Pakistan (NP), Western Ghats (WG), Eastern Ghats (EG), and Tibetan Plateau (TP). These patterns appear during the pre-monsoon season, evolve during the ISM season, and disappear during the post-monsoon season. We compare obtained results with wind fields, temperature, and pressure networks in this region derived from re-analysis data provided by the National Center for Environmental Prediction and National Center for Atmospheric Research (NCEP/NCAR). The areas of Eastern Ghats, Western Ghats, and Tibetan Plateau were previously known as areas that influence the ISM dynamics. These patterns occur because of the intricate topography of this region. The Western Ghats pattern, specifically, the Kerala region, is commonly used by climatologists for the prediction of the onset of the ISM (Pai and Nair, 2009). However, North Pakistan has not widely been considered as an important region in the analysis of the ISM. We have identified a pattern in North Pakistan that plays an important role in the extreme rainfall organization during the ISM, because it is strongly influenced by winter westerlies during the pre-monsoon season and monsoon season, which strongly influence the ISM variability. Accordingly, this pattern may serve as a marker of winter westerlies. We suggest that our obtained spatial patterns are important meteorological features that may be useful for prediction of the Indian Summer Monsoon.