On the mechanism of the Gent–McWilliams instability of a columnar vortex in stratified rotating fluids

In stably stratified and rotating fluids, an axisymmetric columnar vortex can be unstable to a special instability with an azimuthal wavenumber $m=1$ which bends and slices the vortex into pancake vortices (Gent & McWilliams Geophys. Astrophys. Fluid Dyn., vol. 35 (1–4), 1986, pp. 209–233). This bending instability, called the ‘Gent–McWilliams instability’ herein, is distinct from the shear, centrifugal or radiative instabilities. The goals of the paper are to better understand the origin and properties of this instability and to explain why it operates only in stratified rotating fluids. Both numerical and asymptotic stability analyses of several velocity profiles have been performed for wide ranges of Froude number $\mathit{Fr}_{h}={\it\Omega}_{0}/N$ and Rossby number $\mathit{Ro}=2{\it\Omega}_{0}/f$ , where ${\it\Omega}_{0}$ is the angular velocity on the vortex axis, $N$ the Brunt–Väisälä frequency and $f$ the Coriolis parameter. Numerical analyses restricted to the centrifugally stable range show that the maximum growth rate of the Gent–McWilliams instability increases with $\mathit{Ro}$ and is independent of $\mathit{Fr}_{h}$ for $\mathit{Fr}_{h}\leqslant 1$ . In contrast, when $\mathit{Fr}_{h}>1$ , the maximum growth rate decreases dramatically with $\mathit{Fr}_{h}$ . Long axial wavelength asymptotic analyses for isolated vortices prove that the Gent–McWilliams instability is due to the destabilization of the long-wavelength bending mode by a critical layer at the radius $r_{c}$ where the angular velocity ${\it\Omega}$ is equal to the frequency ${\it\omega}$ : ${\it\Omega}(r_{c})={\it\omega}$ . A necessary and sufficient instability condition valid for long wavelengths, finite Rossby number and $\mathit{Fr}_{h}\leqslant 1$ is that the derivative of the vertical vorticity of the basic vortex is positive at $r_{c}$ : ${\it\zeta}^{\prime }(r_{c})>0$ . Such a critical layer $r_{c}$ exists for finite Rossby and Froude numbers because the real part of the frequency of the long-wavelength bending mode is positive instead of being negative as in a homogeneous non-rotating fluid ( $\mathit{Ro}=\mathit{Fr}_{h}=\infty$ ). When $\mathit{Fr}_{h}>1$ , the instability condition ${\it\zeta}^{\prime }(r_{c})>0$ is necessary but not sufficient because the destabilizing effect of the critical layer $r_{c}$ is strongly reduced by a second stabilizing critical layer $r_{c2}$ existing at the radius where the angular velocity is equal to the Brunt–Väisälä frequency. For