## 1. Introduction

Autonomous, buoyancy-propelled underwater gliders have become an important component of modern ocean observation systems (Testor et al. 2019). An advantage offered by the glider is their versatility in allowing for the combination of many different sensors on a single autonomous platform. One such recent use of gliders is in the measurement of ocean turbulent dissipation rates, and the resulting ocean mixing (Fer et al. 2014; Palmer et al. 2015; Schultze et al. 2017). Turbulence plays a key role in the transport of momentum, matter, and energy in the global oceans and coastal seas. However, in contrast to more standard ocean parameters such as temperature, salinity, and oxygen, the measurement of parameters that characterize ocean turbulence is relatively difficult. Consequently, datasets that contain directly measured parameters quantifying turbulence are relatively scarce (Waterhouse et al. 2014).

A primary reason for this scarcity is because the traditional method to measure turbulence profiles, i.e., the dissipation rate of turbulent kinetic energy (*ϵ*, in W kg^{−1}) with depth, is to use shear probes on a free-falling profiler dropped from a floating platform (usually a ship). This procedure, however, is labor intensive, and has a number of other drawbacks such as being confined to relatively calm weather, and involving uncertainties in the near-surface region due to possible wake effects of the floating platform. Using gliders as host platform, the labor cost can be significantly reduced, and measurements made possible in both stormy weather (Schultze et al. 2020; Carpenter et al. 2020), as well as in the undisturbed near-surface ocean (Wijesekera et al. 2020). The use of gliders as an observation platform have thus expanded our sampling capabilities of ocean turbulence, and contributed to a reduction in the scarcity of ocean turbulence measurements.

Despite these advantages in using gliders to measure ocean turbulence, the method is not without its challenges. The computation of *ϵ* from velocity shear probes can be shown to be proportional to the water velocity incident on the sensors to the fourth power (Merckelbach et al. 2019). When the glider is used as host platform, it is the glider speed through water *U* that enters the equation for *ϵ* ∝ *U*^{−4}. Since gliders are normally not equipped with devices to measure *U* directly, some type of glider flight model is needed to infer the incident water velocity [e.g., see Fer et al. (2014), St. Laurent and Merrifield (2017), Schultze et al. (2017), and Carpenter et al. (2020) for different examples in the use of such models]. Recently Merckelbach et al. (2019) derived a dynamic glider flight model and assessed the accuracy of this model for estimating the incident water velocity by comparing with the data from velocity sensors added to the glider. A key assumption made in the development of the model is that the glider flies in an ocean water column at rest. This assumption becomes problematic, however, near the ocean surface in the presence of surface gravity waves.

In the present paper, we extent the dynamic flight model of Merckelbach et al. (2019) by considering the effects of surface gravity waves on the motion of the glider itself, on the relative motion of water past the glider induced by these waves, and on estimates of the glider vertical velocity. We use linear wave theory to model the interaction between surface waves and a moving submersed body (i.e., the glider), and quantify the extent to which the presence of surface waves leads to additional errors in the estimation of the dissipation rate of turbulent kinetic energy when using dynamic glider flight models. This allows for an understanding, and recommendations, for using gliders to study ocean turbulence at the near-surface (as in, e.g., diurnal warm and rain layers), as well as in shallow coastal seas.

## 2. Relative glider velocities induced by surface waves

Ocean surface gravity waves propagate on the air–water interface causing an oscillation in elevation of the free surface. Below the surface, water parcels in such a wave field follow (almost closed) orbital trajectories. The amplitudes of displacement and velocity decrease with depth, and furthermore depend on the wave amplitude, the wave period, and the water depth. When sufficiently deep in the water column compared to the wavelength of the surface wave, the orbital motions are negligible, and the presence of surface waves is not expected to affect the flight of the glider. This changes, however, when the glider approaches the free surface, and enters a medium whose motion increasingly becomes affected by the orbital motions of the surface waves.

To model this process, we now consider a body of water with zero mean flow, excited by a train of monochromatic linear surface waves. For a neutrally buoyant submersed body, with a length scale much smaller than the wavelength, it is expected that the relative velocities vanish and the body follows the orbital motion exactly. Under normal flight conditions, however, a glider attains a nonzero buoyancy, corresponding to a volume change of 250 and 440 mL, for a standard and an extended buoyancy pump, respectively. Taken relative to the glider volume, which is typically on the order of 50 L, this amounts to about 0.5% and 1%, respectively. Thus, the departure from a perfectly neutrally buoyant body is small, and therefore, we shall see that relative motions between the water and the glider due to the wave orbitals are also small.

### a. Linear wave theory

*t*is time,

*ζ*the surface elevation,

*a*the wave amplitude,

*k*the wavenumber, and

*ω*the radial frequency with

*ω*= 2

*πT*

^{−1}, where

*T*is the wave period. The

*x*and

*z*coordinates represent the horizontal and vertical directions, respectively, with the wave traveling in the positive

*x*direction, and

*z*directed upward. The linearized problem for an even bottom at depth

*z*= −

*h*, is given by the velocity potential

*x*and

