A series of recent studies has indicated that the component of the bottom drag caused by irregular small-scale topography in the ocean varies non-monotonically with the flow speed. The roughness-induced forcing increases with the speed of relatively slow abyssal currents but, somewhat counterintuitively, starts to decrease when flows are sufficiently swift. This reduction in drag at high speeds leads to the instability of laterally uniform currents, and the resulting evolutionary patterns are explored using numerical and analytical methods. The drag-law instability manifests in the spontaneous emergence of parallel jets, aligned in the direction of the basic flow and separated by relatively quiescent regions. We hypothesize that the mechanisms identified in this investigation could play a role in the dynamics of zonal striations commonly observed in the ocean.

Seafloor roughness profoundly influences the pattern and dynamics of large-scale oceanic flows. However, these kilometre-scale topographic patterns are unresolved by global numerical Earth system models and will remain subgrid for the foreseeable future. To properly represent the effects of small-scale bathymetry in analytical and coarse-resolution numerical models, we develop the stratified ‘sandpaper’ theory of flow–topography interaction. This model, which is based on the multilayer shallow-water framework, extends its barotropic antecedent to stratified flows. The proposed theory is successfully tested on the configuration representing the interaction of a zonal current with a corrugated cross-flow ridge.

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.

Recent studies reveal the dramatic impact of seafloor roughness on the dynamics and stability of broad oceanic flows. These findings motivate the development of parameterizations that concisely represent the effects of small-scale bathymetric patterns in theoretical and coarse-resolution numerical circulation models. The previously reported quasi-geostrophic ‘sandpaper’ theory of flow–topography interaction a priori assumes gentle topographic slopes and weak flows with low Rossby numbers. Since such conditions are often violated in the ocean, we now proceed to formulate a more general model based on shallow-water equations. The new version of the sandpaper model is validated by comparing roughness-resolving and parametric simulations of the flow over a corrugated seamount.

Surface differential heating on a stably stratified fluid body drives an overturning circulation confined to the upper fluid region – here coined stratified horizontal convection (SHC). In this manuscript, we investigate the dynamics of SHC via laboratory experiments, exploring local and global flow properties. By considering the available potential energy of the system, we derive a unique length scale of SHC and introduce the Péclet number $Pe$ that captures both the stabilising effect of stratification and the destabilising effect of the baroclinic adjustment. We found that $Pe$ characterises local and global flow properties, including the fluid transport of the overturning circulation, the available mechanical energy and the flow dimensionality. Our study provides insights into the fluid dynamics of stratified environments that experience horizontal convection, such as lakes, oceans and atmospheres.

An analytical model is developed which explores the impact of irregular sea-floor roughness on large-scale oceanic flows. The previously reported asymptotic ‘sandpaper’ theory of flow-topography interaction represents relatively swift currents and exhibits singular behaviour in the weak flow limit. The present investigation systematically spans a wider parameter space and identifies the principal dissimilarities in the topographic regulation of slow and fast currents. The fast flows are controlled by the Reynolds stresses produced by topographically generated eddies. In contrast, relatively weak flows are more affected by the eddy-induced bottom form drag. The asymptotic models for fast and slow currents are then combined to arrive at a concise description of flow forcing by small-scale topography in homogeneous and multilayer models. The proposed closure is validated by comparing corresponding topography-resolving and parametric simulations.

A multiscale asymptotic theory is formulated for surface gravity waves and currents in finite-depth water with a vegetation canopy that provides a drag force on both flows with known drag coefficients. It assumes that the density is uniform and the depth is uniform pro tem and that the wave frequency is fast compared to the current advective rate. It is a quasi-linear theory in which the wave dynamics is independent of current and drag to leading order but provides perturbative corrections, and in which wave nonlinear interactions are neglected while quadratic wave-averaged wave fluxes and quadratic wave-drag effects are retained. The primary surface wave is modified by drag and current interactions, and the wave-averaged current momentum balance includes a wave-augmented drag force and several vortex forces due to Earth's rotation, current vorticity, Stokes drift and drag-induced wave vorticity. The wave-averaged current equations derived here are a suitable basis for future large-eddy simulation and submesoscale circulation computational models.

