the Creative Commons Attribution 4.0 License.

the Creative Commons Attribution 4.0 License.

# Internal tide energy flux over a ridge measured by a co-located ocean glider and moored acoustic Doppler current profiler

### Barbara Berx

### Gillian M. Damerell

Internal tide energy flux is an important diagnostic for the study of energy pathways in the ocean, from large-scale input by the surface tide to small-scale dissipation by turbulent mixing. Accurate calculation of energy flux requires repeated full-depth measurements of both potential density (*ρ*) and horizontal current velocity (** u**) over at least a tidal cycle and over several weeks to resolve the internal spring–neap cycle. Typically, these observations are made using full-depth oceanographic moorings that are vulnerable to being “fished out” by commercial trawlers when deployed on continental shelves and slopes. Here we test an alternative approach to minimize these risks, with

**measured by a low-frequency acoustic Doppler current profiler (ADCP) moored near the seabed and**

*u**ρ*measured by an autonomous ocean glider holding station by the ADCP. The method is used to measure the semidiurnal internal tide radiating from the Wyville Thomson Ridge in the North Atlantic. The observed energy flux (4.2±0.2 kW m

^{−1}) compares favourably with historic observations and a previous numerical model study.

Error in the energy flux calculation due to imperfect co-location of the glider and ADCP is estimated by subsampling potential density in an idealized internal tide field along pseudorandomly distributed glider paths. The error is considered acceptable (<10 %) if all the glider data are contained within a “watch circle” with a diameter smaller than 1∕8 the mode-1 horizontal wavelength of the internal tide. Energy flux is biased low because the glider samples density with a broad range of phase shifts, resulting in underestimation of vertical isopycnal displacement and available potential energy. The negative bias increases with increasing watch circle diameter. If watch circle diameter is larger than 1∕8 the mode-1 horizontal wavelength, the negative bias is more than 3 % and all realizations within the 95 % confidence interval are underestimates. Over the Wyville Thomson Ridge, where the semidiurnal mode-1 horizontal wavelength is ≈100 km and all the glider dives are within a 5 km diameter watch circle, the observed energy flux is estimated to have a negative bias of only 0.4 % and an error of less than 3 % at the 95 % confidence limit. With typical glider performance, we expect energy flux error due to imperfect co-location to be <10 % in most mid-latitude shelf slope regions.

Internal tides are a ubiquitous hydrodynamic feature over continental shelves and slopes as they are commonly generated at the shelf break by across-slope tidal flows (Baines, 1982; Pingree et al., 1986; Sharples et al., 2007). However, direct measurement of internal tides can be a challenge in these regions due to intense commercial fishing activity leading to an increased risk of oceanographic mooring loss (Sharples et al., 2013). Calculation of internal tide energy flux, a key diagnostic for the understanding of baroclinic energy pathways, requires repeated full-depth measurements of both potential density (*ρ*) and horizontal current velocity (** u**) over at least a tidal cycle (Nash et al., 2005). If an objective is to resolve the internal spring–neap cycle or observe the effect of seasonal changes in stratification on the internal tide field, repeated full-depth measurements over several weeks or months are required. Typically, these measurements are made using a full-depth oceanographic mooring incorporating an acoustic Doppler current profiler (ADCP) and a string of conductivity–temperature loggers (e.g. Hopkins et al., 2014), or a profiling mooring with a CTD (conductivity, temperature, and depth) and acoustic current meter (e.g. Zhao et al., 2012). On continental shelves and slopes, these full-depth moorings are vulnerable to being “fished out” by demersal and pelagic trawling activity.

Hall et al. (2017b) describe a novel alternative approach to minimize these risks, with ** u** measured by a low-frequency ADCP moored near the seabed and

*ρ*measured by an autonomous ocean glider holding station by the ADCP as a “virtual mooring”. Commercial fishing activity on continental shelves and their adjacent slopes is often intense because these regions are highly biologically productive. However, steps can be taken to reduce the risk of ADCP loss, including deploying deeper than 600 m, keeping mooring lines short, or using trawl-resistant frames. Being relatively small, gliders are unlikely to be fished out, and the risk can be further reduced by real-time evasive action in response to vessel proximity guided by the maritime automatic identification system (AIS). However, this alternate approach was not comprehensively tested by Hall et al. (2017b) because of glider navigation and telemetry problems. In this study we test the method using a co-located glider and ADCP dataset from the Wyville Thomson Ridge in the North Atlantic, a topographic feature that previous observations and numerical model studies suggest is an energetic internal tide generator (Sherwin, 1991; Hall et al., 2011). We also estimate the error in the energy flux calculation due to imperfect co-location of the glider and ADCP and find that, with typical glider performance, it is acceptable in most mid-latitude shelf slope regions.

Ocean gliders have previously been used to observe internal waves and internal tides (Rudnick et al., 2013; Rainville et al., 2013; Johnston and Rudnick, 2015; Boettger et al., 2015; Hall et al., 2017a), including the calculation of energy fluxes using current velocity measurements from gliders equipped with ADCPs (Johnston et al., 2013, 2015). However, ADCPs are not routinely integrated with commercially available glider platforms (Seaglider, Slocum, and SeaExplorer), in part due to their higher power requirement. Synergy with moored ADCP data allows accurate calculation of internal tide energetics without the endurance limitations and data analysis complexities of an ADCP-equipped glider (e.g. Todd et al., 2017).

In Sect. 2 the temporal resolution constraints of glider measurements are explained and the observations used in this study described. The calculation of internal tide energy flux from co-located glider and moored ADCP data is fully described in Sect. 3. Observations of the internal tide radiating from the Wyville Thomson Ridge are presented in Sect. 4 and compared with historic observations and a previous numerical model study. In Sect. 5 the error in the energy flux calculation due to imperfect co-location is estimated. Key results are summarized and discussed in Sect. 6.

High temporal resolution is crucial for internal tide observations; Nash et al. (2005) suggest that a minimum of four evenly distributed independent profiles of ** u** and

*ρ*are required per tidal cycle for an unbiased calculation of energy flux. Over continental shelves and upper shelf slopes, this temporal resolution is achievable using gliders. Typical glider vertical velocities are 15–20 cm s

