An official website of the United States government
Many turbulence estimates require fitting expected model forms, typically nonlinear equations, to observations that have been mathematically converted into spectral observations. This database contains synthetic spectral observations used to evaluate several fitting methods, along with the estimated parameters. The impact of degrees of freedom (spectral averaging) was assessed using two datasets with variability generated using two different statistical distributions. The testing was done as part of ATOMIX SCOR working group #160, with support from NSF grant #OCE-2140395 and contributions from national SCOR committees. The ATOMIX wiki has more information about the group's activities. A manuscript has been submitted relating to the above testing, and this repo will be updated once the citation is available.
more »
« less
This content will become publicly available on September 16, 2026
Assessing statistical fitting methods used for estimating turbulence parameters
Abstract Many turbulence estimates require fitting model forms, typically nonlinear expressions, to observations that have been converted into spectra. Choosing a fitting method usually depends on user preference, such as calculation ease under the spectra's presumed statistical nature or reducing computation demands when turbulence quantities must be estimated onboard expendable instruments. Six different methods are assessed by fitting a known model against synthetic spectra with variability generated from two different statistical distributions. The assessment uses an inertial subrange model to estimate the turbulent kinetic energy dissipation rate from velocity spectra. However, the results and conclusions are relevant to fitting other turbulence inertial subrange models that follow a power law where is the spectral slope and contains the sought‐after turbulence parameter. The two most accurate methods require linearizing the spectral observations by taking the logarithm of the wavenumbers and the dependent spectra power density . These methods are less sensitive to outliers and deviations of the observations from a known statistical distribution. Some methods returned that deviated from the prescribed value by more than 50% depending on the number of samples fitted and the level of uncertainty of the spectra. Methods for estimating the spectral slope, , were also assessed to provide recommendations on using this parameter to flag data which deviates from the expected form so that the spectra (or wavenumbers) can be excluded from further analysis.
more »
« less
- Award ID(s):
- 2140395
- PAR ID:
- 10652498
- Publisher / Repository:
- ASLO
- Date Published:
- Journal Name:
- Limnology and Oceanography: Methods
- ISSN:
- 1541-5856
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The inertial subrange of turbulent scales is commonly reflected by a power law signature in ensemble statistics such as the energy spectrum and structure functions – both in theory and from observations. Despite promising findings on the topic of fractal geometries in turbulence, there is no accepted image for the physical flow features corresponding to this statistical signature in the inertial subrange. The present study uses boundary layer turbulence measurements to evaluate the self-similar geometric properties of velocity isosurfaces and investigate their influence on statistics for the velocity signal. The fractal dimension of streamwise velocity isosurfaces, indicating statistical self-similarity in the size of ‘wrinkles’ along each isosurface, is shown to be constant only within the inertial subrange of scales. For the transition between the inertial subrange and production range, it is inferred that the largest wrinkles become increasingly confined by the overall size of large-scale coherent velocity regions such as uniform momentum zones. The self-similarity of isosurfaces yields power-law trends in subsequent one-dimensional statistics. For instance, the theoretical 2/3 power-law exponent for the structure function can be recovered by considering the collective behaviour of numerous isosurface level sets. The results suggest that the physical presence of inertial subrange eddies is manifested in the self-similar wrinkles of isosurfaces.more » « less
-
Abstract Stratification can cause turbulence spectra to deviate from Kolmogorov's isotropicpower law scaling in the universal equilibrium range at high Reynolds numbers. However, a consensus has not been reached with regard to the exact shape of the spectra. Here we propose a shape of the turbulent kinetic energy and temperature spectra in horizontal wavenumber for the equilibrium range that consists of three regimes at small Froude number: the buoyancy subrange, a transition region, and the isotropic inertial subrange through dimensional analysis and substantial revision of previous theoretical approximation. These spectral regimes are confirmed by various observations in the atmospheric boundary layer. The representation of the transition region in direct numerical simulations will require large‐scale separation between the Dougherty‐Ozmidov scale and the Kolmogorov scale for strongly stratified turbulence at high Reynolds numbers, which is still challenging computationally. In addition, we suggest that the failure of Monin‐Obukhov similarity theory in the very stable atmospheric boundary layer is due to the fact that it does not consider the buoyancy scale that characterizes the transition region.more » « less
-
Horizontal and vertical wavenumbers (kx, kz) immediately below the Ozmidov wavenumber (N3/ε)1/2 are spectrally distinct from both isotropic turbulence (kx, kz > 1 cpm) and internal waves as described by the Garrett–Munk (GM) model spectrum (kz < 0.1 cpm). A towed CTD chain, augmented with concurrent Electromagnetic Autonomous Profiling Explorer (EM-APEX) profiling float microstructure measurements and shipboard ADCP surveys, are used to characterize 2D wavenumber (kx, kz) spectra of isopycnal slope, vertical strain, and isopycnal salinity gradient on horizontal wavelengths from 50 m to 250 km and vertical wavelengths of 2–48 m. For kz < 0.1 cpm, 2D spectra of isopycnal slope and vertical strain resemble GM. Integrated over the other wavenumber, the isopycnal slope 1D kx spectrum exhibits a roughly +1/3 slope for kx > 3 × 10−3 cpm, and the vertical strain 1D kz spectrum a −1 slope for kz > 0.1 cpm, consistent with previous 1D measurements, numerical simulations, and anisotropic stratified turbulence theory. Isopycnal salinity gradient 1D kx spectra have a +1 slope for kx > 2 × 10−3 cpm, consistent with nonlocal stirring. Turbulent diapycnal diffusivities inferred in the (i) internal wave subrange using a vertical strain-based finescale parameterization are consistent with those inferred from finescale horizonal wavenumber spectra of (ii) isopycnal slope and (iii) isopycnal salinity gradients using Batchelor model spectra. This suggests that horizontal submesoscale and vertical finescale subranges participate in bridging the forward cascade between weakly nonlinear internal waves and isotropic turbulence, as hypothesized by anisotropic turbulence theory.more » « less
-
A model for the structure function tensor is proposed, incorporating the e↵ect of anisotropy as a linear perturbation to the standard isotropic form. The analysis extends the spectral approach of Ishihara et al. (2002) to physical space based on Kolmogorov’s theory and is valid in the inertial range of turbulence. Previous results for velocity co-spectra are used to obtain estimates of the model coe"cients. Structure functions measured from direct numerical simulations of channel flow and from experimental measurements in turbulent boundary layers are compared with predicted behaviour and reasonable agreement is found. We note that power-law scaling is more evident in the co-spectra than for the mixed structure functions. New observations are made about countergradient correlation between Fourier modes of wall normal and streamwise velocity components for wavenumbers approaching the Kolmogorov scale.more » « less
