Dissipation Rate Estimates from Microstructure and Finescale Internal Wave Observations along the A25 Greenland–Portugal OVIDE Line

Bruno Ferron Laboratoire de Physique des Océans, UBO-IFREMER-CNRS-IRD, Plouzané, France

Search for other papers by Bruno Ferron in
Current site
Google Scholar
PubMed
Close
,
Florian Kokoszka Laboratoire de Physique des Océans, UBO-IFREMER-CNRS-IRD, Plouzané, France

Search for other papers by Florian Kokoszka in
Current site
Google Scholar
PubMed
Close
,
Herlé Mercier Laboratoire de Physique des Océans, UBO-IFREMER-CNRS-IRD, Plouzané, France

Search for other papers by Herlé Mercier in
Current site
Google Scholar
PubMed
Close
, and
Pascale Lherminier Laboratoire de Physique des Océans, UBO-IFREMER-CNRS-IRD, Plouzané, France

Search for other papers by Pascale Lherminier in
Current site
Google Scholar
PubMed
Close
Full access

We are aware of a technical issue preventing figures and tables from showing in some newly published articles in the full-text HTML view.
While we are resolving the problem, please use the online PDF version of these articles to view figures and tables.

Abstract

A total of 96 finestructure and 30 microstructure full-depth vertical profiles were collected along the A25 Greenland–Portugal Observatoire de la Variabilité Interannuelle et Décennale en Atlantique Nord (OVIDE) hydrographic line in 2008. The microstructure of the horizontal velocity was used to calculate turbulent kinetic energy dissipation rates εvmp, where vmp refers to the vertical microstructure profiler. The lowest dissipation values (εvmp < 0.5 × 10−10 W kg−1) are found below 2000 m in the Iberian Abyssal Plain and in the center of the Irminger basin; the largest values (>5 × 10−10 W kg−1) are found in the main thermocline, around the Reykjanes Ridge, and in a 1000-m-thick layer above the bottom near 48°N. The finestructure of density was used to estimate isopycnal strain and that of the lowered acoustic Doppler current profiler to estimate the vertical shear of horizontal velocities. Strain and shear were used to estimate dissipation rates εG03 (Gregg et al.) associated with the internal wave field. The shear-to-strain ratio correction term of the finescale parameterization εG03 brings the fine- and microscale estimates of the dissipation rate into better agreement as Polzin et al. found. The latitude/buoyancy frequency term slightly improves the parameterization for weakly stratified waters. Correction term εG03 is consistent with εvmp within a factor of 4.5 over 95% of the profiles. This good consistency suggests that most of the turbulent activity recorded in this dataset is due to the internal wave field. The canonical globally averaged diffusivity value of order 10−4 m2 s−1 needed to maintain the global abyssal stratification (Munk) is only reached on the flank of the Reykjanes Ridge and in the region around 48°N.

Corresponding author address: Bruno Ferron, Laboratoire de Physique des Océans, IFREMER, Centre de Brest, ZI Pointe du Diable, CS 10070, 29280 Plouzané, France. E-mail: bruno.ferron@ifremer.fr

Abstract

A total of 96 finestructure and 30 microstructure full-depth vertical profiles were collected along the A25 Greenland–Portugal Observatoire de la Variabilité Interannuelle et Décennale en Atlantique Nord (OVIDE) hydrographic line in 2008. The microstructure of the horizontal velocity was used to calculate turbulent kinetic energy dissipation rates εvmp, where vmp refers to the vertical microstructure profiler. The lowest dissipation values (εvmp < 0.5 × 10−10 W kg−1) are found below 2000 m in the Iberian Abyssal Plain and in the center of the Irminger basin; the largest values (>5 × 10−10 W kg−1) are found in the main thermocline, around the Reykjanes Ridge, and in a 1000-m-thick layer above the bottom near 48°N. The finestructure of density was used to estimate isopycnal strain and that of the lowered acoustic Doppler current profiler to estimate the vertical shear of horizontal velocities. Strain and shear were used to estimate dissipation rates εG03 (Gregg et al.) associated with the internal wave field. The shear-to-strain ratio correction term of the finescale parameterization εG03 brings the fine- and microscale estimates of the dissipation rate into better agreement as Polzin et al. found. The latitude/buoyancy frequency term slightly improves the parameterization for weakly stratified waters. Correction term εG03 is consistent with εvmp within a factor of 4.5 over 95% of the profiles. This good consistency suggests that most of the turbulent activity recorded in this dataset is due to the internal wave field. The canonical globally averaged diffusivity value of order 10−4 m2 s−1 needed to maintain the global abyssal stratification (Munk) is only reached on the flank of the Reykjanes Ridge and in the region around 48°N.

Corresponding author address: Bruno Ferron, Laboratoire de Physique des Océans, IFREMER, Centre de Brest, ZI Pointe du Diable, CS 10070, 29280 Plouzané, France. E-mail: bruno.ferron@ifremer.fr

1. Introduction

Turbulent mixing plays a key role in the global water mass budget of the ocean and contributes to the maintenance of the stratification. Since turbulent mixing is a key component of the global overturning circulation (Munk and Wunsch 1998) and because of its strong spatial heterogeneity and temporal intermittence, continuous efforts to measure its distribution and magnitude and to understand the underlying mechanisms need to be pursued. Turbulent mixing is usually quantified locally from microstructure data obtained from specific instruments with O(1 cm) resolution (e.g., Oakey 1982) and from the finestructure data obtained with the widely used conductivity–temperature–depth (CTD) and lowered acoustic Doppler current profilers (LADCP). The finestructure is either used to quantify Thorpe scales (Thorpe 1977; Dillon 1982) or to diagnose the intensity of the internal wave field (Gregg 1989). At larger scales, turbulent mixing can also be quantified from tracer release experiments and from tracer budgets (e.g., Ledwell et al. 1993, 2000; Hogg et al. 1982). During the last decade, an increasing number of studies used the finescale properties of the internal wave field to quantify turbulent mixing due to the internal wave environment (e.g., Hibiya and Nagasawa 2004; Naveira Garabato et al. 2004a,b; Walter et al. 2005; Kunze et al. 2006; Lauderdale et al. 2008; Stöber et al. 2008; MacKinnon et al. 2008; Daae et al. 2009; Tian et al. 2009). Other studies have compared and validated the local approaches based on microstructure, which contains mixing from different mechanisms, and finescale parameterizations (Gregg 1989; Polzin et al. 1995; Gregg et al. 2003; Hibiya et al. 2012; Frants et al. 2013; Sheen et al. 2013; Waterman et al. 2013) based on the internal wave field, which is presumed to be a major contributor of ocean interior mixing.

In this paper we present full-depth dissipation rate profiles obtained with a vertical microstructure profiler along the A25 Greenland–Portugal Observatoire de la Variabilité Interannuelle et Décennale en Atlantique Nord (OVIDE) hydrographic line and compare them to finescale parameterizations of the dissipation rate based on internal wave–wave interaction theory. We tested several parameterizations of the dissipation rate, which depends on finescale observations of internal wave velocity shear and strain, by comparing them with dissipation data. Gregg (1989) proposed a shear-based parameterization for the turbulent dissipation rate:
e1
where ε0 = 7 × 10−10 W kg−1, and N/N0 and being the buoyancy frequency and the shear variance, respectively, scaled by their respective values from the Garrett and Munk (GM) model (Garrett and Munk 1975; Cairns and Williams 1976). This parameterization has been subsequently refined by several authors. To better represent regions that depart from the canonical GM description, Polzin et al. (1995) incorporated a frequency correction h leading to a new parameterization, defined as
e2
where
eq1
and is the shear (Vz)-to-strain (ζz) variance ratio. Gregg et al. (2003) added a term j dependent on latitude to improve the consistency with the dissipation rates found in the equatorial region that are weaker than at midlatitudes. Gregg et al.’s (2003) parameterization reads
e3
where
eq2
where f is the inertial frequency and . Following Wijesekera et al. (1993), internal wave shear variance ratios in those parameterizations can be replaced by internal wave strain variance ratios (e.g., Kunze et al. 2006), leading to
e4
where
eq3
In the following we present evaluations of εG89, εP95, and εG03, and compare them with dissipation rates measured with a microstructure profiler along the OVIDE transect that spans both a large latitudinal band (40°–59°N; Fig. 1) and a large range of buoyancy frequencies from 1.6 × 10−4 s−1 to 3.2 × 10−3 s−1.
Fig. 1.
Fig. 1.

Position of CTD–LADCP stations (yellow and red dots) and VMP stations (red dots) along the A25 OVIDE hydrographic line. The blue arrow indicates a station (48°02′N, 20°51′W) with enhanced dissipation rates at depth (topography from Smith and Sandwell 1997).

Citation: Journal of Atmospheric and Oceanic Technology 31, 11; 10.1175/JTECH-D-14-00036.1

2. Dataset

The OVIDE project consisted of repeating the A25 Greenland–Portugal hydrographic line every 2 years from 2002 to 2010 (Mercier et al. 2014). Each survey was composed of nearly 100 CTD–LADCP profiles (Fig. 1, yellow and red dots). Measurements lasted 3 weeks usually in June/July. During 2008, an autonomous vertical microstructure profiler (VMP; www.rocklandscientific.com) able to dive to 6000 m was deployed at 30 CTD locations (Fig. 1, red dots). The geographical distribution of the VMP deployments was not as uniform as that of the CTD. In particular, no VMP stations were done on either the southeastern flank of the Reykjanes Ridge because of heavy sea conditions or on the Greenland continental slope because of a technical problem with the instrument. The VMP was deployed 20–30 min before the CTD cast began and about half a nautical mile away from the CTD station location to avoid any collision. The VMP was retrieved after the end of the CTD upcast. The VMP was programmed before deployment to stop its profile 70 m above the bottom, while the CTD–LADCP profiles were stopped 10–15 m above the bottom.

