pith. sign in

arxiv: 2606.27597 · v1 · pith:JEYGTZKQnew · submitted 2026-06-25 · ⚛️ physics.flu-dyn · physics.ao-ph

Toward a Universal Framework for the Internal Gravity Wave Spectrum

Leticia Fabre-Lima (1) , Jeffrey Early (2) , Miles A. Sundermeyer (1) ((1) School for Marine Science , Technology University of Massachusetts Dartmouth , (2) NorthWest Research Associates) This is my paper

Pith reviewed 2026-06-29 00:27 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn physics.ao-ph
keywords internal gravity wavesGarrett-Munk spectrumvertical modeshorizontal wavenumbernon-hydrostatic dynamicsoceanic stratificationSturm-Liouville problem
0
0 comments X

The pith

A generalized framework formulates the internal gravity wave spectrum in horizontal wavenumber space using non-hydrostatic vertical modes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The classic Garrett-Munk spectrum for oceanic internal gravity waves is limited by its hydrostatic assumptions and idealized stratification. This paper establishes a generalized framework that uses numerically solved non-hydrostatic vertical modes formulated in horizontal wavenumber space. Energetic orthogonality of the modes dictates this choice of space over frequency space. The resulting model, paired with a generalized spectral function, gives expected energy distributions in depth, frequency, and wavenumber that account for boundary effects and non-hydrostatic dynamics.

Core claim

The formulation in horizontal wavenumber-vertical mode space with vertical modes from the fixed-K Sturm-Liouville problem, combined with a generalized spectral function, yields expected distributions of horizontal kinetic, vertical kinetic, and potential energy as functions of depth, frequency, and horizontal wavenumber.

What carries the argument

Non-hydrostatic vertical modes obtained from the fixed-K Sturm-Liouville problem in horizontal wavenumber space, which serve as the basis for the spectrum and distinguish hydrostatic and non-hydrostatic regimes via the deformation radius.

If this is right

  • Energy distributions show departures from GM theory associated with boundary effects and non-hydrostatic dynamics.
  • Vertical variance and high-frequency vertical kinetic energy receive improved representation.
  • Observed features of horizontal wavenumber spectra are reproduced.
  • Arbitrary stratification profiles and multiple turning depths are accommodated through numerical mode computation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The wavenumber-space approach enables straightforward extraction of spectra for direct comparison against observational or simulation data in wavenumber form.
  • Numerical solution of the Sturm-Liouville problem for any given profile suggests the framework can be applied to measured ocean density data without analytic simplification.
  • The deformation radius per mode supplies a practical criterion for switching between hydrostatic and non-hydrostatic approximations in layered models.

Load-bearing premise

Energetic orthogonality among wave modes requires that the formulation be cast in horizontal wavenumber space rather than frequency space.

What would settle it

Direct computation of energy orthogonality in a frequency-space version for non-hydrostatic modes, or mismatch between predicted and observed vertical kinetic energy distributions in a boundary-influenced non-hydrostatic regime, would falsify the central requirement.

read the original abstract

The Garrett-Munk (GM) spectrum has long provided a canonical model of the oceanic internal gravity wave field. However, it relies on hydrostatic assumptions and idealized stratification that limit its applicability where non-hydrostatic dynamics, vertical boundary effects, or non-monotonic stratification are important. Here we develop a generalized framework for the internal wave spectrum based on non-hydrostatic vertical modes formulated in horizontal wavenumber-vertical mode space. Energetic orthogonality among wave modes requires that such a formulation be cast in horizontal wavenumber space rather than frequency space. In this formulation, the deformation radius associated with each vertical mode provides a proxy for distinguishing hydrostatic and non-hydrostatic regimes. Vertical modes are obtained numerically from the fixed-K Sturm-Liouville problem, allowing arbitrary stratification and multiple turning depths. Combined with a generalized spectral function, the formulation yields expected distributions of horizontal kinetic, vertical kinetic, and potential energy as functions of depth, frequency, and horizontal wavenumber. Example applications illustrate departures from GM theory associated with boundary effects and non-hydrostatic dynamics, including improved representation of vertical variance and high-frequency vertical kinetic energy, while reproducing observed features of horizontal wavenumber spectra.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript develops a generalized framework for the oceanic internal gravity wave spectrum. It solves the fixed-K Sturm-Liouville problem for non-hydrostatic vertical modes at each horizontal wavenumber, then combines these modes with a generalized spectral function to produce expected distributions of horizontal kinetic, vertical kinetic, and potential energy as functions of depth, frequency, and horizontal wavenumber. The approach is positioned as an extension beyond the hydrostatic and idealized-stratification assumptions of the Garrett-Munk (GM) spectrum, with example applications illustrating improved representation of vertical variance, high-frequency vertical kinetic energy, and observed horizontal-wavenumber spectral features.

