Waves transport particles in the direction of wave propagation with the Stokes drift. When the Earth’s rotation is accounted for, waves induce an additional (Eulerian-mean) current that reduces drift and is known as the anti-Stokes drift. This effect is often ignored in oceanic particle-tracking simulations, despite being important. Although different theoretical models exist, they have not been validated by experiments. We conduct laboratory experiments studying the surface drift induced by deep-water waves in a purpose-built rotating wave flume. With rotation, the Lagrangian-mean drift deflects to the right (counterclockwise rotation) and reduces in magnitude. Compared with two existing steady theoretical models, measured drift speed follows a similar trend with wave Ekman number but is larger. The difference is largely explained by unsteadiness on inertial time scales. Our results emphasise the importance of considering unsteadiness when predicting and analysing the transport of floating material by waves.

Ocean submesoscales, flows with characteristic size $10\,\text{m}{-}10\,\text{km}$, are transitional between the larger, rotationally constrained mesoscale and three-dimensional turbulence. In this paper, we present simulations of a submesoscale ocean filament. In our case, the filament is strongly sheared in both vertical and cross-filament directions, and is unstable. Instability indeed dominates the early behaviour with a fast extraction of kinetic energy from the vertically sheared thermal wind. However, the instability that emerges does not exhibit characteristics that match the perhaps expected symmetric or Kelvin–Helmholtz instabilities, and appears to be non-normal in nature. The prominence of the transient response depends on the initial noise, and for large initial noise amplitudes, saturates before symmetric instability normal modes are able to develop. The action of the instability is sufficiently rapid – with energy extraction from the mean flow emerging and peaking within the first inertial period ($\sim\! 18\ \text{h}$) – that the filament does not respond in a geostrophically balanced sense. Instead, at all initial noise levels, it later exhibits vertically sheared near-inertial oscillations with higher amplitude as the initial minimum Richardson number decreases. Horizontal gradients strengthen only briefly as the fronts restratify. These unstable filaments can be generated by strong mixing events at pre-existing stable structures; we also caution against inadvertently triggering this response in idealised studies that start in a very unstable state.

We investigate the motion of weakly negatively buoyant spheres settling in surface gravity waves using laboratory experiments. The trajectories of the settling spheres are tracked over most of the water depth with simultaneous measurements of the background fluid flow. These experiments are conducted in the regime relevant for environmental and geophysical applications where both particle inertia and fluid inertia are important. Using these data, we show that the sphere motion is well described by the kinematic sum of the undisturbed fluid velocity and the particle terminal settling velocity as long as the fluid inertia is not too large. We show how this result can be understood in the context of an ad hoc Maxey–Riley–Gatignol-type equation where the drag on the particle is given by the Schiller–Naumann drag correlation. We also evaluate whether inertial particles experience enhanced settling in waves, finding that measurement uncertainties in the particle terminal settling velocity and the presence of Eulerian-mean flows do not allow the small percentage increase in the settling velocity to be measured. When the fluid inertia becomes large enough, we observe path instabilities caused by particle wake effects in both quiescent and wavy conditions. However, the particle velocity fluctuations associated with the path instabilities are unaffected by the background flow. The minimal influence of the wavy flow on the particle path instabilities is thought to be due to the large-scale separation between the waves and the particle.

Rogue waves are associated with various ocean processes, both at the coast and in the open ocean. In either zone, inhomogeneities in the wave field caused by shoaling, crossing seas or current interactions disturb the wave statistics, increasing the rogue wave probability and magnitude. Such amplification of the frequency of rogue waves and their intensity, i.e. the maximum normalised height, have been attested to in numerical simulations and laboratory studies, in particular for wave–current interactions. In this study, we investigate the effect of the current intensity and direction on rogue wave probability, by analysing long-term observations from the southern North Sea. We observe that the amplification is similar for opposing and following currents. Despite the sea states being dominantly broadbanded and featuring a large directional spread, the anomalous statistics are of the same order of magnitude as those observed in unidirectional laboratory experiments for stationary currents.

Large-eddy simulations are analysed to determine the influence of suspended canopies, such as those formed in macroalgal farms, on ocean mixed layer (OML) deepening and internal wave generation. In the absence of a canopy, we show that Langmuir turbulence, when compared with wind-driven shear turbulence, results in a deeper OML and more pronounced internal waves beneath the OML. Subsequently, we examine simulations with suspended canopies of varying densities located in the OML, in the presence of a background geostrophic current. Intensified turbulence occurs in the shear layer at the canopy’s bottom edge, arising from the interaction between the geostrophic current and canopy drag. Structures resembling Kelvin–Helmholtz (KH) instability emerge as the canopy shear layer interacts with the underlying stratification, radiating internal waves beneath the OML. Both intensified turbulence and lower-frequency motions associated with KH-type structures are critical factors in enhancing mixing. Consequently, the OML depth increases by up to a factor of two compared with cases without a canopy. Denser canopies and stronger geostrophic currents lead to more pronounced KH-type structures and internal waves, stronger turbulence and greater OML deepening. Additionally, vertical nutrient transport is enhanced as the OML deepens due to the presence of the canopy. Considering that the canopy density investigated in this study closely represents offshore macroalgal farms, these findings suggest a mechanism for passive nutrient entrainment conducive to sustainable farming. Overall, this study reveals the intricate interactions between the suspended canopy, turbulent mixing and stratification, underscoring their importance in reshaping OML characteristics.

