## 1. Introduction

In regions of the ocean where temperature and salinity contributions to the density gradient are in opposition, double-diffusive fluxes may equal or exceed those due to mechanical turbulence caused by breaking internal waves (Schmitt 1994; Mack and Schloberlein 1993; St. Laurent and Schmitt 1999). While detailed observations of both phenomena have been made, comprehensive measurements that can identify both processes and distinguish their differing roles in the evolution of the thermocline are rare. Such insights are required in order to advance mixing parameterizations in large scale numerical models.

*S*), 2-m strain (

*γ*=

*N*

^{2}

*N*

^{2}, where

*N*is the buoyancy frequency), 2-m effective strain rate

*N*

^{2}/

*S*

^{2}). They observe an increased likelihood of overturning in the presence of small (5–20-m vertical scale) internal waves, indexed by high (∂

*ŵ*/∂

*z*)

^{2}. High strain (weak stratification) and low 6.4-m Ri also tend to accompany overturns.

In the present work, the focus is broadened to consider double-diffusive phenomena. Intensive depth–time measurements of temperature, salinity, density, and 8-cm microconductivity fluctuations were sampled every 4 minutes during MBL over a 300-m vertical aperture. Estimates of traditional microstructure quantities are formed from these directly measured microconductivity data, and related to the overturns and the background finescale *T–S* fields. Mixing structures are repeatedly sampled as they drift past the Research Platform *FLIP.* The ability to track individual “microstructure” patches in depth/density and time greatly facilitates the effort, as does the ability to examine patch evolution in an isopycnal-following frame.

Centimeter-scale conductivity gradient fluctuations and meter-scale density overturns represent distinctly different phenomena. To oversimplify, overturns initiate turbulent motion, while conductivity microstructure can represent the “aftermath” (AP), or indicate double-diffusive activity. Overturns are restored to gravitational stability within several buoyancy periods (15–30 min in the present observations). The resultant conductivity signals (which are dominated by temperature) can persist much longer. Once formed, passive temperature fluctuations at the 8-cm limit of our probe’s resolution would be diffused away in *τ* = *L*^{2}*K*^{−1}_{T}

While microstructure can result from either mechanical turbulence or double-diffusive processes, large (>2 m) density overturns are not likely to be formed by double diffusion, whose net effect is always to *increase* density stratification.

We approach the present observations with the hypothesis that both mechanical and double-diffusive fluxes are present. The corresponding conductivity signals are termed “turbulent” and “double-diffusive” microstructure.

Turbulent microstructure and overturning should display generally similar relationships with the Ri and (∂*ŵ*/∂*z*)^{2} fields, which index internal wave instabilities. Double-diffusive microstructure, on the other hand, should occur in regions where the density ratio *R*_{ρ} and (equivalently) Turner angle Tu (Ruddick 1983) indicate that the water is unstable to either salt fingering or diffusive convection.

As in AP, we attempt to use dense depth–time sampling to gain a view of the evolution of upper-ocean processes. We anticipate the presence of mixing “patches,” “puffs,” and “wisps,” (Gregg 1987) as well as zones of salt fingering and diffusive regime that accompany thermohaline intrusions. The various structures presumably result from modulations in the various finescale fields, which index shear, convective, salt fingering, and other instabilities. Given the multiplicity of phenomena at MBL, the discussion takes the tone of a survey, rather than an exhaustive account.

Descriptions of the experiment (section 2), and the instruments and methods used (section 3), follow. We then present depth–time maps of the entire 12-day cruise (section 4), and 12-h close-ups (section 5). In section 6, we undertake a statistical summary of the variety of phenomena associated with microstructure. Section 7 compares dissipation rate estimates based on microconductivity observations to those inferred from overturning scales. We summarize and conclude in section 8.

## 2. Experiment

Phase I of the Marine Boundary Layer Experiment was conducted from the Research Platform *FLIP,* 30 km west of Point Argüello, California (35°N, 121°W), during February–March 1995. Data used in this study were collected over the 12-day period yearday 56–68. During the first 10 days, FLIP was secured with a one-point mooring in 1.5 km of water. On 6 March (yearday 66), the mooring line was released, and *FLIP* drifted freely. An eight-beam Doppler sonar system and a repeat-profiling CTD sampled the upper 400 m of the water column. An SBE dual-electrode microconductivity cell was mounted in a special housing below the CTD package. Details of the experiment, and a schematic of the CTD/sonar array, may be found in AP.

High shears and large baroclinic tidal displacements characterized the MBL site throughout the experiment (AP). Rms 6.4-m shear is about 1.4 times the Garrett–Munk 1975 (GM) value, for depths above 250 m. Below this, mean shear increases with depth (rms shear is 3 times the GM value at 400 m), suggesting a deep source of wave energy at the site. Rms baroclinic tidal displacements are typically 20–30 m. Larger (40–50 m) displacements occurring in a limited depth range near 150 m during yearday 60–62 are associated with the passage of baroclinic energy along a tidal ray, presumably stemming from nearby topography.

This region of the California coast is known to contain *T–S* variations (fronts, “squirts,” “jets,” and intrusions) on a variety of time/space scales. These modulate the double-diffusive stability of the water column (e.g., Gregg 1975). Of particular relevance to the present study is the passage of three sloping *T–S* structures ∼50 km wide over yeardays 59–62 (section 4). Layers of double-diffusively unstable water bound these features above and below. In addition, smaller thermohaline intrusions continually advect through the MBL site. Bursts of heightened temperature dissipation rates *χ̂*

## 3. Instruments and methods

### a. Sensors

During MBL, velocity profiles (0–400 m) are obtained from a Doppler sonar mounted on the stern of *FLIP* at a depth of 87 m. Every minute, estimates of velocity (3.2 m vertical resolution) are recorded. Averages are subsequently computed over 4 minutes to facilitate merging with the CTD data. Use of repeat sequence codes (Pinkel and Smith 1992) improves precision of the system.

An automatic CTD profiles temperature, conductivity, and pressure every 4 min from the surface to 415 m. The CTD, a modified Sea-Bird Instruments SBE-9 unit (Fig. 1a), falls at approximately 3.6 m s^{−1}. At this speed, the 24-Hz sample rate corresponds to a measurement every 16 cm. Conductivity data are corrected for thermal inertia effects (Lueck and Picklo 1990; Morison et al. 1994). A cross-spectral technique, described in Anderson (1993), is then employed to match the responses of the temperature and conductivity probes. Data are low-pass filtered with a 2-m cutoff prior to calculating salinity and density profiles.

