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 threemode) models, and numerically, using direct simulations of CHM equation in the Fourier space. The linear theory [1] 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 [2] 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 zonalaveraged 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).

A detailed study of the Charney-Hasegawa-Mima model and its extensions is presented. These simplest nonlinear partial differential equations suggested for both Rossby waves in the atmosphere and also drift waves in a magnetically-confined plasma exhibits some remarkable and nontrivial properties, which in their qualitative form survive in more realistic and complicated models, and as such form a conceptual basis for understanding the turbulence and zonal flow dynamics in real plasma and geophysical systems. Two idealised scenarios of generation of zonal flows by small-scale turbulence are explored: a modulational instability and turbulent cascades. A detailed study of the generation of zonal flows by the modulational instability reveals that the dynamics of this zonal flow generation mechanism differ widely depending on the initial degree of nonlinearity. The jets in the strongly nonlinear case further roll up into Kármán-like vortex streets and saturate, while for the weaker nonlinearities, the growth of the unstable mode reverses and the system oscillates between a dominant jet, which is slightly inclined to the zonal direction, and a dominant primary wave. A numerical proof is provided for an the extra invariant in Rossby and drift wave turbulence -zonostrophy. While the theoretical derivations of this invariant stem from the wave kinetic equation which assumes weak wave amplitudes, it is shown to be relatively well-conserved for higher nonlinearities also. Together with the energy and enstrophy, these three invariants cascade into anisotropic sectors in the k-space as predicted by the Fjørtoft argument. The cascades are characterised by the zonostrophy pushing the energy to the zonal scales. A small scale instability forcing applied to the model has demonstrated the well-known drift wavezonal flow feedback loop. The drift wave turbulence is generated from this primary instability. The zonal flows are then excited by either one of the generation mechanisms, extracting energy from the drift waves as they grow. Eventually the turbulence is completely suppressed and the zonal flows saturate. The turbulence spectrum is shown to diffuse in a manner which has been mathematically predicted. The insights gained from this simple model could provide a basis for equivalent studies in more sophisticated plasma and geophysical fluid dynamics models in an effort to fully understand the zonal flow generation, the turbulent transport suppression and the zonal flow saturation processes in both the plasma and geophysical contexts as well as other wave and turbulence systems where order evolves from chaos.

We study quasi-geostrophic turbulence and plasma drift turbulence within the Charney-Hasegawa-Mima (CHM) model. We focus, theoretically and using numerical simulations, on conservation of zonostrophy and on its role in the formation of the zonal jets. The zonostrophy invariant was first predicted in [1,2] in two special cases -large-scale turbulence and anisotropic turbulence. Papers [1, 2] also predicted that the three invariants, energy, enstrophy and zonostrophy, will cascade anisotropically into non-intersecting sectors in the k-space, so that the energy cascade is "pushed" into the large-scale zonal scales. In the present paper, we consider the scales much less than the Rossby deformation radius and generalise the Fjørtoft argument of [1, 2] to find the directions of the three cascades in this case. For the first time, we demonstrate numerically that zonostrophy is well conserved by the CHM model, and that the energy, enstrophy and zonostrophy cascade as prescribed by the Fjørtoft argument if the nonlinearity is sufficiently weak. Moreover, numerically we observe that zonostrophy is conserved surprisingly well at late times and the triple-cascade picture is rather accurate even if the initial nonlinearity is strong.