## 1. Introduction

The tides are believed to be a significant source of energy to the deep ocean, energizing the mixing between water masses that drives the overturning circulation (Huang 1994; Munk and Wunsch 1998). In the deep ocean, an important mechanism of energizing the mixing is via the creation of internal tides. Two end points of internal tide generation have been considered, generation over rough but not very steep topography (St. Laurent and Garrett 2002) and generation over abrupt topographic features (Llewellyn Smith and Young 2003; St. Laurent et al. 2003). The latter process has been studied in fjords (Klymak and Gregg 2004; Inall et al. 2005) and more recently at Hawaii (Merrifield and Holloway 2002; Rudnick et al. 2003). Interestingly, both fjord studies and the observations at Hawaii indicate that most of the energy lost from the barotropic conversion at abrupt topography radiates away as low-mode internal tides rather than dissipating locally (Klymak and Gregg 2004; Klymak et al. 2006a). Recent numerical modeling supports the inferences from observations (Carter et al. 2008) and offers a simple theoretical explanation (Klymak et al. 2010b) in terms of trapped lee waves (Nakamura et al. 2000; Legg and Klymak 2008; Klymak et al. 2010a).

Because so much of the energy radiates away from abrupt topography, considerable attention has been placed on understanding where that energy eventually dissipates. One hypothesis is that there are dissipative losses when the tides scatter at isolated seamounts (Johnston and Merrifield 2003) or when they impact continental slopes (Nash et al. 2004, 2007; Martini et al. 2010, manuscript submitted to *J. Phys. Oceanogr.*; Kelly et al. 2010). When internal waves impact a slope, the behavior depends on whether the slope is supercritical (*β* < *dh*/*dx*) or subcritical (*β* > *dh*/*dx*) with respect to the incoming waves, where *β* = [(*N*^{2} − *ω*^{2})/(*ω*^{2} > *f*^{2})]^{1/2} is the aspect ratio of energy propagation of internal waves and *dh*/*dx* is the local topographic slope (*N* is the local buoyancy frequency, *ω* is the tidal frequency and *f* is the Coriolis frequency). For subcritical slopes, presumably much of the energy scatters upslope onto the continental shelf, where it likely dissipates in shallow water, whereas for supercritical slopes an undetermined amount of energy will reflect back into the basin. In both cases, there is presently no method for predicting a priori how much energy is dissipated locally, though laboratory experiments indicate that near-critical topographies will have significant nonlinearities (McPhee-Shaw and Kunze 2002).

In this paper, we consider observations of the internal tide impacting the continental slope on the west side of the South China Sea (SCS). The site is directly west of the Luzon Strait, a region of very strong internal tide generation (Niwa and Hibiya 2004; Jan et al. 2007). This leads to large-amplitude nonlinear internal waves propagating across the basin (Zhao et al. 2004; Ramp et al. 2004; Klymak et al. 2006b; Alford et al. 2010) and impacting the continental slope (Lien et al. 2005; Chang et al. 2006; St. Laurent 2008). The exact interplay between these nonlinear waves and the internal tide is still being investigated, but it seems the nonlinear waves are largely derived from the steeper semidiurnal tide rather than the diurnal tide (Farmer et al. 2009; Alford et al. 2010; Buijsman et al. 2010).

This paper examines what happens when this strong internal tide impacts the continental slope, beginning with a description of our mooring and ship-based operations (section 2). The observations indicate that there is significant reflection from the continental slope in addition to substantial local dissipation (section 3). These observations compare favorably to a linear reflection model that indicates that for this stratification and continental shelf depth there should be significant reflection of incoming internal tides (section 4), and we conclude with a discussion of the results (section 5).

## 2. Experiment

We took part in a two-week experiment aboard the *Ocean Researcher I* with observations centered just east of Dong Sha Island (Fig. 1). A line of moorings was deployed along approximately 21°N, with most of the moorings being deployed farther upslope. Here we will concentrate on a mooring deployed near the 1500-m isobath (MP1, Fig. 1b), a ship survey made on the 700-m isobath, and a mooring placed on the Dongsha Plateau in 580 m of water (LR1, Fig. 1b).

The site of mooring MP1 was chosen as part of an array of sensors meant to understand the evolution of solitary waves impacting the slope, the subject of other papers. Here we were interested in the fate of the internal tide and, at sea, decided to occupy a gap in our mooring array. It appears that the east–west profile of the slope at this latitude has a region that is very near critical to the *K*_{1} internal tide (Fig. 2). This was only apparent when shipboard bathymetry became available, and it supports the idea of Cacchione et al. (2002) that internal tides smooth off sharp edges of continental slopes (see section 5). We positioned our ship-based time series directly upslope from this near-critical region with the expectation that the near criticality would drive enhanced mixing at this location.

Comparison of Smith–Sandwell 8.2 bathymetry, used for planning, and the shipboard observed bathymetry (note that newer versions of the bathymetry dataset are substantially improved; see section 5). Also indicated are the ray paths for the dominant semidiurnal (*M*_{2}) and diurnal (*K*_{1}) tidal frequencies. Note that most of the slope is subcritical with respect to *D*_{2} frequencies and supercritical deeper than 800 m with respect to *D*_{1} frequencies. There is a large region of near criticality for *D*_{1} frequencies 650 and 800 m, observed with the ship’s echo sounder.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Comparison of Smith–Sandwell 8.2 bathymetry, used for planning, and the shipboard observed bathymetry (note that newer versions of the bathymetry dataset are substantially improved; see section 5). Also indicated are the ray paths for the dominant semidiurnal (*M*_{2}) and diurnal (*K*_{1}) tidal frequencies. Note that most of the slope is subcritical with respect to *D*_{2} frequencies and supercritical deeper than 800 m with respect to *D*_{1} frequencies. There is a large region of near criticality for *D*_{1} frequencies 650 and 800 m, observed with the ship’s echo sounder.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Comparison of Smith–Sandwell 8.2 bathymetry, used for planning, and the shipboard observed bathymetry (note that newer versions of the bathymetry dataset are substantially improved; see section 5). Also indicated are the ray paths for the dominant semidiurnal (*M*_{2}) and diurnal (*K*_{1}) tidal frequencies. Note that most of the slope is subcritical with respect to *D*_{2} frequencies and supercritical deeper than 800 m with respect to *D*_{1} frequencies. There is a large region of near criticality for *D*_{1} frequencies 650 and 800 m, observed with the ship’s echo sounder.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

The slope mooring, MP1, was located at 20.919°N, 117.895°E. Its main instrument was a McLane Moored Profiler (MMP) that moved between 20- and 1450-m depths 14.2 times a day. The profiler measured temperature, salinity, and pressure with a Falmouth Scientific CTD, and velocity with a Sontek acoustic Doppler velocimeter. The mooring also had an RD Instruments 75-kHz acoustic Doppler profiler nominally moored at 80-m depth, though it experienced substantial blowdown. Data were saved every 80 s from 16-m vertical depth bins.

The plateau mooring, LR1, was located at 21.087°N, 117.347°E in 580 m of water and consisted of a single upward-looking 75-kHz broadband ADCP that was placed on the seafloor, saving data every 90 s from 16-m vertical bins. This instrument had been deployed 4 months previously and was recovered at the end of the MP1 deployment cruise, so there is only an 11-day overlap with the data from MP1.

The shipboard measurements were made with a fast-profiling CTD. This system drops at approximately 3 m s^{−1}, tethered to a Kevlar cable, and transited the 700-m depth every 8.4 min. These observations were supplemented by acoustic Doppler velocities from a 75-kHz ADCP mounted on the ship. This Doppler collected 16-m velocity measurements in narrowband mode to approximately 600 m deep in the water column. The ship transit was made at 2 kt north-northeast along the 700-m isobath.

*T*–

*S*relationship in the SCS is such that temperature dominates the stratification and there are no salinity-compensated intrusions (Fig. 3). Therefore, the overturning analysis was performed using just temperature sensors. Potential temperature profiles were sorted vertically, and the vertical displacements Δ

*z*were used to calculate a dissipation rate,

*N*is the buoyancy frequency calculated from the density profile under the same sorting (Alford and Pinkel 2000; Nash et al. 2007; Klymak et al. 2008). Overturns that did not meet a noise criteria of Δ

_{s}*T*>

*n*

^{−l/2}0.01°C were discarded, where Δ

*T*is the temperature difference across an overturn and

*n*is the number of temperature samples within the overturn; 0.01°C is much greater than the noise level of the thermistors. This screening step made very little difference in the observed dissipation rate because temperature sensors are very precise (as compared to calculating density overturns, which can suffer from salinity spiking) and because the dominant overturns were very large.

*T*–*S* relation at mooring MP1 over the 44-day deployment. Contours are potential density (kg m^{−3}) relative to the sea surface.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

*T*–*S* relation at mooring MP1 over the 44-day deployment. Contours are potential density (kg m^{−3}) relative to the sea surface.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

*T*–*S* relation at mooring MP1 over the 44-day deployment. Contours are potential density (kg m^{−3}) relative to the sea surface.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

## 3. Experimental results

### a. Overview

First, we briefly overview the data collected before proceeding to more quantitative analyses. The offshore mooring, MP1, shows a strong baroclinic tidal signal (Fig. 4a), with a pronounced near-bottom intensification. Isopycnals near the seafloor demonstrate very large vertical excursions, over 200 m during spring tides, and velocities exceed 0.7 m s^{−1}. The strong near-bottom forcing drives large overturning events, some exceeding 200 m in height, allowing us to infer large turbulence dissipations (Fig. 4b). The overturning events demonstrate a clear spring–neap cycle in phase with the barotropic forcing for the region, though perhaps lagged by two days, consistent with propagation from Luzon Strait (as indicated in Fig. 4d).

Observations at mooring MP1. (a) East–west velocity and isopycnals spaced 200 m in the background density profile. Blue is westward (on slope). (b) Dissipation rate estimated via overturns. (c) Depth-integrated dissipation rate, raw and smoothed by 1 day. (d) Depth-integrated baroclinic westward energy flux smoothed by 24 h: raw (black), sum of flux based on bandpassed *K*_{1} and *M*_{2} components (blue), and flux based on *K*_{1} components (cyan). The square of the tidal velocity at the Luzon straits smoothed by 24 h and lagged by two days is shown for comparison.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Observations at mooring MP1. (a) East–west velocity and isopycnals spaced 200 m in the background density profile. Blue is westward (on slope). (b) Dissipation rate estimated via overturns. (c) Depth-integrated dissipation rate, raw and smoothed by 1 day. (d) Depth-integrated baroclinic westward energy flux smoothed by 24 h: raw (black), sum of flux based on bandpassed *K*_{1} and *M*_{2} components (blue), and flux based on *K*_{1} components (cyan). The square of the tidal velocity at the Luzon straits smoothed by 24 h and lagged by two days is shown for comparison.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Observations at mooring MP1. (a) East–west velocity and isopycnals spaced 200 m in the background density profile. Blue is westward (on slope). (b) Dissipation rate estimated via overturns. (c) Depth-integrated dissipation rate, raw and smoothed by 1 day. (d) Depth-integrated baroclinic westward energy flux smoothed by 24 h: raw (black), sum of flux based on bandpassed *K*_{1} and *M*_{2} components (blue), and flux based on *K*_{1} components (cyan). The square of the tidal velocity at the Luzon straits smoothed by 24 h and lagged by two days is shown for comparison.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

The structure of the internal tide at this site is very reminiscent of the observations at Hawaii and the Oregon shelf (Fig. 5). On-slope flow causes the isopycnals to steepen and eventually break, and the appearance of large overturns tends to be phase locked to the tide. There is also a phase change to the tide between 600- and 1000-m depths. Although this is close to the zero crossing of the mode-1 tide in the open basin (though not the local mode 1; Fig. 6), it is also near the depth of the shelfbreak and is perhaps indicative of a beam like structure at this depth (see below). The upward phase propagation near the seafloor indicates downward energy propagation rather than a simple mode-1 response.

Four days from MP1. Contoured isopycnals are from every 200 m in the background profile.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Four days from MP1. Contoured isopycnals are from every 200 m in the background profile.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Four days from MP1. Contoured isopycnals are from every 200 m in the background profile.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

(a) Mean stratification at mooring MP1 and from 120°E, the open-basin 2005 cruise (“deep” from Klymak et al. 2006b). (b) Mode-1 and mode-2 vertical mode for MP1 (gray) and deep CTD (black). The deep-water linear mode speed is 2.6 m s^{−1}, and the slope speed is 2.1 m s^{−1}. The 500-m mode-1 shape is also indicated (thick black line).

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

(a) Mean stratification at mooring MP1 and from 120°E, the open-basin 2005 cruise (“deep” from Klymak et al. 2006b). (b) Mode-1 and mode-2 vertical mode for MP1 (gray) and deep CTD (black). The deep-water linear mode speed is 2.6 m s^{−1}, and the slope speed is 2.1 m s^{−1}. The 500-m mode-1 shape is also indicated (thick black line).

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

(a) Mean stratification at mooring MP1 and from 120°E, the open-basin 2005 cruise (“deep” from Klymak et al. 2006b). (b) Mode-1 and mode-2 vertical mode for MP1 (gray) and deep CTD (black). The deep-water linear mode speed is 2.6 m s^{−1}, and the slope speed is 2.1 m s^{−1}. The 500-m mode-1 shape is also indicated (thick black line).

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

The character of the tide farther up the slope, near the 700-m isobath, is demonstrated in the ship-based time series (Fig. 7). The diurnal tide dominates at this location, with downward propagating phase. Solitary wave events are resolved by the high-frequency sampling at 0300 and 2000 UTC on day 125 of 2007. Strong dissipation is observed near the seafloor, associated with the onshore phase of the near-bottom baroclinic currents rather than the solitary waves.

On-slope velocity and dissipation observed at the shipboard time series at 700-m depth, just upslope of the critical point. This is during the last flood tide in Fig. 4 (6 May is yearday 125). Contours are isotherms (°C).

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

On-slope velocity and dissipation observed at the shipboard time series at 700-m depth, just upslope of the critical point. This is during the last flood tide in Fig. 4 (6 May is yearday 125). Contours are isotherms (°C).

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

On-slope velocity and dissipation observed at the shipboard time series at 700-m depth, just upslope of the critical point. This is during the last flood tide in Fig. 4 (6 May is yearday 125). Contours are isotherms (°C).

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

### b. Decomposition of signals

Our goal is to understand energy fluxes over the shelfbreak region. To calculate these it is helpful to decompose the signals both in frequency and depth. We use the usual methods for these decompositions; first vertical means are removed, so that *u*′(*z*, *t*) = *u*(*z*, *t*) − 〈*u*(*z*, *t*)〉* _{z}*. The perturbation signals at both moorings are very much concentrated at the diurnal and semidiurnal bands (

*D*

_{1}and

*D*

_{2}; Figs. 8, 9), so these are isolated with second-order Butterworth filters with cutoff frequencies for Δ

*f*(

*D*

_{1}) = [0.71–1.4] cycles per day and Δ

*f*(

*D*

_{2}) = [1.43–3.3] cycles per day (indicated in Fig. 8). The

*D*

_{2}bandpass picks up some of the 2

*D*

_{1}harmonic, so some of the nonlinearity of the tide is included in the filtering.

Spectra of mean velocity between 144 and 224 m on MP1; comparison between the Doppler and the MMP. The main peaks are well captured by the coarser sampling of the MMP, but the high-frequency nonlinear harmonics are underestimated. Gray shading indicates the bandpass frequencies used in the paper.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Spectra of mean velocity between 144 and 224 m on MP1; comparison between the Doppler and the MMP. The main peaks are well captured by the coarser sampling of the MMP, but the high-frequency nonlinear harmonics are underestimated. Gray shading indicates the bandpass frequencies used in the paper.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Spectra of mean velocity between 144 and 224 m on MP1; comparison between the Doppler and the MMP. The main peaks are well captured by the coarser sampling of the MMP, but the high-frequency nonlinear harmonics are underestimated. Gray shading indicates the bandpass frequencies used in the paper.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Kinetic energy spectra at the plateau mooring LR1. Component vertical-mode spectra are shown as gray lines.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Kinetic energy spectra at the plateau mooring LR1. Component vertical-mode spectra are shown as gray lines.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Kinetic energy spectra at the plateau mooring LR1. Component vertical-mode spectra are shown as gray lines.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

It is also worth noting that the large sampling interval of the MMP mooring MP1 means that it misses some high-frequency energy (Figs. 8, 10). Some of this high-frequency energy is contained in solitary waves that pass through the region at a semidiurnal interval (Ramp et al. 2004; Alford et al. 2010). The sampling from the MMP does not capture these signals very clearly. Regardless, the kinetic energy observed by each method is not very different when averaged over a day (Fig. 10, thin lines). If we assume energy flux is proportional to energy density, as is the case for linearly propagating waves, the MMP measurements are capturing almost 80% of the total flux.

Comparison of the velocity signal from the MMP and Doppler. The mean velocity between 144 and 224 m was squared as a rough proxy for energy flux. One-day running means of both signals are shown as thin lines. Based on this, the coarser sampling of the MMP may underestimate energy flux by approximately 20%.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Comparison of the velocity signal from the MMP and Doppler. The mean velocity between 144 and 224 m was squared as a rough proxy for energy flux. One-day running means of both signals are shown as thin lines. Based on this, the coarser sampling of the MMP may underestimate energy flux by approximately 20%.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Comparison of the velocity signal from the MMP and Doppler. The mean velocity between 144 and 224 m was squared as a rough proxy for energy flux. One-day running means of both signals are shown as thin lines. Based on this, the coarser sampling of the MMP may underestimate energy flux by approximately 20%.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

*d*Ψ/

*d*

*z*(

*z*= 0) =

*d*Ψ/

*dz*(

*z*= −

*H*) = 0 and therefore has orthogonal solutions that decompose the flow into discrete vertical modes (Kundu 1990). These modes will be modified in the presence of mean shear, a complication we ignore here because mean flows are relatively weak.

*c*are the group and phase speed for each mode. In a rotating ocean, the horizontally propagating signals are Poincaré waves and the phase and group speeds are

_{e}*c*and vertical structure Ψ(

_{e}*z*) are determined numerically from observed profiles of

*N*(

*z*) and projected onto the perturbation signals,

There is a relatively large caveat in decomposing flows into vertical modes over sloping topography, because technically the decomposition is only valid on a flat-bottom ocean. The problem can be appreciated by considering the shape of vertical modes in the deep basin (Fig. 6, black lines), at MP1 (gray lines), and on the plateau (thick gray line). Clearly, signals identified as mode 1 in the deep basin are modified in space as they propagate onto the plateau. Of particular note is the comparison of mode 2 in the basin and mode 1 on the plateau; if the transition to the plateau is abrupt enough, a mode-2 wave in the basin may strongly excite a mode-1 response on the plateau. As we see below, this scattering of the incoming energy from deep modes into shallow modes can significantly complicate the interpretation of signals at any of our observation sites.

### c. Net fluxes

*u*and

*p*to get

*u*′(

*z*,

*t*) and

*p*′(

*z*,

*t*), where

*D*

_{2}tidal periods for the results shown here. We decompose the perturbations into frequency windows and vertical modes to isolate the various components of the energy fluxes (Fig. 11).

Decomposition of tidal fluxes at MP1 (continental slope). Positive fluxes are on slope (i.e., propagating from east to west): (a) all frequencies decomposed into vertical modes, (b) bandpassed at *D*_{1} [gray shading is total flux from (a)], and (c) bandpassed at *D*_{2}.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Decomposition of tidal fluxes at MP1 (continental slope). Positive fluxes are on slope (i.e., propagating from east to west): (a) all frequencies decomposed into vertical modes, (b) bandpassed at *D*_{1} [gray shading is total flux from (a)], and (c) bandpassed at *D*_{2}.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Decomposition of tidal fluxes at MP1 (continental slope). Positive fluxes are on slope (i.e., propagating from east to west): (a) all frequencies decomposed into vertical modes, (b) bandpassed at *D*_{1} [gray shading is total flux from (a)], and (c) bandpassed at *D*_{2}.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

The internal wave fluxes impacting the continental slope at MP1 are relatively complicated (Fig. 11). In general, the fluxes are on slope; mode 1 dominates, but mode 2 is also important, particularly during diurnal spring tides when the two fluxes are approximately equivalent. Most of the energy is contained in the diurnal and semidiurnal tides (dashed line in Fig. 11c and Table 1), with diurnal mode 1 dominating. The semidiurnal tide is almost all mode 1 and is relatively constant through the mooring deployment, with slight peaks out of phase with the diurnal spring tides (Fig. 11c, solid line). The decomposition is summarized in Table 1.

MP1 and LR1 modal and frequency mean energy fluxes for a 28-day period; positive is shoreward. Energy flux for LR1 is estimated from energy density and is an upper bound.

The energy flux farther onshore on Dongsha Plateau can be estimated from the Doppler mooring LR1 at 580-m depth (Figs. 12, 9). The energy here is mostly mode 1, though the diurnal component has the energy inverted, with more mode 3 than modes 2 or 1 (Fig. 9). Much of the energy is either diurnal or semidiurnal; however, higher harmonics also play a large role, indicative of the nonlinear effects on the plateau. To calculate an energy flux, we assume that the energy is all propagating shoreward and therefore approximate the tidal-period mean flux as *F _{n}* =

*c*=

_{gn}E_{n}*c*Ke2

_{e}*ω*(

*ω*

^{2}−

*f*

^{2})

^{1/2}/(

*ω*

^{2}+

*f*

^{2}) for each mode

*n*and bandpassed frequency

*ω*, where Ke is the local kinetic energy evaluated from the Doppler and

*c*is the linear eigenspeed (Chang et al. 2006).

_{e}Energy flux estimated from standing energy at onshore mooring LR1. Note that the deployment only overlaps with MP1 from day 115 to when it was pulled from the water on day 128.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Energy flux estimated from standing energy at onshore mooring LR1. Note that the deployment only overlaps with MP1 from day 115 to when it was pulled from the water on day 128.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Energy flux estimated from standing energy at onshore mooring LR1. Note that the deployment only overlaps with MP1 from day 115 to when it was pulled from the water on day 128.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

As noted below, this will overpredict the net energy flux propagating shoreward because reflected waves will reduce the net flux but enhance the energy density. In addition, any noninternal wave energy will not be filtered out in the same way that it is if we can use both pressure and velocity perturbations. Therefore, this estimate is a crude upper bound. If we do not do any frequency filtering, we still get very large values averaging almost 20 kW m^{−1}, which is not very realistic and is not presented. If we isolate the diurnal and semidiurnal components, we get values that are higher than observed at the deeper mooring (Fig. 12, Table 1), particularly for *D*_{2}, so we consider this calculation an upper bound on the possible fluxes.

### d. Evidence for reflected energy at the continental slope

*F*to the energy density

*E*for horizontal plane mode-1 and mode-2 waves. For a horizontally progressive wave,

*F*=

*c*if averaged over a number of periods. If there is a horizontally standing wave (100% reflection), the net flux will be zero,

_{g}E*F*=

*F*−

_{o}*F*, but there will be standing energy. In the absence of rotation, the standing energy is simply

_{o}*E*= 2

*F*/

