the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
A turbulence data reduction scheme for autonomous and expendable profiling floats
Kenneth G. Hughes
James N. Moum
Daniel L. Rudnick
Autonomous and expendable profilingfloat arrays such as those deployed in the Argo Program require the transmission of reliable data from remote sites. However, existing satellite data transfer rates preclude complete transmission of rapidly sampled turbulence measurements. It is therefore necessary to reduce turbulence data on board. Here we propose a scheme for onboard data reduction and test it with existing turbulence data obtained with a modified SOLOII profiling float. First, voltage spectra are derived from shear probe and fastthermistor signals. Then, we focus on a fixedfrequency band that we know to be unaffected by vibrations and that approximately corresponds to a wavenumber band of 5–25 cpm. Over the fixedfrequency band, we make simple power law fits that – after calibration and correction in postprocessing – yield values for the turbulent kinetic energy dissipation rate ϵ and thermalvariance dissipation rate χ. With roughly 1 m vertical segments, this scheme reduces the necessary data transfer volume 300fold to approximately 2.5 kB for every 100 m of a profile (when profiling at 0.2 m s^{−1}). As a test, we apply our scheme to a dataset comprising 650 profiles and compare its output to that from our standard turbulenceprocessing algorithm. For ϵ, values from the two approaches agree within a factor of 2 87 % of the time; for χ, they agree 78 % of the time. These levels of agreement are greater than or comparable to that between the ϵ and χ values derived from two shear probes and two fast thermistors, respectively, on the same profiler.
Measurements of oceanic turbulence have been made since the 1950s using platforms and sensors of various shapes and sizes (Lueck et al., 2002). Complete resolution of the turbulence requires measuring temperature and velocity gradients at millimetertocentimeter scale. Hence, sampling turbulence is data intensive. Whereas conventional profiling measurements of temperature, conductivity, and pressure are typically sampled at 1 Hz (e.g., Argo floats; Roemmich et al., 2019a), a turbulence profile involves sampling multiple sensors at 100–1000 Hz. A relatively minimal requirement of five separate signals sampled at 100 Hz and recorded at 16bit resolution equates to 1 kB s^{−1} or 500 kB per 100 m of profiling range at 0.2 m s^{−1} profiling speed. For floats, this is not a trivial volume of data. For example, transmitting only 3 kB of data from a Deep SOLO float takes 100–200 s (Roemmich et al., 2019b). Extended surfacings are also presented with the danger of surface vessels and vandals. Ultimately, raw turbulence profiles are 2–3 orders of magnitude too large to transmit in a reasonable amount of time.
One approach to reducing turbulence data is given by Rainville et al. (2017), who use it for multimonth glider missions. On board the glider, they calculate spectra of raw voltage signals reported by the shear probes and fast thermistors and then band average each of these spectra into 12 bins. After transmission, these binned values are calibrated and fit to model spectra. In other words, they (i) postpone calibration and (ii) minimize the data manipulation and processing that happens on board. These two strategies are shared by our scheme (and also shared by the reduction scheme developed for χpods; Becherer and Moum, 2017). Otherwise, however, our scheme differs from that of Rainville et al. (2017), as it does not suit our scientific goals of measuring turbulence over the upper ∼ 120 m at high vertical resolution (e.g., ∼ 1 m) and as frequently as possible (see Sect. 2). Two shear probes and two fast thermistors would produce $(\mathrm{2}+\mathrm{2})\times \mathrm{12}=\mathrm{48}$ spectral values per segment. Even without considering the other profiling quantities, the spectral values could add up to >20 kB per dive for our scenario. We would be spending as much time transmitting the data as actually measuring the ocean.
Our scheme is for profiling instruments that contain shear probes and, optionally, fast thermistors (Sect. 2). First, we document the necessary calibration details (Sect. 3). Next, we compress raw shear voltages by way of simple power law fits and show how ϵ is derived from these fits in postprocessing (Sect. 4). A test of the scheme employing 650 profiles demonstrates that little accuracy is sacrificed in return for a large reduction in data volume (Sect. 5). Knowing ϵ from the shear probe measurements makes possible a similar method for deriving χ from fastthermistor measurements (Sects. 6 and 7). Adapting the scheme to a different profiler requires minimal modification (Sect. 8). For our setup, the scheme reduces the dataset size by a factor of ∼ 300: only 2.5 kB for each 100 m of a profile (Sect. 9).
We intend our data reduction scheme to be sufficiently general so as to be portable to all vertical turbulence profilers that contain shear probes. It can also be used with gliders if a measure of flow speed past the sensors is available (e.g., Greenan et al., 2001; Merckelbach and Carpenter, 2021). In a general sense, some of the values specified herein ought to be considered variables (Sect. 8). However, we do have a particular platform for which we are developing the scheme, the Flippin' χSOLO (FCS), and the values used here are chosen for the objective of detailed upperocean profiling.
FCS is a descendant of the SOLOII profiling float (Roemmich et al., 2004) with the addition of a turbulence package and extra functionality. The turbulence package includes two shear probes (Osborn, 1974) to measure smallscale velocity gradients from which ϵ is computed, two fast thermistors to measure smallscale temperature fluctuations from which χ is computed, and a pressure sensor from which profiling speed is derived. FCS also includes a threeaxis accelerometer that is used to measure the surface wave field (although with a method not described in this paper). Accelerometer data from when FCS is profiling are not used in our reduction scheme. FCS and its measurements are described more completely in the companion paper (Moum et al., 2023).
To reverse profiling direction, FCS adjusts buoyancy and flips via internal shifting of the battery pack. This causes the turbulence sensors to always point into undisturbed fluid. Flipping therefore permits profiling on both descent and ascent, including sampling of the upper 5 m on the ascent.
As a prototype, a standard (nonflipping) SOLO float with a modified χpod (Moum and Nash, 2009) attached was deployed in the Bay of Bengal to measure the suppression of turbulence by salinity stratification (Shroyer et al., 2016). This instrument – named χSOLO – did not have shear probes and therefore could not have provided estimates of ϵ within mixed layers. (Values of ϵ can be approximated from χ, but only if there is stratification.) Nevertheless, χSOLO's success motivated the development of the FCS units with flipping capabilities and fully integrated turbulence packages. These new instruments retained the χ or C in their name despite their ability to also directly measure ϵ from shear probes.
Two FCS units were vetted over 4 days in May 2019 off the Oregon coast. During this period, each unit profiled from the surface to ∼ 120 m and back at a typical speed of 0.2 m s^{−1}. Adding time at the surface, each dive cycle took ∼ 30 min, and we obtained 650 profiles in total. In this 2019 experiment, one of the shear probes on one of the two units malfunctioned. Hence, the dataset for this paper contains approximately 25 % less shear data than fastthermistor data.
The core of our data reduction scheme uses power law fits of voltage spectra that are calculated on board and subsequently converted to meaningful turbulence quantities in postprocessing. Additional voltage quantities are also recorded to determine temperature, pressure, and profiling speed.
3.1 Nomenclature and conventions

All quantities measured by FCS that are discussed in this paper are sampled at 100 Hz.

All voltage spectra are frequency spectra and are denoted by Ψ_{x}(f) (where x is a label, such as s for shear) with units of V^{2} Hz^{−1}.

