Data from: A method of separating linear internal wave and vortical mode energies using shipboard ADCP velocity measurements
Data files
Dec 17, 2025 version files 7.67 MB
-
README.md
17.79 KB
-
SupportingData.zip
7.66 MB
Abstract
We present a method to quantify total, horizontal kinetic and available potential energies of linear internal waves (IW) and vortical mode (VM) using only two-dimensional (2D) (depth, along-track distance) measurements of horizontal velocity, such as those commonly taken by oceanic shipboard ADCP (SADCP). Previous IW and VM energy decomposition methods (Bühler et al. 2014, 2017, and their extensions) require both velocity and buoyancy measurements. Applying Helmholtz decomposition, 2D horizontal kinetic energy wavenumber (kx , kz) spectra are projected onto divergent Kdiv and rotational Krot components. IW total energy spectrum is EIW = 2Kdiv . VM total energy is EVM=1/Bu[(1+Bu)Krot-Kdiv], where Bu is Burger number with N and f buoyancy and inertial frequencies, kh horizontal wavenumber magnitude and kz vertical wavenumber. IW and VM horizontal kinetic energy (K) and available potential energy (P) can be inferred from EVM and EIW as functions of Bu. The proposed method, derived directly from IW and VM theoretical polarization relations, is demonstrated using two sets of velocity and density data. EVM derived by this new method agrees with results computed using the Bühler et al. (2014, 2017) method at Bu~O(1) and within a factor of ~ 2–3 elsewhere, confirming that IW and VM energy can be separated using only velocity data. At Bu ≪ O (0.1), EVM is dominated by PVM, with KVM/PVM=Bu, and using K alone to extract EVM through the proposed method is challenging due to inherent uncertainty in spectral measurements. This method could be applied to global SADCP datasets to separate upper ocean IW and VM energy contributions in different dynamical regions at horizontal scales O (100m) – O (100km) and vertical scales O (10m) – O (100m).
Dataset DOI: 10.5061/dryad.70rxwdcc9
Description of the data and file structure
Data to accompany the article "A Method of Separating Linear Internal Wave and Vortical Mode Energies Using
Shipboard ADCP Velocity Measurements" (Journal of Atmospheric and Oceanic Technology, DOI: 10.1175/JTECH-D-25-0076.1) by Anda Vladoiu and Ren-Chieh Lien.
File: SupportingData.zip
Description: The .zip archive contains .mat processed data files used in creating the figures in the article. Please see article for detailed data description and processing methods.
Descriptions of each .mat file, corresponding to distinct figures, are below, with units in square brackets. Each file contains structures pertaining to variables plotted in each figure panel, named in alphabetical order.
Contact Anda Vladoiu (avladoiu@uw.edu) for additional information.
file: Fig1.mat
Fig. 1. Model (kh, kz) spectra for IW based on GM IW model spectra: (a) KIWdiv, (b) KIWrot, (c) PIW, and (d) EIW=KIWdiv +KIWrot +PIW, and for VM based on mean, smoothed spectrum from Vladoiu et al. (2024): (f) KVMrot , (f) PVM, and (g) EVM=KVMdiv+KVMrot+PVM, where KVMdiv=0. (h) Total inferred Kdiv=KIWdiv+KVMdiv , (i) total Krot=KIWrot+KVMrot , and (j) the total energy of VM derived from the K energy decomposition method E*= (1/Bu)[(1 + Bu)Krot - Kdiv]. (k) Kdiv^^~ and (l) Krot^^~ are realizations of Kdiv and Krot spectra following chi2 distribution with 20 DOF. (m) EVM~~~ is the VM total energy derived from the K energy decomposition method using Kdiv^^~ and Krot^^~.
Structure "Fig1a" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"KdivIW" = log10 internal wave divergent component of kinetic energy 2D spectrum [m^2/s^2/cpm^2]
Structure "Fig1b" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"KrotIW" = log10 internal wave rotational component of kinetic energy 2D spectrum [m^2/s^2/cpm^2]
Structure "Fig1c" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"PIW" = log10 internal wave potential energy 2D spectrum [m^2/s^2/cpm^2]
Structure "Fig1d" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"EIW" = log10 internal wave total energy 2D spectrum [m^2/s^2/cpm^2]
Structure "Fig1e" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"KrotVM" = log10 vortical mode rotational component of kinetic energy 2D spectrum [m^2/s^2/cpm^2]
Structure "Fig1f" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"PVM" = log10 vortical mode potential energy 2D spectrum [m^2/s^2/cpm^2]
Structure "Fig1g" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"EVM" = log10 vortical mode total energy 2D spectrum [m^2/s^2/cpm^2]
Structure "Fig1h" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"Kdiv" = log10 total divergent kinetic energy 2D spectrum [m^2/s^2/cpm^2]
Structure "Fig1i" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"Krot" = log10 total rotational kinetic energy 2D spectrum [m^2/s^2/cpm^2]
Structure "Fig1j" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"EVMstar" = log10 vortical mode total energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
Structure "Fig1k" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"Kdivtilde" = log10 total divergent kinetic energy 2D spectrum with chi-squared uncertainty [m^2/s^2/cpm^2]
Structure "Fig1l" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"Krottilde" = log10 total rotational kinetic energy 2D spectrum with chi-squared uncertainty [m^2/s^2/cpm^2]
Structure "Fig1j" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"EVMstartilde" = log10 vortical mode total energy 2D spectrum, inferred from the K energy decomposition method with chi-squared uncertainty [m^2/s^2/cpm^2]
file: Fig2.mat
Fig. 2. Dependence on Bu of spectra shown in Fig. 1 for (a) total IW and VM energy, (b) IW and VM horizontal kinetic energy, and (c) IW and VM available potential energy. Curves represent spectra averaged over Burger number bins.
Structure "Fig2a" contains variables:
"Bu" = Burger number
"EIW" = log10 internal wave total energy 2D spectrum [m^2/s^2/cpm^2]
"EVM" = log10 vortical mode total energy 2D spectrum [m^2/s^2/cpm^2]
"EVMstar" = log10 vortical mode total energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
"EVMstartilde" = log10 vortical mode total energy 2D spectrum, inferred from the K energy decomposition method with chi-squared uncertainty [m^2/s^2/cpm^2]
"BuPDF" = Burger number
"PDF" = PDF for all resolved (kh,kz)
"PDF_positiveEVMstartilde" = PDF for (kh,kz) where EVMstartilde > 0
Structure "Fig2b" contains variables:
"Bu" = Burger number
"KIW" = log10 internal wave kinetic energy 2D spectrum [m^2/s^2/cpm^2]
"KVM" = log10 vortical mode kinetic energy 2D spectrum [m^2/s^2/cpm^2]
"KVMstar" = log10 vortical mode kinetic energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
"KVMstartilde" = log10 vortical mode kinetic energy 2D spectrum, inferred from the K energy decomposition method with chi-squared uncertainty [m^2/s^2/cpm^2]
Structure "Fig2c" contains variables:
"Bu" = Burger number
"PIW" = log10 internal wave potential energy 2D spectrum [m^2/s^2/cpm^2]
"PVM" = log10 vortical mode potential energy 2D spectrum [m^2/s^2/cpm^2]
"PVMstar" = log10 vortical mode potential energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
"PVMstartilde" = log10 vortical mode potential energy 2D spectrum, inferred from the K energy decomposition method with chi-squared uncertainty [m^2/s^2/cpm^2]
file: Fig3.mat
Fig. 3. FIG. 3. Ratio between Burger number bin-averaged VM energy (Bu) from Kdiv^^~ and Krot^^~ with chi2 spectral uncertainty and VM energy (Bu) from Kdiv and Krot without the uncertainty for DOF between 4 and 120, assuming each 2D spectral estimate has 4 degrees of freedom.
Structure "Fig3" contains variables:
"Bu" = Burger number
"ratio" = ratio of 2D spectra EVMstartilde/EVMstar
"DOF" = degrees of freedom
file: Fig4.mat
Fig. 4. Location, date, instruments providing velocity and density measurements, depth range, resolved horizontal and vertical wavelengths, number of sections, mean N and f, and trajectories (colors from blue to yellow indicate the time elapsed, and black circles indicate the beginning of each section) for the (a) Sargasso Sea and (b) eastern North Pacific datasets.
Structure "Fig4a" contains variables:
"lon" = longitude [degW] for Sargasso Sea sections
"lat" = latitude [degN] for Sargasso Sea sections
Structure "Fig4b" contains variables:
"lon" = longitude [degW] for eastern North Pacific sections
"lat" = latitude [degN] for eastern North Pacific sections
"time" = Matlab time [days]
file: Fig5.mat
Fig. 5. Mean eastern North Pacific (kh, kz) spectrum for (a) EVM derived from the K-P energy decomposition using both velocity and density measurements and (b) EVM* derived from the K energy decomposition using velocity measurements only. (c) Ratio of the two VM energy estimates and (d) their scatterplot.
Structure "Fig5a" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"EVM" = log10 vortical mode total energy 2D spectrum [m^2/s^2/cpm^2]
Structure "Fig5b" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"EVMstar" = log10 vortical mode total energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
Structure "Fig5c" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"ratio" = log10 ratio of vortical mode total energy 2D spectrum inferred from the K-P energy decomposition method to the vortical mode total energy 2D spectrum inferred from the K energy decomposition method
Structure "Fig5d" contains variables:
"EVM" = vortical mode total energy 2D spectrum [m^2/s^2/cpm^2]
"EVMstar" = vortical mode total energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
"Bu" = log10 Burger number
file: Fig6.mat
Fig. 6. Dependence on Bu of mean eastern North Pacific (kh, kz) spectrum for (a) total IW and VM energy, (b) IW and VM horizontal kinetic energy, and (c) IW and VM available potential energy. Curves represent spectra averaged over Burger number bins.
Structure "Fig6a" contains variables:
"Bu" = Burger number
"EIW" = log10 internal wave total energy 2D spectrum [m^2/s^2/cpm^2]
"EVM" = log10 vortical mode total energy 2D spectrum [m^2/s^2/cpm^2]
"EVMstar" = log10 vortical mode total energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
"BuPDF" = Burger number
"PDF" = PDF for all resolved (kh,kz)
"PDF_positiveEVMstartilde" = PDF for (kh,kz) where EVMstartilde > 0
Structure "Fig6b" contains variables:
"Bu" = Burger number
"KIW" = log10 internal wave kinetic energy 2D spectrum [m^2/s^2/cpm^2]
"KVM" = log10 vortical mode kinetic energy 2D spectrum [m^2/s^2/cpm^2]
"KVMstar" = log10 vortical mode kinetic energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
Structure "Fig6c" contains variables:
"Bu" = Burger number
"PIW" = log10 internal wave potential energy 2D spectrum [m^2/s^2/cpm^2]
"PVM" = log10 vortical mode potential energy 2D spectrum [m^2/s^2/cpm^2]
"PVMstar" = log10 vortical mode potential energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
file: Fig7.mat
Fig. 7. Mean Sargasso Sea (kh, kz) spectrum for (a) EVM derived from the K-P energy decomposition using both velocity and density measurements and (b) EVM* derived from the K energy decomposition using velocity measurements only. (c) Ratio of the two VM energy estimates and (d) their scatterplot.
Structure "Fig7a" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"EVM" = log10 vortical mode total energy 2D spectrum [m^2/s^2/cpm^2]
Structure "Fig7b" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"EVMstar" = log10 vortical mode total energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
Structure "Fig7c" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"kz" = vertical wavenumber [cpm]
"ratio" = log10 ratio of vortical mode total energy 2D spectrum inferred from the K-P energy decomposition method to the vortical mode total energy 2D spectrum inferred from the K energy decomposition method
Structure "Fig7d" contains variables:
"EVM" = vortical mode total energy 2D spectrum [m^2/s^2/cpm^2]
"EVMstar" = vortical mode total energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
"Bu" = log10 Burger number
file: Fig8.mat
Fig. 8. Dependence on Bu of Sargasso Sea (kh, kz) spectrum for (a) total IW and VM energy, (b) IW and VM horizontal kinetic energy, and (c) IW and VM available potential energy. Curves represent spectra averaged over Burger number bins.
Structure "Fig8a" contains variables:
"Bu" = Burger number
"EIW" = log10 internal wave total energy 2D spectrum [m^2/s^2/cpm^2]
"EVM" = log10 vortical mode total energy 2D spectrum [m^2/s^2/cpm^2]
"EVMstar" = log10 vortical mode total energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
"BuPDF" = Burger number
"PDF" = PDF for all resolved (kh,kz)
"PDF_positiveEVMstartilde" = PDF for (kh,kz) where EVMstartilde > 0
Structure "Fig8b" contains variables:
"Bu" = Burger number
"KIW" = log10 internal wave kinetic energy 2D spectrum [m^2/s^2/cpm^2]
"KVM" = log10 vortical mode kinetic energy 2D spectrum [m^2/s^2/cpm^2]
"KVMstar" = log10 vortical mode kinetic energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
Structure "Fig8c" contains variables:
"Bu" = Burger number
"PIW" = log10 internal wave potential energy 2D spectrum [m^2/s^2/cpm^2]
"PVM" = log10 vortical mode potential energy 2D spectrum [m^2/s^2/cpm^2]
"PVMstar" = log10 vortical mode potential energy 2D spectrum, inferred from the K energy decomposition method [m^2/s^2/cpm^2]
file: Fig9.mat
Fig. 9. PDFs of IW to VM spectral ratios for (a) total energy, (b) horizontal kinetic energy, and (c) available potential energy. Spectra were computed from GO-SHIP P02W SADCP velocity data collected during May 2022, between 75- and 225-m depth, and averaged between 0.5<Bu<1.5, horizontal wavelengths 3.3–50 km, and vertical wavelengths 20–150 m.
Structure "Fig9a" contains variables:
"PDFx" = log10 internal wave to vortical mode total energy ratio
"PDFy" = PDF
Structure "Fig9b" contains variables:
"PDFx" = log10 internal wave to vortical mode kinetic energy ratio
"PDFy" = PDF
Structure "Fig9c" contains variables:
"PDFx" = log10 internal wave to vortical mode potential energy ratio
"PDFy" = PDF
file: Fig10.mat
Fig. 10. IW to VM spectral ratios for (a) total energy, (b) horizontal kinetic energy, and (c) available potential energy. Spectra were computed from GO-SHIP P02W SADCP velocity data collected during May 2022, between 75- and 225-m depth, and averaged between 0.5<Bu<1.5, horizontal wavelengths 3.3–50 km, and vertical wavelengths 20–150 m.
Structure "Fig10a" contains variables:
"lon" = longitude [degE] for GO-SHIP P02W sections
"lat" = latitude [degN] for GO-SHIP P02W sections
"ratio" = log10 internal wave to vortical mode total energy ratio
Structure "Fig10b" contains variables:
"lon" = longitude [degE] for GO-SHIP P02W sections
"lat" = latitude [degN] for GO-SHIP P02W sections
"ratio" = log10 internal wave to vortical mode kinetic energy ratio
Structure "Fig10c" contains variables:
"lon" = longitude [degE] for GO-SHIP P02W sections
"lat" = latitude [degN] for GO-SHIP P02W sections
"ratio" = log10 internal wave to vortical mode potential energy ratio
file: FigA1.mat
Fig. A1. IW phase speed C as a function of resolved along-track horizontal wavenumber kx, assuming horizontal isotropy, for the (a) eastern North Pacific, (b) Sargasso Sea , and (c) GO-SHIP P02W measurements, for different vertical wavelengths.
Structure "FigA1a" contains variables:
"kx" = along-track horizontal wavenumber [cpm]
"lambdaz16m" = internal wave phase speed for vertical wavelength lambdaz=16 m [m/s]
"lambdaz50m" = internal wave phase speed for vertical wavelength lambdaz=50 m [m/s]
"lambdaz100m" = internal wave phase speed for vertical wavelength lambdaz=100 m [m/s]
"lambdaz300m" = internal wave phase speed for vertical wavelength lambdaz=300 m [m/s]
"lambdaz500m" = internal wave phase speed for vertical wavelength lambdaz=500 m [m/s]
"lambdaz1000m" = internal wave phase speed for vertical wavelength lambdaz=1000 m [m/s]
Structure "FigA1b" contains variables:
"kx" = along-track horizontal wavenumber [cpm]
"lambdaz20m" = internal wave phase speed for vertical wavelength lambdaz=20 m [m/s]
"lambdaz40m" = internal wave phase speed for vertical wavelength lambdaz=40 m [m/s]
"lambdaz62m" = internal wave phase speed for vertical wavelength lambdaz=62 m [m/s]
"lambdaz300m" = internal wave phase speed for vertical wavelength lambdaz=300 m [m/s]
"lambdaz500m" = internal wave phase speed for vertical wavelength lambdaz=500 m [m/s]
"lambdaz1000m" = internal wave phase speed for vertical wavelength lambdaz=1000 m [m/s]
Structure "FigA1c" contains variables:
"kx" = along-track horizontal wavenumber [cpm]
"lambdaz20m" = internal wave phase speed for vertical wavelength lambdaz=20 m [m/s]
"lambdaz80m" = internal wave phase speed for vertical wavelength lambdaz=80 m [m/s]
"lambdaz150m" = internal wave phase speed for vertical wavelength lambdaz=150 m [m/s]
"lambdaz300m" = internal wave phase speed for vertical wavelength lambdaz=300 m [m/s]
"lambdaz500m" = internal wave phase speed for vertical wavelength lambdaz=500 m [m/s]
"lambdaz1000m" = internal wave phase speed for vertical wavelength lambdaz=1000 m [m/s]
file: FigB1.mat
Fig. B1. (a),(b) Demonstration of 1D kx horizontal wavenumber spectrum to 1D kh horizontal wavenumber magnitude spectrum conversion assuming horizontal isotropy.
Structure "FigB1a" contains variables:
"kx" = along-track horizontal wavenumber [cpm]
"Phi" = 1D one-sided wavenumber kx spectrum [cpm^-1]
Structure "FigB1b" contains variables:
"kh" = horizontal wavenumber magnitude [cpm]
"Phikh" = horizontal wavenumber magnitude wavenumber spectrum [cpm^-1]
"Phikxkh" = kx spectrum converted to kh spectrum [cpm^-1]
"ratio" = ratio of Phikxkh to Phikh
