The stirring and mixing of heat and momentum in the ocean surface boundary layer (OSBL) are dominated by 1 to 10 km fluid flows – too small to be resolved in global and regional ocean models. Instead, these processes are parametrized. Two main parametrizations include vertical mixing by surface-forced metre-scale turbulence and overturning by kilometre-scale submesoscale frontal flows and instabilities. In present models, these distinct parametrizations are implemented in tandem, yet ignore meaningful interactions between these two scales that may influence net turbulent fluxes. Using a large-eddy simulation of frontal spin down resolving processes at both scales, this work diagnoses submesoscale and surface-forced turbulence impacts that are the foundation of OSBL parametrizations, following a traditional understanding of these flows. It is shown that frontal circulations act to suppress the vertical buoyancy flux by surface forced turbulence, and that this suppression is not represented by traditional boundary layer turbulence theory. A main result of this work is that current OSBL parametrizations excessively mix buoyancy and overestimate turbulence dissipation rates in the presence of lateral flows. These interactions have a direct influence on the upper ocean potential vorticity and energy budgets with implications for global upper ocean budgets and circulation.

The recent work of Siegelman & Young (Proc. Natl Acad. Sci. USA, vol. 120, issue 44, 2023, e2308018120) revealed two extreme states reached by the evolution of unforced and weakly damped two-dimensional turbulence above random rough topography, separated by a critical kinetic energy $E_\#$. The low- and high-energy solutions correspond to topographically locked and roaming vortices, surrounded by non-uniform and homogeneous background potential vorticity (PV), respectively. However, we found that these phenomena are restricted to the particular intermediate length scale where the energy was initially injected into the system. Through simulations initialized by injecting the energy at larger and smaller length scales, we found that the long-term state of topographic turbulence is also dependent on the initial length scale and thus the initial enstrophy. If the initial length scale is comparable to the domain size, the long-term flow field resembles the minimum-enstrophy state proposed by Bretherton & Haidvogel (J. Fluid Mech., vol. 78, issue 1, 1976, pp. 129–154), with very few topographically locked vortices; the long-term enstrophy is quite close to the minimum value, especially when the energy is no larger than $E_\#$. As the initial length scale becomes smaller, more vortices nucleate and become more mobile; the long-term enstrophy increasingly deviates from the minimum value. Simultaneously, the background PV tends to homogenization, even if the energy is below $E_\#$. These results complement the phenomenology of topographic turbulence documented by Siegelman & Young, by showing that the minimum-enstrophy and background PV homogenization states can be adequately approached by large- and small-scale initial fields, respectively, with relatively arbitrary energy.

Three-dimensional direct numerical simulations of rotating Rayleigh–Bénard convection in the planar geometry with no-slip top and bottom and periodic lateral boundary conditions are performed for a broad parameter range with the Rayleigh number spanning in $5\times 10^{6}\leq Ra \leq 5\times 10^{13}$, Ekman number within $5\times 10^{-9}\leq Ek \leq 5\times 10^{-5}$ and Prandtl number $Pr=1$. The thermal and Ekman boundary layer (BL) statistics, temperature drop within the thermal BL, interior temperature gradient and scaling behaviours of the heat and momentum transports (reflected in the Nusselt $Nu$ and Reynolds numbers $Re$) as well as the convective length scale are investigated across various flow regimes. The global and local momentum transports are examined via the $Re$ scaling derived from the classical theoretical balances of viscous–Archimedean–Coriolis (VAC) and Coriolis–inertial–Archimedean (CIA) forces. The VAC-based $Re$ scaling is shown to agree well with the data in the cellular and columnar regimes, where the characteristic convective length scales as the onset length scale ${\sim } Ek^{1/3}$, while the CIA-based $Re$ scaling and the inertia length scale $\sim (ReEk)^{1/2}$ work well in the geostrophic turbulence regime for $Ek\leq 1.5\times 10^{-8}$. The examinations of $Nu$, global and local $Re$, and convective length scale as well as the temperature drop within the thermal BL and its thickness scaling behaviours, indicate that for extreme parameters of $Ek\leq 1.5\times 10^{-8}$ and $80\lesssim RaEk^{4/3}\lesssim 200$, we have reached the diffusion-free geostrophic turbulence regime.