Physical spectra of shear and temperature gradient are wavenumber spectra and are denoted by Φ_{x}(k) with units of s^{−2} cpm^{−1} and K^{2} m^{−2} cpm^{−1}, respectively. Figure 4 is an exception in which shear spectra are frequency spectra: Φ_{s}(f).

Wavenumber k has the unit cycles per meter (cpm). Expressions quoted from other papers may differ by factors of 2π for wavenumbers in radians per meter.

The Kraichnan model spectrum Φ_{Kr} primarily depends on the dissipation rates of turbulent kinetic energy and temperature variance (ϵ and χ), but it also depends on the molecular viscosity ν and molecular thermal diffusivity D_{T}. For brevity, we write ${\mathrm{\Phi}}_{\mathrm{Kr}}(k,\mathit{\u03f5},\mathit{\chi})$ rather than the more complete ${\mathrm{\Phi}}_{\mathrm{Kr}}(k,\mathit{\u03f5},\mathit{\chi},\mathit{\nu},{D}_{\mathrm{T}})$. Similarly, the Nasmyth spectrum is written as Φ_{Na}(k,ϵ) rather than ${\mathrm{\Phi}}_{\mathrm{Na}}(k,\mathit{\u03f5},\mathit{\nu})$. In cases where the arguments are unambiguous or unimportant, we simply write Φ_{Na} and Φ_{Kr}.

A pair of angle brackets, 〈⋅〉, denotes the mean value over a segment of length: N_{seg}=512 data points. This equates to ∼ 1 m at our nominal profiling speed of 0.2 m s^{−1}.

