Andrew D. Gilbert

The shallow-water equations are widely used to model interactions between horizontal shear flows and (rotating) gravity waves in thin planetary atmospheres. Their extension to allow for interactions with magnetic fields – the equations of shallow-water magnetohydrodynamics (SWMHD) – is often used to model waves and instabilities in thin stratified layers in stellar and planetary atmospheres, in the perfectly conducting limit. Here we consider how magnetic diffusion should be added to the equations of SWMHD. This is crucial for an accurate balance between advection and diffusion in the induction equation, and hence for modelling instabilities and turbulence. For the straightforward choice of Laplacian diffusion, we explain how fundamental mathematical and physical inconsistencies arise in the equations of SWMHD, and show that unphysical dynamo action can result. We then derive a physically consistent magnetic diffusion term by performing an asymptotic analysis of the three-dimensional equations of magnetohydrodynamics in the thin-layer limit, giving the resulting diffusion term explicitly in both planar and spherical coordinates. We show how this magnetic diffusion term, which allows for a horizontally varying diffusivity, is consistent with the standard shallow-water solenoidal constraint, and leads to negative semidefinite Ohmic dissipation. We also establish a basic type of antidynamo theorem.

Zonal flows are mean flows in the east–west direction, which are ubiquitous on planets, and can be formed through ‘zonostrophic instability’: within turbulence or random waves, a weak large-scale zonal flow can grow exponentially to become prominent. In this paper, we study the statistical behaviour of the zonostrophic instability and the effect of magnetic fields. We use a stochastic white noise forcing to drive random waves, and study the growth of a mean flow in this random system. The dispersion relation for the growth rate of the expectation of the mean flow is derived, and properties of the instability are discussed. In the limits of weak and strong magnetic diffusivity, the dispersion relation reduces to manageable expressions, which provide clear insights into the effect of the magnetic field and scaling laws for the threshold of instability. The magnetic field mainly plays a stabilising role and thus impedes the formation of the zonal flow, but under certain conditions it can also have destabilising effects. Numerical simulation of the stochastic flow is performed to confirm the theory. Results indicate that the magnetic field can significantly increase the randomness of the zonal flow. It is found that the zonal flow of an individual realisation may behave very differently from the expectation. For weak magnetic diffusivity and moderate magnetic field strengths, this leads to considerable variation of the outcome, that is whether zonostrophic instability takes place or not in individual realisations.

This paper studies the instability of two-dimensional magnetohydrodynamic systems on a sphere using analytical methods. The underlying flow consists of a zonal differential rotation and a toroidal magnetic field is present. Semicircle rules that prescribe the possible domain of the wave velocity in the complex plane for general flow and field profiles are derived. The paper then sets out an analytical study of the ‘clamshell instability’, which features field lines on the two hemispheres tilting in opposite directions (Cally, Sol. Phys., vol. 199, 2001, pp. 231–249). An asymptotic solution for the instability problem is derived for the limit of weak shear of the zonal flow, via the method of matched asymptotic expansions. It is shown that when the zonal flow is solid body rotation, there exists a neutral mode that tilts the magnetic field lines, referred to as the ‘tilting mode’. A weak shear of the zonal flow excites the critical layer of the tilting mode, which reverses the tilting direction to form the clamshell pattern and induces the instability. The asymptotic solution provides insights into properties of the instability for a range of flow and field profiles. A remarkable feature is that the magnetic field affects the instability only through its local behaviour in the critical layer.

In this paper, the instability of shallow-water shear flow with a sheared parallel magnetic field is studied. Waves propagating in such magnetic shear flows encounter critical levels where the phase velocity relative to the basic flow, in the critical layer, is found to play a complementary role in the instability.

A numerical investigation is presented of kinematic dynamo action in a dynamically driven fluid flow. The model isolates basic dynamo processes relevant to field generation in the Solar tachocline. The horizontal plane layer geometry adopted is chosen as the local representation of a differentially rotating spherical fluid shell at co-latitude ϑ; the unit vectors [xcirc ], ŷ and zˆ point east, north and vertically upwards respectively. Relative to axes moving easterly with the local bulk motion of the fluid the rotation vector Ω lies in the (y, z)-plane inclined at an angle ϑ to the z-axis, while the base of the layer moves with constant velocity in the x-direction. An Ekman layer is formed on the lower boundary characterized by a strong localized spiralling shear flow. This basic state is destabilized by a convective instability through uniform heating at the base of the layer, or by a purely hydrodynamic instability of the Ekman layer shear flow. The onset of instability is characterized by a horizontal wave vector inclined at some angle ∈ to the x-axis. Such motion is two-dimensional, dependent only on two spatial coordinates together with time. It is supposed that this two-dimensionality persists into the various fully nonlinear regimes in which we study large magnetic Reynolds number kinematic dynamo action.When the Ekman layer flow is destabilized hydrodynamically, the fluid flow that results is steady in an appropriately chosen moving frame, and takes the form of a row of cat's eyes. Kinematic magnetic field growth is characterized by modes of two types. One is akin to the Ponomarenko dynamo mechanism and located close to some closed stream surface; the other appears to be associated with stagnation points and heteroclinic separatrices.When the Ekman layer flow is destabilized thermally, the well-developed convective instability far from onset is characterized by a flow that is intrinsically time-dependent in the sense that it is unsteady in any moving frame. The magnetic field is concentrated in magnetic sheets situated around the convective cells in regions where chaotic particle paths are likely to exist; evidence for fast dynamo action is obtained. The presence of the Ekman layer close to the bottom boundary breaks the up–down symmetry of the layer and localizes the magnetic field near the lower boundary.