a. Microstructure data

The VMP is a 3-m-long instrument that weights 150 kg. Microstructure sensors are located at the bottom of the instrument, so that only data from the downcast are useful. On the OVIDE cruise, the typical downcast fall rate of the instrument was 0.7 m s−1, a velocity high enough to apply Taylor’s frozen turbulence hypothesis, which allows a conversion from time to space using a smoothed instrument fall speed profile.

The VMP carried two Rockland Scientific fast thermistors (FP-07), one Sea-bird Electronics (SBE) microconductivity sensor (SBE-7), and two Rockland Scientific velocity shear probes (SPM-38). Microstructure channels were sampled at 512 Hz. A low-pass analog filter with a cutoff frequency of 165 Hz was applied on all microstructure channels to avoid aliasing.

Data from the two shear probes were used to calculate dissipation rates. The first 10 m of the dissipation profiles were discarded, since data were frequently contaminated by ship wake. In this study, fast thermistors and microconductivity data were only used to visually check the turbulence signal given by the shear probes. Shear data were first edited to remove spikes. Assuming isotropy, the dissipation rate of turbulent kinetic energy per unit mass reads ε = 7.5νuz2〉, where ν is the viscosity of seawater and uz is the vertical shear of the horizontal velocity u. The velocity shear variance 〈uz2〉 was estimated by a spectral analysis of the velocity shear signals. The velocity shear spectra were calculated using half-overlapping 512-point segments. Each segment was linearly detrended and windowed with a Hanning taper function prior to FFT calculation. Spectral coherence between the velocity shear data and the three accelerometer data was used to remove vibrational contamination (Levine and Lueck 1999; Goodman et al. 2006). Because of their finite size, velocity shear probes tend to average the smallest turbulent eddies, thereby underestimating shear variance at high wavenumbers (Oakey 1982). Velocity shear spectra were corrected for velocity shear probe averaging using a single-pole filter for the wavenumber response of the probe with a cutoff wavenumber of 48 cycles per meter (cpm; Macoun and Lueck 2004). To obtain velocity shear variance, shear spectra were integrated up to a cutoff wavenumber kc found iteratively. The procedure followed that described by Moum et al. (1995). For the first iteration, kc was set to 15 cpm and gave a first dissipation rate estimate from which the Kolmogorov wavenumber was estimated. For the next iteration, kc was increased to ks or, to avoid noise contamination at high wavenumbers, to a maximum of 50 cpm. The Kolmogorov wavenumber was then updated. The iterative process was stopped when the relative difference between kc and ks was lower than 0.001%, or when kc was equal to 50 cpm or after 10 iterations. At the end of the iteration process, if kc < ks, then velocity shear variance loss was accounted for assuming Nasmyth (1970) turbulence spectrum form. The background noise on the dissipation rate profile was set to 5 × 10−11 W kg−1. It was estimated from quiescent profiles of the deep ocean located in the Iberian Abyssal Plain. Following Gregg (1999) criteria, we observed that for dissipation rates lower than 5 × 10−11 W kg−1, spectra of microstructure velocity shear tend to depart from Nasmyth spectra. Examples of microstructure velocity shear spectra calculated at two different depths are provided as an illustration for a profile located at 42°35′N, 15°27′W (Fig. 2). While the spectrum followed the Nasmyth spectral shape at a depth of 4200 m and for a dissipation rate of 7.2 × 10−11 W kg−1 (Fig. 2, gray curves), it clearly departed from the Nasmyth spectral shape at a depth of 4270 m for a dissipation rate of 2.6 × 10−11 W kg−1 (Fig. 2, black curves). Note that for those two examples, power spectral density (PSD) rise was due to the presence of noise for wavenumbers larger than 13 cpm. This does not influence the dissipation rate estimate though, since Kolmogorov wavenumbers (viscous cutoff wavenumber) that set the upper bound of the spectral integration in the dissipation rate calculation are lower than 13 cpm (Fig. 2, vertical black and gray lines).

Fig. 2.
Fig. 2.

Example of an average of four PSD of microstructure velocity shear (solid) observed with the VMP for the profile located at 42°35′N, 15°27′W and at a depth of 4200 (gray) and 4270 m (black). Dashed curves are matching Nasmyth (1970) spectra for ε = 7.2 × 10−11 W kg−1 (gray) and for ε = 2.6 × 10−11 W kg−1 (black). Vertical lines give the corresponding Kolmogorov wavenumber.

Citation: Journal of Atmospheric and Oceanic Technology 31, 11; 10.1175/JTECH-D-14-00036.1

b. Finestructure data

On the OVIDE 2008 cruise, one upward-looking 300-kHz and one downward-looking 150-kHz RD Instruments workhorse LADCPs mounted on the CTD carousel provided observations of vertical shear of horizontal currents with 16-m vertical resolution. LADCP data were processed with the velocity-inversion method initially developed by Visbeck (2002).

To compute the velocity shear variance, each velocity profile was split into 640-m segments with 75% overlap. Each segment was linearly detrended prior to FFT calculation. Buoyancy-normalized shear variance was calculated as the integral, from the lowest wavenumber kmin = 2π/640 rad m−1 up to a maximum wavenumber kmax determined for each segment, of the velocity spectra S(u) multiplied by a transfer function T(k)(see discussion below): , where kmax is the minimum of knoise and kc, knoise is the wavenumber above which noise dominates the signal (see discussion below), and kc is the cutoff wavenumber for which the integrated velocity shear variance equals 0.66 , where is the segment-averaged buoyancy frequency. The GM velocity shear spectra were also integrated from kmin to kmax. Finding kc for which (i.e., the variance of the GM velocity spectra integrated up to its cutoff wavenumber of 2π/10 m−1) takes into account the fact that the cutoff wavenumber decreases with increasing energy level in the shear (Polzin et al. 1995; Gregg et al. 2003).

Polzin et al. (2002) showed that spectral corrections T(k) need to be applied to compensate for loss of velocity variance at a wavelength shorter than about 200 m because of the various steps of the acquisition and processing methods of the LADCP data. While Polzin et al. (2002) processed their LADCP data using the shear method (Firing and Gordon 1990; Fischer and Visbeck 1993), Thurnherr (2012) updated the spectral corrections for the velocity-inversion method and showed that it naturally produces less damping at short wavelengths. With the velocity-inversion method, the shear spectra need to be corrected at least for range averaging and depth binning. A supplementary term arises when data are prebinned into superensembles before the inversion, as done for our dataset. When instrument tilting is important, an additional term also applies. Considering the weak averaged tilt (mean tilt angle of 1.8°) due to our carousel configuration, no tilt correction was applied in this study. As Polzin et al. (2002) suggested, beam separation effects were ignored. Spectral correction for our data processing then takes the form of the transfer function , where k is the vertical wavenumber and Δz = 16 m represents the ADCP transmit pulse length, range gate length, superensemble preaveraging interval, and bin length of the output grid. An illustration of an ensemble-averaged raw (uncorrected) velocity shear spectrum scaled by the buoyancy frequency and its corrected version is given for an ensemble of profiles that were taken away from steep bathymetry, so that we expected a relatively flat spectral shape close to that of the GM model (Fig. 3). The ensemble was made of 589 spectra in the depth range 96–2656 m, from longitude 10° to 25°W and from latitude 40° to 52°N. In this example, the raw spectrum falls off by 60% from wavenumbers 0.04 rad m−1 (wavelengths of 160 m) to wavenumbers 0.1 rad m−1 (wavelength of 60 m). Such a damping is typical of the velocity-inversion method (see Thurnherr 2012, his Fig. 1). After correction, the spectral shape is relatively flat up to wavenumbers close to 0.1 rad m−1, above which noise dominates.

Fig. 3.
Fig. 3.

Comparison of ensemble-averaged buoyancy-normalized shear spectra without (raw) and with the use of the transfer function T(k) (see text). The ensemble was made of 589 spectra in the depth range 96–2656 m and from longitude 10°–25°W and latitude 40°–52°N.

Citation: Journal of Atmospheric and Oceanic Technology 31, 11; 10.1175/JTECH-D-14-00036.1

Ensemble-averaged shear spectra clearly show a depth-dependent contamination by noise at high wavenumbers for which the rms noise Srms strongly dominates the signal (Fig. 4a). To estimate the increase in the single-ping error Serr depth, a fit to the high wavenumber part of each spectrum is done, which provides the rms noise Srms. The number of averaged pings Npings into each superensemble varies as a function of depth, since the carousel was stopped more frequently to close Niskin bottles when getting to shallower depth during the upcast (Table 1). The single-ping error [Serr = Srms(Npings)1/2] increased with depth from 0.06 m s−1 above 1400 m to 0.11 m s−1 deeper than 3900 m (Table 1). This is consistent with the single-ping error profile provided by the LADCP software for each station. This single-ping error increase with depth was likely due to a decrease in the number of scatterers with depth. To be consistent with the variable observed noise levels, a depth-dependent wavenumber knoise was chosen in the variance computation (Table 1).

Fig. 4.
Fig. 4.

(a) Buoyancy-normalized shear spectra by depth range (solid) and corresponding GM spectra (dashed). The single-ping error in velocity (dashed–dotted) is given in the legend (cm s−1) for each depth range, taking into account the depth-dependent number of pings (see Table 1). Red triangles located at knoise give the upper limit of spectral integration for each depth range. (b) Strain spectra by depth range (thick lines) and corresponding GM spectra (thin lines). Vertical dashed lines denote the lower and maximum upper limits to determine the strain variances.