^{−1}, so a complete dive cycle to 1000 m can take as little as 3 h. This yields eight profiles (four dives) per semidiurnal (≈12 h) tidal cycle, but near the surface and seabed the descending and ascending profiles converge in time so the number of independent samples is halved to four. For diurnal (≈24 h) internal tides, 16 profiles per cycle are possible. In shallower water the temporal resolution of glider measurements increases further; 40 min full-depth dives are achievable over a 200 m deep shelf break, yielding 36 profiles per semidiurnal tidal cycle. The depth-limiting factor of the methodology is the range of the ADCP. In narrowband mode, 75 kHz ADCPs have a maximum range of around 600 m (dependent on environmental conditions) so multiple ADCPs or additional current meters on the mooring line are required for sites between 600 and 1000 m deep.

The observations used to test the method were collected from the northern flank of the Wyville Thomson Ridge (WTR) in the North Atlantic (Fig. 1a). A Kongsberg Seaglider (SG613; Eriksen et al., 2001) was deployed from NRV *Alliance* between 2 and 5 June 2017 during the fourth Marine Autonomous Systems in Support of Marine Observations mission (MASSMO4). The glider was navigated from the deeper waters of the Faroe–Shetland Channel (FSC) to the WTR and held station for 40 h by a short oceanographic mooring, deployed 5 days previously from MRV *Scotia* (Fig. 1b). The mooring was sited close to the 800 m isobath and instrumented with an upwards-looking 75 kHz RDI Long Ranger ADCP at approximately 722 m and an Aanderaa Seaguard acoustic current meter at 784 m, yielding observations of horizontal current velocity over 78 % of the water column. When on-station by the ADCP mooring, the glider made repeated 2 h dives to 700 m or the seabed, whichever was shallower. This yielded approximately 12 profiles (six independent samples near the surface and seabed; 12 independent samples at mid-depth) per semidiurnal tidal cycle.

Glider location at the surface, before and after each dive, was given by GPS position. Subsurface sample locations were approximated by linearly interpolating surface latitude and longitude onto sample time. When on-station, the glider stayed within 2.5 km of the mooring and the mean horizontal distance between temporally coincident glider and ADCP measurements was 1.3 km. This spatial scattering of the glider data is small compared to the semidiurnal mode-1 horizontal wavelength over the WTR (≈100 km, calculated from the observed buoyancy frequency profile) and so the glider data are initially considered a fixed-point time series with no spatial–temporal aliasing.

As the glider was on-station for only 40 h, the co-located time series is not long enough to resolve the internal spring–neap cycle. As a result, *M*_{2} harmonic fits
to the glider and mooring data (Sect. 3) are contaminated with *S*_{2} variability. To acknowledge this, we refer to the estimated *M*_{2} component of the co-located time series as *D*_{2} following Alford et al. (2011). The comparative numerical model (Sect. 4.1) only includes the *M*_{2} tidal constituent so we refer to model diagnostics as *M*_{2}.

## 2.1 Data processing