Microconductivity data are obtained from a Sea-Bird Instruments SBE 7-02 microconductivity sensor (Fig. 1b). The probe consists of two electrodes separated by approximately 1 mm. It is mounted 80 cm below the main CTD package in a special protected unit. At the 3.6 m s^{−1} fall speed, a 96-Hz sample rate yields data every 4 cm. Prior to digitizing, the analog signals are passed through two antialiasing filters (4 pole Bessel with a time delay of 19.1 msec). These provide a minimum of 20-dB signal reduction for frequencies above the Nyquist.

Ten sample microconductivity profiles from 100 to 400 m near yearday 64.1 (Fig. 2a) demonstrate the probe’s ability to detect large- and small-scale signals. Fifty-meter close-ups of microconductivity (Fig. 2b) and microconductivity gradient (Fig. 2c) data surrounding the energetic “patch” (Fig. 2a) show irregular features, steps, and inversions over a 30-m depth range. Numerous 2-m overturns in density (not shown) accompany the patch. Since the region is stable to double diffusion, the patch presumably results from mechanical turbulence.

To investigate the response of the microconductivity probe, the conductivity gradient spectrum is estimated from 10 profiles inside the “patch” (Fig. 2d, heavy black line). A 512-point (20 m) section of each profile, centered on the instantaneous location of the patch, is tapered with a triangle window and Fourier transformed. The 10 resulting spectral estimates are averaged together, and smoothed over Δ(log_{10}*k*) = 0.2 before plotting. The spectrum is expressed as an equivalent temperature gradient spectrum *φ̂*_{TG}(*k*)*μ* by *b* ≡ ∂*T*/∂*μ* = 9.16 [°C/(S m^{−1})].

Significant attenuation is observed for *k* > 6 cpm, due to the finite sampling rate of the conductivity probe [based on measurements in a stratified tank, Meagher et al. (1982) estimate that the inherent spatial response of the SBE 7-02 probe is uniform to 100 cpm]. It is possible to correct the spectrum for degradation associated with finite sampling. Formally, the inverse of the antialiasing filter applied in the instrument’s A/D processing, expressed in terms of wavenumber, is the appropriate correction. However, *T*(*k*) = sinc^{2}(*πk*/*k*_{N}) [where sinc(*x*) = sin(*x*)/*x,* and *k*_{N} = 13.3 cpm is the Nyquist wavenumber] is a good approximation to the Bessel filter (Fig. 3). The sinc^{2} filter is the appropriate spectral correction if each sample represents a uniformly weighted average over an interval 1/(2*k*_{N}). When the sinc^{2} correction is applied to the wavenumber spectrum, a spectral slope consistent with +⅓ (Fig. 2d, dotted line) is evident for 2 cpm > *k* > 9 cpm. This is consistent with the inertial subrange of the scalar dissipation spectrum, suggesting that the probe is in fact detecting turbulent motions.

*σ*

^{2}

_{N}

^{−10}(S m

^{−1})

^{2}is estimated. Assuming that noise is white in microconductivity, noise in the gradient should exhibit a

*k*

^{+2}slope:

Revisiting the active turbulent patch (Fig. 2d), the modeled noise spectrum is found to be well below the observed spectral levels for *k* < 9 cpm, increasing until it approaches the measured spectral level near 9 cpm (*S*_{N} = 0.62 × 10^{−4} °C^{2} cpm^{−1} at 9 cpm). For *k* > 9 cpm, the corrected spectrum begins to rise above *k*^{+1/3}, indicating that noise may be significant at the highest wavenumbers. Still, in the most active regions such as this one (which dominate contributions to *χ̂**Ĉ*), noise appears to be a minor effect.

### b. *χ̂* and *Ĉ*

*χ̂*

Temperature-fluctuation dissipation rate, *χ̂,**Ĉ,* can be estimated from the microconductivity profiles. We emphasize that the 4-cm sampling of our data is inadequate to resolve most of the variance in *χ**C.* These quantities are underestimated by a dissipation-rate-dependent factor of 6–100 (appendix). Use of the caret (*χ̂,**Ĉ*) distinguishes estimates from the fully resolved quantities. We optimistically assume that space/time patterns in *χ̂**Ĉ* are not significantly different those of the true *χ**C.*

*μ*(

*z*), the Cox number estimate is given by

Overbars represent 1-m vertical averages centered about *z*_{i}. The factor of 3 accounts for the contribution from horizontal gradients, under the assumption (suspicious at these scales) that the conductivity fluctuations are isotropic.

In regions of small background mean gradient the denominator in (2) may be unduly influenced by measurement noise. Very large, spurious, *Ĉ* estimates can result. To minimize this influence on Cox number variability, a threshold denominator *D*_{min} = 10^{−10} (S m^{−1})^{2} is enforced. Occurrences of denominator *D* < *D*_{min} are replaced by *D*_{min} before Cox number is computed. The appearance of the depth–time maps is not influenced by the exact choice of *D*_{min}.

*χ̂,*

*χ̂*

*Ĉ*it is assumed that conductivity variability is dominated by temperature at the scales of interest; that is, ∂

*T*/∂

*z*=

*b*∂

*μ*/∂

*z.*

### c. Semi-Lagrangian reference frame

Interpretation of Eulerian records is complicated by the vertical advection of the thermocline by internal waves. We choose to examine the time evolution of *χ̂**Ĉ* in an isopycnal following, or semi-Lagrangian (s-L), reference frame. To compute s-L quantities, a set of reference isopycnals, separated on average by one meter, is formed using the SBE-9 density profiles. Linear interpolation is used to map the Eulerian observations onto isopycnal surfaces. “Depth–time” s-L maps are presented as functions of the mean depth of each reference isopycnal, and time. On these plots, the cross-isopycnal migration rate is given by the apparent slope of isosurfaces.

### d. Ri and (∂*ŵ*/∂*z*)^{2}

*N*

^{2}/

*S*

^{2}is a measure of the potential for Kelvin–Helmholtz instability. Parallel, stratified shear flow is stable to small perturbations for Ri > ¼ (Miles 1961; Howard 1961). Following AP, 6.4-m Richardson number estimates are formed using combined CTD and sonar data. Shear-squared profiles (

*S*

^{2}) are computed from 6.4-m vertically differenced velocity. Profiles of buoyancy frequency

*N*are formed by differencing the density profiles after “re-sorting” the overturns to enforce monotonicity. These are then smoothed with a 6.4-m Butterworth filter. The 6.4-m Ri is computed from these instantaneous

*S*

^{2}and

*N*

^{2}profiles. Here ∂

*ŵ*/∂

*z*is computed from CTD data as in AP. Briefly, between any two time series of isopycnal depth

*z*(

*t, ρ*

_{i}) and

*z*(

*t, ρ*

_{j}), the normalized depth difference is given by

*z*

*z*

*ρ*

_{i}) +

