Whistler mode wave growth and propagation in the prenoon magnetosphere

Pitch-angle scattering of electrons can limit the stably trapped particle flux in the magnetosphere and precipitate energetic electrons into the ionosphere. Whistler-mode waves generated by a temperature anisotropy can mediate this pitch-angle scattering over a wide range of radial distances and latitudes, but in order to correctly predict the phase-space diffusion, it is important to characterise the whistler-mode wave distributions that result from the instability. We use previously-published observations of number density, pitch-angle anisotropy and phase space density to model the plasma in the quiet pre-noon magnetosphere (defined as periods when AE<100nT). We investigate the global propagation and growth of whistler-mode waves by studying millions of growing ray paths and demonstrate that the wave distribution at any one location is a superposition of many waves at different points along their trajectories and with different histories. We show that for observed electron plasma properties, very few raypaths undergo magnetospheric reflection, most rays grow and decay within 30 degrees of the magnetic equator. The frequency range of the wave distribution at large L can be adequately described by the solutions of the local dispersion relation, but the range of wavenormal angle is different. The wave distribution is asymmetric with respect to the wavenormal angle. The numerical results suggest that it is important to determine the variation of magnetospheric parameters as a function of latitude, as well as local time and L-shell.

Non-ergodicity of inviscid two-dimensional flow on a beta-plane and on the surface of a rotating sphere

It is shown that, for a sufficiently large value of β, two-dimensional flow on a doubly-periodic beta-plane cannot be ergodic (phase-space filling) on the phase-space surface of constant energy and enstrophy. A corresponding result holds for flow on the surface of a rotating sphere, for a sufficiently rapid rotation rate Ω. This implies that the higher-order, non-quadratic invariants are exerting a significant influence on the statistical evolution of the flow. The proof relies on the existence of a finite-amplitude Liapunov stability theorem for zonally symmetric basic states with a non-vanishing absolute-vorticity gradient. When the domain size is much larger than the size of a typical eddy, then a sufficient condition for non-ergodicity is that the wave steepness ε < 1, where ε = 2[surd radical]2Z/βU in the planar case and $\epsilon = 2^{\frac{1}{4}} a^{\frac{5}{2}}Z^{\frac{7}{4}}/\Omega U^{\frac{5}{2}}$ in the spherical case, and where Z is the enstrophy, U the r.m.s. velocity, and a the radius of the sphere. This result may help to explain why numerical simulations of unforced beta-plane turbulence (in which ε decreases in time) seem to evolve into a non-ergodic regime at large scales.

Comments on 'Balance and the slow quasimanifold: some explicit results'

The concept of slow vortical dynamics and its role in theoretical understanding is central to geophysical fluid dynamics. It leads, for example, to “potential vorticity thinking” (Hoskins et al. 1985). Mathematically, one imagines an invariant manifold within the phase space of solutions, called the slow manifold (Leith 1980; Lorenz 1980), to which the dynamics are constrained. Whether this slow manifold truly exists has been a major subject of inquiry over the past 20 years. It has become clear that an exact slow manifold is an exceptional case, restricted to steady or perhaps temporally periodic flows (Warn 1997). Thus the concept of a “fuzzy slow manifold” (Warn and Ménard 1986) has been suggested. The idea is that nearly slow dynamics will occur in a stochastic layer about the putative slow manifold. The natural question then is, how thick is this layer? In a recent paper, Ford et al. (2000) argue that Lighthill emission—the spontaneous emission of freely propagating acoustic waves by unsteady vortical flows—is applicable to the problem of balance, with the Mach number Ma replaced by the Froude number F, and that it is a fundamental mechanism for this fuzziness. They consider the rotating shallow-water equations and find emission of inertia–gravity waves at O(F2). This is rather surprising at first sight, because several studies of balanced dynamics with the rotating shallow-water equations have gone beyond second order in F, and found only an exponentially small unbalanced component (Warn and Ménard 1986; Lorenz and Krishnamurthy 1987; Bokhove and Shepherd 1996; Wirosoetisno and Shepherd 2000). We have no technical objection to the analysis of Ford et al. (2000), but wish to point out that it depends crucially on R 1, where R is the Rossby number. This condition requires the ratio of the characteristic length scale of the flow L to the Rossby deformation radius LR to go to zero in the limit F → 0. This is the low Froude number scaling of Charney (1963), which, while originally designed for the Tropics, has been argued to be also relevant to mesoscale dynamics (Riley et al. 1981). If L/LR is fixed, however, then F → 0 implies R → 0, which is the standard quasigeostrophic scaling of Charney (1948; see, e.g., Pedlosky 1987). In this limit there is reason to expect the fuzziness of the slow manifold to be “exponentially thin,” and balance to be much more accurate than is consistent with (algebraic) Lighthill emission.

Implications of model error for numerical climate prediction

