Iterative spectral integration algorithm

Spectral integration is an iterative procedure because the bandwidth required to estimate the variance of shear depends on the rate of dissipation, the level of electronic noise in the shear-probe signal, on the wavenumber of spurious signals that were not removed by the Goodman algorithm, and on the wavenumber resolution of the shear probe. The rate of dissipation is estimated using

[math]\displaystyle{ \begin{equation} \varepsilon = \frac{15}{2} \nu \int^{k_u}_{k_0} S(k) dk \end{equation} }[/math]

where [math]\displaystyle{ S(k) }[/math] is the dimensional shear spectrum. The upper limit of integration, [math]\displaystyle{ k_u }[/math], is et by the smallest of a number of criteria.

Electronic noise usually manifests itself as a minimum in the spectrum with real shear at wavenumbers smaller than this minimum and electronic noise at higher wavenumbers. The minimum may be found by fitting a low-order polynomial to the spectrum in log-log space. The wavenumber of the spectral minimum sets one of the limits on [math]\displaystyle{ k_u }[/math].

Another upper limit is [math]\displaystyle{ k_u = 150 cpm }[/math] that is imposed by the spatial resolution of the shear probe. The cut-off frequency, [math]\displaystyle{ f_A }[/math], of the anti-aliasing filter used by the shear-probe sampler sets another upper limit of spectral integration, namely [math]\displaystyle{ k_u \leq f_A/U }[/math]. The user may impose an upper limit to exclude wavenumbers that have contaminations that are not correctable. All of these limits are fixed for a particular data set.

The wavenumber range of the spectrum of shear depends on the rate of dissipation. The spectrum broadens in proportion to [math]\displaystyle{ \epsilon^{1/4} }[/math] and the spectrum rises in proportion to [math]\displaystyle{ \epsilon^{3/4} }[/math]. Thus, the fraction of the shear variance that is resolved at any particular wavenumber depends on the rate of dissipation.

However, the non-dimensional spectrum [math]\displaystyle{ G(\hat k) = \frac{L^2_k}{(\varepsilon \nu^5)^{1/4}} S(k) }[/math]

where [math]\displaystyle{ \hat k=kL_k }[/math] and [math]\displaystyle{ L_k=(\nu^3/\varepsilon)^{1/4} }[/math] s the Kolmogorov length, is expected to be universal and independent of the rate of dissipation. There are several analytic models of the non-dimensional spectrum that are based on approximations of empirically derived spectra such as Wolk et al (2002) who provides an approximation to the Nasmyth spectrum, and Lueck (2021b) who provides an improved approximation to the Nasmyth spectrum and to a non-dimensional spectrum derived from 14,000 dimensional spectra. These spectral approximations also provide the fraction of the shear variance that is resolved at any particular non-dimensional wavenumber. For all of these analytic approximations, 95% of the variance of shear is resolved at a non-dimensional wavenumber of approximately 0.12. There is little motivation for integrating the spectrum beyond this wavenumber. This sets another upper limit to spectral integration once the rate of dissipation is known, or at least known approximately.

The upper limit of spectral integration is nearly always higher than 10 cpm. For example, when [math]\displaystyle{ \varepsilon=1\times 10^{-10} }[/math] [math]\displaystyle{ W }[/math] [math]\displaystyle{ kg^{-1} }[/math], this limit resolves a little more than 95% of the shear variance. Less for larger dissipation rates. Thus, if we integrate a shear spectrum to 10 cpm this variance must be related to the total variance if the spectrum follows a universal form. In fact, if the measured spectrum follows the universal form exactly, then the ratio of the “true” dissipation rate, [math]\displaystyle{ \varepsilon }[/math], to the rate derived from integration to 10 cpm, [math]\displaystyle{ \varepsilon_{10} }[/math], is given by

[math]\displaystyle{ \frac{\varepsilon}{\varepsilon_{10}}= \sqrt{1+a\varepsilon_{10}} }[/math]

where [math]\displaystyle{ a= 1.08 \times 10^9 }[/math] Thus, integrating the spectrum to 10 cpm provides the first (and quite rough) estimate of the rate of dissipation and this estimate can be used to set another upper limit to spectral integration. That is, at what wavenumber is 95% of the spectrum resolved if the dissipation rate equals this initial estimate?

The spectrum is then integrated to the lowest of these upper limits to provide an improved estimate of the rate of dissipation. It is an improved estimate because it will nearly always span a range of wavenumbers that is considerably wider than 10 cpm, which improves its statistical reliability. This improved estimate is then used to determine the fraction of the shear variance that is resolved by the upper limit of integration. The estimate is divided by this fraction to form a new estimate of [math]\displaystyle{ \varepsilon }[/math] This process is then repeated -- find the new fraction resolved, adjust the estimate upwards using this fraction, and so on. This iteration converges quickly because the wavenumber is non-dimensionalized by [math]\displaystyle{ \varepsilon^{1/4} }[/math]. If the fraction of the variance resolved by a particular [math]\displaystyle{ k_u }[/math] exceeds 50%, the estimate of the rate of dissipation converges to within 1% of its ultimate value in two or fewer iterations.

For very high rates of dissipation, such as [math]\displaystyle{ \varepsilon=1\times 10^{-5} }[/math] [math]\displaystyle{ W }[/math] [math]\displaystyle{ kg^{-1} }[/math], the shear-probe cannot fully resolve the spectrum of shear. We recommend to estimate the rate by fitting to the spectrum in the inertial subrange.

In this range, the spectrum rises in proportion to [math]\displaystyle{ \varepsilon^{2/3}k^{1/3} }[/math] and, thus, its level provides an estimate of the rate of dissipation. The inertial subrange is confined to wavenumbers smaller than [math]\displaystyle{ k(\nu^3/\varepsilon)^{1/4}=0.02 }[/math], and thus will usually use fewer spectral points than the method of spectral integration. This reduces the statistical reliability of the dissipation estimate but it does avoid the bias introduced by not fully resolving the spectrum of shear.

[math]\displaystyle{ }[/math]