*z*

*ρ*

_{j})]/2 is the mean depth of the pair, and

*z*

*z*

*ρ*

_{i}) −

*z*

*ρ*

_{j}) is the mean separation. The effective strain rate is then defined by

*ŵ*/∂

*z,*and the true strain rate or vertical divergence, ∂

*w*/∂

*z,*since advection of horizontal structure can contribute significantly to ∂

*ŵ*/∂

*z.*

Depth–time maps of ∂*ŵ*/∂*z* reveal wavelike features (AP, Fig. 13). Individual crests within well-defined groups can often be seen. The effective strain rate magnitude, (∂*ŵ*/∂*z*)^{2}, highlights the “envelopes” of these groups and is strongly correlated with the location of overturning events (AP).

### e. Turner angle

To estimate Tu, temperature and salinity are computed on isopycnals. Differences are computed between isopycnals whose mean separation is 4 m. Here Tu is then calculated using Eq. (6).

## 4. Twelve-day maps of scalar dissipation and Cox number

### a. General features

Semi-Lagrangian depth–time maps of _{10}(*χ̂*_{10}(*Ĉ*) (Fig. 4b) are presented for the 12-day MBL record. Since signals above 87 m (the bottom of *FLIP*) are potentially contaminated by *FLIP*’s wake, we restrict our study to a 100–400-m range. Prior to plotting, the logarithm of each quantity is averaged over 1 h in time and 15 m in depth. Black vertical bars represent periods where the CTD and/or the microconductivity probe was out of service.

Here *χ̂**χ̂**Ĉ* with double-diffusive processes. Shorter-lived patches (5–20 m thick, duration several hours), such as the one at yearday 64.1, depth 200 m (shown in Fig. 2) appear to be of mechanical origin.

The scalar dissipation *χ̂**Ĉ* by the mean conductivity gradient squared [Eq. (2)]. Thus, *χ̂*

Heightened *χ̂**Ĉ* levels are present at all depths beginning at yearday 59. Subsequently, three features (beginning at about 100, 220, and 280 m at yearday 59) appear to migrate downward. The middle, most intense one may be followed for a day and a half, when it disappears near 300 m. The top one disappears and reappears several times, but may be tracked until yearday 64. Each extends approximately 20 m in the vertical and appears to migrate downward at approximately 50 m day^{−1}. These features are associated with the advection of sloping, double-diffusively unstable lenses of water past the experimental site. The interpretation of observed temporal evolution as advected horizontal structure is discussed in section 4b. The double-diffusive stability of the lenses is discussed in section 4c.

### b. Advection

Advection affects the observed temporal variability significantly in these moored observations. A progressive vector diagram for the 12 days (Fig. 5) provides a useful index of the flow of water past *FLIP.*

While moored (yeardays 59–66), the current trajectory reflects the superposition of strong westward flow and weaker inertial motions (evident as daily “wiggles” in the progressive vector). During this moored period, no loops are present, indicating that inertial motions are not strong enough to cause the same water to pass *FLIP* twice. Nevertheless, the time coordinate of our depth–time maps indicates a mix of intrinsic and advective variability. We can model the persistent structures as “frozen” in time,^{1} and we can use observed advective velocities to estimate lateral scales. Using a nominal value of 20 cm s^{−1}, one day corresponds to 17 km.

### c. *T–S* structure

Coincident with the large, slanted *χ̂,**Ĉ* features beginning at yearday 59, three sloping frontal features pass underneath the experiment site, evidenced by distinct changes in *T* and *S* along isopycnals. The observed heightened microconductivity levels appear due to double-diffusive processes occurring at the interfaces.

We investigate this hypothesis by examining records of isopycnal temperature anomaly (Fig. 6a). Density is primarily a function of temperature and salinity alone, given the limited depth range of these isopycnal surfaces. An analogous plot of isopycnal salinity anomaly therefore shows identical structure.

Temperature variation along an isopycnal is not sensitive to internal wave straining. Variability in Fig. 6a is therefore due to the horizontal advection of water of varying *T–S* characteristics past the experimental site. Prior to yearday 59, warmer, saltier water (relative to the 12-day mean) overlies colder, fresher. The water undergoes a sharp transition to colder, fresher (yearday 59), and back again (yearday 59.5). Above 230 m, all isopycnals experience the change simultaneously; below this, the transitions slope downward.

The layered structures are reminiscent of interleaving layers typically observed near fronts. An idealized interpretation of the layered structures in terms of frontal interleaving is possible. A vertical boundary between cold, fresh water (to the left in Fig. 6) and warm, salty water (to the right) is encountered at yearday 59.5. In accordance with the laboratory results of Ruddick and Turner (1979), a “nose” of cold, fresh water penetrates into the warmer water to the right. This interleaving process is assumed to have occurred at a previous time. The resultant interface–nose structure then sweeps past *FLIP.*^{2}

Assuming all observed time changes are due to advection, the cold layer slopes downward to the southwest at about (50 m day^{−1})/(17 km day^{−1}) = 0.3% or 0.26°. The protruding layer’s downward slope away from the interface is consistent with higher buoyancy fluxes due to salt fingering at the top of the layer than by the diffusive regime at the bottom. With a net buoyancy flux out of the layer, it sinks into the surrounding fluid as it advances (Stern 1967). This behavior is observed in laboratory experiments with salt and sugar (Turner 1978), and in the ocean (Gregg 1980; Joyce et al. 1978; Gregg and McKenzie 1979; Ruddick 1992; Anderson and Pinkel 1994).

*T–S*gradient has been estimated by Ruddick and Turner (1979). They consider the potential energy in a layer, before and after the

*T–S*gradients have been “run down” by double-diffusive fluxes. This vertical scale is given by their Eq. (11)

*n*=

*αF*

_{T}/

*βF*

_{s}is the flux ratio, expressed in terms of density (

*n*= 0.56 for

*T–S*fingers). Along an isopycnal,

*β*|Δ

*S*| =

*α*|Δ

*T*|. Substituting into (7) and expressing the density gradient in terms of

*N*

^{2},

*T*= 0.4°C,

*N*= 2.53 cph, and

*α*= 1.5 × 10

^{−4}°C

^{−1},

*g*= 9.8 m s

^{−2}yields

*H*≈ 20 m, consistent with the observed vertical scale of the lenses. Based upon their downward slope away from the vertical interface, and their vertical scale, the layers appear to be double-diffusively driven.

The protruding cold, fresh layers in Fig. 6b (e.g., yearday 60, 250 m) are associated with the potential for salt fingering (blue, Tu > 45°) on top, and diffusive regime (red, Tu < −45) below. The opposite holds for hot, salty layers. Regions where *χ̂*^{−9} (°C)^{2} s^{−1} are contoured with black lines in both 6a and 6b. It is clear from Figs. 4, 6 that the protruding lenses are associated with large *χ̂*

Almost all of the persistent regions of elevated *χ̂**χ̂*

Of course, |Tu| > 45° only implies potential for double diffusion. The occurrence of mechanical turbulence is presumably unaffected by double-diffusive phenomena. Patches of turbulence often occur in double-diffusively unstable regions, as appears to be the case in some of the shorter-lived patches apparent in the 12-h maps (next section). But in the 12-day records, the patterns in Tu and *χ̂*

## 5. Twelve-hour maps of scalar dissipation and Cox number

The relationship between the microstructure and the fine-scale fields can be examined in more detail in 12-h maps. Three pairs of (*χ̂,**Ĉ*) maps are presented for yearday 61.5–62.0, centered on one of the slanted *T–S* lenses. The first (Fig. 7) is in an Eulerian frame. The last two (Figs. 8, 9) are in a s-L frame. In these, regions are contoured where extreme values of the finescale fields Ri, (∂*ŵ*/∂*z*)^{2}, and Tu occur. Viewing the three presentations together, one can gain an appreciation of the variety of instability processes that appear to be co-occurring.

### a. Eulerian representation of *χ̂* and *Ĉ* fields

*χ̂*

Twelve-hour Eulerian records of _{10}(*χ̂*)_{10}(*Ĉ*) are presented in Fig. 7. Prior to plotting, the logarithm of each quantity has been averaged over 4 m in depth and 12 min in time. Isopycnal surfaces, separated in the (12-day) mean by 20 m, are plotted as thin black lines. Locations of density overturns >2 m (AP) are marked with black dots.

The most prominent feature in both fields extends from the start of the record to the finish and is associated with one of the *T–S* lenses (Figs. 4 and 5). Midway through the record, a second, shallower, lens appears and eventually converges with the first.

Elsewhere, both *χ̂**Ĉ* display the “patchy” nature of microstructure. Thin (<10 m) layers of elevated scalar dissipation, *χ̂,**χ̂**Ĉ* (e.g., yearday 61.77, 150 m). Features in both fields typically last several hours, implying length scales of several kilometers assuming advection at 20 cm s^{−1}.

In general, layers of *χ̂**Ĉ* demonstrate a strong tendency to follow isopycnal surfaces. Even rapid displacements are often tracked by both microstructure and overturning, as in the feature beginning at yearday 61.75, depth 300 m. Over the course of two hours, the high *Ĉ* feature remains between two isopycnals, which complete two up–down cycles of about 20 m amplitude. High *χ̂*

Observed 2-m overturns display a tendency to lie within high *χ̂**Ĉ* regions. Since salt fingering is not expected to generate large density overturns, mechanical turbulence is likely present in these regions. Note even the double-diffusively unstable lens discussed in the last section displays overturning, implying the simultaneous presence of double-diffusive and mechanical turbulent microstructure.

Better correspondence with overturning is seen in the *Ĉ* field than in *χ̂,**N* scaling of the two quantities. Here *χ̂**Ĉ* has this tendency removed. It has been found that large overturns are more likely in weakly stratified regions (AP). This emphasizes the importance of comparing quantities that scale similarly with *N.*

A few overturns are seen in low Cox number regions (e.g., 150 m and yearday 61.95, 150 m and yearday 61.72, 170 m and yearday 61.8). These may be sporadic, isolated events in which the overturn occurs below *FLIP* but the aftermath is “swept away.”

### b. Mechanical turbulence: Comparison with Ri and (*∂ŵ*/∂*z*)^{2}

Semi-Lagrangian quantities *χ̂*_{s-L}*Ĉ*_{s-L} for the same period and depth range are plotted in Figs. 8a,b. Observed 2-m overturns are overplotted as in Fig. 7. Regions where 6.4-m Ri < 1 are contoured with heavy gray lines. Heavy black lines correspond to regions where (*∂ŵ*/∂*z*)^{2} exceeds twice its mean value over the 12-h, 300-m domain.

Horizontal traces (where *χ̂**Ĉ* layers remain on specific isopycnals) dominate the record. However, some cross-isopycnal migration is also present (e.g., yearday 61.6, 200 m). Without replenishment, turbulent patches are restricted to remain on isopycnals. In some circumstances, microstructure appears to follow disturbances as they propagate past isopycnals. Such features are reminiscent of a spilling breaker, wherein continually generated whitewater follows the passage of the surface bore.

Strong correspondence is observed between both (*Ĉ,* *χ̂**ŵ*/∂*z*)^{2} (heavy black lines). Apparently, microstructure patches, like overturns (AP), are more likely where (∂*ŵ*/∂*z*)^{2} is high.

Here *χ̂**χ̂**χ̂*

Better agreement is seen between *Ĉ* and 6.4-m Ri. Low Richardson number regions often display high Cox number (e.g., yearday 61.55, 350 m; yearday 61.9, 175 m). However, examples exist of low 6.4-m Ri regions that persist for hours without a *Ĉ* signature (yearday 61.8, 180 m). Analogous persistent zones of low 6.4-m Ri that demonstrate no 2-m overturns are documented in AP and are also seen in these records. Evidently, instability does not always develop when 6.4-m Ri is low.

Likewise, *Ĉ* may also be large where the 6.4-m Ri > 1. The region near 200 m, yearday 61.55 is an example, as is the feature that migrates from 240 to 280 m near yearday 61.8. High (∂*ŵ*/∂*z*)^{2} contours and 2-m overturns track the *Ĉ* signatures, indicating that instability of small-scale internal waves may have led to the microstructure. As discussed in AP, shear at unresolved horizontal and vertical scales causes the “true” Ri to be lower than that reported here. Alternatively, such features may be persistent “patches” [shown by Gregg et al. (1986) to be poorly correlated with low 10-m Ri] or double-diffusive structures (for which Ri is irrelevant). The Turner angle for this region (next section) is examined to isolate the latter possibility.

### c. Double diffusion: Comparison with Turner angle

The largest *Ĉ* and *χ̂**χ̂**Ĉ* fields are replotted in Fig. 9. Black lines now enclose regions unstable to salt fingering (Tu > 45°). Gray lines correspond to the diffusive regime (Tu < −45°). This 12-h period is coincident with the passage of the *T–S* lenses; therefore, much double-diffusive activity is anticipated. Indeed, double-diffusively unstable layers 10–25 m thick pervade the record below 220 m. The sloping features evident in the 12-day maps (Fig. 6) are visible here, though the apparent aspect ratio is altered. High *χ̂**Ĉ* levels are associated with most layers. However, some layers exhibit no microstructure (300 m, yearday 61.5 and 320 m, yearday 61.9).

Here *χ̂**Ĉ.* Quite strong correspondence is seen with the lenses, as expected from section 4. For example, the Tu > 45° layer near 230 m, yearday 61.6 bifurcates. Associated microstructure tracks the split. Several regions of Tu > 45° exhibit high *χ̂**Ĉ* (indicating that the feature occurred in a high gradient region). Examples are at (230 m, yearday 61.5–61.7), and (330 m, yearday 61.8–62). The remaining high *Ĉ* regions appear better correlated with Ri and (∂*ŵ*/∂*z*)^{2} (Fig. 8b). The suggestion is that *Ĉ* is a better indicator than *χ̂**χ̂,**Ĉ*), while overturns, low Ri, and high (∂*ŵ*/∂*z*)^{2} all occur in weak stratification (AP).

