The relationship between boundary currents generated by tidal flow over topography and the radiated internal wave power is examined in two-dimensional numerical simulations of a uniformly stratified fluid. The radiated power PIW and kinetic energy density of the boundary currents are computed as a function of the internal wave slope SIW and the criticality parameter ε (ratio of the maximum topographic slope to SIW). Both SIW and ε are varied two orders of magnitude about unity by changing the tidal frequency, stratification, or topographic shape and slope. We consider cases where the hydrostatic approximation is valid (SIW ≪ 1), as well as test theoretical predictions for models of the deep ocean where the beam slope diverges and the hydrostatic approximation fails. We confirm that resonant boundary currents characterized by large kinetic energy densities form over critical topography (ε = 1). However, we find that this resonance phenomenon does not extend to the power radiated by internal waves that propagate away from the topography. Further, by directly comparing the kinetic energy density to the energy flux of the generated internal waves, we find that the more easily measured kinetic energy density cannot be used as a proxy to characterize the conversion of tidal energy to radiated internal wave power. Whether the hydrostatic approximation is valid or fails, our measurements of the radiated power can be described as PIW = Ptidef (ε, shape)/SIW, where Ptide is the effective tidal power that interacts with the topography, and π/8 < f (ε, shape) < π/4 is bounded below by the theoretical prediction of Bell [“Topographically generated internal waves in the open ocean,” J. Geophys. Res. 80, 320–327 (1975)] for ε → 0 and above by Llewellyn Smith and Young [“Tidal conversion at a very steep ridge,” J. Fluid Mech. 495, 175–191 (2003)] for ε → ∞.
I. INTRODUCTION
Internal waves play an important role in the global ocean energy budget by providing mixing necessary to maintain the thermohaline circulation.1 Approximately half of the internal wave energy in the ocean is produced by tidal flow over topography, resulting in waves at the tidal frequency, known as internal tides. Conversion of the barotropic tide into internal tides occurs at continental slopes,2–4 sea-mounts,5,6 ridges,7,8 and rough topography in the deep ocean.9,10 An understanding of the efficiency of this conversion is needed to understand the energy budget of the oceans. Therefore, it is useful to measure the total power PIW converted from barotropic tidal motions over topography into radiated internal waves,
where
We are interested in the conversion of tidal energy into internal wave power, which we characterize as the ratio of the radiated internal wave power PIW to the effective tidal power Ptide,
In addition to radiating internal waves, barotropic tidal motions over bottom topography produce boundary layer currents. These currents can vary with the topographic steepness, which is characterized by the criticality parameter
where Stopo is the maximum topographic slope and
A number of past studies (see Sec. II) have examined the conversion of barotropic tidal motions to internal tides by two-dimensional topography in the hydrostatic approximation, which assumes that horizontal length scales far exceed vertical length scales. This approximation is appropriate for the strongly stratified shallow ocean, where the local buoyancy frequency N(z) is much greater than the tidal frequency ω. However, King et al.11 found that in the deep ocean the buoyancy frequency can become comparable to or even less than the semi-diurnal lunar tidal frequency ωM2 = 1.4052 × 10−4 rad s−1 for depths greater than about 4 km. As an example, measurements from the Venezuelan Basin (see Fig. 1) show that below about 4 km depth the local buoyancy frequency is less than ωM2.
The local buoyancy frequency N(z) (black circles, black axes on the left and bottom) decreases by more than two orders of magnitude from the shallow ocean down to 5 km depth in the Venezuelan Basin in the Caribbean. The data are taken from cast A22-316N151-4-10 from the World Ocean Circulation Experiment, latitude 13.17° N, longitude 66.00° W, and analyzed (averaging over 50 m depths) by the method described by King et al.11 For depths greater than 4 km, the local buoyancy frequency becomes less than the semi-diurnal lunar tidal frequency ωM2 (the black dashed line indicates this turning depth). The internal wave angle of propagation,
The local buoyancy frequency N(z) (black circles, black axes on the left and bottom) decreases by more than two orders of magnitude from the shallow ocean down to 5 km depth in the Venezuelan Basin in the Caribbean. The data are taken from cast A22-316N151-4-10 from the World Ocean Circulation Experiment, latitude 13.17° N, longitude 66.00° W, and analyzed (averaging over 50 m depths) by the method described by King et al.11 For depths greater than 4 km, the local buoyancy frequency becomes less than the semi-diurnal lunar tidal frequency ωM2 (the black dashed line indicates this turning depth). The internal wave angle of propagation,
To investigate the relation of the radiated internal wave power to the corresponding boundary current strength, and also the consequences of weak stratification where the hydrostatic approximation fails, we conduct two-dimensional numerical simulations of the Navier-Stokes equations for tidal flow of a uniformly stratified fluid (N = const) over triangular, Gaussian, knife-edge, and semi-circular topographies. We vary the ratio of the tidal frequency ω to the buoyancy frequency N from values comparable to unity, where nonhydrostatic effects are expected to be significant, to values an order of magnitude smaller, where our measurements may be directly compared with prior hydrostatic theory. We examine the scalings and interplay of the radiated internal wave power and the boundary currents by independently varying the tidal frequency, buoyancy frequency, topographic slope, and topographic shape.
In Sec. II, we summarize past studies of resonant boundary currents and internal wave generation by tidal flow over two-dimensional topography. Our numerical methods and the details of the parameter variation are described in Sec. III. The results are presented in Sec. IV and are discussed in Sec. V.
II. BACKGROUND
A. Generation of boundary currents and internal tides by critical topography
Cacchione et al.2 observed that the topographic slopes off the coasts of California and New Jersey are approximately equal to the local slopes of internal waves, that is, the topography is critical, ε = 1. They suggested that this critical topography resonantly generates intense boundary layer currents that can carry sediment, thus shaping the continental slopes. The resonant generation of boundary currents by critical topography was examined in experimental studies of tidal flow (over a two-dimensional slope) by Zhang et al.,12 who found that at resonance (ε = 1) the boundary current velocity was more than an order of magnitude larger than the maximum tidal velocity. At resonance, the measured maximum boundary current velocity was found to vary as L,4/3 where L was the length of the topography's near-critical region. For sufficiently large L, the boundary layer currents became unstable and formed Kelvin-Helmholtz billows. The increase in the resonant current with increasing slope length was also observed in numerical studies of Gayen and Sarkar,13 who also observed that 3D instabilities could occur along the continental slope, resulting in spanwise wave breaking and turbulent mixing.14,15
Critical topography has also been identified as being particularly effective in converting tidal motions to internal tides by Lamb.16 In a review of internal tide generation, Garrett and Kunze17 argued that intense internal waves would be radiated by slopes near criticality, and these internal waves would lead to local mixing; however, Garrett and Kunze17 also point out that most of the energy flux is expected to reside in low vertical modes. Griffiths and Grimshaw18 used a modal decomposition method to argue that there is a resonant forcing of internal waves by critical topography along a continental slope. In similar experimental studies, Lim et al.19 found that intense internal wave beams were radiated only when a critical slope was present along their continental shelf/slope system. Zhang et al.12 and Gayen and Sarkar13 both also suggested that critical topography leads, in addition to the resonant boundary currents described above, to intense internal waves that propagate away from the continental slope. However, no past work has examined the relationship between the radiated internal wave power and the boundary current intensity as a function of the relevant parameters of the problem.
B. Internal tide generation: Past theoretical studies
Garrett and Kunze17 have reviewed theoretical approaches to understanding internal tide production. Baines20 used the method of characteristics, which is broadly applicable yet difficult to implement. Analytical predictions for the internal wave power generated by inviscid tidal flow in two dimensions have been obtained assuming the hydrostatic approximation (ω ≪ N), in the limits of weak topography by Bell21 and steep topography by Llewellyn Smith and Young.22 Bell21 predicted that in the limit ε → 0 in an infinitely deep ocean
This prediction was extended by Llewellyn Smith and Young23 and Khatiwala24 to account for finite-depth effects, which were shown to reduce the radiated power. Even though these linear, inviscid theoretical results become singular for critical (ε = 1) and supercritical topography (ε > 1), Balmforth et al.25 showed that the power radiated by subcritical topography (ε < 1) in the deep ocean is proportional to H2 and smoothly increases only modestly as the topography steepens from Bell's limit of ε ≪ 1 towards the critical case of ε = 1.
To circumvent the singularities that arise for critical and supercritical topography (ε ⩾ 1), Llewellyn Smith and Young22 used a Green's function approach26 to compute the conversion of the barotropic tide into the internal tide for knife-edge topography (where ε → ∞)
Pétrélis et al.27 used the Green's function method to perform hydrostatic calculations that showed the power smoothly varies from Bell's21 prediction for subcritical topography (Eq. (5)) towards Llewellyn Smith and Young's22 result for supercritical topography (Eq. (6)). The tidal conversion in the extreme cases of ε → ∞ predicted by Llewellyn Smith and Young22 and ε → 0 described by Bell21 only differ by a factor of two. The radiated power shows a similar modest increase as the ratio of the topographic height to the fluid depth approaches unity,22,28 gaining only a factor of two as this ratio increases from 0 to 0.92 when the scaling with H2 is taken into account.
These analytical predictions were shown to agree with hydrostatic numerical simulations by di Lorenzo et al.,29 despite the singularities that arise in linear, inviscid analytical calculations for critical and supercritical topography. di Lorenzo et al.29 also highlighted the importance of resolving the topography, as smoothing the topography acted to reduce estimates of the radiated power. This point was emphasized by Zilberman et al.,30 whose model yielded estimates of tidal conversion for the Mid-Atlantic Ridge that decreased by as much as 20% when the topography was smoothed.
A Green's function26 approach has also been used to predict the power radiated by tidal flow over more complicated topographic structures. An infinite series of evenly spaced knife edges was considered (for an inviscid fluid) by Nycander,31 who found that the radiated internal wave power depends strongly upon the height and spacing of the ridges, and diverges for particular parameter values. Similarly, Balmforth and Peacock32 calculated for periodic Gaussian, sinusoidal, and sawtooth topographies the dependence of the power on the height and spacing of the peaks. These two studies of internal wave power radiated by periodic topography showed that estimates of the average radiated power are independent of the height, which is in contrast to the H2 dependence found for isolated topographic features. The requirement for the topography to be spatially periodic was removed by the studies of Echeverri and Peacock,33 where the Green's function method was extended to arbitrary two-dimensional topography.
Griffiths and Grimshaw,18 using modal decomposition, developed a linear, hydrostatic model that enabled them to consider internal tide generation by arbitrary three-dimensional topography for general background stratifications. Zarroug et al.34 instead accounted for nonuniform stratifications using homogenization theory, which was shown to be more appropriate than the WKB approximation for low vertical modes (long vertical wavelength variations). Since the WKB approximation uses the buoyancy frequency at the bottom boundary, the method may underestimate the global tidal conversion because local mixing over the topography can reduce the buoyancy frequency.31
C. Internal tide generation: Past experimental and numerical studies
Internal tide generation by two-dimensional topography has also been the focus of a number of experimental and numerical studies. In addition to the aforementioned analytical studies, Khatiwala24 performed nonhydrostatic numerical simulations of internal tide generation by truncated sine or Gaussian ridges, focusing on finite-depth effects and the variation of the power with the criticality parameter ε. Khatiwala24 found that finite-depth effects can greatly reduce the radiated power for topography that is wide along the direction of the tidal flow (compared to the wavelength of a mode-1 internal wave), and that the radiated power can saturate for certain supercritical topography. The generation of harmonics and wave-breaking by nonlinear beam-beam interactions were simulated for supercritical topography by Lamb.16 Gostiaux and Dauxois35 experimentally examined the generation of internal tides by a continental shelf break, which was shown to be analogous to internal wave generation by an oscillating cylinder. Along the same lines, Zhang et al.36 oscillated a cylinder horizontally to experimentally simulate the generation of internal waves by supercritical topography. Their results agreed well with the viscous theoretical predictions of Hurley and Keady37 in the linear regime, while higher harmonics were generated for stronger forcing. Internal tide generation by knife edge and Gaussian topography was experimentally shown to agree with prior theory by Peacock et al.38 when viscous dissipation, which smoothed the experimental wavefields and suppressed instabilities, was accounted for in the predictions. The modal structure of internal tides generated by Gaussian topography was experimentally and numerically shown by Echeverri et al.39 to agree with linear theory only for small tidal excursions (∼1% of the topographic width), while higher modes were strongly affected for larger tidal excursions. Qian et al.40 performed numerical simulations that showed that the presence of a strong pycnocline above ridges could greatly enhance the radiated internal wave power, even if the stratification in the vicinity of the topography was weak. Lim et al.19 performed experiments of tidal motions over a continental shelf/slope system, which produced internal wave beams only when a critical slope was present, while boluses could form for other parameter values.
D. Nonhydrostatic effects and turning depths
The work discussed above has either been restricted to or primarily focused on situations where the hydrostatic approximation is valid (ω ≪ N). Approximately half of the internal wave energy in the ocean is at the semi-diurnal lunar tidal frequency1,41 ωM2, which is typically much less than characteristic values of the buoyancy frequency in the shallow ocean. In such cases, internal waves propagate nearly horizontally (small SIW) and have small vertical length scales compared to their horizontal extent. Therefore, applying the hydrostatic approximation to study internal tide generation by topography in the shallow ocean or mid-ocean ridges is valid.
The buoyancy frequency, however, decreases by two orders of magnitude or more from the shallow ocean to the abyss (see Fig. 1). King et al.11 analyzed temperature, pressure, and salinity data from the World Ocean Circulation Experiment (WOCE) for thousands of locations throughout the oceans, and they found that the local buoyancy frequency becomes smaller than ωM2 for many locations in the deep ocean. Internal waves reflect from turning depths where N(z) = ωM2, and become evanescent (exponentially damped) below.42,43 The internal wave slope SIW diverges at a turning depth, where the direction of propagation becomes vertical (
The analytical predictions of Bell21 and Llewellyn Smith and Young,22 given by Eqs. (5) and (6), were derived under the hydrostatic approximation, where the
For ω ≪ N, the nonhydrostatic expressions (7) and (8) reduce to the hydrostatic predictions of Bell21 and Llewellyn Smith and Young22 in Eqs. (5) and (6). As in the hydrostatic limit, the power radiated by infinitely steep topography (ε → ∞) is predicted to be twice that predicted for infinitely shallow topography (ε → 0).
For both hydrostatic and nonhydrostatic cases, the radiated internal wave power is predicted to be proportional to the effective tidal power Ptide (see Eq. (2)) that interacts with the topography. In the hydrostatic limit, the power is predicted to be proportional to N/ω, while in the nonhydrostatic case the power is predicted to be proportional to the inverse beam slope (1/SIW). Hydrostatic and nonhydrostatic predictions differ particularly for the weak stratifications near a turning depth. As N → ωM2, hydrostatic theory predicts the internal wave power to be comparable to the tidal power, while nonhydrostatic theory predicts zero radiated power in this case. For the data in Fig. 1, hydrostatic theory overestimates the radiated power by only 3% at a depth of 2 km, but the error increases to 65% at 3 km depth and diverges at the turning depth, 4 km. Therefore, while a topographic ridge with a height H ∼ 1 km would typically be expected to efficiently convert barotropic tidal motions to internal tides, the radiated power could be significantly decreased in locations in the abyssal oceans where N ∼ ωM2.
III. METHODS
We simulate internal wave generation by tidal flow over two-dimensional topography in a uniformly stratified fluid by solving the Navier-Stokes equations in the Boussinesq approximation. The code is CDP 2.4,44 a parallel, finite-volume-based solver (the subgrid modeling is disabled). Second-order accuracy in space and time is achieved using a fractional-step time-marching scheme and multiple implicit schemes for the spatial operators.45 The equations are solved for the density ρ, pressure p, and velocity
where ρ0 = 1000 kg/m3, ν = 0.01 m2/s (unless otherwise specified) is the kinematic viscosity, and D = 2 × 10−5 m2/s is the salt diffusivity, resulting in a Schmidt number of ν/D = 500. We also simulate cases with ν = 0.003 and 0.03 m2/s and find that the normalized power
The computational domain, generated with Pointwise Gridgen, spans −40 < x < 40 km and 0 < z < 8 km. Different topographic shapes (described below) are centered at x = 0 km with the bottom at z = 0 km. The structured grid is composed of rectangular control volumes with smoothly varying spatial resolution. The horizontal resolution along the topography changes from Δx = 10 m in the center to Δx = 200 m for |x| > 35 km. The boundary layer is resolved using Δz = 0.8 m for the first 200 m above the topography, and then the resolution is monotonically changed to Δz = 15 m for z > 2 km. The domain has approximately 106 control volumes. No-slip boundary conditions are enforced along the upper and lower boundaries, while periodic boundary conditions are used in the x-direction. To mimic an infinitely deep ocean and to minimize finite-size effects, we apply a Rayleigh damping force proportional to the deviation from the tidal velocity (
Simulations are conducted for several symmetric topographies of height H = 0.5 km: multiple triangular-shaped ridges, a Gaussian ridge, a knife edge, and a semi-circular topography; see Table I for the range of parameters for each topographic shape. The triangular ridges have constant slope flanks (0.2 ⩽ Stopo ⩽ 11.25) connected at the peak by a circular cap of radius r = 62.5 m and to the bottom by a smooth parabolic curve of height 40 m (see Fig. 2(a)). The Gaussian ridge has max (Stopo) = 1. The knife edge has thickness 16 m and the semi-circle has a radius of 0.5 km.
(a) The radiated internal wave power is obtained by integrating the energy flux normal to a cross-section (dashed line) that is a distance δ from the topography (gray). The triangular topography here has slope Stopo = 1 and δ = 100 m. (b) The power PIW radiated through a cross-section a distance δ away from the topography reveals the decay of the internal waves as they propagate away from the topography. (Triangular topography with SIW = Stopo = 1 and ε = 1.)
(a) The radiated internal wave power is obtained by integrating the energy flux normal to a cross-section (dashed line) that is a distance δ from the topography (gray). The triangular topography here has slope Stopo = 1 and δ = 100 m. (b) The power PIW radiated through a cross-section a distance δ away from the topography reveals the decay of the internal waves as they propagate away from the topography. (Triangular topography with SIW = Stopo = 1 and ε = 1.)
Numerical simulation parameters. Here, Stopo is the maximum topographic slope, SIW is the slope of the internal wave beam, ε = Stopo/SIW is the criticality parameter, ω is the tidal frequency, and N is the buoyancy frequency.
Topographic shape . | Stopo . | 1/SIW . | ε . | ω (10−4 rad/s) . | N (10−4 rad/s) . |
---|---|---|---|---|---|
Triangular | 0.2–11.25 | 0.18–10 | 2 | 0.154–1.526 | 1.55 |
Triangular | 0.63–11.25 | 0.18–3.2 | 2 | 1.4052 | 1.11–3.64 |
Triangular | 0.2–2 | 0.1–1 | 0.2 | 1.096–1.542 | 1.55 |
Triangular | 1 | 0.1–10 | 0.1–10 | 0.154–1.542 | 1.55 |
Triangular | $1/\sqrt{3}$ | 0.17–17 | 0.1–10 | 0.089–1.527 | 1.55 |
Triangular | 2 | 0.1–5 | 0.2–10 | 0.304–1.542 | 1.55 |
Gaussian | 1 | 0.2–10 | 0.2–10 | 0.154–1.520 | 1.55 |
Knife edge | ∞ | 0.14–10 | ∞ | 0.154–1.535 | 1.55 |
Semicircle | ∞ | 0.2–10 | ∞ | 0.154–1.520 | 1.55 |
Topographic shape . | Stopo . | 1/SIW . | ε . | ω (10−4 rad/s) . | N (10−4 rad/s) . |
---|---|---|---|---|---|
Triangular | 0.2–11.25 | 0.18–10 | 2 | 0.154–1.526 | 1.55 |
Triangular | 0.63–11.25 | 0.18–3.2 | 2 | 1.4052 | 1.11–3.64 |
Triangular | 0.2–2 | 0.1–1 | 0.2 | 1.096–1.542 | 1.55 |
Triangular | 1 | 0.1–10 | 0.1–10 | 0.154–1.542 | 1.55 |
Triangular | $1/\sqrt{3}$ | 0.17–17 | 0.1–10 | 0.089–1.527 | 1.55 |
Triangular | 2 | 0.1–5 | 0.2–10 | 0.304–1.542 | 1.55 |
Gaussian | 1 | 0.2–10 | 0.2–10 | 0.154–1.520 | 1.55 |
Knife edge | ∞ | 0.14–10 | ∞ | 0.154–1.535 | 1.55 |
Semicircle | ∞ | 0.2–10 | ∞ | 0.154–1.520 | 1.55 |
The internal wave power radiated to the far field is obtained for the triangular, semi-circle, and knife edge topographies by integrating the energy flux along paths with the same shape as the topography; for the Gaussian ridge, we use a triangular-shaped path with flanks of slope 1. The qualitative behavior of the radiated power is insensitive to the choice of path. The integration path for obtaining the energy flux is located at a distance δ from the topography, as illustrated by Fig. 2(a). The radiated power measured through cross-sections of increasing δ for triangular topography with SIW = Stopo = 1 and ε = 1 are shown in Fig. 2(b). The radiated power decays rapidly inside the laminar boundary layer (δ/H < 0.1), and then transitions to a gradual decay as the internal wave beam propagates away from the topography (δ/H > 0.1). We therefore choose a value of δ = 100 m, which corresponds to δ/H = 0.2 as shown by the dashed red line in Fig. 2(b), which is outside of the boundary layer yet close enough to the topography so that viscous decay is negligible.
IV. RESULTS
Figure 3 illustrates that resonant boundary currents form over critical topography, but they do not enhance the internal wave power radiated to the far-field. Here, we compare the kinetic energy density and the energy flux in the direction of internal wave propagation (
Strong boundary currents form for critical topography (ε = 1), (c) and (d), while the radiated internal waves are stronger for supercritical topography (ε = 2), (e) and (f). Subcritical topography
Strong boundary currents form for critical topography (ε = 1), (c) and (d), while the radiated internal waves are stronger for supercritical topography (ε = 2), (e) and (f). Subcritical topography
A. Relation of boundary currents to radiated wave power
We examine the connection between the boundary currents and the tidal conversion by computing both the integrated kinetic energy density in the boundary currents, as in Zhang et al.,12 and the power radiated to the far-field. The spatially- and tidally averaged kinetic energy density
(a) Resonant boundary currents form for tidal flow over triangular topography (Stopo = 1) near criticality (ε ≈ 1), as revealed by the peak in
(a) Resonant boundary currents form for tidal flow over triangular topography (Stopo = 1) near criticality (ε ≈ 1), as revealed by the peak in
The resonant peak for the boundary current kinetic energy density is enhanced (diminished) by decreasing (increasing) the fluid viscosity, which acts to increase (decrease) the Reynolds number of the flow. The resonant peak also increases as the topographic length is increased. In the laboratory experiments of Zhang et al.,12 which were at higher Reynolds numbers than we simulated here, the kinetic energy density was found to be approximately 60 times the maximum tidal kinetic energy density. These experimental observations were complemented by numerical studies by Gayen and Sarkar,13 who focused on the scaling with the slope length. The Reynolds numbers of oceanic flows far exceed those achieved in these past studies as well as the work presented here, indicating that the resonant boundary currents can be significantly stronger in the ocean.2
While the boundary currents are strongest near ε = 1, the internal wave power monotonically increases with the criticality parameter, as Fig. 4(b) illustrates. We multiply
Direct measurements of the energy flux of internal wave beams are rare in laboratory experiments and oceanic measurements because of the difficulty in determining the pressure field. As an alternative, the kinetic energy density is often used to characterize the strength of radiated internal waves. Figure 5 compares the kinetic energy density and energy flux for the critical and supercritical topography shown in Figs. 3(c)–3(f). Values of the kinetic energy density for the top internal wave beams (near σ/H = 2.2) are higher than for the beams that initially propagate downward and reflect from the bottom boundary (peaked near σ/H = 0.6). Further, the kinetic energy density is largest for the internal wave beam generated by critical topography. However, the behavior of the energy flux is different. The peak value of the energy flux corresponds to the internal wave beam generated by supercritical topography that reflects off of the bottom boundary (peak in Fig. 5(b) near σ/H = 0.8). While this internal wave beam yields the largest values of energy flux and radiated power, the kinetic energy density only has a modest peak compared to the others that did not reflect from the bottom boundary. The kinetic energy density in the area between the peaks of the internal wave beams (1 < σ/H < 2) is very similar for critical and supercritical topography, although the energy flux in the supercritical case is higher over this entire range. Finally, the largest peak in the kinetic energy density (σ/H ≈ 2.2) of critical topography does correspond to a maximum in the energy flux, but this is not the case for supercritical topography. We therefore conclude that the kinetic energy density should not be used as a proxy for the energy flux to characterize the conversion of tidal motions to internal waves, as the two quantities do not have a consistent relationship.
A comparison of the tidally averaged (a) kinetic energy density and (b) energy flux in the direction of the internal wave beams for the critical topographies (black circles) and supercritical topographies (squares) shown in Figs. 3(c)–3(f). The measurements are taken along a cross-section from (x/H = −2, z/H = 0) to (x/H = 0, z/H = 2), which is shown as the dashed line in panel (c) of Fig. 3.
A comparison of the tidally averaged (a) kinetic energy density and (b) energy flux in the direction of the internal wave beams for the critical topographies (black circles) and supercritical topographies (squares) shown in Figs. 3(c)–3(f). The measurements are taken along a cross-section from (x/H = −2, z/H = 0) to (x/H = 0, z/H = 2), which is shown as the dashed line in panel (c) of Fig. 3.
B. Radiated power dependence on beam slope
Many of the prior studies described in Sec. II focused on the scaling of the radiated power with the criticality parameter, topographic height, and the ratio of the topographic height to the total fluid depth. In the deep ocean, the internal wave beam slope increases rapidly with depth, particularly in the vicinity of a turning depth where it diverges (cf. Fig. 1). We now examine the scaling of the radiated power with beam slope for cases where the hydrostatic approximation is valid and for cases where the approximation fails.
The normalized power
(a) The normalized radiated power
(a) The normalized radiated power
We first compare our measurements to the hydrostatic theories of Bell21 (Eq. (5)) and Llewellyn Smith and Young22 (Eq. (6)), which are also shown in Fig. 6(a). As expected, the theories fail to predict the power for small inverse beam slopes (1/SIW < 1), which correspond to cases with weak stratification where ω → N. Our measurements indicate that the power is inversely proportional to the beam slope for all SIW. While
The nonhydrostatic extensions of the theories of Bell21 and Llewellyn Smith and Young,22 given by Eqs. (7) and (8), each predict P ∝ 1/SIW for all SIW, in accord with our measurements in Fig. 6(a). Analogous to the hydrostatic simulations of Pétrélis et al.,27 the nonhydrostatic prediction of Bell's theory agrees well with our simulations for small ε, while the nonhydrostatic version of Llewellyn Smith and Young's theory yields accurate predictions at large ε, as Fig. 6(b) illustrates. For supercritical topography with very steep beam slopes (1/SIW < 0.3), we find that the radiated power begins to decrease owing to interactions between the internal waves generated on the left- and right-hand sides of the topography.
C. Dependence on topographic shape
Thus far we have considered the internal tidal power radiated by triangular topography with constant slope flanks. This is certainly an idealization of rough ocean topography, which is characterized by many different slope angles and length scales. To extend beyond such simple topography, we compare the results for triangular topography (slope Stopo = 1) to the results for Gaussian, knife-edge, and semi-circular ridges in Fig. 7, thus including cases where the slope smoothly varies over a fixed range (Gaussian), is infinitely steep (knife edge), or spans 0 < Stopo < ∞ (semi-circle). For internal waves with a shallow slope (large 1/SIW), the radiated power in all cases approaches the nonhydrostatic prediction of Llewellyn Smith and Young,22 given by Eq. (8). For steep beam slopes (small 1/SIW), the radiated power slowly tends toward the prediction of the nonhydrostatic extension of the Bell theory. Thus, the nonhydrostatic theories are capable of predicting the radiated power for topography with varying slope angles, in addition to the simple triangular topography discussed above.
The normalized power
The normalized power
The results for Gaussian and triangular ridges with the same maximum slope (Stopo = 1) are very similar, indicating that using the maximum topographic slope to determine the criticality parameter is appropriate. For shallow beams, the knife-edge and semi-circular topographies yield results similar to those for the Gaussian and triangular ridges, as all of the topographic shapes are supercritical. The decrease in power with increasing beam slope, though, is slower for the knife-edge and semi-circular topographies owing to their steeper topographic slopes.
One would expect the power radiated by the knife-edge to closely agree with
V. DISCUSSION
Our simulations show that critical topography generates resonant boundary currents, as found in prior studies.2,12–15 However, this resonance phenomenon does not extend to radiated internal waves that propagate away from the topography. Prior studies have found that this resonance behavior can produce boundary currents that are strong enough to undergo shear instability, leading to turbulence and local mixing,12–15 which can actually reduce the radiated internal wave power rather than enhancing it. Therefore, the turbulence produced by unstable resonant boundary currents over critical topography in the ocean may actually decrease the internal wave generation, in contrast to the notion that critical topography generates intense internal waves.
We have also found that the kinetic energy density and radiated internal wave power are largely decoupled. Decreasing (increasing) the viscosity leads to an increase (decrease) in the kinetic energy density of the resonant boundary currents, while the radiated internal wave power remains almost unchanged. Even though the viscous effects in our simulations are stronger than in the ocean, we believe that the decoupling of boundary current kinetic energy density and the radiated power extends beyond our parameter values. Indeed, our measurements of the radiated power agree quantitatively with the inviscid theoretical predictions of Bell21 and Llewellyn Smith and Young22 for both small and large values of the criticality parameter. Viscous effects are minor in our simulations because we examine the radiated power in the region immediately beyond the laminar boundary layer, whereas experimental studies of internal wave beams farther away from the topography found that viscous effects had to be taken into account.19,38,39 Future studies could bridge the gap between our focus on varying parameters over wide ranges and the studies that have examined higher Reynolds numbers where the boundary currents become turbulent.12–15,46
In summary, the power radiated by internal waves produced by tidal flow over two-dimensional ridges in a uniformly stratified fluid can be described by
where π/8 ⩽ f (ε, shape) ⩽ π/4 depends upon the details of the topography and is bounded by the nonhydrostatic versions of the linear, inviscid theories of Bell21 and Llewellyn Smith and Young.22 From Eq. (12), we see that variations in the criticality parameter ε can only change the power by a factor of two, whereas the power scales as the inverse of the beam slope. This scaling with the beam slope holds for cases where the hydrostatic approximation is appropriate as well as situations where it fails. The beam slope becomes large in the deep ocean, as observed in the Bay of Biscay by Pingree and New,47 and even diverges if there is a turning depth.42,43 The nonhydrostatic effects that arise from weak stratifications reduce the radiated power below hydrostatic estimates; this should be considered in making global estimates of the internal wave power radiated by deep ocean topography.
The present work has examined two-dimensional laminar tidal flow over smooth topography in uniformly stratified fluids. However, in the ocean the topography is rough and three-dimensional, and the fluid is not uniformly stratified.11 Holloway and Merrifield48 performed three-dimensional hydrostatic simulations of tidal motions over steep topographic features and found that the tidal conversion strongly increased when the length of the topography orthogonal to the tidal flow increased, which forced the flow to cross isobaths rather than flowing around the topography. While two-dimensional studies can serve as a foundation for understanding tidal conversion by topographic features that efficiently convert tidal motions to radiated internal waves, the effects of three-dimensionality cannot be neglected. Tidal conversion has been examined in three-dimensional simulations of the effects of corrugations in a continental slope, for both cross- and along-slope tidal forcing.49,50 Studies of tidal flow over a half sphere revealed both a mean flow and an oscillatory flow (at twice the tidal frequency) orthogonal to the tidal direction.51 Along similar lines, King et al.52 found that tidal motions over a three-dimensional Gaussian mountain produced internal wave beams (at twice the tidal frequency) that propagated orthogonal to the tidal direction, but only when the topography was supercritical. Gayen and Sarkar14 observed three-dimensional instabilities produced by resonant boundary currents over two-dimensional critical topography. It will be interesting in the future to search for resonant boundary currents in tidal flow over three-dimensional topography.
ACKNOWLEDGMENTS
We thank Likun Zhang, Bruce Rodenborn, Philip Morrison, and Robert Moser for helpful discussions. The computations were done at the Texas Advanced Computing Center. The research was supported by the Office of Naval Research (ONR) MURI Grant No. N000141110701.