### Refine

#### Year of publication

#### Document Type

- Article (36)
- Preprint (9)
- Monograph/Edited Volume (3)

#### Institute

- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (48) (remove)

The Voyager 2 Photopolarimeter experiment has yielded the highest resolved data of Saturn's rings, exhibiting a wide variety of features. The B-ring region between 105000 km and 110000 km distance from Saturn has been investigated. It has a high matter density and contains no significance features visible by eye. Analysis with statistical methods has let us to the detection of two significant events. These features are correlated with the inner 3:2 resonances of the F-ring shepherd satellites Pandora and Prometheus, and may be evidence of large ring paricles caught in the corotation resonances.

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.

The dynamics of noisy bistable systems is analyzed by means of Lyapunov exponents and measures of complexity. We consider both the classical Kramers problem with additive white noise and the case when the barrier fluctuates due to additional external colored noise. In case of additive noise we calculate the Lyapunov exponents and all measures of complexity analytically as functions of the noise intensity resp. the mean escape time. For the problem of fluctuating barrier the usual description of the dynamics with the mean escape time is not sufficient. The application of the concept of measures of complexity allows to describe the structures of motion in more detail. Most complexity measures sign the value of correlation time at which the phenomenon of resonant activation occurs with an extremum.

Contents: 1 Introduction 2 Experiment 3 Data 4 Symbolic dynamics 4.1 Symbolic dynamics as a tool for data analysis 4.2 2-symbols coding 4.3 3-symbols coding 5 Measures of complexity 5.1 Word statistics 5.2 Shannon entropy 6 Testing for stationarity 6.1 Stationarity 6.2 Time series of cycle durations 6.3 Chi-square test 7 Control parameters in the production of rhythms 8 Analysis of relative phases 9 Discussion 10 Outlook

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.

In the modern industrialized countries every year several hundred thousands of people die due to the sudden cardiac death. The individual risk for this sudden cardiac death cannot be defined precisely by common available, non-invasive diagnostic tools like Holter-monitoring, highly amplified ECG and traditional linear analysis of heart rate variability (HRV). Therefore, we apply some rather unconventional methods of nonlinear dynamics to analyse the HRV. Especially, some complexity measures that are basing on symbolic dynamics as well as a new measure, the renormalized entropy, detect some abnormalities in the HRV of several patients who have been classified in the low risk group by traditional methods. A combination of these complexity measures with the parameters in the frequency domain seems to be a promising way to get a more precise definition of the individual risk. These findings have to be validated by a representative number of patients.

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 have used techniques of nonlinear dynamics to compare a special model for the reversals of the Earth's magnetic field with the observational data. Although this model is rather simple, there is no essential difference to the data by means of well-known characteristics, such as correlation function and probability distribution. Applying methods of symbolic dynamics we have found that the considered model is not able to describe the dynamical properties of the observed process. These significant differences are expressed by algorithmic complexity and Renyi information.

Two deterministic processes leading to roughening interfaces are considered. It is shown that the dynamics of linear perturbations of turbulent regimes in coupled map lattices is governed by a discrete version of the Kardar-Parisi-Zhang equation. The asymptotic scaling behavior of the perturbation field is investigated in the case of large lattices. Secondly, the dynamics of an order-disorder interface is modelled with a simple two-dimensional coupled map lattice, possesing a turbulent and a laminar state. It is demonstrated, that in some range of parameters the spreading of the turbulent state is accompanied by kinetic roughening of the interface.

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.

Objectives: Scoring sleep visually based on polysomnography is an important but time-consuming element of sleep medicine. Where-as computer software assists human experts in the assignment of sleep stages to polysomnogram epochs, their performance is usually insufficient. This study evaluates the possibility to fully automatize sleep staging considering the reliability of the sleep stages available from human expert sleep scorers. Methods: We obtain features from EEG, ECG and respiratory signals of polysomnograms from ten healthy subjects. Using the sleep stages provided by three human experts, we evaluate the performance of linear discriminant analysis on the entire polysomnogram and:only on epochs where the three experts agree in their-sleep stage scoring. Results: We show that in polysomnogram intervals, to which all three scorers assign the same sleep stage, our algorithm achieves 90% accuracy. This high rate of agreement with the human experts is accomplished with only a small set of three frequency features from the EEG. We increase-the performance to 93% by including ECG and respiration features. In contrast, on intervals of ambiguous sleep stage, the sleep stage classification obtained from our algorithm, agrees with the human consensus scorer in approximately 61%. Conclusions: These findings suggest that machine classification is highly consistent with human sleep staging and that error in the algorithm's assignments is rather a problem of lack of well-defined criteria for human experts to judge certain polysomnogram epochs than an insufficiency of computational procedures