*z*directions, are denoted by

*u*

_{x}and

*u*

_{z}, and can be found by differentiation of the velocity potential with respect to

*x*and

*z*, respectively. The oscillating pressure perturbation produced by the wave motion

*ω*

^{2}=

*gk*tanh(

*kh*), and Bernoulli’s equation,

*g*is the acceleration due to gravity, and

*ρ*is the density of water, which is assumed to be constant. This yields the set of equations

### b. Wave-forced glider flight modeling

*m*

_{g}the glider mass,

**V**the glider velocity (relative to the seabed),

**F**the external force vector (composed of drag, lift, and buoyancy forces), and

**F**

_{a}a force vector resulting from the acceleration of the water around the body. One reason for an nonzero

**F**

_{a}to arise is when water around the body accelerates in response to an accelerating submerged body in a flow, that, in the absence of the body, would be steady. This effect can be accounted for by an added mass matrix (Newman 1977), and has been included in the dynamic glider flight model by Merckelbach et al. (2019). Another, complementary situation that gives rise to a nonzero vector

**F**

_{a}, is when the motion of the body is steady, while the flow (in the absence of the body) is unsteady. Such is the case for an ocean surface wave field, since the wave orbitals result in an accelerating flow in the water column below. The presence of the glider in such a flow causes a disturbance that contributes to

**F**

_{a}. Newman (1977, paragraph 4.17) treats this problem in detail for the one-dimensional case, where he makes “the fundamental assumption … that the non-uniformity of the flow is slowly varying, relative to the length scale of the body.” In the present context, this translates to an assumption that the wavelength is much longer than the length of the glider.

*ξ*,

*η*), where the

*ξ*axis points along the principal axis toward the nose of the glider, and the

*η*axis points upward, perpendicular to the principal axis (Fig. 1). Newman’s solution for

**F**

_{a}can then be written as

*m*

_{11}and

*m*

_{22}the added mass coefficients (see also Merckelbach et al. 2019). In this expression the first term on the right-hand side of (6) represents the effect due to an accelerating glider, and has the appearance of a force due to an added mass. The second term on the right-hand side of (6), which is not considered by Merckelbach et al. (2019), represents the effects of the accelerating fluid. This force is sometimes referred to as due to “effective volume.”

**V**into a water velocity due to wave motion

**u**and a relative glider velocity through the water

**v**, we can substitute (6) into (5). From (3) it follows that

**F**for an (

*x*,

*z*) frame of reference as

*θ*and the angle of attack

*α*are defined in Fig. 1. The forces due to buoyancy

*F*

_{B}, gravity

*F*

_{g}, lift

*F*

_{L}, and drag

*F*

_{D}, are (Merckelbach et al. 2010, 2019)

*ϵ*

_{c}is the coefficient of compressibility,

*P*the water pressure,

*α*

_{T}, the thermal expansion coefficient of the glider, Θ the water temperature, Θ

_{0}a reference temperature,

*S*the total surface area of the wings,

*α*, the lift coefficient

*C*

_{L}is assumed to be linear in

*α*, with proportionality coefficient

*A*, the lift angle coefficient.

^{1}

**v**. To solve this equation for

**v**, we perturb it relative to the steady-state solution. We write

**v**=

**v**

_{0}+

**v**

_{1}, and

**F**=

**F**

_{0}+

**F**

_{1}, where subscripts 0 and 1 denote the steady-state solution and the oscillating component due to the waves, respectively. Substituting the perturbations into (7), we can separate the result into a steady-state part and an oscillating part, respectively, as

**F**depends on

**v**through

*U*

^{2}and

*α*. The steady-state parameters can be represented in an (

*x*,

*z*) frame of reference as

**F**

_{0}to be expressed in

**v**

_{0}. Consequently,

*υ*

_{x0}and

*υ*

_{z0}can be determined from (11a).

*U*and

*α*, we can derive the linearized approximations (see appendix for details)

**F**

_{0}from the result, we obtain

**f**are given by

*ξ*,

*η*) coordinates, (15) reads

**(**

*g**t*) vector on the right hand side describes the forcing of the system by the surface wave motion through its dependence on

### c. Application to glider flight in surface waves

The natural response of the glider due to an external disturbance can be investigated by solving the eigenvalue problem for the homogeneous form of (22), that is, for **g**(*t*) = 0. In this case, we find two negative eigenvalues, indicating that the glider dynamic behavior returns to the steady-state conditions after a small perturbation. The longest response time, computed from the reciprocal of the smallest eigenvalue, is about 13 s for nominal flight conditions. The forcing term (24) introduces two more time scales: the wave period and a time scale *υ*_{z}) and the length scale of the vertical structure of the wave field. For the wave periods considered herein, the latter time scale, however, is long compared to the wave period, so that it is permissible to simplify the mathematical problem by regarding the depth dependency to be quasi constant in time. This means that the spatial gradients in the pressure anomalies in (24) can be evaluated for constant *z* during a wave period.

**g**(

*t*) as

**v**