Significance. If the central derivations hold, the framework offers a more flexible, numerically realizable model capable of handling arbitrary stratification, multiple turning depths, and non-hydrostatic regimes. This could be valuable for ocean modeling and wave-energy studies where boundary effects or non-hydrostatic dynamics matter, while still recovering key observed spectral properties.

major comments (1)
  1. [Abstract] Abstract (paragraph beginning 'Energetic orthogonality among wave modes requires...'): the assertion that energetic orthogonality among wave modes necessitates formulation in horizontal-wavenumber space rather than frequency space is presented as a prerequisite but lacks an explicit derivation, orthogonality integral, or reference to the inner-product definition that would demonstrate non-orthogonality in frequency space. Because this choice determines the entire modal basis and spectral construction, a self-contained justification (or citation to a prior derivation) is required.
minor comments (2)
  1. The abstract refers to 'a generalized spectral function' whose explicit functional form, free parameters, and reduction to the GM limit are not stated; the main text should supply the precise definition (including any depth or wavenumber dependence) so that the claimed energy distributions can be reproduced.
  2. The example applications are described only qualitatively ('improved representation of vertical variance...'); quantitative metrics (e.g., RMS error against observations, comparison tables, or specific stratification profiles) should be added to allow assessment of the magnitude of the reported departures from GM theory.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and recommendation of minor revision. We address the single major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract (paragraph beginning 'Energetic orthogonality among wave modes requires...'): the assertion that energetic orthogonality among wave modes necessitates formulation in horizontal-wavenumber space rather than frequency space is presented as a prerequisite but lacks an explicit derivation, orthogonality integral, or reference to the inner-product definition that would demonstrate non-orthogonality in frequency space. Because this choice determines the entire modal basis and spectral construction, a self-contained justification (or citation to a prior derivation) is required.

    Authors: We agree that the abstract presents the requirement without a self-contained derivation. The orthogonality follows directly from the fixed-K Sturm-Liouville eigenvalue problem whose eigenfunctions satisfy the standard weighted inner product ∫ φ_m φ_n N^{2} dz = 0 for m eq n at fixed horizontal wavenumber K. In frequency space the vertical structure operator is no longer self-adjoint in the same way once the dispersion relation is imposed, so the same inner product does not yield orthogonality. We have added a concise statement of this inner-product definition and a one-sentence derivation to the abstract, together with a pointer to the full development in Section 2 of the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The derivation proceeds from the standard fixed-K Sturm-Liouville eigenvalue problem (a linear operator on arbitrary stratification) to obtain vertical modes, then superposes them with an externally supplied generalized spectral function to produce the energy distributions. The orthogonality argument is a physical justification for the choice of wavenumber basis rather than a self-referential definition of the spectrum itself. No parameter is fitted to a subset of the target data and then relabeled a prediction, no uniqueness theorem is imported from prior self-work, and no ansatz is smuggled via citation. The resulting spectra are therefore not equivalent to the inputs by construction; they remain falsifiable against independent observations and the GM benchmark.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based solely on the abstract; the central addition is the numerical vertical-mode solution and the choice of wavenumber space, with the generalized spectral function left unspecified.

axioms (1)
  • domain assumption Energetic orthogonality among wave modes requires formulation in horizontal wavenumber space rather than frequency space
    Explicitly stated as a prerequisite for the framework in the abstract.

pith-pipeline@v0.9.1-grok · 5759 in / 1229 out tokens · 22196 ms · 2026-06-29T00:27:04.363484+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