Salt fingering and diffusive structures associated with the passage of the lenses appear to make up a large portion of the microstructure in this 12-h record, suggesting that double-diffusive fluxes may be important in depth-averaged budgets. Indeed, estimates of kinetic energy dissipation rate, ε, based on the assumption that all observed *χ̂*

### d. Double-diffusive processes and overturning

Salt fingering generates stabilizing density fluxes, and therefore cannot produce large overturns in density. On the other hand, mechanical turbulence is generated by density and shear fluctuations, and should therefore occur approximately independently of the *T–S* characteristics of the water column.

As a simple test of the independence of overturning and double-diffusive processes, the probability density of Turner angle (Fig. 10, black line) is computed for all observations between 200 and 350 m, yeardays 59–68 (the same range used in other PDFs computed in section 6, and in AP). The probability density is then computed using only data exhibiting 2-m overturns (dotted line).^{3} The technique allows statistical evaluation of the degree to which overturning is influenced by the value of local Tu.

The two distributions differ, but not greatly. The largest differences in the two PDFs appear in the double-diffusively stable regime, |Tu| < 45°. The two populations demonstrate nearly identical percentages of the water occupied by the various diffusive regimes (Table 1), suggesting that large-scale overturning occurs independently of double-diffusive processes.

### e. Profiles and spectra inside double-diffusive regions

