Institut für Physik und Astronomie
Refine
Year of publication
Document Type
- Article (212)
- Preprint (9)
- Postprint (8)
- Monograph/Edited Volume (4)
- Other (1)
Keywords
- Amazon rainforest (2)
- Complex networks (2)
- channel (2)
- complex networks (2)
- diffusion (2)
- droughts (2)
- prediction (2)
- space-dependent diffusivity (2)
- 3D medical image analysis (1)
- African climate (1)
Institute
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.
Dynamical organization of connection weights is studied in scale-free networks of chaotic oscillators, where the coupling strength of a node from its neighbors develops adaptively according to the local synchronization property between the node and its neighbors. We find that when complete synchronization is achieved, the coupling strength becomes weighted and correlated with the topology due to a hierarchical transition to synchronization in heterogeneous networks. Importantly, such an adaptive process enhances significantly the synchronizability of the networks, which could have meaningful implications in the manipulation of dynamical networks
We introduce a modified dynamical optimization coupling scheme to enhance the synchronizability in the scale- free networks as well as to keep uniform and converging intensities during the transition to synchronization. Further, the size of networks that can be synchronizable exceeds by several orders of magnitude the size of unweighted networks.
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
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.
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 show that two dynamical invariants, the second order Renyi entropy and the correlation dimension, can be estimated from recurrence plots (RPs) with arbitrary embedding dimension and delay. This fact is interesting as these quantities are even invariant if no embedding is used. This is an important advantage of RPs compared to other techniques of nonlinear data analysis. These estimates for the correlation dimension and entropy are robust and, moreover, can be obtained at a low numerical cost. We exemplify our results for the Rossler system, the funnel attractor and the Mackey-Glass system. In the last part of the paper we estimate dynamical invariants for data from some fluid dynamical experiments and confirm previous evidence for low dimensional chaos in this experimental system. (C) 2004 American Institute of Physics
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.
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.
Standard time and frequency parameters of heart rate variability (HRV) describe only linear and periodic behaviour, whereas more complex relationships cannot be recognised. A method that may be capable of assessing more complex properties is the non-linear measure of 'renormalised entropy.' A new concept of the method, RE(AR), has been developed, based on a non-linear renormalisation of autoregressive spectral distributions. To test the hypothesis that renormalised entropy may improve the result of high-risk stratification after myocardial infarction, it is applied to a clinical pilot study (41 subjects) and to prospective data of the St George's Hospital post- infarction database (572 patients). The study shows that the new RE(AR) method is more reproducible and more stable in time than a previously introduced method (p<0.001). Moreover, the results of the study confirm the hypothesis that on average, the survivors have negative values of RE(AR) (-0.11+/-0.18), whereas the non-survivors have positive values (0.03+/-0.22, p<0.01). Further, the study shows that the combination of an HRV triangular index and RE(AR) leads to a better prediction of sudden arrhythmic death than standard measurements of HRV. In summary, the new RE(AR) method is an independent measure in HRV analysis that may be suitable for risk stratification in patients after myocardial infarction.
Interacting human activities underlie the patterns of many social, technological, and economic phenomena. Here we present clear empirical evidence from Short Message correspondence that observed human actions are the result of the interplay of three basic ingredients: Poisson initiation of tasks and decision making for task execution in individual humans as well as interaction among individuals. This interplay leads to new types of interevent time distribution, neither completely Poisson nor power-law, but a bimodal combination of them. We show that the events can be separated into independent bursts which are generated by frequent mutual interactions in short times following random initiations of communications in longer times by the individuals. We introduce a minimal model of two interacting priority queues incorporating the three basic ingredients which fits well the distributions using the parameters extracted from the empirical data. The model can also embrace a range of realistic social interacting systems such as e-mail and letter communications when taking the time scale of processing into account. Our findings provide insight into various human activities both at the individual and network level. Our analysis and modeling of bimodal activity in human communication from the viewpoint of the interplay between processes of different time scales is likely to shed light on bimodal phenomena in other complex systems, such as interevent times in earthquakes, rainfall, forest fire, and economic systems, etc.
Experimental evidence of anomalous phase synchronization in two diffusively coupled Chua oscillators
(2006)
We study the transition to phase synchronization in two diffusively coupled, nonidentical Chua oscillators. In the experiments, depending on the used parameterization, we observe several distinct routes to phase synchronization, including states of either in-phase, out-of-phase, or antiphase synchronization, which may be intersected by an intermediate desynchronization regime with large fluctuations of the frequency difference. Furthermore, we report the first experimental evidence of an anomalous transition to phase synchronization, which is characterized by an initial enlargement of the natural frequency difference with coupling strength. This results in a maximal frequency disorder at intermediate coupling levels, whereas usual phase synchronization via monotonic decrease in frequency difference sets in only for larger coupling values. All experimental results are supported by numerical simulations of two coupled Chua models
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
We give evidence of frequency entrainment of dominant peaks in the chaotic spectra of two coupled chaotic nonautonomous oscillators. At variance with the autonomous case, the phenomenon is here characterized by the vanishing of a previously positive Lyapunov exponent in the spectrum, which takes place for a broad range of the coupling strength parameter. Such a state is studied also for the case of chaotic oscillators with ill-defined phases due to the absence of a unique center of rotation. Different phase synchronization indicators are used to circumvent this difficulty
We study frequency selectivity in noise-induced subthreshold signal processing in a system with many noise- supported stochastic attractors which are created due to slow variable diffusion between identical excitable elements. Such a coupling provides coexisting of several average periods distinct from that of an isolated oscillator and several phase relations between elements. We show that the response of the coupled elements under different noise levels can be significantly enhanced or reduced by forcing some elements in resonance with these new frequencies which correspond to appropriate phase relations
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.
In this paper we present an approach to recover the dynamics from recurrences of a system and then generate (multivariate) twin surrogate (TS) trajectories. In contrast to other approaches, such as the linear-like surrogates, this technique produces surrogates which correspond to an independent copy of the underlying system, i. e. they induce a trajectory of the underlying system visiting the attractor in a different way. We show that these surrogates are well suited to test for complex synchronization, which makes it possible to systematically assess the reliability of synchronization analyses. We then apply the TS to study binocular fixational movements and find strong indications that the fixational movements of the left and right eye are phase synchronized. This result indicates that there might be one centre only in the brain that produces the fixational movements in both eyes or a close link between two centres.
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.
How do diverse dynamical patterns arise from the topology of complex networks? We study synchronization dynamics in the cortical brain network of the cat, which displays a hierarchically clustered organization, by modeling each node (cortical area) with a subnetwork of interacting excitable neurons. We find that in the biologically plausible regime the dynamics exhibits a hierarchical modular organization, in particular, revealing functional clusters coinciding with the anatomical communities at different scales. Our results provide insights into the relationship between network topology and functional organization of complex brain networks.
We study synchronization behavior in networks of coupled chaotic oscillators with heterogeneous connection degrees. Our focus is on regimes away from the complete synchronization state, when the coupling is not strong enough, when the oscillators are under the influence of noise or when the oscillators are nonidentical. We have found a hierarchical organization of the synchronization behavior with respect to the collective dynamics of the network. Oscillators with more connections (hubs) are synchronized more closely by the collective dynamics and constitute the dynamical core of the network. The numerical observation of this hierarchical synchronization is supported with an analysis based on a mean field approximation and the master stability function. (C) 2006 American Institute of Physics
Higher resilience to climatic disturbances in tropical vegetation exposed to more variable rainfall
(2019)
With ongoing global warming, the amount and frequency of precipitation in the tropics is projected to change substantially. While it has been shown that tropical forests and savannahs are sustained within the same intermediate mean annual precipitation range, the mechanisms that lead to the resilience of these ecosystems are still not fully understood. In particular, the long-term impact of rainfall variability on resilience is as yet unclear. Here we present observational evidence that both tropical forest and savannah exposed to a higher rainfall variability-in particular on interannual scales-during their long-term past are overall more resilient against climatic disturbances. Based on precipitation and tree cover data in the Brazilian Amazon basin, we constructed potential landscapes that enable us to systematically measure the resilience of the different ecosystems. Additionally, we infer that shifts from forest to savannah due to decreasing precipitation in the future are more likely to occur in regions with a precursory lower rainfall variability. Long-term rainfall variability thus needs to be taken into account in resilience analyses and projections of vegetation response to climate change.
We investigate a network of influences connected to global mean temperature. Considering various climatic factors known to influence global mean temperature, we evaluate not only the impacts of these factors on temperature but also the directed dependencies among the factors themselves. Based on an existing recurrence-based connectivity measure, we propose a new and more general measure that quantifies the level of dependence between two time series based on joint recurrences at a chosen time delay. The measures estimated in the analysis are tested for statistical significance using twin surrogates. We find, in accordance with earlier studies, the major drivers for global mean temperature to be greenhouse gases, ENSO, volcanic activity, and solar irradiance. We further uncover a feedback between temperature and ENSO. Our results demonstrate the need to involve multiple, delayed interactions within the drivers of temperature in order to develop a more thorough picture of global temperature variations.
Recurrence plots have recently been recognized as a powerful tool for the analysis of data. Not only the visualization of structures of the time series but also the possibility to estimate invariants from them and the possibility to analyze non-stationary data sets are remarkable. However, the question of how much information is encoded in such a two-dimensional and binary representation has not been discussed so far. In this Letter we show that-under some conditions-it is possible to reconstruct an attractor from the recurrence plot, at least topologically. This means that all relevant dynamical information is contained in the plot. (C) 2004 Elsevier B.V. All rights reserved
Human comment is studied using data from 'tianya' which is one of the most popular on-line social systems in China. We found that the time interval between two consecutive comments on the same topic, called inter-event time, follows a power-law distribution. This result shows that there is no characteristic decay time on a topic. It allows for very long periods without comments that separate bursts of intensive comments. Furthermore, the frequency of a different ID commenting on a topic also follows a power-law distribution. It indicates that there are some "hubs" in the topic who lead the direction of the public opinion. Based on the personal comments habit, a model is introduced to explain these phenomena. The numerical simulations of the model fit well with the empirical results. Our findings are helpful for discovering regular patterns of human behavior in on-line society and the evolution of the public opinion on the virtual as well as real society.
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.
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 show many versatile phase synchronous configurations that emerge in an array of coupled chaotic elements due to the presence of a periodic stimulus. Then, we explain the relevance of these configurations to the understanding of how information about such a. stimulus is transmitted from one side to the other in this array. The stimulus actively creates the ways to be transmitted, by making the chaotic elements to phase synchronize
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
We study Hamiltonian chaos generated by the dynamics of passive tracers moving in a two-dimensional fluid flow and describe the complex structure formed in a chaotic layer that separates a vortex region from the shear flow. The stable and unstable manifolds of unstable periodic orbits are computed. It is shown that their intersections in the Poincare map as an invariant set of homoclinic points constitute the backbone of the chaotic layer. Special attention is paid to the finite time properties of the chaotic layer. In particular, finite time Lyapunov exponents are computed and a scaling law of the variance of their distribution is derived. Additionally, the box counting dimension as an effective dimension to characterize the fractal properties of the layer is estimated for different duration times of simulation. Its behavior in the asymptotic time limit is discussed. By computing the Lyapunov exponents and by applying methods of symbolic dynamics, the formation of the layer as a function of the external forcing strength, which in turn represents the perturbation of the originally integrable system, is characterized. In particular, it is shown that the capture of KAM tori by the layer has a remarkable influence on the averaged Lyapunov exponents. (C) 2004 Elsevier Ltd. All rights reserved
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
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
Locking-based frequency measurement and synchronization of chaotic oscillators with complex dynamics
(2002)
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.
We study several algorithms to simulate bone mass loss in two-dimensional and three-dimensional computed tomography bone images. The aim is to extrapolate and predict the bone loss, to provide test objects for newly developed structural measures, and to understand the physical mechanisms behind the bone alteration. Our bone model approach differs from those already reported in the literature by two features. First, we work with original bone images, obtained by computed tomography (CT); second, we use structural measures of complexity to evaluate bone resorption and to compare it with the data provided by CT. This gives us the possibility to test algorithms of bone resorption by comparing their results with experimentally found dependencies of structural measures of complexity, as well as to show efficiency of the complexity measures in the analysis of bone models. For two-dimensional images we suggest two algorithms, a threshold algorithm and a virtual slicing algorithm. The threshold algorithm simulates bone resorption on a boundary between bone and marrow, representing an activity of osteoclasts. The virtual slicing algorithm uses a distribution of the bone material between several virtually created slices to achieve statistically correct results, when the bone-marrow transition is not clearly defined. These algorithms have been tested for original CT 10 mm thick vertebral slices and for simulated 10 mm thick slices constructed from ten I mm thick slices. For three-dimensional data, we suggest a variation of the threshold algorithm and apply it to bone images. The results of modeling have been compared with CT images using structural measures of complexity in two- and three-dimensions. This comparison has confirmed credibility of a virtual slicing modeling algorithm for two-dimensional data and a threshold algorithm for three-dimensional data