Bottom Boundary Layer Research Articles

Tide-Induced Mixing in the Bottom Boundary Layer in the Western East China Sea

Tide-Induced Mixing in the Bottom Boundary Layer in the Western East China Sea

Analytical Eddy Viscosity Model for Turbulent Wave Boundary Layers: Application to Suspended Sediment Concentrations over Wave Ripples

Turbulence related to flow oscillations near the seabed, in the wave bottom boundary layer (WBBL), is the phenomenon responsible for the suspension and transport of sediments. The vertical distribution of turbulent eddy viscosity within the WBBL is a key parameter that determines the vertical distribution of suspended sediments. For practical coastal engineering applications, the most used method to parameterize turbulence consists in specifying the shape of the one-dimensional-vertical (1DV) profile of eddy viscosity. Different empirical models have been proposed for the vertical variation of eddy viscosity in the WBBL. In this study, we consider the exponential-type profile, which was validated and calibrated by direct numerical simulation (DNS) and experimental data for turbulent channel and open-channel flows, respectively. This model is generalized to the WBBL, and the period-averaged eddy viscosity is calibrated by a two-equation baseline (BSL) k-ω model for different conditions. This model, together with a β-function (where β is the inverse of the turbulent Schmidt number), is used in modeling suspended sediment concentration (SSC) profiles over wave ripples, where field and laboratory measurements of SSC show two kinds of concentration profiles depending on grain particles size. Our study shows that the convection–diffusion equation, for SSC in WBBLs over sand ripples with an upward convection term, reverts to the classical advection–diffusion equation (ADE) with an “apparent” sediment diffusivity εs*=α εs related to the sediment diffusivity εs by an additional parameter α associated with the convective sediment entrainment process over sand ripples, which is defined by two equations. In the first, α depends on the relative importance of upward convection related to coherent vortex shedding and downward settling of sediments. When the convective transfer is very small, above low-steepness ripples, α≈1. In the second, α depends on the relative importance of coherent vortex shedding and random turbulence. When random turbulence is more important than coherent vortex shedding, α≈1, and “apparent” sediment diffusivity reverts to the classical sediment diffusivity εs*≈ εs. Comparisons with experimental data show that the proposed method allows a good description of both SSC for fine and coarse sand and “apparent” sediment diffusivity εs*.

Coupling between Abyssal Boundary Layers and the Interior Ocean in the Absence of Along-Slope Variations

Abstract To close the overturning circulation, dense bottom water must upwell via turbulent mixing. Recent studies have identified thin bottom boundary layers (BLs) as locations of intense upwelling, yet it remains unclear how they interact with and shape the large-scale circulation of the abyssal ocean. The current understanding of this BL–interior coupling is shaped by 1D theory, suggesting that variations in locally produced BL transport generate exchange with the interior and thus a global circulation. Until now, however, this picture has been based on a 1D theory that fails to capture the local evolution in even highly idealized 2D geometries. The present work applies BL theory to revised 1D dynamics, which more naturally generalizes to two and three dimensions. The BL is assumed to be in quasi-equilibrium between the upwelling of dense water and the convergence of downward buoyancy fluxes. The BL transport, for which explicit formulas are presented, exerts an influence on the interior by modifying the bottom boundary condition. In 1D, this BL transport is independent of the interior evolution, but in 2D the BL and interior are fully coupled. Once interior variables and the bottom slope are allowed to vary in the horizontal, the resulting convergences and divergences in the BL transport exchange mass with the interior. This framework allows for the analysis of previously inaccessible problems such as the BL–interior coupling in the presence of an exponential interior stratification, laying the foundation for developing a full theory for the abyssal circulation.

Enhanced Biogeochemical Cycling Along the U.S. West Coast Shelf

