We present a simulation-based study of the effect of a passing wave packet on underlying fully developed turbulence. We propose a novel wave-phase-resolved simulation method inspired by Helmholtz decomposition to directly couple the turbulence simulation with instantaneous wave orbital motions without wave-phase averaging. We also introduce a boundary condition treatment for the turbulence at the wave surface, which allows the turbulence simulation to be conducted in a rectangular domain while retaining the wave-phase effect. The results obtained from the proposed method reveal considerable variations in turbulence statistics, including the enstrophy and Reynolds normal stresses, during wave packet passage. Most changes occur rapidly when the narrow bandwidth around the wave packet core passes. Further analyses of the energy spectra indicate that the enhancement of turbulence occurs across a wide range of scales, with the near-surface small-scale motions experiencing the most significant intensification. Meanwhile, large-scale motions with scales comparable to the boundary layer depth are also enhanced. The mechanisms underlying the Reynolds normal stress variation at different length scales are related to the energy transfer from the wave orbital straining to turbulence through production, the pressure–strain effect, the pressure diffusion and the wave advection. By assessing the turbulence statistics and dynamics impacted by a wave packet in detail, this study provides an improved understanding of the response of a developed turbulent flow to a transient wave field. The proposed simulation method also proves to be a promising phase-resolved approach for efficiently modelling the wave effect on turbulence.

The particle trajectories in irrotational, incompressible and inviscid deep-water surface gravity waves are open, leading to a net drift in the direction of wave propagation commonly referred to as the Stokes drift, which is responsible for catalysing surface wave-induced mixing in the ocean and transporting marine debris. A balance between phase-averaged momentum density, kinetic energy density and vorticity for irrotational, monochromatic and spatially periodic two-dimensional water waves is derived by working directly within the Lagrangian reference frame, which tracks particle trajectories as a function of their labels and time. This balance should be expected as all three of these quantities are conserved following particles in this system. Vorticity in particular is always conserved along particles in two-dimensional inviscid flow, and as such even in its absence it is the value of the vorticity that fundamentally sets the drift, which in the Lagrangian frame is identified as the phase-averaged momentum density of the system. A relationship between the drift and the geometric mean water level of particles is found at the surface, which highlights connections between the geometry and dynamics. Finally, an example of an initially quiescent fluid driven by a wavelike pressure disturbance is considered, showing how the net momentum and energy from the surface pressure disturbance transfer to the wave field, and recognizing the source of the mean Lagrangian drift as the net momentum required to generate an irrotational surface wave by any conservative force.

As most mathematically justifiable Lagrangian coherent structure detection methods rely on spatial derivatives, their applicability to sparse trajectory data has been limited. For experimental fluid dynamicists and natural scientists working with Lagrangian trajectory data via passive tracers in unsteady flows (e.g. Lagrangian particle tracking or ocean buoys), obtaining material measures of fluid rotation or stretching is an active topic of research. To facilitate frame-indifferent investigations in unsteady and sparsely sampled flows, we present a novel approach to quantify fluid stretching and rotation via relative Lagrangian velocities. This technique provides a formal objective extension of quasi-objective metrics to unsteady flows by accounting for mean flow behaviour. For extremely sparse experimental data, fluid structures may be significantly undersampled and the mean flow behaviour becomes difficult to quantify. We provide a means to maintain the accuracy of our novel sparse flow diagnostics in extremely sparse sampling scenarios, such as ocean buoy data and Lagrangian particle tracking. We use data from multiple numerical and experimental flows to show that our methods can identify structures beyond existing limits of sparse, frame-indifferent diagnostics and exhibit improved interpretability over common frame-dependent diagnostics.

We consider two-dimensional flows above topography, revisiting the selective decay (or minimum enstrophy) hypothesis of Bretherton and Haidvogel. We derive a ‘condensed branch’ of solutions to the variational problem where a domain-scale condensate coexists with a flow at the (smaller) scale of the topography. The condensate arises through a supercritical bifurcation as the conserved energy of the initial condition exceeds a threshold value, a prediction that we quantitatively validate using direct numerical simulations. We then consider the forced–dissipative case, showing how weak forcing and dissipation select a single dissipative state out of the continuum of solutions to the energy-conserving system predicted by selective decay. As the forcing strength increases, the condensate arises through a supercritical bifurcation for topographic-scale forcing and through a subcritical bifurcation for domain-scale forcing, both predictions being quantitatively validated by direct numerical simulations. This method provides a way of determining the equilibrated state of forced–dissipative flows based on variational approaches to the associated energy-conserving system, such as the statistical mechanics of two-dimensional flows or selective decay.