The incidence of cardiovascular diseases increases with the growth of the human population and an aging society, leading to very high expenses in the public health system. Therefore, it is challenging to develop sophisticated methods in order to improve medical diagnostics. The question whether the normal heart rate is chaotic or not is an attempt to elucidate the underlying mechanisms of cardiovascular dynamics and therefore a highly controversial topical challenge. In this contribution we demonstrate that linear and nonlinear parameters allow us to separate completely the data sets of the three groups provided for this controversial topic in nonlinear dynamics. The question whether these time series are chaotic or not cannot be answered satisfactorily without investigating the underlying mechanisms leading to them. We give an example of the dominant influence of respiration on heart beat dynamics, which shows that observed fluctuations can be mostly explained by respiratory modulations of heart rate and blood pressure (coefficient of determination: 96%). Therefore, we recommend reformulating the following initial question: "Is the normal heart rate chaotic?" We rather ask the following: " Is the normal heart rate 'chaotic' due to respiration?"

Aims. Sunspot distribution in the northern and southern solar hemispheres exibit striking synchronous behaviour on the scale of a Schwabe cycle. However, sometimes the bilateral symmetry of the Butterfly diagram relative to the solar equatorial plane breaks down. The investigation of this phenomenon is important to explaining the almost-periodic behaviour of solar cycles. Methods. We use cross-recurrence plots for the study of the time-varying phase asymmetry of the northern and southern hemisphere and compare our results with the latitudinal distribution of the sunspots. Results. We observe a long-term persistence of phase leading in one of the hemispheres, which lasts almost 4 solar cycles and probably corresponds to the Gleissberg cycle. Long-term variations in the hemispheric-leading do not demonstrate clear periodicity but are strongly anti-correlated with the long-term variations in the magnetic equator.

Sensory information entering the nervous system follows independent paths of processing such that specific features are individually detected. However, sensory perception, awareness, and cognition emerge from the combination of information. Here we have analyzed the corticocortical network of the cat, looking for the anatomical substrate which permits the simultaneous segregation and integration of information in the brain. We find that cortical communications are mainly governed by three topological factors of the underlying network: (i) a large density of connections, (ii) segregation of cortical areas into clusters, and (iii) the presence of highly connected hubs aiding the multisensory processing and integration. Statistical analysis of the shortest paths reveals that, while information is highly accessible to all cortical areas, the complexity of cortical information processing may arise from the rich and intricate alternative paths in which areas can influence each other.

The response of scale-free networks with community structure to external stimuli is studied. By disturbing some nodes with different strategies, it is shown that the robustness of this kind of network can be enhanced due to the existence of communities in the networks. Some of the response patterns are found to coincide with topological communities. We show that such phenomena also occur in the cat brain network which is an example of a scale-free like network with community structure. Our results provide insights into the relationship between network topology and the functional organization in complex networks from another viewpoint.

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.

Experimental evidences point Out the participation of nonsynaptic mechanisms (e.g., fluctuations in extracellular tons) in epileptiform bursting and spreading depression (SD). During these abnormal oscillatory patterns, it is observed an increase of extracellular potassium concentration [K+](o) and a decrease of extracellular calcium concentration [Ca2+](o) which raises the neuronal excitability. However, whether the high [K+](o) triggers and propagates these abnormal neuronal activities or plays a secondary role into this process is unclear. To better understand the influence of extracellular potassium dynamics in these oscillatory patterns, the experimental conditions of high [K+](o) and zero [Ca2+](o) were replicated in an extended Golomb model where we added important regulatory mechanisms of ion concentration as Na+-K+ pump, ion diffusion and glial buffering. Within these Conditions, simulations of the cell model exhibit seizure-like discharges (ictal bursting). The SD was elicited by the interruption of the Na+- K+ pump activity, mimicking the effect of cellular hypoxia (an experimental protocol to elicit SD, the hypoxia-induced SD). We used the bifurcation theory and the fast-slow method to analyze the interference of K+ dynamics in the cellular excitability. This analysis indicates that the system loses its stability at a high [K+](o), transiting to an elevated state of neuronal excitability. Effects of high [K+](o), are observed in different stages of ictal bursting and SD. In the initial stage, the increase of [K+](o) creates favorable conditions to trigger both oscillatory patterns. During the neuronal activity, a continuous growth of [K+](o) by outward K+ flow depresses K+ Currents in a positive feedback way. At the last stage, due to the depression of K+ currents, the Na+-K+ pump is the main mechanism in the end of neuronal activity. Thus, this work suggests that [K+](o) dynamics may play a fundamental role in these abnormal oscillatory 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.

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.