Inertio-Gravity Wave Band in Geophysical Flows

• The inertio-gravity wave band is defined as the frequency interval (f < |ω| < N) arising from the interplay between Coriolis and buoyancy forces in rotating, stratified fluids.
• It leverages linearized hydrodynamic equations to characterize energy cascades, modal properties, and spectral distributions in oceanic and atmospheric systems.
• Applications include turbulence modeling and boundary-layer analysis, with experimental and observational data validating key scale-invariance and isotropization phenomena.

The inertio-gravity wave (IGW) band denotes the characteristic frequency interval and associated spectral domain in rotating, stably stratified fluids where wave propagation results from interplay between the Coriolis and buoyancy restoring forces. This band arises naturally from the governing linearized hydrodynamic equations, spans a well-defined region in frequency–wavenumber space, and is of central importance to large-scale geophysical dynamics, mesoscale energy cascades, turbulence, and boundary-layer phenomena in the ocean, atmosphere, and planetary interiors.

1. Linear Theory and Definition of the Inertio-Gravity Band

In a rotating, stably stratified Boussinesq fluid, the governing linearized equations for small-amplitude perturbations and their associated pressure, velocity, and buoyancy fields yield a fundamental dispersion relation for plane waves of the form

$\omega^2 = \frac{N^2 k_h^2 + f^2 m^2}{k_h^2 + m^2}$

where $k_h = \sqrt{k_x^2 + k_y^2}$ (horizontal wavenumber), $m$ (vertical wavenumber), $N$ (Brunt–Väisälä frequency), and $f$ (vertical component of the Coriolis parameter; $f=2\Omega\sin\phi$ at latitude $\phi$) (Haren, 2023, Shavit et al., 4 Jan 2026, Heifetz et al., 2021). The admissible frequency band for propagating inertio-gravity waves is bounded from below by $f$ (inertial oscillations) and from above by $N$ (pure internal gravity waves), i.e.,

$f < |\omega| < N$

This band exists in both unbounded and confined domains, with specific boundary-induced modifications where relevant (Verdìère et al., 2024, Vidal et al., 2024).

In the simplified rotating shallow-water system, the inertio-gravity band is explicit in

$k_h = \sqrt{k_x^2 + k_y^2}$0

with $k_h = \sqrt{k_x^2 + k_y^2}$1 the non-rotating shallow-water gravity-wave speed, so that $k_h = \sqrt{k_x^2 + k_y^2}$2 (Heifetz et al., 2021). The general result is a spectral gap below $k_h = \sqrt{k_x^2 + k_y^2}$3 (no propagating waves), a super-inertial region $k_h = \sqrt{k_x^2 + k_y^2}$4, and an upper cutoff for internal waves at $k_h = \sqrt{k_x^2 + k_y^2}$5 in stratified settings.

2. Geometry in Wavenumber Space and Modal Properties

For the three-dimensional Boussinesq system,

$k_h = \sqrt{k_x^2 + k_y^2}$6

the IGW band at fixed frequency corresponds in $k_h = \sqrt{k_x^2 + k_y^2}$7-space to a double cone (the “constant-frequency cone”), comprising all wave vectors for which the above relation is satisfied given $k_h = \sqrt{k_x^2 + k_y^2}$8 (Savva et al., 2020, Kafiabad et al., 2019). In spherical coordinates $k_h = \sqrt{k_x^2 + k_y^2}$9,

$m$0

This geometric characterization is critical for understanding IGW energy transfer, scattering, and cascade phenomena: all nontrivial energy redistribution by nonlinear or stochastic processes acts along or within these surfaces, and the restriction to the band is enforced by the underlying dispersion relation and restoring forces.

3. Impact of Boundaries, Stratification, and Rotation

In bounded domains, surface and interior (bulk) spectrum are distinguished by the properties of the governing differential operators. For a smooth, compact domain $m$1, the linearized pressure satisfies a Poincaré-type elliptic–hyperbolic operator,

$m$2

with rigidity (impenetrable boundary) enforced as $m$3 (Verdìère et al., 2024).

Surface-trapped (Kelvin-like) waves arise when $m$4, with spectrum determined by a pseudo-differential boundary operator $m$5, emerging as a consequence of boundary-induced ellipticity. Bulk IGW modes exist only for $m$6, i.e., within the inertio-gravity band (Verdìère et al., 2024, Vidal et al., 2024). In ellipsoidal geometries, the interior spectrum in $m$7 becomes pure-point: eigenvalues accumulate densely and correspond to polynomially regular eigenfunctions, with spectral density shaped by the geometry and rotation axis tilt.

In layered or nonuniform stratification, transmission and reflection phenomena at interfaces (with $m$8 varying rapidly) cause partial amplitude loss for IGWs crossing from stratified to homogeneous regions, with observed energy transmission coefficients $m$9 (amplitude ratio $N$0) at sharp interfaces (Haren, 2023).