The spontaneous formation of zonal jets is a distinctive feature of geostrophic turbulence with the phenomenon witnessed in numerous numerical studies. In such systems, strong rotation anisotropises the spectral evolution of the energy density such that zonal modes are favoured. In physical space, this manifests as eddies zonally elongating and forming into zonal jets. In the presence of large scale dissipation, the flow may reach statistical stationarity such that the zonal structure persists in the zonal and time mean, and is supported by a flux of eddy momentum. What is unclear is how the excitation of Rossby waves arranges the underlying eddy momentum stresses to support the mean flow structures. To study this, we examine a steady-state flow in the so-called ‘zonostrophic’ regime, in which characteristic scales of geostrophic turbulence are well separated and there are several alternating zonal jets that have formed spontaneously. We apply a geometric eddy ellipse formulation, in which momentum fluxes are cast as ellipses that encode information about the magnitude and direction of flux; the latter is described using the tilt angle. With the aid of a zonal filter, it is revealed that the scales responsible for providing the momentum fluxes associated with the jet structure are much smaller than the characteristic scales identified, and occupy a region of the energy spectrum that has been typically associated with isotropic dynamics.

The parameterization of fluxes associated with representing unresolved dynamics in turbulent flows, especially in the atmosphere and ocean (which have a vast range of scales), remains a challenging task. This is especially true for Earth system models including complex biogeochemistry and requiring very long simulations. The problem of representing the dependence of the mean flux of a passive tracer in terms of the mean has a very long history; in this study, we take a somewhat different approach. We use a formalism showing that the mean flux will be a functional of the mean gradients, a formalism that can be used to calculate the structure of the functional which is non-local in both space and time. Two-dimensional turbulent simulations are used to explore the weight of nearby (in space or time) gradients. We also use stochastic velocities and iterated maps to show that the results are similar. The functional formalism provides an understanding of when non-locality needs to be considered and when a local eddy diffusivity can be a reasonably good approximation. Furthermore, the formalism provides guidance for the development of data-driven parameterizations.

The geostrophic turbulence in rapidly rotating thermal convection exhibits characteristics shared by many highly turbulent geophysical and astrophysical flows. In this regime, the convective length and velocity scales and heat flux are all diffusion-free, i.e. independent of the viscosity and thermal diffusivity. Our direct numerical simulations (DNS) of rotating Rayleigh–Bénard convection in domains with no-slip top and bottom and periodic lateral boundary conditions for a fluid with the Prandtl number $Pr=1$ and extreme buoyancy and rotation parameters (the Rayleigh number up to $Ra=3\times 10^{13}$ and the Ekman number down to $Ek=5\times 10^{-9}$) indeed demonstrate all these diffusion-free scaling relations, in particular, that the dimensionless convective heat transport scales with the supercriticality parameter $\widetilde {Ra}\equiv Ra\, Ek^{4/3}$ as $Nu-1\propto \widetilde {Ra}^{3/2}$, where $Nu$ is the Nusselt number. We further derive and verify in the DNS that with the decreasing $\widetilde {Ra}$, the geostrophic turbulence regime undergoes a transition into another geostrophic regime, the convective heat transport in this regime is characterized by a very steep $\widetilde {Ra}$-dependence, $Nu-1\propto \widetilde {Ra}^{3}$.

There are numerous examples of long-lived, surface-intensified anticyclones over submarine depressions and troughs in the ocean. These often co-exist with a large-scale cyclonic circulation. The latter is predicted by existing barotropic theory but the anticyclone is not. We extend one such theory, which minimizes enstrophy while conserving energy, to two fluid layers. This yields a bottom-intensified flow with cyclonic circulation over a depression. The solution is steady, an enstrophy minimum and stable. When the Lagrange multiplier, $\lambda$, is near zero, the total potential vorticity (PV) becomes homogenized, in both layers. For positive $\lambda$, the surface PV is anticyclonic and strongest at intermediate energies. In quasi-geostrophic numerical simulations with a random initial perturbation PV, the bottom-intensified cyclonic flow always emerges. Vortices evolve independently in the layers and vortex mergers are asymmetric over the depression; cyclones are preferentially strained out at depth while only anticyclones merge at the surface. Both asymmetries are linked to the topographic flow. The deep cyclones feed the bottom-intensified cyclonic circulation while the asymmetry at the surface is only apparent after that circulation has spun up. The result of the surface merger asymmetry is often a lone anticyclone above the depression. This occurs primarily at intermediate energies, when the surface PV predicted by the theory is strongest. Similar results obtain in a full complexity ocean model but with a more pronounced asymmetry in surface vortex mergers and, with bottom friction, significant bottom flow beneath the central anticyclone.