Wave turbulence provides a conceptual framework for weakly nonlinear interactions in dispersive media. Dating from five decades ago, applications of wave turbulence theory to oceanic internal waves assigned a leading-order role to interactions characterized by a near equivalence between the group velocity of high-frequency internal waves with the phase velocity of near-inertial waves. This scale-separated interaction leads to a Fokker–Planck (generalized diffusion) equation. More recently, starting four decades ago, this scale-separated paradigm has been investigated using ray tracing methods. These ray methods characterize spectral transport of energy by counting the amplitude and net velocity of wave packets in phase space past a high-wavenumber gate prior to ‘breaking’. This explicitly advective characterization is based on an intuitive assignment and lacks theoretical underpinning. When one takes an estimate of the net spectral drift from the wave turbulence derivation and makes the corresponding assessment, one obtains a prediction of spectral transport that is an order of magnitude larger than either observations or reported ray tracing estimates. Motivated by this contradiction, we report two parallel derivations for transport equations describing the refraction of high-frequency internal waves in a sea of random inertial waves. The first uses standard wave turbulence techniques and the second is an ensemble-averaged packet transport equation characterized by the dispersion of wave packets about a mean drift in the spectral domain. The ensemble-averaged transport equation for ray tracing differs in that it contains the intuitively motivated advective term. We conclude that the aforementioned contradiction between theory, numerics and observations needs to be taken at face value and present a pathway for resolving this contradiction.

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

We study the spectral energy transfer due to wave–triad interactions in the Garrett–Munk spectrum of internal gravity waves based on a numerical evaluation of the collision integral in the wave kinetic equation. Our numerical evaluation builds on the reduction of the collision integral on the resonant manifold for a horizontally isotropic spectrum. We evaluate directly the downscale energy flux available for ocean mixing, whose value is in close agreement with the finescale parameterization. We further decompose the energy transfer into contributions from different mechanisms, including local interactions and three types of non-local interactions, namely parametric subharmonic instability, elastic scattering (ES) and induced diffusion (ID). Through analysis on the role of each mechanism, we resolve two long-standing paradoxes regarding the mechanism for forward cascade in frequency and zero ID flux for the GM76 spectrum. In addition, our analysis estimates the contribution of each mechanism to the energy transfer in each spectral direction, and reveals new understanding of the importance of local interactions and ES in the energy transfer.

The refraction of surface gravity waves by currents leads to spatial modulations in the wave field and, in particular, in the significant wave height. We examine this phenomenon in the case of waves scattered by a localised current feature, assuming (i) the smallness of the ratio between current velocity and wave group speed, and (ii) a swell-like, highly directional wave spectrum. We apply matched asymptotics to the equation governing the conservation of wave action in the four-dimensional position–wavenumber space. The resulting explicit formulas show that the modulations in wave action and significant wave height past the localised current are controlled by the vorticity of the current integrated along the primary direction of the swell. We assess the asymptotic predictions against numerical simulations using WAVEWATCH III for a Gaussian vortex. We also consider vortex dipoles to demonstrate the possibility of ‘vortex cloaking’ whereby certain currents have (asymptotically) no impact on the significant wave height. We discuss the role of the ratio of the two small parameters characterising assumptions (i) and (ii) above, and show that caustics are significant only for unrealistically large values of this ratio, corresponding to unrealistically narrow directional spectra.

Marine debris pollution is a growing problem impacting aquatic ecosystems, coastal recreation and human society. Beaches are known to be a sink for debris, and beaching needs to be accounted for in marine debris mass balances, but the process of buoyant debris beaching is not yet sufficiently well understood in order for it to be included in coastal models. We develop a simplified model for buoyant marine debris transport in bore-driven swash (where swash refers to the area that the water wets the beach with each incoming wave). We validate the model with laboratory experiments and use the combined results from the model and experiments to understand the parameters that are important for dictating particle beaching. The most relevant parameters are the particle inertia and the initial conditions with which debris particles enter the swash zone.

