Search for: All records

Abstract We describe a process called “squeeze dispersion” in which the squeezing of oceanic tracer gradients by waves, eddies, and bathymetric flow modulates diapycnal diffusion by centimeter to meter‐scale turbulence. Due to squeeze dispersion, the effective diapycnal diffusivity of oceanic tracers is different and typically greater than the average “local” diffusivity, especially when local diffusivity correlates with squeezing. We develop a theory to quantify the effects of squeeze dispersion on diapycnal oceanic transport, finding formulas that connect density‐averaged tracer flux, locally measured diffusivity, large‐scale oceanic strain, the thickness‐weighted average buoyancy gradient, and the effective diffusivity of oceanic tracers. We use this effective diffusivity to interpret observations of abyssal flow through the Samoan Passage reported by Alford et al. (2013,https://doi.org/10.1002/grl.50684) and find that squeezing modulates diapycnal tracer dispersion by factors between 0.5 and 3. 

Abstract The Indo-Pacific Ocean appears exponentially stratified between 1- and 3-km depth with a decay scale on the order of 1 km. In his celebrated paper “Abyssal recipes,” W. Munk proposed a theoretical explanation of these observations by suggesting a pointwise buoyancy balance between the upwelling of cold water and the downward diffusion of heat. Assuming a constant upwelling velocity w and turbulent diffusivity κ , the model yields an exponential stratification whose decay scale is consistent with observations if κ ∼ 10 −4 m 2 s −1 . Over time, much effort has been made to reconcile Munk’s ideas with evidence of vertical variability in κ , but comparably little emphasis has been placed on the even stronger evidence that w decays toward the surface. In particular, the basin-averaged w nearly vanishes at 1-km depth in the Indo-Pacific. In light of this evidence, we consider a variable-coefficient, basin-averaged analog of Munk’s budget, which we verify against a hierarchy of numerical models ranging from an idealized basin-and-channel configuration to a coarse global ocean simulation. Study of the budget reveals that the decay of basin-averaged w requires a concurrent decay in basin-averaged κ to produce an exponential-like stratification. As such, the frequently cited value of 10 −4 m 2 s −1 is representative only of the bottom of the middepths, whereas κ must be much smaller above. The decay of mixing in the vertical is as important to the stratification as its magnitude . Significance Statement Using a combination of theory and numerical simulations, it is argued that the observed magnitude and shape of the global ocean stratification and overturning circulation appear to demand that turbulent mixing increases quasi-exponentially toward the ocean bottom. Climate models must therefore prescribe such a vertical profile of turbulent mixing in order to properly represent the heat and carbon uptake accomplished by the global overturning circulation on centennial and longer time scales. 

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. 

Abstract ABSTRACT: The abyssal overturning circulation is thought to be primarily driven by small-scale turbulent mixing. Diagnosed watermass transformations are dominated by rough topography “hotspots”, where the bottom-enhancement of mixing causes the diffusive buoyancy flux to diverge, driving widespread downwelling in the interior—only to be overwhelmed by an even stronger up-welling in a thin Bottom Boundary Layer (BBL). These watermass transformations are significantly underestimated by one-dimensional (1D) sloping boundary layer solutions, suggesting the importance of three-dimensional physics. Here, we use a hierarchy of models to generalize this 1D boundary layer approach to three-dimensional eddying flows over realistically rough topography. When applied to the Mid-Atlantic Ridge in the Brazil Basin, the idealized simulation results are roughly consistent with available observations. Integral buoyancy budgets isolate the physical processes that contribute to realistically strong BBL upwelling. The downwards diffusion of buoyancy is primarily balanced by upwelling along the sloping canyon sidewalls and the surrounding abyssal hills. These flows are strengthened by the restratifying effects of submesoscale baroclinic eddies and by the blocking of along-ridge thermal wind within the canyon. Major topographic sills block along-thalweg flows from restratifying the canyon trough, resulting in the continual erosion of the trough’s stratification. We propose simple modifications to the 1D boundary layer model which approximate each of these three-dimensional effects. These results provide local dynamical insights into mixing-driven abyssal overturning, but a complete theory will also require the non-local coupling to the basin-scale circulation. 

Abstract Turbulent mixing across density surfaces transforms abyssal ocean waters into lighter waters and is vital to close the deepest branches of the global overturning circulation. Over the last 20 years, mixing rates inferred from in situ microstructure profilers and tracer release experiments (TREs) have provided valuable insights in the connection between small-scale mixing and large-scale ocean circulation. Problematically, estimates based on TREs consistently exceed those from collocated in situ microstructure measurements. These differences have been attributed to a low bias in the microstructure estimates that can miss strong, but rare, mixing events. Here we demonstrate that TRE estimates can suffer from a high bias, because of the approximations generally made to interpret the data. We first derive formulas to estimate mixing from the temporal growth of the second moment of a tracer patch by extending Taylor’s celebrated formula to account for both density stratification and variations in mixing rates. The formulas are validated with tracers released in numerical simulations of turbulent flows and then used to discuss biases in the interpretation of TREs based estimates and how to possibly overcome them. 

Abstract The emerging view of the abyssal circulation is that it is associated with bottom-enhanced mixing, which results in downwelling in the stratified ocean interior and upwelling in a bottom boundary layer along the insulating and sloping seafloor. In the limit of slowly varying vertical stratification and topography, however, boundary layer theory predicts that these upslope and downslope flows largely compensate, such that net water mass transformations along the slope are vanishingly small. Using a planetary geostrophic circulation model that resolves both the boundary layer dynamics and the large-scale overturning in an idealized basin with bottom-enhanced mixing along a midocean ridge, we show that vertical variations in stratification become sufficiently large at equilibrium to reduce the degree of compensation along the midocean ridge flanks. The resulting large net transformations are similar to estimates for the abyssal ocean and span the vertical extent of the ridge. These results suggest that boundary flows generated by mixing play a crucial role in setting the global ocean stratification and overturning circulation, requiring a revision of abyssal ocean theories. 

Abstract Closing the overturning circulation of bottom water requires abyssal transformation to lighter densities and upwelling. Where and how buoyancy is gained and water is transported upward remain topics of debate, not least because the available observations generally show downward-increasing turbulence levels in the abyss, apparently implying mean vertical turbulent buoyancy-flux divergence (densification). Here, we synthesize available observations indicating that bottom water is made less dense and upwelled in fracture zone valleys on the flanks of slow-spreading midocean ridges, which cover more than one-half of the seafloor area in some regions. The fracture zones are filled almost completely with water flowing up-valley and gaining buoyancy. Locally, valley water is transformed to lighter densities both in thin boundary layers that are in contact with the seafloor, where the buoyancy flux must vanish to match the no-flux boundary condition, and in thicker layers associated with downward-decreasing turbulence levels below interior maxima associated with hydraulic overflows and critical-layer interactions. Integrated across the valley, the turbulent buoyancy fluxes show maxima near the sidewall crests, consistent with net convergence below, with little sensitivity of this pattern to the vertical structure of the turbulence profiles, which implies that buoyancy flux convergence in the layers with downward-decreasing turbulence levels dominates over the divergence elsewhere, accounting for the net transformation to lighter densities in fracture zone valleys. We conclude that fracture zone topography likely exerts a controlling influence on the transformation and upwelling of bottom water in many areas of the global ocean.