The quasigeostrophic shallow-water system on the mid-latitude $\beta$ plane with weak, small-scale turbulent forcing is explored in the limit of large energy. Forcing is weak in the sense that the energy input rate relative to the energy of the flow is very small, of the order of $10^{-5}$–$10^{-10}$, and the potential vorticity assumes an approximate staircase structure. The flow has large energy in the sense that the jet spacing is equal to the domain width so that no further jet mergers can occur. Quasi-stationary numerical experiments, in which the energy grows linearly, reveal late-time quasi-steady, translating solutions comprising a single jet and vortex dipole, with details of the jet-vortex configuration depending on the deformation radius. At a smaller deformation radius the jet may traverse the entire domain in the $y$ direction one or more times, giving a jet orientation that is predominantly north–south, rather than the usual east–west orientation characteristic of $\beta$-plane jets at lower energy. In these meandering cases, a mode number is proposed that quantifies the degree of meandering relative to the vortices. Besides the steadily translating solutions, topological changes in the jet-vortex structure are identified that occur via a transient interaction of a meandering jet with a vortex. At high energy, these give rise to apparently periodic solutions of the system; at low energy, before a single, domain-wide jet is established, they indicate that jet merger may occur through more complicated processes than the simple merging of neighbouring jets.

The transport induced by ocean mesoscale eddies remains unresolved in most state-of-the-art climate models and needs to be parametrized instead. The natural scale separation between the forcing and the emergent turbulent flow calls for a diffusive parametrization, where the eddy-induced fluxes are related to the large-scale gradients by a diffusion tensor. The standard parametrization scheme in climate modelling consists in adopting the Gent–McWilliams/Redi (GM/R) form for the diffusion tensor, initially put forward based on physical intuition and educated guesses before being put on firm analytical footing using a thickness-weighted average (TWA). In the present contribution, we provide a direct derivation of this diffusion tensor from the quasi-geostrophic (QG) dynamics of a horizontally homogeneous three-dimensional patch of ocean hosting a large-scale vertically sheared zonal flow on the $\beta$ plane. The derivation hinges on the identification of a useful cross-invariant defined as the product of the buoyancy and QG potential vorticity fluctuations. While less general than the TWA approach, the present QG framework leads to rigorous constraints on the diffusion tensor. First, there is no diapycnal diffusivity arising in the QG GM/R tensor for low viscosity and small-scale diffusivities. The diffusion tensor then involves only two vertically dependent coefficients, namely the GM transport coefficient $K_{GM}(z)$ and the Redi diffusivity $K_R(z)$. Second, as identified already by previous authors, the vertical structures of the two coefficients are related by the so-called Taylor–Bretherton relation. Finally, while the two coefficients differ generically in the interior of the water column, we show that they are equal to one another near the surface and near the bottom of the domain for low-enough dissipative coefficients. We illustrate these findings by simulating numerically the QG dynamics of a horizontally homogeneous patch of ocean hosting a vertically sheared zonal current resembling the Antarctic Circumpolar Current.

Surface buoyancy gradients over a quasigeostrophic fluid permit the existence of surface-trapped Rossby waves. The interplay of these Rossby waves with surface quasigeostrophic turbulence results in latitudinally inhomogeneous mixing that, under certain conditions, culminates in a surface buoyancy staircase: a meridional buoyancy profile consisting of mixed zones punctuated by sharp buoyancy gradients, with eastward jets centred at the sharp gradients and weaker westward flows in between. In this article, we investigate the emergence of this buoyancy staircase limit in surface quasigeostrophic turbulence and we examine the dependence of the resulting dynamics on the vertical stratification. Over decreasing stratification ($\mathrm {d}N(z)/\mathrm {d}z\leq 0$, where $N(z)$ is the buoyancy frequency), we obtain flows with a longer interaction range (than in uniform stratification) and highly dispersive Rossby waves. In the staircase limit, we find straight jets that are perturbed by eastward propagating along jet waves, similar to two-dimensional barotropic $\beta$-plane turbulence. In contrast, over increasing stratification ($\mathrm {d}N(z)/\mathrm {d}z\geq 0$), we obtain flows with shorter interaction range and weakly dispersive Rossby waves. In the staircase limit, we find sinuous jets with large latitudinal meanders whose shape evolves in time due to the westward propagation of weakly dispersive along jet waves. These along jet waves have larger amplitudes over increasing stratification than over decreasing stratification, and, as a result, the ratio of domain-averaged zonal to meridional speeds is two to three times smaller over increasing stratification than over decreasing stratification. Finally, we find that, for a given Rhines wavenumber, jets over increasing stratification are closer together than jets over decreasing stratification.