Numerical climate models constitute the best available tools to tackle the problem of climate prediction. Two assumptions lie at the heart of their suitability: (1) a climate attractor exists, and (2) the numerical climate model's attractor lies on the actual climate attractor, or at least on the projection of the climate attractor on the model's phase space. In this contribution, the Lorenz '63 system is used both as a prototype system and as an imperfect model to investigate the implications of the second assumption. By comparing results drawn from the Lorenz '63 system and from numerical weather and climate models, the implications of using imperfect models for the prediction of weather and climate are discussed. It is shown that the imperfect model's orbit and the system's orbit are essentially different, purely due to model error and not to sensitivity to initial conditions. Furthermore, if a model is a perfect model, then the attractor, reconstructed by sampling a collection of initialised model orbits (forecast orbits), will be invariant to forecast lead time. This conclusion provides an alternative method for the assessment of climate models.

Reconnection at the high-latitude magnetopause during northward interplanetary magnetic field conditions

The Polar spacecraft had a prolonged encounter with the high-latitude dayside magnetopause on May 29, 1996. This encounter with the magnetopause occurred when the interplanetary magnetic field was directed northward. From the three-dimensional electron and ion distribution functions measured by the Hydra instrument, it has been possible to identify nearly all of the distinct boundary layer regions associated with high-latitude reconnection. The regions that have been identified are (1) the cusp; (2) the magnetopause current layer; (3) magnetosheath field lines that have interconnected in only the Northern Hemisphere; (4) magnetosheath field lines that have interconnected in only the Southern Hemisphere; (5) magnetosheath field lines that have interconnected in both the Northern and Southern Hemispheres; (6) magnetosheath that is disconnected from the terrestrial magnetic field; and (7) high-latitude plasma sheet field lines that are participating in magnetosheath reconnection. Reconnection over this time period was occurring at high latitudes over a broad local-time extent, interconnecting the magnetosheath and lobe and/or plasma sheet field lines in both the Northern and Southern Hemispheres. Newly closed boundary layer field lines were observed as reconnection occur-red first at high latitudes in one hemisphere and then later in the other. These observations establish the location of magnetopause reconnection during these northward interplanetary magnetic field conditions as being at high latitudes, poleward of the cusp, and further reinforce the general interpretation of electron and ion phase space density signatures as indicators of magnetic reconnection and boundary layer formation.

The problem of two fixed centers: bifurcation diagram for positive energies

We give a comprehensive analysis of the Euler-Jacobi problem of motion in the field of two fixed centers with arbitrary relative strength and for positive values of the energy. These systems represent nontrivial examples of integrable dynamics and are analysed from the point of view of the energy-momentum mapping from the phase space to the space of the integration constants. In this setting, we describe the structure of the scattering trajectories in phase space and derive an explicit description of the bifurcation diagram, i.e., the set of critical value of the energy-momentum map.

Advancements in decadal climate predictability: the role of nonoceanic drivers

We review recent progress in understanding the role of sea ice, land surface, stratosphere, and aerosols in decadal-scale predictability and discuss the perspectives for improving the predictive capabilities of current Earth system models (ESMs). These constituents have received relatively little attention because their contribution to the slow climatic manifold is controversial in comparison to that of the large heat capacity of the oceans. Furthermore, their initialization as well as their representation in state-of-the-art climate models remains a challenge. Numerous extraoceanic processes that could be active over the decadal range are proposed. Potential predictability associated with the aforementioned, poorly represented, and scarcely observed constituents of the climate system has been primarily inspected through numerical simulations performed under idealized experimental settings. The impact, however, on practical decadal predictions, conducted with realistically initialized full-fledged climate models, is still largely unexploited. Enhancing initial-value predictability through an improved model initialization appears to be a viable option for land surface, sea ice, and, marginally, the stratosphere. Similarly, capturing future aerosol emission storylines might lead to an improved representation of both global and regional short-term climatic changes. In addition to these factors, a key role on the overall predictive ability of ESMs is expected to be played by an accurate representation of processes associated with specific components of the climate system. These act as “signal carriers,” transferring across the climatic phase space the information associated with the initial state and boundary forcings, and dynamically bridging different (otherwise unconnected) subsystems. Through this mechanism, Earth system components trigger low-frequency variability modes, thus extending the predictability beyond the seasonal scale.

Dynamic portrait of the planetary 2/1 mean-motion resonance - I. Systems with a more massive outer planet

We analyse the global structure of the phase space of the planar planetary 2/1 mean-motion resonance in cases where the outer planet is more massive than its inner companion. Inside the resonant domain, we show the existence of two families of periodic orbits, one associated to the librational motion of resonant angle (sigma-family) and the other related to the circulatory motion of the difference in longitudes of pericentre (Delta pi-family). The well-known apsidal corotation resonances (ACR) appear as intersections between both families. A complex web of secondary resonances is also detected for low eccentricities, whose strengths and positions are dependent on the individual masses and spatial scale of the system. The construction of dynamical maps for various values of the total angular momentum shows the evolution of the families of stable motion with the eccentricities, identifying possible configurations suitable for exoplanetary systems. For low-moderate eccentricities, several different stable modes exist outside the ACR. For larger eccentricities, however, all stable solutions are associated to oscillations around the stationary solutions. Finally, we present a possible link between these stable families and the process of resonance capture, identifying the most probable routes from the secular region to the resonant domain, and discussing how the final resonant configuration may be affected by the extension of the chaotic layer around the resonance region.

