🌀 nonlinear sciences

Instability of exact finite-amplitude acoustic-gravity waves

This paper demonstrates that exact finite-amplitude acoustic-gravity waves constructed by Godin are linearly unstable above a critical steepness, a threshold that depends on stratification and can vanish due to resonance in isothermal atmospheres, thereby limiting the validity of such solutions to specific shallow layers or wedges.

Original authors: Christian Puntini

Published 2026-10-27

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Christian Puntini

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The atmosphere and the ocean are not static backdrops; they are fluid layers constantly in motion, carrying energy across vast distances in the form of waves. Some of these waves are gentle ripples, while others are massive, invisible pulses that travel through the air and water, carrying heat and momentum. Scientists have long studied how these waves behave, particularly when they grow large enough to become "finite-amplitude" waves, where the water or air itself moves in complex, looping paths rather than just bobbing up and down. A key question in fluid dynamics is whether these large, exact wave patterns can sustain themselves or if they are destined to break down into turbulence. Understanding this is crucial because the stability of such waves determines how energy is distributed in our weather systems and oceans, and whether these massive motions can persist long enough to be observed or to influence the climate.

In a recent study, researchers investigated a specific family of these large waves, known as acoustic-gravity waves, which were recently described by a scientist named Godin. These are unique because, even though the fluid they travel through is compressible (meaning it can be squeezed, like air), the motion of the fluid parcels within the wave follows a very specific, rigid pattern where the density and pressure of each tiny packet of fluid remain constant as it moves. This creates a wave that looks and moves like a classic water wave, but exists in a compressible medium like the atmosphere. The central question was whether these exact, mathematically perfect waves are stable, or if they are inherently prone to falling apart. The researchers found that these waves are indeed unstable, but only if they become steep enough. Specifically, they discovered a critical threshold: if the wave's height relative to its width exceeds a certain limit, the wave will begin to grow exponentially in its internal disturbances, leading to a breakdown of the smooth motion.

The study revealed that this breaking point is not a single fixed number but depends on the environment. In a perfectly uniform fluid, the wave becomes unstable once its steepness reaches a value of one-third. This confirms a previous finding for simple water waves but extends it to the more complex, compressible waves found in the atmosphere. However, the researchers found that the presence of stable layers in the atmosphere, where temperature changes with height, can actually lower this threshold. In some cases, a very slight amount of stability in the air can make the wave unstable at a steepness slightly less than one-third. This means that in the real world, these waves might break down more easily than previously thought if the atmospheric conditions are just right.

Perhaps the most striking discovery involves a phenomenon called resonance. The researchers identified that for certain wavelengths, the wave becomes unstable even at very small steepness, essentially vanishing the safety margin that usually protects it. This happens when the natural "bouncing" frequency of the air parcels matches a specific multiple of the wave's own speed. For the fast-moving waves observed in the upper atmosphere, known as airglow waves, this resonance condition occurs at a wavelength of approximately 5.5 times the atmospheric scale height. When the researchers applied their findings to a specific wave observed by Hecht and colleagues, they found that the atmospheric layer where the wave was detected sat right on this resonance. In this specific layer, the critical steepness required for instability drops to about 0.1, which is much lower than the standard one-third. However, because the published data does not determine the actual steepness of the wave at that altitude, it remains unknown whether the wave was steep enough to be unstable or how rapidly it might have been growing.

The study also clarified where this instability lives within the wave. For waves traveling on the surface of the ocean, the instability is confined to a layer just beneath the free surface. For waves traveling along a sloping beach, the instability is concentrated in a wedge shape right at the shoreline. The researchers showed that uniform currents, or steady winds blowing in the same direction as the wave, do not change the stability of the wave; they only shift the speed at which the wave is observed. However, if the wind or current changes speed across the wave (a sheared current), the mathematical tools used in this study no longer apply, and the behavior of the wave remains an open question.

By combining advanced mathematical techniques with the specific properties of Godin's waves, the researchers provided a clear map of when and where these waves will break down. They demonstrated that while these waves can exist as exact solutions in theory, the real world—with its varying temperatures and specific wavelengths—often pushes them into a state of instability. The findings offer a new lens through which to view the turbulent upper atmosphere, suggesting that the dramatic, fast-moving waves seen in airglow data might be the result of these very specific instability mechanisms, growing and breaking in a matter of minutes. The work does not claim to solve every mystery of fluid motion, but it firmly establishes the conditions under which these elegant, large-scale waves lose their stability and begin to churn.

Technical Summary: Instability of Exact Finite-Amplitude Acoustic-Gravity Waves

Problem Statement
This paper investigates the linear stability of exact finite-amplitude wave motions in compressible fluids under a uniform gravity field, specifically the class of solutions constructed by Godin (2015). These motions are kinematically identical to Gerstner's trochoidal waves but exist in compressible media where pressure, density, and entropy are materially conserved ( $DP/Dt = D\rho/Dt = DS/Dt = 0$ ). While Godin demonstrated the existence of these exact solutions for surface waves, edge waves, and waves on currents, their stability properties remained unanalyzed. The central question is whether these finite-amplitude exact solutions are linearly stable or unstable, and how compressibility and stratification modify the stability thresholds known for incompressible Gerstner waves.

Methodology
The author employs the short-wavelength instability method (also known as the WKB or geometric optics method) for compressible flows, originally developed by Lifschitz & Hameiri (1991) and Lebovitz & Lifschitz (1992). This approach analyzes the evolution of localized, high-frequency perturbations advected by the basic flow.

