Rotational and magnetic instability in the diffusive tachocline

Abstract

In this article we study the linear instability of the two-dimensional strongly stratified model for global MHD in the diffusive solar tachocline. Gilman and Fox [Gilman, P.A. and Fox, P., Joint instability of the latitudinal differential rotation and toroidal magnetic fields below the solar convection zone. Astrophys. J., 1997, 484, 439–454] showed that for ideal MHD, the observed surface differential rotation becomes more unstable than is predicted by Watson's [Watson, M., Shear instability of differential rotation in stars. Geophys. Astrophys. Fluid Dyn., 1981, 16, 285–298] nonmagnetic analysis. They showed that the solar differential rotation is unstable for essentially all reasonable values of the differential rotation in the presence of an antisymmetric toroidal field. They found that for the broad field case B φ∼sinθcosθ, θ being the co-latitude, instability occurs only for the azimuthal m = 1 mode, and concluded that modes which are symmetric (meridional flow in the same direction) about the equator onset at lower field strengths than the antisymmetric modes. We study the effect of viscosity and magnetic diffusivity in the strongly stably stratified case where diffusion is primarily along the level surfaces. We show that antisymmetric modes are now strongly preferred over symmetric modes, and that diffusion can sometimes be destabilising. Even solid body rotation can be destabilised through the action of magnetic field. In addition, we show that when diffusion is present, instability can occur when the longitudinal wavenumber m = 2.