Dynamic scaling of bred vectors in spatially extended chaotic systems

We unfold a profound relationship between the dynamics of finite-size perturbations in spatially extended chaotic systems and the universality class of Kardar-Parisi-Zhang (KPZ). We show how this relationship can be exploited to obtain a complete theoretical description of the bred vectors dynamics. The existence of characteristic length/time scales, the spatial extent of spatial correlations and how to time it, and the role of the breeding amplitude are all analyzed in the light of our theory. Implications to weather forecasting based on ensembles of initial conditions are also discussed.

Four-dimensional data assimilation and balanced dynamics

The problem of spurious excitation of gravity waves in the context of four-dimensional data assimilation is investigated using a simple model of balanced dynamics. The model admits a chaotic vortical mode coupled to a comparatively fast gravity wave mode, and can be initialized such that the model evolves on a so-called slow manifold, where the fast motion is suppressed. Identical twin assimilation experiments are performed, comparing the extended and ensemble Kalman filters (EKF and EnKF, respectively). The EKF uses a tangent linear model (TLM) to estimate the evolution of forecast error statistics in time, whereas the EnKF uses the statistics of an ensemble of nonlinear model integrations. Specifically, the case is examined where the true state is balanced, but observation errors project onto all degrees of freedom, including the fast modes. It is shown that the EKF and EnKF will assimilate observations in a balanced way only if certain assumptions hold, and that, outside of ideal cases (i.e., with very frequent observations), dynamical balance can easily be lost in the assimilation. For the EKF, the repeated adjustment of the covariances by the assimilation of observations can easily unbalance the TLM, and destroy the assumptions on which balanced assimilation rests. It is shown that an important factor is the choice of initial forecast error covariance matrix. A balance-constrained EKF is described and compared to the standard EKF, and shown to offer significant improvement for observation frequencies where balance in the standard EKF is lost. The EnKF is advantageous in that balance in the error covariances relies only on a balanced forecast ensemble, and that the analysis step is an ensemble-mean operation. Numerical experiments show that the EnKF may be preferable to the EKF in terms of balance, though its validity is limited by ensemble size. It is also found that overobserving can lead to a more unbalanced forecast ensemble and thus to an unbalanced analysis.

Averaging, slaving and balance dynamics in a simple atmospheric model

We report numerical results from a study of balance dynamics using a simple model of atmospheric motion that is designed to help address the question of why balance dynamics is so stable. The non-autonomous Hamiltonian model has a chaotic slow degree of freedom (representing vortical modes) coupled to one or two linear fast oscillators (representing inertia-gravity waves). The system is said to be balanced when the fast and slow degrees of freedom are separated. We find adiabatic invariants that drift slowly in time. This drift is consistent with a random-walk behaviour at a speed which qualitatively scales, even for modest time scale separations, as the upper bound given by Neishtadt’s and Nekhoroshev’s theorems. Moreover, a similar type of scaling is observed for solutions obtained using a singular perturbation (‘slaving’) technique in resonant cases where Nekhoroshev’s theorem does not apply. We present evidence that the smaller Lyapunov exponents of the system scale exponentially as well. The results suggest that the observed stability of nearly-slow motion is a consequence of the approximate adiabatic invariance of the fast motion.

Chaotic mixing and transport in Rossby-wave critical layers

A simple, dynamically consistent model of mixing and transport in Rossby-wave critical layers is obtained from the well-known Stewartson–Warn–Warn (SWW) solution of Rossby-wave critical-layer theory. The SWW solution is thought to be a useful conceptual model of Rossby-wave breaking in the stratosphere. Chaotic advection in the model is a consequence of the interaction between a stationary and a transient Rossby wave. Mixing and transport are characterized separately with a number of quantitative diagnostics (e.g. mean-square dispersion, lobe dynamics, and spectral moments), and with particular emphasis on the dynamics of the tracer field itself. The parameter dependences of the diagnostics are examined: transport tends to increase monotonically with increasing perturbation amplitude whereas mixing does not. The robustness of the results is investigated by stochastically perturbing the transient-wave phase speed. The two-wave chaotic advection model is contrasted with a stochastic single-wave model. It is shown that the effects of chaotic advection cannot be captured by stochasticity alone.