_{o}*c*. Rotation adds an along-flux dependence to the energy density (Nash et al. 2004; Alford and Zhao 2007; Martini et al. 2007) that requires some care to derive. If

_{g}*F*is the incoming flux and

_{o}*F*is the reflected flux, then the energy density is

_{r}*x*is the distance from the reflection point and

*k*is the horizontal wavelength of the internal wave. If

*F*= −

_{r}*αF*, then the observed energy flux is

_{o}*F*=

*F*(1 −

_{o}*α*) and we can write

*α*. There is some uncertainty inherent in choosing exactly how far MP1 is from the loosely defined reflection point. Here we use

*x*= 30 km, and note that a ±10 km change in this has an approximately 20% error to the reflection estimate.

This estimate is relatively crude, because *E* will include any noise in the measurements whereas *F* tends to filter out uncorrelated noise between the density and velocity measurements. However, the vertical-mode fits and the frequency filtering serve to reduce noise. It also should be borne in mind that this analysis does not distinguish between reflected energy and energy that is radiated offshore from barotropic conversion at the topographic break; we estimate barotropic conversion below.

These caveats in mind, it appears that there is considerable reflection of the incoming internal tide energy at this site (Fig. 13a), with the flux-weighted average ⅓ of the incoming average for the mode-1 diurnal signal (Fig. 13c, dashed line). In contrast, the semidiurnal energy flux has weak reflection (Fig. 13i), consistent with the subcriticality of the slope to the semidiurnal tidal frequencies. Interestingly, the mode-2 diurnal reflection has only a weak reflection (Fig. 13f), a surprising feature that we reconcile below.

