Synthetic spectra for assessing statistical fitting methods used to estimate ocean turbulence

Summary

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.

NetCDF-4 datasets

The repository contains the following datasets in grouped NetCDF-4 format. 

SyntheticSpectra.nc that contains:

3200 synthetic spectra with the variability drawn from lognormal distribution

3200 synthetic spectra with the variability drawn from chi2 distribution

LogNormal.nc that contains fitting results for lognormal dataset with four groups.

EPSI_DELTA_025 with the estimated (one paramter) epsilon using a decadal range of 0.25 to fit: 

10 samples in each spectrum

100 samples in each spectrum

EPSI_DELTA_1 with the estimated (one paramter) epsilon using a decadal range of 1 to fit: 

10 samples in each spectrum

100 samples in each spectrum

SLOPE_DELTA_025 with the (two-parameter) estimated spectral slope using a decadal range of 0.25 to fit: 

10 samples in each spectrum

100 samples in each spectrum

SLOPE_DELTA_1 with the (two-parameter) estimated spectral slope using a decadal range of 0.25 to fit: 

10 samples in each spectrum

100 samples in each spectrum

Chi.nc that contains fitting results for chi_d^2 dataset in the same format as LogNormal.nc

The parameters estimated include the sought turbulence (epsilon) quantity and the spectral slope (beta_1). 

Loading NetCDF files

Matlab

load_netcdf_data('LogNormal.nc') which can be found in this repository in https://github.com/SCOR-ATOMIX/ShearProbes_BenchmarkComparison

JuliaLang

using NCDatasets
NCDataset("SyntheticSpectra.nc")