The small-scale structure within the intrusive, double-diffusive regions is of interest. Close-up profiles of two cold, flesh intrusions are presented in Figs. 11, 12. In the first (yearday 61.75, 310 m), elevated microstructure occurs in the salt-fingering-favorable region above the intrusion. In the second (yearday 61.8, 265 m), the more intense activity is in the diffusive-favorable region below. These displays are direct analogs of Fig. 2 (mechanical turbulence), spanning 40 min in time.

Whereas the conductivity profiles in Fig. 2a decrease monotonically on large scales (implying double-diffusive stability), the passage of a cold, fresh “nose” is evident in Fig. 11a near 310 m. Salt fingering is expected when warm, salty water overlies cold, fresh at the top of the layer. Irregular features are seen here, on the strong gradient above the cool nose (11b). Large fluctuations in microconductivity gradient (11c) are present above a relatively quiescent region in the center of the intrusion. Weaker fluctuations below the layer (where cold, fresh overlies hot, salty) appear to be associated with the diffusive regime.

The characteristic staircase signature of double diffusion is absent. Indeed, more “steps” appear in the “turbulent” patch (Fig. 2) than in Fig. 11. However, the spectrum (computed as in Fig. 2d) exhibits slope steeper than *k*^{+1/3} (11d), as is expected for spectra of salt fingering fluctuations (Gargett and Schmitt 1982). We caution that these spectra are pushing the probe’s vertical resolution limit.

Figure 12 also shows a cold, fresh intrusion. Again, activity is detectable above and below the intrusive lens. However, in this case the stronger, thicker signal appears below, in the diffusive regime. Here *φ*_{TG} for this region exhibits a slope closer to *k*^{+1/3}. Perhaps convective motions within the diffusive cell have become turbulent.

In spite of the differences in spectral slope, no obvious differences are discernible by eye between the fluctuations within mechanically turbulent, salt fingering, and diffusive regime patches. As noted by Linden (1971), weak background turbulence may disrupt the double-diffusive structures.

## 6. *C* and *χ*: A statistical summary

Following the approach taken by AP, general trends in the depth–time maps are summarized using probability density functions (PDFs). These are formed using the same data as in AP; that is, the range 200–350 m, over a 9-day period from yearday 59–68. PDFs of *χ̂**Ĉ* are formed in the Eulerian frame, using all 334 012 observations in this depth–time range. Then, dependence upon other variables is examined by forming conditionally sampled PDFs, using subsets of data that satisfy certain conditions in these other variables. These include situations where overturning is observed, unstable Turner angle, low 6.4-m Ri, and high (∂*ŵ*/∂*z*)^{2}. Table 2 summarizes the various PDFs computed.

Heavy black lines (Fig. 13) represent full-population probability densities of *χ̂**Ĉ.* The PDF of _{10}*χ̂*^{−10} °C^{2} s^{−1}. Lognormal *χ* distributions are often observed (e.g., Gregg et al. 1986). Recall that *χ̂**χ.**χ̂*

The PDF of log_{10}*Ĉ* displays a central tendency, but is skewed. The true Cox number often has a lognormal distribution (J. Sherman 1998, personal communication), though it is not strictly expected to (Davis 1996). Were *Ĉ* fully resolved, the PDF would be shifted to the right. Since the right side of the distribution would be shifted more than the left, accurate *C* estimates might appear more nearly lognormal than the underresolved quantity, *Ĉ.*

Subsampled PDFs (light lines) are shifted relative to the general population PDFs. The overturning population (black line) demonstrates a very slightly shifted *χ̂**χ̂**χ̂**χ̂**χ̂*

By contrast, the PDF of Cox number, considering overturns only, displays a marked right shift relative to the overall population PDF. This is consistent with visual conclusions drawn from the maps.