Dynamic portrait of the planetary 2/1 mean-motion resonance - II. Systems with a more massive inner planet

This paper presents the second part in our study of the global structure of the planar phase space of the planetary three-body problem, when both planets lie in the vicinity of a 2/1 mean-motion resonance. While Paper I was devoted to cases where the outer planet is the more massive body, the present work is devoted to the cases where the more massive body is the inner planet. As before, outside the well-known Apsidal Corotation Resonances (ACR), the phase space shows a complex picture marked by the presence of several distinct regimes of resonant and non-resonant motion, crossed by families of periodic orbits and separated by chaotic zones. When the chosen values of the integrals of motion lead to symmetric ACR, the global dynamics are generally similar to the structure presented in Paper I. However, for asymmetric ACR the resonant phase space is strikingly different and shows a galore of distinct dynamical states. This structure is shown with the help of dynamical maps constructed on two different representative planes, one centred on the unstable symmetric ACR and the other on the stable asymmetric equilibrium solution. Although the study described in the work may be applied to any mass ratio, we present a detailed analysis for mass values similar to the Jupiter-Saturn case. Results give a global view of the different dynamical states available to resonant planets with these characteristics. Some of these dynamical paths could have marked the evolution of the giant planets of our Solar system, assuming they suffered a temporary capture in the 2/1 resonance during the latest stages of the formation of our Solar system.

Semiclassical evolution of dissipative Markovian systems

A semiclassical approximation for an evolving density operator, driven by a `closed` Hamiltonian operator and `open` Markovian Lindblad operators, is obtained. The theory is based on the chord function, i.e. the Fourier transform of the Wigner function. It reduces to an exact solution of the Lindblad master equation if the Hamiltonian operator is a quadratic function and the Lindblad operators are linear functions of positions and momenta. Initially, the semiclassical formulae for the case of Hermitian Lindblad operators are reinterpreted in terms of a (real) double phase space, generated by an appropriate classical double Hamiltonian. An extra `open` term is added to the double Hamiltonian by the non-Hermitian part of the Lindblad operators in the general case of dissipative Markovian evolution. The particular case of generic Hamiltonian operators, but linear dissipative Lindblad operators, is studied in more detail. A Liouville-type equivariance still holds for the corresponding classical evolution in double phase space, but the centre subspace, which supports the Wigner function, is compressed, along with expansion of its conjugate subspace, which supports the chord function. Decoherence narrows the relevant region of double phase space to the neighbourhood of a caustic for both the Wigner function and the chord function. This difficulty is avoided by a propagator in a mixed representation, so that a further `small-chord` approximation leads to a simple generalization of the quadratic theory for evolving Wigner functions.

Blowout bifurcation and spatial mode excitation in the bubbling transition to turbulence

The transition to turbulence (spatio-temporal chaos) in a wide class of spatially extended dynamical system is due to the loss of transversal stability of a chaotic attractor lying on a homogeneous manifold (in the Fourier phase space of the system) causing spatial mode excitation Since the latter manifests as intermittent spikes this has been called a bubbling transition We present numerical evidences that this transition occurs due to the so called blowout bifurcation whereby the attractor as a whole loses transversal stability and becomes a chaotic saddle We used a nonlinear three-wave interacting model with spatial diffusion as an example of this transition (C) 2010 Elsevier B V All rights reserved

Robust tori in a double-waved Hamiltonian model

A Hamiltonian system perturbed by two waves with particular wave numbers can present robust tori, which are barriers created by the vanishing of the perturbed Hamiltonian at some defined positions. When robust tori exist, any trajectory in phase space passing close to them is blocked by emergent invariant curves that prevent the chaotic transport. Our results indicate that the considered particular solution for the two waves Hamiltonian model shows plenty of robust tori blocking radial transport. (C) 2010 Elsevier B.V. All rights reserved.

The non-twist standard map with robust tori

The non-twist standard map occurs frequently in many fields of science specially in modelling the dynamics of the magnetic field lines in tokamaks. Robust tori, dynamical barriers that impede the radial transport among different regions of the phase space, are introduced in the non-twist standard map in a conservative fashion. The resulting non-twist standard map with robust tori is an improved model to study transport barriers in plasmas confined in tokamaks.

Clustering and diffusion in a symplectic map lattice with non-local coupling

We study a symplectic chain with a non-local form of coupling by means of a standard map lattice where the interaction strength decreases with the lattice distance as a power-law, in Such a way that one can pass continuously from a local (nearest-neighbor) to a global (mean-field) type of coupling. We investigate the formation of map clusters, or spatially coherent structures generated by the system dynamics. Such clusters are found to be related to stickiness of chaotic phase-space trajectories near periodic island remnants, and also to the behavior of the diffusion coefficient. An approximate two-dimensional map is derived to explain some of the features of this connection. (C) 2008 Elsevier Ltd. All rights reserved.

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.