To calculate spectra for a given 512element voltage segment, we first remove the linear trend, then we use three halfoverlapping, Hammingwindowed, 256element subsegments (i.e., N_{fft}=256; N_{overlap}=128).
In general, the values of N_{seg} and N_{fft} are variables. Our choices are based on the 100 Hz sampling (∼ 500 cpm) and the goals of FCS, which include obtaining highverticalresolution turbulence data, especially near the surface. For different turbulence profilers or different scientific goals, longer segments and/or more overlapping subsegments may be more appropriate (see Sect. 8).
We do not pursue the possibility of using accelerometers to decontaminate spectra (e.g., Levine and Lueck, 1999); three subsegments is too few for this to work well. Rather, we focus on a frequency band that we know to be unaffected by vibration.
3.2 Shear calibration
The voltage reported by the shear probe V_{s} is linearly proportional to shear:
where W is the flow speed past the sensor. The overall engineering calibration α includes the seawater density ρ, the analog circuit gain G_{s} (equal to 1 for FCS circuitry), the probe sensitivity S_{s} (∼ $\mathrm{0.25}\times {\mathrm{10}}^{\mathrm{3}}$ V m^{2} N^{−1}), and the differentiator time constant T_{s} (∼ 1 s). The linearity in V_{s} admits a simple link between the physical and voltage spectra:
where ${H}_{\mathrm{s}}^{\mathrm{2}}\left(k\right)$ is the transfer function that accounts for (i) spatial averaging by the shear probe of highwavenumber motions and (ii) analog and digital filtering of the raw voltage signal (see Appendix A). Note also the use above of the following relation:
3.3 Temperature and temperature gradient calibration
Two voltage signals are recorded for each fast thermistor. V_{T} is the voltage output directly related to T, and V_{Tt} is the differentiated output, which improves resolution at high frequencies (≳ 10 Hz). Temperature is related to V_{T} through a quadratic calibration:
where C_{1T}, C_{2T}, and C_{3T} are coefficients determined from laboratory calibrations. Equation (6) is technically an approximation because it contains 〈V_{T}〉^{2} and not $\langle {V}_{T}^{\mathrm{2}}\rangle $, but over 5 s time scales this changes 〈T〉 by ≲ 0.001 ^{∘}C (estimated from the 2019 dataset).
The gradient of the temperature calibration is
Consequently, the smallscale vertical temperature gradient T_{z} is linearly proportional to the differentiated voltage V_{Tt}. To demonstrate, we first rewrite T_{z} in terms of more directly measured quantities:
The first quantity on the righthand side is Eq. (7), the last is $\mathrm{1}/W$, and the second is
where C_{Tt} is the gain of the analog differentiator.
Rewriting Eq. (8), the aforementioned linear relationship between T_{z} and V_{Tt} becomes
The relationship between physical and voltage spectra is therefore
Again, we have invoked Eq. (4), and the transfer function ${H}_{T\phantom{\rule{0.125em}{0ex}}t}^{\mathrm{2}}\left(k\right)$ is defined in Appendix A.
3.4 Pressure and profilingvelocity calibration
Pressure has a linear calibration:
In our usage, the coefficients C_{1P} and C_{2P} are recorded in units of psi and psi V^{−1}, respectively, and calibrated under total pressure. Subtracting atmospheric pressure makes P=0 at the sea surface. The constant C_{1P} must account for the vertical position of the pressure sensor on the instrument relative to the shear probes and thermistors. Hence, C_{1P} differs between upcasts and downcasts. For the reduced dataset, we record the last pressure voltage in each segment. For example, with N_{seg}=512, we save the 512th, 1024th, and 1536th values of V_{P} for the first three segments. The average pressures in the second and third segments are 0.5(V_{P}(1024)+V_{P}(512)) and 0.5(V_{P}(1536)+V_{P}(1024)), respectively. The average pressure in the first segment is found by extrapolation.
The flow speed past the sensors, denoted W, is derived from the rate of change of the pressure voltages just described:
where Δt is the sampling period (here 0.01 s), and ΔP is V_{P}(1024)−V_{P}(512) and V_{P}(1536)−V_{P}(1024) for the second and third segments, respectively. Extrapolation is again used for the first segment.
No smoothing is necessary before calculating ΔV_{P} because its magnitude is so much larger than the quantization of the signal (this being the limiting factor for precision of pressure recorded by FCS). In physical units, P is precise to 0.003 dbar, which is 𝒪(300) times smaller than ΔP.
Wave orbitals can introduce variability when W is small (≲ 0.15 m s^{−1}). As a diagnostic we calculate and record the minimum value of W for each segment. This also helps to identify the beginning and end of profiles, as shown in Appendix B. In standard processing, we would derive W(t) from the pressure signal low passed at 2 Hz. To avoid the need to lowpass filter the signal on board, we instead make 10 estimates of W(t) per segment and take the minimum of these:
where t_{i}=1, 51, 101, …, 501. Even with this sampling of every 50th element, which follows from subsampling a 100 Hz signal at 2 Hz, ΔV_{P}(t_{i}) is large enough that smoothing is unnecessary.
In this paper, we immediately discard all segments in which W_{min}<0.05 m s^{−1}. This threshold is reached only at the top and bottom of profiles, if at all. Note, however, that this does not imply that a segment with W_{min}>0.05 m s^{−1} is trustworthy. Even segments with W_{min} closer to 0.15 m s^{−1} should be treated with particular caution. Signs that a segment is questionable are that 〈W〉 and W_{min} differ by ∼ 20 % or more and that spectral fit scores are low (see Sects. 4.3 and 6.3). These two issues often cooccur because of the nonlinear relationship between shear and profiling speed (Eq. 1). Lowfrequency variations in W ultimately lead to spectra that are redder than expected and hence have low fit scores. Such segments should be discarded. There is not a simple way to correct the spectra given the nonlinearity.
In this section, we are ultimately going to fit measured spectra to an inertial subrange model that does not necessarily apply at the relevant frequencies or wavenumbers. We will elaborate as we go, but we want to emphasize in advance that measured spectra do not need to conform to an inertial subrange model for us to obtain accurate values of ϵ. The inertial subrange is merely a convenient starting point.
4.1 Summarizing Nasmyth spectra with ${f}^{\mathrm{1}/\mathrm{3}}$ fits
Shear measurements ideally capture both the inertial and viscous subranges and hence use a wide band of the measured spectrum to derive values for ϵ. In practice, noise and sensor resolution limit how well the true environmental spectrum is resolved. Conventional workarounds exploit the Nasmyth model spectrum Φ_{Na}(k,ϵ) (Nasmyth, 1970; Oakey, 1982). One approach is to iterate toward a solution in which the integral of Φ_{Na} over a specific wavenumber band matches that of the measured spectrum ${\mathrm{\Phi}}_{{u}_{z}}$ (e.g., Moum et al., 1995). Another is to find the best fit of ${\mathrm{\Phi}}_{{u}_{z}}$ to Φ_{Na} by using maximumlikelihood estimation together with a model of the expected statistical distribution of the spectral coefficients being fitted (e.g., Bluteau et al., 2016).
Here we develop a new and simpler twostage approach to fitting shear spectra to Φ_{Na}. In the first stage, we use an ${f}^{\mathrm{1}/\mathrm{3}}$ power law fit over a fixedfrequency range of f_{l} to f_{h} = 1–5 Hz, where ${f}^{\mathrm{1}/\mathrm{3}}$ follows from the assumption that we are fitting over the inertial subrange. In the second stage, we correct for this ofteninvalid assumption.
In the inertial subrange, shear spectra are proportional to ${k}^{\mathrm{1}/\mathrm{3}}$ and hence also to ${f}^{\mathrm{1}/\mathrm{3}}$, since f=Wk. With N_{fft}=256 and 100 Hz sampling (Sect. 3.1), spectral coefficients are separated by frequency increments of 100 Hz $/$ 256 = 0.39 Hz, so there are 10 coefficients between 1 and 5 Hz. (Our processing code will actually use bounding frequencies of 0.98 and 4.88 Hz, as these are halfinteger multiples of 0.39 Hz, but for brevity we will write these as 1 and 5 Hz throughout.)
Our choice of f_{l}=1 Hz is dictated by a requirement that we avoid lowfrequency contamination induced by (i) advection by wave orbital motion and (ii) pitch and roll motions of the profiler. Together, these dominate below 0.3 Hz. Setting f_{l}=0.5 Hz would add only one more spectral coefficient. Our choice of f_{h}=5 Hz is a tradeoff between maximizing the bandwidth of the fit and minimizing how much the measured spectra are subject to either noise or viscous roll off. Other profilers may benefit from different frequency bounds (see Sect. 8).
Our inertial subrange assumption is often false. Indeed, “assumption” is perhaps a misnomer, as we do not expect it to be true; we know that viscous roll off will often occur at frequencies lower than 5 Hz (25 cpm for a nominal value of W=0.2 m s^{−1}). However, because there exists an analytical expression for the viscous roll off, we are able to derive an exact expression that quantifies how much ϵ is underestimated. This is the second stage of our approach. We derive an expression for the correction function F_{Na} in such a way that it can be calculated in postprocessing. The benefits of this approach are that (i) we can fit uncalibrated (i.e., voltage) spectra, and (ii) it simplifies the actual onboard fitting routine (Sect. 4.2).
The full Nasmyth spectrum and its inertial range approximation are as follows (Lueck, 2013):
where $\mathit{\eta}=({\mathit{\nu}}^{\mathrm{3}}/\mathit{\u03f5}{)}^{\mathrm{1}/\mathrm{4}}$ is the Kolmogorov length scale.
Let ϵ_{init} denote the initial value of ϵ that comes from fitting a measured spectrum to the approximate form in Eq. (16) using the simple power law fitting method in Appendix C rather than fitting to the full form in Eq. (15). As noted earlier, the fit will be over the f_{l}–f_{h} = 1–5 Hz range which, given our nominal value of W=0.2 m s^{−1}, equates to k_{l}–k_{h} = 5–25 cpm.
Consider two contrasting examples of low and high turbulence with $\mathit{\u03f5}=\mathrm{1}\times {\mathrm{10}}^{\mathrm{9}}$ and $\mathrm{1}\times {\mathrm{10}}^{\mathrm{6}}$ W kg^{−1}, respectively (Fig. 1a). For now, assume the measured spectrum to be fit is itself a Nasmyth spectrum. For $\mathit{\u03f5}={\mathrm{10}}^{\mathrm{6}}$ W kg^{−1}, the ${f}^{\mathrm{1}/\mathrm{3}}$ fit lies on top of Φ_{Na}. Conversely, the ${f}^{\mathrm{1}/\mathrm{3}}$ fit for the smaller ϵ value is seemingly meaningless: the ${f}^{\mathrm{1}/\mathrm{3}}$ fit (dashed line) does not even match the sign of the slope of Φ_{Na}. Worse yet, naively inverting this initial fit produces the underestimate ${\mathit{\u03f5}}_{\mathrm{init}}=\mathrm{1.2}\times {\mathrm{10}}^{\mathrm{10}}$ W kg^{−1}, which is 8 times smaller than the true value of ϵ. However, by adjusting by a factor of $\mathrm{1}/{F}_{\mathrm{Na}}$, defined in the following paragraph, the fit (dotted line) looks like a hypothetical extrapolation of the inertial subrange. Equivalently, ϵ_{init} is corrected to the true value of ϵ as
In our example, $\mathrm{1}\times {\mathrm{10}}^{\mathrm{9}}$ W kg^{−1} = $\mathrm{1.2}\times {\mathrm{10}}^{\mathrm{10}}$ W kg^{−1} $/$ ${\mathrm{0.238}}^{\mathrm{3}/\mathrm{2}}$. The value of 0.238 is the solution to an implicit equation, derived below, that depends on ϵ_{init} and W. For clarity, our demonstration starts by assuming that we know ϵ rather than ϵ_{init}.
Nasmyth spectra can be flattened to unity over the inertial subrange with the normalization $\mathrm{8.05}{k}^{\mathrm{1}/\mathrm{3}}{\mathit{\u03f5}}^{\mathrm{2}/\mathrm{3}}$ (Fig. 1b). Values of F_{Na} are based on the mean of these flattened spectra over the wavenumber range k_{l}–k_{h} ($={f}_{\mathrm{l}}/W$–${f}_{\mathrm{h}}/W$):
To remove the dependence of the true value of ϵ, we substitute Eq. (17) to produce an implicit function for F_{Na}, which can be solved numerically:
Note how the two forms of F_{Na} (Eqs. 18 and 19) are defined with different arguments. For our example, F_{Na} ($\mathit{\u03f5}=\mathrm{1}\times {\mathrm{10}}^{\mathrm{9}}$ W kg^{−1}) = ${F}_{\mathrm{Na}}({\mathit{\u03f5}}_{\mathrm{init}}=\mathrm{1.2}\times {\mathrm{10}}^{\mathrm{10}}$ W kg^{−1}) = 0.238. Hereafter, we use the latter: F_{Na}(ϵ_{init}).
With f_{l} and f_{h} fixed, F_{Na} is a function of three variables: ϵ_{init}, W, and ν. F_{Na} is closer to 1 (less of a correction) for larger values of ϵ_{init} (Fig. 1). It is also closer to 1 for higher values of W (Fig. 2), since k_{l} and k_{h} decrease with increasing W (i.e., k_{l}–k_{h} move closer to the inertial subrange).
To simplify calculations in the upcoming section, we make one final change to Eq. (19) using the following substitution:
Therefore,
Think of this substitution in Eq. (20) as inverting the conventional way that ${H}_{\mathrm{s}}^{\mathrm{2}}\left(k\right)$ is invoked. Usually, a measured shear spectrum is amplified at high wavenumbers by $\mathrm{1}/{H}_{\mathrm{s}}^{\mathrm{2}}\left(k\right)$ and is then fit to the model spectrum Φ_{Na}. Here, instead of amplifying the measured spectrum, we reduce the model spectrum. With this latter approach, ${H}_{\mathrm{s}}^{\mathrm{2}}\left(k\right)$ is calculated and applied only during the postprocessing stage. (It changes F_{Na} by only ∼ 5 %, since we fit over relatively low wavenumbers.)
4.2 Obtaining ϵ_{init} from shear voltage spectra
Since ϵ can be reconstructed from ϵ_{init}, we require an expression linking ϵ_{init} to the shear voltage spectrum Ψ_{s}. Equating Eqs. (3) and (16) gives
where we have left out ${H}_{\mathrm{s}}^{\mathrm{2}}\left(k\right)$, since it has been incorporated into F_{Na}. Rearranging and substituting $k=f/W$ gives
Then, to solve for ϵ_{init}, we use a leastsquares fit (see Appendix C):
where the sums are understood to be over the range f_{l}–f_{h}. The quantities α, W, and ϵ_{init} are calculated in postprocessing.
4.3 Quality control of the shear spectral fits
Measured shear spectra are often quality controlled either by manual visual inspection or, more objectively, by quantifying the level of mismatch between them and their associated model. Possible mismatch quantities include the mean absolute deviation or the variance of the ratio ${\mathrm{\Phi}}_{{u}_{z}}/{\mathrm{\Phi}}_{\mathrm{Na}}$ (e.g., Ruddick et al., 2000; Bluteau et al., 2016). We cannot calculate such quantities with our reduced scheme because we do not know what each spectrum should look like until we calculate its ϵ value in the postprocessing stage. (Recall that Φ_{Na} is a function of ϵ.) By this stage, we have lost information about the spectral shape through the summing operation in Eq. (24).
To retain at least some information about the shape of each voltage spectrum, we will split the 1–5 Hz range and compute two fits rather than one. Doing so allows for a firstorder check that the spectrum over the 1–5 Hz range approximately follows the expected shape.
Mathematically, there is nothing special about our choice f_{l}–f_{h}=1–5 Hz. In theory, we can split the 1–5 Hz range into two (1–3 and 3–5 Hz) and obtain a value of ϵ_{init} for each. These values will differ, but so will the associated values of F_{Na}. For a measured spectrum that conforms to a Nasmyth spectrum, the two values of ϵ calculated with Eq. (17) will not differ (Fig. 3). We therefore calculate on board the sums in Eq. 24 over both f_{l}–f_{m} and f_{m}–f_{h}, where the middle frequency f_{m}=3 Hz. (In our code, f_{m} is actually 7.5×0.39 Hz = 2.93 Hz for the reason given in Sect. 4.1.) Hence, for each spectrum, we are able to postprocess to recover two independent estimates of ϵ, denoted ϵ_{l–m} and ϵ_{m–h}. The mean of these two provides a single, final value for ϵ, and their ratio quantifies the match of a measured spectrum to a Nasmyth spectrum over the range f_{l}–f_{h}:
The best possible fit score is 1; the lower the score, the poorer the fit. The example spectra in Fig. 4 show that a high fit score does not necessarily imply small residuals. Rather, fits with high scores are typically those with random residuals, meaning that a given measured spectral coefficient is just as likely to be above the fit as it is to be below it. Fits with low fit scores are typically those with autocorrelated residuals, meaning that the sign and/or magnitude of a residual is correlated with that of its neighbors. In practice, we expect a range of ϵ fit scores: instantaneous and unaveraged spectra differ from the Nasmyth spectrum because they are derived from a limited sampling of a statistical process and because of nonstationarity, anisotropy, and inhomogeneity of the turbulence.
When ϵ is small (≲ 10^{−9} W kg^{−1}), the fit score may be consistently low if spectral coefficients in the f_{m}–f_{h} range are affected by noise, and consequently ϵ_{m–h}≫ϵ_{l–m}. For such cases, we choose to use only the lowerfrequency fit. We would rather have a more accurate estimate of ϵ and forgo the fit score than have a biasedhigh ϵ value with a biasedlow fit score. (Either way, the small values of ϵ in question will have minimal effect on any averages, given that turbulence distributions have high kurtosis, so high values dominate means.) Specifically,
where the threshold is equivalent to $k<\mathrm{0.1}/\mathit{\eta}$, with η being the Kolmogorov length scale estimated from ϵ_{l–m}. For reference, Φ_{Na} peaks at $k=\mathrm{0.026}/\mathit{\eta}$ and rolls off to 11 % of its maximum by $k=\mathrm{0.1}/\mathit{\eta}$ (see Eq. 15).
To test the accuracy of the shear reduction scheme described in the previous section, we apply it retrospectively to the dataset from the 2019 test cruise (Sect. 2). We compare the results to those obtained with the standard processing scheme. This standard scheme (Appendix D) features a more sophisticated despiking routine than that used for our reduced scheme, which employs a 3standarddeviations threshold filter (Appendix E).
A profilebyprofile comparison of the two schemes is shown in Fig. 5. The comparison is then extended to all 650 profiles (> 77 000 segments of shear), where we find that ϵ from the reduced scheme (${\mathit{\u03f5}}_{\mathrm{init}}/{F}_{\mathrm{Na}}^{\mathrm{3}/\mathrm{2}}$) is within a factor of 2 of that from the standard scheme 87 % of the time over the full range of measured values, ${\mathrm{10}}^{\mathrm{10}}<\mathit{\u03f5}<{\mathrm{10}}^{\mathrm{4}}$ W kg^{−1} (Fig. 6a–b). For comparison, in only 72 % do we obtain a factorof2 agreement between the two independent values of ϵ measured on the unit with two working shear probes (not shown). Further, to obtain this 87 % agreement, we clearly need the correction function F_{Na}: Fig. 6c shows that the uncorrected values ϵ_{init} only have 1 : 1 agreement with ϵ from the standard scheme if $\mathit{\u03f5}\gtrsim {\mathrm{10}}^{\mathrm{7}}$ W kg^{−1}. For the lowest values of ϵ, the ratio is closer to 1 : 30.
To demonstrate the ability of the ϵ fit score to characterize spectra, we show twodimensional histograms of nondimensionalized spectral coefficients from all 77 000 measured shear spectra separated into three classes based on their ϵ fit score: 0.67–1.00, 0.33–0.67, and 0.00–0.33. Only the lowestscoring class fails to collapse to the Nasmyth spectrum (Fig. 7c, f).
The scheme to reduce fastthermistor data to enable measurement of χ is much like the scheme to reduce shear data. As in Sect. 4, we first show how we summarize a model spectrum in terms of a power law fit and a correction factor. (In this case, the correction factor partly depends on the values of ϵ calculated in Sect. 4.) Then we derive the implementation in terms of voltages and calculate a spectralfit metric.
6.1 Summarizing Kraichnan spectra with f^{1} fits
Here we take the Kraichnan spectrum Φ_{Kr} (Kraichnan, 1968) as our model; for its lowwavenumber approximation, we use the viscous–convective subrange, which scales as k^{+1}. In units of K^{2} m^{−2} cpm^{−1}, Φ_{Kr} and its approximation are as follows (e.g., Peterson and Fer, 2014):
where the Batchelor length scale ${\mathit{\lambda}}_{B}=(\mathit{\nu}{D}_{\mathrm{T}}^{\mathrm{2}}/\mathit{\u03f5}{)}^{\mathrm{1}/\mathrm{4}}$, and q is a constant taken to be 5.26. This expression does not include a ${k}^{+\mathrm{1}/\mathrm{3}}$ inertial–convective subrange, which we ignore here, as it increases by less than 1 % the integral of the temperature gradient spectrum from k=0 to k=∞ and therefore has negligible effect on our results.
A fit against Eq. (29) can be rearranged to give χ_{init}, which is related to χ through the correction function F_{Kr} as
F_{Kr} is not raised to a power like F_{Na} (Eq. 17). For small values of k, Φ_{Kr}∝χ, whereas ${\mathrm{\Phi}}_{\mathrm{Na}}\propto {\mathit{\u03f5}}^{\mathrm{2}/\mathrm{3}}$.
The derivation of F_{Kr} is equivalent to F_{Na}. We therefore present only the result:
Note that ${F}_{\mathrm{Kr}}(\mathit{\u03f5},{\mathit{\chi}}_{\mathrm{init}},W)$ depends on the underestimate χ_{init}, whereas it depends on the “true” or “corrected” value of ϵ as calculated in Sect. 4.
6.2 Obtaining χ_{init} from fastthermistor voltage spectra
Like we did for ϵ_{init} in Sect. 4.2, we derive the expression for χ_{init} in three steps. First, equate the righthand sides of Eqs. (11) and (29) (excluding the transfer function ${H}_{T\phantom{\rule{0.125em}{0ex}}t}^{\mathrm{2}}\left(k\right)$, which is incorporated into Eq. 31):
Then, rearrange while substituting $k=f/W$ to get
Finally, solve for χ_{init} using a leastsquares fit (Appendix C):
6.3 Quality control of the temperature gradient spectral fits
The approach to quality controlling the fastthermistor data is the same as that for shear (Sect. 4.3). That is, we fit Ψ_{Tt} over f_{l}–f_{m} and f_{m}–f_{h} (1–3 and 3–5 Hz). This ultimately provides two estimates of χ for each spectrum, which are combined as follows:
We do not apply a low χ threshold equivalent to Eq. (27).
Profiles of χ from the reduced scheme compare well to the standard processing, albeit with a small bias in one direction for low values and in the other direction for high values (Fig. 8). Across all values, the two approaches agree within a factor of 2 78 % of the time (Fig. 9). By comparison, 82 % of segments exhibit a factorof2 agreement between χ values from the two fast thermistors on the same unit.
Compared to shear spectra, nondimensionalized temperature gradient spectra have lower fit scores. Especially for the lowest fit scores, the measured temperature gradient spectra tend to be too high at lower frequencies and too low at frequencies near f_{h} (Fig. 10a–c). As frequency increases beyond f_{h}, the effects of noise and thermalresponse corrections (Appendix A) begin to dominate.
There are three reasons for the poorer fits to temperature gradient spectra compared to that for shear. First, shapes of temperature gradient spectra are often more variable; the best choice for the nondimensional spectral model can be debated (e.g., Sanchez et al., 2011). Second, the temperature gradient fits depend on ϵ, so uncertainties in ϵ propagate into the calculation of χ. Third, for our 2019 experiment, the recorded temperature gradient signals were occasionally affected by digitization noise as a consequence of sampling mixed layers. (Shear signals were not affected by digitization noise.)
8.1 Setting the scheme's parameters
Our scheme requires a few userdefined parameters: f_{l}, f_{h}, N_{seg}, and N_{fft}. For this paper, we based these partly on the profiling speed and scientific goals of FCS. For a different profiler, we suggest the following:

Choose f_{h} based on a typical profiling speed such that ${k}_{\mathrm{h}}={f}_{\mathrm{h}}/W\approx \mathrm{25}$ cpm for a nominal profiling speed W. For a wide range of ϵ values, 25 cpm is close to, or beyond, the peak of the Nasmyth spectrum (Fig. 1). As an example, if FCS profiled at ∼ 0.5 m s^{−1}, we would consider setting f_{h}≈12 Hz.

Ensure that there are no known issues such as vibrations that are likely to adversely affect spectral coefficients within the f_{l}–f_{h} range. If, however, there are known issues within the desired frequency range, then an alternative approach (one that we did not test) is to use accelerometer signals to correct spectra that are contaminated by vibrations (Levine and Lueck, 1999; Goodman et al., 2006).

Define f_{l} and f_{h} separately for shear and temperature gradient if appropriate. Although we set them to be equal here, this is not necessary.

Use more than two fitting bands if desired. We use only two bands (1–3 and 3–5 Hz) so as to minimize the file size to be transmitted, but there is nothing preventing there being three or more bands (e.g., adding a 5–7 Hz band). Indeed, this would enable improved estimates of the fit scores (Eqs. 26 and 36) and more flexibility to discard noiseaffected bands as in Eq. (27). If file size is less of an issue such that it is possible to send back fit values for many more than two bands, then the Rainville et al. (2017) scheme outlined in Sect. 1 may be a better choice than ours.