The glider was equipped with a standard Sea-Bird Electronics conductivity–temperature (CT) sail sampling at 0.2 Hz and the data processed using the UEA Seaglider Toolbox (https://bitbucket.org/bastienqueste/uea-seaglider-toolbox, last access: 9 February 2017) following Queste (2014). Conductivity data were corrected for thermal hysteresis following Garau et al. (2011) and the Seaglider flight model regressed using a method adapted from Frajka-Williams et al. (2011). As the CT sail was unpumped, salinity samples were flagged when the glider's speed was less than 10 cm s^{−1} or it was within 8 m of apogee^{1}. Temperature–salinity profiles from descents and ascents were independently averaged (median value) in 5 m depth bins, typically with 4–5 samples per bin. Sample time was averaged into the same bins to allow accurate temporal analysis at all depths. Absolute Salinity (*S*_{A}), Conservative Temperature (Θ), and potential density (*ρ*) in each bin were calculated using the TEOS-10 equation of state (IOC et al., 2010).

The 75 kHz ADCP was configured in narrowband mode with 10 m bins and 24 pings per 20 min ensemble. The ADCP data were processed using Marine Scotland Science's standard protocols, including correction for magnetic declination and quality assurance based on error velocity, vertical velocity, and percentage good ping thresholds. The ADCP data were then linearly upsampled onto the same Δ5 m depth levels as the glider data. The acoustic current meter was configured with a 10 min sampling interval and linearly downsampled onto the same 20 min sampling interval as the ADCP. Good velocity data were recovered for all depth levels between 85 and 705 m, as well as 780–785 m. In addition to the ADCP and current meter measurements, horizontal velocity was inferred from GPS position and the Seaglider flight model using a dive-average current method (DAC; Eriksen et al., 2001; Frajka-Williams et al., 2011). DAC was only calculated for dives deeper than 500 m so that values were representative of the majority of the water column. All velocities were transformed into along-slope and across-slope components. We take the northern flank of the WTR to be orientated exactly northwest–southeast so along-slope (*u*) is positive southeast and across-slope (*v*) is positive northeast (down-slope).

The full 3-day glider time series of Conservative Temperature and Absolute Salinity is shown in Fig. 2a. A semidiurnal internal tide is evident as a vertical oscillation of the main pycnocline (centred around 550 m) with an amplitude up to 50 m and a period of ≈12 h. Temporally coincident ADCP and current meter measurements (Fig. 2b, c) show dominant semidiurnal periodicity and a reversal of baroclinic current velocity across the main pycnocline, characteristic of a low-mode internal tide. Mode-1 horizontal velocity, calculated from the observed buoyancy frequency profile, reverses at approximately 505 m, slightly above the pycnocline.

Following Kunze et al. (2002) and Nash et al. (2005), internal tide energy flux is calculated as $\mathit{F}=\u2329{\mathit{u}}_{\mathrm{bc}}^{\prime}{p}^{\prime}\u232a$. The method requires repeated full-depth measurements of *ρ* and ** u** over at least a tidal cycle in order to determine pressure perturbation (

*p*

^{′}) and baroclinic velocity (${\mathit{u}}_{\mathrm{bc}}^{\prime}$), respectively.

## 3.1 Pressure perturbation

For the 40 h window when the glider was on-station by the ADCP mooring, potential density anomaly is calculated by subtracting the window-mean density profile from measured potential density,

Before subtraction, $\stackrel{\mathrm{\u203e}}{\mathit{\rho}}\left(z\right)$ is smoothed with a 50 m gaussian tapered running mean (*σ*=10 m) to yield a suitable background density profile. Vertical isopycnal displacement is then calculated as

To separate *D*_{2} internal tide variability from other physical processes, *M*_{2} tidal period (*T*=12.42 h) harmonics are fit to *ξ* on each Δ5 m depth level following Emery and Thomson (2001). This analysis is only applied to depth levels between 10 and 675 m; near-surface and near-bottom bins are excluded because of high numbers of flagged samples and reduced temporal resolution due to the glider going into apogee above 700 m. To obtain a full-depth time series, the *D*_{2} component of *ξ* is linearly extrapolated assuming *ξ*=0 at the surface (*z*=0) and bottom ($z=-H$, where *H* is water depth). Buoyancy frequency squared, ${N}^{\mathrm{2}}=-g/{\mathit{\rho}}_{\mathrm{0}}(\partial \stackrel{\mathrm{\u203e}}{\mathit{\rho}}/\partial z)$, is also linearly extrapolated, assuming ${N}^{\mathrm{2}}={\mathrm{10}}^{-\mathrm{6}}$ s^{−2} at the surface and bottom. Pressure perturbation is then calculated by integrating the hydrostatic equation from the surface,

where ${p}_{\mathrm{surf}}^{\prime}$ is pressure perturbation at the surface due to the internal tide, determined by applying the baroclinicity condition for pressure,

Figure 3 shows potential density (Fig. 3c and d) and the *D*_{2} component of vertical isopycnal displacement (Fig. 3e and f) for the 40 h analysis window. The amplitudes and phases of the *D*_{2} component of *ξ* are shown in Fig. 4a and b.

## 3.2 Baroclinic velocity

For the same 40 h window, horizontal velocity perturbation is calculated,

where $\stackrel{\mathrm{\u203e}}{\mathit{u}}\left(z\right)$ is the window-mean horizontal velocity profile. There are three spatial gaps in the time series: above 85 m, between the ADCP and current meter (705–780 m including blanking distance), and from the current meter to the seabed (785–800 m). To obtain a full-depth time series, *u*^{′} is linearly interpolated between the ADCP and current meter, and extrapolated to the surface and the bottom using a nearest neighbour method. Baroclinic velocity is then calculated,

where ${\mathit{u}}_{\mathrm{bt}}^{\prime}$ is barotropic velocity, assumed here to equal the depth-mean velocity, calculated as

The *D*_{2} components of ${\mathit{u}}_{\mathrm{bc}}^{\prime}$ and ${\mathit{u}}_{\mathrm{bt}}^{\prime}$ are extracted using the same harmonic analysis method applied to *ξ*. Figure 3 shows barotropic (Fig. 3a and b) and baroclinic (Fig. 3c and d) velocities and the *D*_{2} components of barotropic (Fig. 3a and b) and baroclinic (Fig. 3e and f) velocities for the 40 h analysis window. The amplitudes and phases of the *D*_{2} component of ${\mathit{u}}_{\mathrm{bc}}^{\prime}$ are shown in Fig. 4a and b.

## 3.3 Internal tide energetics

Profiles of internal tide energy flux, available potential energy (APE), and horizontal kinetic energy (HKE) are calculated as

and

where 〈⋅〉 denotes an average (mean) over an integer number of *M*_{2} cycles and *ρ*_{0}=1028 kg m^{−3} is a reference density.

Maximum *D*_{2} vertical isopycnal displacement is 42 m and occurs at 565 m (Fig. 4a), within the main pycnocline. This is comparable with historic observations of a semidiurnal internal tide over the northern flank of the WTR. Sherwin (1991) analysed CTD data from a 17 h repeat station (30 min between casts) that was 6.7 km east of the mooring (Fig. 1b) and determined maximum *D*_{2} vertical isopycnal displacement to be 37 at 580 m, again within the pycnocline. Here, almost all of the APE is contained within the pycnocline (Fig. 4c) because maximum *ξ* occurs at a similar depth to maximum *N*^{2} ($\mathrm{4.9}\times {\mathrm{10}}^{-\mathrm{5}}$ s^{−2} at 525 m). *D*_{2} baroclinic velocity is maximum (≈20 cm s^{−1}) near-bottom (Fig. 4a), as is HKE (Fig. 4c). Depth-integrated HKE and APE are 5.5 and 1.7 kJ m^{−2}, respectively.

Both the across- and along-slope components of *D*_{2} internal tide energy flux are maximum (≈12.7 kW m^{−2}) near-bottom and go to zero at the depth of maximum *N*^{2} (Fig. 4d), characteristic of a low-mode internal tide with a pycnocline in the lower half of the water column. Depth-integrated energy flux magnitude is 4.2 kW m^{−1}, directed almost due east (7^{∘} anticlockwise from east). In comparison, Sherwin (1991) estimated the *D*_{2} mode-1 internal tide energy flux at the nearby CTD repeat station to be 4.7 kW m^{−1}, but was unable to diagnose the direction.

## 4.1 Model comparison

In Fig. 5 the observations are compared with the regional tide model described by Hall et al. (2011). The model is a configuration of the Princeton Ocean Model (POM; Blumberg and Mellor, 1987) for the FSC and WTR region, initiated with typical late-summer stratification and forced at the boundaries with *M*_{2} barotropic velocities (see Hall et al., 2011, for full details). Maximum *N*^{2} in the model is slightly higher than observed (Table 1) but the vertical distribution of stratification is similar; the main pycnocline is between 500 and 600 m. *M*_{2} internal tide generation occurs within the model domain, driven by barotropic tidal currents across isobaths, and is diagnosed as positive barotropic-to-baroclinic energy conversion (Fig. 5b). The northern flank of the WTR is an area of energetic internal tide generation, up to 4 W m^{−2}, and radiates an internal tide into the southern FSC. Modelled internal tide energy fluxes are spatially variable, but >5 kW m^{−1} at some locations (Fig. 5a). The mooring was located east of the most energetic generation and up-slope of the largest energy fluxes.

For direct comparison, the model output is interpolated onto the exact location of the mooring (Table 1). The modelled *M*_{2} internal tide energy flux is 6 %–7 % larger than the observed *D*_{2} energy flux, but within 10^{∘} of its direction. Maximum modelled vertical isopycnal displacement is 41 m (slightly smaller than observed) but is compensated by the higher maximum *N*^{2} and results in modelled APE being 30 % larger than observed; modelled HKE is 10 % smaller than observed.

## 4.2 Surface tidal ellipses

As well as measuring potential density by the ADCP mooring, the glider is used to infer a second estimate of barotropic velocity. Harmonic analysis is used to extract the *D*_{2} component of DAC velocity (all dives deeper than 500 m) and compared to the *D*_{2} component of ${\mathit{u}}_{\mathrm{bt}}^{\prime}$ from the ADCP and current meter. Barotropic velocity is highest in the across-slope direction (maximum 15 cm s^{−1}; Fig. 3a) and there is a very close match between the DAC and ADCP estimates (rms difference is 0.8 cm s^{−1}). In the along-slope direction, where barotropic velocity is lower (maximum 0.5 cm s^{−1}; Fig. 3b), the DAC estimate lags the ADCP estimate by 35 min but their amplitudes closely match (rms difference is 1.2 cm s^{−1}). The resulting surface tidal ellipses have similar semi-major axis lengths and phases (Table 1) but the DAC estimate is less eccentric (more circular) and rotated 3^{∘} anticlockwise (Fig. 5c). Compared with *M*_{2} surface ellipses from the regional tide model described by Hall et al. (2011), both observational estimates are less eccentric and have shorter semi-major axes (Table 1; Fig. 5c). However, the inclination of observed and modelled ellipses are comparable, with their semi-major axes orientated across-slope. This is the orientation required to generate an energetic internal tide at the WTR.

The separation of spatial and temporal variability is a common problem when interpreting glider data due to their slow speed (Rudnick and Cole, 2011) and imperfect positioning. In this context, the inability of the glider to perfectly hold station by the ADCP mooring leads to error in the calculation of internal tide energy flux (Sect. 3) due to mis-sampling of the spatially and temporally varying density field. An understanding of this error is important for both mission planning and interpretation of results. Other missions along the European continental slope (e.g. Hall et al., 2017a) have shown that a glider operating as a virtual mooring by repeatedly diving to 1000 m around a fixed station can maintain a “watch circle” with a diameter of approximately 5 km, i.e. all dives start and end within 2.5 km of the target location. The ability to do this is dependent on environmental conditions, particularly tidal and slope currents, but the lower limit is effectively set by the glide angle; a steep 45^{∘} glide angle will result in around 2 km horizontal travel over a complete dive cycle to 1000 m.

The size of the energy flux error is related to the length scale of the sampling cloud (*d*, the diameter of the watch circle) and the horizontal wavelength of internal tide being measured (*λ*). If *d*≪*λ* we can consider the glider data a fixed-point time series with no spatial–temporal aliasing, and so the error will be small. As *d* increases the glider will increasingly sample density at the wrong phase of the internal tide and so the error will increase because the measured pressure perturbation (${p}_{\mathrm{Glider}}^{\prime}$) will deviate from the pressure perturbation at the ADCP (${p}_{\mathrm{ADCP}}^{\prime}$), located at the centre of the watch circle. If *d*≃*λ* the glider will sample density at random phases of the internal tide and so ${p}_{\mathrm{Glider}}^{\prime}$ and ${p}_{\mathrm{ADCP}}^{\prime}$ will be uncorrelated.

Here we use a Monte Carlo approach to estimate the energy flux error. Potential density in an idealized internal tide field is subsampled along pseudorandomly distributed glider paths contained within watch circles of varying diameters. The “true” depth-integrated energy flux at the ADCP, ${\mathit{F}}_{\mathrm{true}}={\int}_{-H}^{\mathrm{0}}\u2329{\mathit{u}}_{\mathrm{bc}}^{\prime}{p}_{\mathrm{ADCP}}^{\prime}\u232a\phantom{\rule{0.125em}{0ex}}\mathrm{d}z$, is then compared with the “observed” depth-integrated energy flux, ${\mathit{F}}_{\mathrm{obs}}={\int}_{-H}^{\mathrm{0}}\u2329{\mathit{u}}_{\mathrm{bc}}^{\prime}{p}_{\mathrm{Glider}}^{\prime}\u232a\phantom{\rule{0.125em}{0ex}}\mathrm{d}z$. In both equations ${\mathit{u}}_{\mathrm{bc}}^{\prime}$ is baroclinic velocity at the ADCP. An idealized *M*_{2} multi-mode internal tide field is created for a 1000 m deep water column with uniform stratification (Appendix A). The mode-1 horizontal wavelength (*λ*) is 80 km and mode-1 vertical isopycnal displacement is 50 m, typical of mid-latitude shelf slope regions. Glider sampling is modelled as a group of twelve 1000 m dives (denoted here as a twelve-dive scenario) over 37 h (≈3 *M*_{2} cycles), within a watch circle of diameter *d*. Each dive is 2 h 50 min long, with 15 min at the surface between dives. Horizontal distance travelled during each dive cycle is between 1.5 and 4 km (typical of real glider missions), but there is no surface drift. The glider's path during each dive is determined by randomly selecting a start position within the watch circle then randomly selecting an end position 1.5–4 km away, but still within the watch circle. The start position of the following dive is the same as the end position. Potential density is linearly interpolated onto this pseudorandom glider path and the resulting density “observations” analysed using the method described in Sect. 3.1 to yield ${p}_{\mathrm{Glider}}^{\prime}$.

Nine cases are investigated, with *d* ranging from *λ*∕32 (2.5 km) to *λ*∕4 (20 km), and for each case 5000 different twelve-dive scenarios are simulated. A different random set of baroclinic mode phases is used for each scenario. Example pseudorandomly distributed glider paths for four cases are shown in Fig. 6. Energy flux relative error is defined as ${\mathit{F}}_{\mathrm{err}}=({\mathit{F}}_{\mathrm{obs}}-{\mathit{F}}_{\mathrm{true}})/{\mathit{F}}_{\mathrm{true}}$ so positive error indicates an overestimation and negative error indicates an underestimation. Similarly, APE relative error is defined as ${\mathrm{APE}}_{\mathrm{err}}=({\mathrm{APE}}_{\mathrm{obs}}-{\mathrm{APE}}_{\mathrm{true}})/{\mathrm{APE}}_{\mathrm{true}}$, where APE_{true} is “true” depth-integrated APE (calculated from *ξ*_{ADCP}) and APE_{obs} is “observed” depth-integrated APE (calculated from *ξ*_{Glider}).

## 5.1 Single-scenario example

A single twelve-dive scenario for the $d=\mathit{\lambda}/\mathrm{4}$ case is shown in Fig. 7 to highlight the impact of mis-sampling density on observed energy flux and APE. This is an extreme example, with all the glider dives 6–10 km from the ADCP (Fig. 7d), and features near-bottom internal tide intensification similar to that observed on the northern flank of the WTR. In this example, the error in measured density is maximum in the lower half of the water column (where *ξ*_{ADCP} is up to 80 m; Fig. 7a, b); the resulting *ξ*_{Glider} underestimates *ξ*_{ADCP} by up to 20 m and leads by up to 40 min. Observed energy flux and APE underestimate true energy flux and APE over the majority of the water column (Fig. 7c); depth-integrated observed energy flux and APE underestimate depth-integrated true energy flux and APE by 772 W m^{−1} (${\mathit{F}}_{\mathrm{err}}=-\mathrm{0.09}$) and 615 J m^{−2} (${\mathrm{APE}}_{\mathrm{err}}=-\mathrm{0.2}$), respectively.

## 5.2 Energy flux error

Histograms of *F*_{err} (0.005 wide bins) for four watch circle diameter cases are shown in Fig. 8a. The peaked distribution for the $d=\mathrm{1}/\mathrm{32}\mathit{\lambda}$ case broadens with increasing watch circle diameter as well as becoming biased towards negative error. The negative bias results from two related mechanisms. Firstly, the amplitude of *ξ*_{Glider} (and therefore ${p}_{\mathrm{Glider}}^{\prime}$) is typically underestimated for large watch circles because the glider samples density with a broad range of phase shifts, causing spectral smearing and poor harmonic fits to *ξ*. Secondly, maximum energy flux occurs when *p*^{′} and ${\mathit{u}}_{\mathrm{bc}}^{\prime}$ are exactly in phase so any error in the phase of ${p}_{\mathrm{Glider}}^{\prime}$, positive or negative, will also result in a negative bias.

*F*_{err} distributions for all nine watch circle diameter cases are shown in Fig. 8c, including the 99 % and 95 % confidence limits and the bias (median value). As watch circle diameter increases, the width of the confidence intervals increases and the bias becomes progressively more negative. For the $d=\mathrm{1}/\mathrm{32}\mathit{\lambda}$ case, *F*_{err} is ±0.04 at the 99 % limit and the bias is near zero (−0.002). For the $d=\mathrm{1}/\mathrm{4}\mathit{\lambda}$ case at the other extreme, *F*_{err} is 0 to −0.31 at the 99 % limit and the bias is −0.1.

## 5.3 APE error

Histograms of APE_{err} for four watch circle diameter cases are shown in Fig. 8b. Compared with *F*_{err}, the distributions are broader and with a more negative bias for small watch circles. The broader distribution is explained by the error in *ξ*_{Glider} being squared in Eq. (9). The negative bias is explained by the first mechanism described in Sect. 5.2. APE_{err} distributions for all nine watch circle diameter cases are shown in Fig. 8d. Similar to *F*_{err}, the width of the confidence intervals increases and the bias becomes progressively more negative as watch circle diameter increases. For $d=\mathrm{1}/\mathrm{32}\mathit{\lambda}$, APE_{err} is ±0.08 at the 99 % limit and the bias is only −0.005. For $d=\mathrm{1}/\mathrm{4}\mathit{\lambda}$, APE_{err} is 0.02 to −0.33 at the 99 % limit and the bias is −0.08. Unlike *F*_{err}, the bias converges towards a constant value for very large watch circles.

A novel approach to measuring internal tide energy flux using a co-located ocean glider and moored ADCP is tested using a dataset collected from the WTR in the North Atlantic. Gliders cannot perfectly hold station, even when operating as a virtual mooring, so error in the energy flux calculation due to imperfect co-location of the glider and ADCP is estimated by subsampling potential density in an idealized internal tide field along pseudorandomly distributed glider paths. If we consider the maximum acceptable energy flux error to be 0.1 (10 %), all the glider data must be contained within a watch circle with a diameter smaller than 1∕8 the mode-1 horizontal wavelength of the internal tide. Energy flux is biased low and the negative bias increases with increasing watch circle diameter. If watch circle diameter is larger than 1∕8 the mode-1 horizontal wavelength, the negative bias is more than −0.03 (3 %) and all realizations within the 95 % confidence interval are underestimates. When on-station over the WTR, the glider stayed within 2.5 km of the mooring so watch circle diameter *d*=5 km. The local *D*_{2} mode-1 horizontal wavelength *λ*≈100 km so the $d/\mathit{\lambda}=\mathrm{0.05}$ case (Table 2) is the most appropriate for the observations presented here. The observed energy flux is estimated to have a negative bias of only −0.004 (0.4 %) and an error of less than ±0.03 (3 %) at the 95 % confidence limit. This estimate does not include the effect of internal tide advection by the barotropic tide (Stephenson et al., 2016), which can lead to an additional negative bias if barotropic velocity amplitude is of a similar size to baroclinic phase speed. Over the WTR, *D*_{2} mode-1 phase speed is ≈2.2 m s^{−1} and barotropic velocity amplitude is <0.2 m s^{−1} so we expect this effect to be negligible for our observations.

At mid-latitudes, *D*_{2} mode-1 horizontal wavelength for a 1000 m deep water column is typically in the range 40-160 km. The results presented here suggest energy flux error due to imperfect co-location can be reduced to an acceptable level (10 %) if the glider maintains a 5 to 20 km diameter watch circle. In the absence of strong tidal and slope currents, a well-trimmed glider diving to 1000 m with a relatively steep glide angle can usually maintain a watch circle with a diameter of 5 km or less, so energy flux error will typically be <10 %. Where horizontal wavelengths are shorter, for example at lower latitudes or in shallower and less stratified water columns, a smaller watch circle will be required to maintain an acceptable level of error. In shallower water, smaller watch circles are generally achievable because horizontal travel over a complete dive cycle scales with dive depth. Diurnal internal tides have longer horizontal wavelengths so larger watch circles are acceptable. For mission planning, the mode-1 horizontal wavelength of a tidal frequency *ω* can be estimated, $\mathit{\lambda}=\mathrm{2}\mathit{\pi}{c}_{\mathrm{1}}/\sqrt{{\mathit{\omega}}^{\mathrm{2}}-{f}^{\mathrm{2}}}$, where *f* is the inertial frequency and ${c}_{\mathrm{1}}=NH/\mathit{\pi}$ is an approximation of mode-1 eigenspeed. If the assumption of uniform stratification is not appropriate, *c*_{1} can be calculated by solving the boundary value problem for a given *N*(*z*) (Gill, 1982). Table 2 can then be used to estimate the energy flux bias and error that can be expected for a given value of *d*∕*λ*.

Including the above estimate of error due to imperfect co-location, the observed depth-integrated *D*_{2} internal tide energy flux over the northern flank of the WTR is 4.2±0.2 kW m^{−1}. This is considerably larger than previous internal tide observations over the southeastern bank of the FSC: 0.2 kW m^{−1} (90 km northeast of the WTR; Hall et al., 2011) and 0.4–0.6 kW m^{−1} (105 km northeast of the WTR; Hall et al., 2017b), but small compared with some deep-ocean ridges, e.g. the Hawaiian Ridge (up to 33 kW m^{−1}; Lee et al., 2006) and Luzon Strait (up to 41 kW m^{−1}; Alford et al., 2011). More comparable to the WTR is the Mendocino Escarpment, where a ridge is orientated perpendicular to the continental slope and the observed energy flux is 7 kW m^{−1} (Althaus et al., 2003).

The 40 h co-located time series presented here is not long enough to resolve the internal spring–neap cycle. Peak neap tide occurred on yearday 153^{2}, 1 day before the majority of the co-located time series. Assuming the internal tide is generated locally at the WTR, the surface and internal spring–neap cycles will be in phase. The observed *D*_{2} energy flux is therefore representative of neap internal tide and so an underestimate of the true *M*_{2} internal tide. This may somewhat explain the slight underestimate compared to the *M*_{2}-only regional tide model. Interestingly, the CTD time series used by Sherwin (1991) was recorded 2 days after peak spring tide so is more representative of spring internal tide. The fact that two observational estimates of *D*_{2} vertical isopycnal displacement, 6.7 km apart and at different phases of the internal spring–neap cycle, are so similar implies that there are compensating spatial gradients in internal tide magnitude. The regional tide model shows the possible extent of these gradients and suggests that accurate siting of moorings is crucial for repeated long-term observations.

For future experiments, spatial gaps in the time series can be minimized with conductivity–temperature loggers and additional current meters on the mooring line. We have also shown that glider-inferred DAC can provide an accurate estimate of tidal current velocity that could be used to constrain barotropic velocity in the absence of full-depth data coverage by ADCPs and current meters. However, the major limitation of the dataset presented here is the short length of the co-located time series. Future glider missions will hold station by an ADCP mooring for several weeks to resolve the internal spring–neap cycle. Calculating *D*_{2} internal tide energetics in a 36 h moving window will yield a time-varying energy flux that can be related to seasonal changes in stratification, advection by mesoscale eddies, spatial and temporal patterns in internal tide-driven turbulent mixing, and the resulting biogeochemical response.

The Seaglider data were processed using the UEA Seaglider Toolbox (https://bitbucket.org/bastienqueste/uea-seaglider-toolbox, last access: 9 February 2017) and are available from the British Oceanographic Data Centre (https://doi.org/10.5285/9373933d-48c1-5a37-e053-6c86abc0e213; Wynn et al., 2019). The ADCP and acoustic current meter data are available from Marine Scotland (https://doi.org/10.7489/12217-1; Berx et al., 2019). Data analysis code is available on request from the corresponding author.

An idealized *M*_{2} multi-mode internal tide field is created for a 1000 m deep water column with uniform stratification (${N}^{\mathrm{2}}=\mathrm{6.1}\times {\mathrm{10}}^{-\mathrm{6}}$ s^{−2}). Horizontal current velocity, $\mathit{u}=(u,v)$, and vertical isopycnal displacement, *ξ*, are defined by summing the first 10 baroclinic modes,

and

where *u*_{n} and *ϕ*_{n} are the velocity amplitude and the phase of the *n*th baroclinic mode, respectively; $\mathit{\omega}=\mathrm{1.41}\times {\mathrm{10}}^{-\mathrm{4}}$ s^{−1} is the *M*_{2} frequency; and $f=\mathrm{1.26}\times {\mathrm{10}}^{-\mathrm{4}}$ s^{−1} is the inertial frequency at 60^{∘} N. *A*_{n}(*z*) and *B*_{n}(*z*) are the vertical structures of horizontal current velocity and vertical isopycnal displacement for each baroclinic mode, and are equivalent to cos (*n**π**z*∕*H*) and sin (*n**π**z*∕*H*), respectively, where *n* is mode number. Horizontal wavenumber ${k}_{n}=\sqrt{{\mathit{\omega}}^{\mathrm{2}}-{f}^{\mathrm{2}}}/{c}_{n}$, where ${c}_{n}=NH/n\mathit{\pi}$, is an approximation of mode eigenspeed (Gill, 1982). Velocity amplitude decays with mode number, ${u}_{n}={u}_{\mathrm{1}}{e}^{-\mathrm{0.5}(n-\mathrm{1})}$, where *u*_{1} is the mode-1 velocity amplitude. This decay rate results in a well-defined internal tide beam if velocity phase is approximately equal for each baroclinic mode. However, a different random set of baroclinic mode phases (*ϕ*_{n}) is used for each scenario simulated so internal tide beams are only apparent in a subset of scenarios. *u*_{1}=0.28 m s^{−1} yields a mode-1 vertical isopycnal displacement amplitude of 50 m but energy flux error and APE error are not sensitive to absolute amplitude. The time-varying potential density field is then

where $\stackrel{\mathrm{\u203e}}{\mathit{\rho}}\left(z\right)$ is a background density profile with a vertical gradient equivalent to *N*^{2}. Barotropic velocity (${\mathit{u}}_{\mathrm{bt}}^{\prime}$) and residual flow ($\stackrel{\mathrm{\u203e}}{\mathit{u}}$) are both zero so ${\mathit{u}}_{\mathrm{bc}}^{\prime}=\mathit{u}$.

RH led the glider mission, analysed the co-located glider and ADCP dataset, and developed the method for estimating glider sampling error. BB lead the mooring deployment and processing of the ADCP data. GD processed and quality-controlled the glider data. The paper was written by RH with input from the other authors.

The authors declare that they have no conflict of interest.

This article is part of the special issue “Developments in the science and history of tides (OS/ACP/HGSS/NPG/SE inter-journal SI)”. It is not associated with a conference.

SG613 is owned and maintained by the UEA Marine Support Facility. The glider and ADCP mooring were deployed as part of the fourth Marine Autonomous Systems in Support of Marine Observations mission (MASSMO4; funded primarily by the Defence Science and Technology Laboratory) and the Marine Scotland Science Offshore Monitoring Programme. The cooperation of the captain and crew of NRV *Alliance* (CMRE, Centre for Maritime Research and Experimentation) and MRV *Scotia* (Marine Scotland) are gratefully acknowledged. The glider data were processed by Gillian Damerell, the ADCP data were processed by Helen Smith and Barbara Berx, and the acoustic current meter data were processed by Jennifer Hindson and Helen Smith. Assistance with glider piloting was provided by the UEA Glider Group.

This paper was edited by Mattias Green and reviewed by two anonymous referees.

Alford, M. H., MacKinnon, J. A., Nash, J. D., Simmons, H., Pickering, A., Klymak, J. M., Pinkel, R., Sun, O., Rainville, L., Musgrave, R., Beitzel, T., Fu, K.-H., and Lu, C.-W.: Energy flux and dissipation in Luzon Strait: Two tales of two ridges, J. Phys. Oceanogr., 41, 2211–2222, https://doi.org/10.1175/JPO-D-11-073.1, 2011. a, b

Althaus, A. M., Kunze, E., and Sanford, T. B.: Internal tide radiation from Mendocino Escarpment, J. Phys. Oceanogr., 33, 1510–1527, 2003. a

Baines, P. G.: On internal tide generation models, Deep-Sea Res., 29, 307–338, 1982. a

Berx, B., Hindson, J., and Smith, H.: Moored data from NWZ-E monitoring site in the Faroe-Shetland Channel, Marine Scotland, UK, https://doi.org/10.7489/12217-1, 2019. a

Blumberg, A. F. and Mellor, G. L.: A description of a three-dimensional coastal ocean circulation model, in: Three-Dimensional Coastal Ocean Models, Vol. 4, edited by: Heaps, N. S., American Geophysical Union, Washington, DC, 1–16, 1987. a

Boettger, D., Robertson, R., and Rainville, L.: Characterizing the semidiurnal internal tide off Tasmania using glider data, J. Geophys. Res.-Oceans, 120, 3730–3746, https://doi.org/10.1002/2015JC010711, 2015. a

Emery, W. J. and Thomson, R. E.: Data Analysis Methods in Physical Oceanography, Elsevier, Amsterdam, 2 Edn., 654 pp., 2001. a

Eriksen, C. C., Osse, T. J., Light, R. D., Wen, T., Lehman, T. W., Sabin, P. J., Ballard, J. W., and Chiodi, A. M.: Seaglider: a long-range autonomous underwater vehicle for oceanographic research, IEEE J. Oceanic Eng., 26, 424–436, https://doi.org/10.1109/48.972073, 2001. a, b

Frajka-Williams, E., Eriksen, C. C., Rhines, P. B., and Harcourt, R. R.: Determining vertical water velocities from Seaglider, J. Atmos. Ocean. Tech., 28, 1641–1656, https://doi.org/10.1175/2011JTECHO830.1, 2011. a, b

Garau, B., Ruiz, S., Zhang, W. G., Pascual, A., Heslop, E., Kerfoot, J., and Tintoré, J.: Thermal lag correction on Slocum CTD glider data, J. Atmos. Ocean. Tech., 28, 1065–1071, https://doi.org/10.1175/JTECH-D-10-05030.1, 2011. a

Gill, A. E.: Atmosphere-Ocean Dynamics, Academic Press, 662 pp., 1982. a, b

Hall, R. A., Huthnance, J. M., and Williams, R. G.: Internal tides, nonlinear internal wave trains, and mixing in the Faroe-Shetland Channel, J. Geophys. Res., 116, C03008, https://doi.org/10.1029/2010JC006213, 2011. a, b, c, d, e, f

Hall, R. A., Aslam, T., and Huvenne, V. A. I.: Partly standing internal tides in a dendritic submarine canyon observed by an ocean glider, Deep-Sea Res. Pt. I, 126, 73–84, https://doi.org/10.1016/j.dsr.2017.05.015, 2017a. a, b

Hall, R. A., Berx, B., and Inall, M. E.: Observing internal tides in high-risk regions using co-located ocean gliders and moored ADCPs, Oceanography, 30, 51–52, https://doi.org/10.5670/oceanog.2017.220, 2017b. a, b, c

Hopkins, J. E., Stephenson, G. R., Green, J. A. M., Inall, M. E., and Palmer, M. R.: Storms modify baroclinic energy fluxes in a seasonally stratified shelf sea: Inertial-tidal interaction, J. Geophys. Res.-Oceans, 119, 6863–6883, https://doi.org/10.1002/2014JC010011, 2014. a

IOC, SCOR, and IAPSO: The international thermodynamic equation of seawater – 2010: Calculation and use of thermodynamics properties, in: Intergovernmental Oceanographic Commission, Manuals and Guides, 56, p. 196, UNESCO, 2010. a

Johnston, T. M. S. and Rudnick, D. L.: Trapped diurnal internal tides, propagating semidiurnal internal tides, and mixing estimates in the California Current System from sustained glider observations, 2006–2012, Deep-Sea Res. Pt. II, 112, 61–78, https://doi.org/10.1016/j.dsr2.2014.03.009, 2015. a

Johnston, T. M. S., Rudnick, D. L., Alford, M. H., Pickering, A., and Simmons, H. J.: Internal tide energy fluxes in the South China Sea from density and velocity measurements by gliders, J. Geophys. Res.-Oceans, 118, 3939–3949, https://doi.org/10.1002/jgrc.20311, 2013. a

Johnston, T. M. S., Rudnick, D. L., and Kelly, S. M.: Standing internal tides in the Tasman Sea observed by gliders, J. Phys. Oceanogr., 45, 2715–2737, https://doi.org/10.1175/JPO-D-15-0038.1, 2015. a

Kunze, E., Rosenfeld, L. K., Carter, G. S., and Gregg, M. C.: Internal waves in Monterey Submarine Canyon, J. Phys. Oceanogr., 32, 1890–1913, 2002. a

Lee, C. M., Kunze, E., Sanford, T. B., Nash, J. D., Merrifield, M. A., and Holloway, P. E.: Internal tides and turbulence along the 3000-m isobath of the Hawaiian Ridge, J. Phys. Oceanogr., 36, 1165–1182, 2006. a

Nash, J. D., Alford, M. H., and Kunze, E.: Estimating internal wave energy fluxes in the ocean, J. Atmos. Ocean. Tech., 22, 1551–1570, 2005. a, b, c

Pingree, R. D., Mardell, G. T., and New, A. L.: Propagation of internal tides from the upper slopes of the Bay of Biscay, Nature, 321, 154–158, 1986. a

Queste, B. Y.: Hydrographic observations of oxygen and related physical variables in the North Sea and Western Ross Sea Polynya, PhD thesis, School of Environmental Sciences, University of East Anglia, 215 pp., 2014. a

Rainville, L., Lee, C. M., Rudnick, D. L., and Yang, K.-C.: Propagation of internal tides generated near Luzon Strait: Observations from autonomous gliders, J. Geophys. Res., 118, 4125–4138, https://doi.org/10.1002/jgrc.20293, 2013. a

Rudnick, D. L. and Cole, S. T.: On sampling the ocean using underwater gliders, J. Geophys. Res., 116, C08010, https://doi.org/10.1029/2010JC006849, 2011. a

Rudnick, D. L., Johnston, T. M. S., and Sherman, J. T.: High-frequency internal waves near the Luzion Strait observed by underwater gliders, J. Geophys. Res., 118, 1–11, https://doi.org/10.1002/jgrc.20083, 2013. a

Sharples, J., Tweddle, J. F., Green, J. A. M., Palmer, M. R., Kim, Y.-N., Hickman, A. E., Holligan, P. M., Moore, C. M., Rippeth, T. P., Simpson, J. H., and Krivtsov, V.: Spring-neap modulation of internal tide mixing and vertical nitrate fluxes at a shelf edge, Limnol. Oceanogr., 52, 1735–1747, 2007. a

Sharples, J., Ellis, J. R., Nolan, G., and Scott, B. E.: Fishing and the oceanography of stratified shelf seas, Prog. Oceanogr., 117, 130–139, https://doi.org/10.1016/j.pocean.2013.06.014, 2013. a

Sherwin, T. J.: Evidence of a deep internal tide in the Faeroe-Shetland channel, in: Tidal Hydrodynamics, edited by: Parker, B. B., John Wiley & Sons, New York, 469–488, 1991. a, b, c, d, e

Stephenson, G. R., Green, J. A. M., and Inall, M. E.: Systematic bias in baroclinic energy estimates in shelf seas, J. Phys. Oceanogr., 46, 2851–2862, https://doi.org/10.1175/JPO-D-15-0215.1, 2016. a

Todd, R. E., Rudnick, D. L., Sherman, J. T., Owens, W. B., and George, L.: Absolute velocity estimates from autonomous underwater gliders equipped with Doppler current profilers, J. Atmos. Ocean. Tech., 34, 309–333, https://doi.org/10.1175/JTECH-D-16-0156.1, 2017. a

Wynn, R. B., Wihsgott, J. U., Palmer, M. R., Lichtman, I. D., Miller, P., Goult, S., Nencioli, F., Loveday, B. R., Jones, S., Inall, M. E., Dumont, E., Venables, E., Jones, O., Risch, D., Hall, R. A., Cauchy, P., Pierpoint, C., Doran, J., Mowat, R., and Damerell, G. M.: MASSMO 4 project ocean glider and autonomous surface vehicle data, British Oceanographic Data Centre, National Oceanography Centre, NERC, UK, https://doi.org/10.5285/9373933d-48c1-5a37-e053-6c86abc0e213, 2019. a

Zhao, Z., Alford, M. H., Lien, R.-C., Gregg, M. C., and Carter, G. S.: Internal tides and mixing in a submarine canyon with time-varying stratification, J. Phys. Oceanogr., 42, 2121–2142, https://doi.org/10.1175/JPO-D-12-045.1, 2012. a

^{1}

Apogee is the phase of the dive between descent and ascent, when the glider pitches upwards and increases its buoyancy. Flow through the conductivity cell is unpredictable during this phase and so salinity spikes are common.

^{2}

We refer to time using yearday, defined as decimal days since midnight on 31 December 2016 (e.g. noon on 31 January 2017 is yearday 30.5).

fished-outby commercial trawlers. We show that autonomous ocean gliders and acoustic Doppler current profilers can be used together to accurately measure the amount of energy carried by internal tides.