Citation: Journal of Atmospheric and Oceanic Technology 31, 11; 10.1175/JTECH-D-14-00036.1

Table 1.

Number of pings and single-ping error by depth range. Given the number of pings, the single-ping error was estimated as a fit to the corresponding shear spectrum (Fig. 4).

Table 1.

A SBE 911 CTD was mounted on the carousel. CTD data were processed to provide 1-m vertical resolution density profiles consistent with World Ocean Circulation Experiment (WOCE) standards. Strain defined as the vertical derivative of the isopycnal displacement ζ was estimated on the same vertical segments as the velocity shear. For a given downcast CTD profile, ζ were calculated from density fluctuations as , where ρHP is the density profile filtered by a high-pass Butterworth filter with a cutoff wavelength of 200 m, is the segment-averaged buoyancy frequency, ρ0 is a reference density, and g is the standard gravity. The filter cutoff frequency was chosen to avoid contamination by background stratification (Kunze et al. 2006). The use of a local filter does not make any assumption on the shape of the density profile on the segment. Fitting a polynomial (e.g., quadratic fit for Kunze et al. 2006; linear fit for Lauderdale et al. 2008) to remove the low-frequency part from the density profile can produce in the thermocline increased high-frequency density fluctuations if the degree of the polynomial is not consistent with the shape of the density profile. Strain variances were computed as the integral of the isopycnal displacement spectra S(ζ) from to , and , where = 2π/160 rad m−1 and is the minimum of 2π/10 rad m−1 and a cutoff wavenumber for which the integrated strain variance equals 0.22 (i.e., the variance of the GM strain spectra integrated up to its cutoff wavenumber of 2π/10 rad m−1). The GM strain spectra were also integrated from to . Ensemble-averaged strain spectra exhibit levels that are 40% (upper water column) to 100% (at depth) larger than their corresponding GM levels (Fig. 4b).

Since wavenumber ranges for shear and strain variances are not equal, the shear-to-strain ratio is calculated as , where both shear terms of the numerator are evaluated on the range [kmin; kmax], both strain terms of the denominator are evaluated on the range [kζmin; kζmax], and .

3. Results

a. Dissipation rates from microstructure

The distribution of dissipation rates was first inferred from microstructure velocity shear data (Fig. 5a) and averaged on the same grid as the finestructure estimates (160-m-long bins; εbin = ∑εi/Nbin, where εbin is the linear average of the Nbin individual microstructure dissipation rate estimates εi present in the bin). The weakest dissipation rate is found below 2000 m in the Iberian Abyssal Plain, to the north of the Azores–Biscay Rise and in the center of the Irminger basin with values mostly below 0.5 × 10−10 W kg−1. Above 2000 m, dissipation rates remain weak on the transect, usually around 1–5 × 10−10 W kg−1 with a few peaks in the range 10–50 × 10−10 W kg−1. The most intense values are found above the northwestern flank of the Reykjanes Ridge with values in the range of 5–100 × 10−10 W kg−1 and at 48°N, 21°W (blue arrow in Fig. 1) from 3400 m down to the seafloor at 4400 m with dissipation rates ranging from 10 to 100 × 10−10 W kg−1.

Fig. 5.
Fig. 5.

(a) Dissipation rate profiles εvmp (W kg−1) calculated from microstructure velocity shear data. (b) Corresponding diffusivity profiles Kvmp (m2 s−1). Bathymetry, as derived from CTD casts, is shaded in gray. Note the nonlinear color scales.

Citation: Journal of Atmospheric and Oceanic Technology 31, 11; 10.1175/JTECH-D-14-00036.1

To quantify the associated vertical diffusivity K, an equilibrium between production of turbulent kinetic energy by vertical velocity shear, increase in potential energy (mixing) and dissipation rate is often assumed leading to K = ΓεN−2 (Osborn 1980) with a constant mixing efficiency Γ = 0.2. Shih et al. (2005) and more recently Bouffard and Boegman (2013) pointed to a dependency of the diffusivity upon the turbulent intensity parameter Reb = ε/(νN2). Depending on the value of Reb, four regimes for the diffusivity were identified: molecular (Reb < 1.7), buoyancy-controlled (1.7 < Reb < 8.5), transitional (8.5 < Reb < 400) and energetic (Reb > 400). The boundaries between the various regimes are given for a Prandtl number of 7 and are supported by field data (Bouffard and Boegman 2013). In the transitional regime, the Osborn relation with Γ = 0.2 can be applied whereas in the energetic regime K = 4ν(Reb)1/2 is more appropriate. In the buoyancy-controlled (energetic) regime the use of the Osborn relation with Γ = 0.2 underestimates (overestimates) the diffusivity. In this study, Reb was derived from microstructure data (Fig. 6) and used to calculate the diffusivity. Most of the section belongs to the transitional regime for which the Osborn relation is valid. The energetic regime is found at a depth mainly at 21°W and on the northwestern flank of the Reykjanes Ridge.

Fig. 6.
Fig. 6.

Turbulence intensity parameter Reb = ε/(νN2). The color scale is nonlinear and emphasizes four mixing regimes (bounded by 1.7, 8.5, 400).

Citation: Journal of Atmospheric and Oceanic Technology 31, 11; 10.1175/JTECH-D-14-00036.1

The vertical diffusivities were usually below 0.5 × 10−4 m2 s−1 on the transect, with values well below 0.11 × 10−4 m2 s−1 in the upper 2000 m of the Iberian Abyssal Plain and at depth in the Irminger Basin (Fig. 5b). The two major regions with enhanced diffusivity concerned the bottom of the station at 21°W and 48°N where K reaches 5 × 10−4 m2 s−1 to 50 × 10−4 m2 s−1 and the flank of the Reykjanes Ridge with moderate values ranging from 1 to 5 × 10−4 m2 s−1. Note that for those two regions, using the Osborn relation with Γ = 0.2 as done in most of previous studies would lead to a diffusivity larger by a factor of 3–4 along the Reykjanes Ridge and by a factor 4–8 at 21°W and 48°N.

b. Comparison between finestructure and microstructure estimates of the dissipation rate

Using the procedure described in section 2b, finescale estimates εG89, εP95, and εG03 are calculated using vertical profiles of velocity and density. The comparison of those estimates with the microstructure dissipation rate εvmp allows us to quantify how the different terms entering the parameterizations improve the consistency with εvmp (Fig. 7). Since the dissipation rate varies on several orders of magnitude and in order to compare equally the under- and overestimation of εvmp, we defined the ratio Rε = log10(εfinestructure/εvmp), whose mean and standard deviation calculated for each of the tested parameterization is a measure of the bias between the parameterization and εvmp and of the dispersion of the estimates, respectively (Table 2).

Fig. 7.
Fig. 7.

Dissipation rate εfinestructure as a function of εvmp. (a)–(d) Color scale denotes log10 (N2) with εfinestructure= εG89, εP95, εG03, εG03, respectively. Solid line denotes perfect agreement with microstructure, dashed lines represent a factor of 3, and dotted lines represent a factor of 10.

Citation: Journal of Atmospheric and Oceanic Technology 31, 11; 10.1175/JTECH-D-14-00036.1

Table 2.

Averaged (〈 Rε〉) and standard deviation (σ) of the magnitude of the ratio between the finestructure and the microstructure estimates of ε with Rε = log10(εfinestructure/εmicrostructure). Columns F(m) gives the cumulative distribution P(1/m < εfinestructure/εmicrostructure < m) with P denoting the probability that εfinestructure/εmicrostructure lies in the interval (1/m, m). Dissipation ε*G89 (ε*ζ) is the same as εG89 (εζ) but restricted to the region where 〈ω〉 > 1.1f. εG03cons uses conservative fixed wavenumber bands to estimate shear and strain variances (see text for details).

Table 2.

The simplest parameterization, εG89, only agrees on average within an order of magnitude with εvmp (Fig. 7a) and clearly overestimates εvmp. The mean Rε of 0.17 (Table 2) corresponds to an averaged overestimation by 50% (100.17 = 1.5). Weakly stratified regions are mostly responsible for this overestimation (Fig. 7a, dark blue dots). The exclusion of these regions where the buoyancy frequency N2 falls below 10−7 s−2 improves the results and leads to an overestimation of 30%. The standard deviation of the ratio is high (0.52) when all segments are considered and slightly decreases to 0.44 when weakly stratified regions are omitted. Taking the frequency modulation term h(Rω) into account as done in the parameterization εP95 significantly improves the consistency of the parameterization with εvmp (Fig. 7b). The term h(Rω) is lower than unity over most of the section (Figs. 8a,b). This strongly reduces the positive bias observed with εG89. The lowest values of h(Rω) are found in weakly stratified regions where h(Rω) is close to 0.2. Most of the values larger than one are found in the upper 2000 m east of 21°W (distance to Portugal < 1500 km). The mean of Rε falls close to zero (−0.02) and its standard deviation decreases to a value of 0.39 with all segments taken into account. The added term j(f/N) in the parameterization εG03 does not have a large impact on our dataset. Gregg et al. (2003) showed its strong impact at low latitudes. But it has also a significant influence in the case of a weak stratification. Indeed, this term decreases the estimated dissipation rate by a factor of 2–3 in the weakly stratified regions (Fig. 8c) of the Iberian Abyssal Plains. Elsewhere on the section, the term j(f/N) tends to decrease the dissipation rates below 1500 m and to increase it above, with a maximum increase of 20%–30% in the stratified region of the North Atlantic Drift between 48°N, 21°W and the Reykjanes Ridge. Overall, the contribution of j(f/N) slightly reduces the scatter in εG03 versus εvmp (Fig. 7c). The mean of Rε is −0.05, slightly larger than for εP95, and its standard deviation is slightly improved with a value of 0.37 when all segments are taken into account (Table 2). Consistently, the cumulative distribution F(m) = P(1/m < εfinestructure/εvmp < m), with P being the probability that εfinestructure/εvmp lies in the interval [1/m, m], gives better results for εP95 than for εG89 (Table 2). The cumulative distributions of εP95/εvmp and εG03/εvmp are similar. The fact that the standard deviation is nonetheless better for εG03 than for εP95 is associated with the fewer extreme points in the εG03 parameterization. The probability density function (PDF) of the ratio εG03/εvmp is symmetric about unity (Fig. 9). For 3% of the segments, the ratio εG03/εvmp presents extreme values larger than 7 or weaker than (4% for the ratio εP95/εvmp) .

