Mixing-Driven Mean Flows and Submesoscale Eddies over Mid-Ocean Ridge Flanks and Fracture Zone Canyons

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

Bottom Boundary Potential Vorticity Injection from an Oscillating Flow: A PV Pump

Oceanic boundary currents over the continental slope exhibit variability with a range of time scales. Numerical studies of steady, along-slope currents over a sloping bathymetry have shown that cross-slope Ekman transport can advect buoyancy surfaces in a bottom boundary layer (BBL) so as to produce vertically sheared geostrophic flows that bring the total flow to rest: a process known as buoyancy shutdown of Ekman transport or Ekman arrest. This study considers the generation and evolution of near-bottom flows due to a barotropic, oscillating, and laterally sheared flow over a slope. The sensitivity of the boundary circulation to changes in oscillation frequency ω, background flow amplitude, bottom slope, and background stratification is explored. When ω/f ≪ 1, where f is the Coriolis frequency, oscillations allow the system to escape from the steady buoyancy shutdown scenario. The BBL is responsible for generating a secondary overturning circulation that produces vertical velocities that, combined with the potential vorticity (PV) anomalies of the imposed barotropic flow, give rise to a time-mean, rectified, vertical eddy PV flux into the ocean interior: a “PV pump.” In these idealized simulations, the PV anomalies in the BBL make a secondary contribution to the time-averaged PV flux. Numerical results show the domain-averaged eddy PV flux increases nonlinearly with ω with a peak near the inertial frequency, followed by a sharp decay for ω/f > 1. Different physical mechanisms are discussed that could give rise to the temporal variability of boundary currents

A simple physics-based improvement to the positive degree day model

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

Oceanic Bottom Boundary Layers and Abyssal Overturning Circulation

The vast amount of carbon and heat exchange between the abyssal and upper ocean and subsequently the atmosphere is paced by the abyssal overturning circulation. A key component of the abyssal overturning circulation is the formation and consumption of the densest water mass on Earth, Antarctic Bottom Water (AABW), namely the conversion of North Atlantic Deep Water (NADW) to AABW and the consumption of AABW via small-scale diapycnal mixing. Yet, this pathway of AABW spanning thousands of kilometers has not been successfully reproduced in large-scale general circulation models (GCM). What is missing is essentially the understanding and resolution of small-scale physics involved in converting deep and bottom waters from one density class to another, the water mass transformation (WMT). In this thesis, we focus on small-scale (both in the horizontal and vertical directions) dynamics near the BBL, where enhanced shear, mixing and turbulence exist to facilitate effective WMT above the seafloor. From high-resolution ocean glider observations around Antarctica, we find that a portion of Lower Circumpolar Deep Water, a branch of NADW, becomes lighter via mixing with light shelf water over the continental slope and shelf, instead of being converted into dense AABW under sea ice. This mixing is likely induced by submesoscale symmetric instability coming from a strong boundary current interacting with the sloping topography in the BBL. We then consider how to sustain the consumption of AABW in the global mid-ocean ridge system. Using numerical models, we show that submesoscale baroclinic eddies are crucial to maintaining strong stratification over the flanks of the mid-ocean ridges and thus permitting effective WMT. Lastly, we consider the interaction between external mean flows and stratified BBL over sloping topography. With the large-scale turbulence resolved in a large-eddy simulation model, we propose a new theoretical framework to describe the evolution of the BBL and the Eulerian advection of its associated stratification when external barotropic flows are present. This new framework can be used to parameterize bottom friction, important for closing the kinetic energy budget of the global ocean. We further extend this interaction to a horizontally-sheared and temporally-oscillating external mean flow and explore the response of the BBL and the BBL-interior mass exchange with simple turbulent parameterizations. Using a combination of different approaches, we confirm that the long-overlooked oceanic BBL is the key location for closing the abyssal overturning circulation. More importantly, without appropriate techniques to tackle the currently unresolved small-scale processes, they will likely remain a narrow bottleneck in understanding the abyssal overturning circulation.</p

Bottom Boundary Potential Vorticity Injection from an Oscillating Flow: A PV Pump