AbstractContinental margins play an essential role in global ocean biogeochemistry and the carbon cycle; however, global assessments of this role remain highly uncertain. This uncertainty arises from the large variability over a broad range of temporal and spatial scales of the processes that characterize these environments. High‐resolution simulations with ocean biogeochemical models have emerged as essential tools to advance biogeochemical assessments at regional scales. Here, we examine the processes and balances for carbon, oxygen, and nitrogen cycles along the U.S. West Coast in an 11‐year hindcast simulation with a submesoscale‐permitting oceanic circulation coupled to a biogeochemical model. We describe and quantify the biogeochemical cycles on the continental shelf, and their connection to the broader regional context encompassing the California Current System. On the shelf, coastal and wind stress curl upwelling drive a vigorous overturning circulation that supports biogeochemical rates and fluxes that are approximately twice as large as offshore. Exchanges with the proximate sediments, submesoscale shelf currents, bottom boundary layer transport, and intensified cross‐shelf export of shelf‐produced materials further impact coastal and open‐ocean balances. While regional variability prevents extrapolation of our results to global margins, our approach provides a powerful tool to identify the dominant dynamics in different shelf setting and quantify their large‐scale consequences.

Eddy viscosity profiles and currents characterization of the bottom Ekman layer for the Western Gulf of Mexico continental shelf

By analyzing a 4-year time series of currents and temperature from a set of moorings deployed over the Western Gulf of Mexico shelf, typical bottom Ekman layer parameters as the buoyancy frequency and bottom stress were estimated. Mean vertical structure and variance of the near-bottom currents are presented and the near-neutral stratification assumption is discussed considering what is expected from a turbulent bottom Ekman layer. Constant values for eddy viscosity (K) are obtained here from measured bottom Ekman spirals, which are further used to reproduce theoretical K profiles and the expected variability at each location. The results offer the first regional characterization of the bottom Ekman boundary layer, and are intended to serve in the calculation of mesoscale cross-shelf transport, as well as a reference to the turbulent parameters necessary in the development of benthic biogeochemical models.

Internal solitary waves enhancing turbulent mixing in the bottom boundary layer of continental slope

The importance of internal solitary waves (ISWs) to the energy and material transport in the upper ocean has been confirmed, but how ISWs affect turbulent mixing in the bottom boundary layer is lack of direct field observation. To reveal the characteristics of ISWs and the turbulent mixing they cause in the bottom boundary layer, a 36-h fine-scale observation of ISWs and their influence on turbulent mixing was conducted at the South China Sea slope (water depth: 659 m). During the observation, a group of ISWs passed by with horizontal and vertical velocities of 0.2 and 0.04 m/s at 0.1 m above bottom (mab) respectively, resulting in strong velocity shear exceeding 0.5 m/s per meter. When the ISWs passed through the slope, although they did not deform and break near the seabed, they increased the bottom shear stress, turbulent kinetic energy production rate and turbulent kinetic energy dissipation rate ( ε ) at 0.6 mab several times. The ISWs induced ε reached O (10 −6 ) W/kg, which was one order of magnitude higher than the previously observed ε in the upper ocean of the South China Sea slope. At the same observation site, the internal tides induced ε in the bottom boundary layer was O (10 −5 ) W/kg, which indicated that the internal tides played a more important role than the ISWs in enhancing bottom turbulent mixing. • Direct observation reveals strong internal solitary waves near seabed of slope. • Solitary waves cause strong velocity (0.2 m/s) near seabed at depth of 659 m. • Internal solitary waves increase bottom layer turbulent mixing several times.

Internal Hydraulic Transition and Turbulent Mixing Observed in the Kuroshio over the I-Lan Ridge off Northeastern Taiwan