Linearization and WKB Ansatz: Small perturbations to velocity, density, and pressure are superimposed on Godin's exact solution. The perturbations are assumed to have a rapidly oscillating phase $\Phi$ and slowly varying amplitudes.
Reduction to ODEs: By considering the limit of vanishing wavelength, the partial differential equations governing the perturbations reduce to a system of ordinary differential equations (ODEs) along the particle trajectories of the basic flow. The wavevector evolves according to the basic flow's velocity gradient, while the amplitudes evolve according to a linear system driven by the strain and vorticity of the wave.
Floquet Analysis: Since the basic flow is periodic in time (along particle paths), the resulting amplitude system is a linear ODE with periodic coefficients. The stability is determined by the Floquet multipliers of the monodromy matrix. The system is shown to be mathematically equivalent to the amplitude system analyzed by Puntini (2027) for mountain waves, allowing the direct application of previous stability criteria.
Parameter Mapping: The stability of Godin's waves is mapped to a two-parameter space defined by the wave steepness $\mathfrak{e} = e^{ka_2}$ and a stratification parameter $\nu^2 = (1-\mathfrak{e}^2)\nu_0^2$ , where $\nu_0^2 = \cos^2\alpha N_0^2 / \Omega^2$ . Here, $N_0$ is the buoyancy frequency of the hydrostatic reference state, $\alpha$ is the inclination of the wave plane, and $\Omega$ is the intrinsic orbital frequency.

Key Contributions and Results

General Stability Criterion: The stability of any fluid level in Godin's wave field is determined by the Hill discriminant $\Delta(\mathfrak{e}, \nu^2)$ derived in Puntini (2027). A level is exponentially unstable if and only if $|\Delta| > 2$ .
Isentropic Fluids (Recovery of Leblanc's Threshold): For isentropic fluids (including homogeneous incompressible liquids and adiabatic atmospheres) where the reference buoyancy frequency $N_0 = 0$ , the stratification parameter vanishes ( $\nu^2 = 0$ ). The analysis recovers the critical steepness threshold found by Leblanc (2004) for Gerstner waves: the waves are unstable if and only if $\mathfrak{e} > 1/3$ . This result extends Leblanc's finding to compressible surface and edge waves.
Effect of Weak Stratification: For weakly stable stratification, the critical steepness threshold is lowered. The paper derives an explicit expansion showing the threshold becomes $\mathfrak{e}_c \approx 1/3 - \frac{32}{81}\nu_0^2$ . Conversely, unstable stratification ( $N_0^2 < 0$ ) lowers the threshold further, allowing instability at arbitrarily small steepnesses.
Resonance in Isothermal Atmospheres: In an isothermal atmosphere, the stratification parameter $\nu_0$ $ν_{0}$ is constant across the wave field. The stability diagram reveals "resonance tongues" issuing from values where $\nu_0 = n/2$ $ν_{0} = n /2$ ( $n \in \mathbb{N}$ $n \in N$ ).
- The principal subharmonic resonance occurs at $\nu_0 = 1/2$ , corresponding to horizontal wavelengths $\lambda \approx 5.5 H / \cos\alpha$ (where $H$ is the pressure scale height).
- At this resonance, the critical steepness vanishes, meaning the wave is unstable even for infinitesimal amplitudes.
- Harmonic resonances ( $\nu_0 = 1, 2, \dots$ ) also create regions of instability.
Role of Currents: Uniform currents (both longitudinal and transverse) do not affect the stability criterion; they only Doppler-shift the observed frequency. However, sheared transverse currents destroy the invariant wavevector direction required for the current analysis, leaving that case open.
Spatial Localization of Instability:
- For surface waves in isentropic fluids, instability is confined to a layer beneath the free surface where the steepness exceeds $1/3$ .
- For edge waves, the unstable region forms a wedge bounded by the shoreline, the free surface, and a critical steepness contour.

Physical Estimates and Application
The paper applies these theoretical results to the "fast airglow wave" observed by Hecht et al. (2002), which Godin (2015) had previously interpreted using these exact solutions.

Using lidar temperature and wind profiles from Hecht et al., the author identifies an altitude range (approx. 99–104 km) where the local intrinsic frequency and wind conditions are consistent with Godin's solution.
In this specific layer, the stratification parameter places the wave near the principal subharmonic resonance ( $\nu_0 \approx 0.5$ ).
Consequently, the critical steepness in this layer is significantly reduced (median $\approx 0.09$ ) compared to the isentropic value of $1/3$ .
While the published data do not definitively determine the wave steepness at this altitude, the author notes that if the steepness were $\mathfrak{e} \approx 0.2$ , the layer would be unstable with e-folding times of 6–9 minutes, which is shorter than the observation duration.

Significance
The paper establishes that Godin's exact finite-amplitude solutions, while mathematically elegant, are generally linearly unstable above a critical steepness. It successfully bridges the gap between incompressible Gerstner wave theory and compressible acoustic-gravity wave theory. The primary significance lies in:

Extending the known instability threshold of $1/3$ from incompressible liquids to compressible fluids.
Demonstrating that weak stratification can lower this threshold, making large-amplitude waves more prone to instability.
Identifying that in the upper atmosphere (isothermal regime), resonance effects can render waves unstable at very small amplitudes, potentially explaining the rapid breakdown or dissipation of observed acoustic-gravity waves.

The work relies entirely on the mathematical machinery of short-wavelength analysis and Floquet theory, providing rigorous sufficient conditions for instability without proposing new experimental setups or future applications beyond the analysis of existing observational data.

More like this