_{1}further. We can, however, deduce from (26) that the amplitude of the relative velocity perturbation is proportional to the wave amplitude, and it exhibits a depth dependency, which, in the deep-water wave limit, is given by the factor

*e*

^{kz}. The behavior of the solution is explored graphically in Fig. 2, with results shown for waves for a range of wave periods, and an amplitude of 1 m in water depths of 40 and 1000 m. Furthermore, nominal flight conditions are assumed, for which the relevant parameter settings are listed in Table 1. Figure 2 shows the amplitude of the perturbed velocity in the

*ξ*direction (Fig. 2a), and

*η*direction (Fig. 2b) as a function of depth for a number of different wave periods. In general, we find that the amplitude of

*υ*

_{ξ1}is larger than that of

*υ*

_{η1}, and that the amplitudes of both components are largest near the water surface, where also shorter wave periods result in larger relative velocity perturbations. The vertical structure of the solution is shaped by the linear combination of the horizontal and vertical pressure gradients, so that a rapid decay in amplitude is observed with depth, but less so for larger wave periods and shallower water depths, facts that follow from linear wave theory in (4).

The amplitude of the perturbed velocity in (a) *ξ* direction and (b) *η* direction, as function of depth for water depths of 40 m (solid lines) and 1000 m (dashed lines), for a range of wave periods. The wave amplitude is set to 1 m, whereas nominal flight parameters values are assumed (pitch angle of 25° and a variable displacement achieved by the buoyancy pump of 250 mL). Further flight parameter settings are listed in Table 1. (c) The ratio of the amplitudes of the relative velocity perturbation and the wave orbital velocity.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

The amplitude of the perturbed velocity in (a) *ξ* direction and (b) *η* direction, as function of depth for water depths of 40 m (solid lines) and 1000 m (dashed lines), for a range of wave periods. The wave amplitude is set to 1 m, whereas nominal flight parameters values are assumed (pitch angle of 25° and a variable displacement achieved by the buoyancy pump of 250 mL). Further flight parameter settings are listed in Table 1. (c) The ratio of the amplitudes of the relative velocity perturbation and the wave orbital velocity.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

The amplitude of the perturbed velocity in (a) *ξ* direction and (b) *η* direction, as function of depth for water depths of 40 m (solid lines) and 1000 m (dashed lines), for a range of wave periods. The wave amplitude is set to 1 m, whereas nominal flight parameters values are assumed (pitch angle of 25° and a variable displacement achieved by the buoyancy pump of 250 mL). Further flight parameter settings are listed in Table 1. (c) The ratio of the amplitudes of the relative velocity perturbation and the wave orbital velocity.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

Flight model parameter settings.

The ratios of the amplitudes of the relative velocity perturbation **v**_{1} and the wave orbital velocity **u** are shown in Fig. 2c. The key deduction that can be made from this panel is that the magnitude of the wave-induced relative velocity of the glider is less than 1% of the orbital velocity magnitude. This means that the glider position is significantly influenced by the orbital motion of the waves, and the glider’s perturbed motion due to the surface waves follows, to a good approximation, the wave orbitals. We conclude that because the nominal glider speed through water is significantly larger than the relative velocity perturbations, the glider dynamics problem can be treated as if waves were absent. The implications of this small relative velocity perturbation for the estimation of the dissipation rate of turbulent kinetic energy as measured from microstructure sensor data, are further discussed in section 4a.

## 3. Surface wave effects in the estimation of vertical velocity

In the application of glider flight models to estimate the turbulent dissipation rate from microstructure sensors, a number of studies have utilized the vertical velocity of the glider as an input to the flight model (Schultze et al. 2017; St. Laurent and Merrifield 2017; Scheifele et al. 2018; Moldotsov et al. 2020). The determination of vertical glider velocities is arrived at through differentiation of the glider pressure measurement. In this section we show that alterations in the pressure measurement can arise from the presence of surface waves, and that these can affect the calculation of vertical glider velocity.

*z*

_{g}(

*t*) is the vertical position of the glider relative to the still water level, and, for the sake of simplicity, any variability in the in situ density is ignored. The vertical position of the glider can be decomposed into

*z*

_{0}, due to the wave motion, giving

*z*′ = 0, rather than to the instantaneous water surface elevation

*z*′ =

*ζ*. This can readily be verified if we set

*z*

_{0}= 0 in (31).

*z*

_{0}=

*z*

_{0}, yielding

*K*) and the pressure anomaly due to the vertical orbital excursion (second term in

*K*). For circular orbitals in the deep water limit, both pressure contributions counteract each other and cancel. Indeed, for circular orbitals at

*kh*≫ 1, it can be readily shown that

*K*→ 0. However, when waves enter shallow water, the orbital motion is forced to become more horizontal, rendering the orbitals elliptic. While the pressure anomaly remains unaffected, the vertical excursion length is reduced, and both effects become unbalanced, resulting in

*K*> 0.

*K*on depth and

*kh*. For deep-water waves,

*K*= 0 for all depth, as argued above. For intermediate and shallow water waves,

*K*is larger than zero and increases with depth. The glider depth that would be inferred from dividing the observed pressure by −

