Magnetic induction in magnetohydrodynamic fluids at magnetic Reynolds number (Rm) less than~1 has long been known to cause magnetic drag. Here, we show that when $\mathrm{Rm} \gg 1$ and the fluid is in a hydrodynamic-dominated regime in which the magnetic energy is much smaller than the kinetic energy, induction due to a mean shear flow leads to a magnetic eddy viscosity. The magnetic viscosity is derived from simple physical arguments, where a coherent response due to shear flow builds up in the magnetic field until decorrelated by turbulent motion. The dynamic viscosity coefficient is approximately $(B_p^2/2\mu_0) \tau_{\rm corr}$, the poloidal magnetic energy density multiplied by the correlation time. We confirm the magnetic eddy viscosity through numerical simulations of two-dimensional incompressible magnetohydrodynamics. We also consider the three-dimensional case, and in cylindrical or spherical geometry, theoretical considerations similarly point to a nonzero viscosity whenever there is differential rotation. Hence, these results serve as a dynamical generalization of Ferraro's law of isorotation. The magnetic eddy viscosity leads to transport of angular momentum and may be of importance to zonal flows in astrophysical domains such as the interior of some gas giants.

Zonal flows in rotating systems have been previously shown to be suppressed by the imposition of a background magnetic field aligned with the direction of rotation. Understanding the physics behind the suppression may be important in systems found in astrophysical fluid dynamics, such as stellar interiors. However, the mechanism of suppression has not yet been explained. In the idealized setting of a magnetized beta plane, we provide a theoretical explanation that shows how magnetic fluctuations directly counteract the growth of weak zonal flows. Two distinct calculations yield consistent conclusions. The first, which is simpler and more physically transparent, extends the Kelvin-Orr shearing wave to include magnetic fields and shows that weak, long-wavelength shear flow organizes magnetic fluctuations to absorb energy from the mean flow. The second calculation, based on the quasilinear, statistical CE2 framework, is valid for arbitrary wavelength zonal flow and predicts a self-consistent growth rate of the zonal flow. We find that a background magnetic field suppresses zonal flow if the bare Alfven frequency is comparable to or larger than the bare Rossby frequency. However, suppression can occur for even smaller magnetic fields if the resistivity is sufficiently small enough to allow sizable magnetic fluctuations. Our calculations reproduce the $\eta/B_0^2 = \text{const.}$ scaling that describes the boundary of zonation, as found in previous work, and we explicitly link this scaling to the amplitude of magnetic fluctuations.

Zonal jets and non-zonal large-scale flows are often present in forced-dissipative barotropic turbulence on a beta-plane. The dynamics underlying the formation of both zonal and non-zonal coherent structures is investigated in this work within the statistical framework of Stochastic Structural Stability Theory (S3T). Previous S3T studies have shown that the homogeneous turbulent state undergoes a bifurcation at a critical parameter and becomes inhomogeneous with the emergence of zonal and/or large-scale non-zonal flows and that these statistical predictions of S3T are reflected in direct numerical simulations. In this paper, we study the dynamics underlying the S3T statistical instability of the homogeneous state as a function of parameters. It is shown that for weak planetary vorticity gradient, $\beta$, both zonal jets and non-zonal large-scale structures form from upgradient momentum fluxes due to shearing of the eddies by the emerging infinitesimal flow. For large $\beta$, the dynamics of the S3T instability differs for zonal and non-zonal flows but in both the destabilizing vorticity fluxes decrease with increasing $\beta$. Shearing of the eddies by the mean flow continues to be the mechanism for the emergence of zonal jets while non-zonal large-scale flows emerge from resonant and near resonant triad interactions between the large-scale flow and the stochastically forced eddies. The relation between the formation of large-scale structure through modulational instability and the S3T instability of the homogeneous state is also investigated and it is shown that the modulational instability results are subsumed by the S3T results.

Stochastic Structural Stability Theory (S3T) provides analytical methods for understanding the emergence and equilibration of jets from the turbulence in planetary atmospheres based on the dynamics of the statistical mean state of the turbulence closed at second order. Predictions for formation and equilibration of turbulent jets made using S3T are critically compared with results of simulations made using the associated quasi-linear and nonlinear models. S3T predicts the observed bifurcation behavior associated with the emergence of jets, their equilibration and their breakdown as a function of parameters. Quantitative differences in bifurcation parameter values between predictions of S3T and results of nonlinear simulations are traced to modification of the eddy spectrum which results from two processes: nonlinear eddy-eddy interactions and formation of discrete non-zonal structures. Remarkably, these non-zonal structures, which substantially modify the turbulence spectrum, are found to arise from S3T instability. Formation as linear instabilities and equilibration at finite amplitude of multiple equilibria for identical parameter values in the form of jets with distinct meridional wavenumbers is verified as is the existence at equilibrium of finite amplitude non-zonal structures in the form of nonlinearly modified Rossby waves. When zonal jets and nonlinearly modified Rossby waves coexist at finite amplitude the jet structure is generally found to dominate even if it is linearly less unstable. The physical reality of the manifold of S3T jets and non-zonal structures is underscored by the existence in nonlinear simulations of jet structure at subcritical S3T parameter values which are identified with stable S3T jet modes excited by turbulent fluctuations.

Highly stratified shear layers are rendered unstable even at high stratifications by Holmboe instabilities when the density stratification is concentrated in a small region of the shear layer. These instabilities may cause mixing in highly stratified environments. However these instabilities occur in limited bands in the parameter space. We perform Generalized Stability analysis of the two dimensional perturbation dynamics of an inviscid Boussinesq stratified shear layer and show that Holmboe instabilities at high Richardson numbers can be excited by their adjoints at amplitudes that are orders of magnitude larger than by introducing initially the unstable mode itself. We also determine the optimal growth that is obtained for parameters for which there is no instability. We find that there is potential for large transient growth regardless of whether the background flow is exponentially stable or not and that the characteristic structure of the Holmboe instability asymptotically emerges as a persistent quasi-mode for parameter values for which the flow is stable.