On Hamiltonian balanced dynamics and the slowest invariant manifold

The concept of a slowest invariant manifold is investigated for the five-component model of Lorenz under conservative dynamics. It is shown that Lorenz's model is a two-degree-of-freedom canonical Hamiltonian system, consisting of a nonlinear vorticity-triad oscillator coupled to a linear gravity wave oscillator, whose solutions consist of regular and chaotic orbits. When either the Rossby number or the rotational Froude number is small, there is a formal separation of timescales, and one can speak of fast and slow motion. In the same regime, the coupling is weak, and the Kolmogorov–Arnold-Moser theorem is shown to apply. The chaotic orbits are inherently unbalanced and are confined to regions sandwiched between invariant tori consisting of quasi-periodic regular orbits. The regular orbits generally contain free fast motion, but a slowest invariant manifold may be geometrically defined as the set of all slow cores of invariant tori (defined by zero fast action) that are smoothly related to such cores in the uncoupled system. This slowest invariant manifold is not global; in fact, its structure is fractal; but it is of nearly full measure in the limit of weak coupling. It is also nonlinearly stable. As the coupling increases, the slowest invariant manifold shrinks until it disappears altogether. The results clarify previous definitions of a slowest invariant manifold and highlight the ambiguity in the definition of “slowness.” An asymptotic procedure, analogous to standard initialization techniques, is found to yield nonzero free fast motion even when the core solutions contain none. A hierarchy of Hamiltonian balanced models preserving the symmetries in the original low-order model is formulated; these models are compared with classic balanced models, asymptotically initialized solutions of the full system and the slowest invariant manifold defined by the core solutions. The analysis suggests that for sufficiently small Rossby or rotational Froude numbers, a stable slowest invariant manifold can be defined for this system, which has zero free gravity wave activity, but it cannot be defined everywhere. The implications of the results for more complex systems are discussed.

Complex dynamics in simple systems with periodic parameter oscillations

We study systems with periodically oscillating parameters that can give way to complex periodic or nonperiodic orbits. Performing the long time limit, we can define ergodic averages such as Lyapunov exponents, where a negative maximal Lyapunov exponent corresponds to a stable periodic orbit. By this, extremely complicated periodic orbits composed of contracting and expanding phases appear in a natural way. Employing the technique of ϵ-uncertain points, we find that values of the control parameters supporting such periodic motion are densely embedded in a set of values for which the motion is chaotic. When a tiny amount of noise is coupled to the system, dynamics with positive and with negative nontrivial Lyapunov exponents are indistinguishable. We discuss two physical systems, an oscillatory flow inside a duct and a dripping faucet with variable water supply, where such a mechanism seems to be responsible for a complicated alternation of laminar and turbulent phases.

Dynamics of two planets in co-orbital motion

We study the stability regions and families of periodic orbits of two planets locked in a co-orbital configuration. We consider different ratios of planetary masses and orbital eccentricities; we also assume that both planets share the same orbital plane. Initially, we perform numerical simulations over a grid of osculating initial conditions to map the regions of stable/chaotic motion and identify equilibrium solutions. These results are later analysed in more detail using a semi-analytical model. Apart from the well-known quasi-satellite orbits and the classical equilibrium Lagrangian points L(4) and L(5), we also find a new regime of asymmetric periodic solutions. For low eccentricities these are located at (delta lambda, delta pi) = (+/- 60 degrees, -/+ 120 degrees), where delta lambda is the difference in mean longitudes and delta pi is the difference in longitudes of pericentre. The position of these anti-Lagrangian solutions changes with the mass ratio and the orbital eccentricities and are found for eccentricities as high as similar to 0.7. Finally, we also applied a slow mass variation to one of the planets and analysed its effect on an initially asymmetric periodic orbit. We found that the resonant solution is preserved as long as the mass variation is adiabatic, with practically no change in the equilibrium values of the angles.

Chaotic transport in reversed shear tokamaks