Abstract Complex small-scale processes and energetic turbulence are observed at a sill located on the I-Lan Ridge that spans across the strong Kuroshio off Taiwan. The current speed above the sill is strong (1.5 m s−1) and unsteady (±0.5 m s−1) due to the Kuroshio being modulated by the semidiurnal tide. Above the sill crest, isothermal domes, with vertical scales of ∼20 and ∼50 m during the low and high tides, respectively, are generated by turbulent mixing as a result of shear instability in the bottom boundary layer. Tidally modulated hydraulic character modifies the small-scale processes occurring on the leeward side of the sill. Criticality analysis, performed by solving the Taylor–Goldstein equation, suggests that the observed lee waves and intermediate layer sandwiched by two free shear layers are related to the mode-1 and mode-2 critical control between the sill crest and immediate lee, respectively. Around high tide, lee waves are advected further downstream, and only mode-1 critical control can occur, leading to a warm water depression. The shear instabilities ensuing from the hydraulic transition processes continuously mediate flow kinetic energy to turbulence such that the status of marginal instability where the Richardson number converges at approximately 0.25 is reached. The resultant eddy diffusivity Kρ is concentrated from O(10−4) to O(10−3) m2 s−1 and has a maximum value of 10 m2 s−1. The sill on the western flank of the Kuroshio is a hotspot for energetic mixing of Kuroshio waters and South China Sea waters.

Effects of Nonlinearity on Velocity, Acceleration and Pressure Gradient in Free-Stream Zone of Solitary Wave over Horizontal Bed—An Experimental Study

For solitary waves on a horizontal bed, the study deals experimentally with the high ratio of wave height (H0) to still water depth (h0) that amplifies the wave nonlinearity. The value of H0/h0 tested in a wave flume ranges from 0.050 to 0.550, indicating the shift from a quasi-linear solitary wave to a highly nonlinear one. A high-speed particle image velocimetry (HSPIV) and a flow visualization technique of particle-trajectory tracking method are utilized to measure velocity fields and identify near-bed flow structures. The unsteady free-stream velocities with equal magnitude take place in a free-stream zone, FSZ). The FSZ underlies the internal flow zone, over which the external free surface of solitary wave exists and is situated beyond the boundary layer. The spatio-temporal variation of free-stream velocity, moving in phase with the free surface elevation, characterizes the pattern of pressure gradient in the FSZ and thus dominates the behavior of boundary layer flow. Accordingly, nonlinear effects on the time series as well as the maximum values of horizontal velocity, particle acceleration, and pressure gradient in the FSZs of solitary waves are presented. Before, at, and after the wave crest’s intersection with a given measurement location, favorable, zero, and adverse pressure gradients occur in the FSZ, respectively. For H0/h0 = 0.179, 0.363, and 0.550, the values of the dimensionless maximum free-stream velocity are about 3.10, 5.32, and 6.20 times that (= 0.0473) for H0/h0 = 0.050; and the corresponding values of the dimensionless maximum adverse pressure gradient are about 5.74, 14.54 and 19.84 times that (= 0.0061) for H0/h0 = 0.050. This evidence highlights the nonlinear effect on the kinematic and hydrodynamic features of solitary waves. Finally, the effect of nonlinearity on the relationship between the dimensionless time for the maximum adverse pressure gradient in the FSZ and that for the incipient flow reversal in the bottom boundary layer is explored for the first time. It is found that the incipient flow reversal takes place immediately after the maximum adverse pressure gradient, together with a decrease in the dimensionless time for flow reversal if H0/h0 increases. The fact accentuates the nonlinear effect on the incipient flow reversal right above the bed.

Adaptation of the SANTOSS transport formula for 3D nearshore models: Application to cross-shore sandbar migration

The ability to model bar migration in the nearshore zone is critically important for predicting coastal morphology. We present a cost-effective wave-averaged 3D model (CROCO) with parameterization of wave-related bedload transport (adapted SANTOSS), and evaluate/calibrate its performance in comparison with data measured in a two-dimensional wave channel. In this model, suspended-load transport (both in and above the wave bottom boundary layer) is resolved by the flow model. We assume that if the tuning parameters are not robust to varying forcing conditions, the parameterization can be considered a failure. The present work is therefore a first step towards a numerical model capable of predicting the onshore and offshore migration of a sandbar under storm and post-storm conditions using the same set of parameters. One way to achieve this goal is to treat the effects of waves and currents in separate formulations to avoid conflicts and redundancy. We show that the implementation and adaptation of SANTOSS in CROCO can achieve this goal, provided that the hydrodynamics are accurately reproduced and the separation between the wave-related bedload and current-related suspended load formulations is respected.