Oceanic boundary currents over the continental slope exhibit variability with a range of time scales. Numerical studies of steady, along-slope currents over a sloping bathymetry have shown that cross-slope Ekman transport can advect buoyancy surfaces in a bottom boundary layer (BBL) so as to produce vertically sheared geostrophic flows that bring the total flow to rest: a process known as buoyancy shutdown of Ekman transport or Ekman arrest. This study considers the generation and evolution of near-bottom flows due to a barotropic, oscillating, and laterally sheared flow over a slope. The sensitivity of the boundary circulation to changes in oscillation frequency ω, background flow amplitude, bottom slope, and background stratification is explored. When ω/f ≪ 1, where f is the Coriolis frequency, oscillations allow the system to escape from the steady buoyancy shutdown scenario. The BBL is responsible for generating a secondary overturning circulation that produces vertical velocities that, combined with the potential vorticity (PV) anomalies of the imposed barotropic flow, give rise to a time-mean, rectified, vertical eddy PV flux into the ocean interior: a “PV pump.” In these idealized simulations, the PV anomalies in the BBL make a secondary contribution to the time-averaged PV flux. Numerical results show the domain-averaged eddy PV flux increases nonlinearly with ω with a peak near the inertial frequency, followed by a sharp decay for ω/f > 1. Different physical mechanisms are discussed that could give rise to the temporal variability of boundary currents

LMG 14-11: Cruise Report - ChinStrAP: Changes in Stratification at the Antarctic Peninsula

The goal of the ChinStrAP (Changes in Stratification at the Antarctic Peninsula) project is to assess the role of mesoscale and submesoscale variability on water mass transformation and exchange across the continental shelf and slope in southern Drake Passage. Specifically we plan to: 1. Assess the influence of flow-topography interactions on the frequency and characteristics of mesoscale eddy generation along the southern boundary of the ACC. We will be sampling a region where the strong southern ACC front (SACCF) and the ACC’s southern boundary run along the continental slope before interacting with the Shackleton Fracture Zone (SFZ) and deflecting northward. This is a known region of eddy generation, as observed from remotely-sensed sea surface height (SSH), sea surface temperature (SST) and ocean color observations. We hope to obtain the in situ observations necessary to determine the mechanisms by which these eddies are formed and how they contribute to cross-shelf exchange. 2. Explore the interactions between surface wind and buoyancy forcing on mixed layer depth variability and its implications for ventilation of the deep ocean. Southern Drake Passage is a key location where deep isopycnals rise towards the surface across the ACC and outcrop allowing direct exchange with atmospheric temperatures and gases. This process is critical to the equilibration of dissolved gas concentrations with atmospheric values and thus influences large-scale characteristics of Earth’s climate. A large number of recent studies have pointed to the strong influence of submesoscale processes, both due to surface forcing and stirring by mesoscale eddies, on rapid changes in mixed layer depth. These changes come about through dynamic instabilities related to lateral gradients in mixed layer properties. The ACC is a location where strong lateral fronts align with strong westerly winds. This situation is similar to western boundary currents, however the internal density structure is considerably different in the ACC. Our measurement strategy should allow us to capture the evolution of these dynamical processes. 3. Carry out an XBT/XCTD transect across Drake Passage on the southbound leg as a contribution to Scripps High Resolution XBT/XCTD observing line (WOCE AX22). This information will provide larger-scale context for the data collected from the three gliders. 4. Determine with high spatially- and temporally-resolved measurements the characteristic internal variability in this regional current system. This work will help to better interpret long-term observations in this same location (e.g., LTER monitoring, the AX22 high resolution XBT/XCTD line). The observations will also help validate high frequency variability in numerical models as they push towards resolving key dynamical processes at the continental shelf break

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

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 (&lt; 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.First author draf

The Evolution and Arrest of a Turbulent Stratified Oceanic Bottom Boundary Layer over a Slope: Downslope Regime

The dynamics of a stratified oceanic bottom boundary layer (BBL) over an insulating, sloping surface depend critically on the intersection of density surfaces with the bottom. For an imposed along-slope flow, the cross-slope Ekman transport advects density surfaces and generates a near-bottom geostrophic thermal wind shear that opposes the background flow. A limiting case occurs when a momentum balance is achieved between the Coriolis force and a restoring buoyancy force in response to the displacement of stratified fluid over the slope: this is known as Ekman arrest. However, the turbulent characteristics that accompany this adjustment have received less attention. We present two estimates to characterize the state of the BBL based on the mixed layer thickness: H_a and H_L. The former characterizes the steady Ekman arrested state, and the latter characterizes a relaminarized state. The derivation of H_L makes use of a newly defined slope Obukhov length L_s that characterizes the relative importance of shear production and cross-slope buoyancy advection. The value of H_a can be combined with the temporally evolving depth of the mixed layer H to form a nondimensional variable H/H_a that provides a similarity prediction of the BBL evolution across different turbulent regimes. The length scale L_s can also be used to obtain an expression for the wall stress when the BBL relaminarizes. We validate these relationships using output from a suite of three-dimensional large-eddy simulations. We conclude that the BBL reaches the relaminarized state before the steady Ekman arrested state. Calculating H/H_a and H/H_L from measurements will provide information on the stage of oceanic BBL development being observed. These diagnostics may also help to improve numerical parameterizations of stratified BBL dynamics over sloping topography