Choose N_{seg} and N_{fft} based on scientific goals and, possibly, any logistical constraints; the data reduction scheme is agnostic to these numbers. For example, at the expense of vertical resolution, we could halve the file size of our transmitted dataset by doubling N_{seg} from 512 to 1024.

Reasonable choices for N_{fft} are ${N}_{\mathrm{seg}}/\mathrm{2}$ or ${N}_{\mathrm{seg}}/\mathrm{4}$, which correspond to three or seven halfoverlapping subsegments, respectively. There is little to be gained by dividing a segment into even more subsegments so as to produce smoother spectra before fitting. As Ruddick et al. (2000) note, the choice is analogous to fitting a line to 20 points at once or first clumping them into groups of, say, five and then fitting the four averaged points.
8.2 Evaluating the reduced data
One step that cannot be automated is the heuristic evaluation of the reduced turbulence data after they have been converted from voltage quantities to physical ones. For this evaluation, we recommend looking into multiple quantities. First consider the fit scores (Eqs. 26 and 36). We recommend discarding any ϵ or χ values with an associated fit score lower than 0.33. Note, however, that these scores are not a perfect measure of fit. They should be used together with other quality control checks such as comparing

W and W_{min} (Eqs. 13 and 14) to check whether the profiling speed is constant over a segment;

ϵ values from the two shear probes; and