Many applications in geosciences require solving inverse problems to estimate the state of a physical system. Data assimilation provides a strong framework to do so when the system is partially observed and its underlying dynamics are known to some extent. In the variational flavor, it can be seen as an optimal control problem where initial conditions are the control parameters. Such problems are often ill-posed, regularization may be needed using explicit prior knowledge to enforce a satisfying solution. In this work, we propose to use a deep prior, a neural architecture that generates potential solutions and acts as implicit regularization. The architecture is trained in a fully-unsupervised manner using the variational data assimilation cost so that gradients are backpropagated through the dynamical model and then through the neural network. To demonstrate its use, we set a twin experiment using a shallow-water toy model, where we test various variational assimilation algorithms on an ocean-like circulation estimation.

This study explores the impact of small-scale variability in the bottom relief on the dynamics and evolution of broad baroclinic flows in the ocean. The analytical model presented here generalizes the previously reported barotropic ‘sandpaper’ theory of flow–topography interaction to density-stratified systems. The multiscale asymptotic analysis leads to an explicit representation of the large-scale effects of irregular bottom roughness. The utility of the multiscale model is demonstrated by applying it to the problem of topography-induced spin-down of an axisymmetric vortex. We find that bathymetry affects vortices by suppressing circulation in their deep regions. As a result, vortices located above rough topography tend to be more stable than their flat-bottom counterparts. The multiscale theory is validated by comparing corresponding topography-resolving and parametric simulations.

Maas (J. Fluid Mech., vol. 684, 2011, pp. 5–24) showed that, for an oscillating two-dimensional barotropic tide flowing over sub-critical topography of compact support, some topographic forms existed that produced non-radiating baroclinic disturbances. The problem is related to ‘stealth’ and ‘cloaking’ problems. Here Maas's result is derived using a simpler approach, not involving complicated mappings, but formally restricted to perturbation topography. Wider results come from the discussion of nearly compact support topographic disturbances provided by Schwartz functions with weak high-wavenumber radiation and by exploiting both a known functional equation formulation and Fourier methods. The problem is extended to disturbances on uniform slopes. A variety of non-radiating topographies can be found, although they are mathematically delicate and unlikely to be found in nature. Topography with weak radiation at high wavenumber is a much wider class of structures. Application of these solutions would lie with the ability to estimate dissipation over and near the topography from motions observed at a distance.

This study examines the impact of small-scale irregular topographic features on the dynamics and evolution of large-scale barotropic flows in the ocean. A multiscale theory is developed, which makes it possible to represent large-scale effects of the bottom roughness without explicitly resolving small-scale variability. The analytical model reveals that the key mechanism of topographic control involves the generation of a small-scale eddy field associated with considerable Reynolds stresses. These eddy stresses are inversely proportional to the large-scale velocity and adversely affect mean circulation patterns. The multiscale model is applied to the problem of topography-induced spin-down of a large circularly symmetric vortex and is validated by corresponding topography-resolving simulations. The small-scale bathymetry chosen for this configuration conforms to the Goff–Jordan statistical spectrum. While the multiscale model formally assumes a substantial separation between the scales of interacting flow components, it is remarkably accurate even when scale separation is virtually non-existent.

Lagrangian averaging has been shown to be more effective than the Eulerian mean in separating waves from slow dynamics in two time scale flows. It also appears in many reduced models that capture the wave feedback on the slow flow. Its calculation, however, requires tracking particles in time, which imposes several difficulties in grid-based numerical simulations or estimation from fixed-point measurements. To circumvent these difficulties, we propose a grid-based iterative method to calculate the Lagrangian mean without tracking particles in time, which also reduces computation, memory footprint and communication between processors in parallelised numerical models. To assess the accuracy of this method several examples are examined and discussed. We also explore an application of this method in the context of shallow-water equations by quantifying the validity of wave-averaged geostrophic balance – a modified form of geostrophic balance accounting for the effect of strong waves on slow dynamics.

When a fluid stream in a conduit splits in order to pass around an obstruction, it is possible that one branch will be critically controlled while the other remains not so. This is apparently the situation in Pacific Ocean abyssal circulation, where most of the northward flow of Antarctic bottom water passes through the Samoan Passage, where it is hydraulically controlled, while the remainder is diverted around the Manihiki Plateau and is not controlled. These observations raise a number of questions concerning the dynamics necessary to support such a regime in the steady state, the nature of upstream influence and the usefulness of rotating hydraulic theory to predict the partitioning of volume transport between the two paths, which assumes the controlled branch is inviscid. Through the use of a theory for constant potential vorticity flow and accompanying numerical model, we show that a steady-state regime similar to what is observed is dynamically possible provided that sufficient bottom friction is present in the uncontrolled branch. In this case, the upstream influence that typically exists for rotating channel flow is transformed into influence into how the flow is partitioned. As a result, the partitioning of volume flux can still be reasonably well predicted with an inviscid theory that exploits the lack of upstream influence.