*ρg*, as is common oceanographic procedure, is given by

*z*

_{0}augmented with a modulated surface elevation signal. This signal vanishes for deep-water waves, but

*increases*with depth for shallow and intermediate water waves (Fig. 3).

(a) The behavior of *K* = cosh[*k*(*h* + *z*_{0})]/cosh(*kh*) – sinh[*k*(*h* + *z*_{0})]/sinh(*kh*) as function of the relative depth *z*/*h* for values of *kh* = {1.2, 2.1, ∞}. The depths correspond to a monochromatic wave with a wave period of 9 s. (b) Wave-induced error in the nondimensional vertical glider velocity observation, as function of depth and wave period for *a* = 1 m, *h* = 40 m, and *w*_{g} = 10 cm s^{−1}.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

(a) The behavior of *K* = cosh[*k*(*h* + *z*_{0})]/cosh(*kh*) – sinh[*k*(*h* + *z*_{0})]/sinh(*kh*) as function of the relative depth *z*/*h* for values of *kh* = {1.2, 2.1, ∞}. The depths correspond to a monochromatic wave with a wave period of 9 s. (b) Wave-induced error in the nondimensional vertical glider velocity observation, as function of depth and wave period for *a* = 1 m, *h* = 40 m, and *w*_{g} = 10 cm s^{−1}.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

