Coherent jets with most of the kinetic energy of the flow are common in atmospheric turbulence. In the gaseous planets these jets are maintained by incoherent turbulence excited by small-scale convection. Large-scale coherent waves are sometimes observed to coexist with the jets; a prominent example is Saturn's hexagonal North polar jet (NPJ). The mechanism responsible for forming and maintaining such a turbulent state remains elusive. The coherent planetary-scale component of the turbulence arises and is maintained by interaction with the incoherent small-scale turbulence component. Theoretical understanding of the dynamics of the jet/wave/turbulence coexistence regime is gained by employing a statistical state dynamics (SSD) model. Here, a second-order closure implementation of a two-layer beta-plane SSD is used to develop a theory that accounts for the structure and dynamics of the NPJ. Asymptotic analysis of the SSD equilibrium in the weak jet damping limit predicts a universal jet structure in agreement with NPJ observations. This asymptotic theory also predicts the wavenumber (six) of the prominent jet perturbation. Analysis with this model of the jet/wave/turbulence regime dynamics reveals that jet formation is controlled by the effective value of $\beta$; the required value of this parameter for correspondence with observation is obtained. As this is a robust prediction it is taken as an indirect observation of a deep poleward sloping stable layer beneath the NPJ. The slope required is obtained from observations of NPJ structure as is the small-scale turbulence excitation required to maintain the jet. The observed jet structure is then predicted by the theory as is the wave-six disturbance. This wave, which is identified with the least stable mode of the equilibrated jet, is shown to be primarily responsible for equilibrating the jet with the observed structure and amplitude.

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.