Three-Wave Nonlinear Tides: Resonant Energy Transfer

• Three-wave nonlinear tides are resonant interactions where a parent tidal mode exchanges energy with two daughter modes via triadic coupling.
• They explain phenomena such as parametric subharmonic instability in ocean internal tides and analogous mode coupling in stars and compact binary systems.
• Diagnostic tools like bispectrum analysis and mode decomposition reveal the phase-coherent energy cascades and influence of dissipation on these tidal processes.

Three-wave nonlinear tides are tidal interactions in which a parent tide or tidally forced oscillation exchanges energy with two daughter waves or modes through a triad. In oceanography, this framework underlies parametric subharmonic instability of internal tides and the subsequent transfer of energy toward smaller, more dissipative scales. In stellar and compact-object astrophysics, the analogous language is three-mode coupling: a tidally forced parent mode destabilizes daughter modes, modifies dissipation, and alters orbital and gravitational-wave evolution. Across these settings, the common structure is resonant or near-resonant matching among three components, nonlinear phase coherence, and energy transfer that can be either sharply resonant or broadened by detuning, dissipation, or multibranch propagation physics (Frajka-Williams et al., 2014, Weinberg et al., 2011, Pitre et al., 10 Jun 2025).

1. Resonant structure of three-wave tidal coupling

The canonical triad conditions are conservation of frequency and wavenumber. In hydrodynamic notation these are often written as

$$k_1 = k_2 + k_3,\qquad \omega_1 = \omega_2 + \omega_3,$$

while in internal-wave and turbulence formalisms one equivalently writes