The transfer of momentum between the atmosphere and oceanic motions affected by the Earth's rotation occurs through the thin surface Ekman layer. The exchange depends on the surface wind stress, which produces the Ekman pumping of fluid to the ocean upper layer. The Ekman pumping mainly depends on: (i) the curl of the wind stress and (ii) the advection of vorticity due to the Ekman transport. The wind stress is usually parametrised as a quadratic function of the relative speed between the wind and the ocean currents, providing a feedback mechanism between the two fluids. Under steady and spatially uniform wind conditions over mesoscale vortices, the first mechanism generates vertical motions that induce the vortex decay (top drag), while the second promotes the horizontal advection of vorticity in the Ekman transport direction. This study examines the nonlinear effects of both mechanisms in cyclonic and anticyclonic vortices. The analyses consist of simple analytical approximations and nonlinear numerical simulations of quasi-two-dimensional vortices. When considering only the top drag mechanism, it is found that anticyclones decay faster than cyclones. By considering only the vorticity-advection effect, the vortices acquire horizontal momentum and drift; furthermore, anticyclones are reinforced while cyclones are weakened. The joint action of both mechanisms and the possible consequences on vertical transport properties are also discussed.

In this paper we study the self-induced low-Reynolds-number flow generated by a cylinder immersed in a stratified fluid. In the low Péclet limit, where the Péclet number is the ratio of the radius of the cylinder and the Phillips length scale, the flow is captured by a set of linear equations obtained by linearising the governing equations with respect to the prescribed far field conditions. We specifically focus on the low Péclet regime and develop a Green's function approach to solve the linearised equations governing the flow over the cylinder. We cross check our analytical solution against numerical solution of the nonlinear equations to obtain the range of the Péclet numbers for which the linear solution is valid. We then take advantage of the analytical solution to find explicit far-field decay rates of the flow. Our detailed analysis points out that the streamfunction and the velocity field decays algebraically in the far field. Intriguingly, this algebraic decay of the flow is much slower when compared with the exponential decay of the flow generated by a slow moving cylinder in the homogeneous Stokes regime, in the absence of stratification. Consequently, the flow generated by a cylinder in the stratified Stokes regime will have a larger domain of influence when compared with the flow generated by a cylinder in the homogeneous Stokes regime.

While large-scale seafloor features (e.g. continental slopes, mid-ocean ridges) help shape the broad patterns of ocean circulation, small-scale rough topography (e.g. seamounts, abyssal hills) can also impact the large-scale dynamics in two significant ways. First, they impact the momentum budget through flow steering, flow blocking and drag force, non-local momentum transfer via the generation of radiating internal waves and momentum dissipation by topographically induced turbulence. Second, they impact the density budget. Turbulence induced by flow–topography interactions facilitates ocean overturning circulation by upwelling dense, deep waters which form in polar oceans and sink to the abyss. Rough topography and its associated dynamics are of subgrid scale in Earth systems models (ESMs) and need parameterization for the foreseeable future. A parameterization of the intertwined impacts of flow–topography interactions on momentum and density budgets is currently non-existent, although its importance is well established. Radko (J. Fluid Mech., vol. 961, 2023, A24) provides a novel analytical model for representing the impact of rough seafloor on larger-scale flows, spanning slow to fast flow-speed regimes in homogeneous and multilayer models. The proposed ‘closure’ is in remarkable agreement with high-resolution numerical simulations and provides a crucial step forward in parameterizing the impact of rough topography in coarse-resolution ocean models. Extending the model to the full Navier–Stokes equations, linking it to the density budget and considering more realistic ocean topography are three critical future steps towards implementing Radko's theory in operational ESMs.

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.

We conduct three-dimensional direct numerical simulations to investigate the mixing, entrainment and energy budgets of gravity currents emerging from two-layer stratified locks. Depending on the density and layer thickness ratios, we find that either the upper layer or lower layer fluid can propagate faster, and that the density structure of the overall gravity current can range from strongly stratified to near-complete mixing. We furthermore observe that intermediate values of the density ratio can maximise mixing between the gravity current layers. Based on the vorticity budget, we propose a theoretical model for predicting the overall gravity current height, along with the front velocity of the two layers, for situations in which the lower layer moves faster than the upper layer. The model identifies the role of the height and thickness ratios in determining the velocity structure of the current, and it clarifies the dynamics of the ambient counter-current. A detailed analysis of the energy budget quantifies the conversion of potential into kinetic energy as a function of the governing parameters, along with the energy transfer between the different layers of the gravity current and the ambient fluid. Depending on the values of the density and layer thickness ratios, we find that the lower lock layer can gain or lose energy, whereas the upper layer always loses energy.

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.