Zonal jets and inertia–gravity waves are ubiquitous on planets such as Earth, Jupiter and Saturn. Motivated by the modification of energy flux of balanced flow by inertia–gravity waves, this paper studies the impact of near-inertial waves (NIWs) on zonal jets on a $\beta$-plane. Using a two-dimensional quasi-geostrophic and NIW coupled system on a $\beta$-plane (Xie & Vanneste, J. Fluid Mech., vol. 774, 2015, pp. 143–169), we find NIWs catalytically impact several features of zonal jets. The NIWs inhibit jet formation due to the waves’ catalytic induction of downscale mean energy flux. As the strength of NIWs increases, a critical point exists beyond which zonal jets are annihilated. The jet spacing is captured by the Rhines scale $L\sim \sqrt {U/\beta }$ with $U$ estimated from the upscale energy flux induced by the mean flow alone, which again shows that the NIWs’ impact is catalytic. Also, the temporal asymmetry of NIWs leads to the spatial asymmetry of jet dynamics. The jet profiles are asymmetric with a stronger shear on the left flank. And similar to the left turning of vortex dipole under the impact of NIWs, the NIW-modified jets migrate poleward. The NIWs also show a catalytic role in jet migration: the net momentum flux directly induced by NIWs is of secondary importance in the zonal mean momentum dynamics and impedes jet migration, while the advective effect of NIW-modified mean flow dominates the jet migration velocity.

The spontaneous emergence of structure is a ubiquitous process observed in fluid and plasma turbulence. These structures typically manifest as flows which remain coherent over a range of spatial and temporal scales, resisting statistically homogeneous description. This work conducts a computational and theoretical study of coherence in turbulent flows in the stochastically forced barotropic $\beta$-plane quasi-geostrophic equations. These equations serve as a prototypical two-dimensional model for turbulent flows in Jovian atmospheres, and can also be extended to study flows in magnetically confined fusion plasmas. First, analysis of direct numerical simulations demonstrates that a significant fraction of the flow energy is organized into coherent large-scale Rossby wave eigenmodes, comparable with the total energy in the zonal flows. A characterization is given for Rossby wave eigenmodes as nearly integrable perturbations to zonal flow Lagrangian trajectories, linking finite-dimensional deterministic Hamiltonian chaos in the plane to a laminar-to-turbulent flow transition. Poincaré section analysis reveals that Lagrangian flows induced by the zonal flows plus large-scale waves exhibit localized chaotic regions bounded by invariant tori, manifesting as Rossby wave breaking in the absence of critical layers. It is argued that the surviving invariant tori organize the large-scale flows into a single temporally and zonally varying laminar flow, suggesting a form of self-organization and wave stability that can account for the resilience of the observed large-amplitude Rossby waves.

The scattering of three-dimensional inertia-gravity waves by a turbulent geostrophic flow leads to the redistribution of their action through what is approximately a diffusion process in wavevector space. The corresponding diffusivity tensor was obtained by Kafiabad et al. (J. Fluid Mech., vol. 869, 2019, R7) under the assumption of a time-independent geostrophic flow. We relax this assumption to examine how the weak diffusion of wave action across constant-frequency cones that results from the slow time dependence of the geostrophic flow affects the distribution of wave energy. We find that the stationary wave-energy spectrum that arises from a single-frequency wave forcing is localised within a thin boundary layer around the constant-frequency cone, with a thickness controlled by the acceleration spectrum of the geostrophic flow. We obtain an explicit analytic formula for the wave-energy spectrum which shows good agreement with the results of a high-resolution simulation of the Boussinesq equations.