4. Energy Cascades, Wave Turbulence, and Spectral Distributions

The inertio-gravity wave band forms the backbone of kinetic energy transfer in geophysical flows, mediating the forward cascade from large to small scales via both wave–wave resonances and wave–flow scattering. Kinetic equations derived under random phase assumptions and weak turbulence closures yield steady-state turbulent spectra, notably:

• The scale-invariant solution $N$1 with $N$2, $N$3, which transforms directly into the near-universal Garrett–Munk spectrum observed in the ocean, $N$4 (Shavit et al., 4 Jan 2026).
• Spectra on the constant-frequency cone exhibit a $N$5 power-law in wavenumber, as predicted by both diffusion approximations and kinetic equations for IGW–turbulence interactions (Kafiabad et al., 2019, Savva et al., 2020). Laboratory and observational data find $N$6 at large (rotational) scales transitioning to $N$7 and then $N$8 (wave-dominated and turbulent IGW band) at smaller and intermediate scales (Rodda et al., 2020).

These results are robust to moderate variations in $N$9, $f$0, and stratification profile, provided nonlinearity is weak–moderate and the Rossby number remains small.

5. Scattering, Isotropization, and Forward Cascade Mechanisms

Scattering of IGWs by background geostrophic turbulence redistributes wave-action within the cone of constant frequency, generating both radial (cascade) and azimuthal (isotropization) transfer in wavenumber space (Savva et al., 2020, Kafiabad et al., 2019). The key mechanism is summarized by the kinetic equation

$f$1

with $f$2 determined by the flow energy spectrum and resonance conditions. Same- and cross-nappe scattering (vertical propagation direction) are both realized, resulting in horizontal isotropization and a forward energy cascade to high $f$3. The steady-state solution in the WKBJ limit is $f$4 (Savva et al., 2020).

Validation via direct numerical simulation of the three-dimensional Boussinesq equations confirms the dominance of this process in the regime of subcritical Rossby and Froude numbers and the restriction of transfer to the IGW band (Kafiabad et al., 2019, Savva et al., 2020).

6. Experimental and Observational Evidence

Laboratory rotating annulus experiments with stably stratified, differentially heated fluids manifest clear IGW bands: energy spectra transition from a $f$5 vortex-dominated regime at the largest scales to an $f$6 slope in the IGW-dominated mesoscale, where gravity waves carry the majority ($f$770–90%) of kinetic energy. A two-step Helmholtz and linear–wave decomposition confirms that the boundaries of the IGW band coincide with Doppler-shifted $f$8 and $f$9 (Rodda et al., 2020).

Deep-sea observations in stratified–homogeneous layer transitions (Western Mediterranean) find near-inertial IGW amplitudes reduced by a factor of $f=2\Omega\sin\phi$01.3 at rapid changes in stratification, and polarization spectra bounded by theoretical IGW frequency limits dictated by local values of $f=2\Omega\sin\phi$1, $f=2\Omega\sin\phi$2, and non-traditional Coriolis parameter $f=2\Omega\sin\phi$3 (Haren, 2023).

7. Relevance, Extensions, and Outstanding Issues

The inertio-gravity wave band is a fundamental organizing structure in geophysical fluid dynamics. It governs the spectrum of fast Poincaré waves, determines the partition between balanced (geostrophic) and unbalanced (wave) motions, underlies oceanic energy spectra (Garrett–Munk), and shapes how mesoscale and submesoscale turbulence interacts with large-scale forcing and boundary layers (Heifetz et al., 2021, Shavit et al., 4 Jan 2026, Verdìère et al., 2024, Vidal et al., 2024, Haren, 2023).

Key physical features of the band:

Regime	Frequency interval	Dominant restoring force
Inertial	$f=2\Omega\sin\phi$4	Coriolis (rotation)
SGW/IGW Band	$f=2\Omega\sin\phi$5	Mixed (rot.+buoyancy)
Gravity wave limit	$f=2\Omega\sin\phi$6	Buoyancy (stratification)

A plausible implication is that boundary-induced surface-wave branches and spectral discretization in confined domains further modulate the canonical IGW band, affecting energy localization and parametric instability thresholds (Verdìère et al., 2024, Vidal et al., 2024). The role of non-traditional Coriolis terms ($f=2\Omega\sin\phi$7) and fine-scale stratification remains an active area with direct importance for ocean mixing, wave reflection/refraction, and internal wave energetics (Haren, 2023).

The IGW band thus provides a rigorous, predictive framework for interpreting spectral observations, laboratory results, and high-resolution numerical simulation in rotating stratified fluids across a vast range of geophysical and astrophysical settings.