The propagation of waves from a vertical uplift of a slender rectangular fault in a sea of constant depth is discussed, accounting for water compressibility, gravity and seabed elasticity. The compressed water column results in the generation of acoustic–gravity waves that travel at the speed of sound in water. Acoustic–gravity waves are found to terminate after a finite time, with the decay time most influenced by seabed rigidity, which is in contrast to the rigid stationary-phase model where signals persist indefinitely. At certain frequencies acoustic–gravity waves couple with the elastic seabed and travel at the shear velocity (speed of sound in an elastic solid). Improved estimates of the critical frequencies are derived. Moreover, besides the usual tsunami, a second – very small amplitude – surface wave mode travelling at the speed of sound arises under certain frequencies. We derive the cut-off frequency for this mode. The acoustic modes possess a frequency spectrum which depends on the time evolution and spatial properties of the rupture. We find that appropriate filtering of the acoustic–gravity wave signal can reveal characteristic peaks that encode information on the fault's geometry and dynamics.

Buoyant material such as microplastics accumulate near the ocean surface in regions with convergent surface currents where they can be harmful to marine life. Here, we use large eddy simulations to investigate the transport and accumulation of buoyant material in a turbulent ocean mixed layer under combined wind and convection forcing. We model non-inertial buoyant particles with a combination of buoyant tracers and Lagrangian surface particles, which allows us to explore a wide range of particle buoyancies. Surface cooling drives convection, and under this regime persistent convective vortices form that trap buoyant particles, leading to large concentrations. Despite their small size, the convective vortices exhibit a bias towards cyclonic vorticity that has not been reported previously. Based on an analysis of Lagrangian trajectories, the average time that a particle spends inside a convective vortex is long enough for planetary vorticity to become important and further vortex stretching causes an exponential increase in vorticity. When wind forcing is included, there is a transition from convective cells to longitudinal wind rolls with three distinct flow patterns observed under weak, moderate and strong wind forcing. For sufficiently weak winds, convective vortices survive but are less effective at trapping buoyant material. Under strong wind forcing, convective vortices no longer exist, but some clustering occurs in regions of high speed associated with longitudinal wind rolls. We quantify the degree of clustering using the Gini coefficient and find that clustering is strongly influenced by the relative size of the friction and convective velocities and the particle buoyancy.

The role of the Lagrangian mean flow, or drift, in modulating the geometry, kinematics and dynamics of rotational and irrotational deep-water surface gravity waves is examined. A general theory for permanent progressive waves on an arbitrary vertically sheared steady Lagrangian mean flow is derived in the Lagrangian reference frame and mapped to the Eulerian frame. A Lagrangian viewpoint offers tremendous flexibility due to the particle labelling freedom and allows us to reveal how key physical wave behaviour arises from a kinematic constraint on the vorticity of the fluid, inter alia the nonlinear correction to the phase speed of irrotational finite amplitude waves, the free surface geometry and velocity in the Eulerian frame, and the connection between the Lagrangian drift and the Benjamin–Feir instability. To complement and illustrate our theory, a small laboratory experiment demonstrates how a specially tailored sheared mean flow can almost completely attenuate the Benjamin–Feir instability, in qualitative agreement with the theory. The application of these results to problems in remote sensing and ocean wave modelling is discussed. We provide an answer to a long-standing question: remote sensing techniques based on observing current-induced shifts in the wave dispersion will measure the Lagrangian, not the Eulerian, mean current.

Bifurcation analysis of confined salt-finger convection using single-mode equations obtained from a severely truncated Fourier expansion in the horizontal is performed. Strongly nonlinear staircase-like solutions having, respectively, one (S1), two (S2) and three (S3) regions of mixed salinity in the vertical direction are computed using numerical continuation, and their stability properties are determined. Near onset, the one-layer S1 solution is stable and corresponds to maximum salinity transport among the three solutions. The S2 and S3 solutions are unstable but exert an influence on the statistics observed in direct numerical simulations (DNS) in larger two-dimensional (2-D) domains. Secondary bifurcations of S1 lead either to tilted-finger (TF1) or to travelling wave (TW1) solutions, both accompanied by the spontaneous generation of large-scale shear, a process favoured for lower density ratios and Prandtl numbers ($Pr$). These states at low $Pr$ are associated, respectively, with two-layer and three-layer staircase-like salinity profiles in the mean. States breaking reflection symmetry in the midplane are also computed. In two dimensions and for low $Pr$, the DNS results favour direction-reversing tilted fingers resembling the pulsating wave state observed in other systems. Two-layer and three-layer mean salinity profiles corresponding to reversing tilted fingers and TW1 are observed in 2-D DNS averaged over time. The single-mode solutions close to the high wavenumber onset are in an excellent agreement with 2-D DNS in small horizontal domains and compare well with 3-D DNS.

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.