Energy density compared to energy flux at MP1: (a)–(c) mode-1 bandpassed at *D*_{1} energy fluxes, (d)–(f) mode-2 *D*_{1} energy, and (g)–(i) mode-1 *D*_{2} energy flux. Shown are (top) the fluxes and the standing energy, (middle) the ratio of these two (for a purely progressive wave, these two quantities should be equal) and (bottom) the reflection coefficient [Eq. (9); solid line] and the flux-weighted mean of that coefficient (dashed line). There is significantly more standing mode-1 *D*_{1} energy than energy flux, indicating strong reflection [(c)], but mode 2 has almost no reflection [(f)]. The energy in the *D*_{2} band also appears to be largely progressive [(i)].

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Energy density compared to energy flux at MP1: (a)–(c) mode-1 bandpassed at *D*_{1} energy fluxes, (d)–(f) mode-2 *D*_{1} energy, and (g)–(i) mode-1 *D*_{2} energy flux. Shown are (top) the fluxes and the standing energy, (middle) the ratio of these two (for a purely progressive wave, these two quantities should be equal) and (bottom) the reflection coefficient [Eq. (9); solid line] and the flux-weighted mean of that coefficient (dashed line). There is significantly more standing mode-1 *D*_{1} energy than energy flux, indicating strong reflection [(c)], but mode 2 has almost no reflection [(f)]. The energy in the *D*_{2} band also appears to be largely progressive [(i)].

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Energy density compared to energy flux at MP1: (a)–(c) mode-1 bandpassed at *D*_{1} energy fluxes, (d)–(f) mode-2 *D*_{1} energy, and (g)–(i) mode-1 *D*_{2} energy flux. Shown are (top) the fluxes and the standing energy, (middle) the ratio of these two (for a purely progressive wave, these two quantities should be equal) and (bottom) the reflection coefficient [Eq. (9); solid line] and the flux-weighted mean of that coefficient (dashed line). There is significantly more standing mode-1 *D*_{1} energy than energy flux, indicating strong reflection [(c)], but mode 2 has almost no reflection [(f)]. The energy in the *D*_{2} band also appears to be largely progressive [(i)].

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