Periodic boundary layer separation and lateral intrusions observed above a sloping lakebed

AbstractLateral intrusions flowing from the boundaries of lakes and oceans might influence basin‐scale distributions of heat, chemical solutes, sediments, and organisms. Here, observations of numerous intrusions from multiple locations in a small lake were used to examine the relation of intrusions to a lakewide internal seiche, and to estimate mean exchange between sloping bottom boundary layers (BBLs) and the lake interior. BBL waters were observed periodically separating from the bed and flowing laterally offshore. Boundary layer separation was initiated where strong downslope flows encountered upslope‐propagating temperature fronts. These separation events propagated upslope, coherent with the vertically propagating internal seiche, resembling previously reported cases of upslope‐propagating internal bores. Following separation, jets flowed offshore immediately above temperature fronts, transporting boundary layer water at least 61 m into the interior. Jet propagation could be traced upwards, from the bed at 9.4 m depth to the base of the surface mixed layer. Mean offshore transport was most intense within the metalimnion (4–7 m depth). Intruding intermediate‐temperature thermocline water was likely supplied both by river inputs and by boundary mixing of cold and warm waters. These findings suggest that periodic jets, resulting from seiche‐induced boundary layer separation, transported significant quantities of water from the boundary layer toward the stratified lake interior, with possible implications for distributions of sediments, solutes, and organisms.

A new process-based, wave-resolving, 2DH circulation model for the evolution of natural sand bars: The role of nearbed dynamics and suspended sediment transport

We study the migration of natural sand bars that evolve in the nearshore by means of a new process-based wave-resolving, 2DH circulation model. In order to perform reliable and accurate computations, the robust Nonlinear Shallow Water Equations (NSWEs) hydro-morphodynamic solver of Brocchini et al. (2001) and Postacchini et al. (2012) implements a detailed description of the Bottom Boundary Layer (BBL) dynamics and a new predictor for the Suspended Sediment Transport (SST) based on the solution of a Depth-Averaged Advection-Diffusion Equation (DAADE) for the sediment concentration. The robustness and accuracy of the enhanced model are validated against literature theoretical, experimental, and numerical results, all comparisons highlighting good performances and clarifying the role of both BBL and SST contributions, the former one having a larger positive influence than the latter one on the results. Both original and enhanced models are, then, used to predict the evolution of the sand bar system that characterizes the nearshore of Senigallia (AN). The analysis leverages the field observations collected at such a site by means of the Sena Gallica Speculator video-monitoring system. Modeling of the storm-forced sand bar migration patterns reveals that: 1) the enhanced model can adequately reproduce the seaward migration of the sand bars of the system; 2) the process of shoreline retreat in coincidence with the generation of a new-born bar is well described; 3) inclusion of the BBL improves quantitative prediction of the bar crest migration; 4) the SST, beyond improving the prediction of the bar crest location, induces some smoothing of the bar profile, in line with the literature findings of SST being a stabilizing factor for the bar emergence. • Sand bar evolution in the nearshore zone is modeled using a process-based 2DH wave resolving NSWE model. • The model includes BBL dynamics and predicts the SST based on the depth-averaged advection-diffusion equation. • The model is validated against theoretical, experimental, and numerical results. • The model is used to predict the sandbar migration in the nearshore of Senigallia (AN) during a storm. • The developed model can accurately reproduce the sandbar evolution in both field and laboratory conditions.

Sensing of ocean mixing via forward scattering acoustic scintillation

