1. Introduction
Interest in internal tides has been revived by recognition that up to one-third of the surface tide dissipation occurs in the deep ocean (Egbert and Ray 2001; Kantha and Tierney 1997) and that this tidal energy could contribute significantly to the roughly 2 TW of abyssal mixing thought to be needed to close the upwelling branch of the meridional thermohaline circulation (Munk and Wunsch 1998); the remainder presumably comes from wind forcing into near-inertial motions (e.g., Alford 2003) and into the mean flow (Wunsch 1998). Assimilating satellite altimetry data into a barotropic tide model, Egbert and Ray (2001) found that 0.8 TW of the surface tide was lost at ridges, trenches, and island chains in the deep ocean. Two-layer numerical simulations (Simmons et al. 2004) show that 75% of the deep-ocean surface tide loss enters the low-mode wave field through just 12 regions of abrupt topography. Satellite altimetry also reveals low-mode internal tides radiating from these same topographic features (Ray and Mitchum 1997). While most tide–topography interaction theory is posed in the weak-amplitude limit with subcritical slopes and small heights (Bell 1975; St. Laurent and Garrett 2002; Llewellyn Smith and Young 2002), features like the Hawaiian Ridge are largely characterized by steep supercritical slopes and topography that occupy a significant fraction of the water column. For these cases, more complicated theories (Baines 1974, 1982; St. Laurent et al. 2003; Llewellyn Smith and Young 2003) or numerical modeling (Holloway and Merrifield 1999; Merrifield et al. 2001; Johnston and Merrifield 2003; Simmons et al. 2004) are required.
The Hawaiian Ridge is the site of the most intense surface tide loss per unit area in the global deep ocean (Egbert and Ray 2001; Simmons et al. 2004). Numerical modeling (Merrifield et al. 2001) and in situ measurements (Lee et al. 2006; Rudnick et al. 2003) find particularly intense internal tide energy fluxes radiating from French Frigate Shoals and Kaena Ridge in the Kauai Channel west of Oahu.
Few observations detail the structure of the deep-ocean internal tide where it is generated [for a review see Huthnance (1989)]. In the Bay of Biscay, Pingree and New (1989, 1991) found internal tide beams radiating from the continental shelf break. Holloway et al. (2001) found that the internal tide on the Australian northwest shelf is dominated by a vertically standing low mode that is dissipated close to its generation site. In Monterey Canyon, Petruncio et al. (1998) found that internal tides generated at the canyon mouth propagate either upcanyon or form horizontally standing modes, depending on the stratification. Kunze et al. (2002) computed baroclinic energy fluxes in the canyon and found them to be upcanyon and to decay onshore. At the Mendocino Escarpment, Althaus et al. (2003) found beams of energy radiating from the 1500-m-deep ridge; strong energy flux divergence characterized the generation site. In each of these studies, the observations were generally consistent with energy flow along internal tide characteristics, in accord with linear wave propagation. Highly nonlinear generation has also been observed in association with strong topographic flow over relatively shallow bathymetry (Chereskin 1983; Hibiya 1988; Loder et al. 1992; Farmer and Smith 1980) and from river plumes (Nash and Moum 2005).
In this paper, we describe the formation and structure of the internal tide across the southern half of Kaena Ridge (Fig. 1). Repeated full-depth profiles of velocity and density are used to examine the barotropic-to-baroclinic tidal conversion, the semidiurnal kinetic and potential energies, phase propagation, and energy flux. Low-vertical-wavenumber waves (modes 1–10) resolved by the sampling array carry most of the energy flux. Turbulent dissipation associated with higher vertical wavenumbers is presented by Klymak et al. (2006) and Lee et al. (2006).
This paper is organized as follows. In section 2, we describe the observations and analysis methods. The structure of the velocity and vertical displacements is described in section 3. In section 4, we present observations of the energetics. In section 5, we reconcile the generation and propagation pathways using two simple models. This is followed by conclusions in section 6.
2. Observations
During a spring tide in October 2000, a five-station section across Kaena Ridge was sampled from ridge crest (at about 1000-m depth) to the 3000-m isobath (Fig. 1). To obtain snapshots of the wave field that resolve the semidiurnal internal tide, stations 1–5 were sampled sequentially with expendable current profilers (XCPs; Sanford et al. 1993) for 21 h so that each station was occupied seven times with a 3-h repeat interval. Thirty-five XCPs were deployed with the ship under way, each returning baroclinic horizontal velocity (u, υ) and temperature T with 2-m vertical and 4.5-km horizontal resolution.1 Velocity profiles were made absolute using 30-min-averaged shipboard ADCP profiles (RD Instruments, 150 kHz; GPS referenced) over the upper 150 m.
At station 1, the 3000-m water depth exceeds the range of XCPs, so six full-depth profiles of u, υ, T, salinity S, and microscale shear were collected using the absolute velocity profiler (AVP; Sanford et al. 1978, 1985) during the 15-h period immediately prior to the 21-h XCP surveys. AVP velocity profiles are collected during both descent and ascent, yielding 12 profiles of u and υ. In the upper 2000 m, XCP and AVP data have statistically similar semidiurnal phase and amplitude of u, υ, and T, allowing us to merge the data after applying a 24.8-h offset (i.e., Figs. 2e and 3e); XCP data from station 1 are not presented in this paper.
In the following we use a coordinate system aligned with the bathymetry (and our transect) so that cross ridge (x and u) is 37.4° east of north and along ridge (y and υ) is 37.4° north of west.
a. Methods
The temporal evolution and spatial structure of the baroclinic perturbation fields of velocity, vertical displacement, and pressure (u′, ξ, and p′) are used to diagnose the dynamics of the semidiurnal internal wave field. From these, the energy flux FE = 〈u′p′〉 is computed as the covariance of the wave-induced pressure p′ and velocity u′. By definition, p′ and u′ are baroclinic fluctuations, so have vanishing temporal and vertical averages. Our methods for computing these fluctuations, outlined below, follow Kunze et al. (2002), Althaus et al. (2003), and Nash et al. (2005).
Only temperature and horizontal velocity were measured by the XCPs. Four CTD profiles obtained before and after the XCP surveys had an unambiguous T−S relation and were used to infer S from T. To test its validity, the relationship was applied to the AVP T data, for which S was also measured but not used in the derived S(T) relation. Based on this comparison, an rms error of 0.013 psu in S and 0.010 kg m−3 in potential density is associated with this approximation. Based on the mean stratification, this corresponds to errors in vertical displacement of ξerr = 12 m at z = 1000 m and ξerr = 30 m at z = 1500 m.

