Refine
Year of publication
Document Type
- Article (285)
- Monograph/Edited Volume (8)
- Postprint (7)
- Preprint (6)
- Other (1)
Is part of the Bibliography
- yes (307) (remove)
Keywords
- Complex networks (5)
- Event synchronization (4)
- precipitation (3)
- synchronization (3)
- Amazon rainforest (2)
- Extreme rainfall (2)
- Synchronization (2)
- channel (2)
- classification (2)
- climate networks (2)
Institute
- Institut für Physik und Astronomie (227)
- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (45)
- Institut für Geowissenschaften (21)
- Department Psychologie (14)
- Institut für Biochemie und Biologie (4)
- Department Linguistik (3)
- Extern (3)
- Institut für Informatik und Computational Science (2)
- Mathematisch-Naturwissenschaftliche Fakultät (2)
- Department Sport- und Gesundheitswissenschaften (1)
Recurrence plots exhibit line structures which represent typical behaviour of the investigated system. The local slope of these line structures is connected with a specific transformation of the time scales of different segments of the phase-space trajectory. This provides us a better understanding of the structures occurring in recurrence plots. The relationship between the time-scales and line structures are of practical importance in cross recurrence plots. Using this relationship within cross recurrence plots, the time-scales of differently sampled or time- transformed measurements can be adjusted. An application to geophysical measurements illustrates the capability of this method for the adjustment of time-scales in different measurements. (C) 2005 Elsevier B.V. All rights reserved
Three-dimensional quantification of structures in trabecular bone using measures of complexity
(2009)
The study of pathological changes of bone is an important task in diagnostic procedures of patients with metabolic bone diseases such as osteoporosis as well as in monitoring the health state of astronauts during long-term space flights. The recent availability of high-resolution three-dimensional (3D) imaging of bone challenges the development of data analysis techniques able to assess changes of the 3D microarchitecture of trabecular bone. We introduce an approach based on spatial geometrical properties and define structural measures of complexity for 3D image analysis. These measures evaluate different aspects of organization and complexity of 3D structures, such as complexity of its surface or shape variability. We apply these measures to 3D data acquired by high-resolution microcomputed tomography (mu CT) from human proximal tibiae and lumbar vertebrae at different stages of osteoporotic bone loss. The outcome is compared to the results of conventional static histomorphometry and exhibits clear relationships between the analyzed geometrical features of trabecular bone and loss of bone density, but also indicate that the measures reveal additional information about the structural composition of bone, which were not revealed by the static histomorphometry. Finally, we have studied the dependency of the developed measures of complexity on the spatial resolution of the mu CT data sets.
The rescaling of geological data series to a geological reference time series is of major interest in many investigations. For example, geophysical borehole data should be correlated to a given data series whose time scale is known in order to achieve an age-depth function or the sedimentation rate for the borehole data. Usually this synchronization is performed visually and by hand. Instead of using this wiggle matching by eye, we present the application of cross recurrence plots for such tasks. Using this method, the synchronization and rescaling of geological data to a given time scale is much easier and faster than by hand.
Recurrence-plot-based time series analysis is widely used to study changes and transitions in the dynamics of a system or temporal deviations from its overall dynamical regime. However, most studies do not discuss the significance of the detected variations in the recurrence quantification measures. In this letter we propose a novel method to add a confidence measure to the recurrence quantification analysis. We show how this approach can be used to study significant changes in dynamical systems due to a change in control parameters, chaos-order as well as chaos-chaos transitions. Finally we study and discuss climate transitions by analysing a marine proxy record for past sea surface temperature. This paper is dedicated to the 25th anniversary of the introduction of recurrence plots.
Climatic changes are of major importance in landslide generation in the Argentine Andes. Increased humidity as a potential influential factor was inferred from the temporal clustering of landslide deposits during a period of significantly wetter climate, 30,000 years ago. A change in seasonality was tested by comparing past (inferred from annual-layered lake deposits, 30,000 years old) and modern (present-day observations) precipitation changes. Quantitative analysis of cross recurrence plots were developed to compare the influence of the El Nino/Southern Oscillation (ENSO) on present and past rainfall variations. This analysis has shown the stronger influence of NE trades in the location of landslide deposits in the intra-andean basin and valleys, what caused a higher contrast between summer and winter rainfall and an increasing of precipitation in La Nina years. This is believed to reduce thresholds for landslide generation in the arid to semiarid intra-andean basins and valleys.
Higher variability in rainfall and river discharge could be of major importance in landslide generation in the north-western Argentine Andes. Annual layered (varved) deposits of a landslide dammed lake in the Santa Maria Basin (26°S, 66°W) with an age of 30,000 14C years provide an archive of precipitation variability during this time. The comparison of these data with present-day rainfall observations tests the hypothesis that increased rainfall variability played a major role in landslide generation. A potential cause of such variability is the El Niño/ Southern Oscillation (ENSO). The causal link between ENSO and local rainfall is quantified by using a new method of nonlinear data analysis, the quantitative analysis of cross recurrence plots (CRP). This method seeks similarities in the dynamics of two different processes, such as an ocean-atmosphere oscillation and local rainfall. Our analysis reveals significant similarities in the statistics of both modern and palaeo-precipitation data. The similarities in the data suggest that an ENSO-like influence on local rainfall was present at around 30,000 14C years ago. Increased rainfall, which was inferred from a lake balance modeling in a previous study, together with ENSO-like cyclicities could help to explain the clustering of landslides at around 30,000 14C years ago.
Recurrence-plot-based measures of complexity and its application to heart-rate-variability data
(2002)
The knowledge of transitions between regular, laminar or chaotic behavior is essential to understand the underlying mechanisms behind complex systems. While several linear approaches are often insufficient to describe such processes, there are several nonlinear methods which however require rather long time observations. To overcome these difficulties, we propose measures of complexity based on vertical structures in recurrence plots and apply them to the logistic map as well as to heart rate variability data. For the logistic map these measures enable us not only to detect transitions between chaotic and periodic states, but also to identify laminar states, i.e. chaos-chaos transitions. The traditional recurrence quantification analysis fails to detect the latter transitions. Applying our new measures to the heart rate variability data, we are able to detect and quantify the laminar phases before a life-threatening cardiac arrhythmia occurs thereby facilitating a prediction of such an event. Our findings could be of importance for the therapy of malignant cardiac arrhythmias.
Phase synchronization is an important phenomenon that occurs in a wide variety of complex oscillatory processes. Measuring phase synchronization can therefore help to gain fundamental insight into nature. In this Letter we point out that synchronization analysis techniques can detect spurious synchronization, if they are fed with a superposition of signals such as in electroencephalography or magnetoencephalography data. We show how techniques from blind source separation can help to nevertheless measure the true synchronization and avoid such pitfalls
We present experimental and numerical evidence of synchronization of burst events in two different modulated CO2 lasers. Bursts appear randomly in each laser as trains of large amplitude spikes intercalated by a small amplitude chaotic regime. Experimental data and model show the frequency locking of bursts in a suitable interval of coupling strength. We explain the mechanism of this phenomenon and demonstrate the inhibitory properties of the implemented coupling.
We analyze synchronization between two interacting populations of different phase oscillators. For the important case of asymmetric coupling functions, we find a much richer dynamical behavior compared to that of symmetrically coupled populations of identical oscillators. It includes three types of bistabilities, higher order entrainment and the existence of states with unusual stability properties. All possible routes to synchronization of the populations are presented and some stability boundaries are obtained analytically. The impact of these findings for neuroscience is discussed.
When we fixate our gaze on a stable object, our eyes move continuously with extremely small involuntary and autonomic movements, that even we are unaware of during their occurrence. One of the roles of these fixational eye movements is to prevent the adaptation of the visual system to continuous illumination and inhibit fading of the image. These random, small movements are restricted at long time scales so as to keep the target at the centre of the field of view. In addition, the synchronisation properties between both eyes are related to binocular coordination in order to provide stereopsis. We investigated the roles of different time scale behaviours, especially how they are expressed in the different spatial directions (vertical versus horizontal). We also tested the synchronisation between both eyes. Results show different scaling behaviour between horizontal and vertical movements. When the small ballistic movements, i.e., microsaccades, are removed, the scaling behaviour in both axes becomes similar. Our findings suggest that microsaccades enhance the persistence at short time scales mostly in the horizontal component and much less in the vertical component. We also applied the phase synchronisation decay method to study the synchronisation between six combinations of binocular fixational eye movement components. We found that the vertical-vertical components of right and left eyes are significantly more synchronised than the horizontal-horizontal components. These differences may be due to the need for continuously moving the eyes in the horizontal plane in order to match the stereoscopic image for different viewing distances.
Heterogeneity in the degree (connectivity) distribution has been shown to suppress synchronization in networks of symmetrically coupled oscillators with uniform coupling strength (unweighted coupling). Here we uncover a condition for enhanced synchronization in weighted networks with asymmetric coupling. We show that, in the optimum regime, synchronizability is solely determined by the average degree and does not depend on the system size and the details of the degree distribution. In scale-free networks, where the average degree may increase with heterogeneity, synchronizability is drastically enhanced and may become positively correlated with heterogeneity, while the overall cost involved in the network coupling is significantly reduced as compared to the case of unwcighted coupling
Many complex networks display strong heterogeneity in the degree (connectivity) distribution. Heterogeneity in the degree distribution often reduces the average distance between nodes but, paradoxically, may suppress synchronization in networks of oscillators coupled symmetrically with uniform coupling strength. Here we offer a solution to this apparent paradox. Our analysis is partially based on the identification of a diffusive process underlying the communication between oscillators and reveals a striking relation between this process and the condition for the linear stability of the synchronized states. We show that, for a given degree distribution, the maximum synchronizability is achieved when the network of couplings is weighted and directed and the overall cost involved in the couplings is minimum. This enhanced synchronizability is solely determined by the mean degree and does not depend on the degree distribution and system size. Numerical verification of the main results is provided for representative classes of small-world and scale-free networks
A new globally uniform Lagrangian transport scheme for large ensembles of passive tracer particles is presented and applied to wind data from a coupled atmosphere-ocean climate model that includes interactive dynamical feedback with stratospheric chemistry. This feedback from the chemistry is found to enhance large-scale meridional air mass exchange in the northern winter stratosphere as well as intrusion of stratospheric air into the troposphere, where both effects are due to a weakened polar vortex.
We study phase synchronization effects in a chain of nonidentical chaotic oscillators with a type-I intermittent behavior. Two types of parameter distribution, linear and random, are considered. The typical phenomena are the onset and existence of global (all-to-all) and cluster (partial) synchronization with increase of coupling. Increase of coupling strength can also lead to desynchronization phenomena, i.e., global or cluster synchronization is changed into a regime where synchronization is intermittent with incoherent states. Then a regime of a fully incoherent nonsynchronous state (spatiotemporal intermittency) appears. Synchronization-desynchronization transitions with increase of coupling are also demonstrated for a system resembling an intermittent one: a chain of coupled maps replicating the spiking behavior of neurobiological networks
We demonstrate the existence of phase synchronization of two chaotic rotators. Contrary to phase synchronization of chaotic oscillators, here the Lyapunov exponents corresponding to both phases remain positive even in the synchronous regime. Such frequency locked dynamics with different ratios of frequencies are studied for driven continuous-time rotators and for discrete circle maps. We show that this transition to phase synchronization occurs via a crisis transition to a band-structured attractor.
Attractor-repeller collision and eyelet intermittency at the transition to phase synchronization
(1997)
The chaotically driven circle map is considered as the simplest model ofphase synchronization of a chaotic continuous-time oscillator by external periodic force. The phase dynamics is analyzed via phase-locking regions of the periodic cycles embedded in the strange attractor. It is shown that full synchronization, where all the periodic cycles are phase locked, disappears via the attractor-repeller collision. Beyond the transition an intermittent regime with exponentially rare phase slips, resulting from the trajectory's hits on an eyelet, is observed.
Frequent and intense rainfall events demand innovative techniques to better predict the extreme rainfall dynamics. This task requires essentially the assessment of the basic types of atmospheric processes that trigger extreme rainfall, and then to examine the differences between those processes, which may help to identify key patterns to improve predictive algorithms. We employ tools from network theory to compare the spatial features of extreme rainfall over the Japanese archipelago and surrounding areas caused by two atmospheric processes: the Baiu front, which occurs mainly in June and July (JJ), and the tropical storms from August to November (ASON). We infer from complex networks of satellite-derived rainfall data, which are based on the nonlinear correlation measure of event synchronization. We compare the spatial scales involved in both systems and identify different regions which receive rainfall due to the large spatial scale of the Baiu and tropical storm systems. We observed that the spatial scales involved in the Baiu driven rainfall extremes, including the synoptic processes behind the frontal development, are larger than tropical storms, which even have long tracks during extratropical transitions. We further delineate regions of coherent rainfall during the two seasons based on network communities, identifying the horizontal (east-west) rainfall bands during JJ over the Japanese archipelago, while during ASON these bands align with the island arc of Japan.
Reconciling the paths of extreme rainfall with those of typhoons remains difficult despite advanced forecasting techniques. We use complex networks defined by a nonlinear synchronization measure termed event synchronization to track extreme rainfall over the Japanese islands. Directed networks objectively record patterns of heavy rain brought by frontal storms and typhoons but mask out contributions of local convective storms. We propose a radial rank method to show that paths of extreme rainfall in the typhoon season (August-November, ASON) follow the overall southwest-northeast motion of typhoons and mean rainfall gradient of Japan. The associated eye-of-the-typhoon tracks deviate notably and may thus distort estimates of heavy typhoon rainfall. We mainly found that the lower spread of rainfall tracks in ASON may enable better hindcasting than for westerly-fed frontal storms in June and July.
The 300 year record of the yearly sunspot numbers and numerically generated trajectory of the solar inertial motion (SIM) were subjects of a synchronization analysis. Phase synchronization of the sunspot cycle and a fast component of the SIM have been found and confirmed with statistical significance in three epochs (1727-1757, 1802-1832 and 1863-1922) of the entire 1700-1997 record. This result can be considered as a quantitative support for the hypothesis that there is a weak interaction of gravity and solar activity.
We study possible interrelations between the 300-year record of the yearly sunspot numbers and the solar inertial motion (SIM) using the recently developed technique of synchronization analysis. Phase synchronization of the sunspot cycle and the SIM is found and statistically confirmed in three epochs (1734-1790, 1855-1875 and 1907-1960) of the whole period 1700-2000. These results give quantitative support to the hypothesis that there is a weak interaction between the solar activity and the SIM.
Analysis of blood pressure dynamics in male and female rats using the continuous wavelet transform
(2009)
We study gender-related particularities in cardiovascular responses to stress and nitric oxide (NO) deficiency in rats using HR, mean arterial pressure (MAP) and a proposed wavelet-based approach. Blood pressure dynamics is analyzed: (1) under control conditions, (2) during immobilization stress and recovery and (3) during nitric oxide blockade by N-G-nitro-L-arginine-methyl ester (L-NAME). We show that cardiovascular sensitivity to stress and NO deficiency depends upon gender. Actually, in females the chronotropic effect of stress is more pronounced, while the pressor effect is weakened compared with males. We conclude that females demonstrate more favorable patterns of cardiovascular responses to stress and more effective NO control of cardiovascular activity than males.
We show a scenario of a two-frequeney torus breakdown, in which a global bifurcation occurs due to the collision of a quasi-periodic torus T-2 with saddle points, creating a heteroclinic saddle connection. We analyze the geometry of this torus-saddle collision by showing the local dynamics and the invariant manifolds (global dynamics) of the saddle points. Moreover, we present detailed evidences of a heteroclinic saddle-focus orbit responsible for the type- if intermittency induced by this global bifurcation. We also characterize this transition to chaos by measuring the Lyapunov exponents and the scaling laws.
This paper treats a problem of reconstructing ordinary differential equation from a single analytic time series with observational noise. We suppose that the noise is Gaussian (white). The investigation is presented in terms of classical theory of dynamical systems and modern time series analysis. We restrict our considerations on time series obtained as a numerical analytic solution of autonomous ordinary differential equation, solved with respect to the highest derivative and with polynomial right-hand side. In case of an approximate numerical solution with a rather small error, we propose a geometrical basis and a mathematical algorithm to reconstruct a low-order and low-power polynomial differential equation. To reduce the noise the given time series is smoothed at every point by moving polynomial averages using the least-squares method. Then a specific form of the least-squares method is applied to reconstruct the polynomial right-hand side of the unknown equation. We demonstrate for monotonous, periodic and chaotic solutions that this technique is very efficient
We study the dynamics of the excitable Fitz Hugh-Nagumo system under external noisy driving. Noise activates the system producing a sequence of pulses. The coherence of these noise-induced oscillations is shown to be maximal for a certain noise amplitude. This new effect of coherence resonance is explained by different noise dependencies of the activation and the excursion times. A simple one-dimensional model based on the Langevin dynamics is proposed for the quantitative description of this phenomenon.
We demonstrate the occurrence of regimes with singular continuous (fractal) Fourier spectra in autonomous dissipative dynamical systems. The particular example in an ODE system at the accumulation points of bifurcation sequences associated to the creation of complicated homoclinic orbits. Two different machanisms responsible for the appearance of such spectra are proposed. In the first case when the geometry of the attractor is symbolically represented by the Thue-Morse sequence, both the continuous-time process and its descrete Poincaré map have singular power spectra. The other mechanism owes to the logarithmic divergence of the first return times near the saddle point; here the Poincaré map possesses the discrete spectrum, while the continuous-time process displays the singular one. A method is presented for computing the multifractal characteristics of the singular continuous spectra with the help of the usual Fourier analysis technique.
The transition from fully synchronized behavior to two-cluster dynamics is investigated for a system of N globally coupled chaotic oscillators by means of a model of two coupled logistic maps. An uneven distribution of oscillators between the two clusters causes an asymmetry to arise in the coupling of the model system. While the transverse period-doubling bifurcation remains essentially unaffected by this asymmetry, the transverse pitchfork bifurcation is turned into a saddle-node bifurcation followed by a transcritical riddling bifurcation in which a periodic orbit embedded in the synchronized chaotic state loses its transverse stability. We show that the transcritical riddling transition is always hard. For this, we study the sequence of bifurcations that the asynchronous point cycles produced in the saddle-node bifurcation undergo, and show how the manifolds of these cycles control the magnitude of asynchronous bursts. In the case where the system involves two subpopulations of oscillators with a small mismatch of the parameters, the transcritical riddling will be replaced by two subsequent saddle-node bifurcations, or the saddle cycle involved in the transverse destabilization of the synchronized chaotic state may smoothly shift away from the synchronization manifold. In this way, the transcritical riddling bifurcation is substituted by a symmetry-breaking bifurcation, which is accompanied by the destruction of a thin invariant region around the symmetrical chaotic state.
We propose a technique to calculate large-scale dimension densities in both higher-dimensional spatio-temporal systems and low-dimensional systems from only a few data points, where known methods usually have an unsatisfactory scaling behavior. This is mainly due to boundary and finite size effects. With our rather simple method we normalize boundary effects and get a significant correction of the dimension estimate. This straightforward approach is basing on rather general assumptions. So even weak coherent structures obtained from small spatial couplings can be detected with this method, what is impossible by using the Lyapunov-dimension density. We demonstrate the efficiency of our technique for coupled logistic maps, coupled tent maps, the Lorenz-attractor and the Roessler-attractor.
In this work, we reanalyze the heart rate variability (HRV) data from the 2002 Computers in Cardiology (CiC) Challenge using the concept of large-scale dimension densities and additionally apply this technique to data of healthy persons and of patients with cardiac diseases. The large-scale dimension density (LASDID) is estimated from the time series using a normalized Grassberger-Procaccia algorithm, which leads to a suitable correction of systematic errors produced by boundary effects in the rather large scales of a system. This way, it is possible to analyze rather short, nonstationary, and unfiltered data, such as HRV. Moreover, this method allows us to analyze short parts of the data and to look for differences between day and night. The circadian changes in the dimension density enable us to distinguish almost completely between real data and computer-generated data from the CiC 2002 challenge using only one parameter. In the second part we analyzed the data of 15 patients with atrial fibrillation (AF), 15 patients with congestive heart failure (CHF), 15 elderly healthy subjects (EH), as well as 18 young and healthy persons (YH). With our method we are able to separate completely the AF (rho(mu)(ls)=0.97 +/- 0.02) group from the others and, especially during daytime, the CHF patients show significant differences from the young and elderly healthy volunteers (CHF, 0.65 +/- 0.13; EH, 0.54 +/- 0.05; YH, 0.57 +/- 0.05; p < 0.05 for both comparisons). Moreover, for the CHF patients we find no circadian changes in rho(mu)(ls) (day, 0.65 +/- 0.13; night, 0.66 +/- 0.12; n.s.) in contrast to healthy controls (day, 0.54 +/- 0.05; night, 0.61 +/- 0.05; p=0.002). Correlation analysis showed no statistical significant relation between standard HRV and circadian LASDID, demonstrating a possibly independent application of our method for clinical risk stratification
Identifying causal relations from observational data sets has posed great challenges in data-driven causality inference studies. One of the successful approaches to detect direct coupling in the information theory framework is transfer entropy. However, the core of entropy-based tools lies on the probability estimation of the underlying variables. Herewe propose a data-driven approach for causality inference that incorporates recurrence plot features into the framework of information theory. We define it as the recurrence measure of conditional dependence (RMCD), and we present some applications. The RMCD quantifies the causal dependence between two processes based on joint recurrence patterns between the past of the possible driver and present of the potentially driven, excepting the contribution of the contemporaneous past of the driven variable. Finally, it can unveil the time scale of the influence of the sea-surface temperature of the Pacific Ocean on the precipitation in the Amazonia during recent major droughts.
The analysis of baroreflex sensitivity (BRS) and heart rate variability (HRV) leads to additional insights into patients' prognosis after cardiovascular events. The following study was performed to assess the differences in the post-operative recovery of autonomic regulation after mitral valve (MV) and aortic valve (AV) surgery with a heart lung machine. Among the 43 consecutive male patients enrolled in a prospective study, 26 underwent isolated AV surgery and 17 isolated MV surgery. Blood pressure as well as ECG signals were recorded the day before, 24 hours after and one week after surgery. BRS was calculated according to the dual sequence method, and HRV was calculated using standard linear as well as nonlinear parameters. There were no major differences between the two groups in the pre-operative values. At 24 hours a comparable depression of HRV and BRS in both groups was observed, while at 7 days there was partial recovery in AV patients, which was absent in MV patients: p(AV versus MV) < 0.001. While the response of the autonomic system to surgery is similar in AV and MV patients, there is obviously a decreased ability to recover in MV patients, probably attributed to traumatic lesions of the autonomic nervous system by opening the atria. Ongoing research is required for further clarification of the pathophysiology of this phenomenon and to establish strategies to restore autonomic function.
Non-linear time series analysis of precipitation events using regional climate networks for Germany
(2016)
Synchronous occurrences of heavy rainfall events and the study of their relation in time and space are of large socio-economical relevance, for instance for the agricultural and insurance sectors, but also for the general well-being of the population. In this study, the spatial synchronization structure is analyzed as a regional climate network constructed from precipitation event series. The similarity between event series is determined by the number of synchronous occurrences. We propose a novel standardization of this number that results in synchronization scores which are not biased by the number of events in the respective time series. Additionally, we introduce a new version of the network measure directionality that measures the spatial directionality of weighted links by also taking account of the effects of the spatial embedding of the network. This measure provides an estimate of heavy precipitation isochrones by pointing out directions along which rainfall events synchronize. We propose a climatological interpretation of this measure in terms of propagating fronts or event traces and confirm it for Germany by comparing our results to known atmospheric circulation patterns.
The investigation of foetal reaction to internal and external conditions and stimuli is an important tool in the characterization of the developing neural integration of the foetus. An interesting example of this is the study of the interrelationship between the foetal and the maternal heart rate. Recent studies have shown a certain likelihood of occasional heart rate synchronization between mother and foetus. In the case of respiratory-induced heart rate changes, the comparison with maternal surrogates suggests that the evidence for detected synchronization is largely statistical and does not result from physiological interaction. Rather, they simply reflect a stochastic, temporary stability of two independent oscillators with time-variant frequencies. We reanalysed three datasets from that study for a more local consideration. Epochs of assumed synchronization associated with short-term regulation of the foetal heart rate were selected and compared with synchronization resulting from white noise instead of the foetal signal. Using data-driven modelling analysis, it was possible to identify the consistent influence of the heartbeat duration of maternal beats preceding the foetal beats during epochs of synchronization. These maternal beats occurred approximately one maternal respiratory cycle prior to the affected foetal beat. A similar effect could not be found in the epochs without synchronization. Simulations based on the fitted models led to a higher likelihood of synchronization in the data segments with assumed foetal-maternal interaction than in the segment without such assumed interaction. We conclude that the data-driven model-based analysis can be a useful tool for the identification of synchronization.
We present a new method to detect phase as well as generalized synchronization in a wide class of complex systems. It is based on the recurrences of the system's trajectory to the neighborhood of a former state in phase space. We illustrate the applicability of the algorithm for the paradigmatic chaotic Rossler system in the funnel regime and for noisy data, where other methods to detect phase synchronization fail. Furthermore, we demonstrate for electrochemical experiments that the method can easily detect phase and generalized synchronization in non-phase- coherent and even non-stationary time series
The method of twin surrogates has been introduced to test for phase synchronization of complex systems in the case of passive experiments. In this paper we derive new analytical expressions for the number of twins depending on the size of the neighborhood, as well as on the length of the trajectory. This allows us to determine the optimal parameters for the generation of twin surrogates. Furthermore, we determine the quality of the twin surrogates with respect to several linear and nonlinear statistics depending on the parameters of the method. In the second part of the paper we perform a hypothesis test for phase synchronization in the case of experimental data from fixational eye movements. These miniature eye movements have been shown to play a central role in neural information processing underlying the perception of static visual scenes. The high number of data sets (21 subjects and 30 trials per person) allows us to compare the generated twin surrogates with the "natural" surrogates that correspond to the different trials. We show that the generated twin surrogates reproduce very well all linear and nonlinear characteristics of the underlying experimental system. The synchronization analysis of fixational eye movements by means of twin surrogates reveals that the synchronization between the left and right eye is significant, indicating that either the centers in the brain stem generating fixational eye movements are closely linked, or, alternatively that there is only one center controlling both eyes.
We propose a new approach to calculate recurrence plots of multivariate time series, based on joint recurrences in phase space. This new method allows to estimate dynamical invariants of the whole system, like the joint Renyi entropy of second order. We use this entropy measure to quantitatively study in detail the phase synchronization of two bidirectionally coupled chaotic systems and identify different types of transitions to chaotic phase synchronization in dependence on the coupling strength and the frequency mismatch. By means of this analysis we find several new phenomena, such a chaos-period-chaos transition to phase synchronization for rather large coupling strengths. (C) 2004 Elsevier B.V. All rights reserved
We investigate the relationship between precipitation and runoff data from a small forested catchment in the Harz mountains (Germany). For this purpose, we develop a conceptual model including memory effects to predict the runoff signal using the precipitation data as input. An enhanced variant of the model also includes air temperature as input variable. We show in terms of correlation functions that this model describes main dynamical properties of the runoff, especially the delay between rain event and runoff response as the annual persistence in the runoff data.
Chen et al. [Phys. Rev. E 61, 2559 (2000)] recently proposed an extension of the concept of phase for discrete chaotic systems. Using the newly introduced definition of phase they studied the dynamics of coupled map lattices and compared these dynamics with phase synchronization of coupled continuous-time chaotic systems. In this paper we illustrate by two simple counterexamples that the angle variable introduced by Chen et al. fails to satisfy the basic requirements to the proper phase. Furthermore, we argue that an extension of the notion of phase synchronization to generic discrete maps is doubtful.
In this article we review the application of the synchronization theory to the analysis of multivariate biological signals. We address the problem of phase estimation from data and detection and quantification of weak interaction, as well as quantification of the direction of coupling. We discuss the potentials as well as limitations and misinterpretations of the approach
We study synchronization transitions in a system of two coupled self-sustained chaotic oscillators. We demonstrate that with the increase of coupling strength the system first undergoes the transition to phase synchronization. With a further increase of coupling, a new synchronous regime is observed, where the states of two oscillators are nearly identical, but one system lags in time to the other. We describe thisregime as a state with correlated amplitudes and a constant phase shift. These transitions are traced in the Lyapunov spectrum.
Locking-based frequency measurement and synchronization of chaotic oscillators with complex dynamics
(2002)
During life bones constantly adapt their structure to their mechanical environment via a mechanically controlled process called bone remodeling. For trabecular bone, this process modifies the thickness of each trabecula leading occasionally to full resorption. We describe the irreversible dynamics of the trabecular thickness distribution (TTD) by means of a Markov chain discrete in space and time. By using thickness data from adult patients, we derive the transition probabilities in the chain. This allows a quantification, in terms of geometrical quantities, of the control of bone remodeling and thus to determine the evolution of the TTD with age.
Changes in trabecular bone composition during development of osteoporosis are used as a model for bone loss in microgravity conditions during a space flight. Symbolic dynamics and measures of complexity are proposed and applied to assess quantitatively the structural composition of bone tissue from 3D data sets of human tibia bone biopsies acquired by a micro-CT scanner. In order to justify the newly proposed approach, the measures of complexity of the bone architecture were compared with the results of traditional 2D bone histomorphometry. The proposed technique is able to quantify the structural loss of the bone tissue and may help to diagnose and to monitor changes in bone structure of patients on Earth as well as of the space-flying personnel. © 2005 Elsevier Ltd. All rights reserved
In this paper we apply symbolic transformations as a visualisation technique for analysing rhythm production. It is shown that qualitative information can be extracted from the experimental data. This approach may provide new insights into the organisation of temporal order by the brain on different levels of description. A simple phenomenological model for the explanation of the observed phenomena is proposed.
We investigate the cognitive control in polyrhythmic hand movements as a model paradigm for bimanual coordination. Using a symbolic coding of the recorded time series, we demonstrate the existence of qualitative transitions induced by experimental manipulation of the tempo. A nonlinear model with delayed feedback control is proposed, which accounts for these dynamical transitions in terms of bifurcations resulting from variation of the external control parameter. Furthermore, it is shown that transitions can also be observed due to fluctuations in the timing control level. We conclude that the complexity of coordinated bimanual movements results from interactions between nonlinear control mechanisms with delayed feedback and stochastic timing components.
Graphical models applying partial coherence to multivariate time series are a powerful tool to distinguish direct and indirect interdependencies in multivariate linear systems. We carry over the concept of graphical models and partialization analysis to phase signals of nonlinear synchronizing systems. This procedure leads to the partial phase synchronization index which generalizes a bivariate phase synchronization index to the multivariate case and reveals the coupling structure in multivariate synchronizing systems by differentiating direct and indirect interactions. This ensures that no false positive conclusions are drawn concerning the interaction structure in multivariate synchronizing systems. By application to the paradigmatic model of a coupled chaotic Roessler system, the power of the partial phase synchronization index is demonstrated
In the recent past, recurrence quantification analysis (RQA) has gained an increasing interest in various research areas. The complexity measures the RQA provides have been useful in describing and analysing a broad range of data. It is known to be rather robust to noise and nonstationarities. Yet, one key question in empirical research concerns the confidence bounds of measured data. In the present Letter we suggest a method for estimating the confidence bounds of recurrence-based complexity measures. We study the applicability of the suggested method with model and real- life data.
The EEG is one of the most commonly used tools in brain research. Though of high relevance in research, the data obtained is very noisy and nonstationary. In the present article we investigate the applicability of a nonlinear data analysis method, the recurrence quantification analysis (RQA), to Such data. The method solely rests on the natural property of recurrence which is a phenomenon inherent to complex systems, such as the brain. We show that this method is indeed suitable for the analysis of EEG data and that it might improve contemporary EEG analysis.
Understanding the functional dynamics of the mammalian brain is one of the central aims of modern neuroscience. Mathematical modeling and computational simulations of neural networks can help in this quest. In recent publications, a multilevel model has been presented to simulate the resting-state dynamics of the cortico-cortical connectivity of the mammalian brain. In the present work we investigate how much of the dynamical behavior of the multilevel model can be reproduced by a strongly simplified model. We find that replacing each cortical area by a single Rulkov map recreates the patterns of dynamical correlation of the multilevel model, while the outcome of other models and setups mainly depends on the local network properties, e. g. the input degree of each vertex. In general, we find that a simple simulation whose dynamics depends on the global topology of the whole network is far from trivial. A systematic analysis of different dynamical models and coupling setups is required.
Using quantities of symbolic dynamics, such as mutual information, Shannon information and algorithmic complexity, we have searched for interrelations of spikes emitted simultaneously at different frequencies during the impulsive phase of a flare event. As the spikes are related to the flare energy release and are interpreted as emissions originating at different sites having different magnetic field strengths, any relation in frequency is interpretated as a relation in space. This approach is appropriate to characterize such spatio-temporal patterns, whereas the popular estimate of fractal dimensions can be applied to low-dimensional systems only. Depending on the energy release and emission processes, two types of fragmentation are possible: a scenario of global organization (spikes are emitted in a succession of similar events by the same system) or a scenario of local organization (many systems triggered by an initial event). Mutual information which is a generalization of correlation indicates a relation in frequency beyond the bandwidth of individual spikes. The scans in the spectrograms with large mutual information also show a low level of Shannon information and algorithmic complexity, indicating that the simultaneous appearance of spikes at other frequencies is not a completely stochastic phenomenon (white noise). It may be caused by a nonlinear deterministic system or by a Markov process. By means of mutual information we find a memory over frequency intervals up to 60 MHz. Shannon information and algorithmic complexity concern the mbox{whole} frequency region, i.e. the global source region. A global organization is also apparent in quasi-periodic changes of the Shannon information and algorithmic complexity in the range of 2 - 8 seconds. The finding is compatible with a scenario of local organization in which the information of one event spreads spatially and triggers further events at different places. The region is not an ensemble of independently flashing sources, each representing a system that cascades in energy after an initial trigger. On the contrary, there is a causal connection between the sources at any time. The analysis of the four spike events suggests that the structure in frequency is not stochastic but a process in which spikes at nearby locations are simultaneously triggered by a common exciter.
We investigate the characteristics of time-delay systems in the presence of Gaussian noise. We show that the delay time embedded in the time series of time-delay system with constant delay cannot be estimated in the presence noise for appropriate values of noise intensity thereby forbidding any possibility of phase space reconstruction. We also demonstrate the existence of complete synchronization between two independent identical time-delay systems driven by a common noise without explicitly establishing any external coupling between them.
We use an index of chaotic synchronization based on the averaged coherence function for the quantitative analysis of the process of the complete synchronization loss in unidirectionally coupled oscillators and maps. We demonstrate that this value manifests different stages of the synchronization breaking. It is invariant to time delay and insensitive to small noise and distortions, which can influence the accessible signals at measurements. Peculiarities of the synchronization destruction in maps and oscillators are investigated
Climate change perception
(2018)
Climate change and variability have created widespread risks for farmers’ food and livelihood security in the Himalayas. However, the extent of impacts experienced and perceived by farmers varies, as there is substantial diversity in the demographic, social, and economic conditions. Therefore, it is essential to understand how farmers with different resource-endowment and household characteristics perceive climatic risks. This study aims to analyze how farmer types perceive climate change processes and its impacts to gain insight into locally differentiated concerns by farming communities. The present study is based in the Uttarakhand state of Indian Western Himalayas. We examine farmer perceptions of climate change and how perceived impacts differ across farmer types. Primary household interviews with farming households (n = 241) were done in Chakrata and Bhikiyasian tehsil in Uttarakhand, India. In addition, annual and seasonal patterns of historical data of temperature (1951–2013) and precipitation (1901–2013) were analyzed to estimate trends and validate farmers’ perception. Using statistical methods farmer typology was constructed, and five unique farmer types are identified. Majority of respondents across all farmer types noticed a decrease in summer and winter precipitation and an increase in summer temperature. Whereas the perceptions of impacts of climate change diverged across farmer types, as specific farmer types exclusively experienced few impacts. Impact of climatic risks on household food security and income was significantly perceived stronger by low-resource-endowed subsistence farmers, whereas the landless farmer type exclusively felt impacts on the communities social bond. This deeper understanding of the differentiated perception of impacts has strong implications for agricultural and development policymaking, highlighting the need for providing flexible adaptation options rather than specific solutions to avoid inequalities in fulfilling the needs of the heterogeneous farming communities.
Acoustic emission signals generated during high speed cutting of steel are investigated. The data are represen ted in time-folded form. Several methods from linear and nonlinear data analysis based on time- and frequency- domain are applied to the data and reveal signatures of the observed acoustic emission signal. These investiga tions are necessary for modeling the cutting process by means of differential equations.
Estimation of parameters and unobserved components for nonlinear systems from noisy time series
(2002)
We study the problem of simultaneous estimation of parameters and unobserved states from noisy data of nonlinear time-continuous systems, including the case of additive stochastic forcing. We propose a solution by adapting the recently developed statistical method of unscented Kalman filtering to this problem. Due to its recursive and derivative-free structure, this method minimizes the cost function in a computationally efficient and robust way. It is found that parameters as well as unobserved components can be estimated with high accuracy, including confidence bands, from heavily noise-corrupted data.
Charged dust grains in circumplanetary environments experience, beyond various deterministic forces, also stochastic perturbations caused, by fluctuations of the magnetic field, the charge of the grains, by chaotic rotation of aspherical grains, etc. Here we investigate the dynamics of a dust population in a circular orbit around a planet which is perturbed by a stochastic planetary magnetic field B', modeled by an isotropically Gaussian white noise. The resulting perturbation equations give rise to a modified diffusion of the inclinations i and eccentricities e. The diffusion coefficient is found to be D proportional to w^2 O /n^2 , where the gyrofrequency, the Kepler frequency, and the synodic frequency are denoted by w , O, and n, respectively. This behavior has been checked against numerical simulations. We have chosen dust grains (1 m in radius) ejected from Jupiter's satellite Europa in circular equatorial orbits around Jupiter and integrated numerically their trajectories over their typical lifetimes (100 years). The particles were exposed to a Gaussian fluctuating magnetic field B' with the same statistical properties as in the analytical treatment. These simulations have confirmed the analytical results. The theoretical studies showed the statistical properties of B' to be of decisive importance. To estimate them, we analyzed the magnetic field data obtained by the Galileo spacecraft magnetometer at Jupiter and found almost Gaussian fluctuations of about 5% of the mean field and exponentially decaying correlations. This results in a diffusion of orbital inclinations and eccentricities of the dust grains of about ten percent over the lifetime of the particles. For smaller dusty motes or for close-in particles (e.g., in Jovian gossamer rings) stochastics might well dominate the dynamics.
We construct a new RC phase shift network based Chua's circuit, which exhibits a period-doubling bifurcation route to chaos. Using coupled versions of such a phase-shift network based Chua's oscillators, we describe a new method for achieving complete synchronization (CS), approximate lag synchronization (LS), and approximate anticipating synchronization (AS) without delay or parameter mismatch. Employing the Pecora and Carroll approach, chaos synchronization is achieved in coupled chaotic oscillators, where the drive system variables control the response system. As a result, AS or LS or CS is demonstrated without using a variable delay line both experimentally and numerically.
Synchronization transitions in coupled time-delay electronic circuits with a threshold nonlinearity
(2011)
Experimental observations of typical kinds of synchronization transitions are reported in unidirectionally coupled time-delay electronic circuits with a threshold nonlinearity and two time delays, namely feedback delay tau(1) and coupling delay tau(2). We have observed transitions from anticipatory to lag via complete synchronization and their inverse counterparts with excitatory and inhibitory couplings, respectively, as a function of the coupling delay tau(2). The anticipating and lag times depend on the difference between the feedback and the coupling delays. A single stability condition for all the different types of synchronization is found to be valid as the stability condition is independent of both the delays. Further, the existence of different kinds of synchronizations observed experimentally is corroborated by numerical simulations and from the changes in the Lyapunov exponents of the coupled time-delay systems.
Correlations, as observed between the concentrations of metabolites in a biological sample, may be used to gain additional information about the physiological state of a given tissue. in this mini-review, we discuss the integration of these observed correlations into metabolomic networks and their relationships with the underlying biochemical pathways
Biochemical and genetic regulatory systems that involve low concentrations of molecules are inherently noisy. This intrinsic stochasticity, has received considerable interest recently, leading to new insights about the sources and consequences of noise in complex systems of genetic regulation. However, most prior work was devoted to the reduction of fluctuation and the robustness of cellular function with respect to intrinsic noise. Here, we focus on several scenarios in which the inherent molecular fluctuations are not merely a nuisance, but act constructively and bring about qualitative changes in the dynamics of the system. It will be demonstrated that in many typical situations biochemical and genetic regulatory systems may utilize intrinsic noise to their advantage. (C) 2002 Elsevier Ireland Ltd. All rights reserved
This paper employs a complex network approach to determine the topology and evolution of the network of extreme precipitation that governs the organization of extreme rainfall before, during, and after the Indian Summer Monsoon (ISM) season. We construct networks of extreme rainfall events during the ISM (June-September), post-monsoon (October-December), and pre-monsoon (March-May) periods from satellite-derived (Tropical Rainfall Measurement Mission, TRMM) and rain-gauge interpolated (Asian Precipitation Highly Resolved Observational Data Integration Towards the Evaluation of Water Resources, APHRODITE) data sets. The structure of the networks is determined by the level of synchronization of extreme rainfall events between different grid cells throughout the Indian subcontinent. Through the analysis of various complex-network metrics, we describe typical repetitive patterns in North Pakistan (NP), the Eastern Ghats (EG), and the Tibetan Plateau (TP). These patterns appear during the pre-monsoon season, evolve during the ISM, and disappear during the post-monsoon season. These are important meteorological features that need further attention and that may be useful in ISM timing and strength prediction.
Forecasting the onset and withdrawal of the Indian summer monsoon is crucial for the life and prosperity of more than one billion inhabitants of the Indian subcontinent. However, accurate prediction of monsoon timing remains a challenge, despite numerous efforts. Here we present a method for prediction of monsoon timing based on a critical transition precursor. We identify geographic regions-tipping elements of the monsoon-and use them as observation locations for predicting onset and withdrawal dates. Unlike most predictability methods, our approach does not rely on precipitation analysis but on air temperature and relative humidity, which are well represented both in models and observations. The proposed method allows to predict onset 2 weeks earlier and withdrawal dates 1.5 months earlier than existing methods. In addition, it enables to correctly forecast monsoon duration for some anomalous years, often associated with El Nino-Southern Oscillation.
We study the effects of parametric noise on a lattice network, which is locally modeled by a two-dimensional Rulkov map. We conclude that at some intermediate noise intensity, parametric noise can induce ordered circular patterns, which indicates the appearance of spatiotemporal coherence resonance in the studied lattice. With the observation of coherence-like manner in linear spatial cross-correlation, the coherence phenomena can be analyzed quantitatively.
We report the identification of global phase synchronization (GPS) in a linear array of unidirectionally coupled Mackey-Glass time-delay systems exhibiting highly non-phase-coherent chaotic attractors with complex topological structure. In particular, we show that the dynamical organization of all the coupled time-delay systems in the array to form GPS is achieved by sequential synchronization as a function of the coupling strength. Further, the asynchronous ones in the array with respect to the main sequentially synchronized cluster organize themselves to form clusters before they achieve synchronization with the main cluster. We have confirmed these results by estimating instantaneous phases including phase difference, average phase, average frequency, frequency ratio, and their differences from suitably transformed phase coherent attractors after using a nonlinear transformation of the original non-phase-coherent attractors. The results are further corroborated using two other independent approaches based on recurrence analysis and the concept of localized sets from the original non-phase-coherent attractors directly without explicitly introducing the measure of phase.
We study prebifurcation fluctuation amplification in nonlinear oscillators subject to bifurcations of spontaneous symmetry breaking which are manifest in the doubling of stable equilibrium states. Our theoretical estimates of both the linear growth and the nonlinear saturation of the fluctuations are in good agreement with our results from numerical simulations. We show that in the saturation mode, the fluctuation variance is proportional to the standard deviation of the external noise, whereas in the linear mode, the fluctuation variance is proportional to the noise variance. It is demonstrated that the phenomenon of prebifurcation noise amplification is more pronounced in the case of a slow transition through the bifurcation point. The amplification of fluctuations in this case makes it easier to form a symmetric probability of the final equilibrium states. In contrast, for a fast transition through the bifurcation point, the effect of amplification is much less pronounced. Under backward and forward passages through the bifurcation point, a loop of noise-dependent hysteresis emerges here. We find that for a fast transition of the nonlinear oscillator through the bifurcation point, the probability symmetry of the final equilibrium states is destroyed