Modulational instability of Rossby and drift waves and generation of zonal jets

We study the modulational instability of geophysical Rossby and plasma drift waves within the Charney-Hasegawa-Mima (CHM) model both theoretically, using truncated (four-mode and three-mode) models, and numerically, using direct simulations of CHM equation in the Fourier space. The linear theory predicts instability for any amplitude of the primary wave. For strong primary waves the most unstable modes are perpendicular to the primary wave, which correspond to generation of a zonal flow if the primary wave is purely meridional. For weaker waves, the maximum growth occurs for off-zonal inclined modulations. For very weak primary waves the unstable waves are close to being in three-wave resonance with the primary wave. The nonlinear theory predicts that the zonal flows generated by the linear instability experience pinching into narrow zonal jets. Our numerical simulations confirm the theoretical predictions of the linear theory as well as of the nonlinear pinching. We find that, for strong primary waves, these narrow zonal jets further roll up into Karman-like vortex streets. On the other hand, for weak primary waves, the growth of the unstable mode reverses and the system oscillates between a dominant jet and a dominate primary wave. The 2D vortex streets appear to be more stable than purely 1D zonal jets, and their zonal-averaged speed can reach amplitudes much stronger than is allowed by the Rayleigh-Kuo instability criterion for the 1D case. We find that the truncation models work well for both the linear stage and and often even for the medium-term nonlinear behavior. In the long term, the system transitions to turbulence helped by the vortex-pairing instability (for strong waves) and by the resonant wave-wave interactions (for weak waves).