Parameterizing mesoscale eddies in ocean circulation models remains an open problem due to the ambiguity with separating the eddies from large-scale flow, so that their interplay is consistent with the resolving skill of the employed non-eddy-resolving model. One way to address the issue is by using recently formulated dynamically filtered eddies. These eddies are obtained as the field errors of fitting some given reference ocean circulation into the employed coarse-grid ocean model. The main strengths are (i) no explicit spatio-temporal filter is needed for separating the large-scale and eddy flow components, (ii) the eddies are dynamically translated into the error-correcting forcing that perfectly augments the coarse-grid model towards reproducing the reference circulation. We uncovered physical properties of the eddies by interpreting involved nonlinear eddy/large-scale interactions via the classical flux-gradient relation. We described the eddies in terms of their full, space–time dependent transport tensor, which was made unique by constraining it to be the same for the potential vorticity, momentum and buoyancy fluxes. Both diffusive and advective parts of the transport tensor were found to be significant. The diffusive tensor component is characterised by polar eigenvalues and is further decomposed into isotropic and filamentation components. The latter component completely dominates, therefore, it should be taken into account by eddy parameterizations, which is not yet the case. We also showed that spatial inhomogeneities of the transport tensor components are important. Comparing these properties with those obtained for more common, locally filtered eddies revealed that they are distinctly different.

The meridional component of the earth's rotation is often neglected in geophysical contexts. This is referred to as the ‘traditional approximation’ and is justified by the typically small vertical velocity and aspect ratio of such problems. Ocean fronts are regions of strong horizontal buoyancy gradient and are associated with strong vertical transport of tracers and nutrients. Given these comparatively large vertical velocities, non-traditional rotation may play a role in governing frontal dynamics. Here the effects of non-traditional rotation on a front in turbulent thermal wind balance are considered using an asymptotic approach. Solutions are presented for a general horizontal buoyancy profile and examined in the simple case of a straight front. Non-traditional effects are found to depend strongly on the direction of the front and may lead to the generation of jets and the modification of the frontal circulation and vertical transport.

We present a comprehensive study of the mesoscale eddy forcing in the ocean by proposing spatially local filtering of the high-resolution double-gyre ocean circulation solution into its large- and small-scale (eddy) components. The large-scale component is dominated by the mid-latitude gyres, the western boundary currents and their highly transient eastward jet extension; the eddy component is concentrated around the eastward jet and strongly interacts with it. The proposed decomposition method achieves flow filtering based on the local spatial correlations. This is different from the existing decomposition methods, e.g. classical Reynolds decomposition and moving-average (spatial) filtering with a constant filter size based on the first baroclinic Rossby deformation radius. Next, we characterize the dynamical impacts of the resulting eddy forcing on the large-scale flow in terms of their mutual time-lagged spatial correlations, formulated as product integral characteristics. Its temporal statistics uncover robust causality between the eddy forcing and the induced large-scale potential vorticity anomalies – referred to as the eddy backscatter. The results also prove the significance of the transient eddy forcing and the time lag dependence of the eddy backscatter. We argue that these properties are to be considered by eddy parametrization schemes. We further used the decomposed eddy fields to augment a coarse-resolution ocean model. The augmented solution statistically reproduces the missing eastward jet extension, enhances the eddy activities around it and recovers the essential large-scale low-frequency variability. This justifies a reduced-order statistical emulation of the eddies – an emerging methodology for including eddy effects in non-eddy-resolving ocean models.

The problem of defining oceanic mesoscale eddies remains generally unresolved because there is no unique local spatio-temporal filter that can be used for extracting the eddies, and it is unclear what part of the eddy field cannot be actually resolved and needs to be parameterized in a coarse-grid model. We propose using of the coarse-grid model itself for reconstructing dynamically unresolved eddies, which are actually field errors on the top of the dynamically resolved, large-scale reference flow. The novelties and strengths of the approach are that (i) no spatio-temporal filtering is ever needed, (ii) field errors are dynamically translated into the error-correcting forcing and (iii) the latter exactly augments the coarse-grid model solution towards the reference flow. After implementation of the proposed approach, we study statistical properties of the field errors, show their robustness and reveal their significant differences from the locally filtered eddies. We argue that dynamical effects of unresolved eddies can be ultimately parameterized by emulating field errors and closing them on the dynamically resolved flow. So far, our results are limited to the quasigeostrophic approximation, but this serves as a proof of concept and starting point for the follow-up extension into the primitive equations, which are used routinely in the comprehensive oceanic general circulation models.