(a) The behavior of *K* = cosh[*k*(*h* + *z*_{0})]/cosh(*kh*) – sinh[*k*(*h* + *z*_{0})]/sinh(*kh*) as function of the relative depth *z*/*h* for values of *kh* = {1.2, 2.1, ∞}. The depths correspond to a monochromatic wave with a wave period of 9 s. (b) Wave-induced error in the nondimensional vertical glider velocity observation, as function of depth and wave period for *a* = 1 m, *h* = 40 m, and *w*_{g} = 10 cm s^{−1}.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

A common use of the glider depth measurement is to map sensor parameters, such as temperature and salinity, onto a vertical coordinate. Since the time-averaged vertical coordinate of a water parcel that follows a wave orbital, is equal to *z*_{0}, this is also the preferred vertical coordinate to use in mapping water properties. Therefore, we consider the wave-dependent pressure (depth) anomaly to be a pressure (depth) error. To illustrate this error, Fig. 3a shows the specific example of a wave with 9 s period. In this case, near-bed depth errors can be as large as 30% of the surface wave amplitude in 40 m water depth, increasing to even 50% in water depths of 20 m; a potentially significant fraction of the total depth, depending on wave amplitudes.

*U*= 0.4 m s

^{−1}these terms can safely be neglected for waves that are long compared to the length of the glider. Thus, the observed vertical glider velocity contains an oscillating error term with an amplitude of

*aωK*. As an example, consider a coastal sea with a water depth of 40 m, and a wave amplitude of 1 m. The amplitude of the wave dependent term in the observed vertical velocity, relative to a nominal vertical glider speed of

*w*

_{g}= 0.10 m s

^{−1}, is shown in Fig. 3b. For small wave periods, the relative wave contribution is small. However, under windy conditions when wave periods easily reach 9 s, it can be seen that the wave-induced error can be even larger than the nominal vertical glider velocity.

## 4. Discussion

### a. Implications for turbulent dissipation rate measurements

In section 2, we examined how surface waves can induce an oscillatory flow relative to the glider’s steady flight. The analysis presented assumed idealized monochromatic waves of unit amplitude. The amplitudes of the oscillatory relative flow components were shown to be linearly proportional to the wave height and were a maximum near the surface (Fig. 2). We now demonstrate the implications of these findings in a practical context; we examine the extremes in wave-induced relative velocities that a glider is likely to encounter in realistic ocean settings.

Observed and modeled wind-wave data are often condensed into two representative parameters: a peak period *T*_{p}, where the wave spectrum reaches maximum spectral density, and a significant wave height *H*_{s}, the mean of the highest third of the observed wave heights (with wave height defined as twice the amplitude *a*). Together they serve to define a spectrum of waves that spans a range of frequencies *f* (in Hz), and amplitudes. To account for this range, we expand the monochromatic analysis of section 2 to account for a spectrum of frequencies that is characteristic of realistic ocean wave fields. This is done by considering the surface wave field as input, and the wave-induced relative glider velocities as outputs, to a system with a transfer function described by (27). The response of the relative velocities to an oscillating wave forcing—in signal processing normally referred to as *gain*—can readily be computed, and is shown for the near-surface region (evaluated at *z*_{0} = −1 m) in Fig. 4a as a function of *f*. The gain is seen to approach zero in the limit of *f* → 0, as expected. For the frequency range where waves still can be considered long compared to the glider length (wave periods > 3 s) the gain increases nearly linearly with increasing frequency. Therefore, the shapes of the wave and wave-induced relative velocity spectra are not congruent, rather the wave-induced relative velocity spectrum is biased toward higher frequencies.

(a) Gain of the *υ*_{ξ1} signal as function of frequency *f*. (b) *υ*_{ξ1} spectra (solid lines) for the two wave conditions listed in the legends, assuming a PM-wave spectrum. The dashed lines show the corresponding scaled PM spectra. (c) Estimated probability density functions of the instantaneous velocity *υ*_{ξ1} for the two wave conditions. The shaded regions delimit the velocity ranges comprising of 90% of the observations, whereas the dashed vertical lines denote the *T* = *T*_{p} and *H* = *H*_{s}.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

(a) Gain of the *υ*_{ξ1} signal as function of frequency *f*. (b) *υ*_{ξ1} spectra (solid lines) for the two wave conditions listed in the legends, assuming a PM-wave spectrum. The dashed lines show the corresponding scaled PM spectra. (c) Estimated probability density functions of the instantaneous velocity *υ*_{ξ1} for the two wave conditions. The shaded regions delimit the velocity ranges comprising of 90% of the observations, whereas the dashed vertical lines denote the *T* = *T*_{p} and *H* = *H*_{s}.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

(a) Gain of the *υ*_{ξ1} signal as function of frequency *f*. (b) *υ*_{ξ1} spectra (solid lines) for the two wave conditions listed in the legends, assuming a PM-wave spectrum. The dashed lines show the corresponding scaled PM spectra. (c) Estimated probability density functions of the instantaneous velocity *υ*_{ξ1} for the two wave conditions. The shaded regions delimit the velocity ranges comprising of 90% of the observations, whereas the dashed vertical lines denote the *T* = *T*_{p} and *H* = *H*_{s}.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

In computing the turbulent dissipation rate from microstructure measurements, it is important to quantify the *instantaneous* velocities (*υ*_{ξ1}) a glider may encounter. To gain insight into this, we determine probability density functions for *υ*_{ξ1}, computed from time series constructed for a spectrum of *υ*_{ξ1} by assuming a random phase model. [Since freak waves and other extreme wave events (e.g., Teutsch et al. 2020) are in this context statistically insignificant, they are not considered].

We illustrate this process here by adopting the empirical Pierson–Moskowitz (PM) energy spectrum for ocean waves (e.g., Tucker 1991), and computing the corresponding spectrum of *υ*_{ξ1}. Figure 4b shows the power spectra for *υ*_{ξ1} at *z*_{0} = −1 m, resulting from two PM spectra with *T*_{p} = 6 and 9 s. For reference, the PM wave input spectra, scaled by the power gains for the respective peak frequencies, are indicated by the dashed lines. Although the power spectra for *υ*_{ξ1} have roughly the same shape as the PM spectra, it is clearly seen that the velocity spectra have a stronger high-frequency content. Hence, the peak frequency and peak value of the wave-induced relative velocity spectrum are about 5% higher than their wave spectrum equivalents.

*υ*

_{ξ1}are shown in Fig. 4c. The instantaneous velocities have Gaussian-like distributions with a zero mean. The shaded regions denote the velocities that can be expected to occur in 90% of the observations in a sea state described by the PM spectra. The vertical dashed lines mark the amplitudes of

*υ*

_{ξ1}for monochromatic wave trains with

*T*=

*T*

_{p}and

*H*=

*H*

_{s}, and nearly coincide with the edges of the corresponding shaded region. Hence, for a sea state that can be described by a PM spectrum, roughly 95% of the wave-induced relative velocities

*υ*

_{ξ1}the glider encounters are smaller than those computed with a monochromatic wave with the PM spectral parameters

*T*

_{p}and

*H*

_{s}, denoted by

*H*

_{s}and

*T*

_{p}, are processed into probability distributions (as above), and presented in Fig. 5b (German Bight) and Fig. 5c (Bay of Biscay) using shades of red as a measure of the probability density. The PM spectrum is a one-parameter model, taking the wind speed as input, and models a fully developed sea state. A property of the PM spectrum is that there is an assumed relationship between

*T*

_{p}and

*H*

_{s}, as shown by the blue lines in Figs. 5b and 5c. This is clearly not the case for the datasets presented, which show a much broader distribution. This discrepancy can be attributed to both the model and the measurement data reflecting a sea state that is not fully developed due to time varying wind conditions and limited fetch. In addition, in shallow water (German Bight) bathymetry may also play a role. A multiparameter JONSWAP-type spectrum (e.g., Tucker 1991) may provide a more accurate description for the spectra of the wave data used here. However, because the simpler PM spectrum has a significantly broader frequency range than the JONSWAP spectrum, the probability density function of the relative vertical velocity, resulting from a JONSWAP spectrum is narrower, and therefore (38) provides a conservative estimate.

Site details of used wave data.

Near-surface wave-induced relative velocities. (a) Maximum amplitude of relative velocities in *ξ* and *η* directions for the German Bight (GB) and Bay of Biscay (BB), computed along the solid and dashed lines in the other panels. (b),(c) Estimated probability density maps of the wave climate in the German Bight (shallow water, 40 m water depth) and Bay of Biscay (deep water, 4700 m water depth), respectively, where the dashed black lines indicate the extremities of the observations, and the solid black lines (part of) the 90th percentile.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

Near-surface wave-induced relative velocities. (a) Maximum amplitude of relative velocities in *ξ* and *η* directions for the German Bight (GB) and Bay of Biscay (BB), computed along the solid and dashed lines in the other panels. (b),(c) Estimated probability density maps of the wave climate in the German Bight (shallow water, 40 m water depth) and Bay of Biscay (deep water, 4700 m water depth), respectively, where the dashed black lines indicate the extremities of the observations, and the solid black lines (part of) the 90th percentile.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

Near-surface wave-induced relative velocities. (a) Maximum amplitude of relative velocities in *ξ* and *η* directions for the German Bight (GB) and Bay of Biscay (BB), computed along the solid and dashed lines in the other panels. (b),(c) Estimated probability density maps of the wave climate in the German Bight (shallow water, 40 m water depth) and Bay of Biscay (deep water, 4700 m water depth), respectively, where the dashed black lines indicate the extremities of the observations, and the solid black lines (part of) the 90th percentile.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

The dashed lines in Figs. 5b and 5c roughly delimit the region of observed pairs of peak frequency and significant wave height for the respective datasets, as determined by eye. The solid lines mark the upper 95th percentiles. The amplitudes of the wave-induced relative velocities are computed along the dashed and solid lines, again evaluated near the surface at *z*_{0} = −1 m and presented in Fig. 5a, maintaining the same line style convention. We observe that *υ*_{η1} remains below 3 mm s^{−1}, even during the most extreme sea states (red and green curves). The component incident to the sensors, *υ*_{ξ1}, remains below 1 and 0.8 cm s^{−1} for the deep and shallow sites, respectively, with the 95th-percentile values dropping further to 7 and 5 mm s^{−1}, respectively. Based on these results, we find that a conservative upper bound on the wave-induced along glider velocities in the most extreme sea states is 1 cm s^{−1}. This translates into an error in the computation of turbulence dissipation rates, based on the uncertainty in *U*, that remains below 10% for a typical glider speed of 40 cm s^{−1}. This is small compared to the roughly factor-of-2 uncertainty that is commonly cited for dissipation measurements from shear probes (Moum et al. 1995), and is comparable to the errors associated with the glider flight model itself (Merckelbach et al. 2019).

Finally, we note that this analysis neglects the possibility of wave-breaking effects that are certainly present in many sea states that a glider is likely to encounter. This will affect the velocities in the very near surface (i.e., upper couple of meters), and this is a limitation to the analysis presented herein. This limitation should be kept in mind when glider flight in the near surface region is considered.

### b. Depth projection of measured parameters

A glider mission was run in the German Bight of the North Sea during the passage of a summer storm on 10 August 2019. A significant wave height of 4.5 m and a peak period of 9.3 s for the glider location (54°40.8′N, 6°43.8′E) were computed using a preoperational wave model (Staneva et al. 2014), run for the German Bight/southern North Sea region. At this location the water depth is ~34 m, so that waves at this peak period have *kh* ≈ 1.6. The near-bed horizontal orbital velocity amplitude decay [i.e., Eq. (4a)] amounts to approximately 30%, indicating that the influence of waves with 9 s period reached the sea bed. We therefore expect that near the sea bed, the depth of the glider, computed from the CTD pressure sensor, is at least partially corrupted with a wave-induced signal.

To remove this signal from the inferred depth, we apply a low-pass filter with, in this case, a cutoff frequency of 1/15 Hz to the observed (i.e., apparent) depth position *z*_{p}, as computed directly from the CTD pressure data. This was found to provide a reasonable estimate of the glider vertical position *z*_{p}, to the depth with wave effects removed *z*_{g}, the glider height above the sea bed. This is done using a glider-mounted, downward-looking Teledyne RDI Explorer 600 kHz Doppler velocity log (DVL), as described in Merckelbach et al. (2019). Subtracting this value of *z*_{g} from the water depth *z*_{0} gives an estimate of the vertical glider excursion due to wave orbital motions, and is shown by the orange curves in Fig. 6. The corresponding amplitude envelope of the vertical orbital excursion length is indicated by the orange shaded region.

Depth anomaly using CTD pressure data (blue) and depth anomaly using DVL altitude data (orange) for four consecutive profiles during 2043–2112 UTC 10 Aug 2019. The dashed line indicates the depth where the blue and orange envelopes coincide.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

Depth anomaly using CTD pressure data (blue) and depth anomaly using DVL altitude data (orange) for four consecutive profiles during 2043–2112 UTC 10 Aug 2019. The dashed line indicates the depth where the blue and orange envelopes coincide.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

Depth anomaly using CTD pressure data (blue) and depth anomaly using DVL altitude data (orange) for four consecutive profiles during 2043–2112 UTC 10 Aug 2019. The dashed line indicates the depth where the blue and orange envelopes coincide.

Citation: Journal of Atmospheric and Oceanic Technology 38, 7; 10.1175/JTECH-D-20-0206.1

The anomalies of both depth estimates appear to be reasonably well delimited by the respective envelopes, bearing in mind that the waves are not monochromatic, but have dominant frequencies spread around a band that is centered on a peak period of 9.3 s. This spread in frequency gives rise to wave groups, the effect of which is apparent in the measured data. The envelopes for the pressure and the altitude based anomalies coincide at a depth of approximately 23.5 m. It is also at this depth that the pressure (blue curves) and altitude (orange curves) anomalies coincide. Furthermore, it is noted that, both anomalies are seen to have the same phase in vertical space, which is consistent with both anomalies having the same phase in time, see Eqs. (33) and (36). The consistency of this analysis suggests that the preferred measure of glider depth is *z*_{0}, which can be estimated by applying a low-pass filter on the pressure-derived depth, using a cutoff frequency small enough to remove the wave-induced anomalies. Based on this recommendation, measured parameters, such as temperature and salinity, are projected to a depth coordinate that is not perturbed by the influence of surface waves.

## 5. Conclusions

In this work, we quantified the effects of a surface wave field on the flight characteristics of ocean gliders. Two primary effects were identified: (i) the orbital motion of the surface waves results in a flow past the glider that deviates from the prediction of standard flight models, and (ii) glider vertical positions, and therefore also vertical velocity estimates, can be altered by wave-induced pressure and depth perturbations to the measured pressure signal. These effects were quantified by using a linear surface wave field as a forcing term in a glider flight model that accounts for both added mass and “effective volume” terms, thus extending the Merckelbach et al. (2019) analysis to include surface wave effects. The results show that the wave-induced velocities relative to the glider are generally small, due to the small net (positive or negative) buoyancy of the glider. In other words, the glider’s perturbed motion follows closely the wave’s orbital motion, and experiences little relative motion resulting from the wave forcing. By applying the extended flight model to wave conditions in two coastal study regions, characteristic of a shelf sea and a deep ocean basin, we conclude that the wave-induced relative velocities remain less than 1 cm s^{−1} for the most extreme sea states. When using gliders as a platform for turbulence microstructure observations, the associated error in the dissipation rate of turbulent kinetic energy remains under 10%, and is considered insignificant for this type of measurement. The glider platform thus appears ideal for the study of near-surface ocean turbulence in the presence of nonbreaking surface waves.

On the other hand, surface wave effects can have a significant influence on the estimated depth, and vertical velocity, of the glider. This is only expected to occur in shallow water wave fields, and is due to an incomplete cancellation of pressure anomalies arising from both the vertical excursion of the glider in the water column, and the pressure perturbation of the wave orbital acceleration. Observations taken in the shallow North Sea during a summer storm confirm the prediction of an increasing amplitude of this effect with depth. We propose that an effective remedy to this effect is to apply a low-pass filter to the pressure data. The filter should be chosen with a cutoff frequency low enough to remove the dominant frequency in the wave signal. The filtered depth reading then corresponds to a depth that is representative of an undisturbed water column, and provides a sensible coordinate to map the observations to.

## Acknowledgments

We acknowledge financial support from the Helmholtz Association through the PACES II program, as well as the German Research Foundation (DFG). This paper is a contribution to the projects M6 and T2 of the Collaborative Research Centre TRR 181, “Energy Transfers in Atmosphere and Ocean” funded through DFG Grant 274762653.

## Data availability statement

All datasets used herein are publicly available at the following sources: wave buoy data platform Gascogne Buoy, WMO6200001 at https://www.emodnet-physics.eu/Map/platinfo/piroosplot.aspx?platformid=1043915, model data German Bight at https://codm.hzg.de/codm/ dataset “wave model for the German Bight,” and glider data North Sea 2019 at https://doi.org/10.5281/zenodo.4740627.

## APPENDIX

### Linearization of *U*^{2} and *α*

*U*

^{2}and

*α*. The perturbation approach employed in this work, favors the use of velocities

*υ*

_{x}and

*υ*

_{z}. We write

*υ*

_{i}=

*υ*

_{i0}+

*υ*

_{i1}, with

*i*∈ {

*x*,

*z*}. The indices 0 and 1 denote the steady-state and perturbed velocity, respectively. To express

*U*

^{2}in steady-state parameters and the perturbed parameters

*υ*

_{x1}and

*υ*

_{z1}, we write

*υ*

_{i1}≪

*υ*

_{i0}with

*i*∈ {

*x*,

*z*}. The angle of attack is written as

*α*=

*α*

_{0}+

*α*

_{1}and

*α*

_{1}≪

*α*

_{0}. To express

*α*

_{1}in terms of perturbed velocities, we write

## REFERENCES

Carpenter, J. R., A. Rodrigues, L. K. P. Schultze, L. M. Merckelbach, N. Suzuki, B. Baschek, and L. Umlauf, 2020: Shear instability and turbulence in a submesoscale front following a storm.

,*Geophys. Res. Lett.***47**, e2020GL090365, https://doi.org/10.1029/2020GL090365.Dingemans, M. W., 1997:

*Water Wave Propagation over Uneven Bottoms: Part 1—Linear Wave Propagation*. Advanced Series on Ocean Engineering, Vol. 13. World Scientific, 471 pp.Fer, I., A. K. Peterson, and J. E. Ullgren, 2014: Microstructure measurements from and underwater glider in the turbulent Faroe Bank channel overflow.

,*J. Atmos. Oceanic Technol.***31**, 1128–1150, https://doi.org/10.1175/JTECH-D-13-00221.1.Kundu, P., I. Cohen, and H. Hu, 2004:

*Fluid Mechanics.*2nd ed. Elsevier, 759 pp.Merckelbach, L., D. Smeed, and G. Griffiths, 2010: Vertical velocities from underwater gliders.

,*J. Atmos. Oceanic Technol.***27**, 547–563, https://doi.org/10.1175/2009JTECHO710.1.Merckelbach, L., A. Berger, G. Krahmann, M. Dengler, and J. R. Carpenter, 2019: A dynamic flight model for Slocum gliders and implications for turbulence microstructure measurements.

,*J. Atmos. Oceanic Technol.***36**, 281–296, https://doi.org/10.1175/JTECH-D-18-0168.1.Moldotsov, S., A. Annis, R. Amon, and P. Perez-Brunius, 2020: Turbulent mixing in a loop-current eddy from glider-based microstructure observations.

,*Geophys. Res. Lett.***47**, e2020GL088033, https://doi.org/10.1029/2020GL088033.Moum, J., M. Gregg, R. Lien, and M. Carr, 1995: Comparison of turbulence kinetic energy dissipation rate estimates from two ocean microstructure profilers.

,*J. Atmos. Oceanic Technol.***12**, 346–366, https://doi.org/10.1175/1520-0426(1995)012<0346:COTKED>2.0.CO;2.Newman, J., 1977:

*Marine Hydrodynamics.*MIT Press, 402 pp.Palmer, M., G. Stephenson, M. Inall, C. Balfour, A. Düsterhus, and J. Green, 2015: Turbulence and mixing by internal waves in the Celtic Sea determined from ocean glider microstructure measurements.

,*J. Mar. Syst.***144**, 57–69, https://doi.org/10.1016/j.jmarsys.2014.11.005.Scheifele, B., S. Waterman, L. Merckelbach, and J. R. Carpenter, 2018: Measuring the dissipation rate of turbulent kinetic energy in strongly stratified, low-energy environments: A case study from the Arctic Ocean.

,*J. Geophys. Res. Oceans***123**, 5459–5480, https://doi.org/10.1029/2017JC013731.Schultze, L. K. P., L. M. Merckelbach, and J. R. Carpenter, 2017: Turbulence and mixing in a shallow stratified shelf sea from underwater gliders.

,*J. Geophys. Res. Oceans***122**, 9092–9109, https://doi.org/10.1002/2017JC012872.Schultze, L. K. P., L. M. Merckelbach, and J. R. Carpenter, 2020: Storm-induced turbulence alters shelf sea vertical fluxes.

,*Limnol. Oceanogr. Lett.***5**, 264–270, https://doi.org/10.1002/lol2.10139.St. Laurent, L., and S. Merrifield, 2017: Measurements of near-surface turbulence and mixing from autonomous ocean gliders.

,*Oceanography***30**(2), 116–125, https://doi.org/10.5670/oceanog.2017.231.Staneva, J., A. Behrens, and N. Groll, 2014: Recent advances in wave modelling for the North Sea and German Bight.

,*Küste***81**, 233–254.Staneva, J., A. Behrens, and K. Wahle, 2015: Wave modelling for the German Bight coastal-ocean predicting system.

,*J. Phys. Conf. Ser.***633**, 012117, https://doi.org/10.1088/1742-6596/633/1/012117.Testor, P., and Coauthors, 2019: OceanGliders: A component of the integrated GOOS.

,*Front. Mar. Sci.***6**, 422, https://doi.org/10.3389/fmars.2019.00422.Teutsch, I., R. Weisse, J. Moeller, and O. Krueger, 2020: A statistical analysis of rogue waves in the southern North Sea.

,*Nat. Hazards Earth Syst. Sci.***20**, 2665–2680, https://doi.org/10.5194/nhess-20-2665-2020.Tucker, M. J., 1991:

*Waves in Ocean Engineering: Measurement, Analysis and Interpretation*. Marine Science, Ellis Horwood, 431 pp.Waterhouse, A., and Coauthors, 2014: Global patterns of diapycnal mixing from measurements of the turbulent dissipation rate.

,*J. Phys. Oceanogr.***44**, 1854–1872, https://doi.org/10.1175/JPO-D-13-0104.1.Wijesekera, H., D. Wang, and E. Jarosz, 2020: Dynamics of the diurnal warm layer: Surface jet, high-frequency internal waves, and mixing.

,*J. Phys. Oceanogr.***50**, 2053–2070, https://doi.org/10.1175/JPO-D-19-0285.1.

^{1}

In Merckelbach et al. (2019) the lift angle coefficient is denoted by *a*. In this work the symbol *a* is reserved for the wave amplitude.