Fig. 8.
Fig. 8.

(a) Shear-to-strain ratio Rω, (b) frequency modulation term h(Rω), and (c) latitude/buoyancy frequency modulation term j(f/N) used in the finescale parameterizations εP95 and εG03.

Citation: Journal of Atmospheric and Oceanic Technology 31, 11; 10.1175/JTECH-D-14-00036.1

Fig. 9.
Fig. 9.

PDF of the ratio εG03/εvmp.

Citation: Journal of Atmospheric and Oceanic Technology 31, 11; 10.1175/JTECH-D-14-00036.1

The regional distribution of the ratio εG03/εvmp does not provide a clear picture (Fig. 10). Some specific features stand out: εG03 tends to overestimate εvmp in the upper 1000 m of the stratified waters south of the Reykjanes Ridge and to underestimate εvmp in the upper 1500 m of the weakly stratified waters of the Irminger basin; in the deepest part of the transect, εG03 overestimates εvmp in the poorly stratified waters of the Iberian Abyssal Plain (z > 4500 m) and in the deepest part of the Irminger Basin, while it underestimates εvmp from 21° to 26°W and on the northwestern flank of the Reykjanes Ridge.

Fig. 10.
Fig. 10.

Distribution of the ratio εG03/εvmp.

Citation: Journal of Atmospheric and Oceanic Technology 31, 11; 10.1175/JTECH-D-14-00036.1

To test the sensitivity of εG03 to the spectral ranges chosen to calculate the shear and strain variances, another estimate εG03cons is presented that also uses Gregg et al. (2003) finescale parameterization with “conservative” choices taken from Kunze et al. (2006): segments are 320 m long, strain (shear) variances are calculated using a [30–150] m ([150–320] m) wavenumber range. The new estimate εG03cons exhibits a larger dispersion than εG03, especially for observed dissipation rates εvmp < 1 × 10−10 W kg−1and buoyancy frequencies N < 1.3 × 10−3 s−1 (Fig. 7d; Table 2). A closer inspection of the results shows that most of the increase in the dispersion comes from the shear variance estimate being restricted to a shorter wavenumber range. Indeed, only two spectral components are used in the range [150–320] m to calculate the shear variance for εG03cons, which produces a larger uncertainty than for εG03. The moderate reduction in the range of wavenumbers for the strain variance only slightly increases the dispersion but not significantly (εG03cons dispersion σ increased by 0.02 compared with εG03).

c. Dissipation rates and diffusivity from finestructure on the full transect

To get a complete view of the distribution of the dissipation rate on the OVIDE transect, the parameterization εG03 is applied to all the CTD–LADCP casts and the diffusivity is derived following Bouffard and Boegman (2013; Fig. 11). For buoyancy frequencies weaker than 1.4 s−1 (potential temperature θ < 4.5°C), the dissipation rate usually reaches weak values, εG03 < 10−10 W kg−1. High values of the dissipation rate are found above the bottom at some locations as the microstructure data suggested. The southeastern flank of the Reykjanes Ridge, for which no microstructure data were available because of a heavy sea state, shows the same moderate intensification as that found on the northwestern flank. A moderate intensification is also visible along the Greenland slope.

Fig. 11.
Fig. 11.

(a) Dissipation rate profiles εG03 (W kg−1) calculated from the finestructure of CTD–LADCP data. (b) Corresponding diffusivity profiles KG03 (m2 s−1). Bathymetry, as derived from CTD casts, is shaded in gray. The nonlinear color scales are the same as in Fig. 5.

Citation: Journal of Atmospheric and Oceanic Technology 31, 11; 10.1175/JTECH-D-14-00036.1

Since the diffusivity scales as εG03N−2, the diffusivity distribution is somehow opposite to that of the dissipation. The largest values are found at a depth where the stratification is weak. The most intense values are found in a 1000-m-thick layer above the seafloor from the Iberian Plain up to 21°W (from 1 to 10 × 10−4 m2 s−1), around the Reykjanes Ridge (1–5 × 10−4 m2 s−1) and above 1200 m on the Greenland slope (1–5 × 10−4 m2 s−1). Elsewhere on the transect, the diffusivity falls below 5 × 10−5 m2 s−1.

In regions of weakly stratified waters or when only CTD data are available, Kunze et al. (2006) suggested the use of a strain-based parameterization εζ as presented in section 2. The application of such a parameterization with a constant shear-to-strain ratio representative of the section-averaged Rω gives contrasted results. From 10° to 23°W, the estimated εζ below a depth of 2500–3000 m strongly underestimates εvmp by a factor of 3–1000 (Fig. 12). The factor 1000 is reached for the profile at 21°W, where the microstructure showed enhanced levels above the bottom. For comparison, εG03 underestimated εvmp by a moderate factor of 3–5 at the bottom of this profile. Parameterization εζ completely misses the intensification of the dissipation rate at 21°W (cf. Figs. 5a, 12). Note that the region between 10° and 23°W is weakly stratified at depth, which hinders good estimates of the strain from one instantaneous vertical density profile. Moreover, the strain-based parameterization is expected to underestimate dissipation rates there, since this region is biased toward low frequencies (i.e., high values of Rω; Fig. 8a), which means weak signatures on the strain. In contrast, εζ overestimates εvmp by a factor of 2–10 above 1500–2000 m and south of the Reykjanes Ridge (Fig. 12b), where Rω is relatively low (region biased toward higher frequencies). Again, εG03 provides a better estimate despite a similar bias in this region. North of the Reykjanes Ridge εζ and εG03 are similar. The averaged ratio log10(εζ/εvmp) is −0.12 with a standard deviation of 0.48 (Table 2). Those statistics are close to that obtained with εG89 that similarly assumes a constant shear-to-strain ratio [h(Rω) = 1 <=> Rω = RωGM = 3] together with j(f/N) = 1. However, εG89 produces reasonable estimates where εvmp was enhanced at depth. From this comparison, we conclude that the use of εζ based on a single density profile with a constant Rω has to be used with care and may miss important features and produce regionally strong biases even when a regional value of Rω is chosen.

Fig. 12.
Fig. 12.

(a) Dissipation rate profiles εζ (W kg−1) calculated from the finestructure of CTD data with a constant shear-to-strain ratio of 5.2. (b) Distribution of the ratio εζ/εvmp.

Citation: Journal of Atmospheric and Oceanic Technology 31, 11; 10.1175/JTECH-D-14-00036.1

To gain more insight into the frequency content of the internal wave field, an estimate of the wave-field-averaged frequency 〈ω〉 is diagnosed from Rω [see Polzin et al. 1995, Eq. (15)]. For the whole transect, 〈ω〉 is close to the inertial frequency (between f and 2f) apart from a small region located between 14° and 16°N and at 1200–2200-m depth for which 〈ω〉 is nearly 3 times as large as f (Fig. 13a). Hibiya et al. (2012) used another approach to estimate the relative contribution of α, the “high” (2f-N) frequency band, and of β, the “low” (f-2f) frequency band, to the observed shear and strain variances. The ratio β/α measures the spectral distortion from the GM spectra. Applied to our dataset, we found that the high-frequency band is close to the GM level (α < 2) apart from the southeastern flank of the Reykjanes Ridge (6 < α < 10) and in a small region centered at 15°N and 1700 dbar (4 < α < 7). The contribution of the low-frequency band is usually close to the GM level (β < 2) above 2000 dbar. The largest values (β > 6) are found at depth below 3500 dbar to the south of 26°N and on the northwestern flank of the Reykjanes Ridge. Elsewhere, the low-frequency band is only moderately higher than the GM prescription (2 < β < 5). Interestingly, the ratio β/α shows a good consistency with 〈ω〉 (Fig. 13b). The high-frequency band tends to dominate (β/α < ⅓) in a few patches above 2000 m located south of 22°N. At depth, the ratio β/α is also smaller than ⅓ on the southern flank of the Reykjanes Ridge but not on its northern flank. Elsewhere, the low-frequency band mostly dominates (β/α > 4). Note that 〈ω〉 and β/α are indirect estimates of the averaged frequency content of the wave field subject to uncertainties. Several criteria based on either α, β, β/α, or 〈ω〉 have been tested to define a validity threshold above which εζ remains consistent with εvmp, thus avoiding the large discrepancy found at depth in the weakly stratified waters. The only acceptable criterion found is to consider the use of εζ only when 〈ω〉 > 1.1f, with f being the Coriolis parameter (Fig. 13a, white line). Using this criterion, the statistics of the ratio εζ/εvmp are not significantly different from those of εG03/εvmp (Table 2, line ε*ζ). This confirms that the sole isopycnal strain data diagnosed from the CTD cannot be used to quantify dissipation rates in a near-inertial wave field with 〈ω〉 < 1.1f and in the presence of a weakly stratified environment. No criterion was found that improves the statistics of the ratio εG89/εvmp to the level of εG03/εvmp, although restricting the calculation to either 〈ω〉 > 1.1f or β < 10 improves by 4% the cumulative distribution while still showing a mean positive bias (+25%)(Table 2, line ε*G89).

