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.

A Hamiltonian weak-wave model for shallow-water flow

A reduced dynamical model is derived which describes the interaction of weak inertia–gravity waves with nonlinear vortical motion in the context of rotating shallow–water flow. The formal scaling assumptions are (i) that there is a separation in timescales between the vortical motion and the inertia–gravity waves, and (ii) that the divergence is weak compared to the vorticity. The model is Hamiltonian, and possesses conservation laws analogous to those in the shallow–water equations. Unlike the shallow–water equations, the energy invariant is quadratic. Nonlinear stability theorems are derived for this system, and its linear eigenvalue properties are investigated in the context of some simple basic flows.

DK-1D: a drift-kinetic simulation tool for modelling the shear Alfvén wave and its interaction with collisionless plasma

We present a highly accurate tool for the simulation of shear Alfven waves (SAW) in collisionless plasma. SAW are important in space plasma environments because for small perpendicular scale lengths they can support an electric field parallel to the ambient magnetic field. Electrons can be accelerated by the parallel electric field and these waves have been implicated as the source of vibrant auroral displays. However, the parallel electric field carried by SAW is small in comparison to the perpendicular electric field of the wave, making it difficult to measure directly in the laboratory, or by satellites in the near-Earth plasma environment. In this paper, we present a simulation code that provides a means to study in detail the SAW-particle interaction in both space and laboratory plasma. Using idealised, small-amplitude propagating waves with a single perpendicular wavenumber, the simulation code accurately reproduces the damping rates and parallel electric field amplitudes predicted by linear theory for varying temperatures and perpendicular scale lengths. We present a rigorous kinetic derivation of the parallel electric field strength for small-amplitude SAW and show that commonly-used inertial and kinetic approximations are valid except for where the ratio of thermal to Alfv\'{e}n speed is between 0.7 and 1.0. We also present nonlinear simulations of large-amplitude waves and show that in cases of strong damping, the damping rates and parallel electric field strength deviate from linear predictions when wave energies are greater than only a few percent of the plasma kinetic energy, a situation which is often observed in the magnetosphere. The drift-kinetic code provides reliable, testable predictions of the parallel electric field strength which can be investigated directly in the laboratory, and will help to bridge the gap between studies of SAW in man-made and naturally occuring plasma.

On Rossby waves modified by weak sinusoidal shear

The spatial structure of beta-plane Rossby waves in a sinusoidal basic zonal flow U 0cos(γ,y) is determined analytically in the (stable) asymptotic limit of weak shear, U 0γ2 0/β≈1. The propagating neutral normal modes are found to take their greatest amplitude in the region of maximum westerly flow, while their most rapid phase variation is achieved in the region of maximum easterly flow. These results are shown to be consistent with what is obtained by ray-tracing methods in the limit of small meridional disturbance wavelength.