turbulent features in successive profiles.
The last point is most applicable for a vertical profiler cycling rapidly – for example, twice per hour for FCS. In this case, the profiler is nominally sampling the same vertical fragment of the ocean on a timescale comparable to that over which turbulence evolves. In our experience, many turbulent patches extend over 5–10 profiles.
Recall, also, that all uncertainty in ϵ propagates into the calculation of χ (Sect. 7). If ϵ for a given segment cannot be trusted, neither can χ.
We have developed a data reduction scheme applicable to vertical profiling of turbulence variables in which each ∼ 5 s segment is distilled to 12 quantities (Scheme 1). In postprocessing, we reconstruct estimates of ϵ and χ, associated quality control metrics, and other quantities such as the temperature and profiling speed. The raw data that go into the 12 quantities are seven different voltages (V_{P}; V_{T} and V_{Tt} for each thermistor; and V_{s} for each shear probe). Hence, for each 512element segment, we effectively reduce the data by a factor of $\mathrm{512}\times \mathrm{7}/\mathrm{12}\approx \mathrm{300}$.
This reduction compresses the output data file size for each dive from megabytes to kilobytes. For example, the total amount of data per dive (two profiles) can be estimated assuming our nominal dive depth and profiling velocity of 120 m and 0.2 m s^{−1}. Each dive creates $\mathrm{12}\times \mathrm{2}\times \mathrm{120}$ m $/$ (0.2 m s${}^{\mathrm{1}}\times \mathrm{5.12}\phantom{\rule{0.125em}{0ex}}\mathrm{s})\approx \mathrm{2800}$ quantities. Transmitting each quantity as a 16bit float or integer equates to approximately 6 kB per dive. This can be reduced by onethird if the spectralfit metrics are suitably scaled logarithmically and then transmitted as 8bit integers.
One luxury we lose is the ability to inspect the raw signals. Typically this would help to (i) cultivate faith in the data, (ii) flag which segments to discard, and (iii) inform workarounds such as filtering out potential narrowband vibrations in shear spectra. Our scheme accounts for this constraint in two ways. First, we fit spectra over relatively low frequencies (1–5 Hz) that are unlikely to be affected by noise or vibration. Second, we reduce the data in a way that uses as little arithmetic as possible. Obviously, we cannot reverseengineer the raw signals, but by making the onboard calculations simple, we give ourselves the best chance to later fix or identify any unforeseen issues.
Although the onboard reduction eliminates possibilities in how we process turbulence data, it opens up possibilities in how we obtain turbulence data. By visualizing how turbulence evolves over successive dives in nearreal time, we can concentrate on regions of interest by adapting the dive schedule to profile more frequently or to different depths. If instead we encounter quiescent periods, we might consider profiling less frequently, thereby conserving battery life. Our ultimate objective is to treat FCS floats as expendable.
Voltage signals from shear probes and thermistors are a smoothed representation of the true environmental signal. If the smoothing is a spatial effect, it is described by a transfer function H^{2}(k). If the smoothing is a temporal effect, it is more natural to use H^{2}(f). We can use these interchangeably because f=Wk, and therefore H^{2}(f)=H^{2}(Wk). For FCS, there are three components to the transfer function for each sensor:
where we have used the following shorthand: SP is shear probe, FT is fast thermistor, AA is antialiasing, and D is digital. We describe each of these in turn.
Shear probes built and calibrated by the Ocean Mixing Group are very close in dimension to those examined by Ninnis (1984), who measured their wavenumber response and represented it as
where a_{0}=1.000, ${a}_{\mathrm{1}}=\mathrm{0.164}$, ${a}_{\mathrm{2}}=\mathrm{4.537}$, a_{3}=5.503, ${a}_{\mathrm{4}}=\mathrm{1.804}$, and k_{0}=170 cpm.
Temporal averaging of temperature at high frequencies due to the thermal response of the fast thermistor is modeled using a doublepole filter:
where the cutoff frequency f_{c}=30 Hz. This comes from Nash et al. (1999), who measured the frequency response for two different thermistors on an instrument profiling at 0.3 m s^{−1} and found cutoff frequencies of 25.1 and 36.7 Hz (see their Fig. A2). The 30 Hz value is the approximate mean of these two values.
Raw shear and thermistor voltage signals are both subject to two filters. First, an analog antialiasing filter (twopole Butterworth) with an f_{c}=40 Hz cutoff:
After the analog signal is antialiased, it is digitized at 400 Hz. Before subsampling to the final 100 Hz output, the signal is digitally filtered. For the 2019 FCS cruise, the signal was convolved with a symmetric 29element kernel in which the first 15 elements were
This is a sinc kernel but with negative values set to zero. (We are currently investigating better choices for future implementations). The filter has a halfpower (−3 db) point at 25 Hz.
Early in our processing routine, we partition the raw voltage signals into 512element segments. In order to discard the segments in which FCS was not profiling, we need robust (yet simple) criteria that demarcate the start and end of a profile. For the start, we search for the first three consecutive segments in which W_{min}>0.05 m s^{−1}. For the end, we swap the inequality.
A drawback of this approach is the appearance of a quantity in physical units (0.05 m s^{−1}). This is the one instance where we hard code a calibration coefficient in the onboard software rather than applying it in postprocessing. Fortunately, the relevant coefficient can be approximated as constant: C_{2P} = 76.7 psi V^{−1} (barring a redesign of the circuitry or the use of a different brand or model of pressure sensor). For the two units already built, C_{2P}=76.81 psi V^{−1} and 76.53 psi V^{−1}. By comparison, among the four shear probes on the two units, the calibration coefficients vary by 30 %.
At least for the initial implementation of our scheme, we do not include an algorithm to detect the surface to within centimeters. Doing so would let us work backward to put our uppermost depth bin as close to the surface as possible. However, we expect that this could be a fragile part of the scheme. Further, FCS lacks a microconductivity sensor, which is likely the sensor best suited for identifying the air–sea interface (e.g., Ward et al., 2014).
Without surface detection, the depths of the uppermost bins will be realized randomly. In the worst cases, we would discard the top ∼ 1 m (5 s at ∼ 0.2 m s^{−1}). To alleviate this, we may use halfoverlapping bins near the surface. The exact implementation will be determined later in the development.
In this paper, we use power law fits to derive turbulence quantities: ${\mathrm{\Psi}}_{\mathrm{s}}={A}_{\mathit{\u03f5}}{f}^{\mathrm{1}/\mathrm{3}}$ and Ψ_{Tt}=A_{χ}f^{1}, where A_{ϵ} and A_{χ} are substitutes for the expressions in Eqs. (23) and (33). With only a single parameter for each fit, implementing a leastsquares fit is easy.
Assume we are fitting the vector Ψ_{i} to the function $A{f}_{i}^{n}$ where n is either $\mathrm{1}/\mathrm{3}$ or 1. The sum of squared residuals is therefore
The minimum with respect to A is where the derivative is zero:
Hence,
We had originally intended to find A by following Becherer and Moum (2017), who were fitting ${f}^{\mathrm{1}/\mathrm{3}}$ spectra. Their simpler method, $A=\sum \left({\mathrm{\Psi}}_{i}/{f}_{i}^{n}\right)$, is equivalent to a leastsquares fit, except the quantity minimized is the sum of the squares of the adjusted residuals, where adjusted means divided by f^{n}. Differences can be ignored when $n=\mathrm{1}/\mathrm{3}$ but not when n=1.
The standard processing of FCS turbulence data differs from the reduced scheme in three ways. First, raw data are despiked differently (Appendix E). Second, the 100 Hz raw voltage signals are calibrated into physical quantities right away. Hence, means and spectra are calculated in physical units and not in voltage units. Third, the integration of spectra occurs over a variable wavenumber band, which is found iteratively.
When integrating shear spectra (after correction; see Appendix A) to find ϵ, we follow the approach used for the Chameleon profiler (Moum et al., 1995). A first estimate of ϵ is made by integrating over k=4–10 cpm. This value provides a first estimate of the Kolmogorov wavenumber ${k}_{\mathrm{s}}=(\mathit{\u03f5}/{\mathit{\nu}}^{\mathrm{3}}{)}^{\mathrm{1}/\mathrm{4}}/\mathrm{2}\mathit{\pi}$. (The lower limit for Chameleon is 2 cpm, but we increase this for FCS given its slower profiling speed and hence the possibility of contamination by waves at lower wavenumbers.) The upper integral limit is then set to 0.5k_{s} (with a minimum of 10 cpm and a maximum of 45 cpm). The Nasmyth spectra (Eq. 15) is integrated over the same wavenumber range. If the measured and Nasmyth integrals are within 1 %, then ϵ is set to be equal to the integral of the Nasmyth spectrum over all k. Otherwise, ϵ and k_{s} are adjusted iteratively until the two integrals agree.
A similar but noniterative approach is used for integrating T_{z} spectra to find χ. The model spectrum is the Kraichnan spectrum (Eq. 28), and again, the lower limit of integration is 4 cpm. The upper limit is the Batchelor wavenumber ${k}_{b}=(\mathit{\u03f5}/\mathit{\nu}{D}_{\mathrm{T}}^{\mathrm{2}}{)}^{\mathrm{1}/\mathrm{4}}/\mathrm{2}\mathit{\pi}$ (with a maximum defined by kW=15 Hz).
To properly despike the raw output of a shear probe requires several steps. Lueck et al. (2018) describe a process in which the signal is highpassed, then rectified, and then lowpassed to derive a measure of the local variance. A value is defined as a spike if it is more than 8 times (or a similar threshold) higher than the local variance. Spikes are replaced with an average based on surrounding points. This process is then repeated on the new signal and so on until no spikes are identified.
In our standard processing of FCS data, we use the Lueck et al. (2018) despiking routine. For our data reduction scheme, we use an approach that is easier to implement and quicker to compute, albeit one that is less precise. For each 512element segment of data, a spike is defined as any data point larger than 3 standard deviations from the mean. These spikes are replaced by the mean of the remaining values in the segment.
Our MATLAB implementation of the processing code is available from https://doi.org/10.5281/zenodo.7644701 (Hughes and Vutukur, 2023) or https://github.com/OceanMixingGroup/flippinchisolo (last access: 17 February 2023).
Raw and processed data for the 2019 experiment are available at https://doi.org/10.5281/zenodo.5719505 (Hughes, 2022) or https://kghughes.com/data (last access: 17 February 2023).
KGH designed the reduction scheme and led the writing of the paper. All authors contributed to the final version. JNM and DLR lead the development of the FCS profiler on which much of the paper is based.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This work was funded by the Office of Naval Research under grant nos. N000141712700 (OSU) and N000141712762 (SIO) and continued as part of the ARCTERX (Island Arc Turbulent Eddy Regional Exchange) program under grant nos. N000142112878 (OSU), N000142112762, and N000142112747 (SIO). Engineers who contributed to the design and construction of FCS and its sensors are Craig Van Appledorn, Kerry Latham, Pavan Vutukur, and Mark Borgerson (all from OSU) and Ben Reineman, Kyle Grindley, and Jeff Sherman from SIO. Aurélie Moulin executed initial turbulence processing, and Emily Shroyer provided many helpful comments on early drafts. Thanks to the reviewers, including Cynthia Bluteau and Toshiyuki Hibiya.
This research has been supported by the Office of Naval Research (grant nos. N000141712700, N000141712762, N000142112878, N000142112762, N000142112747).
This paper was edited by Katsuro Katsumata and reviewed by Toshiyuki Hibiya and Cynthia Bluteau Toshiyuki Hibiya and one anonymous referee.
Becherer, J. and Moum, J. N.: An efficient scheme for onboard reduction of moored χpod data, J. Atmos. Ocean. Tech., 34, 2533–2546, https://doi.org/10.1175/JTECHD170118.1, 2017. a, b
Bluteau, C. E., Jones, N. L., and Ivey, G. N.: Estimating turbulent dissipation from microstructure shear measurements using maximum likelihood spectral fitting over the inertial and viscous subranges, J. Atmos. Ocean. Tech., 33, 713–722, https://doi.org/10.1175/JTECHD150218.1, 2016. a, b
Goodman, L., Levine, E. R., and Lueck, R. G.: On measuring the terms of the turbulent kinetic energy budget from an AUV, J. Atmos. Ocean. Tech., 23, 977–990, https://doi.org/10.1175/JTECH1889.1, 2006. a
Greenan, B. J. W., Oakey, N. S., and Dobson, F. W.: Estimates of dissipation in the ocean mixed layer using a quasihorizontal microstructure profiler, J. Phys. Oceanogr., 31, 992–1004, https://doi.org/10.1175/15200485(2001)031<0992:EODITO>2.0.CO;2, 2001. a
Hughes, K.: A turbulence data reduction scheme for autonomous and expendable profiling floats: datasets, Zenodo [data set], https://doi.org/10.5281/zenodo.5719505, 2022. a
Hughes, K. and Vutukur, P.: OceanMixingGroup/flippinchisolo: Version 1, Version v1, Zenodo [code], https://doi.org/10.5281/zenodo.7644701, 2023. a
Kraichnan, R. H.: Smallscale structure of a scalar field convected by turbulence, Phys. Fluids, 11, 945–953, https://doi.org/10.1063/1.1692063, 1968. a
Levine, E. R. and Lueck, R. G.: Turbulence measurement from an autonomous underwater vehicle, J. Atmos. Ocean. Tech., 16, 1533–1544, https://doi.org/10.1175/15200426(1999)016<1533:TMFAAU>2.0.CO;2, 1999. a, b
Lueck, R.: Calculating the rate of dissipation of turbulent kinetic energy, Tech. rep., Rockland Scientific, https://rocklandscientific.com/support/knowledgebase/technicalnotes/ (last access: 2 August 2018), 2013. a
Lueck, R., Murowinski, E., and McMillan, J.: A guide to data processing: ODAS Matlab library v4.3., Tech. rep., Rockland Scientific, https://rocklandscientific.com/support/knowledgebase/technicalnotes/ (last access: 6 August 2021), 2018. a, b
Lueck, R. G., Wolk, F., and Yamazaki, H.: Oceanic velocity microstructure measurements in the 20th century, J. Oceanogr., 58, 153–174, https://doi.org/10.1023/A:1015837020019, 2002. a
Merckelbach, L. M. and Carpenter, J. R.: Ocean glider flight in the presence of surface waves, J. Atmos. Ocean. Tech., 38, 1265–1275, https://doi.org/10.1175/JTECHD200206.1, 2021. a
Moum, J. N. and Nash, J. D.: Mixing measurements on an equatorial ocean mooring, J. Atmos. Ocean. Tech., 26, 317–336, https://doi.org/10.1175/2008JTECHO617.1, 2009. a
Moum, J. N., Gregg, M. C., Lien, R. C., and Carr, M. E.: Comparison of turbulence kinetic energy dissipation rate estimates from two ocean microstructure profilers, J. Atmos. Ocean. Tech., 12, 346–366, https://doi.org/10.1175/15200426(1995)012<0346:COTKED>2.0.CO;2, 1995. a, b
Moum, J. N., Rudnick, D. L., Shroyer, E. L., Hughes, K. G., Reineman, B. D., Grindley, K., Sherman, J., Vutukur, P., Van Appledorn, C., Latham, K., Moulin, A. J., and Johnston, T. M. S.: Flippin' χSOLO, an upper ocean autonomous turbulence profiling float, J. Atmos. Ocean. Tech., in review, 2023. a
Nash, J. D., Caldwell, D. R., Zelman, M. J., and Moum, J. N.: A thermocouple probe for highspeed temperature measurement in the ocean, J. Atmos. Ocean. Tech., 16, 1474–1482, https://doi.org/10.1175/15200426(1999)016<1474:ATPFHS>2.0.CO;2, 1999. a
Nasmyth, P. W.: Oceanic turbulence, PhD thesis, Univ. British Columbia, https://doi.org/10.14288/1.0302459, 1970. a
Ninnis, R.: The effects of spatial averaging on airfoil probe measurements of oceanic velocity microstructure, PhD thesis, University of British Columbia, https://doi.org/10.14288/1.0053131, 1984. a
Oakey, N. S.: Determination of the rate of dissipation of turbulent energy from simultaneous temperature and velocity shear microstructure measurements, J. Phys. Oceanogr., 12, 256–271, https://doi.org/10.1175/15200485(1982)012<0256:DOTROD>2.0.CO;2, 1982. a
Osborn, T. R.: Vertical profiling of velocity microstructure, J. Phys. Oceanogr., 4, 109–115, https://doi.org/10.1175/15200485(1974)004<0109:VPOVM>2.0.CO;2, 1974. a
Peterson, A. K. and Fer, I.: Dissipation measurements using temperature microstructure from an underwater glider, Methods Oceanogr., 10, 44–69, https://doi.org/10.1016/j.mio.2014.05.002, 2014. a
Rainville, L., Gobat, J. I., Lee, C., and Shilling, G.: Multimonth dissipation estimates using microstructure from autonomous underwater gliders, Oceanography, 30, 49–50, https://doi.org/10.5670/oceanog.2017.219, 2017. a, b, c
Roemmich, D., Riser, S., Davis, R., and Desaubies, Y.: Autonomous profiling floats: Workhorse for broadscale ocean observations, Mar. Technol. Soc. J., 38, 21–29, https://doi.org/10.4031/002533204787522802, 2004. a
Roemmich, D., Alford, M. H., Claustre, H., Johnson, K., King, B., Moum, J., Oke, P., Owens, W. B., Pouliquen, S., Purkey, S., Scanderbeg, M., Suga, T., Wijffels, S., Zilberman, N., Bakker, D., Baringer, M., Belbeoch, M., Bittig, H. C., Boss, E., Calil, P., Carse, F., Carval, T., Chai, F., Conchubhair, D., d'Ortenzio, F., Dall'Olmo, G., Desbruyeres, D., Fennel, K., Fer, I., Ferrari, R., Forget, G., Freeland, H., Fujiki, T., Gehlen, M., Greenan, B., Hallberg, R., Hibiya, T., Hosoda, S., Jayne, S., Jochum, M., Johnson, G. C., Kang, K., Kolodziejczyk, N., Körtzinger, A., Traon, P.Y. L., Lenn, Y.D., Maze, G., Mork, K. A., Morris, T., Nagai, T., Nash, J., Garabato, A. N., Olsen, A., Pattabhi, R. R., Prakash, S., Riser, S., Schmechtig, C., Schmid, C., Shroyer, E., Sterl, A., Sutton, P., Talley, L., Tanhua, T., Thierry, V., Thomalla, S., Toole, J., Troisi, A., Trull, T. W., Turton, J., VelezBelchi, P. J., Walczowski, W., Wang, H., Wanninkhof, R., Waterhouse, A. F., Waterman, S., Watson, A., Wilson, C., Wong, A. P. S., Xu, J., and Yasuda, I.: On the future of Argo: A global, fulldepth, multidisciplinary array, Front. Mar. Sci., 6, 439, https://doi.org/10.3389/fmars.2019.00439, 2019a. a
Roemmich, D., Sherman, J. T., Davis, R. E., Grindley, K., McClune, M., Parker, C. J., Black, D. N., Zilberman, N., Purkey, S. G., Sutton, P. J. H., and Gilson, J.: Deep SOLO: A FullDepth Profiling Float for the Argo Program, J Atmos. Ocean. Tech., 36, 1967–1981, https://doi.org/10.1175/JTECHD190066.1, 2019b. a
Ruddick, B., Anis, A., and Thompson, K.: Maximum likelihood spectral fitting: The Batchelor Spectrum, J. Atmos. Ocean. Tech., 17, 1541–1555, https://doi.org/10.1175/15200426(2000)017<1541:MLSFTB>2.0.CO;2, 2000. a, b
Sanchez, X., Roget, E., Planella, J., and Forcat, F.: Smallscale spectrum of a scalar field in water: The Batchelor and Kraichnan models, J. Phys. Oceanogr., 41, 2155–2167, https://doi.org/10.1175/JPOD11025.1, 2011. a
Shroyer, E. L., Rudnick, D. L., Farrar, J. T., Lim, B., Venayagamoorthy, S. K., St. Laurent, L. C., Garanaik, A., and Moum, J. N.: Modification of upperocean temperature structure by subsurface mixing in the presence of strong salinity stratification, Oceanography, 29, 62–71, https://doi.org/10.5670/oceanog.2016.39, 2016. a
Ward, B., Fristedt, T., Callaghan, A. H., Sutherland, G., Sanchez, X., Vialard, J., and ten Doeschate, A.: The Air–Sea Interaction Profiler (ASIP): An autonomous upwardly rising profiler for microstructure measurements in the upper ocean, J. Atmos. Ocean. Tech., 31, 2246–2267, https://doi.org/10.1175/JTECHD1400010.1, 2014. a
 Abstract
 Introduction
 The Flippin' χSOLO (FCS)
 Conversion of measured voltages to physical units
 Reduction of shear data
 Test of the reduction scheme for ϵ
 Reduction of fastthermistor data
 Test of the reduction scheme for χ
 Recommendations
 Conclusions
 Appendix A: Transfer functions for FCS sensors
 Appendix B: Identifying the start and end of a profile
 Appendix C: Leastsquares fitting of power laws
 Appendix D: Standard processing of FCS turbulence measurements
 Appendix E: Identifying and removing noise and spikes in the shear signals
 Code availability
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References
 Abstract
 Introduction
 The Flippin' χSOLO (FCS)
 Conversion of measured voltages to physical units
 Reduction of shear data
 Test of the reduction scheme for ϵ
 Reduction of fastthermistor data
 Test of the reduction scheme for χ
 Recommendations
 Conclusions
 Appendix A: Transfer functions for FCS sensors
 Appendix B: Identifying the start and end of a profile
 Appendix C: Leastsquares fitting of power laws
 Appendix D: Standard processing of FCS turbulence measurements
 Appendix E: Identifying and removing noise and spikes in the shear signals
 Code availability
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References