As expected from the maps, PDFs of *χ̂* and *Ĉ**χ̂**Ĉ,* again because *Ĉ* is less sensitive than *χ̂*

For both *χ̂**Ĉ,* the PDF of the diffusive subsample is right shifted slightly relative to that of salt fingering. Based on this result and close-ups of individual regions (Figs. 11, 12), both the fingering and diffusive forms of double diffusion appear able to generate strong conductivity fluctuations.

Considering only regions of 6.4-m Ri, the distribution is slightly left shifted relative to the general population. The poor correspondence of Ri with persistent mechanical turbulence patches (Gregg et al. 1986), and with double-diffusive structures, is responsible. In addition, the different dependence upon stratification of *χ̂,**Ĉ,* where this problem is avoided, differences are seen in the expected sense: the low Ri population exhibits a higher *Ĉ* distribution.

The strongest shift in both curves is seen in the subsample exhibiting (∂*ŵ*/∂*z*)^{2} > 2 × 10^{−5} s^{−2}. Evidently, (∂*ŵ*/∂*z*)^{2} (which indexes the presence of small-scale wave groups) is a robust indicator of microstructure, as it is of overturns (AP).

## 7. Depth-averaged dissipation rate

Microconductivity fluctuations in MBL are due to a combination of turbulent and double-diffusive processes. Time series of a depth-averaged “effective dissipation rate” (ε_{μ}) can be computed by assuming that all fluctuations are due to mechanical turbulence. We anticipate that ε_{μ} will be a gross overestimate of true dissipation if double-diffusive processes, rather than turbulence, produce the fluctuations. By contrast, the dissipation rate inferred from observed density overturns (ε_{T}) is not affected by double-diffusive processes. By comparing these two different estimates of dissipation rate, the relative importance of the two sources of microstructure can be assessed.

The ε_{μ} estimate assumes that all fluctuations are due to turbulence. Assuming that the conductivity spectrum takes the Batchelor form and that the dissipation rate is simply related to potential energy production through a constant mixing efficiency Γ = 0.2 (Oakey 1982), ε_{μ} may be estimated in terms of observed *χ̂.*

We emphasize that ε_{μ} is an approximate measure based on numerous assumptions, some of which do not strictly hold in the present case. For example, the cruise-averaged *R*_{ρ}, rather than its in situ value, is used in the calculation. Given the variability of *R*_{ρ} (Fig. 6), this assumption is clearly suspect. The hope in computing ε_{μ} is that these and other errors do not bias the depth and time averages significantly.

_{T}is computed from the overturning (Thorpe) scale

*L*

_{T}. The method is based on Dillon’s (1982) observation that the Ozmidov scale

*L*

_{O}=

*N*

^{3}

*L*

_{T}. The technique is described in AP. Briefly, for each 4-min CTD profile, the mean dissipation rate is evaluated by computing the volume average of the dissipation produced by the overturns:

_{i}= 0.64

*L*

^{2}

_{T}

*N*

^{3}is the Thorpe-inferred dissipation rate accomplished by the

*i*th overturn. Here

*L*

_{P,i}is the vertical extent (patch size) of the

*i*th overturn, and

*H*is the depth range sampled. The method produces estimates that demonstrate the same magnitude and dependence upon shear and stratification as the Gregg (1989) parameterization (AP), which seems to be quite robust (e.g., Polzin et al. 1995). Since no “ground truth” dissipation rate measurements are available for MBL, we settle upon ε

_{T}as a reasonable proxy.

Time series of ε_{μ} and ε_{T} (Fig. 14a) show extremely strong correspondence in spite of the numerous assumptions made. Here ε_{μ} nearly always exceeds ε_{T}, but during most of the record the difference is less than 40% of ε_{T}. No ad hoc scaling of either series has been made. The implication is that during regions of close agreement, observed microstructure is mainly turbulent and not double diffusive. This agreement increases our confidence that the microconductivity probe is indeed responding to turbulent signals.

Poor agreement (ε_{μ} ≫ ε_{T}) is expected if observed *χ̂*_{T}/ε_{μ} is plotted in 14b as a thick gray line. Also plotted are percentages of the water column that are subject to salt fingering *p*_{SF}, dotted line), diffusive regime (*p*_{DR}, dash–dot line), and either double-diffusive process (*p*_{DD} = *p*_{SF} + *p*_{DR}, solid black line). The stable percentage *p*_{stable} = 1 − *p*_{DD} is plotted with a heavy black line. Comparing these series to Figs. 4 and 6, the sloping lenses (yearday 59–62) occupy large fractions (∼40%) of the sampled water. The greatest differences between ε_{μ} and ε_{T} (yearday 59–60.5, yearday 67–68) occur when much of the water column is double-diffusively unstable (Fig. 14b). The best agreement (yearday 64–65) occurs when less than 10% of the water column is double-diffusively unstable. The agreement in this region is interpreted as evidence that both the Thorpe method and the microconductivity method are reasonable estimates of depth-averaged dissipation rate, provided that double-diffusive processes do not dominate. Large overestimates in the lens regions suggest that buoyancy fluxes supported by double-diffusive processes can be significant.

## 8. Summary and conclusions

Eight-centimeter scale microconductivity fluctuations have been monitored over a 12-day, 300-m domain in the California Current. The site demonstrates frequent overturns, large tidal displacements, and strong shears. These increase with depth below 250 m, implicating a deep or bottom energy source.

Between yearday 59 and 62, three sloping thermohaline intrusions drift past. Their vertical scale and sense of slope are consistent with laboratory models of double diffusion (Ruddick and Turner 1979). Elevated *χ̂**Ĉ* are observed at the boundaries of these features, in regions where |Tu| > 45. The lenses and associated microstructure appear to be double-diffusively driven.