Convection occurs ubiquitously on and in rotating geophysical and astrophysical bodies. Prior spherical shell studies have shown that the convection dynamics in polar regions can differ significantly from the lower latitude, equatorial dynamics. Yet most spherical shell convective scaling laws use globally-averaged quantities that erase latitudinal differences in the physics. Here we quantify those latitudinal differences by analysing spherical shell simulations in terms of their regionalized convective heat-transfer properties. This is done by measuring local Nusselt numbers in two specific, latitudinally separate, portions of the shell, the polar and the equatorial regions, $Nu_p$ and $Nu_e$, respectively. In rotating spherical shells, convection first sets in outside the tangent cylinder such that equatorial heat transfer dominates at small and moderate supercriticalities. We show that the buoyancy forcing, parameterized by the Rayleigh number $Ra$, must exceed the critical equatorial forcing by a factor of ${\approx }20$ to trigger polar convection within the tangent cylinder. Once triggered, $Nu_p$ increases with $Ra$ much faster than does $Nu_e$. The equatorial and polar heat fluxes then tend to become comparable at sufficiently high $Ra$. Comparisons between the polar convection data and Cartesian numerical simulations reveal quantitative agreement between the two geometries in terms of heat transfer and averaged bulk temperature gradient. This agreement indicates that rotating spherical shell convection dynamics is accessible both through spherical simulations and via reduced investigatory pathways, be they theoretical, numerical or experimental.

We numerically and theoretically investigate the Boussinesq Eady model, where a rapidly rotating density-stratified layer of fluid is subject to a meridional temperature gradient in thermal wind balance with a uniform vertically sheared zonal flow. Through a suite of numerical simulations, we show that the transport properties of the resulting turbulent flow are governed by quasigeostrophic (QG) dynamics in the rapidly rotating strongly stratified regime. The ‘vortex gas’ scaling predictions put forward in the context of the two-layer QG model carry over to this fully three-dimensional system: the functional dependence of the meridional flux on the control parameters is the same, the two adjustable parameters entering the theory taking slightly different values. In line with the QG prediction, the meridional heat flux is depth-independent. The vertical heat flux is such that turbulence transports buoyancy along isopycnals, except in narrow layers near the top and bottom boundaries, the thickness of which decreases as the diffusivities go to zero. The emergent (re)stratification is set by a simple balance between the vertical heat flux and diffusion along the vertical direction. Overall, this study demonstrates how the vortex-gas scaling theory can be adapted to quantitatively predict the magnitude and vertical structure of the meridional and vertical heat fluxes, and of the emergent stratification, without additional fitting parameters.

The late-time behaviour of freely evolving quasi-geostrophic flows with initial characteristic length scale $L_0$ larger than or equal to the deformation radius $L_D$, $L_0/L_D \geq 1$, is studied. At late time the flows are dominated by large multi-level vortices consisting of ascending terraces of well-mixed potential vorticity (PV), i.e. PV staircases. We examine how the number of mixed PV levels depends on the initial conditions, in particular $L_0/L_D$. For sufficiently large values of $L_0/L_D \approx 5$, a complete staircase with regular steps forms, but as $L_0/L_D$ decreases, the staircase becomes more irregular, with fewer mixed levels and the appearance of a large step centred on zero PV, corresponding to large regions of near-zero PV separating the multi-level vortices. This occurs because weak PV features in the initial field with scales smaller than $L_D$ undergo filamentation and are coarse-grained away or homogenised. For all values of $L_0/L_D$ considered, inverse cascades of potential energy commence at sufficiently late times. The onset of these cascades, even when the flow is initialised well within the ‘asymptotic model’ (AM) regime, suggests that the AM regime is not self-consistent: when potential vorticity fronts are well-resolved, frontal dynamics eventually drive ongoing flow evolution.