Hydrodynamic consequences of using simpler geometric shapes to represent coral canopies are examined through a laboratory study. A canopy composed of cylinders is compared with a canopy composed of 3-D-printed, scaled down coral heads in a recirculating flume. Vertical velocity profiles are measured at four horizontal locations for each canopy type, and mean velocity and turbulence statistics are compared both within and above the canopy. A narrow, defined wake on the scale of the canopy element is present behind the cylinder canopy elements that is absent in the coral canopy. There is also a peak in shear stress at the top of the cylinder canopy, likely due to the sharp edge at the top of the cylinder. Above the canopy, however, turbulence statistics and friction velocities behave similarly for both canopy types. Therefore, the results indicate we may reasonably get coral reef drag estimates from canopies with simpler geometric surrogates, especially when the mean free-stream and within-canopy flow speeds are matched to environmental conditions.

We have investigated the dynamics of floating tracer in an idealised turbulent quasi-geostrophic ocean by advecting Lagrangian particles in a high-resolution velocity field enhanced by the potential flow associated with vortex stretching. At first order in the Rossby number expansion, this component of the ageostrophic circulation can be derived through a diagnostic equation in terms of the geostrophic velocities. Borrowing methods from the theory of Lagrangian coherent structures, we identify coherent material loops around strong vortex cores using the Lagrangian averaged vorticity deviation (LAVD). Building on studies of clustering in kinematic, stochastic velocity fields, we utilise methods from statistical topography to show that the coherent vortices dominate the distribution of extreme values of the concentration field. We find that the presence of clusters and voids in a coherent vortex depends on more than just the sense of rotation, but also on the full evolution of the vorticity over its lifecycle. We identify the mechanism behind the cluster formation that respects the symmetries of the quasi-geostrophic equations but can be expected to hold robustly in more complicated regimes, due to the simple physical description. The association of cluster formation with vortex stretching implies that LAVD is a particularly relevant metric for floating tracer dynamics. The detection of intense clustering also has implications for reaction rates between ocean-borne flotsam, meaning that our results are relevant to understanding the general risk of floating microplastics and marine biological populations.

The transformation of internal waves on a stepwise underwater obstacle is studied in the linear approximation. The transmission and reflection coefficients are derived for a two-layer fluid. The results are obtained and presented as functions of incident wave wavenumber, density ratio of layers, pycnocline position, and height of the bottom step. Excitation coefficients of evanescent modes are also calculated, and their importance is demonstrated. This allows one to estimate the number of evanescent modes necessary to take into account to attain the required accuracy for the transformation coefficients.

When atmospheric storms pass over the ocean, they resonantly force near-inertial waves (NIWs), internal waves with a frequency close to the local Coriolis frequency $f$. It has long been recognised that the evolution of NIWs is modulated by the ocean's mesoscale eddy field. This can result in NIWs being concentrated into anticyclones which provide an efficient pathway for NIW propagation to depth. Here we analyse the eigenmodes of NIWs in the presence of mesoscale eddies and heavily draw on parallels with quantum mechanics. Whether the eddies are effective at modulating the behaviour of NIWs depends on the wave dispersiveness $\varepsilon ^2 = f\lambda ^2/\varPsi$, where $\lambda$ is the deformation radius and $\varPsi$ is a scaling for the eddy streamfunction. If $\varepsilon \gg 1$, NIWs are strongly dispersive, and the waves are only weakly affected by the eddies. We calculate the perturbations away from a uniform wave field and the frequency shift away from $f$. If $\varepsilon \ll 1$, NIWs are weakly dispersive, and the wave evolution is strongly modulated by the eddy field. In this weakly dispersive limit, the Wentzel–Kramers–Brillouin approximation, from which ray tracing emerges, is a valid description of the NIW evolution even if the large-scale atmospheric forcing apparently violates the requisite assumption of a scale separation between the waves and the eddies. The large-scale forcing excites many wave modes, each of which varies on a short spatial scale and is amenable to asymptotic analysis analogous to the semi-classical analysis of quantum systems. The strong modulation of weakly dispersive NIWs by eddies has the potential to modulate the energy input into NIWs from the wind, but we find that this effect should be small under oceanic conditions.

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.

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.