Close-up maps of *χ̂**Ĉ* over a 12-h period centered on these features reveal that while the most persistent *χ̂,**Ĉ* patches are associated with double-diffusively unstable Tu values, short-term correspondence is also found between *Ĉ* and both 6.4-m Ri and (∂*ŵ*/∂*z*)^{2}: *χ̂**ŵ*/∂*z*)^{2} as well, but not with 6.4-m Ri (because high stratification, persistent mixing patches, and double-diffusive processes, which are poorly or anticorrelated with Ri, all lead to high *χ̂**χ̂**Ĉ.* Since overturning, Ri, and (∂*ŵ*/∂*z*)^{2} are associated with mechanical turbulence, it appears that both mechanical turbulence and double-diffusive processes are active at the site.

This interpretation is supported by comparison of time series of depth-averaged dissipation rate inferred from overturning scales (ε_{T}; Dillon 1982) and from the microconductivity fluctuations (ε_{μ}; appendix) using a model that assumes a Batchelor scalar spectrum and a constant mixing efficiency. Good correspondence is observed between the two dissipation rate estimates during periods when only a small fraction of the water column is double-diffusively unstable. In the lens regions, where much of the water is double-diffusively unstable, ε_{μ} ≫ ε_{T}, consistent with microconductivity variance stemming from double-diffusive processes.

## Acknowledgments

This work was funded by the Office of Naval Research under the Marine Boundary Layer Program. The authors wish to thank Eric Slater, Lloyd Green, Mike Goldin, and Chris Neely of the Ocean Physics Group at SIO for their help in the design, construction, deployment, and operation of the sensors used during MBL. The experiment would not have been possible without the expertise and cooperation of the captain, Tom Golfinos, and crew of FLIP. Discussions with Mike Gregg, Chip Cox, Eric D’Asaro, Myrl Hendershott, and Lawrence Armi were helpful.

## REFERENCES

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

*J. Phys. Oceanogr.,***30,**805–832.Anderson, S. P., 1993: Shear, strain and thermohaline vertical fine structure in the upper ocean. Ph.D. thesis, University of California, San Diego. [Available from Scripps Institution of Oceanography, 9500 Gilman Drive, La Jolla, CA 92093-0208.].

——, and R. Pinkel, 1994: Double diffusively unstable intrusions near an oceanic front: Observations from R/P

*FLIP. Double-Diffusion Convection, Geophys. Monogr.,*No. 94, Amer. Geophys. Union, 195–211.Davis, R. E., 1996: Sampling turbulent dissipation.

*J. Phys. Oceanogr.,***26,**341–358.Dillon, T. M., 1982: Vertical overturns: A comparison of Thorpe and Ozmidov length scales.

*J. Geophys. Res.,***87,**9601–9613.——, and D. R. Caldwell, 1980: The Batchelor spectrum and dissipation in the upper ocean.

*J. Geophys. Res.,***85,**1910–1916.Gargett, A. E., and R. W. Schmitt, 1982: Observations of salt fingers in the central waters of the eastern North Pacific.

*J. Geophys. Res.,***87,**8017–8092.Garrett, C., and W. Munk, 1975: Space–time scales of internal waves:A progress report.

*J. Geophys. Res.,***80,**291–297.Gibson, C. H., and W. H. Schwarz, 1963: The universal equilibrium spectra of turbulent velocity and scalar fields.

*J. Fluid Mech.,***41,**153–167.Gregg, M. C., 1975: Microstructure and intrusions in the California Current.

*J. Phys. Oceanogr.,***5,**253–278.——, 1980: The three-dimensional mapping of a small thermohaline intrusion.

*J. Phys. Oceanogr.,***10,**1468–1492.——, 1987: Diapycnal mixing in the thermocline: A review.

*J. Geophys. Res.,***92,**5249–5286.——, 1989: Scaling turbulent dissipation in the thermocline.

*J. Geophys. Res.,***94,**9686–9698.——, and J. H. MacKenzie, 1979: Thermohaline intrusions lie across isopycnals.

*Nature,***280,**310–311.——, E. A. D’Asaro, T. J. Shay, and N. Larson, 1986: Observations of persistent mixing and near-inertial internal waves.

*J. Phys. Oceanogr.,***16,**856–884.Howard, L. N., 1961: Note on a paper of John W. Miles.

*J. Fluid Mech.,***10,**509–512.Joyce, T. M., W. Zenk, and J. M. Toole, 1978: The anatomy of the Antarctic Polar Front in the Drake Passage.

*J. Geophys. Res.,***83,**6093–6113.Linden, P. F., 1971: Salt fingers in the presence of grid-generated turbulence.

*J. Fluid. Mech.,***49,**611–624.Lueck, R. G., and J. J. Picklo, 1990: Thermal inertia of conductivity cells: Observations with a Sea-Bird cell.

*J. Atmos. Oceanic Technol.,***7,**756–768.Mack, S. A., and H. C. Schloberlein, 1993: Discriminating salt fingering from turbulence-induced microstructure: Analysis of towed temperature–conductivity chain data.

*J. Phys. Oceanogr.,***23,**2073–2106.Meagher, T. B., A. M. Pederson, and M. C. Gregg, 1982: A low-noise conductivity microstructure instrument.

*Oceans ‘82: Conf. Record,*Washington, DC, Marine Technology Society, 283–290.Miles, J. W., 1961: On the stability of heterogeneous shear flows.

*J. Fluid Mech.,***10,**496–508.Morison, J., R. Anderson., N. Larson, and E. D’Asaro, 1994: The correction for thermal-lag effects in Sea-Bird CTD data.

*J. Atmos. Oceanic Technol.,***11,**1151–1164.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,**256–271.Osborn, T. R., and C. S. Cox, 1972: Oceanic fine structure.

*Geophys. Fluid Dyn.,***3,**321–345.Pinkel, R., and J. A. Smith, 1992: Repeat sequence coding for improved precision of Doppler sonar and sodar.

*J. Atmos. Oceanic Technol.,***9,**149–163.Polzin, K. L., J. M. Toole, and R. W. Schmitt, 1995: Finescale parameterizations of turbulent dissipation.

*J. Phys. Oceanogr.,***25,**306–327.Ruddick, B. R., 1983: A practical indicator of the stability of the water column to double-diffusive activity.

*Deep-Sea Res.,***30,**1105–1107.——, 1992: Intrusive mixing in a Mediterranean salt lens: Intrusion slopes and dynamical mechanisms.

*J. Phys. Oceanogr.,***22,**1274–1285.——, and J. S. Turner, 1979: The vertical length scale of double-diffusive intrusions.

*Deep-Sea Res.,***26,**903–913.Schmitt, R. W., 1994: Double diffusion in oceanography.

*Annu. Rev. Fluid Mech.,***26,**255–285.Seim, H. E., and M. C. Gregg, 1994: Detailed observations of a naturally occurring shear instability.

*J. Geophys. Res.,***99,**10 049–10 073.Stern, M. E., 1967: Lateral mixing of water masses.

*Deep-Sea Res.,***14,**747–753.St. Laurent, L., and R. W. Schmitt, 1999: The contribution of salt fingers to vertical mixing in the North Atlantic Tracer Release Experiment.

*J. Phys. Oceanogr.,***29,**1404–1424.Turner, J. S., 1978: Double-diffusive intrusions into a density gradient.

*J. Geophys. Res.,***83,**2887–2901.

## APPENDIX

### Estimating Dissipation Rate from Underresolved Conductivity Gradient Fluctuations

Dissipation rate estimates are often based on underresolved shear measurements. In situations where the measured fluctuations exhibit an inertial subrange, then ε can be inferred from the spectral level in that range. The assumption is that the unresolved scales exhibit predicted spectral forms.

Underresolved scalar gradient data cannot be used in this manner since the spectral level in the inertial subrange is determined by both ε and *χ.* Though ε and *χ* are related through the physics of stratified turbulence, they are independent parameters in the Batchelor spectrum formulation. As a result, neither can be determined precisely without fully resolved spectra.

Nonetheless, approximate ε estimates can be made from the MBL microconductivity data. If the fluctuations are described by the Batchelor spectrum for scalars, then the fraction of total conductivity variance resolved by the probe may be expressed as a function of ε (as ε increases, more variance is shifted to higher wavenumbers, and the fraction at resolvable wavenumbers decreases). Assuming quasi-steady conditions, a constant mixing efficiency and constant density ratio, ε is a simple function of *χ.* The value of ε consistent with both constraints may then be determined. Observed correspondence between ε_{μ} and ε_{T} during double-diffusively stable periods (Fig. 14) implies that in the depth mean, the numerous assumptions made are reasonable (or at least the errors cancel out).

*φ*

_{TG}(

*k*) are plotted, given by

*k** =

*c*

_{∗}Pr

^{1/2}

*k*

_{B}

*ν*

*K*

^{−1}

_{T}

*c*∗ = 0.0968 is a constant (Gibson and Schwarz 1963). The Batchelor subrange

*α*= (2

*q*)

^{1/2}

*k*/

*k*

_{B}. Here

*q*is a universal constant, taken to be 2(3)

^{1/2}(as argued theoretically by Gibson and Schwarz 1963). The parameter

*β*=

*qc*

^{2/3}

_{∗}

Representative Batchelor spectra are plotted in Fig. A1. Each has equal variance (*χ* ≡ 6*K* ^{∞}_{0}*φ*_{TG}(*k*) *dk* = 10^{−8} °C^{2} s^{−1} for all curves). As ε increases, *k*_{B} increases, shifting more variance out to unresolved scales. Integrating out to the probe’s Nyquist wavenumber *k*_{N} = 13.3 cpm yields measured *χ̂**K* ^{kN}_{0}*T*(*k*)*φ*_{TG}(*k*) *dk.* Here *T*(*k*) is the probe response function (section 2b), taken to be *T*(*k*) = sinc^{2}(*πk*/*k*_{N}) (if Sea-Bird’s antialiasing response function (Fig. 3) is used instead, there is little difference). We define *r*(ε) ≡ *χ̂**χ* as the ratio of turbulent temperature gradient variance at wavenumbers less than *k*_{N} (which are resolved by the probe) to the total.

*r*(ε) (Fig. A2) is well fitted by a power law

*r*

*c*

^{p}

^{−5}

^{−0.4}

^{−1}

*χ̂,*

*χ̂*

_{pe}

The density ratio *R*_{ρ} represents the relative contributions of temperature and salinity to density. Strictly, *χ*_{pe} should be estimated using in situ *R*_{ρ}. However, to avoid introducing too much complexity, a constant value of *R*_{ρ} = −1 is used, consistent with the 12-day mean value.

*χ̂*

_{pe}

*χ̂*

*χ̂*

_{pe}

*χ̂*

*c*= 1.99 × 10

^{−5}(if ε is expressed in W kg

^{−1}), and

*p*= −0.4.

*χ*

_{pe}. On average, the buoyancy flux

*J*

_{b}is expressible in terms of kinetic energy dissipation and a mixing efficiency Γ (taken to be 0.2, Oakey 1982):

*J*

_{b}

*J*

_{b}= ½

*χ*

_{pe}(Seim and Gregg 1994). Thus ε is determined in terms of

*χ*

_{pe}:

*χ̂*

_{pe}

*c*

^{−1/(p+1)}

*χ̂*

^{1/(p+1)}

_{pe}

*χ̂*

*p*= 0,

*c*= 1) (A7) reduces to the familiar (A5), which results from the Osborn–Cox (1972) assumptions, with

*K*

_{T}=

*K*

_{ρ}and

*J*

_{b}= Γε.

The method assumes that the Batchelor spectrum obtains. In practice, the Batchelor subrange in temperature gradient spectra is usually observed in high Cox number situations (Dillon and Caldwell 1980). Additional assumptions are made, that is, that reasonable estimates of *χ*_{pe} are given by using in situ *N*^{2} but cruise-mean *R*_{ρ}, that *χ*_{pe} yields the buoyancy flux (i.e., quasi-steady state), and that buoyancy flux and ε are related by a constant mixing efficiency Γ = 0.2. In spite of all these assumptions, remarkable agreement (Fig. 14) is obtained between ε_{μ} (which involves centimeter-scale microconductivity fluctuations) and ε_{T} (which involves 2-m overturns and greater). The tight agreement suggests that, during periods where double-diffusive processes are unimportant, both methods produce reliable depth-averaged dissipation rate estimates.

Double-diffusive percentages.

Conditionally sampled variables.

^{1}

A simple test of the validity of this model is possible. Alford and Pinkel (2000) observe that |∂*ŵ*/∂*z*| is reduced immediately following overturning events, relative to values in advance. Since overturning represents irreversible processes, the observed differences imply that overturning alters the ∂*ŵ*/∂*z* field. Lacking a preferred spatial direction, some variability in ∂*ŵ*/∂*z* is therefore due to intrinsic temporal evolution. The same technique is applied to the *Ĉ* and *χ̂**χ̂**Ĉ* should be higher following overturns. However, we detect no difference in *χ̂**Ĉ* after overturns relative to before them. Apparently, intrinsic temporal evolution in *χ̂,**Ĉ* is more difficult to detect than in ∂*ŵ*/∂*z,* given their longer decay times.

^{2}

That is, the observed layer slope does not represent sinking with time, but rather downward spatial slope away from the interface of a previously formed structure. Were the mean currents reversed, the structure would appear inverted, but the sense of cold, fresh water sloping downward into warmer water would remain.

^{3}

Within a statically unstable region, |Tu| > 90°. Calculation of Tu on isopycnals, which are computed from sorted density profiles and therefore cannot cross, removes this possibility. Here Tu “at” an overturning point represents the relative *T–S* contributions to the background (4-m) stratification that the overturn is overcoming.