Refine
Has Fulltext
- no (289) (remove)
Year of publication
Document Type
- Article (280)
- Monograph/Edited Volume (8)
- Other (1)
Language
- English (289) (remove)
Is part of the Bibliography
- yes (289)
Keywords
- Complex networks (5)
- Event synchronization (4)
- Extreme rainfall (2)
- Synchronization (2)
- synchronization (2)
- 3D medical image analysis (1)
- African climate (1)
- Algebraic geometry (1)
- Amazon rainforest (1)
- Anisotropy (1)
Institute
- Institut für Physik und Astronomie (215)
- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (38)
- Institut für Geowissenschaften (20)
- Department Psychologie (12)
- Institut für Biochemie und Biologie (4)
- Department Linguistik (2)
- Institut für Informatik und Computational Science (2)
- Department Sport- und Gesundheitswissenschaften (1)
- Institut für Umweltwissenschaften und Geographie (1)
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.
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.
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.
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
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.
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.
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
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.
Locking-based frequency measurement and synchronization of chaotic oscillators with complex dynamics
(2002)
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.
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.
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.
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
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 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 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.
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 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.
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.
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
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.
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 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.
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 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.
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 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.
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.
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.
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.
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.
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.
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
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.
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
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.
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.
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.
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
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.
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.
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.
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.
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.
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.
We use the extension of the method of recurrence plots to cross recurrence plots (CRP) which enables a nonlinear analysis of bivariate data. To quantify CRPs, we develop further three measures of complexity mainly basing on diagonal structures in CRPs. The CRP analysis of prototypical model systems with nonlinear interactions demonstrates that this technique enables to find these nonlinear interrelations from bivariate time series, whereas linear correlation tests do not. Applying the CRP analysis to climatological data, we find a complex relationship between rainfall and El Nino data.
In the recent article "Stochastic analysis of recurrence plots with applications to the detection of deterministic signals" (Physica D 237 (2008) 619-629), Rohde et al. stated that the performance of RQA in order to detect deterministic signals would be below traditional and well-known detectors. However, we have concerns about such a general statement. Based on our own studies we cannot confirm their conclusions. Our findings suggest that the measures of complexity provided by RQA are useful detectors outperforming well-known traditional detectors, in particular for the detection of signals of complex systems, with phase differences or signals modified due to the measurement process.
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
We present a modern method used in nonlinear time series analysis to investigate the relation of two oscillating systems with respect to their phases, independently of their amplitudes. We study the difference of the phase dynamics between El Nino/Southern Oscillation (ENSO) and the Indian Monsoon on inter-annual time scales. We identify distinct epochs, especially two intervals of phase coherence, 1886 - 1908 and 1964 - 1980, corroborating earlier findings from a new point of view. A significance test shows that the coherence is very unlikely to be the result of stochastic fluctuations. We also detect so far unknown periods of coupling which are invisible to linear methods. These findings suggest that the decreasing correlation during the last decades might be a typical epoch of the ENSO/ Monsoon system having occurred repeatedly. The high time resolution of the method enables us to present an interpretation of how volcanic radiative forcing could cause the coupling
In this paper, we present a detailed evaluation of cross wavelet analysis of bivariate time series. We develop a statistical test for zero wavelet coherency based on Monte Carlo simulations. If at least one of the two processes considered is Gaussian white noise, an approximative formula for the critical value can be utilized. In a second part, typical pitfalls of wavelet cross spectra and wavelet coherency are discussed. The wavelet cross spectrum appears to be not suitable for significance testing the interrelation between two processes. Instead, one should rather apply wavelet coherency. Furthermore we investigate problems due to multiple testing. Based on these results, we show that coherency between ENSO and NAO is an artefact for most of the time from 1900 to 1995. However, during a distinct period from around 1920 to 1940, significant coherency between the two phenomena occurs
We propose a novel approach based on the fluctuation of similarity to identify regimes of distinct dynamical complexity in short time series. A statistical test is developed to estimate the significance of the identified transitions. Our method is verified by uncovering bifurcation structures in several paradigmatic models, providing more complex transitions compared with traditional Lyapunov exponents. In a real-world situation, we apply this method to identify millennial-scale dynamical transitions in Plio-Pleistocene proxy records of the South Asian summer monsoon system. We infer that many of these transitions are induced by the external forcing of the solar insolation and are also affected by internal forcing on Monsoonal dynamics, i.e., the glaciation cycles of the Northern Hemisphere and the onset of the Walker circulation.
Analysis of spatial and temporal extreme monsoonal rainfall over South Asia using complex networks
(2012)
We present a detailed analysis of summer monsoon rainfall over the Indian peninsular using nonlinear spatial correlations. This analysis is carried out employing the tools of complex networks and a measure of nonlinear correlation for point processes such as rainfall, called event synchronization. This study provides valuable insights into the spatial organization, scales, and structure of the 90th and 94th percentile rainfall events during the Indian summer monsoon (June-September). We furthermore analyse the influence of different critical synoptic atmospheric systems and the impact of the steep Himalayan topography on rainfall patterns. The presented method not only helps us in visualising the structure of the extreme-event rainfall fields, but also identifies the water vapor pathways and decadal-scale moisture sinks over the region. Furthermore a simple scheme based on complex networks is presented to decipher the spatial intricacies and temporal evolution of monsoonal rainfall patterns over the last 6 decades.
Precipitation patterns and extremes are significantly influenced by various climatic factors and large-scale atmospheric circulation patterns. This study uses wavelet coherence analysis to detect significant interannual and interdecadal oscillations in monthly precipitation extremes across India and their teleconnections to three prominent climate indices, namely, Nino 3.4, Pacific Decadal Oscillation, and Indian Ocean Dipole (IOD). Further, partial wavelet coherence analysis is used to estimate the standalone relationship between the climate indices and precipitation after removing the effect of interdependency. The wavelet analysis of monthly precipitation extremes at 30 different locations across India reveals that (a) interannual (2-8 years) and interdecadal (8-32 years) oscillations are statistically significant, and (b) the oscillations vary in both time and space. The results from the partial wavelet coherence analysis reveal that Nino 3.4 and IOD are the significant drivers of Indian precipitation at interannual and interdecadal scales. Intriguingly, the study also confirms that the strength of influence of large-scale atmospheric circulation patterns on Indian precipitation extremes varies with spatial physiography of the region.
We present a study of ocean convection parameterization based on a novel approach which includes both eddy diffusion and advection and consists of a two-dimensional lattice of bistable maps. This approach retains important features of usual grid models and allows to assess the relative roles of diffusion and advection in the spreading of convective cells. For large diffusion our model exhibits a phase transition from convective patterns to a homogeneous state over the entire lattice. In hysteresis experiments we find staircase behavior depending on stability thresholds of local convection patterns. This nonphysical behavior is suspected to induce spurious abrupt changes in the spreading of convection in ocean models. The final steady state of convective cells depends not only on the magnitude of the advective velocity but also on its direction, implying a possible bias in the development of convective patterns. Such bias points to the need for an appropriate choice of grid geometry in ocean modeling
This work focuses on the dynamics of particles in a confined geometry with position-dependent diffusivity, where the confinement is modelled by a periodic channel consisting of unit cells connected by narrow passage ways. We consider three functional forms for the diffusivity, corresponding to the scenarios of a constant (D ₀), as well as a low (D ₘ) and a high (D d) mobility diffusion in cell centre of the longitudinally symmetric cells. Due to the interaction among the diffusivity, channel shape and external force, the system exhibits complex and interesting phenomena. By calculating the probability density function, mean velocity and mean first exit time with the Itô calculus form, we find that in the absence of external forces the diffusivity D d will redistribute particles near the channel wall, while the diffusivity D ₘ will trap them near the cell centre. The superposition of external forces will break their static distributions. Besides, our results demonstrate that for the diffusivity D d, a high dependence on the x coordinate (parallel with the central channel line) will improve the mean velocity of the particles. In contrast, for the diffusivity D ₘ, a weak dependence on the x coordinate will dramatically accelerate the moving speed. In addition, it shows that a large external force can weaken the influences of different diffusivities; inversely, for a small external force, the types of diffusivity affect significantly the particle dynamics. In practice, one can apply these results to achieve a prominent enhancement of the particle transport in two- or three-dimensional channels by modulating the local tracer diffusivity via an engineered gel of varying porosity or by adding a cold tube to cool down the diffusivity along the central line, which may be a relevant effect in engineering applications. Effects of different stochastic calculi in the evaluation of the underlying multiplicative stochastic equation for different physical scenarios are discussed.
We present conditions for the local and global synchronizations in coupled-map networks using the matrix measure approach. In contrast to many existing synchronization conditions, the proposed synchronization criteria do not depend on the solution of the synchronous state and give less limitation on the network connections. Numerical simulations of the coupled quadratic maps demonstrate the potentials of our main results.
In the present paper, two kinds of dynamical complex networks are considered. The first is that elements of every node have different time delays but all nodes in Such networks have the same time-delay vector. The second is that different nodes have different time-delay vectors, and the elements of each node also have different time delays. Corresponding synchronization theorems are established. Numerical examples show the efficiency of the derived theorems.
The effect of additive noise on transitions in nonlinear systems far from equilibrium is studied. It is shown that additive noise in itself can induce a hidden phase transition, which is similar to the transition induced by multiplicative noise in a nonlinear oscillator [P. Landa and A. Zaikin, Phys. Rev. E 54, 3535 (1996)]. Investigation of different nonlinear models that demonstrate phase transitions induced by multiplicative noise shows that the influence of additive noise upon such phase transitions can be crucial: additive noise can either blur such a transition or stabilize noise-induced oscillations.
The results of the theoretical consideration of stochastic resonance in overdamped bistable oscillators are given. These results are founded not on the model of two states as in [McNamara B, Wiesenfeld K. Theory of stochastic resonance. Phys Rev A 1989;39:4854-69], but on splitting of motion into regular and random and the rigorous solution of the Fokker-Planck equation for the random component. We show that this resonance is caused by a change, under the influence of noise, of the system's effective stiffness and damping factor contained in the equation for the regular component. For a certain value of the noise intensity the effective stiffness is minimal, and this fact causes non-monotonic change of the output signal amplitude as the noise intensity changes. It is important that the location of the minimum and its value depend essentially on the signal frequency.
We have discussed some tools from nonlinear dynamics which may help to analyze transient phenomena, such as solar bursts. The structure function known from turbulence theory is an appropriate method to find out some scaling behavior of fluctuations in time. More generally, the wavelet analysis, which is some generalization of the power spectrum, exhibits information on the location as well as the size of hidden characteristic features. Applying both techniques to microwave bursts, we have found some scaling properties that refer to the existence of hierarchic time structures. This is in good accordance with the electric circuit model for describing the flare-particle energization process.
The radiocarbon record that has been extended from 7199 BC to 1891 AD is of fundamental importance to understand century-scale variations of solar activity. We have, therefore, studied how to extract information from dynamic reconstructions of this observational record. Using some rather unusual methods of nonlinear dynamics, we have found that the data are significantly different from linear colored noise and that there is some evidence of nonlinear behavior. The method of recurrence plots exhibits that the grand minima of solar activity are quite different in their recurrence. Most remarkably, it suggests that the recent epoch seems to be similar to the Medieval maximum.
We look for structural properties in the light curve of the dwarf nova SS Cyg by means of techniques from nonlinear dynamics. Applying the popular Grassberger-Procaccia procedure, Cannizzo and Goddings (1988) showed that there is no evidence for a low-dimensional attractor underlying this record. Because there are some hints for order in the light curve, we search for other signatures of deterministic systems. Therefore, we use other methods recently developed in this theory, such as local linear prediction and recurrence maps. Our main findings are: i] the prediction error grows exponentially during outburst phases, but via a power law in the quiescent states, ii] there are some rather regular patterns in this light curve which sometimes recur, but the recurrence is not regular. This leads to the following conclusions: i] The outburst dynamics shows a higher degree of order than the quiescent one. There are some hints for deterministic chaos in the outburst behavior. ii] The light curve is a complex mixture of deterministic and stochastic structures. The analysis presented in this paper shows that methods of nonlinear dynamics can be an efficient tool for the study of complex processes, even if there is no evidence for a low-dimensional attractor.
The application of chaos theory has become popular to understand the nature of various features of solar activity because most of them are far from regular. The usual approach, however, that is basing on finding low- dimensional structures of the underlying processes seems to be successful only in a few exceptional cases, such as in rather coherent phenomena as coronal pulsations. It is important to note that most phenomena in solar radio emission are more complex. We present two kinds of techniques from nonlinear dynamics which can be useful to analyse such phenomena: i] Fragmentation processes observed in solar spike events are studied by means of symbolic dynamics methods. Different measures of complexity calculated from such observations reveal that there is some order in this fragmentation. ii] Bursts are a typical transient phenomenon. To study energization processes causing impulsive microwave bursts, the wavelet analysis is applied. It exhibits structural differences of the pre- and post-impulsive phase in cases where the power spectra of both are not distinct.
We present two different approaches to detect and quantify phase synchronization in the case of coupled non- phase coherent oscillators. The first one is based on the general idea of curvature of an arbitrary curve. The second one is based on recurrences of the trajectory in phase space. We illustrate both methods in the paradigmatic example of the Rossler system in the funnel regime. We show that the second method is applicable even in the case of noisy data. Furthermore, we extend the second approach to the application of chains of coupled systems, which allows us to detect easily clusters of synchronized oscillators. In order to illustrate the applicability of this approach, we show the results of the algorithm applied to experimental data from a population of 64 electrochemical oscillators