$$\mathbf{k} + \mathbf{k'} + \mathbf{k''} = 0,\qquad \omega + \omega' + \omega'' = 0.$$

Each participating component must also satisfy its own dispersion relation. Direct laboratory measurements on gravity-capillary surface waves verified the frequency and wavenumber conditions simultaneously and detected the daughter wave predicted by the triad geometry (Haudin et al., 2016).

When one member of the triad is treated as a pump, the daughter amplitudes obey damped amplitude equations. In one experimentally resolved form,

$$\dot{A}_2 = i I_{13} A_1 A_3^* - \alpha_2 A_2,\qquad \dot{A}_3 = i I_{12} A_1 A_2^* - \alpha_3 A_3,$$

with exponential growth when the forcing exceeds damping. The observed daughter amplitudes satisfy $A_{2,3}\propto e^{t/\tau}$ during the instability stage, and the interaction phase

$\Phi = \varphi_1 - \varphi_2 - \varphi_3$

locks in the stationary regime; in a toroidal hydrodynamic experiment it was measured close to $-\pi/2$, in quantitative agreement with resonant three-wave theory (Novkoski et al., 2023).

A common misconception is that all three-wave tidal interactions are parametric subharmonic instability. PSI is only a special case in which the parent wave transfers energy to two daughters with frequencies near half the parent frequency. More general triadic resonant instability need not produce exact subharmonics; the sloshing–gravity-capillary experiment explicitly distinguished its observed instability from parametric subharmonic instability (Frajka-Williams et al., 2014, Novkoski et al., 2023).

2. Internal tides and parametric subharmonic instability

In the oceanic internal-wave problem, three-wave nonlinear tides are classically represented by PSI of a coherent semidiurnal internal tide. A parent internal tide at frequency $\omega$ transfers energy to two daughter waves with $\omega' \approx \omega'' \approx \omega/2$, subject to

$$\mathbf{k} + \mathbf{k'} + \mathbf{k''} = 0,\qquad \omega + \omega' + \omega'' = 0.$$

The participating waves satisfy the internal-wave dispersion relation

$$\frac{k_1^2 + k_2^2}{k_3^2} = \frac{\omega^2 - f^2}{N^2 - \omega^2}.$$

For an $$\mathbf{k} + \mathbf{k'} + \mathbf{k''} = 0,\qquad \omega + \omega' + \omega'' = 0.$$0 tide, PSI is restricted to latitudes where $$\mathbf{k} + \mathbf{k'} + \mathbf{k''} = 0,\qquad \omega + \omega' + \omega'' = 0.$$1, that is, less than about $$\mathbf{k} + \mathbf{k'} + \mathbf{k''} = 0,\qquad \omega + \omega' + \omega'' = 0.$$2 (Frajka-Williams et al., 2014).

Assuming interaction in a vertical plane, the daughter-wave locus can be derived analytically. For exactly subharmonic daughters, the allowed vertical-wavenumber triads lie on an ellipse in nondimensional wavenumber space,

$$\mathbf{k} + \mathbf{k'} + \mathbf{k''} = 0,\qquad \omega + \omega' + \omega'' = 0.$$3

This geometric construction is not merely kinematic: peaks in the bispectrum aligned with this ellipse are diagnostic of PSI activity (Frajka-Williams et al., 2014).

The bispectrum supplies both identification and quantification. With

$$\mathbf{k} + \mathbf{k'} + \mathbf{k''} = 0,\qquad \omega + \omega' + \omega'' = 0.$$4

and normalized bicoherence,

$$\mathbf{k} + \mathbf{k'} + \mathbf{k''} = 0,\qquad \omega + \omega' + \omega'' = 0.$$5

one can isolate phase-locked triads and compute transfer rates from the nonlinear terms in the energy equations. In a fully nonlinear, non-hydrostatic Boussinesq model, bispectral transfer rates for PSI compared well to model growth rates of the daughter waves. By contrast, bispectra computed from HOME velocity profiles were relatively noisy and the signal was inconclusive, illustrating the observational difficulty of diagnosing triads in situ (Frajka-Williams et al., 2014).

3. From internal tides to internal wave turbulence

The broader oceanographic significance of three-wave nonlinear tides is the pathway from large-scale tidal injection to interior mixing. An estimated $$\mathbf{k} + \mathbf{k'} + \mathbf{k''} = 0,\qquad \omega + \omega' + \omega'' = 0.$$6 of power is required to support interior mixing, and roughly half is believed to come from tidal flow over topography producing internal gravity waves. A large-scale laboratory realization of this pathway used a $$\mathbf{k} + \mathbf{k'} + \mathbf{k''} = 0,\qquad \omega + \omega' + \omega'' = 0.$$7 stratified wavetank, an oscillating idealized ridge, and Reynolds numbers up to $$\mathbf{k} + \mathbf{k'} + \mathbf{k''} = 0,\qquad \omega + \omega' + \omega'' = 0.$$8 to mimic tidally forced internal-wave generation (Taebel et al., 2024).

In that system, Background Oriented Schlieren over the full tank showed the formation of various sets of subharmonics driven by Triadic Resonant Instabilities. At later times, the subharmonics engaged in further interactions and developed a continuum of waves at frequencies up to $$\mathbf{k} + \mathbf{k'} + \mathbf{k''} = 0,\qquad \omega + \omega' + \omega'' = 0.$$9. Fourier decomposition validated the three-wave resonant conditions, and the cascade displayed a backward transfer in frequency together with a forward transfer in vertical wavenumber (Taebel et al., 2024).

The same experiment also sharpened the classification of active triads. Elastic scattering was identified as a relevant nonlocal interaction in the fully evolved state, but the majority of the triads were local and had been historically overlooked. The reported partition was that local triads were the most common, at about $$\dot{A}_2 = i I_{13} A_1 A_3^* - \alpha_2 A_2,\qquad \dot{A}_3 = i I_{12} A_1 A_2^* - \alpha_3 A_3,$$0, while about $$\dot{A}_2 = i I_{13} A_1 A_3^* - \alpha_2 A_2,\qquad \dot{A}_3 = i I_{12} A_1 A_2^* - \alpha_3 A_3,$$1 were nonlocal, especially elastic scattering. Induced diffusion was largely absent because the experiment had limited spatial scale separation. This empirical balance complicates older pictures that emphasized scale-separated nonlocal pathways as the dominant route from tides to dissipation (Taebel et al., 2024).

4. Detuning, dissipation, and multibranch triads

Exact linear resonance is not the only route to strong three-wave transfer. In deep-water gravity waves, a distinct mechanism called precession resonance permits enhanced energy exchange even when the linear detuning is nonzero, $$\dot{A}_2 = i I_{13} A_1 A_3^* - \alpha_2 A_2,\qquad \dot{A}_3 = i I_{12} A_1 A_2^* - \alpha_3 A_3,$$2. The key condition is synchronization between a triad phase-precession frequency and a nonlinear frequency,

$$\dot{A}_2 = i I_{13} A_1 A_3^* - \alpha_2 A_2,\qquad \dot{A}_3 = i I_{12} A_1 A_2^* - \alpha_3 A_3,$$3

with a low-amplitude gravity-wave condition

$$\dot{A}_2 = i I_{13} A_1 A_3^* - \alpha_2 A_2,\qquad \dot{A}_3 = i I_{12} A_1 A_2^* - \alpha_3 A_3,$$4

Model calculations reported transfer efficiencies of up to $$\dot{A}_2 = i I_{13} A_1 A_3^* - \alpha_2 A_2,\qquad \dot{A}_3 = i I_{12} A_1 A_2^* - \alpha_3 A_3,$$5. This suggests that tidal systems may support dynamically important triads even when strict linear resonance fails, provided nonlinear phase dynamics produce the appropriate synchronization (Lucas et al., 2016).

Dissipation broadens triad selectivity in a different way. In capillary-gravity surface waves, experiments found daughter waves that verified the frequency and wavevector sums but not the linear dispersion relation. The mechanism was identified as forced three-wave interaction: the quadratic product of the two mother waves acts as an external forcing at $$\dot{A}_2 = i I_{13} A_1 A_3^* - \alpha_2 A_2,\qquad \dot{A}_3 = i I_{12} A_1 A_2^* - \alpha_3 A_3,$$6, and significant viscous dissipation broadens the bandwidth of the linear transfer function of the free surface. The result is a non-negligible off-dispersion response that remains phase coherent and measurable (Cazaubiel et al., 2019).

A further extension is multibranch coupling. In a torus of fluid, nonlinear three-wave resonant interactions were observed between the gravity-capillary and sloshing branches of the dispersion relation. The mother wave was a sloshing mode and the daughter waves were gravity-capillary modes. A triadic resonance instability was observed, with exponential growth, phase locking, and a cascade of additional waves for stronger forcing. The interaction efficiency was maximal when the group velocity of the sloshing mode matched the phase velocity of the gravity-capillary mode. The authors stated that this two-branch three-wave resonance mechanism is probably not restricted to hydrodynamics. A plausible implication is that tidal systems with several propagation branches or mode families may admit efficient branch-crossing energy transfer even when single-branch resonances are sparse (Novkoski et al., 2023).

5. Three-wave nonlinear tides in stars and compact binaries

In close binary systems, three-wave nonlinear tides are formulated as mode-amplitude dynamics. A representative amplitude equation is

$$\dot{A}_2 = i I_{13} A_1 A_3^* - \alpha_2 A_2,\qquad \dot{A}_3 = i I_{12} A_1 A_2^* - \alpha_3 A_3,$$7

where $$\dot{A}_2 = i I_{13} A_1 A_3^* - \alpha_2 A_2,\qquad \dot{A}_3 = i I_{12} A_1 A_2^* - \alpha_3 A_3,$$8 is the linear tidal forcing, $$\dot{A}_2 = i I_{13} A_1 A_3^* - \alpha_2 A_2,\qquad \dot{A}_3 = i I_{12} A_1 A_2^* - \alpha_3 A_3,$$9 a nonlinear tidal forcing term, and $A_{2,3}\propto e^{t/\tau}$0 the internal three-mode coupling coefficient. For solar-type stars, the linear tidal solution often used in binary evolution studies was found to be unstable over much of the parameter space in which it is employed. Resonantly excited gravity waves are unstable to parametric resonance for companion masses $A_{2,3}\propto e^{t/\tau}$1 at orbital periods $A_{2,3}\propto e^{t/\tau}$2 days, whereas the nearly static equilibrium tide is parametrically stable except for solar binaries with $A_{2,3}\propto e^{t/\tau}$3 days. The same study also found a collective instability in which a single parent drives $A_{2,3}\propto e^{t/\tau}$4 daughter waves as a coherent unit, with growth rates $A_{2,3}\propto e^{t/\tau}$5 times faster than the standard three-wave instability (Weinberg et al., 2011).

For white dwarf binaries, nonlinear dynamical tides are dominated over a broad period range by a global three-mode parametric instability rather than local wave breaking. The threshold energy for this instability is much lower than the local wave-breaking condition, and networks of coupled modes redistribute energy efficiently. A phenomenological consequence is that spin and orbit go in and out of synchronization, producing brief but significant dips in the tidal heating rate; a few percent of systems were predicted to be about ten times dimmer because they reside in such dips. The same framework implies that LISA and TianGO can constrain a white dwarf’s moment of inertia to better than $A_{2,3}\propto e^{t/\tau}$6 for deci-Hz systems (Yu et al., 2020).

For neutron stars in relativistic inspiral, the tidal deformation can be expanded simultaneously in time derivatives and nonlinearity,

$A_{2,3}\propto e^{t/\tau}$7

with the frequency-domain response

$A_{2,3}\propto e^{t/\tau}$8

The nonlinear constant $A_{2,3}\propto e^{t/\tau}$9 lowers the frequency parameter $\Phi = \varphi_1 - \varphi_2 - \varphi_3$0 by as much as $\Phi = \varphi_1 - \varphi_2 - \varphi_3$1 relative to a purely linear estimate. The paper explicitly linked this relativistic result to the three-wave mode-coupling picture emphasized by Yu et al. in Newtonian theory (Pitre et al., 10 Jun 2025).

Observationally, precise transit timing has started to constrain nonlinear tidal dissipation in planet-star systems. For six hot Jupiters around $\Phi = \varphi_1 - \varphi_2 - \varphi_3$2 GK hosts, transit times followed refined linear ephemerides and no orbital decay was detected. In HATS-18, WASP-19, and WASP-43, the analysis rejected a scenario with total dissipation of internal gravity waves, and for WASP-19 and WASP-43 the lower bounds on $\Phi = \varphi_1 - \varphi_2 - \varphi_3$3 exceeded the weakly nonlinear predictions by about an order of magnitude. These results indicate that the most efficient nonlinear dynamical-tide scenarios are not operating in those systems at current epochs (Maciejewski et al., 2024).

6. Diagnostics, interpretation, and recurrent misunderstandings

The observational and experimental diagnostics of three-wave nonlinear tides are fundamentally phase sensitive. Bispectrum and bicoherence are the standard third-order tools for frequency or wavenumber triads; in gravity-wave turbulence, significant bicoherence greater than $\Phi = \varphi_1 - \varphi_2 - \varphi_3$4, and sometimes up to $\Phi = \varphi_1 - \varphi_2 - \varphi_3$5 or higher, was used to identify robust three-wave phase coupling, primarily involving quasi-resonances of waves with second or higher-order harmonics. In the internal-tide context, transfer bispectra linked directly to the nonlinear energy budget and showed which terms were responsible for the exchange. In high-Reynolds-number internal-wave-turbulence experiments, Bispectral Mode Decomposition and Spectral Proper Orthogonal Decomposition extended this logic to fully space-time-resolved triad identification (Aubourg et al., 2017, Frajka-Williams et al., 2014, Taebel et al., 2024).

Several misunderstandings recur across the literature. First, exact linear resonance is not a universal prerequisite: precession resonance, quasi-resonance, and forced interactions show that detuning and dissipation can still permit large transfers (Lucas et al., 2016, Cazaubiel et al., 2019). Second, PSI is not the only relevant triadic process; more general triadic resonance instabilities, branch-crossing interactions, and collective three-mode instabilities all fall within the broader category of three-wave nonlinear tides. Third, the assumption that nonlocal, scale-separated triads dominate is not generally secure: the high-Reynolds-number internal-wave experiment found that most active triads were local. Finally, negative observational results do not imply absence of nonlinear tides; they may instead reflect noisy bispectra, insufficient scale separation, or dissipation levels that shift the system away from the idealized weakly nonlinear limit (Taebel et al., 2024, Frajka-Williams et al., 2014).

Taken together, the modern picture is that three-wave nonlinear tides are a unifying mechanism for energy transfer in oceanic, laboratory, and astrophysical tidal systems. What varies across contexts is not the existence of triads, but the geometry of the resonance manifold, the role of damping and detuning, the accessibility of phase-resolved diagnostics, and whether the observable outcome is ocean mixing, spectral cascade, orbital decay, tidal heating, or a measurable modulation of gravitational waves.