### e. Dissipation on the continental slope

There are very large overturns observed at both MP1 (Fig. 4), and during the shipboard time series (Fig. 7). Overturns of greater than 100 m thick were observed at both sites, with temperature differences exceeding 0.1°C. These overturns are concentrated near the seafloor, and phase locked to the tide, reminiscent of similar features seen at Hawaii (Levine and Boyd 2006; Aucan et al. 2006; Klymak et al. 2008) and on the Oregon Slope (Nash et al. 2007). The strongest dissipation at both locations is associated with the on-slope flow, where the oncoming tide forms a nonlinear shock near the seafloor. Sometimes there are overturning phenomena concentrated at the seafloor associated with off-slope flow (i.e., day 141.5, Fig. 5b).

The time series of the depth-integrated dissipation at MP1 reveals a pronounced spring–neap cycle (Fig. 4d). This cycle correlates with the incoming tidal *D*_{1} flux consistent with a power law *D*_{2} tidal flux (and is in fact slightly inversely scaled, no doubt because the *D*_{2} and *D*_{1} spring–neap cycles are out of phase during these observations) and only scales well with the total tidal flux when *D*_{1} dominates (Figs. 14b,c). This implicates the *D*_{1} internal tide in the local mixing on the continental slope. The scaling *F*^{3/2} was chosen to be consistent with the observations (Levine and Boyd 2006) and theory (Klymak et al. 2010b) that indicates dissipations scales with the barotropic tidal forcing as *F*, which will tend to scale as

Scaling of the 1-day-averaged dissipation with (a) the *D*_{1} energy flux, (b) the *D*_{2} energy flux, and (c) the total energy flux.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Scaling of the 1-day-averaged dissipation with (a) the *D*_{1} energy flux, (b) the *D*_{2} energy flux, and (c) the total energy flux.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Scaling of the 1-day-averaged dissipation with (a) the *D*_{1} energy flux, (b) the *D*_{2} energy flux, and (c) the total energy flux.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

The vertically integrated average dissipation rate at the shipboard occupation at 700 m was 0.105 W m^{−2}, which is about twice the dissipation rate on yearday 125 at MP1 (Fig. 4d). The stratification is such that *N*^{2} is 4 times higher at the shallower depth, so the stratification-scaled energy sink is a factor of 2 higher at the deeper site. This indicates that a dissipation scheme based on a simple halo of stratification-scaled dissipation will not be adequate to describe dissipation at this site.

The mean vertically averaged dissipation at MP1 for the spring–neap cycle from day 125 to day 139 is 7 × 10^{−2} W m^{−2}. Given MP1’s high dissipation signal, we suspect this represents an upper bound for the slope region. If we multiply by the width of the slope, approximately 40 km, we get 3 kW m^{−1} of local dissipation. It is certainly a large dissipation and on the order of the energy fluxes impacting the slope. Without a more general turbulence survey, we are not able to constrain this value further.

### f. Generation at the continental slope

We have no method for separating the local generation at the continental slope from the signals generated from the internal tide. However, we can estimate the order of magnitude of the generation from a simple linear step model (St. Laurent et al. 2003). A similar model for reflection is discussed below, including the details on dealing with stratification, so we do not present the derivation here. Using barotropic *D*_{1} velocities from altimetry (Egbert and Erofeeva 2002) and using a water depth of 2800 m, a ridge height of 700 m, and the observed basin stratification, we get an average seaward generation of approximately 1 kW m^{−1} (Fig. 15). Landward, the predicted generation is only 1/10th of this. The local *D*_{2} tides are very weak (less than 0.01 m s^{−1}), and the local generation is very weak.

(a) Local east–west tide in 2600 m of water. (b) Seaward tidal flux generation estimated from a linear step-topography model of the continental slope.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

(a) Local east–west tide in 2600 m of water. (b) Seaward tidal flux generation estimated from a linear step-topography model of the continental slope.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

(a) Local east–west tide in 2600 m of water. (b) Seaward tidal flux generation estimated from a linear step-topography model of the continental slope.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

## 4. Linear theory for reflection

The observations above indicate that there is considerably more internal tide energy at the MP1 mooring than a progressive wave interpretation of *F* = *c _{g}E* would indicate, and, if we were to take the observations at face value, perhaps 33% of the incoming mode-1

*D*

_{1}energy reflects back into the larger basin and 67% transmits up the slope (Fig. 13). Conversely, very little mode-2 energy is reflected. Here we investigate if this proportion is reasonable by considering a simple linear model of tidal reflection.

