Segmenting datasets

Once the raw observations have been quality-controlled, then you must split the time series into shorter segments by considering:

• Time and length scales of turbulence
• Stationarity of the segment and Taylor's frozen turbulence hypothesis
• Required statistical significance of the resulting spectra (important if you need to remove motion-induced contamination from the spectra)

Considerations

Measurements are typically collected in the following two ways:

• continuously, or in such long bursts that they can be considered continuous
• short bursts that are typically at most 2-3x the expected largest turbulence time scales (e.g., 10 min in ocean environments)

This segmenting step dictates the minimum burst duration when setting up your equipment. The act of chopping a time series into smaller subsets, i.e., segments, is effectively a form of low-pass (box-car) filtering. The length of the segment in time is usually a more important consideration than detrending the time series when estimating [math]\displaystyle{ \varepsilon }[/math] from the inertial subrange of the final spectra.

The shorter the segment, the higher the temporal resolution of the final [math]\displaystyle{ \varepsilon }[/math] time series, and the more likely the segment will be stationary. However, the spectrum's lowest resolved frequency and final resolution are the inverse of the fft-length, i.e., the duration of the signal used to construct the spectrum. Therefore, the segment must remain sufficiently long such that the lowest wavenumber (frequencies) of the inertial subrange are retained by the spectra. This is particularly important when measurement noise drowns the highest wavenumber (frequencies) of the inertial subrange. Thus, using too short segments may inadvertently render the spectra unusable for deriving [math]\displaystyle{ \varepsilon }[/math] from the inertial subrange by virtue of no longer resolving this subrange as shown in (Fig. 3)

Recommendations

A good rule of thumb for tidally-influenced environments is 5 to 15 min segments, but this may be shorter in certain energetic and fast-moving flows (Fig. 1) and longer in less energetic environments (Fig.2). The spectra are often estimated by block averaging numerous spectra (FFT) estimated from smaller chunks of data within each segment. Another strategy is band-averaging spectra in the frequency domain. The fft-length i.e., the duration of data used to estimate each spectrum dictates the final range of frequencies resolved by the spectra.

Fig. 3 provides a guide to how much of the fft-length to theoretically resolve different subrange as a function of the fft-length, speed past the sensor, and [math]\displaystyle{ \varepsilon }[/math]. For instance, an fft-length of 4 s would resolve one decade of the inertial subrange at speeds past the sensor of 0.5 m/s and [math]\displaystyle{ \varepsilon\sim10^{-7} }[/math] W/kg. Longer segments would be required for slower flows or lower [math]\displaystyle{ \varepsilon }[/math]. At [math]\displaystyle{ \varepsilon\approx10^{-9} }[/math] W/kg, one decade of the inertial subrange would be resolved with an fft-length longer than 10s provided the speed was faster than 0.5 m/s.

Because the inertial subrange may be contaminated at the highest wavenumbers by instrument noise, we suggest using longer segments than the minimum. This strategy also enables having a larger number of spectral observations to fit over the inertial subrange given the spectral resolution also depends on the fft-length.

The final segment length may be larger than the fft-length if using block averaging for the spectral computations. Is this explained in the spectral page?

Are the peaks in the MAVS data vortex shedding from the rings. Check the motion sensors onboard?

Fig. 1: Example theoretical velocity spectra for different [math]\displaystyle{ \varepsilon }[/math] with the empirical limit [math]\displaystyle{ \hat{k}L_k\sim0.1 }[/math] denoted by the diamonds ([math]\displaystyle{ \hat{k} }[/math] is in rad/m). The inertial subrange extends to smaller wavenumber [math]\displaystyle{ k }[/math] [cpm] as [math]\displaystyle{ \varepsilon }[/math] increases. The lowest frequency resolved by a spectra is the inverse of the fft-length used when computing the spectra. The colored lines are spectral observations from a dataset with fast speeds and large [math]\displaystyle{ \varepsilon }[/math]. In this example, we used relatively short segments (128s) to estimate the spectra from fft-length of 32 s (2048 samples @ 64 Hz). The impact of turbulence anisotropy is also visible through the flattening of the spectra around 1 cpm. The secondary x-axis show the corresponding frequencies for a range of mean speeds past the sensors

Fig. 2: Same as Fig 1 but for a different dataset with low speeds and low [math]\displaystyle{ \varepsilon }[/math], requiring the use of relatively long segments (1024s) to estimate the spectra from fft-length of 512 s (4096 samples @ 8 Hz).

Fig.3 Contours represent the log of the fft-length required to resolve the non-dimensional wavenumber [rad/m] indicated in each panel's title. The inertial subrange ends at approximately [math]\displaystyle{ \hat{k}L_k\approx0.1 }[/math] (or [math]\displaystyle{ kL_k\approx0.015 }[/math] in cpm), and so panel (c) denotes the fft-length that resolves the end of the inertial subrange i.e., the beginning of the viscous subrange. The fft-length must be at least 10x longer (see b), preferably 50x (panel c) given the low number of spectral observations at the lowest frequencies (wavenumbers)

---

Return to Preparing_quality-controlled_velocities