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.

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.