Classification of matter through topological phases and topological edge states between distinct materials has been a subject of great interest recently. While lattices have been the main setting for these studies, a relatively unexplored realm for this physics is that of continuum fluids. In the typical case of a fluid model with a point spectrum, nontrivial topology and associated edge modes have been observed previously. However, another possibility is that a continuous spectrum can conexist with the point spectrum. Here we demonstrate that a fluid plasma model can harbor nontrivial topology within its continuous spectrum, and that there are boundary modes at the interface between topologically distinct regions. We consider the ideal magnetohydrodynamics (MHD) model. In the presence of magnetic shear, we find nontrivial topology in the {\Alfven} continuum. For strong shear, the Chern number is $\pm 1$, depending on the sign of the shear. If the magnetic shear changes sign within the plasma, a topological phase transition occurs, and bulk-boundary correspondence then suggests a mode localized to the layer of zero magnetic shear. We confirm the existence of this mode numerically. Moreover, this reversed-shear {\Alfven} eigenmode (RSAE) is well known within the magnetic fusion community as it has been observed in several tokamaks. A potential implication of our findings is that the RSAE may be topologically robust to three-dimensional perturbations such as magnetic islands. Moreover, this work raises the broader possibility that fusion devices could be carefully tailored to produce other topological states and potentially useful behavior. In considering how the MHD model might be regularized at small scales, we also consider the electron limit of Hall MHD. We show that the whistler band, which plays an important role in planetary ionospheres, has nontrivial topology.

Direct Statistical Simulation (DSS) solves the equations of motion for the statistics of turbulent flows in place of the traditional route of accumulating statistics by Direct Numerical Simulation (DNS). That low-order statistics usually evolve slowly compared with instantaneous dynamics is one important advantage of DSS. Depending on the symmetry of the problem and the choice of averaging operation, however, DSS is usually more expensive computationally than DNS because even low order statistics typically have higher dimension than the underlying fields. Here we show that it is possible to go much further by using Proper Orthogonal Decomposition (POD) to address the "curse of dimensionality." We apply POD directly to DSS in the form of expansions in the equal-time cumulants to second order (CE2). We explore two averaging operations (zonal and ensemble) and test the approach on two idealized barotropic models on a rotating sphere (a jet that relaxes deterministically towards an unstable profile, and a stochastically-driven flow that spontaneously organizes into jets). Order-of-magnitude savings in computational cost are obtained in the reduced basis, potentially enabling access to parameter regimes beyond the reach of DNS.

Global magnetohydrodynamic (MHD) instabilities are investigated in a computationally tractable two-dimensional model of the solar tachocline. The model's differential rotation yields stability in the absence of a magnetic field, but if a magnetic field is present, a joint instability is observed. We analyze the nonlinear development of the instability via fully nonlinear direct numerical simulation, the generalized quasilinear approximation (GQL), and direct statistical simulation (DSS) based upon low-order expansion in equal-time cumulants. As the magnetic diffusivity is decreased, the nonlinear development of the instability becomes more complicated until eventually a set of parameters are identified that produce a previously unidentified long-term cycle in which energy is transformed from kinetic energy to magnetic energy and back. We find that the periodic transitions, which mimic some aspects of solar variability -- for example, the quasiperiodic seasonal exchange of energy between toroidal field and waves or eddies -- are unable to be reproduced when eddy-scattering processes are excluded from the model.

In this paper we examine the interaction of convection, rotation and mean flows in a thermal annulus. In this system mean flows are driven by correlations induced by rotation leading to non-trivial Reynolds stresses. The mean flows act back on the convective turbulence acting as a barrier to transport. For this system we demonstrate that the Generalised Quasilinear Approximation (GQL) (Marston et al 2016) may provide a much better approximation to the complicated full nonlinear dynamics than the widely used Quasilinear Approximation (QL). This result will enable the construction of more accurate statistical theories for the description of geophysical and astrophysical flows.

Topology sheds new light on the emergence of unidirectional edge waves in a variety of physical systems, from condensed matter to artificial lattices. Waves observed in geophysical flows are also robust to perturbations, which suggests a role for topology. We show a topological origin for two celebrated equatorially trapped waves known as Kelvin and Yanai modes, due to the Earth's rotation that breaks time-reversal symmetry. The non-trivial structure of the bulk Poincar\'e wave modes encoded through the first Chern number of value $2$ guarantees existence for these waves. This invariant demonstrates that ocean and atmospheric waves share fundamental properties with topological insulators, and that topology plays an unexpected role in the Earth climate system.

We review progress that has been made in utilizing one form of Direct Statistical Simulation (DSS) to describe geophysical and astrophysical flows that are anisotropic and inhomogeneous. We first explain the approach, which is based upon a systematic and conservative expansion of the equations of motion for low-order equal-time cumulants. We place the method into context with other statistical procedures. Truncation at second order in the hierarchy of cumulants is equivalent to retaining the interaction between zonal mean flows and eddies. Eddy-eddy interactions appear at higher orders, but care must be taken to keep the higher-order expansions realizable with non-negative probability distribution functions. The strengths and weaknesses of different levels of approximation are assessed with numerical experiments on the fiducial problem of a stochastically forced jet on a spherical surface. The results give an insight into the mechanisms that may control jet spacing and strength, and indicate interesting avenues for future research.

Quasilinear theory is often utilized to approximate the dynamics of fluids exhibiting significant interactions between mean flows and eddies. In this paper we present a generalization of quasilinear theory to include dynamic mode interactions on the large scales. This generalized quasilinear (GQL) approximation is achieved by separating the state variables into large and small zonal scales via a spectral filter rather than by a decomposition into a formal mean and fluctuations. Nonlinear interactions involving only small zonal scales are then removed. The approximation is conservative and allows for scattering of energy between small-scale modes via the large scale (through non-local spectral interactions). We evaluate GQL for the paradigmatic problems of the driving of large-scale jets on a spherical surface and on the beta-plane and show that it is accurate even for a small number of large-scale modes. As this approximation is formally linear in the small zonal scales it allows for the closure of the system and can be utilized in direct statistical simulation schemes that have proved an attractive alternative to direct numerical simulation for many geophysical and astrophysical problems.

In this paper we introduce the concept of Direct Statistical Simulation (DSS) for astrophysical flows. This technique may be appropriate for problems in astrophysical fluids where the instantaneous dynamics of the flows are of secondary importance to their statistical properties. We give examples of such problems including mixing and transport in planets, stars and disks. The method is described for a general set of evolution equations, before we consider the specific case of a spectral method optimised for problems on a spherical surface. The method is illustrated for the simplest non-trivial example of hydrodynamics and MHD on a rotating spherical surface. We then discuss possible extensions of the method both in terms of computational methods and the range of astrophysical problems that are of interest.