Geophysical turbulent flows, characterized by rapid rotation, quantified by small Rossby number, and stable stratification, often self-organize into a collection of coherent vortices, referred to as a vortex gas. The lowest-order asymptotic expansion in Rossby number is quasigeostrophy, which has purely horizontal velocities and cyclone–anticyclone antisymmetry. Ageostrophic effects are important components of many geophysical flows and, as such, these phenomena are not well modelled by quasigeostrophy. The next-order correction in Rossby number, which includes ageostrophic effects, is the so-called balanced dynamics. Balanced dynamics includes ageostrophic vertical velocity and breaks the geostrophic cyclone–anticyclone antisymmetry. Point-vortex solutions are well known in two-dimensional and quasigeostrophic dynamics and are useful for studying the vortex-gas regime of geophysical turbulence. Here, we find point-vortex solutions in fully three-dimensional continuously stratified QG$^{+1}$ dynamics, a particular formulation of balanced dynamics. Simulations of QG$^{+1}$ point vortices show several interesting features not captured by quasigeostrophic point vortices including significant vertical transport on long time scales. The ageostrophic component of QG$^{+1}$ point vortex point-vortex dynamics renders them useful in modelling flows where quasigeostrophy filters out important physical processes.

Relative dispersion experiments are often analysed using theoretical predictions from two- and three-dimensional turbulence. These apply to infinite inertial ranges, assuming the same dispersive behaviour over all scales. With finite inertial ranges, the metrics are less conclusive. We examine this using pair separation probability density functions (PDFs), obtained by integrating a Fokker–Planck equation with different diffusivity profiles. We consider time-based metrics, such as the relative dispersion, and separation-based metrics, such as the finite scale Lyapunov exponent (FSLE). As the latter cannot be calculated from a PDF, we introduce a new measure, the cumulative inverse separation time (CIST), which can. This behaves like the FSLE, but advantageously has analytical solutions in the inertial ranges. This allows the establishment of consistency between the time- and space-based metrics, something which has been lacking previously. We focus on three dispersion regimes: non-local spreading (as in a two-dimensional enstrophy inertial range), Richardson dispersion (as in an energy inertial range) and diffusion (for uncorrelated pair motion). The time-based metrics are more successful with non-local dispersion, as the corresponding PDF applies from the initial time. Richardson dispersion is barely observed, because the self-similar PDF applies only asymptotically in time. In contrast, the separation-based CIST correctly captures the dependencies, even with a short (one decade) inertial range, and is superior to the traditional FSLE at large scales. Nevertheless, it is advantageous to use all measures together, to seek consistent indications of the dispersion.

Based on an assumption of strongly inhomogeneous potential vorticity mixing in quasi-geostrophic $\beta$-plane turbulence, a relation is obtained between the mean spacing of latitudinally meandering zonal jets and the total kinetic energy of the flow. The relation applies to cases where the Rossby deformation length is much smaller than the Rhines scale, in which kinetic energy is concentrated within the jet cores. The relation can be theoretically achieved in the case of perfect mixing between regularly spaced jets with simple meanders, and of negligible kinetic energy in flow structures other than in jets. Incomplete mixing or unevenly spaced jets will result in jets being more widely separated than the estimate, while significant kinetic energy outside the jets will result in jets closer than the estimate. An additional relation, valid under the same assumptions, is obtained between the total kinetic and potential energies. In flows with large-scale dissipation, the two relations provide a means to predict the jet spacing based only on knowledge of the energy input rate of the forcing and dissipation rate, regardless of whether the latter takes the form of frictional or thermal damping. Comparison with direct numerical integrations of the forced system shows broad support for the relations, but differences between the actual and predicted jet spacings arise both from the complex structure of jet meanders and the non-negligible kinetic energy contained in the turbulent background and in coherent vortices lying between the jets.

Inertia-gravity waves in the atmosphere and ocean are transported and refracted by geostrophic turbulent currents. Provided that the wave group velocity is much greater than the speed of geostrophic turbulent currents, kinetic theory can be used to obtain a comprehensive statistical description of the resulting interaction (Savva et al., J. Fluid Mech., vol. 916, 2021, A6). The leading-order process is scattering of wave energy along a surface of constant frequency, $\omega$, in wavenumber space. The constant-$\omega$ surface corresponding to the linear dispersion relation of inertia-gravity waves is a cone extending to arbitrarily high wavenumbers. Thus, wave scattering by geostrophic turbulence results in a cascade of wave energy to high wavenumbers on the surface of the constant-$\omega$ cone. Solution of the kinetic equations shows establishment of a wave kinetic energy spectrum $\sim k_h^{-2}$, where $k_h$ is the horizontal wavenumber.