Data from: Observations of the internal wave to turbulence cascade
Data files
Mar 19, 2026 version files 16.50 MB
-
DA_final.zip
16.49 MB
-
README.md
5.97 KB
Abstract
Shear spectral energy density (shear spectra) is measured across three decades in vertical wavenumber using a Wirewalker profiler equipped with a microstructure instrument and a pulse-coherent Doppler sonar. We identify the features of the canonical ocean vertical shear spectrum, including an internal wave band, an intermediate saturation band whose spectral level is independent of turbulent dissipation, and a three-dimensional turbulence band. The internal wave band and saturation band of the spectrum scale as $\Phi_{IW}\sim\varepsilon\frac{1}{2}N1f{- \frac{1}{2}}$ and $\Phi_{sat}\sim N2 k_z^{-1}$ respectively. These scalings hold despite the deployment location at the head of a La Jolla canyon, a deep canyon incising the shelf off La Jolla, California, where weakly non-linear wave-wave interaction is not the primary physical process driving the forward energy cascade. In La Jolla Canyon, high-amplitude, tidally driven internal waves generate significant strain, resulting in turbulent events that cover a majority of the water column. During these events, we observe shear spectra with energy above the saturation level, which we interpret as a sign of the wave-turbulence transition. Finestructure parameterizations developed to predict mixing from shear spectra in the open ocean thermocline continue to predict average mixing well. The success of finestructure parameterizations implies a rate of downscale energy transfer consistent with the rate predicted from weakly non- linear wave-wave interactions, and suggests that the theoretical framework of the canonical shear spectrum can be used to make useful predictions in shallow, high-energy environments.
Dataset DOI: 10.5061/dryad.8sf7m0d44
Description of the data and file structure
This set of data and scripts is provided alongside the manuscript "Observations of the internal wave to turbulence cascade" to fulfill data availability requirements. Data used to generate figures in that manuscript are included. All code written by Devon Northcott (dnorthco@ucsd.edu). All data is in NetCDF format; NetCDF can be read in Python using the netCDF4 or Xarray packages, or in MATLAB, with instructions on how to open NetCDF files in MATLAB using their high-level NetCDT functions. Data is from a Wirewalker wave-powered profiler deployment in roughly 100m of water at the head of La Jolla Canyon (32.8564 N,-117.2665 W).
Files and variables
File: DA_final.zip
Description:
Data
RawProfile2_20231005-1132.nc: Contains a single wirewalker profile worth of Nortek signature 1000 coherent (HR) mode data. This includes intensity, signal correlation, and radial velocity, as well as sensor orientation information. The data in this file can be used to calculate turbulent dissipation of kinetic energy, epsilon.
Dimensions:
- Time: time from Nortek Signature 1000 UTC.
- Range: Range from the sonar head.
Variables:
- time: Time reported in seconds after 01-01-1970
- pressure: Pressure in dbar recorded by Nortek Signature 1000.
- v_a: Alias velocity in meters per second. This is based on the setup of the Nortek.
- pitch: instrument pitch in degrees. The instrument is mounted at 25 degrees off vertical.
- roll: instrument roll in degrees. The instrument is mounted at 25 degrees off vertical.
- yaw: instrument yaw in degrees.
- ncells: Number of configured measurement cells for each sonar ping
- vel: radial velocity in meters per second.
- corr: signal correlation in percent.
- int: signal intensity in relative dB.
ShearSpectraData.nc: Contains data required to calculate full water column shear spectra, as well as data required to sort and bin these spectra by water column average N^2 of epsilon. Data is from three sensors: incoherent mode Nortek signature 1000 pings processed into mean velocity, coherent (HR) mode pings processed into water column average spectra and epsilon, and Epsilometer shear probe data processed into water column average spectra.
Dimensions:
- depth: Depth in meters
- time: average time for each wirewalker profile
- k_NortekHR: wavenumber in for spectra calculated from Nortek Signature 1000 HR mode data.
- k_Epsi: wavenumber for spectra calculated from Epsilometer shear probe data.
Variables:
- time: average time for each wirewalker profile. Time reported in seconds after 01-01-1970
- velE: eastward current in m/s.
- velN: northward current in m/s.
- depth: depth in meters.
- epsiNortek: water column average turbulent dissipation, calculated by taking the mean of the log of dissipation weighted with a Hamming window over the full water column.
- N2: water column average buoyancy frequency, calculated by taking the mean buoyancy frequency weighted with a Hamming window over the full water column.
- specNortek: average water column shear spectra from Nortek coherent (HR) mode velocity fluctuations. Calculated by taking velocity spectra, subtracting out a fitted noise level, and multiplying by (2pik)^2 to convert to wavenumber. Spectra are then log-transformed and averaged over the water column.
- specEpsi: average water column shear spectra from Epsilometer shear probe data. Calculated by taking the log of the spectra and averaging over the full water column.
- kNortek: wavenumber in cpm corresponding to Nortek-derived spectra.
- kEpsi: wavenumber in cpm corresponding to Epsilometer-derived spectra.
EpsilometerTurbulenceData.nc: Turbulent dissipation of kinetic energy estimated from Epsilometer shear probes. Contains the full record of processed Epsilometer data. This record is slightly shorter than other records from the same wirewalker deployment because of power constraints on the Epsilometer.
Dimensions:
- depth: Depth in meters
- time: average time for each wirewalker profile
Variables:
- time: average time for each wirewalker profile. Time reported in seconds after 01-01-1970
- depth: depth in meters.
- Epsilon: Turbulent dissipation of kinetic energy in W/kg.
- fom: Figure of merit. Used to determine how closely the measured spectra align with a theoretical spectrum. This is used as a data quality flag.
NortekTurbulenceData.nc: Turbulent dissipation of kinetic energy estimated from Nortek signature 1000 coherent (HR) mode velocity measurements. Contains the full record of processed Nortek signature coherent mode data.
Dimensions:
- depth: Depth in meters
- time: average time for each wirewalker profile
Variables:
- time: average time for each wirewalker profile. Time reported in seconds after 01-01-1970
- depth: depth in meters.
- Epsilon: Turbulent dissipation of kinetic energy in W/kg.
- N: Fitted noise floor in the spectra for a given bin. This quantity is in velocity spectrum units, normalized to cycles per radian.
- corr: bin average signal correlation for coherent mode velocity measurements, expressed as a percentage.
- N2: buoyancy frequency in s^-2.
Code/software
Software:
Matlab is required to run these scripts; they were created using version 2024b.
ProcessSingleProfile: Processes raw Nortek signature 1000 coherent (HR) mode to calculate depth-binned turbulent dissipation of kinetic energy (epsilon). Can be run on the data in the RawProfile2_20231005-1132.nc netCDF file.
PlotShearSpectraFromNCs.m: Uses data from ShearSpectraData.nc to generate shear spectra binned by water column average dissipation.
Utility: Various scripts that are required to run the two data processing scripts above.