Ocean mixing process observations show complex spatio-temporal patterns as well as energy and dissipation levels that vary over six or more orders of magnitude, creating a challenging sampling problem. Acoustic imaging of ocean turbulent structures is informative, but extraction of quantities like spectra and dissipation rates of kinetic energy, ε, and temperature variance, Χ, requires geometric sampling constraints. Here, a remote sensing approach is proposed in the spirit of Di Iorio and Farmer (1994), and Duda et al. (1988), based on weak multiple-forward scattering physics. In the correct scattering regime, observations of the wavenumber spectrum of acoustic intensity obtained on a vertical array gives a direct measure of the range average turbulence index of the refraction spectrum in a plane transverse to the acoustic path. In this talk, the theory will be reviewed and some illustrative simulations will be presented; importantly, one-way acoustic transmission measures a combination of both ε and Χ. Experimental arrangements (frequency/range) where the theory is expected to apply will be outlined with an eye towards observing bottom boundary layer mixing. With a suitable vertical array, one could resolve scales from the smallest internal waves (∼several meters), through the transition regime, and down to scales of turbulence (∼tens of cm).

Impacts of Instrumented Bottom Frame on Flow and Turbulence Measurements

Abstract A series of laboratory experiments are carried out to demonstrate the impacts of instrumented bottom frame legs on flow and turbulence. The magnitudes of vertical velocity, turbulent kinetic energy, dissipation, and shear stress induced by the frame legs depend on several factors, including the diameter and number of the frame legs, distances between the legs and the observational location, and the magnitude of the incoming flow and its direction with respect to the layout of the frame. In situ observations were carried out near the mouth of the Yellow River using two acoustic Doppler velocimeters mounted on a bottom frame. The estimated vertical velocity and turbulent kinetic energy dissipation rate show a significant asymmetry with flood and ebb tidal flows. This asymmetry can be partly explained by the influences of the bottom frame legs. Finally, the design and deployment principles of bottom frames are discussed for the purpose of reducing the impacts of the frame legs. Significance Statement Instrumented bottom frames are widely used for observations in the oceanic bottom boundary layer and above. However, the impacts of frame legs on the observed flow and turbulence have rarely been investigated. A series of laboratory experiments demonstrate that frame legs can induce vertical flow and enhanced turbulence, and the magnitudes of these influences vary with the size and layout of the frame legs and the magnitude and direction of the background flow. The results of the laboratory experiments can partially explain an “asymmetry” behavior of the vertical flow and turbulent kinetic energy with the flood and ebb tidal flows, derived from in situ observations near the mouth of the Yellow River.

Influence of Stratification and Bottom Boundary Layer on the Classical Ekman Model

A depth understanding of the different processes of water movements produced by the wind surface stress yields a better description and improvement of the marine food chain and ecosystem. The classical Ekman model proposes a hypothetical ocean, excluding the influence of continents and the Coriolis force. It also assumes infinite depth and a constant vertical eddy viscosity. The current study aims to understand how the vertical velocity profile is affected by the variation of the eddy viscosity coefficient (kz) and the consideration of a finite depth. The study uses an ideal analytical model with the Ekman classical model as a starting point. It has been demonstrated that, for a very stratified profile, when the depth is not considered infinity, the Ekman transport tends to a direction smaller than 80°. It differs from the classical Ekman model, which proposes an approximated angle equal to 90°. Considering the modified model, it was also found that the surface current deviation is smaller than 40°, which differs from the 45° proposed by the classical model. In addition, it was determined that for ocean depths smaller than 180 m, the maximum velocity does not occur at the water surface, as in the classical model, but at deeper levels.

Evolving mantle convection from bottom up to top down

When it comes to convection, what goes up must come down. Or is it, what goes down must come up? The truth is it depends. Although convection must be mass balanced, there is no reason that it must be force balanced: the positive and negative buoyancy forces driving convection up and down, respectively, do not necessarily need to be balanced. The balance, or imbalance, all depends on the top and bottom boundary layers. Thus, convection in Earth's mantle depends on the temperature differences across the core-mantle boundary below and the lithosphere-asthenosphere boundary above. Convective asymmetry predominated by positive buoyancy, or bottom-up convection, would be driven by plume ascent, whereas if it were predominated by negative buoyancy, or top-down convection, it would be driven by plate subduction. Symmetric convection would balance plume ascent and plate subduction. Is mantle convection on Earth balanced, dominantly top down or bottom up, or time dependent?