35 extracted references · 12 canonical work pages

  1. [1]

    G., 1975 a : Introduction to collection of papers on oceanic internal waves

    Briscoe, M. G., 1975 a : Introduction to collection of papers on oceanic internal waves. J.\ Geophys.\ Res., 80 (3), 289--290

    1975

  2. [2]

    G., 1975 b : Preliminary results from the trimoored internal wave experiment (iwex)

    Briscoe, M. G., 1975 b : Preliminary results from the trimoored internal wave experiment (iwex). J.\ Geophys.\ Res., 80 (27), 3872--3884

    1975

  3. [3]

    Chelton, D. B., R. A. DeSzoeke, M. G. Schlax, K. El Naggar, and N. Siwertz, 1998: Geographical variability of the first baroclinic rossby radius of deformation. J.\ Phys.\ Oceanogr., 28 (3), 433--460

    1998

  4. [4]

    Dashen, R., W. H. Munk, and K. M. Watson, 1979: Linear internal waves. Sound transmission through a fluctuating ocean , S. Flattee, Ed., Cambridge University Press, 44--54

    1979

  5. [6]

    a is \"a l \

    Desaubies, Y. J. F., 1975: A linear theory of internal wave spectra and coherences near the v \"a is \"a l \"a frequency. J.\ Geophys.\ Res., 80, 895–899, doi:10.1080/03091927308236112

    work page doi:10.1080/03091927308236112 1975

  6. [7]

    Early, J. J., M. P. Lelong, and K. S. Smith, 2020: Fast and accurate computation of vertical modes. J.\ Adv.\ Model.\ Earth Syst., 12 (2), e2019MS001\,939

    2020

  7. [8]

    C., 1985: Some characteristics of internal gravity waves in the equatorial pacific

    Eriksen, C. C., 1985: Some characteristics of internal gravity waves in the equatorial pacific. Journal of Geophysical Research: Oceans, 90 (C4), 7243--7255

    1985

  8. [9]

    Munk, 1972: Space-time scales of internal waves

    Garrett, C., and W. Munk, 1972: Space-time scales of internal waves. Geophys. Fluid Dyn., 3, 225–264, doi:10.1080/03091927208236082

    work page doi:10.1080/03091927208236082 1972

  9. [10]

    Munk, 1975: Space-time scales of internal waves: A progress report

    Garrett, C., and W. Munk, 1975: Space-time scales of internal waves: A progress report. J.\ Geophys.\ Res., 80 (3), 291--297

    1975

  10. [11]

    Munk, 1979: Internal waves in the ocean

    Garrett, C., and W. Munk, 1979: Internal waves in the ocean. Annu.\ Rev.\ Fluid Mech., 11, 339--369, doi:10.1146/annurev.fl.11.010179.002011

    work page doi:10.1146/annurev.fl.11.010179.002011 1979

  11. [12]

    C., 1989: Scaling turbulent dissipation in the thermocline

    Gregg, M. C., 1989: Scaling turbulent dissipation in the thermocline. J.\ Geophys.\ Res.: Oceans, 94, 9686–9698, doi:10.1029/JC094iC07p09686

    work page doi:10.1029/jc094ic07p09686 1989

  12. [13]

    induced diffusion

    Henyey, F. S., and N. Pomphrey, 1983: Eikonal description of internal wave interactions: A non-diffusive picture of “induced diffusion”. Dynamics of atmospheres and oceans, 7 (4), 189--219

    1983

  13. [14]

    Henyey, F. S., J. Wright, and S. M. Flatté, 1986: Energy and action flow through the internal wave field: An eikonal approach. J.\ Geophys.\ Res.: Oceans, 91, 8487–8495, doi:10.1029/JC091iC07p08487

    work page doi:10.1029/jc091ic07p08487 1986

  14. [15]

    Johnson, C. L., C. S. Cox, and B. Gallagher, 1978: The separation of wave-induced and intrusive oceanic finestructure. J.\ Phys.\ Oceanogr., 8 (5), 846--860

    1978

  15. [16]

    H., and G

    K \"a se, R. H., and G. Siedler, 1980: Internal wave kinematics in the upper tropical atlantic. Oceanography and Surface Layer Meteorology in the B/C Scale, Elsevier, 161--189

    1980

  16. [17]

    J., 1973: Profile of an isopycnal surface in the main thermocline of the sargasso sea

    Katz, E. J., 1973: Profile of an isopycnal surface in the main thermocline of the sargasso sea. Journal of Physical Oceanography, 3 (4), 448--457

    1973

  17. [18]

    J., and M

    Katz, E. J., and M. G. Briscoe, 1979: Vertical coherence of the internal wave field from towed sensors. J.\ Phys.\ Oceanogr., 9 (3), 518--530

    1979

  18. [19]

    M., 2016: The vertical mode decomposition of surface and internal tides in the presence of a free surface and arbitrary topography

    Kelly, S. M., 2016: The vertical mode decomposition of surface and internal tides in the presence of a free surface and arbitrary topography. J.\ Phys.\ Oceanogr., 46 (12), 3777--3788

    2016

  19. [20]

    Le Boyer, A., and M. H. Alford, 2021: Variability and sources of the internal wave continuum examined from global moored velocity records. J.\ Phys.\ Oceanogr., 51 (9), 2807--2823

    2021

  20. [21]

    D., 1983: Internal waves in the ocean: A review

    Levine, M. D., 1983: Internal waves in the ocean: A review. Rev. Geophys., 21 (5), 1206--1216

    1983

  21. [22]

    and Kelly, K

    Levine, M. D., 2002: A modification of the garrett–munk internal wave spectrum. J.\ Phys.\ Oceanogr., 32, 3166–3181, doi:10.1175/1520-0485(2002)032<3166:AMOTGM>2.0.CO;2

    work page doi:10.1175/1520-0485(2002)032 2002

  22. [23]

    Levine, M. D., J. D. Irish, T. E. Ewart, and S. A. Reynolds, 1986: Simultaneous spatial and temporal measurements of the internal wave field during mate. Journal of Geophysical Research: Oceans, 91 (C8), 9709--9719, doi:https://doi.org/10.1029/JC091iC08p09709

    work page doi:10.1029/jc091ic08p09709 1986

  23. [24]

    A., and Coauthors, 2019: World ocean atlas 2018, volume 5: Density

    Locarnini, R. A., and Coauthors, 2019: World ocean atlas 2018, volume 5: Density. Tech. rep., NOAA

    2019

  24. [25]

    H., and P

    McComas, C. H., and P. M \"u ller, 1981: The dynamic balance of internal waves. J.\ Phys.\ Oceanogr., 11, 970--986, doi:10.1175/1520-0485(1981)011<0970:TDBOIW>2.0.CO;2

    work page doi:10.1175/1520-0485(1981)011 1981

  25. [26]

    M\"uller, P., D. J. Olbers, and J. Willebrand, 1978: The iwex spectrum. Journal of Geophysical Research: Oceans, 83 (C1), 479--500, doi:https://doi.org/10.1029/JC083iC01p00479

    work page doi:10.1029/jc083ic01p00479 1978

  26. [27]

    Evolution of Physical Oceanography , B

    Munk, W., 1981: Internal waves and small-scale processes. Evolution of Physical Oceanography , B. A. Stommel, and C. Wunsch, Eds., The MIT Press, 265

    1981

  27. [28]

    Goldin, J

    Pinkel, R., M. Goldin, J. Smith, O. Sun, A. Aja, M. Bui, and T. Hughen, 2011: The wirewalker: A vertically profiling instrument carrier powered by ocean waves. J.\ Atmos.\ Oceanic Technol., 28 (3), 426--435

    2011

  28. [29]

    Nguyen, J

    Pinkel, R., S. Nguyen, J. A. Smith, A. J. Lucas, B. D. Reineman, and A. F. Waterhouse, 2023: Vertical momentum transport by internal gravity waves above the equatorial undercurrent at 140°w. Geophys.\ Res.\ Lett., 50, n/a -- n/a, doi:10.1029/2022GL101630

    work page doi:10.1029/2022gl101630 2023

  29. [30]

    L., and Y

    Polzin, K. L., and Y. V. Lvov, 2011: Toward regional characterizations of the oceanic internal wavefield. Rev. Geophys., 49 (4)

    2011

  30. [31]

    Polzin, K. L., A. C. Naveira Garabato, T. N. Huussen, B. M. Sloyan, and S. Waterman, 2014: Finescale parameterizations of turbulent dissipation. J.\ Geophys.\ Res., 119 (2), 1383–1419

    2014

  31. [32]

    Polzin, K. L., J. M. Toole, and R. W. Schmitt, 1995: Finescale parameterizations of turbulent dissipation. J.\ Phys.\ Oceanogr., 25, 306--328, doi:10.1029/JC080i003p00291

    work page doi:10.1029/jc080i003p00291 1995

  32. [33]

    Hibiya, and A

    Takahashi, A., T. Hibiya, and A. C. Naveira Garabato, 2021: Influence of the distortion of vertical wavenumber spectra on estimates of turbulent dissipation using the finescale parameterization: Eikonal calculations. J.\ Phys.\ Oceanogr., 51, 1723–1733, doi:10.1175/JPO-D-20-0196.1

    work page doi:10.1175/jpo-d-20-0196.1 2021

  33. [34]

    J.\ Phys.\ Oceanogr., 6 (4), 471--485

    Wunsch, C., 1976: Geographical variability of the internal wave field: A search for sources and sinks. J.\ Phys.\ Oceanogr., 6 (4), 471--485

    1976

  34. [35]

    Webb, 1979: The climatology of deep ocean internal waves

    Wunsch, C., and S. Webb, 1979: The climatology of deep ocean internal waves. J.\ Phys.\ Oceanogr., 9 (2), 235--243

    1979

  35. [36]

    Zheng, B., A. J. Lucas, R. Pinkel, and A. Le Boyer, 2022: Fine-scale velocity measurement on the wirewalker wave-powered profiler. J.\ Atmos.\ Oceanic Technol., 39 (2), 133--147

    2022