For tokamak models using simplified geometries and reversed shear plasma profiles, we have numerically investigated how the onset of Lagrangian chaos at the plasma edge may affect the plasma confinement in two distinct but closely related problems. Firstly, we have considered the motion of particles in drift waves in the presence of an equilibrium radial electric field with shear. We have shown that the radial particle transport caused by this motion is selective in phase space, being determined by the resonant drift waves and depending on the parameters of both the resonant waves and the electric field profile. Moreover, we have shown that an additional transport barrier may be created at the plasma edge by increasing the electric field. In the second place, we have studied escape patterns and magnetic footprints of chaotic magnetic field lines in the region near a tokamak wall, when there are resonant modes due to the action of an ergodic magnetic limiter. A non-monotonic safety factor profile has been used in the analysis of field line topology in a region of negative magnetic shear. We have observed that, if internal modes are perturbed, the distributions of field line connection lengths and magnetic footprints exhibit spatially localized escape channels. For typical physical parameters of a fusion plasma, the two Lagrangian chaotic processes considered in this work can be effective in usual conditions so as to influence plasma confinement. The reversed shear effects discussed in this work may also contribute to evaluate the transport barrier relevance in advanced confinement scenarios in future tokamak experiments.

Nonlinear Schrodinger equation with chaotic, random, and nonperiodic nonlinearity

In this Letter we deal with a nonlinear Schrodinger equation with chaotic, random, and nonperiodic cubic nonlinearity. Our goal is to study the soliton evolution, with the strength of the nonlinearity perturbed in the space and time coordinates and to check its robustness under these conditions. Here we show that the chaotic perturbation is more effective in destroying the soliton behavior, when compared with random or nonperiodic perturbation. For a real system, the perturbation can be related to, e.g., impurities in crystalline structures, or coupling to a thermal reservoir which, on the average, enhances the nonlinearity. We also discuss the relevance of such random perturbations to the dynamics of Bose-Einstein condensates and their collective excitations and transport. (C) 2010 Elsevier B.V. All rights reserved.

A note on learning chaotic sunspot equilibrium

In this paper we prove convergence to chaotic sunspot equilibrium through two learning rules used in the bounded rationality literature. The rst one shows the convergence of the actual dynamics generated by simple adaptive learning rules to a probability distribution that is close to the stationary measure of the sunspot equilibrium; since this stationary measure is absolutely continuous it results in a robust convergence to the stochastic equilibrium. The second one is based on the E-stability criterion for testing stability of rational expectations equilibrium, we show that the conditional probability distribution de ned by the sunspot equilibrium is expectational stable under a reasonable updating rule of this parameter. We also report some numerical simulations of the processes proposed.

DYNAMICS AT INFINITY of A CUBIC CHUA'S SYSTEM

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Dynamics at infinity and the existence of singularly degenerate heteroclinic cycles in the Lorenz system

In this paper, by using the Poincare compactification in R(3) we make a global analysis of the Lorenz system, including the complete description of its dynamic behavior on the sphere at infinity. Combining analytical and numerical techniques we show that for the parameter value b = 0 the system presents an infinite set of singularly degenerate heteroclinic cycles, which consist of invariant sets formed by a line of equilibria together with heteroclinic orbits connecting two of the equilibria. The dynamical consequences related to the existence of such cycles are discussed. In particular a possibly new mechanism behind the creation of Lorenz-like chaotic attractors, consisting of the change in the stability index of the saddle at the origin as the parameter b crosses the null value, is proposed. Based on the knowledge of this mechanism we have numerically found chaotic attractors for the Lorenz system in the case of small b > 0, so nearby the singularly degenerate heteroclinic cycles.

Bifurcation analysis of a new Lorenz-like chaotic system

In this paper we study the local codimension one and two bifurcations which occur in a family of three-dimensional vector fields depending on three parameters. An equivalent family, depending on five parameters, was recently proposed as a new chaotic system with a Lorenz-like butterfly shaped attractor and was studied mainly from a numerical point of view, for particular values of the parameters, for which computational evidences of the chaotic attractor was shown. In order to contribute to the understand of this new system we present an analytical study and the bifurcation diagrams of an equivalent three parameter system, showing the qualitative changes in the dynamics of its solutions, for different values of the parameters. (C) 2007 Elsevier Ltd. All rights reserved.

Global dynamics of the Rikitake system

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

The Orbital Dynamics of Synchronous Satellites: Irregular Motions in the 2:1 Resonance

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)