Particle trajectories in an eastern boundary current using a regional ocean model at two horizontal resolutions

Lagrangian particle tracking (LPT) models are used to study the transport and dispersal of marine organisms. In LPT studies, the accuracy of the circulation is essential for nearshore habitats of Eastern Boundary Current (EBC) regions that are areas of high productivity and economically important fisheries. We used the California Current System as an example of an EBC region, specifically the Oregon coast located in the northern California Current System because it has distinct upwelling and downwelling regimes and variable shelf width. More specifically, we developed and applied a LPT model to compare and contrast particle drift patterns during the spring transition as it is an important period for spawning. We contrasted years (2016–18) using Regional Ocean Modeling System (ROMS) with different horizontal spatial resolutions (2 km, 250 m). Lagrangian particles experience stronger downward velocities and displacements to greater depths in the 250 m ROMS simulations that used a finer resolution bathymetry. Consequently, retention along the Oregon coast increases in the 250 m ROMS compared to the 2 km ROMS. After 10 days, 37%–83% of particles forced with the 2 km ROMS remain in the model domain, compared to 61%–86% of particles remaining when using the 250 m ROMS. Particles in the 250 m ROMS are advected to depth at specific times and locations for each simulated year, coinciding with the location and timing of a strong and shallow alongshore undercurrent that is not present in the 2 km ROMS. Additionally, ageostrophic dynamics close to shore, in the bottom boundary layer, and around headlands emerge in the 250 m resolution model, while they are at best poorly resolved in the 2 km resolution case. We conclude that the higher horizontal model resolution and bathymetry used in the 250 m ROMS generates well-resolved mesoscale and submesoscale features (e.g., surface, subsurface, and nearshore jet) that vary annually. These physical features are significantly different than those modeled by the 2 km model and may be responsible for these differences in particle dispersal. These results have implications for modeling the dispersal, growth, and development of coastal organisms with dispersing early life stages. • We compare ROMS of different horizontal resolutions in particle tracking models. • The higher resolution ROMS (250 m) better resolves mesoscale and submesoscale features. • The resolved mesoscale and submesoscale features may be due to higher resolution and bathymetry. • An increase in particle retention is linked to more resolved mesoscale features.

Insights into the Mixing Efficiency of Submesoscale Centrifugal–Symmetric Instabilities

Abstract Submesoscale processes provide a pathway for energy to transfer from the balanced circulation to turbulent dissipation. One class of submesoscale phenomena that has been shown to be particularly effective at removing energy from the balanced flow is centrifugal–symmetric instabilities (CSIs), which grow via geostrophic shear production. CSIs have been observed to generate significant mixing in both the surface boundary layer and bottom boundary layer flows along bathymetry, where they have been implicated in the mixing and water mass transformation of Antarctic Bottom Water. However, the mixing efficiency (i.e., the fraction of the energy extracted from the flow used to irreversibly mix the fluid) of these instabilities remains uncertain, making estimates of mixing and energy dissipation due to CSI difficult. In this work we use large-eddy simulations to investigate the mixing efficiency of CSIs in the submesoscale range. We find that centrifugally dominated CSIs (i.e., CSI mostly driven by horizontal shear production) tend to have a higher mixing efficiency than symmetrically dominated ones (i.e., driven by vertical shear production). The mixing efficiency associated with CSIs can therefore alternately be significantly higher or significantly lower than the canonical value used by most studies. These results can be understood in light of recent work on stratified turbulence, whereby CSIs control the background state of the flow in which smaller-scale secondary overturning instabilities develop, thus actively modifying the characteristics of mixing by Kelvin–Helmholtz instabilities. Our results also suggest that it may be possible to predict the mixing efficiency with more readily measurable parameters (viz., the Richardson and Rossby numbers), which would allow for parameterization of this effect.

