Search

Slowly evolving stratified flow over rough topography is subject to substantial drag due to internal motions, but often numerical simulations are carried out at resolutions where this “wave” drag must be parameterized. Here we highlight the importance of internal drag from topography with scales that cannot radiate internal waves, but may be highly nonlinear, and we propose a simple parameterization of this drag that has a minimum of fit parameters compared to existing schemes. The parameterization smoothly transitions from a quadratic drag law () for low Nh/u0 (linear wave dynamics) to a linear drag law () for high Nh/u0 flows (nonlinear blocking and hydraulic dynamics), where N is the stratification, h is the height of the topography, and u0 is the near-bottom velocity; the parameterization does not have a dependence on Coriolis frequency. Simulations carried out in a channel with synthetic bathymetry and steady body forcing indicate that this parameterization accurately predicts drag across a broad range of forcing parameters when the effect of reduced near-bottom mixing is taken into account by reducing the effective height of the topography. The parameterization is also tested in simulations of wind-driven channel flows that generate mesoscale eddy fields, a setup where the downstream transport is sensitive to the bottom drag parameterization and its effect on the eddies. In these simulations, the parameterization replicates the effect of rough bathymetry on the eddies. If extrapolated globally, the subinertial topographic scales can account for 2.7 TW of work done on the low-frequency circulation, an important sink that is redistributed to mixing in the open ocean.

The large-scale middepth circulation and eddy diffusivities in the southeast Pacific Ocean and Scotia Sea sectors between 110° and 45°W of the Antarctic Circumpolar Current (ACC) are described based on a subsurface quasi-isobaric RAFOS-float-based Lagrangian dataset. These RAFOS float data were collected during the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES). The mean flow, adjusted to a common 1400-m depth, shows the presence of jets in the time-averaged sense with speeds of 6 cm s−1 in the southeast Pacific Ocean and upward of 13 cm s−1 in the Scotia Sea. These jets appear to be locked to topography in the Scotia Sea but, aside from negotiating a seamount chain, are mostly free of local topographic constraints in the southeast Pacific Ocean. The eddy kinetic energy (EKE) is higher than the mean kinetic energy everywhere in the sampled domain by about 50%. The magnitude of the EKE increases drastically (by a factor of 2 or more) as the current crosses over the Hero and Shackleton fracture zones into the Scotia Sea. The meridional isopycnal stirring shows lateral and vertical variations with local eddy diffusivities as high as 2800 ± 600 m2 s−1 at 700 m decreasing to 990 ± 200 m2 s−1 at 1800 m in the southeast Pacific Ocean. However, the cross-ACC diffusivity in the southeast Pacific Ocean is significantly lower, with values of 690 ± 150 and 1000 ± 200 m2 s−1 at shallow and deep levels, respectively, due to the action of jets. The cross-ACC diffusivity in the Scotia Sea is about 1200 ± 500 m2 s−1.

As part of the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES), 210 subsurface floats were deployed west of the Drake Passage on two targeted density surfaces. Absolute (single particle) diffusivities are calculated for the floats. The focus is on the meridional component, which is less affected by the mean shear. The diffusivities are estimated in several ways, including a novel method based on the probability density function of the meridional displacements. This allows the determination of the range of possible lateral diffusivities, as well as the period over which the spreading can be said to be diffusive. The method is applied to the float data and to synthetic trajectories generated with the Massachusetts Institute of Technology General Circulation Model (MITgcm). Because of ballasting problems, many of the floats did not remain on their targeted density surface. However, the float temperature records suggest that most occupied a small range of densities, so the floats were grouped together for the analysis. The latter focuses on a subset of 109 of the floats, launched near 105°W. The different methods yield a consistent estimate for the diffusivity of 800 ± 200 m2 s−1. The same calculations were made with model particles deployed on 20 different density surfaces and the result for the particles deployed on the neutral density surface γ = 27.7 surface was the same within the errors. The model was then used to map the variation of the diffusivity in the vertical, near the core of the Antarctic Circumpolar Current (ACC). The results suggest mixing is intensified at middepths, between 1500 and 2000 m, consistent with several previous studies.