Profiles of isopycnal displacement, perturbation velocity and pressure are used to estimate the baroclinic energy flux FE = 〈u′p′〉, baroclinic horizontal kinetic energy density HKE = 〈u′2 + υ′2〉/2, and available potential energy density APE = N 2〈ξ2〉/2. The angle brackets represent depth-dependent temporal averages, as computed from vertical profiles of semidiurnal amplitude and phase (described below).
Monte Carlo–type simulations are used to estimate how a broadband Garrett–Munk (GM) wave field alters results obtained from short-duration, discretely sampled time series. Following Nash et al. (2005), we generate an ensemble of 100 independent realizations of randomly phased, GM-contaminated wave fields with the observed semidiurnal energetics. Each realization is then sampled as were the observations; standard deviations may then be computed from the ensemble of these regenerated samplings. These error estimates account for broadband aliasing but do not account for aliasing by narrowband diurnal and near-inertial signals. Cospectral analysis of u and p′ from a long, rapidly sampled time series (Rainville and Pinkel 2006) indicates that the semidiurnal waveband dominates the energy flux; diurnal fluctuations contribute to a lesser extent, accounting for ∼15% of the total energy flux.
3. Structure of velocity and isopycnal displacement
Time series of perturbation velocity and vertical displacement constructed from individual XCP and AVP profiles are shown in Figs. 2 and 3. Stations 5 and 4 (Figs. 2a,b and 3a,b) represent observations on the ridge crest, while stations 3, 2, and 1 (Figs. 2c–e and 3c–e) are measurements over the steep ridge flank.
Amplitude and phase were computed over 50-m vertical intervals for each cross-ridge velocity, along-ridge velocity, vertical displacement, and pressure perturbation. The semidiurnal vertical velocity fluctuation w′ was estimated from the fits as dξ/dt. Semidiurnal fits capture roughly 70% of the variance of the cross-ridge perturbation velocity (i.e., ≈ 2.5
). Semidiurnal fits to ξ(z, t) capture >90% of the isopycnal displacement variance off ridge (stations 1 and 2), but only 30% over the ridge crest (stations 4 and 5), where the signal is weak.
Time-averaged, subtidal currents
a. Phase propagation of cross-ridge velocity
Raw time series of cross-ridge baroclinic velocity u′ and its semidiurnal phase ϕ(u) and amplitude uo are shown in Fig. 2. The velocity is strongest and exhibits the simplest phase propagation at the shallowest stations, becoming progressively weaker and more complicated moving farther off ridge.
Semidiurnal barotropic velocity amplitudes were computed by applying harmonic analysis (5) to the depth-mean observations. The amplitude of the barotropic tide’s cross-ridge velocity uobt increases from 0.04 m s−1 at 3000 m to 0.19 m s−1 at the ridge crest (Fig. 2, thick vertical gray lines in rightmost panels). Barotropic velocity amplitudes exceed baroclinic over the ridge crest, while the reverse is true on the ridge flanks (station 1). Barotropic phase was estimated as ϕ(u)bt = 254 ± 3°, where ± represents variability between stations. In comparison, the phase of the barotropic tide from the TPX0.5 model (Egbert 1997) is 262 ± 2°, with amplitudes ranging from 0.05 m s−1 at station 1 to 0.18 m s−1 at station 5.
Over the ridge crest (stations 4 and 5; Figs. 2a,b), the semidiurnal phase propagation of upslope (cross ridge) velocity is downward and increases approximately linearly with depth, indicating upward energy propagation with a vertical wavelength of 1000 m. The semidiurnal amplitude is relatively constant with depth at ∼0.16 m s−1.
Over the steep slopes (stations 2 and 3; Figs. 2c,d), downward phase (upward energy) propagation is observed in the upper 600–700 m, similar to that over the ridge crest. This is mirrored by upward phase (downward energy) propagation below 700 m. Farthest off ridge (station 1, Fig. 2e), the phase structure becomes more complex. As with stations 2 and 3, downward (upward) phase propagation is observed in the upper (lower) part of the water column. However, there are also hints of vertically standing behavior, with a roughly 180° phase shift between upper and lower water column velocities. The semidiurnal velocity amplitude uo has significant vertical structure, exhibiting both surface and bottom intensification with range from 0.02 to 0.15 m s−1.
b. Phase propagation of vertical displacement
In contrast to the baroclinic velocity, which is strongest and exhibits the clearest phase propagation over the ridge crest, vertical displacements are only significant over the steep flanks. There, the direction of phase propagation for ξ is qualitatively similar to that of u′—downward above 1000 m and upward below 1000 m. In section 5, we suggest that a superposition of at least two internal waves is needed to achieve symmetry over the ridge crest. One of the consequences of such a superposition is that vertical displacements are minimized over the ridge, as observed.
Over the slope (stations 1 and 2), the observed vertical displacements are of similar amplitude to ξbt. However, in the region of largest ξ (station 1, z > 2000 m), ϕ(ξ) does not match ϕ(ξ)bt, but lags ϕ(u) by ∼90° (see Fig. 2). This suggests that those displacements are associated with a baroclinic wave, not the barotropic forcing. At station 3, the computed ξbt far exceeds the observed vertical displacements. Because the bottom slope changes rapidly at this station (see Fig. 6), this may indicate that our model of constant bottom slope is inappropriate for estimating ξbt near the slope break. This dissimilarity may also indicate that the observed displacements represent the baroclinic response, not barotropic heaving.
c. Model comparisons
We compare our observations with the model results of Merrifield and Holloway (2002). There is qualitative agreement between the observed and modeled velocity and vertical displacement (right-hand panels of Figs. 2 and 3). Both model and observations indicate an approximate linear increase in phase with depth at the ridge crest (stations 5 and 4; Figs. 2a,b and 3a,b). Farther from the ridge, the bifurcation in phase propagation (station 3; Figs. 2c and 3c) is less apparent in the model than the observations. The most notable disagreement occurs for vertical displacement at station 1 for which the observations at 2500 m are roughly 3 times those modeled.
We emphasize that the scale on each phase plot (720°) highlights gross similarities but conceals differences. Closer scrutiny of Figs. 2 and 3 indicates that differences in phase between model and observations may exceed 90°. Nevertheless, most vertical trends (and many of the smaller features) are similar in the model and observations. This is noteworthy, considering that the model represents a simplified deterministic process and our observations represent not much more than a single tidal cycle.
d. Modal content
The barotropic plus first 10 flat-bottom baroclinic modes were computed using the stratification at each station. Semidiurnal harmonic fits were projected onto these modes to determine the dominant vertical scales of u′, υ′, p′, and 〈u′p′〉 (Fig. 4). Over the ridge (stations 4 and 5), the up-ridge and along-ridge velocities (u′, υ′) are strong and dominated by the lowest modes, while the pressure fluctuations are comparatively weak. In contrast, pressure fluctuations over the slope are strong (stations 1 and 2), despite weak velocities. The correlation 〈u′p′〉 is much stronger over the slope than over the ridge. Only on the flank (station 1) does mode 1 dominate the energy flux, with almost all being carried in modes 1 and 2. Lee et al. (2006) find similar modal distributions at all of the 3000-m-isobath stations.
4. Energetics
In the following, we examine the structure of the horizontal and vertical energy flux, as well as the kinetic and potential energy. Each are derived from semidiurnal harmonic fits to perturbation quantities (u′, υ′, ξ, and p′), as discussed in sections 2 and 3. Errors represent standard deviations computed from 100 randomly phased, GM-contaminated synthetic time series with energetics and sampling consistent with the observations (Nash et al. 2005).
a. Energy flux
The cross-ridge and vertical components of the semidiurnal energy flux (〈u′p′〉, 〈w′p′〉) are shown in Fig. 5. Over the ridge (stations 4 and 5) the energy-flux density and its depth average are small. The depth-integrated flux increases monotonically over the slope (stations 1–3), where it forms two distinct beams:
in the upper 500 m, where the energy flux is off ridge and upward;
in the lower 1000 m, where the energy flux is off ridge and downward.
b. Energy
Horizontal kinetic energy HKE, available potential energy APE, and their depth integrals are shown in Fig. 6. The energy density is dominated by HKE over the ridge and by APE over the slope, with the ratio in HKE/APE varying by more than an order of magnitude. This structure contrasts the expected value for a single propagating wave for which HKE/APE is an intrinsic property equal to (ω2 + f 2)/(ω2 − f 2) = 1.34 for M2 at 21.6°N.
Only at the farthest off-ridge station does HKE/APE approach its expected value for a freely propagating wave. This is also the only station where APE and HKE have similar vertical structure.
c. Ridge-top energetics
The near-constant horizontal and vertical phase propagation that spans the ridge crest and slope [i.e., (6) and Fig. 2] suggests that the south-going near-surface beam (Fig. 5) originates on the opposite (north) side of the crest. In addition, the total semidiurnal energy density (APE + HKE, Figs. 4 and 6c) is highest over the ridge crest despite the vanishing energy flux there (Fig. 5).
5. Interpretation
a. Generation
A comparison between (12) and the numerical simulations of Khatiwala (2003) indicates that ΔFE represents an upper bound on internal tide generation. For a ridge width of similar scale to the mode-1 horizontal wavelength, as over Kaena Ridge, Khatiwala’s numerical results agree with (12) to within 10%. However, (12) overpredicts the energy flux in the limits of both longer and shorter wavelength topography. For long-wavelength topography, the horizontal scale of 〈Fbwbt〉 exceeds that of the lowest vertical mode so that much of the transfer remains as part of the barotropic flow. For short-wavelength topography, Baines (1982) overpredicts the induced vertical velocity because of blocking, yielding infinite wbt, 〈Fbwbt〉 and ΔFE in the limit of a knife edge.
Integration of (12) over Kaena Ridge yields ΔFE(x) ∼ 25 kW m−1 per side, as shown in Fig. 7 (top). This exceeds observed 〈u′p′〉 by ∼30%–50%, so is consistent with ΔFE as an upper bound for generation, even after turbulent dissipation is taken into consideration as an energy sink (∫ρϵ dx dz ≃ 1.2 kW m−1 for each side of Kaena Ridge; Klymak et al. 2006).
For Kaena Ridge (Fig. 7), barotropic power Fbwbt is evenly distributed below 300 m—a result of the compensation of the decay of w2bt away from the bottom by the vertical structure of N 2. For a two-dimensional ridge, energy may radiate along one of four internal tide characteristics, given by dz/dx = ±α = ±[(ω2 − f 2)/(N 2 − ω2)]1/2. The pathways for energy flow can be determined by assuming that energy refracts by stratification and reflects from the surface and boundaries according to linear theory (Eriksen 1982). At Kaena Ridge, downward propagation is expected off ridge from the deep supercritical flanks. However, the broad subcritical ridge top leads to the predominance of upward energy propagation in the upper water column (Fig. 7). This is consistent with the observed phase structure of u′ in Figs. 2a and 2b and the intense HKE over the ridge crest. Off ridge, the energy pathways are corroborated by the observed vector energy flux (Fig. 5).
b. The ridge crest as a horizontally standing mode
Barotropic flow over Kaena Ridge simultaneously generates upward ξbt at x ∼ −10 km and downward ξbt at x ∼15 km. Assuming these displacements drive a baroclinic wave spanning the ridge crest, the 25-km separation implies a wave with λx = 50 km. From the dispersion relation, this implies λz = 850 m, similar to that observed in sections 3a and 3b. This suggests that the internal tide’s vertical wavelength over the ridge crest may be set by the width of the ridge (i.e., the separation between generation regions).
Over the ridge crest (stations 5 and 4; Figs. 2a,b), a linear regression of ϕ(u) with depth suggests λz = 1000 m; the average amplitude is 0.16 m s−1.2 We model this velocity structure as a superposition of two linear waves with amplitude uo = 0.08 m s−1 and vertical wavelength m = −2π/1000 m (downward phase propagation) in constant stratification. They have opposite horizontal wavenumbers k = ±αm, so that one propagates north and the other south, and are phased so that the resultant wave field is symmetric about x = 0. The horizontal energy flux density associated with each is 4.2 W m−2, yielding a depth integral of 4.2 kW m−1—about 30% less than that estimated from (10) because only one wavelength is considered here. Neither surface nor bottom boundary conditions can be satisfied by this formulation. However, the resultant cross-ridge standing mode satisfies our above intuition that ξbt is 180° out of phase in the opposing generation regions (separated horizontally by 25 km).
Maxima in APE and 〈w′p′〉 occur at x = (2n − 1)π/(2k); HKE is out of phase with maxima at x = nπ/k (Fig. 8b). As compared with a single propagating wave for which HKE/APE = (ω2 + f 2)/(ω2 − f 2) = 1.34 uniformly, the ratio HKE/APE for this vertically propagating, horizontally standing mode ranges from ∞ over the ridge crest to 0 at x = ±λx/2. Note that the transverse energy flux 〈υ′p′〉tot does not vanish, having extrema at 14-km intervals. This is consistent with the westward component of the observed energy flux at stations 3 and 4 (Fig. 1).
We stress that the standing wave pattern (13) is plausible only where both north- and southbound waves are present (i.e., between the two generation regions—over the ridge crest). In this region, the comparison in Fig. 8 suggests wave interference as a means for generating the observed HKE/APE structure. HKE density is maximum over the ridge crest in both model and observations, with values of 7 and 8–11 J m−3, respectively. APE increases over the slope to ∼6 J m−3 at x = 14 km in both observations and model. The model is not applicable seaward of the generation region (x ∼ −14 km), where most of the wave flux is off ridge (Fig. 5).
Regions of high APE and 〈w′p′〉 coincide with the steep ridge flanks where ξbt and 〈Fbwbt〉 are largest (Fig. 7). In contrast, the superposition produces weak vertical displacements and APE over the relatively flat ridge crest. We hypothesize that it is no coincidence that the dominant vertical wavelength over Kaena Ridge is 1000 m, as this is the required wavelength to create HKE/APE banding that matches the ∼25 km horizontal spacing between generation regions. It therefore appears that the vertical and horizontal wavelengths of the internal tide (λz ∼ 1000 m and λx ∼ 50 km) are tuned to the topography.
6. Conclusions
Phase propagation, energy density, and energy flux 〈u′p′〉 of the semidiurnal internal tide have been examined at Kaena Ridge during a spring tide in October 2000. Over the ridge crest, u′ and p′ are out of phase and yield negligible net depth-integrated energy flux. Over the steep slopes of the ridge flank, the energy flux diverges strongly, forming two beams radiating to the south: one propagates downward near the flank bottom, the other upward near the surface. In these beams, u′ and p′ are in phase, with 90% of the fluctuations contributing to the ∼16 kW m−1 depth-integrated energy flux carried past the 3000-m isobath.
The beam near the flank bottom occurs below a semidiurnal internal wave characteristic emanating from the rim of the ridge crest. According to linear theory, ridge-crest energy cannot penetrate into this region, so the source of this beam must be the nearby steep, supercritical slope. In contrast, phase propagation in the upper 500 m is continuous across the ridge, suggesting that both northern and southern flanks influence the upward-propagating beams. Based on the phase propagation of velocity, the vertical and horizontal wavelengths are 800 m and 44 km, respectively (the corresponding wavelengths based on ξ are 1000 m and 100 km). The source of the discrepancy between λ(u)x and λ(ξ)x has yet to be rationalized but may arise because ξ represents total vertical displacement, not just the baroclinic response.
APE and HKE vary dramatically across the ridge. The highest APE occurs over the steep slopes, whereas the highest HKE occurs over the crest. As a result, the energy ratio HKE/APE varies by more than a factor of 10 across the ridge and only attains the intrinsic value for a freely propagating wave (1.34) at the farthest off-ridge station (Fig. 6).
Because there is no net energy flux over the ridge crest, it is tempting to assume that the wave field emerging from the ridge is generated at the nearest flank. However, the intense ridge-crest HKE indicates that this is not the case and that the wave field has multiple sources. While the three-dimensionality of the ridge adds complexity and produces along-ridge propagation (Rainville and Pinkel 2006), we employ a two-dimensional approach to explain the cross-ridge energy and energy-flux structure. Assuming a symmetric superposition of north- and southbound waves over the ridge crest, each must carry ∼5–6 kW m−1 to account for the modal distribution of energy density (10); inclusion of along-ridge contributions to the energy density would reduce this estimate. Based on linear theory, such a superposition leads to a horizontally standing mode (Nash et al. 2004) with alternating bands of HKE and APE. The spacing between bands is set solely by the wavenumber of the internal tide. We find that the magnitude and expected spacing between bands of maximum HKE and APE (∼14 km) are consistent with that observed.
Barotropic flow over Kaena Ridge simultaneously displaces isopycnals in opposite directions over either ridge flank. The ∼25 km horizontal spacing between generation regions (Fig. 7) sets a horizontal scale of the forcing to generate waves with horizontal wavelength λx ∼50 km. Through the dispersion relation, this implies a vertical wavelength of λz ∼1000 m, as observed.
We thus suggest that the vertical wavelength of baroclinic response is set by the horizontal scale of the feature that generates it. For Kaena Ridge, this wavelength is 1000 m, similar to the water depth. We speculate that Kaena Ridge is a particularly efficient generation site for the semidiurnal internal tide because of this agreement. A much narrower feature may generate a higher wavenumber response and, conversely, a wider feature may generate lower vertical wavenumbers. Since high wavenumbers dissipate more quickly than the low modes (St. Laurent and Nash 2004), the ridge width (which presumably sets the near-field wavenumber content) may be the important parameter in determining the fraction of energy that is dissipated locally as compared with the fraction that is radiated away.
Determination of the barotropic-to-baroclinic conversion rate from steep, large-amplitude topography has had recent theoretical success for knife-edge and top-hat bathymetries (Llewellyn Smith and Young 2003; St. Laurent et al. 2003). For Hawaiian Ridge geometry, these predict energy fluxes of appropriate magnitude and an energy-flux spectrum in agreement with observations (St. Laurent and Nash 2004). A common feature of the solutions is the beamlike character of the baroclinic response, similar to that observed. For example, in the case of the top-hat ridge (St. Laurent et al. 2003), beams emanate from ridge discontinuities along ray paths. These beams form a superposition of mostly upward-propagating energy over the ridge crest, and upward- and downward-propagating waves away from the ridge. Significant baroclinic energy occurs over the ridge top despite zero net energy flux there. However, the vertical wavenumbers generated by top-hat topography depend primarily on the ridge height (δ = ho/H) and only very weakly on the ridge width (L. St. Laurent 2005, personal communication) Such steep-topography models therefore do not predict that the ridge-top vertical wavelength is set by the ridge width.
The power available from the barotropic tide (Baines 1982; 〈Fbwbt〉) represents an upper bound for baroclinic generation. Comparisons with Khatiwala (2003) indicate that it fails in the limits of large- and small-wavelength topography. For Kaena Ridge, it overpredicts the energy flux by 30%–50%.
Klymak et al. (2006) find 1.1–1.3 kW m−1 of turbulent dissipation integrated over the cross-ridge section considered here, amounting to ∼10% of the total flux observed. Klymak et al. (2006) and Levine and Boyd (2006) provide evidence for enhanced mixing over and on the upper flanks of Kaena Ridge, sites with near-critical slopes where linear theory predicts a high wavenumber response and associated wave dissipation.
These observations produce as many questions as answers. For example, what is the fate of baroclinic tidal beams as they reflect from the bottom and highly stratified upper ocean? Does the ridge width set the energy partitioning between high and low wavenumbers? Are high modes simply dissipated, or do they undergo other internal wave transformations? And, most important, does the radiated energy ultimately dissipate in the appropriate locations and with a sufficient mixing efficiency to contribute significantly to global mixing budgets?
Acknowledgments
The authors thank John Dunlap, Art Bartlett, Bob Drever, Eric Boget, Dicky Allison, Alana Althaus, and the captain and crew of the R/V Wecoma for making these measurements possible. Princeton Ocean Model output was kindly provided by Mark Merrifield. Lou St. Laurent shared his code to determine top-hat generation. We thank Luc Rainville and Chris Garrett for their helpful reviews of this manuscript. This research was funded under NSF Grants OCE-98-19537 and OCE-01-36116.
REFERENCES
Alford, M. H., 2003: Energy available for ocean mixing redistributed by long-range propagation of internal waves. Nature, 423 , 159–162.
Althaus, A. M., E. Kunze, and T. B. Sanford, 2003: Internal tide radiation from Mendicino Escarpment. J. Phys. Oceanogr., 33 , 1510–1527.
Baines, P. G., 1974: The generation of internal tides over steep continental slopes. Proc. Roy. Soc. London, 277A , 27–58.
Baines, P. G., 1982: On internal tide generation models. Deep-Sea Res., 29 , 307–338.
Bell, T. H., 1975: Topographically generated internal waves in the open ocean. J. Geophys. Res., 80 , 320–327.
Chereskin, T. K., 1983: Generation of internal waves in Massachusetts Bay. J. Geophys. Res., 88 , 2649–2661.
Egbert, G. D., 1997: Tidal data inversion: Interpolation and inference. Progress in Oceanography, Vol. 40, Pergamon. 81–108.
Egbert, G. D., and R. D. Ray, 2001: Estimates of M2 tidal energy dissipation from TOPEX/Poseidon altimeter data. J. Geophys. Res., 106 , 22475–22502.
Eriksen, C. C., 1982: Observations of internal wave reflection off sloping bottoms. J. Geophys. Res., 87 , 525–538.
Farmer, D. M., and J. D. Smith, 1980: Tidal interaction of stratified flow with a sill in Knight Inlet. Deep-Sea Res., 27A , 239–254.
Hibiya, T., 1988: The generation of internal waves by tidal flow over Stellwagen Bank. J. Geophys. Res., 93 , 533–542.
Holloway, P. E., and M. A. Merrifield, 1999: Internal tide generation by seamounts, ridges and islands. J. Geophys. Res., 104 , 25937–25951.
Holloway, P. E., P. G. Chatwin, and P. Craig, 2001: Internal tide observations from the Australian North West Shelf in summer 1995. J. Phys. Oceanogr., 31 , 1182–1199.
Huthnance, J. M., 1989: Internal tides and waves near the continental shelf edge. Geophys. Astrophys. Fluid Dyn., 48 , 81–105.
Johnston, T. M. S., and M. Merrifield, 2003: Internal tide scattering at seamounts, ridges and islands. J. Geophys. Res., 108 .3180, doi:10.1029/2002JC001528.
Kantha, L. H., and C. C. Tierney, 1997: Global baroclinic tides. Progress in Oceanography, Vol. 40, Pergamon. 163–178.
Khatiwala, S., 2003: Generation of internal tides in an ocean of finite depth: Analytical and numerical calculations. Deep-Sea Res., 50 , 3–21.
Coauthors, 2006: An estimate of tidal energy lost to turbulence at the Hawaiian Ridge. J. Phys. Oceanogr., 36 , 1148–1164.
Kunze, E. L., L. K. Rosenfeld, G. S. Carter, and M. C. Gregg, 2002: Internal waves in Monterey Submarine Canyon. J. Phys. Oceanogr., 32 , 1890–1913.
Lee, C. M., E. Kunze, T. B. Sanford, J. D. Nash, M. A. Merrifield, and P. E. Holloway, 2006: Internal tides and turbulence along the 3000-m isobath of the Hawaiian Ridge. J. Phys. Oceanogr., 36 , 1165–1183.
Levine, M. D., and T. J. Boyd, 2006: Tidally forced internal waves and overturns observed on a slope: Results from HOME. J. Phys. Oceanogr., 36 , 1184–1201.
Llewellyn Smith, S. G., and W. R. Young, 2002: Conversion of the barotropic tide. J. Phys. Oceanogr., 32 , 1554–1566.
Llewellyn Smith, S. G., and W. R. Young, 2003: Tidal conversion at a very steep ridge. J. Fluid Mech., 495 , 175–191.
Loder, J. W., D. Brickman, and E. P. W. Home, 1992: Detailed structure of currents and hydrography on the northern side of Georges Bank. J. Geophys. Res., 97 , 14331–14351.
Merrifield, M. A., and P. E. Holloway, 2002: Model estimates of M2 internal tide energetics at the Hawaiian Ridge. J. Geophys. Res., 107 .3179, doi:10.1029/2001JC000996.
Merrifield, M. A., P. E. Holloway, and T. M. S. Johnston, 2001: Internal tide generation at the Hawaiian Ridge. Geophys. Res. Lett., 28 , 559–562.
Munk, W., and C. Wunsch, 1998: Abyssal recipes II: Energetics of tidal and wind mixing. Deep-Sea Res., 45 , 1977–2010.
Nash, J. D., and J. N. Moum, 2005: River plumes as a source of large amplitude internal waves in the coastal ocean. Nature, 437 .doi:10.1038/nature03936.
Nash, J. D., E. Kunze, J. M. Toole, and R. W. Schmitt, 2004: Internal tide reflection and turbulent mixing on the continental slope. J. Phys. Oceanogr., 34 , 1117–1134.
Nash, J. D., M. H. Alford, and E. Kunze, 2005: Estimating internal wave energy fluxes in the ocean. J. Atmos. Oceanic Technol., 22 , 1551–1570.
Petruncio, E. T., L. K. Rosenfeld, and J. D. Paduan, 1998: Observations of the internal tide in Monterey Canyon. J. Phys. Oceanogr., 28 , 1873–1903.
Pingree, R. D., and A. L. New, 1989: Downward propagation of internal tidal energy into the Bay of Biscay. Deep-Sea Res., 36 , 735–758.
Pingree, R. D., and A. L. New, 1991: Abyssal penetration and bottom reflection of internal tide energy into the Bay of Biscay. J. Phys. Oceanogr., 21 , 28–39.
Rainville, L., and R. Pinkel, 2006: Baroclinic energy flux at the Hawaiian Ridge: Observations from the R/P FLIP. J. Phys. Oceanogr., 36 , 1104–1122.
Ray, R. D., and G. T. Mitchum, 1997: Surface manifestation of internal tides in the deep ocean: Observations from altimetry and island gauges. Progress in Oceanography, Vol. 40, Pergamon. 135–162.
Coauthors, 2003: From tides to mixing along the Hawaiian Ridge. Science, 301 , 355–357.
Sanford, T. B., R. G. Drever, and J. H. Dunlap, 1978: A velocity profiler based on the principles of geomagnetic induction. Deep-Sea Res., 25 , 183–210.
Sanford, T. B., R. G. Drever, and J. H. Dunlap, 1985: An acoustic Doppler and electromagnetic velocity profiler. J. Atmos. Oceanic Technol., 2 , 110–124.
Sanford, T. B., E. A. D’Asaro, E. L. Kunze, J. H. Dunlap, R. G. Drever, M. A. Kennelly, M. D. Prater, and M. S. Horgan, 1993: An XCP user’s guide and reference manual. Applied Physics Laboratory Tech. Rep. APL-UW TR9309, University of Washington, Seattle, WA, 59 pp. [Available from APL, University of Washington, 1013 NE 40th St., Seattle, WA 98105-6698.].
Simmons, H., R. Hallberg, and B. Arbic, 2004: Internal wave generation in a global baroclinic tide model. Deep-Sea Res., 51 , 3043–3068.
St. Laurent, L., and C. Garrett, 2002: The role of internal tides in mixing the deep ocean. J. Phys. Oceanogr., 32 , 2882–2899.
St. Laurent, L., and J. D. Nash, 2004: An examination of the radiative and dissipative properties of deep ocean internal tides. Deep-Sea Res., 51 , 3029–3042.
St. Laurent, L., S. Stringer, C. Garrett, and D. Perrault-Joncas, 2003: The generation of internal tides at abrupt topography. Deep-Sea Res., 50 , 987–1003.
Wunsch, C., 1998: The work done by the wind on the oceanic general circulation. J. Phys. Oceanogr., 28 , 2332–2340.
Bathymetry of Kauai Channel (500-m contour intervals; bathymetry smoothed at 1 km) showing the five XCP stations (1–5) and AVP station (1) used in this study. The depth-integrated semidiurnal baroclinic energy flux computed from these observations is shown for reference.
Citation: Journal of Physical Oceanography 36, 6; 10.1175/JPO2883.1
Time series of upslope velocity fluctuation [u′ = u(z, t) −
Citation: Journal of Physical Oceanography 36, 6; 10.1175/JPO2883.1
As in Fig. 2, but for vertical displacement ξ, semidiurnal phase ϕ(ξ), and amplitude ξo at five XCP stations across Kaena Ridge. Shading in middle panels represents 95% confidence limits on the baroclinic phase. Vertical gray lines in the right-hand panels indicate the phase and amplitude of vertical displacements induced by barotropic flow up a slope [ξbt; Eq. (7)].
Citation: Journal of Physical Oceanography 36, 6; 10.1175/JPO2883.1
Modal decomposition of velocity, pressure, and energy flux as a function of distance across Kaena Ridge. These modal amplitudes represent the projection of semidiurnal fits to u′, υ′, and p′ onto the first 10 flat-bottom vertical modes at each station; 〈u′p′〉 is the product of the amplitudes of u′, p′, and the cosine of the difference in their phases. Because of changes in the water depth, the dimensional wavelengths associated with mode one over the ridge are similar to mode 2 off ridge. Modal structure functions are normalized to unity variance, so modal amplitudes are representative of depth averages not depth integrals. Depth-integrated energy fluxes (the sum of the modal amplitudes of 〈u′p′〉, multiplied by the water depth) are shown in Fig. 5.
Citation: Journal of Physical Oceanography 36, 6; 10.1175/JPO2883.1
Cross-ridge distribution of the upslope (cross ridge) component of (a) energy flux, (b) depth-integrated energy flux with standard error, and (c) coherence between u′ and p′. Black vectors in (a) represent the vertical and horizontal components of the energy flux (〈u′p′〉, 〈w′p′〉); the vertical component has been stretched commensurate with the figure’s vertical exaggeration. Two internal tide ray paths originating from critical topography are indicated with dashed lines. In (c), the absolute value of the depth-averaged energy flux density (dark shading) is compared to the square root of the product of variances
(light shading), where the overbar represents a vertical average; the indicated coherence represents the ratio of these two quantities.
Citation: Journal of Physical Oceanography 36, 6; 10.1175/JPO2883.1
Cross-ridge distribution of (a) available potential energy density (APE, dark shading) and horizontal kinetic energy density (HKE, light shading) in the semidiurnal waveband. (b) The depth-integrated kinetic and potential energies and their ratios (numbers). (c) The associated energy densities. Error bars represent Monte Carlo–derived standard deviations.
Citation: Journal of Physical Oceanography 36, 6; 10.1175/JPO2883.1
(bottom) Distribution of the time-average power associated with barotropic flow over Kaena Ridge (11) using TPX0.5 barotropic currents and 3-km smoothed bathymetry. Arrows indicate pathways for energy flow along internal tide ray paths. Internal tide characteristics emanating from critical bathymetry (dotted lines) are indicated and divide the dynamical regions. (top) An upper bound to the depth-integrated energy flux is calculated by assuming ridge symmetry and integrating (12) from 0 to x.
Citation: Journal of Physical Oceanography 36, 6; 10.1175/JPO2883.1
The structure of an internal tide formed from the superposition of two plane waves: one propagating toward +x, +z and the other toward −x, +z. Each wave has 1000-m vertical wavelength, as observed over the ridge crest (the associated horizontal wavelength λx = 56 km). Average APE (dark shading) and HKE (light shading) are shown by smooth curves in (a); curves are dashed where a two-wave superposition is not valid (|x| > 15). Bars in (a) indicate the observed energy densities shown in Fig. 6c. The resultant current ellipses are shown in (b); vertical exaggeration of current ellipses is consistent with the figure aspect ratio. Arrows represent the distribution of vertical energy flux, and the along-ridge energy flux is indicated with ⊗ (into the page) and ȯ (out of the page). All energetics are independent of depth and periodic across the ridge with period λx/2 except for the cross-ridge energy flux, 〈u′p′〉, which is zero everywhere.
Citation: Journal of Physical Oceanography 36, 6; 10.1175/JPO2883.1
Roughly one-third of the XCP drops during this cruise were affected by pressure-induced thermistor failures over a segment of their depth range. This necessitated some interpolation over bad T data for this analysis.
For stations 4 and 5, λz ∼1000 m is slightly larger than that determined when all stations are used in the regression (6).