A machine learning approach to predicting equilibrium ripple wavelength

Sand ripples are geomorphic features on the seafloor that affect bottom boundary layer dynamics including wave attenuation and sediment transport. We present a new equilibrium ripple predictor using a machine learning approach that outputs a probability distribution of wave-generated equilibrium wavelengths and statistics including an estimate of ripple height, the most probable ripple wavelength, and sediment and flow parameterizations. The Bayesian Optimal Model System (BOMS) is an ensemble machine learning system that combines two machine learning algorithms and two deterministic empirical ripple predictors with a Bayesian meta-learner to produce probabilistic wave-generated equilibrium ripple wavelength estimates in sandy locations. A ten-fold cross validation of BOMS resulted in an adjusted R-squared value of 0.93 and an average root mean square error (RMSE) of 8.0 cm. During both cross validation and testing on three unique field datasets, BOMS provided more accurate wavelength predictions than each individual base model and other common ripple predictors.

Microstructure observations of three-layered mixing on the shelf of the northwestern South China Sea

This study presents microstructure observations on the northwestern South China Sea (nwSCS) shelf, where energetic turbulence drove a high dissipation rate on the order of 10−5 W kg−1 in May 2021. Together with the hydrography and current observations, the results show a typical three-layered structure of dissipation on the shelf and reveal different mixing regimes for different layers. Below the mixed layer, mixing over the shallow shelf displays periodically varying characteristics, indicating low-mode internal tides. Strong mixing occurred in the stably stratified layer, reaching an approximate magnitude of 1 × 10−5 W kg−1. Our analysis suggests that the strong mixing in this layer is attributed to the shear behavior of the low-mode internal waves because the observed dissipation rate increases with increasing vertical shear and stratification as described by the MacKinnon and Gregg parameterization. In the bottom boundary layer, mixing by tidal shear is dependent on the Richardson number and trends with velocity and vertical shear. In addition, our analysis also reveals that mixing continuously contributed to downward heat transport to the deep layers. This study can improve our understanding of the vertical structure of fine-scale mixing and its regimes, and their effects on the shelf of the nwSCS.

Diapycnal displacement, diffusion, and distortion of tracers in the ocean

Abstract Small-scale mixing drives the diabatic upwelling that closes the abyssal ocean overturning circulation. Indirect microstructure measurements of in-situ turbulence suggest that mixing is bottom-enhanced over rough topography, implying downwelling in the interior and stronger upwelling in a sloping bottom boundary layer. Tracer Release Experiments (TREs), in which inert tracers are purposefully released and their dispersion is surveyed over time, have been used to independently infer turbulent diffusivities—but typically provide estimates in excess of microstructure ones. In an attempt to reconcile these differences, Ruan and Ferrari (2021) derived exact tracer-weighted buoyancy moment diagnostics, which we here apply to quasi-realistic simulations. A tracer’s diapycnal displacement rate is exactly twice the tracer-averaged buoyancy velocity, itself a convolution of an asymmetric upwelling/downwelling dipole. The tracer’s diapycnal spreading rate, however, involves both the expected positive contribution from the tracer-averaged in-situ diffusion as well as an additional non-linear diapycnal distortion term, which is caused by correlations between buoyancy and the buoyancy velocity, and can be of either sign. Distortion is generally positive (stretching) due to bottom-enhanced mixing in the stratified interior but negative (contraction) near the bottom. Our simulations suggest that these two effects coincidentally cancel for the Brazil Basin Tracer Release Experiment, resulting in negligible net distortion. By contrast, near-bottom tracers experience leading-order distortion that varies in time. Errors in tracer moments due to realistically sparse sampling are generally small (< 20%), especially compared to the O(1) structural errors due to the omission of distortion effects in inverse models. These results suggest that TREs, although indispensable, should not be treated as “unambiguous” constraints on diapycnal mixing.