Fig. 13.
Fig. 13.

(a) Ratio 〈ω〉/f of the averaged internal wave frequency scaled by the inertial frequency; the white line denotes the threshold 〈ω〉 = 1.1f. (b) Ratio β/α of the low-to-high wavenumber band contribution to the observed shear and strain variances.

Citation: Journal of Atmospheric and Oceanic Technology 31, 11; 10.1175/JTECH-D-14-00036.1

4. Summary and discussion

For our dataset, the largest overestimation of εvmp by the finescale estimate appears in the Iberian Abyssal Plain below a depth of 4400 m (Fig. 10). In this region of weak stratification (N ≈ 3f), both the terms h(Rω) and j(f/N) contribute to improve the consistency of εG03 with εvmp (see the evolution of the dark blue dots from Fig. 7a to Fig. 7c). However, some of the segments remain strongly inconsistent (by a factor of 10). From an analysis of the shear-to-strain ratio Rω, Kunze et al. (2006) argued that weakly stratified waters are associated with velocity profiles that are potentially strongly contaminated by noise, since the number of scatterers there is often low. They reported very high values of Rω ranging from 30 to more than 1000 for N < 4.5 × 10−4 s−1 (see their . Below 4000 m in the Iberian Plain, Rω also exhibits values ranging from 30 to 300 for a buoyancy frequency N < 3.2 × 10−4 s−1 (Fig. 8a). The noise associated with those deep profiles (z > 3900 m) is more intense by a factor of 2 compared to depths z < 2700 m and is clearly visible for wavenumbers larger than 0.07 rad m−1 (Table 1; Fig. 4). At low wavenumbers (k < 0.05 rad m−1), however, the spectral shape of the buoyancy-normalized shear is flat as it is at shallower depths (Fig. 4). This means that if noise exists on velocity profiles for k < 0.05 rad m−1, then it coincidentally has the same k−2 dependence before changing to white noise for k > 0.07 rad m−1. Although one cannot exclude that noise is responsible for the overestimation of the dissipation rate in those regions, one could also argue that the scaling of observed shear by the traditional GM shear may not be appropriate for weakly stratified waters. Indeed, Gerkema and Shrira (2005) and Gerkema and Exarchou (2008) showed that nontraditional effects (effects of the Coriolis terms proportional to the cosine of latitude) can significantly change internal wave properties in weakly stratified waters.

Shallower than 3500 m, the averaged transect value of Rω is 5.2 (standard deviation of 3), nearly twice the GM value of 3. Thus, there is on average more energy in near-inertial frequencies in our dataset than there is in the GM model (there is no peak of energy at tidal frequencies in the GM model). Kunze et al. (2006) found an averaged Rω of 7 for various sections distributed over the Atlantic, Pacific and Indian Oceans. The ratio of the low (f-2f) to high (2f-N) wavenumber band contribution to the observed shear and strain variances confirms that at depth spectra are biased toward low frequencies compared to GM spectra. This is somewhat expected, since GM spectra does not take the tidal forcing into account. Interestingly the Reykjanes Ridge shows a significant cross-ridge contrast, the northwestern flank being biased toward low frequencies, while the southeastern flank is biased toward high frequencies. Over the Emperor seamount chain, Hibiya et al. (2012) also found some regions where β and α dominate at depth (e.g., their profiles H3 and H1, respectively). They attributed the bias toward high frequencies to tidal interactions with topographic features, which also possibly occur on the southeastern flank of the Reykjanes Ridge. Additional temporal measurements from moorings are needed to clearly quantify whether the contrast as diagnosed from β/α and 〈ω〉 is really significant. One can, however, note from altimetric data that there is also a cross-ridge asymmetry in the mesoscale eddy kinetic energy distribution with much higher levels to the northwest than to the southeast of the ridge. Nikurashin and Ferrari (2010) report enhanced near-inertial motions due to lee-wave generation by geostrophic currents, a process that is possibly at play to the north of the ridge.

On most parts of the section, the finescale parameterization εG03 is consistent with the microstructure dissipation rate εvmp. This suggests that most of the turbulent activity recorded in this dataset is due to the internal wave field. Over 75% of the transect, the agreement is within a factor of 2.5 and almost 95% within a factor 4.5. Although we did not average the segments into bins as done in Polzin et al. (1995), which would reduce the scatter, those statistics are consistent with previous published studies (e.g., Polzin et al. 1995, 2002; Gregg et al. 2003). The use of the frequency-based correction term h(Rω) improves the consistency as Polzin et al. (1995) found. The latitude/buoyancy frequency term slightly improves the parameterization, especially for weakly stratified waters. On average, the parameterization εG03 tends to underestimate the microstructure dissipation rate by 11%. Most of the transect takes place away from internal wave generation sites (away from the surface and away from rough topography). The environment is then favorable for the comparison with finescale parameterizations based on the GM model, which is representative of the internal wave field away from boundaries.

Apart from a few locations near the topography, diffusivity remains weak. It is lower than 1 × 10−5 m2 s−1 over 44% of the transect and lower than 5 × 10−5 m2 s−1 over 90% of the transect. The picture emerging from this dataset is consistent with the finescale estimate of the dissipation rate along the A16N WOCE hydrographic line in the eastern North Atlantic for the same latitude range (Kunze et al. 2006, their Fig. 11). Stöber et al. (2008) reported frequent occurrences of enhanced diffusivities larger than 1 × 10−4 m2 s−1 on the Greenland slope off Cape Farewell (Eirik Ridge) for profiles located over isobaths ranging from 2700 up to 1100 m (their Fig. 4e). The part of our transect above the Greenland slope is located 100 km eastward and only exhibits such enhanced diffusivities for slopes shallower than 1150 m despite similar velocity amplitudes of the boundary current (maxima of 0.2–0.3 m s−1), which extends down to the isobath 2500 m. Consistently, Lauderdale et al. (2008) reported bottom-intensified diffusivities frequently beyond 10 × 10−4 m2 s−1 around the Eirik Ridge and a factor of 5–10 less for a transect located on the Greenland slope near ours. Note that Lauderdale et al. (2008) used a strain-based parameterization with a constant “conservative” shear-to-strain ratio of 7.6 [hζ (Rω) = 3 in Eq. (4)]. At 48°N and from 22° to 14°W, Walter et al. (2005) applied the shear variance [Eq. (3)] to a set of LADCP and CTD profiles to calculate dissipation and diffusivity. They found enhanced diffusivities that are systematically in the range of 10–100 × 10−4 m2 s−1 near the bottom and that they gradually decreased to 1 × 10−4 m2 s−1 at a distance of more than 1000 m above the seafloor. Such a strong intensification widely distributed in the eastern basin was partly attributed to the roughness of the bathymetry due to the proximity of the Mid-Atlantic Ridge. They also pointed out that from 40° to 50° of latitude, subinertial wave energy can be trapped (Gerkema and Shrira 2005) and lead to such intense diffusivities. In the region of our transect located in the same eastern basin, we indeed observe an intensification of turbulent diffusivities but 1) it is one to two orders of magnitude less important near the bottom (at most 1–5 × 10−4 m2 s−1; Fig. 11a) than Walter et al. (2005) found, and 2) it is not associated with an intensification of the dissipation rate (Fig. 11b) but results from the weak values of the buoyancy frequency N used to calculated diffusivities from the Osborn (1980) relation. East of 21°W, our microstructure data do not show an intensification of the dissipation rate at depth either. Consistent with our study, Kunze et al. (2006) only find a few spots of elevated diffusivity that hardly reach 1 × 10−4 m2 s−1i in the same deep basin. We believe that part of the discrepancy observed in this region at depth comes from the fact that Walter et al. (2005) did not use the scaling term function of the buoyancy frequency {cosh−1(N/f)/[cosh−1(N0/f30)]} in j(f/N) [Eq. (3)]. Our dataset showed that for weakly stratified waters, this term decreased the finescale estimate of the dissipation rates and diffusivity, and improved the consistency with microstructure observations of the dissipation rates. Thus, our microstructure data and our finestructure shear/strain dataset do not support the idea of a basinwide intensification of the turbulent diffusivities at depth because of wave trapping between 40° and 50°N. The canonical globally averaged value of order 10−4 m2 s−1 needed to maintain the abyssal stratification (Munk 1966) is only reached on the flank of the Reykjanes Ridge and in the region around 21°W. Enhanced dissipation rates and diffusivity above the Reykjanes Ridge are potentially sustained by enhanced levels of semidiurnal tidal energy observed above the Reykjanes Ridge (van Haren 2007) and by breaking of lee waves generated by bottom geostrophic flows above this rough topography (Nikurashin and Ferrari 2011).

Acknowledgments

The microstructure profiler was funded by the French Agence Nationale de la Recherche (ANR) through Grant ANR-JC05_50690 and by the French Institute for Marine Science (IFREMER). B. Ferron and H. Mercier were supported by the French National Center for Scientific Research (CNRS), F. Kokoszka by a doctoral grant from the Ministère de la Rechercher et de l'Enseignement Supérieur, and P. Lherminier by IFREMER.

