Search for: All records

Abstract Diapycnal upwelling along sloping topography has been shown to be an important component of the abyssal overturning circulation. Theoretical studies of mixing-driven upwelling have mostly relied on a time-averaged description of mixing acting on a mean stratification which ignores the intermittency of mixing. Typically, these studies prescribed a time-invariant turbulent diffusivity profile motivated by scenarios where tidal currents encounter gentle topography with small-scale corrugations, leading to subsequent propagation and breaking of internal waves. Here, a different scenario is considered where a tidal current interacts with smooth but steep topography, a case often encountered near continental margins and troughs. The performed nonhydrostatic simulations resolve both the strong oscillatory shear that develops along the steep critical topography and the associated mixing events. Strong diapycnal mixing is observed during the upslope phase of the tidal flow when both the near-boundary stratification and shear are enhanced. During the downslope phase, strong overturning events do develop, but they are associated with weak stratification and less efficient diapycnal mixing. These results highlight that the temporal evolution of both shear and stratification play a key role in setting when diapycnal mixing and water mass transformation occur along steep topography. In contrast, over gentle topography, tidal shears do not become sufficiently large to generate strong local mixing for typical oceanographic parameters. 

Abstract Coastal upwelling systems play a key role in sustaining productive coastal ecosystems in the global ocean by transporting nutrients to surface waters. However, the fundamental mechanisms and pathways responsible for nutrient upwelling are not fully understood, largely due to the historically employed two-dimensional frameworks in which coastal upwelling systems have long been studied. Using both observations and idealized numerical simulations, we identify and quantify two primary routes of nutrient upwelling: the residual circulation, resulting from a significant cancellation between Eulerian-mean and eddy-induced circulations, and along-isopycnal eddy stirring. Our analysis demonstrates that their relative contributions depend on two distinct parameters: 1) the slope Burger number S, defined here as S=αN/f, whereαis the topographic slope angle and Nandfare the buoyancy and Coriolis frequencies, and 2) the surface nutrient uptake rate by biological activities. Specifically, we propose that wind forcing induces isopycnal tilting and surface outcropping, which creates favorable conditions for along-isopycnal nutrient gradients to develop in regions of strong biological activity at the surface. The magnitude of these gradients depends on both the slope Burger number S, which influences the strength of the residual circulation bringing nutrients from depths, and the surface biological uptake rate, which consumes nutrients. Our diagnostics provide insights into the intricate pathways for nutrient upwelling and underscore the significance of eddy stirring in coastal upwelling systems. 

Abstract Small-scale turbulent mixing drives the upwelling of deep water masses in the abyssal ocean as part of the global overturning circulation¹. However, the processes leading to mixing and the pathways through which this upwelling occurs remain insufficiently understood. Recent observational and theoretical work^2–5has suggested that deep-water upwelling may occur along the ocean’s sloping seafloor; however, evidence has, so far, been indirect. Here we show vigorous near-bottom upwelling across isopycnals at a rate of the order of 100 metres per day, coupled with adiabatic exchange of near-boundary and interior fluid. These observations were made using a dye released close to the seafloor within a sloping submarine canyon, and they provide direct evidence of strong, bottom-focused diapycnal upwelling in the deep ocean. This supports previous suggestions that mixing at topographic features, such as canyons, leads to globally significant upwelling^3,6–8. The upwelling rates observed were approximately 10,000 times higher than the global average value required for approximately 30 × 10⁶m³s⁻¹of net upwelling globally⁹. 

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. 

The dissipation of the kinetic energy (KE) associated with oceanic flows is believed to occur primarily in the oceanic bottom boundary layer (BBL), where bottom drag converts the KE from mean flows to heat loss through irreversible mixing at molecular scales. Due to the practical difficulties associated with direct observations on small-scale turbulence close to the seafloor, most up-to-date estimates on bottom drag rely on a simple bulk formula (CdU3) proposed by G.I. Taylor that relates the integrated BBL dissipation rate to a drag coefficient (Cd) as well as a flow magnitude outside of the BBL (U). Using output from several turbulence-resolving direct numerical simulations, it is shown that the true BBL-integrated dissipation rate is approximately 90% of that estimated using the classic bulk formula, applied here to the simplest scenario where a mean flow is present over a flat and hydrodynamically smooth bottom. It is further argued that Taylor’s formula only provides an upper bound estimate and should be applied with caution in the future quantification of BBL dissipation; the performance of the bulk formula depends on the distribution of velocity and shear stress near the bottom, which, in the real ocean, could be disrupted by bottom roughness. 

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 To close the abyssal overturning circulation, dense bottom water has to become lighter by mixing with lighter water above. This diapycnal mixing is strongly enhanced over rough topography in abyssal mixing layers, which span the bottom few hundred meters of the water column. In particular, mixing rates are enhanced over mid-ocean ridge systems, which extend for thousands of kilometers in the global ocean and are thought to be key contributors to the required abyssal water mass transformation. To examine how stratification and thus diabatic transformation is maintained in such abyssal mixing layers, this study explores the circulation driven by bottom-intensified mixing over mid-ocean ridge flanks and within ridge-flank canyons. Idealized numerical experiments show that stratification over the ridge flanks is maintained by submesoscale baroclinic eddies and that stratification within ridge-flank canyons is maintained by mixing-driven mean flows. These restratification processes affect how strong a diabatic buoyancy flux into the abyss can be maintained, and they are essential for maintaining the dipole in water mass transformation that has emerged as the hallmark of a diabatic circulation driven by bottom-intensified mixing. 

ABSTRACT Meltwater is important to understanding glacier health and dynamics. Since melt measurements are uncommon, ice ablation estimates are often based on models including the positive degree day (PDD) model. The PDD estimate is popular since it only requires air temperature as input, but suffers from the lack of physical motivation of an energy-balance model. We present a physics-based alternative to the PDD model that still only takes air/surface temperature as input. The model resembles the PDD model except accounting for time lags in ablation when cold ice needs to be warmed. The model is expressed as a differential equation with a single extra parameter related to the efficiency of heating a near-surface layer of ice. With zero thickness, the model reduces to the PDD model, providing a physical basis for the PDD model. Applying the model to data from Greenland, it improves modestly upon the PDD model, with the main improvement being better prediction of early season melting. This new model is a useful compromise, with some of the physics of more realistic models and the simplicity of a PDD model. The model should improve estimates of meltwater production and help constrain PDD parameters when empirical calibration is challenging. 

Abstract Bottom drag is believed to be one of the key mechanisms that remove kinetic energy from the ocean's general circulation. However, large uncertainty still remains in global estimates of bottom drag dissipation. One significant source of uncertainty comes from the velocity structures near the bottom where the combination of sloping topography and stratification can reduce the mean flow magnitude, and thus the bottom drag dissipation. Using high‐resolution numerical simulations, we demonstrate that previous estimates of bottom drag dissipation are biased high because they neglect velocity shear in the bottom boundary layer. The estimated bottom drag dissipation associated with geostrophic flows over the continental slopes is at least 56% smaller compared with prior estimates made using total velocities outside the near‐bottom layer. The diagnostics suggest the necessity of resolving the bottom boundary layer structures in coarse‐resolution ocean models and observations in order to close the global kinetic energy budget.