An official website of the United States government
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.
more »
« less
Evolution of Bottom Boundary Layers on Three Dimensional Topography—Buoyancy Adjustment and Instabilities
Abstract A current along a sloping bottom gives rise to upwelling, or downwelling Ekman transport within the stratified bottom boundary layer (BBL), also known as the bottom Ekman layer. In 1D models of slope currents, geostrophic vertical shear resulting from horizontal buoyancy gradients within the BBL is predicted to eventually bring the bottom stress to zero, leading to a shutdown, or “arrest,” of the BBL. Using 3D ROMS simulations, we explore how the dynamics of buoyancy adjustment in a current‐ridge encounter problem differs from 1D and 2D temporal spin up problems. We show that in a downwelling BBL, the destruction of the ageostrophic BBL shear, and hence the bottom stress, is accomplished primarily by nonlinear straining effects during the initial topographic encounter. As the current advects along the ridge slopes, the BBL deepens and evolves toward thermal wind balance. The onset of negative potential vorticitymodes of instability and their subsequent dissipation partially offsets the reduction of the BBL dissipation during the ridge‐current interaction. On the upwelling side, although the bottom stress weakens substantially during the encounter, the BBL experiences a horizontal inflectional point instability and separates from the slopes before sustained along‐slope stress reduction can occur. In all our solutions, both the upwelling and downwelling BBLs are in a partially arrested state when the current separates from the ridge slope, characterized by a reduced, but non‐zero bottom stress on the slopes.
more »
« less
- Award ID(s):
- 1751386
- PAR ID:
- 10409218
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Oceans
- Volume:
- 128
- Issue:
- 4
- ISSN:
- 2169-9275
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract This study examines the utility of Eady-type theories as applied to understanding baroclinic instability in coastal flows where depth variations and bottom drag are important. The focus is on the effects of nongeostrophy, boundary dissipation, and bottom slope. The approach compares theoretically derived instability properties against numerical model calculations, for experiments designed to isolate the individual effects and justified to have Eady-like basic states. For the nongeostrophic effect, the theory of Stone (1966) is shown to give reasonable predictions for the most unstable growth rate and wavelength. It is also shown that the growing instability in a fully nonlinear model can be interpreted as boundary-trapped Rossby wave interactions—that is, wave phase locking and westward phase tilt allow waves to be mutually amplified. The analyses demonstrate that both the boundary dissipative and bottom slope effects can be represented by vertical velocities at the lower boundary of the unstable interior, via inducing Ekman pumping and slope-parallel flow, respectively, as proposed by the theories of Williams and Robinson (1974; referred to as the Eady–Ekman problem) and Blumsack and Gierasch (1972). The vertical velocities, characterized by a friction parameter and a slope ratio, modify the bottom wave and thus the scale selection. However, the theories have inherent quantitative limitations. Eady–Ekman neglects boundary layer responses that limit the increase of bottom stress, thereby overestimating the Ekman pumping and growth rate reduction at large drag. Blumsack and Gierasch’s (1972) model ignores slope-induced horizontal shear in the mean flow that tilts the eddies to favor converting energy back to the mean, thus having limited utility over steep slopes.more » « less
-
Abstract Realistic computational simulations in different oceanic basins reveal prevalent prograde mean flows (in the direction of topographic Rossby wave propagation along isobaths; aka topostrophy) on topographic slopes in the deep ocean, consistent with the barotropic theory of eddy-driven mean flows. Attention is focused on the western Mediterranean Sea with strong currents and steep topography. These prograde mean currents induce an opposing bottom drag stress and thus a turbulent boundary layer mean flow in the downhill direction, evidenced by a near-bottom negative mean vertical velocity. The slope-normal profile of diapycnal buoyancy mixing results in downslope mean advection near the bottom (a tendency to locally increase the mean buoyancy) and upslope buoyancy mixing (a tendency to decrease buoyancy) with associated buoyancy fluxes across the mean isopycnal surfaces (diapycnal downwelling). In the upper part of the boundary layer and nearby interior, the diapycnal turbulent buoyancy flux divergence reverses sign (diapycnal upwelling), with upward Eulerian mean buoyancy advection across isopycnal surfaces. These near-slope tendencies abate with further distance from the boundary. An along-isobath mean momentum balance shows an advective acceleration and a bottom-drag retardation of the prograde flow. The eddy buoyancy advection is significant near the slope, and the associated eddy potential energy conversion is negative, consistent with mean vertical shear flow generation for the eddies. This cross-isobath flow structure differs from previous proposals, and a new one-dimensional model is constructed for a topostrophic, stratified, slope bottom boundary layer. The broader issue of the return pathways of the global thermohaline circulation remains open, but the abyssal slope region is likely to play a dominant role.more » « less
-
Turbulence and mixing in a near-bottom convectively driven flow are examined by numerical simulations of a model problem: a statically unstable disturbance at a slope with inclination $$\unicode[STIX]{x1D6FD}$$ in a stable background with buoyancy frequency $$N$$ . The influence of slope angle and initial disturbance amplitude are quantified in a parametric study. The flow evolution involves energy exchange between four energy reservoirs, namely the mean and turbulent components of kinetic energy (KE) and available potential energy (APE). In contrast to the zero-slope case where the mean flow is negligible, the presence of a slope leads to a current that oscillates with $$\unicode[STIX]{x1D714}=N\sin \unicode[STIX]{x1D6FD}$$ and qualitatively changes the subsequent evolution of the initial density disturbance. The frequency, $$N\sin \unicode[STIX]{x1D6FD}$$ , and the initial speed of the current are predicted using linear theory. The energy transfer in the sloping cases is dominated by an oscillatory exchange between mean APE and mean KE with a transfer to turbulence at specific phases. In all simulated cases, the positive buoyancy flux during episodes of convective instability at the zero-velocity phase is the dominant contributor to turbulent kinetic energy (TKE) although the shear production becomes increasingly important with increasing $$\unicode[STIX]{x1D6FD}$$ . Energy that initially resides wholly in mean available potential energy is lost through conversion to turbulence and the subsequent dissipation of TKE and turbulent available potential energy. A key result is that, in contrast to the explosive loss of energy during the initial convective instability in the non-sloping case, the sloping cases exhibit a more gradual energy loss that is sustained over a long time interval. The slope-parallel oscillation introduces a new flow time scale $$T=2\unicode[STIX]{x03C0}/(N\sin \unicode[STIX]{x1D6FD})$$ and, consequently, the fraction of initial APE that is converted to turbulence during convective instability progressively decreases with increasing $$\unicode[STIX]{x1D6FD}$$ . For moderate slopes with $$\unicode[STIX]{x1D6FD}<10^{\circ }$$ , most of the net energy loss takes place during an initial, short ( $$Nt\approx 20$$ ) interval with periodic convective overturns. For steeper slopes, most of the energy loss takes place during a later, long ( $Nt>100$ ) interval when both shear and convective instability occur, and the energy loss rate is approximately constant. The mixing efficiency during the initial period dominated by convectively driven turbulence is found to be substantially higher (exceeds 0.5) than the widely used value of 0.2. The mixing efficiency at long time in the present problem of a convective overturn at a boundary varies between 0.24 and 0.3.more » « less
-
Abstract A three-dimensional inertial model that conserves quasigeostrophic potential vorticity is proposed for wind-driven coastal upwelling along western boundaries. The dominant response to upwelling favorable winds is a surface-intensified baroclinic meridional boundary current with a subsurface countercurrent. The width of the current is not the baroclinic deformation radius but instead scales with the inertial boundary layer thickness while the depth scales as the ratio of the inertial boundary layer thickness to the baroclinic deformation radius. Thus, the boundary current scales depend on the stratification, wind stress, Coriolis parameter, and its meridional variation. In contrast to two-dimensional wind-driven coastal upwelling, the source waters that feed the Ekman upwelling are provided over the depth scale of this baroclinic current through a combination of onshore barotropic flow and from alongshore in the narrow boundary current. Topography forces an additional current whose characteristics depend on the topographic slope and width. For topography wider than the inertial boundary layer thickness the current is bottom intensified, while for narrow topography the current is wave-like in the vertical and trapped over the topography within the inertial boundary layer. An idealized primitive equation numerical model produces a similar baroclinic boundary current whose vertical length scale agrees with the theoretical scaling for both upwelling and downwelling favorable winds.more » « less
