# Browse

## Abstract

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 (*Nh*/*u*
_{0} (linear wave dynamics) to a linear drag law (*Nh*/*u*
_{0} flows (nonlinear blocking and hydraulic dynamics), where *N* is the stratification, *h* is the height of the topography, and *u*
_{0} 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.

## Abstract

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 (*Nh*/*u*
_{0} (linear wave dynamics) to a linear drag law (*Nh*/*u*
_{0} flows (nonlinear blocking and hydraulic dynamics), where *N* is the stratification, *h* is the height of the topography, and *u*
_{0} 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.

## Abstract

As part of the Flow Encountering Abrupt Topography (FLEAT) program, an array of pressure-sensor equipped inverted echo sounders (PIESs) was deployed north of Palau where the westward-flowing North Equatorial Current encounters the southern end of the Kyushu–Palau Ridge in the tropical North Pacific. Capitalizing on concurrent observations from satellite altimetry, FLEAT Spray gliders, and shipboard hydrography, the PIESs’ 10-month duration hourly bottom pressure *p* and round-trip acoustic travel time *τ* records are used to examine the magnitude and predictability of sea level and pycnocline depth changes and to track signal propagations through the array. Sea level and pycnocline depth are found to vary in response to a range of ocean processes, with their magnitude and predictability strongly process dependent. Signals characterized here comprise the barotropic tides, semidiurnal and diurnal internal tides, southeastward-propagating superinertial waves, westward-propagating mesoscale eddies, and a strong signature of sea level increase and pycnocline deepening associated with the region’s relaxation from El Niño to La Niña conditions. The presence of a broad band of superinertial waves just above the inertial frequency was unexpected and the FLEAT observations and output from a numerical model suggest that these waves detected near Palau are forced by remote winds east of the Philippines. The PIES-based estimates of pycnocline displacement are found to have large uncertainties relative to overall variability in pycnocline depth, as localized deep current variations arising from interactions of the large-scale currents with the abrupt topography around Palau have significant travel time variability.

## Abstract

As part of the Flow Encountering Abrupt Topography (FLEAT) program, an array of pressure-sensor equipped inverted echo sounders (PIESs) was deployed north of Palau where the westward-flowing North Equatorial Current encounters the southern end of the Kyushu–Palau Ridge in the tropical North Pacific. Capitalizing on concurrent observations from satellite altimetry, FLEAT Spray gliders, and shipboard hydrography, the PIESs’ 10-month duration hourly bottom pressure *p* and round-trip acoustic travel time *τ* records are used to examine the magnitude and predictability of sea level and pycnocline depth changes and to track signal propagations through the array. Sea level and pycnocline depth are found to vary in response to a range of ocean processes, with their magnitude and predictability strongly process dependent. Signals characterized here comprise the barotropic tides, semidiurnal and diurnal internal tides, southeastward-propagating superinertial waves, westward-propagating mesoscale eddies, and a strong signature of sea level increase and pycnocline deepening associated with the region’s relaxation from El Niño to La Niña conditions. The presence of a broad band of superinertial waves just above the inertial frequency was unexpected and the FLEAT observations and output from a numerical model suggest that these waves detected near Palau are forced by remote winds east of the Philippines. The PIES-based estimates of pycnocline displacement are found to have large uncertainties relative to overall variability in pycnocline depth, as localized deep current variations arising from interactions of the large-scale currents with the abrupt topography around Palau have significant travel time variability.

## Abstract

Hydrographic measurements recently acquired along the thalweg of the Lifamatola Passage combined with historical moored velocity measurements immediately downstream of the sill are used to study the hydraulics, transport, mixing, and entrainment in the dense overflow. The observations suggest that the mean overflow is nearly critical at the mooring site, suggesting that a weir formula may be appropriate for estimating the overflow transport. Our assessment suggests that the weir formulas corresponding to a rectangular, triangular, or parabolic cross section all result in transports very close to the observation, suggesting their potential usage in long-term monitoring of the overflow transport or parameterizing the transport in numerical models. Analyses also suggest that deep signals within the overflow layer are blocked by the shear flow from propagating upstream, whereas the shallow wave modes of the full-depth continuously stratified flow are able to propagate upstream from the Banda Sea into the Maluku Sea. Strong mixing is found immediately downstream of the sill crest, with Thorpe-scale-based estimates of the mean dissipation rate within the overflow up to 1.1 × 10^{−7} W kg^{−1} and the region-averaged diapycnal diffusivity within the downstream overflow in the range of 2.3 × 10^{−3} to 10.1 × 10^{−3} m^{2} s^{−1}. Mixing in the Lifamatola Passage results in 0.6–1.2-Sv (1 Sv ≡ 10^{6} m^{3} s^{−1}) entrainment transport added to the overflow, enhancing the deep-water renewal in the Banda Sea. A bulk diffusivity coefficient estimated in the deep Banda Sea yields 1.6 × 10^{−3} ± 5 × 10^{−4} m^{2} s^{−1}, with an associated downward turbulent heat flux of 9 W m^{−2}.

## Abstract

Hydrographic measurements recently acquired along the thalweg of the Lifamatola Passage combined with historical moored velocity measurements immediately downstream of the sill are used to study the hydraulics, transport, mixing, and entrainment in the dense overflow. The observations suggest that the mean overflow is nearly critical at the mooring site, suggesting that a weir formula may be appropriate for estimating the overflow transport. Our assessment suggests that the weir formulas corresponding to a rectangular, triangular, or parabolic cross section all result in transports very close to the observation, suggesting their potential usage in long-term monitoring of the overflow transport or parameterizing the transport in numerical models. Analyses also suggest that deep signals within the overflow layer are blocked by the shear flow from propagating upstream, whereas the shallow wave modes of the full-depth continuously stratified flow are able to propagate upstream from the Banda Sea into the Maluku Sea. Strong mixing is found immediately downstream of the sill crest, with Thorpe-scale-based estimates of the mean dissipation rate within the overflow up to 1.1 × 10^{−7} W kg^{−1} and the region-averaged diapycnal diffusivity within the downstream overflow in the range of 2.3 × 10^{−3} to 10.1 × 10^{−3} m^{2} s^{−1}. Mixing in the Lifamatola Passage results in 0.6–1.2-Sv (1 Sv ≡ 10^{6} m^{3} s^{−1}) entrainment transport added to the overflow, enhancing the deep-water renewal in the Banda Sea. A bulk diffusivity coefficient estimated in the deep Banda Sea yields 1.6 × 10^{−3} ± 5 × 10^{−4} m^{2} s^{−1}, with an associated downward turbulent heat flux of 9 W m^{−2}.

## Abstract

Ocean lee waves occur on length scales that are smaller than the grid scale of global circulation models (GCMs). Therefore, such models must parameterize the drag associated with launching lee waves. This paper compares the lee wave drag predicted by existing parameterizations with the drag measured in high-resolution nonhydrostatic numerical simulations of a lee wave over periodic sinusoidal bathymetry. The simulations afford a time-varying glimpse at the nonlinear and nonhydrostatic oceanic lee wave spinup process and identify a characteristic time scale to reach steady state. The maximum instantaneous lee wave drag observed during the spinup period is found to be well predicted by linear lee wave theory for all hill heights. In steady state, the simulations demonstrate the applicability of parameterizing the drag based on applying linear theory to the lowest overtopping streamline of the flow (LOTS), as is currently employed in GCMs. However, because existing parameterizations are based only on the height of the LOTS, they implicitly assume hydrostatic flow. For hills tall enough to trap water in their valleys, the simulations identify a set of nonhydrostatic processes that can result in a reduction of the lee wave drag from that given by hydrostatic parameterizations. The simulations suggest implementing a time-dependent nonhydrostatic version of the LOTS-based parameterization of lee wave drag and demonstrate the remarkable applicability of linear lee wave theory to oceanic lee waves.

## Abstract

Ocean lee waves occur on length scales that are smaller than the grid scale of global circulation models (GCMs). Therefore, such models must parameterize the drag associated with launching lee waves. This paper compares the lee wave drag predicted by existing parameterizations with the drag measured in high-resolution nonhydrostatic numerical simulations of a lee wave over periodic sinusoidal bathymetry. The simulations afford a time-varying glimpse at the nonlinear and nonhydrostatic oceanic lee wave spinup process and identify a characteristic time scale to reach steady state. The maximum instantaneous lee wave drag observed during the spinup period is found to be well predicted by linear lee wave theory for all hill heights. In steady state, the simulations demonstrate the applicability of parameterizing the drag based on applying linear theory to the lowest overtopping streamline of the flow (LOTS), as is currently employed in GCMs. However, because existing parameterizations are based only on the height of the LOTS, they implicitly assume hydrostatic flow. For hills tall enough to trap water in their valleys, the simulations identify a set of nonhydrostatic processes that can result in a reduction of the lee wave drag from that given by hydrostatic parameterizations. The simulations suggest implementing a time-dependent nonhydrostatic version of the LOTS-based parameterization of lee wave drag and demonstrate the remarkable applicability of linear lee wave theory to oceanic lee waves.

## Abstract

Turbulent mixing adjacent to the Velasco Reef and Kyushu–Palau Ridge, off northern Palau in the western equatorial Pacific Ocean, is examined using shipboard and moored observations. The study focuses on a 9-day-long, ship-based microstructure and velocity survey, conducted in November–December 2016. Several sections (9–15 km in length) of microstructure, hydrographic, and velocity fields were acquired over and around the reef, where water depths ranged from 50 to 3000 m. Microstructure profiles were collected while steaming slowly either toward or away from the reef, and underway current surveys were conducted along quasi-rectangular boxes with side lengths of 5–10 km. Near the reef, both tidal and subtidal motions were important, while subtidal motions were stronger away from the reef. Vertical shears of currents and mixing were stronger on the northern and eastern flanks of the reef than on the western flanks. High turbulent kinetic energy dissipation rates, 10^{−6}–10^{−4} W kg^{−1}, and large values of eddy diffusivities, 10^{−4}–10^{−2} m^{2} s^{−1}, with strong turbulent heat fluxes, 100–500 W m^{−2}, were found. Currents flowing along the eastern side separated at the northern tip of the reef and generated submesoscale cyclonic vorticity of about 2–4 times the planetary vorticity. The analysis suggests that a torque, imparted by the turbulent bottom stress, generated the cyclonic vorticity at the northern boundary. The northern reef is associated with high vertical transports resulting from both submesoscale flow convergences and energetic mixing. Even though the area around Palau represents a small footprint of the ocean, vertical velocities and mixing rates are several orders magnitude larger than in the open ocean.

## Abstract

Turbulent mixing adjacent to the Velasco Reef and Kyushu–Palau Ridge, off northern Palau in the western equatorial Pacific Ocean, is examined using shipboard and moored observations. The study focuses on a 9-day-long, ship-based microstructure and velocity survey, conducted in November–December 2016. Several sections (9–15 km in length) of microstructure, hydrographic, and velocity fields were acquired over and around the reef, where water depths ranged from 50 to 3000 m. Microstructure profiles were collected while steaming slowly either toward or away from the reef, and underway current surveys were conducted along quasi-rectangular boxes with side lengths of 5–10 km. Near the reef, both tidal and subtidal motions were important, while subtidal motions were stronger away from the reef. Vertical shears of currents and mixing were stronger on the northern and eastern flanks of the reef than on the western flanks. High turbulent kinetic energy dissipation rates, 10^{−6}–10^{−4} W kg^{−1}, and large values of eddy diffusivities, 10^{−4}–10^{−2} m^{2} s^{−1}, with strong turbulent heat fluxes, 100–500 W m^{−2}, were found. Currents flowing along the eastern side separated at the northern tip of the reef and generated submesoscale cyclonic vorticity of about 2–4 times the planetary vorticity. The analysis suggests that a torque, imparted by the turbulent bottom stress, generated the cyclonic vorticity at the northern boundary. The northern reef is associated with high vertical transports resulting from both submesoscale flow convergences and energetic mixing. Even though the area around Palau represents a small footprint of the ocean, vertical velocities and mixing rates are several orders magnitude larger than in the open ocean.

## Abstract

Towed shipboard and moored observations show internal gravity waves over a tall, supercritical submarine ridge that reaches to 1000 m below the ocean surface in the tropical western Pacific north of Palau. The lee-wave or topographic Froude number, *Nh*
_{0}/*U*
_{0} (where *N* is the buoyancy frequency, *h*
_{0} the ridge height, and *U*
_{0} the farfield velocity), ranged between 25 and 140. The waves were generated by a superposition of tidal and low-frequency flows and thus had two distinct energy sources with combined amplitudes of up to 0.2 m s^{−1}. Local breaking of the waves led to enhanced rates of dissipation of turbulent kinetic energy reaching above 10^{−6} W kg^{−1} in the lee of the ridge near topography. Turbulence observations showed a stark contrast between conditions at spring and neap tide. During spring tide, when the tidal flow dominated, turbulence was approximately equally distributed around both sides of the ridge. During neap tide, when the mean flow dominated over tidal oscillations, turbulence was mostly observed on the downstream side of the ridge relative to the mean flow. The drag exerted by the ridge on the flow, estimated to

## Abstract

Towed shipboard and moored observations show internal gravity waves over a tall, supercritical submarine ridge that reaches to 1000 m below the ocean surface in the tropical western Pacific north of Palau. The lee-wave or topographic Froude number, *Nh*
_{0}/*U*
_{0} (where *N* is the buoyancy frequency, *h*
_{0} the ridge height, and *U*
_{0} the farfield velocity), ranged between 25 and 140. The waves were generated by a superposition of tidal and low-frequency flows and thus had two distinct energy sources with combined amplitudes of up to 0.2 m s^{−1}. Local breaking of the waves led to enhanced rates of dissipation of turbulent kinetic energy reaching above 10^{−6} W kg^{−1} in the lee of the ridge near topography. Turbulence observations showed a stark contrast between conditions at spring and neap tide. During spring tide, when the tidal flow dominated, turbulence was approximately equally distributed around both sides of the ridge. During neap tide, when the mean flow dominated over tidal oscillations, turbulence was mostly observed on the downstream side of the ridge relative to the mean flow. The drag exerted by the ridge on the flow, estimated to

## Abstract

The effects of topography on the linear stability of both barotropic vortices and two-layer, baroclinic vortices are examined by considering cylindrical topography and vortices with stepwise relative vorticity profiles in the quasigeostrophic approximation. Four vortex configurations are considered, classified by the number of relative vorticity steps in the horizontal and the number of layers in the vertical: barotropic one-step vortex (Rankine vortex), barotropic two-step vortex, and their two-layer, baroclinic counterparts with the vorticity steps in the upper layer. In the barotropic calculation, the vortex is destabilized by topography having an oppositely signed potential vorticity jump while stabilized by topography of same-signed jump, that is, anticyclones are destabilized by seamounts while stabilized by depressions. Further, topography of appropriate sign and magnitude can excite a mode-1 instability for a two-step vortex, especially relevant for topographic encounters of an otherwise stable vortex. The baroclinic calculation is in general consistent with the barotropic calculation except that the growth rate weakens and, for a two-step vortex, becomes less sensitive to topography (sign and magnitude) as baroclinicity increases. The smaller growth rate for a baroclinic vortex is consistent with previous findings that vortices with sufficient baroclinic structure could cross the topography relatively easily. Nonlinear contour dynamics simulations are conducted to confirm the linear stability analysis and to describe the subsequent evolution.

## Abstract

The effects of topography on the linear stability of both barotropic vortices and two-layer, baroclinic vortices are examined by considering cylindrical topography and vortices with stepwise relative vorticity profiles in the quasigeostrophic approximation. Four vortex configurations are considered, classified by the number of relative vorticity steps in the horizontal and the number of layers in the vertical: barotropic one-step vortex (Rankine vortex), barotropic two-step vortex, and their two-layer, baroclinic counterparts with the vorticity steps in the upper layer. In the barotropic calculation, the vortex is destabilized by topography having an oppositely signed potential vorticity jump while stabilized by topography of same-signed jump, that is, anticyclones are destabilized by seamounts while stabilized by depressions. Further, topography of appropriate sign and magnitude can excite a mode-1 instability for a two-step vortex, especially relevant for topographic encounters of an otherwise stable vortex. The baroclinic calculation is in general consistent with the barotropic calculation except that the growth rate weakens and, for a two-step vortex, becomes less sensitive to topography (sign and magnitude) as baroclinicity increases. The smaller growth rate for a baroclinic vortex is consistent with previous findings that vortices with sufficient baroclinic structure could cross the topography relatively easily. Nonlinear contour dynamics simulations are conducted to confirm the linear stability analysis and to describe the subsequent evolution.

## Abstract

Microstructure measurements in Drake Passage and on the flanks of Kerguelen Plateau find turbulent dissipation rates *ε* on average factors of 2–3 smaller than linear lee-wave generation predictions, as well as a factor of 3 smaller than the predictions of a well-established parameterization based on finescale shear and strain. Here, the possibility that these discrepancies are a result of conservation of wave action *E*/*ω*
_{L} = *E*/|*kU*| is explored. Conservation of wave action will transfer a fraction of the lee-wave radiation back to the mean flow if the waves encounter weakening currents *U*, where the intrinsic or Lagrangian frequency *ω*
_{L} = |*kU*| ↓ |*f*| and *k* the along-stream horizontal wavenumber, where *kU* ≡ **k** ⋅ **V**. The dissipative fraction of power that is lost to turbulence depends on the Doppler shift of the intrinsic frequency between generation and breaking, hence on the topographic height spectrum and bandwidth *N*/*f*. The partition between dissipation and loss to the mean flow is quantified for typical topographic height spectral shapes and *N*/*f* ratios found in the abyssal ocean under the assumption that blocking is local in wavenumber. Although some fraction of lee-wave generation is always dissipated in a rotating fluid, lee waves are not as large a sink for balanced energy or as large a source for turbulence as previously suggested. The dissipative fraction is 0.44–0.56 for topographic spectral slopes and buoyancy frequencies typical of the deep Southern Ocean, insensitive to flow speed *U* and topographic splitting. Lee waves are also an important mechanism for redistributing balanced energy within their generating bottom current.

## Abstract

Microstructure measurements in Drake Passage and on the flanks of Kerguelen Plateau find turbulent dissipation rates *ε* on average factors of 2–3 smaller than linear lee-wave generation predictions, as well as a factor of 3 smaller than the predictions of a well-established parameterization based on finescale shear and strain. Here, the possibility that these discrepancies are a result of conservation of wave action *E*/*ω*
_{L} = *E*/|*kU*| is explored. Conservation of wave action will transfer a fraction of the lee-wave radiation back to the mean flow if the waves encounter weakening currents *U*, where the intrinsic or Lagrangian frequency *ω*
_{L} = |*kU*| ↓ |*f*| and *k* the along-stream horizontal wavenumber, where *kU* ≡ **k** ⋅ **V**. The dissipative fraction of power that is lost to turbulence depends on the Doppler shift of the intrinsic frequency between generation and breaking, hence on the topographic height spectrum and bandwidth *N*/*f*. The partition between dissipation and loss to the mean flow is quantified for typical topographic height spectral shapes and *N*/*f* ratios found in the abyssal ocean under the assumption that blocking is local in wavenumber. Although some fraction of lee-wave generation is always dissipated in a rotating fluid, lee waves are not as large a sink for balanced energy or as large a source for turbulence as previously suggested. The dissipative fraction is 0.44–0.56 for topographic spectral slopes and buoyancy frequencies typical of the deep Southern Ocean, insensitive to flow speed *U* and topographic splitting. Lee waves are also an important mechanism for redistributing balanced energy within their generating bottom current.

## Abstract

In this study, a 2-yr time series of velocity profiles to 1000 m from meridional glider surveys is used to characterize the wake in the lee of a large island in the western tropical North Pacific Ocean, Palau. Surveys were completed along sections to the east and west of the island to capture both upstream and downstream conditions. Objectively mapped in time and space, mean sections of velocity show the incident westward North Equatorial Current accelerating around the island of Palau, increasing from 0.1 to 0.2 m s^{−1} at the surface. Downstream of the island, elevated velocity variability and return flow in the lee are indicative of boundary layer separation. Isolating for periods of depth-average westward flow reveals a length scale in the wake that reflects local details of the topography. Eastward flow is shown to produce an asymmetric wake. Depth-average velocity time series indicate that energetic events (on time scales from weeks to months) are prevalent. These events are associated with mean vorticity values in the wake up to 0.3*f* near the surface and with instantaneous values that can exceed *f* (the local Coriolis frequency) during periods of sustained, anomalously strong westward flow. Thus, ageostrophic effects become important to first order.

## Abstract

In this study, a 2-yr time series of velocity profiles to 1000 m from meridional glider surveys is used to characterize the wake in the lee of a large island in the western tropical North Pacific Ocean, Palau. Surveys were completed along sections to the east and west of the island to capture both upstream and downstream conditions. Objectively mapped in time and space, mean sections of velocity show the incident westward North Equatorial Current accelerating around the island of Palau, increasing from 0.1 to 0.2 m s^{−1} at the surface. Downstream of the island, elevated velocity variability and return flow in the lee are indicative of boundary layer separation. Isolating for periods of depth-average westward flow reveals a length scale in the wake that reflects local details of the topography. Eastward flow is shown to produce an asymmetric wake. Depth-average velocity time series indicate that energetic events (on time scales from weeks to months) are prevalent. These events are associated with mean vorticity values in the wake up to 0.3*f* near the surface and with instantaneous values that can exceed *f* (the local Coriolis frequency) during periods of sustained, anomalously strong westward flow. Thus, ageostrophic effects become important to first order.

## Abstract

Drag and turbulence in steady stratified flows over “abyssal hills” have been parameterized using linear theory and rates of energy cascade due to wave–wave interactions. Linear theory has no drag or energy loss due to large-scale bathymetry because waves with intrinsic frequency less than the Coriolis frequency are evanescent. Numerical work has tested the theory by high passing the topography and estimating the radiation and turbulence. Adding larger-scale bathymetry that would generate evanescent internal waves generates nonlinear and turbulent flow, driving a dissipation approximately twice that of the radiating waves for the topographic spectrum chosen. This drag is linear in the forcing velocity, in contrast to atmospheric parameterizations that have quadratic drag. Simulations containing both small- and large-scale bathymetry have more dissipation than just adding the large- and small-scale dissipations together, so the scales couple. The large-scale turbulence is localized, generally in the lee of large obstacles. Medium-scale regional models partially resolve the “nonpropagating” wavenumbers, leading to the question of whether they need the large-scale energy loss to be parameterized. Varying the resolution of the simulations indicates that if the ratio of gridcell height to width is less than the root-mean-square topographic slope, then the dissipation is overestimated in coarse models (by up to 25%); conversely, it can be underestimated by up to a factor of 2 if the ratio is greater. Most regional simulations are likely in the second regime and should have extra drag added to represent the large-scale bathymetry, and the deficit is at least as large as that parameterized for abyssal hills.

## Abstract

Drag and turbulence in steady stratified flows over “abyssal hills” have been parameterized using linear theory and rates of energy cascade due to wave–wave interactions. Linear theory has no drag or energy loss due to large-scale bathymetry because waves with intrinsic frequency less than the Coriolis frequency are evanescent. Numerical work has tested the theory by high passing the topography and estimating the radiation and turbulence. Adding larger-scale bathymetry that would generate evanescent internal waves generates nonlinear and turbulent flow, driving a dissipation approximately twice that of the radiating waves for the topographic spectrum chosen. This drag is linear in the forcing velocity, in contrast to atmospheric parameterizations that have quadratic drag. Simulations containing both small- and large-scale bathymetry have more dissipation than just adding the large- and small-scale dissipations together, so the scales couple. The large-scale turbulence is localized, generally in the lee of large obstacles. Medium-scale regional models partially resolve the “nonpropagating” wavenumbers, leading to the question of whether they need the large-scale energy loss to be parameterized. Varying the resolution of the simulations indicates that if the ratio of gridcell height to width is less than the root-mean-square topographic slope, then the dissipation is overestimated in coarse models (by up to 25%); conversely, it can be underestimated by up to a factor of 2 if the ratio is greater. Most regional simulations are likely in the second regime and should have extra drag added to represent the large-scale bathymetry, and the deficit is at least as large as that parameterized for abyssal hills.