Refine
Year of publication
Document Type
- Preprint (50)
- Article (43)
- Monograph/Edited Volume (19)
- Doctoral Thesis (1)
Language
- English (113) (remove)
Keywords
- aerosol size distribution (2)
- inversion (2)
- Ill-posed problem (1)
- MHD-equations (1)
- Magnetfeld-Satellit (1)
- Magnetic field measurements (1)
- Magnetische Feldmessungen (1)
- Magnetometer-Kalibrierung (1)
- Multiwavelength LIDAR (1)
- Planetary Rings (1)
Institute
- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (113) (remove)
The modular structure is pervasive in many complex networks of interactions observed in natural, social and technological sciences. Its study sheds light on the relation between the structure and the function of complex systems. Generally speaking, modules are islands of highly connected nodes separated by a relatively small number of links. Every module can have the contributions of links from any node in the network. The challenge is to disentangle these contributions to understand how the modular structure is built. The main problem is that the analysis of a certain partition into modules involves, in principle, as much data as the number of modules times the number of nodes. To confront this challenge, here we first define the contribution matrix, the mathematical object containing all the information about the partition of interest, and then we use truncated singular value decomposition to extract the best representation of this matrix in a plane. The analysis of this projection allows us to scrutinize the skeleton of the modular structure, revealing the structure of individual modules and their interrelations.
During the last glacial period, climate records from the North Atlantic region exhibit a pronounced spectral component corresponding to a period of about 1470 years, which has attracted much attention. This spectral peak is closely related to the recurrence pattern of Dansgaard-Oeschger (DO) events. In previous studies a red noise random process, more precisely a first-order autoregressive (AR1) process, was used to evaluate the statistical significance of this peak, with a reported significance of more than 99%. Here we use a simple mechanistic two-state model of DO events, which itself was derived from a much more sophisticated ocean-atmosphere model of intermediate complexity, to numerically evaluate the spectral properties of random (i.e., solely noise-driven) events. This way we find that the power spectral density of random DO events differs fundamentally from a simple red noise random process. These results question the applicability of linear spectral analysis for estimating the statistical significance of highly non-linear processes such as DO events. More precisely, to enhance our scientific understanding about the trigger of DO events, we must not consider simple "straw men" as, for example, the AR1 random process, but rather test against realistic alternative descriptions.
We have shown that the two-dimensional complex Ginzburg-Landau equation exhibits supertransient chaos in a certain parameter range. Using numerical methods this behavior is found near the transition line separating frozen spiral solutions from turbulence. Supertransient chaos seems to be a common phenomenon in extended spatiotemporal systems. These supertransients are characterized by an average transient lifetime which depends exponentially on the size of the system and are due to an underlying nonattracting chaotic set.
We have numerically studied the bifurcations and transition to chaos in a two-dimensional fluid for varying values of the Reynolds number. These investigations have been motivated by experiments in fluids, where an array of vortices was driven by an electromotive force. In these experiments, successive changes leading to a complex motion of the vortices, due to increased forcing, have been explored [Tabeling, Perrin, and Fauve, J. Fluid Mech. 213, 511 (1990)]. We model this experiment by means of two-dimensional Navier-Stokes equations with a special external forcing, driving a linear chain of eight counter-rotating vortices, imposing stress-free boundary conditions in the vertical direction and periodic boundary conditions in the horizontal direction. As the strength of the forcing or the Reynolds number is raised, the original stationary vortex array becomes unstable and a complex sequence of bifurcations is observed. Several steady states and periodic branches and a period doubling cascade appear on the route to chaos. For increasing values of the Reynolds number, shear flow develops, for which the spatial scale is large compared to the scale of the forcing. Furthermore, we have investigated the influence of the aspect ratio of the container as well as the effect of no-slip boundary conditions at the top and bottom, on the bifurcation scenario.
We have studied the bifurcation structure of the incompressible two-dimensional Navier-Stokes equations with a special external forcing driving an array of 8×8 counterrotating vortices. The study has been motivated by recent experiments with thin layers of electrolytes showing, among other things, the formation of large-scale spatial patterns. As the strength of the forcing or the Reynolds number is raised the original stationary vortex array becomes unstable and a complex sequence of bifurcations is observed. The bifurcations lead to several periodic branches, torus and chaotic solutions, and other stationary solutions. Most remarkable is the appearance of solutions characterized by structures on spatial scales large compared to the scale of the forcing. We also characterize the different dynamic regimes by means of tracers injected into the fluid. Stretching rates and Hausdorff dimensions of convected line elements are calculated to quantify the mixing process. It turns out that for time-periodic velocity fields the mixing can be very effective.
The ill-posed problem of aerosol size distribution determination from a small number of backscatter and extinction measurements was solved successfully with a mollifier method which is advantageous since the ill-posed part is performed on exactly given quantities, the points r where n(r) is evaluated may be freely selected. A new twodimensional model for the troposphere is proposed.
The determination of the atmospheric aerosol size distribution is an inverse illposed problem. The shape and the material composition of the air-carried particles are two substantial model parameters. Present evaluation algorithms only used an approximation with spherical homogeneous particles. In this paper we propose a new numerically efficient recursive algorithm for inhomogeneous multilayered coated and absorbing particles. Numerical results of real existing particles show that the influence of the two parameters on the model is very important and therefore cannot be ignored.
The ill-posed inversion of multiwavelength lidar data by a hybrid method of variable projection
(1999)
The ill-posed problem of aerosol distribution determination from a small number of backscatter and extinction lidar measurements was solved successfully via a hybrid method by a variable dimension of projection with B-Splines. Numerical simulation results with noisy data at different measurement situations show that it is possible to derive a reconstruction of the aerosol distribution only with 4 measurements.
The ill-posed inversion of multiwavelength lidar data by a hybrid method of variable projection
(1999)
Laser beam melt ablation - a contact-free machining process - offers several advantages compared to conventional processing mechanisms: there exists no tool wear and even extremely hard or brittle materials can be processed. During ablation the workpiece is molten by a CO2-laser beam, this melt is then driven out by the impulse of a process gas. The idea behind laser ablation is rather simple, but it has a major limitation in practical applications: with increasing ablation rates surface quality of the workpiece processed declines rapidly. At high ablation rates, depending on the process parameters different periodic-like structures can be observed on the ablated surface. These structures show a dependence on the line energy, which has been identified as a fundamental control parameter. In dependence on this parameter several regimes with different behaviours of the process have been separated. These regimes are distinguishable as well in the surfaces obtained as in the signals gained by the measurement of the process emissions. Further aim is to identify the different modes of the system and reach a deeper understanding of the dynamics of the molten material in order to understand the formation of these surface structures. With this it should be possible to influence the system in the direction of avoiding structure formation even at high ablation rates. Relying on the results on-line monitoring and control of the process should be studied.
We investigate numerically the appearance of heteroclinic behavior in a three-dimensional, buoyancy-driven fluid layer with stress-free top and bottom boundaries, a square horizontal periodicity with a small aspect ratio, and rotation at low to moderate rates about a vertical axis. The Prandtl number is 6.8. If the rotation is not too slow, the skewed-varicose instability leads from stationary rolls to a stationary mixed-mode solution, which in turn loses stability to a heteroclinic cycle formed by unstable roll states and connections between them. The unstable eigenvectors of these roll states are also of the skewed-varicose or mixed-mode type and in some parameter regions skewed-varicose like shearing oscillations as well as square patterns are involved in the cycle. Always present weak noise leads to irregular horizontal translations of the convection pattern and makes the dynamics chaotic, which is verified by calculating Lyapunov exponents. In the nonrotating case, the primary rolls lose, depending on the aspect ratio, stability to traveling waves or a stationary square pattern. We also study the symmetries of the solutions at the intermittent fixed points in the heteroclinic cycle.
We numerically investigate nonlinear asymmetric square patterns in a horizontal convection layer with up-down reflection symmetry. As a novel feature we find the patterns to appear via the skewed varicose instability of rolls. The time-independent nonlinear state is generated by two unstable checkerboard (symmetric square) patterns and their nonlinear interaction. As the bouyancy forces increase, the interacting modes give rise to bifurcations leading to a periodic alternation between a nonequilateral hexagonal pattern and the square pattern or to different kinds of standing oscillations.
In single photon emission computed tomography (SPECT) one is interested in reconstructing the activity distribution f of some radiopharmaceutical. The data gathered suffer from attenuation due to the tissue density µ. Each imaged slice incorporates noisy sample values of the nonlinear attenuated Radon transform (formular at this place in the original abstract) Traditional theory for SPECT reconstruction treats µ as a known parameter. In practical applications, however, µ is not known, but either crudely estimated, determined in costly additional measurements or plainly neglected. We demonstrate that an approximation of both f and µ from SPECT data alone is feasible, leading to quantitatively more accurate SPECT images. The result is based on nonlinear Tikhonov regularization techniques for parameter estimation problems in differential equations combined with Gauss-Newton-CG minimization.
Projection methods based on wavelet functions combine optimal convergence rates with algorithmic efficiency. The proofs in this paper utilize the approximation properties of wavelets and results from the general theory of regularization methods. Moreover, adaptive strategies can be incorporated still leading to optimal convergence rates for the resulting algorithms. The so-called wavelet-vaguelette decompositions enable the realization of especially fast algorithms for certain operators.
The structure of time series and letter sequences is investigated using the concepts of entropy and complexity. First conditional entropy and transinformation are introduced and several generalizations are discussed. Further several measures of complexity are introduced and discussed. The capability of these concepts to describe the structure of time series and letter sequences generated by nonlinear maps, data series from meteorology, astrophysics, cardiology, cognitive psychology and finance is investigated. The relation between the complexity and the predictability of informational strings is discussed. The relation between local order and the predictability of time series is investigated.
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.
Based on the data of the Magion2 subsatellite of the Intercosmos24 satellite, an example of small-scale irregularities of the electron concentration with linear dimensions l ~ 100-300 m in the polar ion- osphere of the morning sector under field-aligned currents at altitudes of 1800-2030 km during the main phase of the magnetic storm of June 13, 1990 is presented. The dependence of the spectral index of the above small-scale irregularities on latitude is determined for the first time. Certain mechanisms of the generation of these small-scale irregularities are also qualitatively discussed.
We report on bifurcation studies for the incompressible Navier-Stokes equations in two space dimensions with periodic boundary conditions and an external forcing of the Kolmogorov type. Fourier representations of velocity and pressure have been used to approximate the original partial differential equations by a finite-dimensional system of ordinary differential equations, which then has been studied by means of bifurcation-analysis techniques. A special route into chaos observed for increasing Reynolds number or strength of the imposed forcing is described. It includes several steady states, traveling waves, modulated traveling waves, periodic and torus solutions, as well as a period-doubling cascade for a torus solution. Lyapunov exponents and Kaplan-Yorke dimensions have been calculated to characterize the chaotic branch. While studying the dynamics of the system in Fourier space, we also have transformed solutions to real space and examined the relation between the different bifurcations in Fourier space and toplogical changes of the streamline portrait. In particular, the time-dependent solutions, such as, e.g., traveling waves, torus, and chaotic solutions, have been characterized by the associated fluid-particle motion (Lagrangian dynamics).
We have studied bifurcation phenomena for the incompressable Navier-Stokes equations in two space dimensions with periodic boundary conditions. Fourier representations of velocity and pressure have been used to transform the original partial differential equations into systems of ordinary differential equations (ODE), to which then numerical methods for the qualitative analysis of systems of ODE have been applied, supplemented by the simulative calculation of solutions for selected initial conditions. Invariant sets, notably steady states, have been traced for varying Reynolds number or strength of the imposed forcing, respectively. A complete bifurcation sequence leading to chaos is described in detail, including the calculation of the Lyapunov exponents that characterize the resulting chaotic branch in the bifurcation diagram.
We have studied the bifurcations in a three-dimensional incompressible magnetofluid with periodic boundary conditions and an external forcing of the Arnold-Beltrami-Childress (ABC) type. Bifurcation-analysis techniques have been applied to explore the qualitative behavior of solution branches. Due to the symmetry of the forcing, the equations are equivariant with respect to a group of transformations isomorphic to the octahedral group, and we have paid special attention to symmetry-breaking effects. As the Reynolds number is increased, the primary nonmagnetic steady state, the ABC flow, loses its stability to a periodic magnetic state, showing the appearance of a generic dynamo effect; the critical value of the Reynolds number for the instability of the ABC flow is decreased compared to the purely hydrodynamic case. The bifurcating magnetic branch in turn is subject to secondary, symmetry-breaking bifurcations. We have traced periodic and quasi- periodic branches until they end up in chaotic states. In particular detail we have analyzed the subgroup symmetries of the bifurcating periodic branches, which are closely related to the spatial structure of the magnetic field.
We report on bifurcation studies for the incompressible magnetohydrodynamic equations in three space dimensions with periodic boundary conditions and a temporally constant external forcing. Fourier reprsentations of velocity, pressure and magnetic field have been used to transform the original partial differential equations into systems of ordinary differential equations (ODE), to which then special numerical methods for the qualitative analysis of systems of ODE have been applied, supplemented by the simulative calculation of solutions for selected initial conditions. In a part of the calculations, in order to reduce the number of modes to be retained, the concept of approximate inertial manifolds has been applied. For varying (incereasing from zero) strength of the imposed forcing, or varying Reynolds number, respectively, time-asymptotic states, notably stable stationary solutions, have been traced. A primary non-magnetic steady state loses, in a Hopf bifurcation, stability to a periodic state with a non-vanishing magnetic field, showing the appearance of a generic dynamo effect. From now on the magnetic field is present for all values of the forcing. The Hopf bifurcation is followed by furhter, symmetry-breaking, bifurcations, leading finally to chaos. We pay particular attention to kinetic and magnetic helicities. The dynamo effect is observed only if the forcing is chosen such that a mean kinetic helicity is generated; otherwise the magnetic field diffuses away, and the time-asymptotic states are non-magnetic, in accordance with traditional kinematic dynamo theory.
The bifurcation behaviour of the 3D magnetohydrodynamic equations has been studied for external forcings of varying degree of helicity. With increasing strength of the forcing a primary non-magnetic steady state loses stability to a magnetic periodic state if the helicity exceeds a threshold value and to different non-magnetic states otherwise.
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.
We employ a spectral decomposition method to analyze synchronization of a non-identical oscillator network. We study the case that a small parameter mismatch of oscillators is characterized by one parameter and phase synchronization is observed. We derive a linearized equation for each eigenmode of the coupling matrix. The parameter mismatch is reflected on inhomogeneous term in the linearized equation. We find that the oscillation of each mode is essentially characterized only by the eigenvalue of the coupling matrix with a suitable normalization. We refer to this property as spectral universality, because it is observed irrespective of network topology. Numerical results in various network topologies show good agreement with those based on linearized equation. This universality is also observed in a system driven by additive independent Gaussian noise.
Recent studies have pointed out the importance of transient synchronization between widely distributed neural assemblies to understand conscious perception. These neural assemblies form intricate networks of neurons and synapses whose detailed map for mammals is still unknown and far from our experimental capabilities. Only in a few cases, for example the C. elegans, we know the complete mapping of the neuronal tissue or its mesoscopic level of description provided by cortical areas. Here we study the process of transient and global synchronization using a simple model of phase-coupled oscillators assigned to cortical areas in the cerebral cat cortex. Our results highlight the impact of the topological connectivity in the developing of synchronization, revealing a transition in the synchronization organization that goes from a modular decentralized coherence to a centralized synchronized regime controlled by a few cortical areas forming a Rich-Club connectivity pattern.
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
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
In one of the data mining techniques, change-point detection is of importance in evaluating time series measured in real world. For decades this technique has been developed as a nonlinear dynamics. We apply the method for detecting the change points, Singular Spectrum Transformation (SST), to the climate time series. To know where the structures of climate data sets change can reveal a climate background. In this paper we discuss the structures of precipitation data in Kenya and Wrangel Island (Arctic land) by using the SST.
The paper presents a method that determines, by standard numerical means, the type of mutual relations of fold and flip bifurcations (configured as a so-called communication area) of a map. Equation systems are developed for the computation of points where a transition between areas of different types occurs. Furthermore, it is shown that saddle area<->spring area transitions can exist which have not yet been considered in the literature. Analytical conditions of that transition are derived.
Contents: I. Algorithms 1. Theoretical Backround 2. Numerical Procedures 3. Graph Representation of the Solutions 4. Applications and Example II. Users' Manual 5. About the Program 6. The Course of a Qualitative Analysis 7. The Model Module 8. Input description 9. Output Description 10. Example 11. Graphics
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 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.
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 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.
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.
In the present work, phenomena in the ionosphere are studied, which are connected with earthquakes (16 events) having a depth of less than 50 km and a magnitude M larger than 4. Analysed are night-time Es-spread effects using data of the vertical sounding station Petropavlovsk- Kanchatsky (φ=53.0°, λ=158.7°) from May 2004 until August 2004 registered every 15 minutes. It is found that the maximum distance of the earthquake from the sounding station, where pre-seismic phenomena are yet observable, depends on the magnitude of the earthquake. Further it is shown that 1-2 days before the earthquakes, in the premidnight hours, the appearance of Es-spread increases. The reliability of this increase amounts to 0.95.
The statistical analysis of the variations of the dayly-mean frequency of the maximum ionospheric electron density foF2 is performed in connection with the occurrence of (more than 60) earthquakes with magnitudes M > 6.0, depths h < 80 km and distances from the vertical sounding station R < 1000 km. For the study, data of the Tokyo sounding station are used, which were registered every hour in the years 1957-1990. It is shown that, on the average, foF2 decreases before the earthquakes. One day before the shock the decrease amounts to about 5 %. The statistical reliability of this phenomenon is obtained to be better than 0.95. Further, the variations of the occurrence probability of the turbulization of the F-layer (F spread) are investigated for (more than 260) earthquakes with M > 5.5, h < 80 km, R < 1000 km. For the analysis, data of the Japanese station Akita from 1969-1990 are used, which were obtained every hour. It is found that before the earthquakes the occurrence probability of F spread decreases. In the week before the event, the decrease has values of more than 10 %. The statistical reliability of this phenomenon is also larger than 0.95. Examining the seismo-ionospheric effects, here periods of time with weak heliogeomagnetic disturbances are considered, the Wolf number is less than 100 and the index ∑ Kp is smaller than 30.
A model of the generation of pulses of local electric fields with characteristic time scales of 1–10 minutes is considered for atmospheric conditions above fracture regions of earthquakes. In the model, it is proposed that aerosols, increased ionization velocity and upstreaming air flows occur at night-time conditions. The pulses of local electric fields cause respective pulses of infrared emissions. But infrared emissions with time scales of 1–10 minutes were not observed up to now experimentally. The authors think, that the considered non-stationary field and radiation effects might be a new-type of applicable earthquake indicators and ask to perform special earth-based and satellite observations of the night-time atmosphere in seismoactive fracture regions.
Contents: 1 Introduction 2 Formation and destruction of sporadic E-layers 3 Temporal variations of parameters of sporadic E-layers during earthquake preparation 3.1 Temporal variations of fbEs with time-scales of a few hours 3.2 Study of fbEs variations with characteristic time-scales of 0.5-3 hours 3.3 Variations of the parameters of sporadic E-layers with characteristic time-scales of 15-45 minutes 3.4 Sporadic E-layer variations with characteristic time-scales of 2-15 minutes 4 On the spatial scales of sporadic E-layer disturbances related to seismic activity 5 Complex experimental researches of the ionosphere, electromagnetic noise and the geomagnetic field 5.1 Ionospheric and electromagnetic phenomena of the Kayraccum earthquake in 1985 5.2 Comparison of anomalies with characteristic time-scales of 2-3 hours for ionospheric E- and F-layers, and temporal behaviour of electromagnetic noise emission intensity 5.3 Night airglow emissions in the E-region before earthquakes and sporadic E-layer variations 6 Physical models of lithosphere-ionosphere links 6.1 Lithosphere-ionosphere links due to AGW 6.2 Electromagnetic models for the lithosphere-ionosphere coupling 6.3 Sporadic E-layers as current generators 7 Discussion and conclusion
This paper deals with the electrical conductivity problem in geophysics. It is formulated as an elliptic boundary value problem of second order for a large class of bounded and unbounded domains. A special boundary condition, the so called "Complete Electrode Model", is used. Poincaré inequalities are formulated and proved in the context of weighted Sobolev spaces, leading to existence and uniqueness statements for the boundary value problem. In addition, a parameter-to-solution operator arising from the inverse conductivity problem in medicine (EIT) and geophysics is investigated mathematically and is shown to be smooth and analytic.
The aim of this paper is to describe an efficient strategy for descritizing ill-posed linear operator equations of the first kind: we consider Tikhonov-Phillips-regularization χ^δ α = (a * a + α I)^-1 A * y ^δ with a finite dimensional approximation A n instead of A. We propose a sparse matrix structure which still leads to optimal convergences rates but requires substantially less scalar products for computing A n compared with standard methods.
Contents: 1 Introduction 1.1 Tikhanov-Phillips Regularization of Ill-Posed Problems 1.2 A Compact Course to Wavelets 2 A Multilevel Iteration for Tikhonov-Phillips Regularization 2.1 Multilevel Splitting 2.2 The Multilevel Iteration 2.3 Multilevel Approach to Cone Beam Reconstuction 3 The use of approximating operators 3.1 Computing approximating families {Ah}
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.
In this paper an analysis of the excitation conditions of mirror waves is done, which propagate parallel to an external magnetic field. There are found analytical expressions for the dispersion relations of the waves in case of different plasma conditions. These relations may be used in future to develop the nonlinear theory of mirror waves. In comparison with former analytical works, in the study the inuence of the magnetic field and nite temperatures of the ions parallel to the magnetic field are taken into account. Application is done for the earth's magnetosheath.
This paper deals with the Mie scattering kernels for multi-spectral data. The kernels may be represented in form of power series. Furthermore, the singular-value spectrum and the degree of ill-posedness in dependence on the refractive index of the particles are numerically approximated. A special hybrid regularization technique allows us to determine via inversion the particle distributions of different types.
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.
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.
Strange nonchaotic attractors typically appear in quasiperiodically driven nonlinear systems. Two methods of their characterization are proposed. The first one is based on the bifurcation analysis of the systems, resulting from periodic approximations of the quasiperiodic forcing. Secondly, we propose th characterize their strangeness by calculating a phase sensitivity exponent, that measures the sensitivity with respect to changes of the phase of the external force. It is shown, that phase sensitivity appears if there is a non-zero probability for positive local Lyapunov exponents to occur.
We demonstrate the occurrence of regimes with singular continuous (fractal) Fourier spectra in autonomous dissipative dynamical systems. The particular example in an ODE system at the accumulation points of bifurcation sequences associated to the creation of complicated homoclinic orbits. Two different machanisms responsible for the appearance of such spectra are proposed. In the first case when the geometry of the attractor is symbolically represented by the Thue-Morse sequence, both the continuous-time process and its descrete Poincaré map have singular power spectra. The other mechanism owes to the logarithmic divergence of the first return times near the saddle point; here the Poincaré map possesses the discrete spectrum, while the continuous-time process displays the singular one. A method is presented for computing the multifractal characteristics of the singular continuous spectra with the help of the usual Fourier analysis technique.
The 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.
The investigation of foetal reaction to internal and external conditions and stimuli is an important tool in the characterization of the developing neural integration of the foetus. An interesting example of this is the study of the interrelationship between the foetal and the maternal heart rate. Recent studies have shown a certain likelihood of occasional heart rate synchronization between mother and foetus. In the case of respiratory-induced heart rate changes, the comparison with maternal surrogates suggests that the evidence for detected synchronization is largely statistical and does not result from physiological interaction. Rather, they simply reflect a stochastic, temporary stability of two independent oscillators with time-variant frequencies. We reanalysed three datasets from that study for a more local consideration. Epochs of assumed synchronization associated with short-term regulation of the foetal heart rate were selected and compared with synchronization resulting from white noise instead of the foetal signal. Using data-driven modelling analysis, it was possible to identify the consistent influence of the heartbeat duration of maternal beats preceding the foetal beats during epochs of synchronization. These maternal beats occurred approximately one maternal respiratory cycle prior to the affected foetal beat. A similar effect could not be found in the epochs without synchronization. Simulations based on the fitted models led to a higher likelihood of synchronization in the data segments with assumed foetal-maternal interaction than in the segment without such assumed interaction. We conclude that the data-driven model-based analysis can be a useful tool for the identification of synchronization.
We 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.
The dynamics of tail-like current sheets under the influence of small-scale plasma turbulence
(1999)
A 2D-magnetohydrodynamic model of current-sheet dynamics caused by anomalous electrical resistivity as result of small-scale plasma turbulence is proposed. The anomalous resistivity is assumed to be proportional to the square of the gradient of the magnetic pressure as may be valid for instance in the case of lower-hybrid-drift turbulence. The initial resistivity pulse is given. Then the temporal and spatial evolution of the magnetic and electric fields, plasma density, pressure, convection and resistivity are considered. The motion of the induced electric field is discussed as indicator of the plasma disturbances. The obtained results found using much improved numerical methods show a magnetic field evolution with x-line formation and plasma acceleration. Besides, in the current sheet, three types of magnetohydrodynamic waves occur, fast magnetoacoustic waves of compression and rarefaction as well as slow magnetoacoustic waves.
The bifurcations in a three-dimensional incompressible, electrically conducting fluid with an external forcing of the Roberts type have been studied numerically. The corresponding flow can serve as a model for the convection in the outer core of the Earth and is realized in an ongoing laboratory experiment aimed at demonstrating a dynamo effect. The symmetry group of the problem has been determined and special attention has been paid to symmetry breaking by the bifurcations. The nonmagnetic, steady Roberts flow loses stability to a steady magnetic state, which in turn is subject to secondary bifurcations. The secondary solution branches have been traced until they end up in chaotic states.
Three-dimensional bouyancy-driven convection in a horizontal fluid layer with stress-free boundary conditions at the top and bottom and periodic boundary conditions in the horizontal directions is investigated by means of numerical simulation and bifurcation-analysis techniques. The aspect ratio is fixed to a value of 2√2 and the Prandtl number to a value of 6.8. Two-dimensional convection rolls are found to be stable up to a Rayleigh number of 17 950, where a Hopf bifurcation leads to traveling waves. These are stable up to a Rayleigh number of 30 000, where a secondary Hopf bifurcation generates modulated traveling waves. We pay particular attention to the symmetries of the solutions and symmetry breaking by the bifurcations.
In the recent past, recurrence quantification analysis (RQA) has gained an increasing interest in various research areas. The complexity measures the RQA provides have been useful in describing and analysing a broad range of data. It is known to be rather robust to noise and nonstationarities. Yet, one key question in empirical research concerns the confidence bounds of measured data. In the present Letter we suggest a method for estimating the confidence bounds of recurrence-based complexity measures. We study the applicability of the suggested method with model and real- life data.
The EEG is one of the most commonly used tools in brain research. Though of high relevance in research, the data obtained is very noisy and nonstationary. In the present article we investigate the applicability of a nonlinear data analysis method, the recurrence quantification analysis (RQA), to Such data. The method solely rests on the natural property of recurrence which is a phenomenon inherent to complex systems, such as the brain. We show that this method is indeed suitable for the analysis of EEG data and that it might improve contemporary EEG analysis.
The present paper is related to the problem of approximating the exact solution to the magnetohydrodynamic equations (MHD). The behaviour of a viscous, incompressible and resistive fluid is exemined for a long period of time. Contents: 1 The magnetohydrodynamic equations 2 Notations and precise functional setting of the problem 3 Existence, uniqueness and regularity results 4 Statement and Proof of the main theorem 5 The approximate inertial manifold 6 Summary
The usage of nonlinear Galerkin methods for the numerical solution of partial differential equations is demonstrated by treating an example. We desribe the implementation of a nonlinear Galerkin method based on an approximate inertial manifold for the 3D magnetohydrodynamic equations and compare its efficiency with the linear Galerkin approximation. Special bifurcation points, time-averaged values of energy and enstrophy as well as Kaplan-Yorke dimensions are calculated for both schemes in order to estimate the number of modes necessary to correctly describe the behavior of the exact solutions.
A numerical bifurcation analysis of the electrically driven plane sheet pinch is presented. The electrical conductivity varies across the sheet such as to allow instability of the quiescent basic state at some critical Hartmann number. The most unstable perturbation is the two-dimensional tearing mode. Restricting the whole problem to two spatial dimensions, this mode is followed up to a time-asymptotic steady state, which proves to be sensitive to three-dimensional perturbations even close to the point where the primary instability sets in. A comprehensive three-dimensional stability analysis of the two-dimensional steady tearing-mode state is performed by varying parameters of the sheet pinch. The instability with respect to three-dimensional perturbations is suppressed by a sufficiently strong magnetic field in the invariant direction of the equilibrium. For a special choice of the system parameters, the unstably perturbed state is followed up in its nonlinear evolution and is found to approach a three-dimensional steady state.
Using quantities of symbolic dynamics, such as mutual information, Shannon information and algorithmic complexity, we have searched for interrelations of spikes emitted simultaneously at different frequencies during the impulsive phase of a flare event. As the spikes are related to the flare energy release and are interpreted as emissions originating at different sites having different magnetic field strengths, any relation in frequency is interpretated as a relation in space. This approach is appropriate to characterize such spatio-temporal patterns, whereas the popular estimate of fractal dimensions can be applied to low-dimensional systems only. Depending on the energy release and emission processes, two types of fragmentation are possible: a scenario of global organization (spikes are emitted in a succession of similar events by the same system) or a scenario of local organization (many systems triggered by an initial event). Mutual information which is a generalization of correlation indicates a relation in frequency beyond the bandwidth of individual spikes. The scans in the spectrograms with large mutual information also show a low level of Shannon information and algorithmic complexity, indicating that the simultaneous appearance of spikes at other frequencies is not a completely stochastic phenomenon (white noise). It may be caused by a nonlinear deterministic system or by a Markov process. By means of mutual information we find a memory over frequency intervals up to 60 MHz. Shannon information and algorithmic complexity concern the mbox{whole} frequency region, i.e. the global source region. A global organization is also apparent in quasi-periodic changes of the Shannon information and algorithmic complexity in the range of 2 - 8 seconds. The finding is compatible with a scenario of local organization in which the information of one event spreads spatially and triggers further events at different places. The region is not an ensemble of independently flashing sources, each representing a system that cascades in energy after an initial trigger. On the contrary, there is a causal connection between the sources at any time. The analysis of the four spike events suggests that the structure in frequency is not stochastic but a process in which spikes at nearby locations are simultaneously triggered by a common exciter.
It is shown that the ff effect of mean-field magnetohydrodynamics, which consists in the generation of a mean electromotive force along the mean magnetic field by turbulently fluctuating parts of velocity and magnetic field, is equivalent to the simultaneous generation of both turbulent and mean-field magnetic helicities, the generation rates being equal in magnitude and opposite in sign. In the particular case of statistically stationary and homogeneous fluctuations this implies that the ff effect can increase the energy in the mean magnetic field only under the condition that also magnetic helicity is accumulated there.
The stability of the quiescent ground state of an incompressible, viscous and electrically conducting fluid sheet, bounded by stress-free parallel planes and driven by an external electric field tangential to the boundaries, is studied numerically. The electrical conductivity varies as cosh–2(x1/a), where x1 is the cross-sheet coordinate and a is the half width of a current layer centered about the midplane of the sheet. For a <~ 0.4L, where L is the distance between the boundary planes, the ground state is unstable to disturbances whose wavelengths parallel to the sheet lie between lower and upper bounds depending on the value of a and on the Hartmann number. Asymmetry of the configuration with respect to the midplane of the sheet, modelled by the addition of an externally imposed constant magnetic field to a symmetric equilibrium field, acts as a stabilizing factor.
The stability of the quiescent ground state of an incompressible viscous fluid sheet bounded by two parallel planes, with an electrical conductivity varying across the sheet, and driven by an external electric field tangential to the boundaries is considered. It is demonstrated that irrespective of the conductivity profile, as magnetic and kinetic Reynolds numbers (based on the Alfvén velocity) are raised from small values, two-dimensional perturbations become unstable first.
We have numerically studied the bifurcation properties of a sheet pinch with impenetrable stress-free boundaries. An incompressible, electrically conducting fluid with spatially and temporally uniform kinematic viscosity and magnetic diffusivity is confined between planes at x1=0 and 1. Periodic boundary conditions are assumed in the x2 and x3 directions and the magnetofluid is driven by an electric field in the x3 direction, prescribed on the boundary planes. There is a stationary basic state with the fluid at rest and a uniform current J=(0,0,J3). Surprisingly, this basic state proves to be stable and apparently to be the only time-asymptotic state, no matter how strong the applied electric field and irrespective of the other control parameters of the system, namely, the magnetic Prandtl number, the spatial periods L2 and L3 in the x2 and x3 directions, and the mean values B¯2 and B¯3 of the magnetic-field components in these directions.
Acoustic emission signals generated during high speed cutting of steel are investigated. The data are represen ted in time-folded form. Several methods from linear and nonlinear data analysis based on time- and frequency- domain are applied to the data and reveal signatures of the observed acoustic emission signal. These investiga tions are necessary for modeling the cutting process by means of differential equations.
Charged dust grains in circumplanetary environments experience, beyond various deterministic forces, also stochastic perturbations caused, by fluctuations of the magnetic field, the charge of the grains, by chaotic rotation of aspherical grains, etc. Here we investigate the dynamics of a dust population in a circular orbit around a planet which is perturbed by a stochastic planetary magnetic field B', modeled by an isotropically Gaussian white noise. The resulting perturbation equations give rise to a modified diffusion of the inclinations i and eccentricities e. The diffusion coefficient is found to be D proportional to w^2 O /n^2 , where the gyrofrequency, the Kepler frequency, and the synodic frequency are denoted by w , O, and n, respectively. This behavior has been checked against numerical simulations. We have chosen dust grains (1 m in radius) ejected from Jupiter's satellite Europa in circular equatorial orbits around Jupiter and integrated numerically their trajectories over their typical lifetimes (100 years). The particles were exposed to a Gaussian fluctuating magnetic field B' with the same statistical properties as in the analytical treatment. These simulations have confirmed the analytical results. The theoretical studies showed the statistical properties of B' to be of decisive importance. To estimate them, we analyzed the magnetic field data obtained by the Galileo spacecraft magnetometer at Jupiter and found almost Gaussian fluctuations of about 5% of the mean field and exponentially decaying correlations. This results in a diffusion of orbital inclinations and eccentricities of the dust grains of about ten percent over the lifetime of the particles. For smaller dusty motes or for close-in particles (e.g., in Jovian gossamer rings) stochastics might well dominate the dynamics.
We analyse the X-ray light curves of compact objects using linear and nonlinear time series analysis methods. A Power Density Spectrum (PDS) describes the overall second order properties of the observed data well. To look beyond we propose the nonlinear Q-statistic to detect an asymmetry of the time series. This allows us to find relevant time scales. This method even grants a subclassification of the known states of X-ray sources.
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 analyze the variability in the x-ray lightcurves of the black hole candidate Cygnus X-1 by linear and nonlinear time series analysis methods. While a linear model describes the overall second order properties of the observed data well, surrogate data analysis reveals a significant deviation from linearity. We discuss the relation between shot noise models usually applied to analyze these data and linear stochastic autoregressive models. We debate statistical and interpretational issues of surrogate data testing for the present context. Finally, we suggest a combination of tools from linear and nonlinear time series analysis methods as a procedure to test the predictions of astrophysical models on observed data.
Basing on recent solar models, the excitation of ion-acoustic turbulence in the weaklycollisional, fully and partially-ionized regions of the solar atmosphere is investigated. Within the frame of hydrodynamics, conditions are found under which the heating of the plasma by ion-acoustic type waves is more effective than the Joule heating. Taking into account wave and Joule heating effects, a nonlinear differential equation is derived, which describes the evolution of nonlinear ion-acoustic waves in the collisional plasma.
The nonlinear interaction of waves excited by the modified two-stream instability (Farley-Buneman instability) is considered. It is found that, during the linear stage of wave growth, the enhanced pressure of the high-frequency part of the waves locally generates a ponderomotive force. This force acts on the plasma particles and redistributes them. Thus an additional electrostatic polarization field occurs, which influences the low-frequency part of the waves. Then, the low-frequency waves also cause a redistribution of the high-frequency waves. In the paper, a self-consistent system of equations is obtained, which describes the nonlinear interaction of the waves. It is shown that the considered mechanism of wave interaction causes a nonlinear stabilization of the high-frequency waves’ growth and a formation of local density structures of the charged particles. The density modifications of the charged particles during the non-linear stage of wave growth and the possible interval of aspect angles of the high-frequency waves are estimated.
Using a special technique of data analysis, we have found out 34 grand minima of solar activity obtained from a 7,700 years long Δ14C record. The method used rests on a proper filtering of the Δ14C record and the extrapolation of verifiable results for the later history back in time. Additionally, we use a method of nonlinear dynamics, the recurrence rate, to back up the results. Our findings are not contradictory to the record of solar maxima resp. minima by Eddy [5], but constitute a considerable extension. Hence, it has become possible to look closer at the validity of models. This way, we have tested several models for solar activity, esp. the model of Barnes et al. [1]. There are hints for that the grand minima might solely be driven by the 209 year period found in the Δ14C record.
Using a special technique of data analysis, we have found out 34 grand minima of solar activity in a 7,700 years long C14 record. The method used rests on a proper filtering of the C14 record and the extrapolation of verifiable results for the later history back in time. Additionally, we have applied a method of nonlinear dynamics, the recurrence rate, to back up the results. Our findings are not contradictory to the record of grand minima by Eddy, but constitute a considerable extension. Hence, it has become possible to look closer at the validity of models. This way, we have tested esp. the model of Barnes et al. There are hints for that the grand minima might solely be driven by the 209--year period found in the C14 record.