To solve this problem we follow the geometric approach suggested by St. Laurent et al. (2003) for baroclinic-to-barotropic tidal conversion.^{1} For supercritical topography, the continental slope is approximated as a step change from the basin depth *H* and to a shelf depth (or plateau depth in our case) *h _{s}* (Fig. 16). The flow is decomposed into an incident wave with velocities

*u*, a reflected wave

_{i}*u*, and a transmitted wave

_{r}*u*. At the shelf break, the velocities match so that

_{t}*u*=

_{t}*u*+

_{i}*u*for

_{r}*z*> −

*h*, and 0 =

*u*+

_{i}*u*for −

_{r}*H*<

*z*< −

*h*, and

*w*=

_{t}*w*+

_{i}*w*. Following St. Laurent et al. (2003), we make a Wentzel–Kramers–Brillouin (WKB) coordinate transform to a constant-stratification ocean. The coordinate transform from real to WKB-stretched coordinates is normalized so that

_{r}Velocity snapshots of linear model for (left) mode-1 and (right) mode-2 incoming waves, with *F*_{1} = 2*F*_{2}. Depth and velocity are WKB scaled,

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Velocity snapshots of linear model for (left) mode-1 and (right) mode-2 incoming waves, with *F*_{1} = 2*F*_{2}. Depth and velocity are WKB scaled,

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Velocity snapshots of linear model for (left) mode-1 and (right) mode-2 incoming waves, with *F*_{1} = 2*F*_{2}. Depth and velocity are WKB scaled,

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

*k*and

_{n}*a*,

_{m}*b*,

_{n}*A*and

_{mn}*B*are given by St. Laurent et al. (2003),

_{nl}*I*is the identity matrix (this notation departs slightly from St. Laurent et al. 2003, who use

_{ml}*δ*= 1 −

*γ*).

### a. Application to South China Sea

If the incoming wave field is largely composed of a single mode, the response of the Dongsha Plateau is relatively straight forward. The fraction of energy reflected back into the basin is a function of the depth of the continental shelf compared to the vertical-mode shapes. The response for mode-1 and mode-2 incoming waves shows that a substantial amount of energy can be reflected back into the basin (Fig. 16), with a standing pattern in the shape of beams that have the vertical extent of the plateau depth. This results in a largely red mode-number spectrum of the reflected energy, analogous to the spectra for barotropic conversion (St. Laurent et al. 2003). The fraction of energy reflected depends on the depth of the plateau, with more reflection for shallower shelves (Fig. 17) and a slight structure in the mode-2 response because of the midwater maximum in the velocity structure. For a mode-1-only incoming flux, a very large fraction of the energy is reflected for a WKB shelf of *F _{r}* ≈ 0.75

*F*.

_{o}Fraction of energy reflected from continental slope as a function of shelf depth (in WKB coordinates) for monochromatic incoming waves.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Fraction of energy reflected from continental slope as a function of shelf depth (in WKB coordinates) for monochromatic incoming waves.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Fraction of energy reflected from continental slope as a function of shelf depth (in WKB coordinates) for monochromatic incoming waves.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

The actual situation in SCS is more complicated, with both mode-1 and mode-2 energy impacting Dongsha Plateau (Fig. 13). In this case, the phase of the incoming waves relative to each other strongly impacts the response. If the waves are in phase at the plateau break such that modes 1 and 2 are positively reinforced at the sea surface, there is very strong transmission onto the plateau and relatively weak reflection (Fig. 18, left). If they are out of phase, strong incoming velocities are felt at the bottom and most of the energy is reflected, with almost none transmitted onto the plateau (Fig. 18, right).

Velocity snapshots of a linear model for mode-1 and mode-2 incoming waves, with *F*_{1} = 2*F*_{2}: (left) waves in phase (=0°) at the sea surface at *x* = 0 and (right) waves out of phase at the sea surface (phase = 180°). Depth and velocity are WKB scaled. Energy fluxes are indicated with text inside the domains.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Velocity snapshots of a linear model for mode-1 and mode-2 incoming waves, with *F*_{1} = 2*F*_{2}: (left) waves in phase (=0°) at the sea surface at *x* = 0 and (right) waves out of phase at the sea surface (phase = 180°). Depth and velocity are WKB scaled. Energy fluxes are indicated with text inside the domains.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Velocity snapshots of a linear model for mode-1 and mode-2 incoming waves, with *F*_{1} = 2*F*_{2}: (left) waves in phase (=0°) at the sea surface at *x* = 0 and (right) waves out of phase at the sea surface (phase = 180°). Depth and velocity are WKB scaled. Energy fluxes are indicated with text inside the domains.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

It is noteworthy that a linear model would have such an asymmetric response. What happens is that the response to an incoming wave is phase locked with the forcing, either perfectly in phase or 180° out of phase. If the mode-1 and mode-2 incident waves are out of phase, their reflections will also be out of phase, and this can lead to positive reinforcement or cancellation of the modes in the response. The waves are linear, but the energy flux is inherently nonlinear.

There is a continuum between these two extremes depending on the exact phase (and amplitude) of the incoming waves (Fig. 19). Interestingly, the amount of mode-2 reflected energy is very low if the waves are nearly in phase (i.e., phase < 50° or > 310°), because the mode-2 wave reflected from the incident mode-1 wave has the opposite phase as the one reflected from the incident mode-2 wave. This compares favorably with the observation above that there is little energy reflected from mode 2 at the MP1 mooring (Fig. 13, middle), indicating that the incoming mode-1 and mode-2 waves may be nearly in phase at the topographic break.

Dependence of fraction of energy reflected on phase between mode-1 and mode-2 waves. The mode-2 wave has half the flux of the mode-1 wave as pictured in Fig. 18. (a) Total fraction reflected as a fraction of energy incident. Thick line is total, and thin lines are cumulative fractions in mode 1 and in modes 1 and 2. (b) Fraction of mode 1 and mode 2 reflected compared to the incident mode-1 and mode-2 energy. The vertical dashed lines are the approximate phase we think is valid at Dongsha Plateau.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Dependence of fraction of energy reflected on phase between mode-1 and mode-2 waves. The mode-2 wave has half the flux of the mode-1 wave as pictured in Fig. 18. (a) Total fraction reflected as a fraction of energy incident. Thick line is total, and thin lines are cumulative fractions in mode 1 and in modes 1 and 2. (b) Fraction of mode 1 and mode 2 reflected compared to the incident mode-1 and mode-2 energy. The vertical dashed lines are the approximate phase we think is valid at Dongsha Plateau.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Dependence of fraction of energy reflected on phase between mode-1 and mode-2 waves. The mode-2 wave has half the flux of the mode-1 wave as pictured in Fig. 18. (a) Total fraction reflected as a fraction of energy incident. Thick line is total, and thin lines are cumulative fractions in mode 1 and in modes 1 and 2. (b) Fraction of mode 1 and mode 2 reflected compared to the incident mode-1 and mode-2 energy. The vertical dashed lines are the approximate phase we think is valid at Dongsha Plateau.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

*x*= 32 km offshore of the reflection point (Fig. 2),

*ϕ*≈ 320°, though slight variations in the distance of the slope break can change this phase substantially (Fig. 21).

Coherence and phase difference between mode-1 and mode-2 east–west velocity at MP1. The diurnal tide (gray box) has a phase difference of almost exactly 180°.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Coherence and phase difference between mode-1 and mode-2 east–west velocity at MP1. The diurnal tide (gray box) has a phase difference of almost exactly 180°.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Coherence and phase difference between mode-1 and mode-2 east–west velocity at MP1. The diurnal tide (gray box) has a phase difference of almost exactly 180°.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Phase of mode-1 and mode-2 velocity signals at *x* = 32 km in the linear model as a function of their phase at the topographic break of the plateau. Thin lines are for *x* = 30 and 34 km. The observed phase at MP1 (*x* ≈ 30 km) is approximately 180° (Fig. 20), so the phase at the topographic break is between 270° and 350°.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Phase of mode-1 and mode-2 velocity signals at *x* = 32 km in the linear model as a function of their phase at the topographic break of the plateau. Thin lines are for *x* = 30 and 34 km. The observed phase at MP1 (*x* ≈ 30 km) is approximately 180° (Fig. 20), so the phase at the topographic break is between 270° and 350°.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Phase of mode-1 and mode-2 velocity signals at *x* = 32 km in the linear model as a function of their phase at the topographic break of the plateau. Thin lines are for *x* = 30 and 34 km. The observed phase at MP1 (*x* ≈ 30 km) is approximately 180° (Fig. 20), so the phase at the topographic break is between 270° and 350°.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