This study examines the influence of flow nonlinearity in western boundary layers upon the turbulent wind-driven ocean gyres. Our analysis involves comparisons between large-scale circulation properties of the linear and nonlinear states, as well as a Lagrangian particle analysis of relevant flow features. We find that the so-called counter-rotating gyre anomalies, which are nonlinear circulation features embedded in the gyres, are consistent in shape with the linear, weakened, wind-curl response created by the geometric wind effect. However, the linear response is far too weak without considering nonlinear effects. Within the western boundary layer lobe of these features, the nonlinear boundary layer has a pivotal impact upon the global circulation. Effects of potential vorticity advection inhibit viscous relative vorticity fluxes through the western boundary. This creates a significant potential vorticity imbalance between the gyres. Consequently, this generates an accumulation of enstrophy downstream in the inertial recirculation zones, which in turn supports the eastward jet. However, within the ocean basin, the growing imbalance is eventually rectified by inter-gyre potential vorticity exchanges owing to nonlinear fluxes. The Lagrangian particle analysis reveals the inter-gyre exchange mechanism, where particles seeded within the western boundary layer migrate between the gyres and weaken the eastward jet extension.

Conservation laws that relate the local time-rate-of-change of the spatial integral of a density function to the divergence of its flux through the boundaries of the integration domain provide integral constraints on the spatio-temporal development of a field. Here we show that a new type of conserved quantity exists that does not require integration over a particular domain but which holds locally at any point in the field. This is derived for the pseudo-energy density of non-divergent Rossby waves where local invariance is obtained for (i) a single plane wave, and (ii) waves produced by an impulsive point source of vorticity. The definition of pseudo-energy used here consists of a conventional kinetic part, as well as an unconventional pseudo-potential part, proposed by Buchwald (Proc. R. Soc. Lond. A, vol. 328, issue 1572, 1972, pp. 37–48). The anisotropic nature of the energy flux that appears in response to the point source further clarifies the role of the beta plane in the observed western intensification of ocean currents.

In the problem of horizontal convection a non-uniform buoyancy, $b_{s}(x,y)$, is imposed on the top surface of a container and all other surfaces are insulating. Horizontal convection produces a net horizontal flux of buoyancy, $\boldsymbol{J}$, defined by vertically and temporally averaging the interior horizontal flux of buoyancy. We show that $\overline{\boldsymbol{J}\boldsymbol{\cdot }\unicode[STIX]{x1D735}b_{s}}=-\unicode[STIX]{x1D705}\langle |\unicode[STIX]{x1D735}b|^{2}\rangle$; the overbar denotes a space–time average over the top surface, angle brackets denote a volume–time average and $\unicode[STIX]{x1D705}$ is the molecular diffusivity of buoyancy $b$. This connection between $\boldsymbol{J}$ and $\unicode[STIX]{x1D705}\langle |\unicode[STIX]{x1D735}b|^{2}\rangle$ justifies the definition of the horizontal-convective Nusselt number, $Nu$, as the ratio of $\unicode[STIX]{x1D705}\langle |\unicode[STIX]{x1D735}b|^{2}\rangle$ to the corresponding quantity produced by molecular diffusion alone. We discuss the advantages of this definition of $Nu$ over other definitions of horizontal-convective Nusselt number. We investigate transient effects and show that $\unicode[STIX]{x1D705}\langle |\unicode[STIX]{x1D735}b|^{2}\rangle$ equilibrates more rapidly than other global averages, such as the averaged kinetic energy and bottom buoyancy. We show that $\unicode[STIX]{x1D705}\langle |\unicode[STIX]{x1D735}b|^{2}\rangle$ is the volume-averaged rate of Boussinesq entropy production within the enclosure. In statistical steady state, the interior entropy production is balanced by a flux through the top surface. This leads to an equivalent ‘surface Nusselt number’, defined as the surface average of vertical buoyancy flux through the top surface times the imposed surface buoyancy $b_{s}(x,y)$. In experimental situations it is easier to evaluate the surface entropy flux, rather than the volume integral of $|\unicode[STIX]{x1D735}b|^{2}$ demanded by $\unicode[STIX]{x1D705}\langle |\unicode[STIX]{x1D735}b|^{2}\rangle$.