REFERENCES

  • Bouffard, D., and Boegman L. , 2013: A diapycnal diffusivity model for stratified environmental flows. Dyn. Atmos. Oceans, 61–62, 1434, doi:10.1016/j.dynatmoce.2013.02.002.

    • Search Google Scholar
    • Export Citation
  • Cairns, J. L., and Williams G. O. , 1976: Internal wave observations from a midwater float, 2. J. Geophys. Res., 81, 19431950, doi:10.1029/JC081i012p01943.

    • Search Google Scholar
    • Export Citation
  • Daae, K. L., Fer I. , and Abrahamsen E. P. , 2009: Mixing on the continental slope of the southern Weddell Sea. J. Geophys. Res., 114, C09018, doi:10.1029/2008JC005259.

    • Search Google Scholar
    • Export Citation
  • Dillon, T. M., 1982: Vertical overturns: A comparison of Thorpe and Ozmidov length scales. J. Geophys. Res., 87, 96019613, doi:10.1029/JC087iC12p09601.

    • Search Google Scholar
    • Export Citation
  • Firing, E., and Gordon R. , 1990: Deep ocean acoustic Doppler current profiling. Proceedings of the IEEE Fourth Working Conference on Current Measurement, IEEE, 192–201, doi:10.1109/CURM.1990.110905.

  • Fischer, J., and Visbeck M. , 1993: Deep velocity profiling with self-contained ADCPs. J. Atmos. Oceanic Technol., 10, 764773, doi:10.1175/1520-0426(1993)010<0764:DVPWSC>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Frants, M., Damerell G. M. , Gille S. T. , Heywood K. J. , MacKinnon J. , and Sprintall J. , 2013: An assessment of density-based finescale methods for estimating diapycnal diffusivity in the Southern Ocean. J. Atmos. Oceanic Technol., 30, 26472661, doi:10.1175/JTECH-D-12-00241.1.

    • Search Google Scholar
    • Export Citation
  • Garrett, C., and Munk W. , 1975: Space–time scales of internal waves: A progress report. J. Geophys. Res., 80, 291297, doi:10.1029/JC080i003p00291.

    • Search Google Scholar
    • Export Citation
  • Gerkema, T., and Shrira V. I. , 2005: Near-inertial waves on the ‘‘non-traditional’’ β plane. J. Geophys. Res., 110, C01003, doi:10.1029/2004JC002519.

    • Search Google Scholar
    • Export Citation
  • Gerkema, T., and Exarchou E. , 2008: Internal-wave properties in weakly stratified layers. J. Mar. Res., 66, 617644, doi:10.1357/002224008787536817.

    • Search Google Scholar
    • Export Citation
  • Goodman, L., Levine E. R. , and Lueck R. G. , 2006: On measuring the terms of the turbulent kinetic energy budget from an AUV. J. Atmos. Oceanic Technol., 23, 977990, doi:10.1175/JTECH1889.1.

    • Search Google Scholar
    • Export Citation
  • Gregg, M. C., 1989: Scaling turbulent dissipation in the thermocline. J. Geophys. Res., 94, 96869698, doi:10.1029/JC094iC07p09686.

  • Gregg, M. C., 1999: Uncertainties and limitations in measuring ϵ and χT. J. Atmos. Oceanic Technol., 16, 14831490, doi:10.1175/1520-0426(1999)016<1483:UALIMA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Gregg, M. C., Sanford T. B. , and Winkel D. P. , 2003: Reduced mixing from the breaking of internal waves in equatorial ocean waters. Nature, 422, 513515, doi:10.1038/nature01507.

    • Search Google Scholar
    • Export Citation
  • Hibiya, T., and Nagasawa M. , 2004: Latitudinal dependence of diapycnal diffusivity in the thermocline estimated using a finescale parameterization. Geophys. Res. Lett., 31, L01301, doi:10.1029/2003GL017998.

    • Search Google Scholar
    • Export Citation
  • Hibiya, T., Furuichi N. , and Robertson R. , 2012: Assessment of fine-scale parameterizations of turbulent dissipation rates near mixing hotspots in the deep ocean. Geophys. Res. Lett., 39, L24601, doi:10.1029/2012GL054068.

    • Search Google Scholar
    • Export Citation
  • Hogg, N., Biscaye P. , Gardner W. , and Schmitz W. J. Jr., 1982: On transport and modification of Antarctic Bottom Water in the Vema Channel. J. Mar. Res., 40 (Suppl.), 231263.

    • Search Google Scholar
    • Export Citation
  • Kunze, E., Firing E. , Hummon J. M. , Chereskin T. K. , and Thurnherr A. M. , 2006: Global abyssal mixing inferred from lowered ADCP shear and CTD strain profiles. J. Phys. Oceanogr., 36, 15531576, doi:10.1175/JPO2926.1.

    • Search Google Scholar
    • Export Citation
  • Lauderdale, J. M., Bacon S. , Naveira Garabato A. C. , and Holiday N. P. , 2008: Intensified turbulent mixing in the boundary current system of southern Greenland. Geophys. Res. Lett., 35, L04611, doi:10.1029/2007GL032785.

    • Search Google Scholar
    • Export Citation
  • Ledwell, J. R., Watson A. J. , and Law C. S. , 1993: Evidence for slow mixing across the pycnocline from an open-ocean tracer release. Nature, 364, 701702, doi:10.1038/364701a0.

    • Search Google Scholar
    • Export Citation
  • Ledwell, J. R., Montgomery E. T. , Polzin K. L. , St. Laurent L. C. , Schmitt R. W. , and Toole J. M. , 2000: Evidence for enhanced mixing over rough topography in the abyssal ocean. Nature, 403, 179182, doi:10.1038/35003164.

    • Search Google Scholar
    • Export Citation
  • Levine, E. R., and Lueck R. G. , 1999: Turbulence measurement from an autonomous underwater vehicle. J. Atmos. Oceanic Technol., 16, 15331544, doi:10.1175/1520-0426(1999)016<1533:TMFAAU>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • MacKinnon, J. A., Johnston T. M. S. , and Pinkel R. , 2008: Strong transport and mixing of deep water through the Southwest Indian Ridge. Nat. Geosci., 1, 755758, doi:10.1038/ngeo340.

    • Search Google Scholar
    • Export Citation
  • Macoun, P., and Lueck R. , 2004: Modeling the spatial response of the airfoil shear probe using different sized probes. J. Atmos. Oceanic Technol., 21, 284297, doi:10.1175/1520-0426(2004)021<0284:MTSROT>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Mercier, H., and Coauthors, 2014: Variability of the meridional overturning circulation at the Greenland–Portugal OVIDE section from 1993 to 2010. Prog. Oceanogr., doi:10.1016/j.pocean.2013.11.001, in press.

    • Search Google Scholar
    • Export Citation
  • Moum, J. N., Gregg M. C. , Lien R. C. , and Carr M. E. , 1995: Comparison of turbulence kinetic energy dissipation rate estimates from two ocean microstructure profilers. J. Atmos. Oceanic Technol., 12, 346366, doi:10.1175/1520-0426(1995)012<0346:COTKED>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Munk, W. H., 1966: Abyssal recipes. Deep-Sea Res. Oceanogr. Abstr.,13, 707730, doi:10.1016/0011-7471(66)90602-4.

  • Munk, W. H., and Wunsch C. , 1998: Abyssal recipes II: Energetics of tidal and wind mixing. Deep-Sea Res. I, 45, 19772010, doi:10.1016/S0967-0637(98)00070-3.

    • Search Google Scholar
    • Export Citation
  • Nasmyth, P. V., 1970: Ocean turbulence. Ph.D. thesis, University of British Columbia, 69 pp.

  • Naveira Garabato, A. C., Oliver K. I. C. , Watson A. J. , and Messias M.-J. , 2004a: Turbulent diapycnal mixing in the Nordic seas. J. Geophys. Res., 109, C12010, doi:10.1029/2004JC002411.

    • Search Google Scholar
    • Export Citation
  • Naveira Garabato, A. C., Polzin K. L. , King B. A. , Heywood K. J. , and Visbeck M. , 2004b: Widespread intense turbulent mixing in the Southern Ocean. Science, 303, 210213, doi:10.1126/science.1090929.

    • Search Google Scholar
    • Export Citation
  • Nikurashin, M., and Ferrari R. , 2010: Radiation and dissipation of internal waves generated by geostrophic motions impinging on small-scale topography: Application to the Southern Ocean. J. Phys. Oceanogr., 40, 20252042, doi:10.1175/2010JPO4315.1.

    • Search Google Scholar
    • Export Citation
  • Nikurashin, M., and Ferrari R. , 2011: Global energy conversion rate from geostrophic flows into internal lee waves in the deep ocean. Geophys. Res. Lett., 38, L08610, doi:10.1029/2011GL046576.

    • Search Google Scholar
    • Export Citation
  • Oakey, N. S., 1982: Determination of the rate of dissipation of turbulent energy from simultaneous temperature and velocity shear microstructure measurements. J. Phys. Oceanogr., 12, 256271, doi:10.1175/1520-0485(1982)012<0256:DOTROD>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Osborn, T. R., 1980: Estimates of the local rate of vertical diffusion from dissipation measurements. J. Phys. Oceanogr., 10, 8389, doi:10.1175/1520-0485(1980)010<0083:EOTLRO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Polzin, K. L., Toole J. M. , and Schmitt R. W. , 1995: Finescale parameterization of turbulent dissipation. J. Phys. Oceanogr., 25, 306328, doi:10.1175/1520-0485(1995)025<0306:FPOTD>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Polzin, K. L., Kunze E. , Hummon J. , and Firing E. , 2002: The finescale response of lowered ADCP velocity profiles. J. Atmos. Oceanic Technol., 19, 205224, doi:10.1175/1520-0426(2002)019<0205:TFROLA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Sheen, K. L., and Coauthors, 2013: Rates and mechanisms of turbulent dissipation and mixing in the Southern Ocean: Results from the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES). J. Geophys. Res. Oceans, 118, 27742792, doi:10.1002/jgrc.20217.

    • Search Google Scholar
    • Export Citation
  • Shih, L. H., Koseff J. R. , Ivey G. N. , and Ferziger J. H. , 2005: Parameterization of turbulent fluxes and scales using homogeneous sheared stably stratified turbulence simulations. J. Fluid Mech., 525, 193214, doi:10.1017/S0022112004002587.

    • Search Google Scholar
    • Export Citation
  • Smith, W. H. F., and Sandwell D. T. , 1997: Global seafloor topography from satellite altimetry and ship depth soundings. Science, 277, 19561962, doi:10.1126/science.277.5334.1956.

    • Search Google Scholar
    • Export Citation
  • Stöber, U., Walter M. , Mertens C. , and Rhein M. , 2008: Mixing estimates from hydrographic measurements in the Deep Western Boundary Current of the North Atlantic. Deep-Sea Res. I, 55, 721736, doi:10.1016/j.dsr.2008.03.006.

    • Search Google Scholar
    • Export Citation
  • Thorpe, S. A., 1977: Turbulence and mixing in a Scottish loch. Philos. Trans. Roy. Soc. London, A286, 125181, doi:10.1098/rsta.1977.0112.

    • Search Google Scholar
    • Export Citation
  • Thurnherr, A. M., 2012: The finescale response of lowered ADCP velocity measurements processed with different methods. J. Atmos. Oceanic Technol., 29, 597600, doi:10.1175/JTECH-D-11-00158.1.

    • Search Google Scholar
    • Export Citation
  • Tian, J., Yang Q. , and Zhao W. , 2009: Enhanced diapycnal mixing in the South China Sea. J. Phys. Oceanogr., 39, 31913203, doi:10.1175/2009JPO3899.1.

    • Search Google Scholar
    • Export Citation
  • van Haren, H., 2007: Inertial and tidal shear variability above Reykjanes Ridge. Deep-Sea Res. I, 54, 856870, doi:10.1016/j.dsr.2007.03.003.

    • Search Google Scholar
    • Export Citation
  • Visbeck, M., 2002: Deep velocity profiling using lowered acoustic Doppler current profilers: Bottom track and inverse solutions. J. Atmos. Oceanic Technol., 19, 794807, doi:10.1175/1520-0426(2002)019<0794:DVPULA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Walter, M., Mertens C. , and Rhein M. , 2005: Mixing estimates from a large-scale hydrographic survey in the North Atlantic. Geophys. Res. Lett.,32, L13605, doi:10.1029/2005GL022471.

  • Waterman, S., Naveira Garabato A. C. , and Polzin K. L. , 2013: Internal waves and turbulence in the Antarctic Circumpolar Current. J. Phys. Oceanogr., 43, 259282, doi:10.1175/JPO-D-11-0194.1.

    • Search Google Scholar
    • Export Citation
  • Wijesekera, H., Padman L. , Dillon T. , Levine M. , Paulson C. , and Pinkel R. , 1993: The application of internal-wave dissipation models to a region of strong forcing. J. Phys. Oceanogr., 23, 269286, doi:10.1175/1520-0485(1993)023<0269:TAOIWD>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