This phase difference between the mode-1 and mode-2 signals at MP1 and the results of the linear model therefore indicate that the fraction of mode-1 energy reflected would be 40% (Fig. 19b) and the fraction of mode-2 energy reflected would be 5%. This is in rough agreement with the observations (Fig. 13), though the predicted reflection for both modes is higher than observed. This fraction is close enough to justify the iterative approach to determining the phase.

If the incoming mode-1 and mode-2 internal tides are 320° out of phase at the reflection point, then we can trace eastward in the basin to find possible points of generation. At a single knife-edge ridge, mode 1 and mode 2 are either in phase or 180° out of phase; for ridges that are the height of those in the Luzon Strait, we would expect the modes to be out of phase. Tracing back (Fig. 22), we find only three locations where the mode-1 or mode-2 internal tides are in or out of phase. One location just west of 119°E is unlikely to be the source of a strong internal tide, leaving the shorter ridge at 120.6°E where the internal tides are predicted to be in phase and 20 km to the east of the taller ridge at 122°E. This second potential generation point is the correct phase to indicate simple knife-edge generation at this location. However, the generation point does not lie exactly at the ridge crest, perhaps because this estimate is very crude; perhaps because of three dimensionality; or perhaps because the reflections, both from the continental slope and from the smaller ridge, complicate the phasing of the internal tide generation. Finally, mode 2 has the potential to be refracted more than mode 1 (Rainville and Pinkel 2006), perhaps changing the phase relation by the time the waves reach the continental slope. The phase response of the complete system will be the subject of further research.

(top) Topography into Luzon Strait. Idealized step topography is indicated as a dashed line. Gray vertical lines are the topographic break and possible generation points. (bottom) Locations where the mode-1 and mode-2 *K*_{1} tides are calculated to be in phase (thin and solid lines intersect) or exactly out of phase (dashed and solid lines intersect). Phases start at the continental slope (117.6°E) and are traced backward to the east.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

(top) Topography into Luzon Strait. Idealized step topography is indicated as a dashed line. Gray vertical lines are the topographic break and possible generation points. (bottom) Locations where the mode-1 and mode-2 *K*_{1} tides are calculated to be in phase (thin and solid lines intersect) or exactly out of phase (dashed and solid lines intersect). Phases start at the continental slope (117.6°E) and are traced backward to the east.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

(top) Topography into Luzon Strait. Idealized step topography is indicated as a dashed line. Gray vertical lines are the topographic break and possible generation points. (bottom) Locations where the mode-1 and mode-2 *K*_{1} tides are calculated to be in phase (thin and solid lines intersect) or exactly out of phase (dashed and solid lines intersect). Phases start at the continental slope (117.6°E) and are traced backward to the east.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