Save
  • Bouffard, D., and Boegman L. , 2013: A diapycnal diffusivity model for stratified environmental flows. Dyn. Atmos. Oceans, 61–62, 1434, doi:10.1016/j.dynatmoce.2013.02.002.

    • Search Google Scholar
    • Export Citation
  • Cairns, J. L., and Williams G. O. , 1976: Internal wave observations from a midwater float, 2. J. Geophys. Res., 81, 19431950, doi:10.1029/JC081i012p01943.

    • Search Google Scholar
    • Export Citation
  • Daae, K. L., Fer I. , and Abrahamsen E. P. , 2009: Mixing on the continental slope of the southern Weddell Sea. J. Geophys. Res., 114, C09018, doi:10.1029/2008JC005259.

    • Search Google Scholar
    • Export Citation
  • Dillon, T. M., 1982: Vertical overturns: A comparison of Thorpe and Ozmidov length scales. J. Geophys. Res., 87, 96019613, doi:10.1029/JC087iC12p09601.

    • Search Google Scholar
    • Export Citation
  • Firing, E., and Gordon R. , 1990: Deep ocean acoustic Doppler current profiling. Proceedings of the IEEE Fourth Working Conference on Current Measurement, IEEE, 192–201, doi:10.1109/CURM.1990.110905.

  • Fischer, J., and Visbeck M. , 1993: Deep velocity profiling with self-contained ADCPs. J. Atmos. Oceanic Technol., 10, 764773, doi:10.1175/1520-0426(1993)010<0764:DVPWSC>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Frants, M., Damerell G. M. , Gille S. T. , Heywood K. J. , MacKinnon J. , and Sprintall J. , 2013: An assessment of density-based finescale methods for estimating diapycnal diffusivity in the Southern Ocean. J. Atmos. Oceanic Technol., 30, 26472661, doi:10.1175/JTECH-D-12-00241.1.

    • Search Google Scholar
    • Export Citation
  • Garrett, C., and Munk W. , 1975: Space–time scales of internal waves: A progress report. J. Geophys. Res., 80, 291297, doi:10.1029/JC080i003p00291.

    • Search Google Scholar
    • Export Citation
  • Gerkema, T., and Shrira V. I. , 2005: Near-inertial waves on the ‘‘non-traditional’’ β plane. J. Geophys. Res., 110, C01003, doi:10.1029/2004JC002519.

    • Search Google Scholar
    • Export Citation
  • Gerkema, T., and Exarchou E. , 2008: Internal-wave properties in weakly stratified layers. J. Mar. Res., 66, 617644, doi:10.1357/002224008787536817.

    • Search Google Scholar
    • Export Citation
  • Goodman, L., Levine E. R. , and Lueck R. G. , 2006: On measuring the terms of the turbulent kinetic energy budget from an AUV. J. Atmos. Oceanic Technol., 23, 977990, doi:10.1175/JTECH1889.1.

    • Search Google Scholar
    • Export Citation
  • Gregg, M. C., 1989: Scaling turbulent dissipation in the thermocline. J. Geophys. Res., 94, 96869698, doi:10.1029/JC094iC07p09686.

  • Gregg, M. C., 1999: Uncertainties and limitations in measuring ϵ and χT. J. Atmos. Oceanic Technol., 16, 14831490, doi:10.1175/1520-0426(1999)016<1483:UALIMA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Gregg, M. C., Sanford T. B. , and Winkel D. P. , 2003: Reduced mixing from the breaking of internal waves in equatorial ocean waters. Nature, 422, 513515, doi:10.1038/nature01507.

    • Search Google Scholar
    • Export Citation
  • Hibiya, T., and Nagasawa M. , 2004: Latitudinal dependence of diapycnal diffusivity in the thermocline estimated using a finescale parameterization. Geophys. Res. Lett., 31, L01301, doi:10.1029/2003GL017998.

    • Search Google Scholar
    • Export Citation
  • Hibiya, T., Furuichi N. , and Robertson R. , 2012: Assessment of fine-scale parameterizations of turbulent dissipation rates near mixing hotspots in the deep ocean. Geophys. Res. Lett., 39, L24601, doi:10.1029/2012GL054068.

    • Search Google Scholar
    • Export Citation
  • Hogg, N., Biscaye P. , Gardner W. , and Schmitz W. J. Jr., 1982: On transport and modification of Antarctic Bottom Water in the Vema Channel. J. Mar. Res., 40 (Suppl.), 231263.

    • Search Google Scholar
    • Export Citation
  • Kunze, E., Firing E. , Hummon J. M. , Chereskin T. K. , and Thurnherr A. M. , 2006: Global abyssal mixing inferred from lowered ADCP shear and CTD strain profiles. J. Phys. Oceanogr., 36, 15531576, doi:10.1175/JPO2926.1.

    • Search Google Scholar
    • Export Citation
  • Lauderdale, J. M., Bacon S. , Naveira Garabato A. C. , and Holiday N. P. , 2008: Intensified turbulent mixing in the boundary current system of southern Greenland. Geophys. Res. Lett., 35, L04611, doi:10.1029/2007GL032785.

    • Search Google Scholar
    • Export Citation
  • Ledwell, J. R., Watson A. J. , and Law C. S. , 1993: Evidence for slow mixing across the pycnocline from an open-ocean tracer release. Nature, 364, 701702, doi:10.1038/364701a0.

    • Search Google Scholar
    • Export Citation
  • Ledwell, J. R., Montgomery E. T. , Polzin K. L. , St. Laurent L. C. , Schmitt R. W. , and Toole J. M. , 2000: Evidence for enhanced mixing over rough topography in the abyssal ocean. Nature, 403, 179182, doi:10.1038/35003164.

    • Search Google Scholar
    • Export Citation
  • Levine, E. R., and Lueck R. G. , 1999: Turbulence measurement from an autonomous underwater vehicle. J. Atmos. Oceanic Technol., 16, 15331544, doi:10.1175/1520-0426(1999)016<1533:TMFAAU>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • MacKinnon, J. A., Johnston T. M. S. , and Pinkel R. , 2008: Strong transport and mixing of deep water through the Southwest Indian Ridge. Nat. Geosci., 1, 755758, doi:10.1038/ngeo340.

    • Search Google Scholar
    • Export Citation
  • Macoun, P., and Lueck R. , 2004: Modeling the spatial response of the airfoil shear probe using different sized probes. J. Atmos. Oceanic Technol., 21, 284297, doi:10.1175/1520-0426(2004)021<0284:MTSROT>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Mercier, H., and Coauthors, 2014: Variability of the meridional overturning circulation at the Greenland–Portugal OVIDE section from 1993 to 2010. Prog. Oceanogr., doi:10.1016/j.pocean.2013.11.001, in press.

    • Search Google Scholar
    • Export Citation
  • Moum, J. N., Gregg M. C. , Lien R. C. , and Carr M. E. , 1995: Comparison of turbulence kinetic energy dissipation rate estimates from two ocean microstructure profilers. J. Atmos. Oceanic Technol., 12, 346366, doi:10.1175/1520-0426(1995)012<0346:COTKED>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Munk, W. H., 1966: Abyssal recipes. Deep-Sea Res. Oceanogr. Abstr.,13, 707730, doi:10.1016/0011-7471(66)90602-4.

  • Munk, W. H., and Wunsch C. , 1998: Abyssal recipes II: Energetics of tidal and wind mixing. Deep-Sea Res. I, 45, 19772010, doi:10.1016/S0967-0637(98)00070-3.

    • Search Google Scholar
    • Export Citation
  • Nasmyth, P. V., 1970: Ocean turbulence. Ph.D. thesis, University of British Columbia, 69 pp.

  • Naveira Garabato, A. C., Oliver K. I. C. , Watson A. J. , and Messias M.-J. , 2004a: Turbulent diapycnal mixing in the Nordic seas. J. Geophys. Res., 109, C12010, doi:10.1029/2004JC002411.

    • Search Google Scholar
    • Export Citation
  • Naveira Garabato, A. C., Polzin K. L. , King B. A. , Heywood K. J. , and Visbeck M. , 2004b: Widespread intense turbulent mixing in the Southern Ocean. Science, 303, 210213, doi:10.1126/science.1090929.

    • Search Google Scholar
    • Export Citation
  • Nikurashin, M., and Ferrari R. , 2010: Radiation and dissipation of internal waves generated by geostrophic motions impinging on small-scale topography: Application to the Southern Ocean. J. Phys. Oceanogr., 40, 20252042, doi:10.1175/2010JPO4315.1.

    • Search Google Scholar
    • Export Citation
  • Nikurashin, M., and Ferrari R. , 2011: Global energy conversion rate from geostrophic flows into internal lee waves in the deep ocean. Geophys. Res. Lett., 38, L08610, doi:10.1029/2011GL046576.

    • Search Google Scholar
    • Export Citation
  • Oakey, N. S., 1982: Determination of the rate of dissipation of turbulent energy from simultaneous temperature and velocity shear microstructure measurements. J. Phys. Oceanogr., 12, 256271, doi:10.1175/1520-0485(1982)012<0256:DOTROD>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Osborn, T. R., 1980: Estimates of the local rate of vertical diffusion from dissipation measurements. J. Phys. Oceanogr., 10, 8389, doi:10.1175/1520-0485(1980)010<0083:EOTLRO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Polzin, K. L., Toole J. M. , and Schmitt R. W. , 1995: Finescale parameterization of turbulent dissipation. J. Phys. Oceanogr., 25, 306328, doi:10.1175/1520-0485(1995)025<0306:FPOTD>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Polzin, K. L., Kunze E. , Hummon J. , and Firing E. , 2002: The finescale response of lowered ADCP velocity profiles. J. Atmos. Oceanic Technol., 19, 205224, doi:10.1175/1520-0426(2002)019<0205:TFROLA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Sheen, K. L., and Coauthors, 2013: Rates and mechanisms of turbulent dissipation and mixing in the Southern Ocean: Results from the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES). J. Geophys. Res. Oceans, 118, 27742792, doi:10.1002/jgrc.20217.

    • Search Google Scholar
    • Export Citation
  • Shih, L. H., Koseff J. R. , Ivey G. N. , and Ferziger J. H. , 2005: Parameterization of turbulent fluxes and scales using homogeneous sheared stably stratified turbulence simulations. J. Fluid Mech., 525, 193214, doi:10.1017/S0022112004002587.

    • Search Google Scholar
    • Export Citation
  • Smith, W. H. F., and Sandwell D. T. , 1997: Global seafloor topography from satellite altimetry and ship depth soundings. Science, 277, 19561962, doi:10.1126/science.277.5334.1956.

    • Search Google Scholar
    • Export Citation
  • Stöber, U., Walter M. , Mertens C. , and Rhein M. , 2008: Mixing estimates from hydrographic measurements in the Deep Western Boundary Current of the North Atlantic. Deep-Sea Res. I, 55, 721736, doi:10.1016/j.dsr.2008.03.006.

    • Search Google Scholar
    • Export Citation
  • Thorpe, S. A., 1977: Turbulence and mixing in a Scottish loch. Philos. Trans. Roy. Soc. London, A286, 125181, doi:10.1098/rsta.1977.0112.

    • Search Google Scholar
    • Export Citation
  • Thurnherr, A. M., 2012: The finescale response of lowered ADCP velocity measurements processed with different methods. J. Atmos. Oceanic Technol., 29, 597600, doi:10.1175/JTECH-D-11-00158.1.

    • Search Google Scholar
    • Export Citation
  • Tian, J., Yang Q. , and Zhao W. , 2009: Enhanced diapycnal mixing in the South China Sea. J. Phys. Oceanogr., 39, 31913203, doi:10.1175/2009JPO3899.1.

    • Search Google Scholar
    • Export Citation
  • van Haren, H., 2007: Inertial and tidal shear variability above Reykjanes Ridge. Deep-Sea Res. I, 54, 856870, doi:10.1016/j.dsr.2007.03.003.

    • Search Google Scholar
    • Export Citation
  • Visbeck, M., 2002: Deep velocity profiling using lowered acoustic Doppler current profilers: Bottom track and inverse solutions. J. Atmos. Oceanic Technol., 19, 794807, doi:10.1175/1520-0426(2002)019<0794:DVPULA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Walter, M., Mertens C. , and Rhein M. , 2005: Mixing estimates from a large-scale hydrographic survey in the North Atlantic. Geophys. Res. Lett.,32, L13605, doi:10.1029/2005GL022471.

  • Waterman, S., Naveira Garabato A. C. , and Polzin K. L. , 2013: Internal waves and turbulence in the Antarctic Circumpolar Current. J. Phys. Oceanogr., 43, 259282, doi:10.1175/JPO-D-11-0194.1.

    • Search Google Scholar
    • Export Citation
  • Wijesekera, H., Padman L. , Dillon T. , Levine M. , Paulson C. , and Pinkel R. , 1993: The application of internal-wave dissipation models to a region of strong forcing. J. Phys. Oceanogr., 23, 269286, doi:10.1175/1520-0485(1993)023<0269:TAOIWD>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Fig. 1.

    Position of CTD–LADCP stations (yellow and red dots) and VMP stations (red dots) along the A25 OVIDE hydrographic line. The blue arrow indicates a station (48°02′N, 20°51′W) with enhanced dissipation rates at depth (topography from Smith and Sandwell 1997).

  • Fig. 2.

    Example of an average of four PSD of microstructure velocity shear (solid) observed with the VMP for the profile located at 42°35′N, 15°27′W and at a depth of 4200 (gray) and 4270 m (black). Dashed curves are matching Nasmyth (1970) spectra for ε = 7.2 × 10−11 W kg−1 (gray) and for ε = 2.6 × 10−11 W kg−1 (black). Vertical lines give the corresponding Kolmogorov wavenumber.

  • Fig. 3.

    Comparison of ensemble-averaged buoyancy-normalized shear spectra without (raw) and with the use of the transfer function T(k) (see text). The ensemble was made of 589 spectra in the depth range 96–2656 m and from longitude 10°–25°W and latitude 40°–52°N.

  • Fig. 4.

    (a) Buoyancy-normalized shear spectra by depth range (solid) and corresponding GM spectra (dashed). The single-ping error in velocity (dashed–dotted) is given in the legend (cm s−1) for each depth range, taking into account the depth-dependent number of pings (see Table 1). Red triangles located at knoise give the upper limit of spectral integration for each depth range. (b) Strain spectra by depth range (thick lines) and corresponding GM spectra (thin lines). Vertical dashed lines denote the lower and maximum upper limits to determine the strain variances.

  • Fig. 5.

    (a) Dissipation rate profiles εvmp (W kg−1) calculated from microstructure velocity shear data. (b) Corresponding diffusivity profiles Kvmp (m2 s−1). Bathymetry, as derived from CTD casts, is shaded in gray. Note the nonlinear color scales.

  • Fig. 6.

    Turbulence intensity parameter Reb = ε/(νN2). The color scale is nonlinear and emphasizes four mixing regimes (bounded by 1.7, 8.5, 400).

  • Fig. 7.

    Dissipation rate εfinestructure as a function of εvmp. (a)–(d) Color scale denotes log10 (N2) with εfinestructure= εG89, εP95, εG03, εG03, respectively. Solid line denotes perfect agreement with microstructure, dashed lines represent a factor of 3, and dotted lines represent a factor of 10.

  • Fig. 8.

    (a) Shear-to-strain ratio Rω, (b) frequency modulation term h(Rω), and (c) latitude/buoyancy frequency modulation term j(f/N) used in the finescale parameterizations εP95 and εG03.

  • Fig. 9.

    PDF of the ratio εG03/εvmp.

  • Fig. 10.

    Distribution of the ratio εG03/εvmp.

  • Fig. 11.

    (a) Dissipation rate profiles εG03 (W kg−1) calculated from the finestructure of CTD–LADCP data. (b) Corresponding diffusivity profiles KG03 (m2 s−1). Bathymetry, as derived from CTD casts, is shaded in gray. The nonlinear color scales are the same as in Fig. 5.

  • Fig. 12.

    (a) Dissipation rate profiles εζ (W kg−1) calculated from the finestructure of CTD data with a constant shear-to-strain ratio of 5.2. (b) Distribution of the ratio εζ/εvmp.

  • Fig. 13.

    (a) Ratio 〈ω〉/f of the averaged internal wave frequency scaled by the inertial frequency; the white line denotes the threshold 〈ω〉 = 1.1f. (b) Ratio β/α of the low-to-high wavenumber band contribution to the observed shear and strain variances.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 572 265 69
PDF Downloads 286 93 9