The linear theory is valid for supercritical topography (*β* > *dh*/*dx*). This approximation appears to be valid for knife-edge generation problems (Di Lorenzo et al. 2006; Klymak et al. 2010b), so we anticipate that it is valid for the scattering problem, though verification will need to wait for appropriate numerical studies. Most of the slope is supercritical to the diurnal tides, but it is subcritical to the semidiurnal tides (Fig. 2). For the *D*_{2} tides, we expect most of the energy to propagate onshore, again as indicated by the observations (Fig. 13).

## 5. Discussion

### a. Slope criticality

There is a near-critical region in the continental slope with respect to the *D*_{1} internal tide, which prompted the location of the ship time series. It is interesting to consider how much of the slope is near critical by comparing the bottom slope to the slope of internal waves of either *D*_{1} or *D*_{2} frequency (Figs. 23a,b). Doing this with the readily available Smith–Sandwell database (Smith and Sandwell 1997) yields poor results, but a high-resolution database from National Taiwan University indicates that there are wide regions of the slope that are near critical to *D*_{1} internal waves. There are narrow regions of near criticality for *D*_{2} waves, but they are not as prevalent as for the dominant *D*_{1} tide.

(a) Slope criticality for *D*_{1} tides. (b) Slope criticality for *D*_{2} tides. Bathymetry is a 1.8-km-resolution product from National Taiwan University. Older versions of Smith–Sandwell (i.e., 8.2 and previous) do not show as coherent near-critical regions.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

(a) Slope criticality for *D*_{1} tides. (b) Slope criticality for *D*_{2} tides. Bathymetry is a 1.8-km-resolution product from National Taiwan University. Older versions of Smith–Sandwell (i.e., 8.2 and previous) do not show as coherent near-critical regions.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

(a) Slope criticality for *D*_{1} tides. (b) Slope criticality for *D*_{2} tides. Bathymetry is a 1.8-km-resolution product from National Taiwan University. Older versions of Smith–Sandwell (i.e., 8.2 and previous) do not show as coherent near-critical regions.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

The hypothesis put forward by Cacchione et al. (2002) is that the shelf breaks of continental slopes are eroded by tidal flows over them until they reach the critical condition. For the slope we are considering, there are diurnal internal tides impacting the slope from the east-southeast, so the hypothesis would be that the slope facing in these directions would be more near critical than slopes facing in other directions. We test this by sweeping through the wave directions, and the normalized distribution of slope criticality (Fig. 24) indicates that the southeast is more likely to have near-critical slopes than slopes facing in the other directions. The histograms are diffuse and this observation could have geological reasons, but the observation is indicative.

Probability density functions for criticality for waves coming from a spread of directions. The slopes with the largest near-critical areas are found for waves coming from the southeast, broadly consistent with an internal tide source near the Luzon Strait.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Probability density functions for criticality for waves coming from a spread of directions. The slopes with the largest near-critical areas are found for waves coming from the southeast, broadly consistent with an internal tide source near the Luzon Strait.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Probability density functions for criticality for waves coming from a spread of directions. The slopes with the largest near-critical areas are found for waves coming from the southeast, broadly consistent with an internal tide source near the Luzon Strait.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

There is much less semidiurnal energy on this continental slope, partly because the semidiurnal internal tides generated at Luzon Strait appear to disintegrate into solitons (Farmer et al. 2009; Alford et al. 2010). Although the solitons have considerable energy and energy flux, it is largely at higher frequencies and there is little opportunity for interactions with the seafloor that might cause localized erosion at the continental break. Certainly, in the shipboard observations, there was little evidence of turbulence near the seafloor associated with the solitons (Fig. 7).

The continental slope of the South China Sea is a passive margin and composed of a very thick layer of sediments (>1 km; Nissen et al. 1995), so the mechanism proposed by Cacchione et al. (2002) is feasible insofar as the internal-wave-induced stress does not need to act on bedrock. However, being a passive margin, there are geological reasons that the slope is more gentle here than at an active margin. In their geological history, the islands of Taiwan and Luzon moved to the east-southeast from this continental slope, so the fact that the gentler slopes are in the direction of rifting is perhaps coincidental. However, Luzon Strait itself is more near critical to *D*_{1} tides (Fig. 25), and this is the “active” margin for this basin. The *D*_{1} criticality for a basin that is very resonant to the *D*_{1} tides and has a very strong *D*_{1} internal tide is indicative that the Cacchione et al. (2002) hypothesis is applicable here.

Criticality in Luzon Strait. Again, the strait is more near critical to *D*_{1} than *D*_{2} tides.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Criticality in Luzon Strait. Again, the strait is more near critical to *D*_{1} than *D*_{2} tides.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Criticality in Luzon Strait. Again, the strait is more near critical to *D*_{1} than *D*_{2} tides.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

### b. Discussion of the energy balance terms

The rough energy budget inferred from the observations is that there is approximately 7 kW m^{−1} of low-mode *D*_{1} energy approaching the continental slope (Fig. 26). Of that, maybe one-third of the mode-1 energy reflects (1.5 kW m^{−1}), whereas most of the mode-2 energy propagates onshore. The sum of the mode-1 and mode-2 energy at the plateau mooring LR1 balances is approximately 6.5 kW m^{−1} propagating onshore.

Sketch of rough diurnal (*D*_{1}) energy budget inferred from the observations. Local generation from a linear model (gray arrow) and the upper bound on the dissipation as extrapolated from MP1 are noted in gray, though they are not very well constrained. The shoreward estimate at MP1 is also an upper bound and only from modes 1 and 2.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Sketch of rough diurnal (*D*_{1}) energy budget inferred from the observations. Local generation from a linear model (gray arrow) and the upper bound on the dissipation as extrapolated from MP1 are noted in gray, though they are not very well constrained. The shoreward estimate at MP1 is also an upper bound and only from modes 1 and 2.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Sketch of rough diurnal (*D*_{1}) energy budget inferred from the observations. Local generation from a linear model (gray arrow) and the upper bound on the dissipation as extrapolated from MP1 are noted in gray, though they are not very well constrained. The shoreward estimate at MP1 is also an upper bound and only from modes 1 and 2.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

This picture is complicated by the high-mode energy flux inferred at LR1 (Table 1), though that energy could be subject to considerable noise or balanced by unresolved high-mode energy at MP1. It is also complicated by the presence of significant dissipation at both MP1 and the shipboard site, the simple extrapolation of which may lead us to expect as much as a 3 kW m^{−1} energy loss.

Similarly, there is the potential for local generation of diurnal energy (Fig. 15) of similar magnitude to the fraction of mode-1 energy we estimate as being reflected. We do not have any way to improve our knowledge of this system from just one mooring offshore. If the phase between the incoming waves and those generated locally were well separated, it might be possible to distinguish the relative energy in each, but that requires better knowledge of the geometry of the system than we have here.

Our budget for the semidiurnal tide is not as well constrained, with 3.5 kW m^{−1} observed offshore at MP1 but 7.5 kW m^{−1} observed onshore at LR1. This discrepancy is unlikely due to local generation because the *D*_{2} tides are quite weak here. Again, it is possible that the energy flux is simply being overestimated at LR1 because we are using the energy density multiplied by group speed rather than a proper estimate of flux.

Even for *D*_{1}, the energy budget presented above is very schematic and would require more data and modeling to make it more robust. The first limitation is that we have made vertical model decompositions. This is necessary for getting a reflection coefficient because we require the group speed of a progressive wave to compare against. We plan to test the robustness of this assumption in a series of numerical experiments where it is possible to separate eastward-propagating signals from westward-propagating signals using complex demodulation. To do the same thing with real data, a pair of moorings spaced significantly apart so that the phase lags could be determined would be necessary. For this experiment, another mooring farther east in the basin would have helped.

Similarly, we are hampered by our assumption of a two-dimensional wave. In the real ocean, the scattering from the three-dimensional shelf break will lead to out of phase signals that will increase *E* but not be readily apparent in the energy flux.

Finally, the fraction of energy reflected is difficult to determine, even for the simple linear step model. If all the incoming energy is mode 1, the answer is unambiguous. If the incoming energy has a broader modal content, the fraction of reflected energy becomes a very complicated function of the modal content, its phase, and the depth of the continental shelf. Evaluating this further is beyond the scope of this work simply because we do not have the appropriate data to compare with.

In terms of where the dissipation of low-mode tidal energy occurs, these observations and the linear model indicate that continental shelves have the potential to dissipate a large fraction of the low-mode energy. However, depending on the relative phase of the modes in the incoming tide, the theoretical results here indicate that, for many supercritical slopes, most of the energy may be reflected. It should be noted that our WKB-scaled plateau depth was exactly half the WKB water depth in the South China Sea, whereas most continental shelves are shallower, *F*_{2}/*F*_{1} = 0.5, almost all the energy would be reflected, with less dependence on the phasing of the incoming waves (Fig. 27). The dependence becomes more complicated for other values of *F*_{2}/*F*_{1} with the result weighted toward the single-mode values (Fig. 17).

Fraction of energy flux reflected from step-bathymetry continental shelves for incoming internal tides with *F*_{2} = 0.5*F*_{1} at different phases from one another. The SCS continental slope is unusually deep.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Fraction of energy flux reflected from step-bathymetry continental shelves for incoming internal tides with *F*_{2} = 0.5*F*_{1} at different phases from one another. The SCS continental slope is unusually deep.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Fraction of energy flux reflected from step-bathymetry continental shelves for incoming internal tides with *F*_{2} = 0.5*F*_{1} at different phases from one another. The SCS continental slope is unusually deep.

Citation: Journal of Physical Oceanography 41, 5; 10.1175/2010JPO4500.1

Of great interest then is what happens to the fraction of energy that reflects to the deep basin. If it radiates all the way back to Luzon Strait, then a lossy standing wave will set up, and this will fundamentally alter the generation mechanisms at the strait. If the dissipation occurs before the reflected energy makes it that far east, the question is, what mechanisms drive the dissipation?

Finally, the observations here are not so well constrained that we have adequately answered how “efficient” the reflection/transmission problem is at this locale. The large dissipations at MP1 and the shipboard survey indicate that a large fraction of the incoming energy may dissipate on the slope. Similar observations made on the Oregon continental slope indicate that there are many near-critical regions where local hydraulics and bores form (Nash et al. 2007), indicating that local roughness may be an important factor in dissipating the tide. Determining the fraction lost more precisely will need to be the subject of future research.

## Acknowledgments

Thanks to the captain, crew, and technical support staff of the *Ocean Researcher I*. Michael Goldin, Tyler Hughen, Oliver Sun, and Andrew Coockson all helped at sea. Thanks also to Harper Simmons for numerical model runs of the region and the other NLIWI PIs. Thanks to Laurence Coogan for discussing the geology of the South China Sea. Comments from two anonymous reviewers were very helpful in clarifying important aspects of the paper. Bruno Voisin pointed out that the reflection decomposition had been done before by Chapman and Hendershott (1981). This work was supported under US-ONR Grants N00014-08-1-0376, N00014-09-1-0219, N00014-05-1-0283, and N00014-05-0284.

## REFERENCES

Alford, M. H., and R. Pinkel, 2000: Observations of overturning in the thermocline: The context of ocean mixing.

,*J. Phys. Oceanogr.***30**, 805–832.Alford, M. H., and Z. Zhao, 2007: Global patterns of low-mode internal-wave propagation. Part II: Group velocity.

,*J. Phys. Oceanogr.***37**, 1849–1858.Alford, M. H., R. Lien, H. Simmons, J. Klymak, S. Ramp, Y. Yang, D. Tang, and M.-H. Chang, 2010: Speed and evolution of nonlinear internal waves transiting the South China Sea.

,*J. Phys. Oceanogr.***40**, 1338–1355.Aucan, J., M. A. Merrifield, D. S. Luther, and P. Flament, 2006: Tidal mixing events on the deep flanks of Kaena Ridge, Hawaii.

,*J. Phys. Oceanogr.***36**, 1202–1219.Buijsman, M. C., Y. Kanarska, and J. C. McWilliams, 2010: On the generation and evolution of nonlinear internal waves in the South China Sea.

,*J. Geophys. Res.***115**, C02012, doi:10.1029/2009JC005275.Cacchione, D. A., L. F. Pratson, and A. S. Ogston, 2002: The shaping of continental slopes by internal tides.

,*Science***296**, 724–727.Carter, G. S., and Coauthors, 2008: Energetics of

*M*2 barotropic-to-baroclinic tidal conversion at the Hawaiian Islands.,*J. Phys. Oceanogr.***38**, 2205–2223.Chang, M.-H., R.-C. Lien, T.-Y. Tang, E. A. D’Asaro, and Y.-J. Yang, 2006: Energy flux of nonlinear internal waves in northern South China Sea.

,*Geophys. Res. Lett.***33**, L03607, doi:10.1029/2005GL025196.Chapman, D., and M. Hendershott, 1981: Scattering of internal waves obliquely incident upon a step change in bottom relief.

,*Deep-Sea Res.***28A**, 1323–1338.Di Lorenzo, E., W. R. Young, and S. Llewellyn Smith, 2006: Numerical and analytical estimates of

*M*2 tidal conversion at steep oceanic ridges.,*J. Phys. Oceanogr.***36**, 1072–1084.Egbert, G. D., and S. Y. Erofeeva, 2002: Efficient inverse modeling of barotropic ocean tides.

,*J. Atmos. Oceanic Technol.***19**, 183–204.Farmer, D., Q. Li, and J. Park, 2009: Internal wave observations in the South China Sea: The role of rotation and non-linearity.

,*Atmos.–Ocean***47**, 267–280.Huang, R., 1994: Thermohaline circulation: Energetics and variability in a single-hemisphere basin model.

,*J. Geophys. Res.***99**, 12 471–12 485.Inall, M., T. Rippeth, C. Griffiths, and P. Wiles, 2005: Evolution and distribution of TKE production and dissipation within stratified flow over topography.

,*Geophys. Res. Lett.***32**, L08607, doi:10.1029/2004GL022289.Jan, S., C. Chern, J. Wang, and S. Chao, 2007: Generation of diurnal K1 internal tide in the Luzon Strait and its influence on surface tide in the South China Sea.

,*J. Geophys. Res.***112**, C06019, doi:10.1029/2006JC004003.Johnston, T. M. S., and M. A. Merrifield, 2003: Internal tide scattering at seamounts, ridges, and islands.

,*J. Geophys. Res.***108**, 3180, doi:10.1029/2002JC001528.Kelly, S. M., J. D. Nash, and E. Kunze, 2010: Internal-tide energy over topography.

,*J. Geophys. Res.***115**, C06014, doi:10.1029/2009JC005618.Klymak, J. M., and M. C. Gregg, 2004: Tidally generated turbulence over the Knight Inlet sill.

,*J. Phys. Oceanogr.***34**, 1135–1151.Klymak, J. M., and Coauthors, 2006a: An estimate of tidal energy lost to turbulence at the Hawaiian Ridge.

,*J. Phys. Oceanogr.***36**, 1148–1164.Klymak, J. M., R. Pinkel, C.-T. Liu, A. K. Liu, and L. David, 2006b: Prototypical solitons in the South China Sea.

,*Geophys. Res. Lett.***33**, L11607, doi:10.1029/2006GL025932.Klymak, J. M., R. Pinkel, and L. Rainville, 2008: Direct breaking of the internal tide near topography: Kaena Ridge, Hawaii.

,*J. Phys. Oceanogr.***38**, 380–399.Klymak, J. M., S. Legg, and R. Pinkel, 2010a: High-mode stationary waves in stratified flow over large obstacles.

,*J. Fluid Mech.***644**, 312–336.Klymak, J. M., S. Legg, and R. Pinkel, 2010b: A simple parameterization of turbulent tidal mixing near supercritical topography.

,*J. Phys. Oceanogr.***40**, 2059–2074.Kundu, P. K., 1990:

*Fluid Mechanics*. Academic Press, 638 pp.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.Legg, S., and J. Klymak, 2008: Internal hydrualic jumps and overturning generated by tidal flow over a tall steep ridge.

,*J. Phys. Oceanogr.***38**, 1949–1964.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.Lien, R.-C., T. Y. Tang, M. H. Chang, and E. A. D’Asaro, 2005: Energy of nonlinear internal waves in the South China Sea.

,*Geophys. Res. Lett.***32**, L05615, doi:10.1029/2004GL022012.Llewellyn Smith, S. G., and W. R. Young, 2003: Tidal conversion at a very steep ridge.

,*J. Fluid Mech.***495**, 175–191.Martini, K. I., M. H. Alford, J. D. Nash, E. Kunze, and M. A. Merrifield, 2007: Diagnosing a partly standing internal wave in Mamala Bay, Oahu.

,*Geophys. Res. Lett.***34**, L17604, doi:10.1029/2007GL029749.McPhee-Shaw, E. E., and E. Kunze, 2002: Boundary layer intrusions from a sloping bottom: A mechanism for generating intermediate nepheloid layers.

,*J. Geophys. Res.***107**, 3050, doi:10.1029/2001JC000801.Merrifield, M. A., and P. E. Holloway, 2002: Model estimate of

*M*2 internal tide energetics at the Hawaiian Ridge.,*J. Geophys. Res.***107**, 3179, doi:10.1029/2001JC000996.Munk, W., and C. Wunsch, 1998: Abyssal recipes II: Energetics of tidal and wind mixing.

,*Deep-Sea Res.***45**, 1977–2010.Nakamura, T., T. Awaji, T. Hatayama, K. Akitomo, T. Takizawa, T. Kono, Y. Kawasaki, and M. Fukasawa, 2000: The generation of large-amplitude unsteady lee waves by subinertial

*K*1 tidal flow: A possible vertical mixing mechanism in the Kuril Straits.,*J. Phys. Oceanogr.***30**, 1601–1621.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. Ocean. Technol.***22**, 1551–1570.Nash, J. D., M. H. Alford, E. Kunze, K. Martini, and S. Kelly, 2007: Hotspots of deep ocean mixing on the Oregon continental slope.

,*Geophys. Res. Lett.***34**, L01605, doi:10.1029/2006GL028170.Nissen, S., D. Hayes, Y. Bochu, Z. Weijun, C. Yongqin, and N. Xiaupin, 1995: Gravity, heat flow, and seismic constraints on the processes of crustal extension: Northern margin of the South China Sea.

,*J. Geophys. Res.***100**, 22 447–22 483.Niwa, Y., and T. Hibiya, 2004: Three-dimensional numerical simulation of M2 internal tides in the East China Sea.

,*J. Geophys. Res.***109**, C04027, doi:10.1029/2003JC001923.Rainville, L., and R. Pinkel, 2006: Observations of the propagation and nonlinear interaction of the internal tide generated at the Hawaiian Ridge.

,*J. Phys. Oceanogr.***36**, 1104–1122.Ramp, S. R., and Coauthors, 2004: Internal solitons in the northeastern South China Sea. Part I: Sources and deep water propagation.

,*IEEE J. Oceanic Eng.***29**, 1157–1181.Rudnick, D. L., and Coauthors, 2003: From tides to mixing along the Hawaiian Ridge.

,*Science***301**, 355–357.Smith, W., and D. Sandwell, 1997: Global sea floor topography from satellite altimetry and ship depth soundings.

,*Science***277**, 1956.St. Laurent, L., 2008: Turbulent dissipation on the margins of the South China Sea.

,*Geophys. Res. Lett.***35**, L23615, doi:10.1029/2008GL035520.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., S. Stringer, C. Garrett, and D. Perrault-Joncas, 2003: The generation of internal tides at abrupt topography.

,*Deep-Sea Res. I***50**, 987–1003, doi:10.1016/S0967-0637(03)00096-7.Thorpe, S. A., 1977: Turbulence and mixing in a Scottish loch.

,*Phil. Trans. Roy. Soc. London***286A**, 125–181.Zhao, Z., V. Klemas, Q. Zheng, and X.-H. Yan, 2004: Remote sensing evidence for baroclinic tide origin of internal solitary waves in the northeastern South China Sea.

,*Geophys. Res. Lett.***31**, L06302, doi:10.1029/2003GL019077.

^{1}

It was pointed out to us during revision that the method described here closely follows Chapman and Hendershott (1981).