REVIEW 3 major objections 3 minor 41 references
Instabilities of internal gravity waves in the two-dimensional Boussinesq system
T0 review · 3 major / 3 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read The paper proves that small-amplitude internal gravity waves in an inviscid stably stratified Boussinesq fluid are linearly unstable, with an explicit linear-in-amplitude growth rate.
desk verdict First rigorous PSI-type instability for inviscid internal waves, with a proof that mostly holds up; the main caveat is a typo in Proposition 3.3 that needs fixing. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the resonant set $R_k$, the curve of Floquet parameters $\mu$ for which the unperturbed linearized operator $L_{\mu,0}$ has a double eigenvalue $\lambda_+ = \lambda^-_k(\mu) = \lambda^+_0(\mu)$ on the invariant subspace $H^1_k$ of functions whose Fourier coefficients are supported on wavevectors proportional to the primary wavevector $k$. Kato's similarity transformation and spectral projectors show that this double eigenvalue persists as an isolated pair for $\epsilon>0$, reducing the problem to the eigenvalues of a $2\times 2$ matrix. The product of the off-diagonal entries of that matrix yields the explicit growth-rate function $e(\mu)$ in equation (2.38), whose sign controls the instability.
What would settle it
Direct numerical diagonalization of the restricted operator $L_{\mu,\epsilon}|_{H^1_k}$ for a concrete case, say $k=(1,1)$, small $\epsilon$, and $\mu$ on the resonant branch $R_k$, should reproduce an unstable eigenvalue whose real part equals $\epsilon\sqrt{e(\mu)}$ at leading order; finding no such eigenvalue, or a leading rate that disagrees, would refute the theorem.
Extended reading notes
Core claim
The paper's central claim is that for every wavevector $k=(m,n)$ with $m,n>0$ and every sufficiently small amplitude $\epsilon$, the operator obtained by linearizing the two-dimensional inviscid Boussinesq equations in vorticity-stream form around the plane wave solution has at least one eigenvalue with strictly positive real part. The sharper statement is that, restricting to the invariant subspace of functions supported on wavevectors proportional to $k$ and choosing the Floquet parameter $\mu$ on the resonant set defined by $\Omega(k)-\Omega(k+\mu)=\Omega(\mu)$, the perturbed eigenvalues take the form $\lambda_+ + iO(\epsilon^2) \pm \epsilon\sqrt{e(\mu)} + O(\epsilon)$, where $e(\mu)$ is given explicitly and is positive in the small-$|\mu|$ and large-$|\mu|$ limits. Because $e(\mu)>0$, the real part of one of the eigenvalues is positive, an instability whose leading growth rate is linear in the amplitude $\epsilon$.
Load-bearing premise
The linear instability is a spectral point of the Floquet-Bloch operator on the torus, and the growing solution $e^{\lambda t} e^{i\mu x} v(x)$ is not square-integrable over $\mathbb{R}^2$, so the theorem does not provide an $L^2$ eigenfunction of the operator on the whole plane.
Editorial extensions
If this is right
- If the theorem is right, every sufficiently small-amplitude internal plane wave in the inviscid two-dimensional Boussinesq system is modulationally unstable, with growth rate proportional to amplitude.
- The explicit formula for $e(\mu)$ makes it possible to predict, for a given primary wave vector $k$, which resonant subharmonic pairs grow fastest.
- The instability persists all along each resonant branch, from small to large Floquet parameters, as the asymptotic formulas (2.40) and (2.41) show positivity in both limits.
- The method extends the rigorous modulational-instability toolkit to a non-Hamiltonian, reversible system with an unbounded perturbation, where classical perturbation theory requiring boundedness does not apply.
- The result rigorously confirms the physical PSI picture: a primary wave of frequency $\Omega(k)$ transfers energy to two secondary waves whose frequencies are, to leading order, half of $\Omega(k)$.
Reading between the lines
- One testable extension is to compute the spectrum of the truncated operator on $H^1_k$ numerically for a fixed $k$ and small $\epsilon$, and check that the leading growth rate matches $\epsilon\sqrt{e(\mu)}$; agreement to order $\epsilon$ would confirm the sharpness of the formula.
- Because the unstable Bloch mode is not square-integrable over $\mathbb{R}^2$, the established instability is spectral; whether nonlinear effects convert it into actual norm growth for localized data remains an open question that the paper does not address.
- The structure of the proof suggests that the resonance condition, rather than the specific form of $\Omega$, drives the instability, so an analogous result may hold for more general stable density profiles with a homogeneous dispersion relation.
- The relation in Proposition 6.1 between $e(\mu)$ and the physical interaction coefficients $I_\pm(\mu)$ indicates that the rigorous growth rate coincides with the weakly-nonlinear triad expansion in the small-Floquet limit, which could be used to calibrate reduced models of internal-wave energy transfer.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the two-dimensional inviscid Boussinesq system and proves a modulational instability result for small-amplitude internal plane waves. The authors linearize around an exact plane-wave solution, use Floquet-Bloch decomposition, restrict to the invariant subspace of harmonics of the primary wavevector, and show that for resonant Floquet parameters a purely imaginary double eigenvalue of the unperturbed operator splits into a pair with nonzero real part. The growth rate is computed explicitly as ϵ√e(µ) with an explicit function e(µ), and its small- and large-Floquet-parameter asymptotics are given. The paper also compares the resulting rate with the physical literature on parametric subharmonic instability and includes a Mathematica code for the algebraic computations.
Significance. If the result is correct, this is the first rigorous treatment of inviscid parametric subharmonic instability for internal gravity waves, a mechanism that is widely used in oceanography and experimental fluid mechanics. The explicit, parameter-free formula for the first-order growth rate and the verification of both small- and large-Floquet-parameter regimes are concrete and falsifiable predictions. The analytic setup, including the characterization of the resonant set in Lemma 2.5 and the spectral-isolation argument in Proposition 2.7, is a nontrivial adaptation of recent Stokes-wave techniques. The availability of machine-checkable code for the entanglement coefficients is a further strength. However, as detailed below, two load-bearing points need correction: a displayed algebra error in Proposition 3.3 that temporarily breaks the proof of Theorem 2, and an overstatement in Theorem 1 about the nature of the instability (spectral versus L²-eigenvalue instability).
major comments (3)
- [Proposition 3.3, Eqs. (3.18)-(3.19)] Equations (3.18) and (3.19) contain an extra factor ϵ on the right-hand side. The coefficients β1, β0, γ1 are ϵ-independent Taylor coefficients defined in (3.16), and the derivation in Section 5 (the displayed chain ending with 'proving formula (3.18)') actually computes (β1−γ1w)ϵ = ϵ/8·..., i.e. β1−γ1w = 1/8·... . As printed, the quantities b1 and b0 defined in the proof of Theorem 2 become O(ϵ), the off-diagonal entries in (3.20) become O(ϵ²), and the eigenvalue splitting ±ϵ√e(µ) in (2.37) does not follow. The final formulas (2.36) confirm that the intended identities are O(1). This is a local but load-bearing typo; it must be corrected and the surrounding computation reconciled.
- [Theorem 1 and Eqs. (2.6)-(2.7)] Theorem 1 claims that the linearized operator Lϵ on L²(R²) has an unstable eigenvalue. What is actually constructed is a Bloch eigenvalue of L_{µ,ϵ} on L²(T²); the corresponding spacetime function h(t,x)=e^{λt}e^{iµx}v(x) is not square-integrable over R², as the text itself notes after (2.6)-(2.7). Thus the result establishes spectral instability (Re σ(L_{µ,ϵ})>0 for some µ, equivalently the spectrum of Lϵ intersects the right half-plane) rather than existence of an L² eigenfunction. This is a standard accepted notion in modulational-instability theory, but Theorem 1 and the abstract should state this explicitly so that the claim is not overread as an L²-instability result.
- [Lemma 3.1] Lemma 3.1 is the main perturbative tool and is stated without proof. Because L0 is unbounded and the introduction emphasizes that the perturbation is not bounded, the analyticity of the spectral projectors P_{µ,ϵ} and the existence of the contour integral for small ϵ require justification. On the invariant subspace H¹_k this is straightforward: in view of (4.15), the jets L^{±k}_1 act on modes nk through the constants k⊥·µ, so L1 is actually bounded on H¹_k. Please add a proof or a precise Kato-type statement verifying the hypotheses, so that the reduction to the 2×2 matrix representation is fully justified.
minor comments (3)
- [Lemma 4.6, Eq. (4.19)] The reversibility identities in (4.19) appear to be missing complex conjugation or an overline; as typeset, the first identity is tautological. Please check the intended relations.
- [Appendix B, Step 2 near 0] For ℓ=−1 the leading coefficient displayed after (B.7) contains the factor m²−2n²; a parenthetical noting that this factor is nonzero for nonzero integers would remove an apparent edge case.
- [General notation] The notation r(ϵⁿ) is defined in the introduction, but in (2.37) the expression ±ϵ√e(µ)+r(ϵ) could be confused with the earlier real-valued remainders r(ϵ²); a brief reminder or a different symbol would improve readability.
Circularity Check
No circularity: the PSI growth rate is derived from the linearized Boussinesq operator with no fitted input, and the comparison with the physics literature is a post-hoc consistency check.
full rationale
The derivation is self-contained. The instability growth rate e(µ) is computed from the linearized Boussinesq operator: Theorem 2's b1(µ) and b0(µ) in (2.36) are obtained by expanding the 2x2 matrix (3.7) using the jets (4.3)-(4.4) and the entanglement coefficients (4.18), with no fitted constants. The resonance set R_k in (2.28) is a mathematical definition, not an assumption of instability; Lemma 2.5 and Proposition 2.7 only locate where the unperturbed eigenvalue has algebraic multiplicity two, and Theorem 2 then proves e(µ)>0 along the relevant parametrizations. The comparison with Dauxois et al. in Section 6 occurs after the derivation and is a consistency check, not an input. Self-citations to [9,10,15] supply technical vocabulary (entanglement coefficients, Kato transformations, the plane-wave ansatz), but the paper proves the lemmas it uses rather than importing the conclusion by citation. The apparent extra factor ϵ in (3.18)-(3.19) is a typographical/correctness issue: the Section 5 algebra and the final O(1) expressions for b1,b0 in (2.36) show the ϵ belongs on the left-hand side of those identities, so a reader must repair the displayed formulas, but this does not make the derivation circular. No load-bearing step reduces to its own input, and no fitted parameter is later relabelled as a prediction.
Assumptions & free parameters
assumptions (5)
- standard math Floquet-Bloch decomposition characterizes the L^2 spectrum of periodic-coefficient operators as the union of spectra of Bloch operators L_{mu,epsilon} on the torus.
- standard math Kato perturbation theory for isolated eigenvalues and spectral projectors applies to the unbounded operator family L_{mu,epsilon}: Y_k -> X_k.
- domain assumption The 2D inviscid Boussinesq system (1.1) is an adequate model for internal gravity wave instability.
- domain assumption The resonance condition (2.27) defines the set R_k of Floquet parameters studied; instability needs only one such parameter with e(mu)>0.
- domain assumption A Bloch eigenvalue with positive real part constitutes linearized instability of the original PDE.
Cite this review
Pith. "Pith review of Instabilities of internal gravity waves in the two-dimensional Boussinesq system." pith.science (2026). https://pith.science/paper/WZYO52QO
@misc{pith2026250710390,
author = {Pith},
title = {Pith review of: Instabilities of internal gravity waves in the two-dimensional Boussinesq system},
year = {2026},
howpublished = {\url{https://pith.science/paper/WZYO52QO}},
note = {Machine review of arXiv:2507.10390}
}
read the original abstract
We consider a two-dimensional, incompressible, inviscid fluid with variable density, subject to the action of gravity. Assuming a stable equilibrium density profile, we adopt the so-called Boussinesq approximation, which neglects density variations in all terms except those involving gravity. This model is widely used in the physical literature to describe internal gravity waves. In this work, we prove a modulational instability result for such a system: specifically, we show that the linearization around a small-amplitude travelling wave admits at least one eigenvalue with positive real part, bifurcating from double eigenvalues of the linear, unperturbed equations. This can be regarded as the first rigorous justification of the Parametric Subharmonic Instability (PSI) of inviscid internal waves, wherein energy is transferred from an initially excited primary wave to two secondary waves with different frequencies. Our approach uses Floquet-Bloch decomposition and Kato's similarity transformations to compute rigorously the perturbed eigenvalues without requiring boundedness of the perturbed operator - differing fundamentally from prior analyses involving viscosity. Notably, the inviscid setting is especially relevant in oceanographic applications, where viscous effects are often negligible.
Reference graph
Works this paper leans on
-
[1]
H. D. Abarbanel, D. D. Holm, J. E. Marsden, and T. S. Ratiu. Nonlinear stability analysis of stratified fluid equi- libria.Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences, 318(1543):349–409, 1986
work page 1986
-
[2]
T. Akylas and C. Kakoutas. Stability of internal gravity wave modes: from triad resonance to broadband instability. Journal of Fluid Mechanics, 961:A22, 2023
work page 2023
-
[3]
J. Bedrossian, R. Bianchini, M. C. Zelati, and M. Dolce. Nonlinear inviscid damping and shear-buoyancy instability in the two-dimensional Boussinesq equations.Comm. Pure Appl. Math., 76(12):3685–3768, 2023
work page 2023
-
[4]
D. Benielli and J. Sommeria. Excitation and breaking of internal gravity waves by parametric instability.Journal of Fluid Mechanics, 374:117–144, 1998
work page 1998
-
[5]
T. B. Benjamin. Internal waves of finite amplitude and permanent form.Journal of Fluid Mechanics, 25(2):241–270, 1966
work page 1966
-
[6]
T. B. Benjamin. Instability of periodic wavetrains in nonlinear dispersive systems.Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences, 299(1456):59–76, 1967
work page 1967
-
[7]
T. B. Benjamin. On the Boussinesq model for two-dimensional wave motions in heterogeneous fluids.Journal of Fluid Mechanics, 165:445–474, 1986
work page 1986
-
[8]
T. B. Benjamin and J. E. Feir. The disintegration of wave trains on deep water Part 1. Theory.Journal of Fluid Mechanics, 27(3):417–430, 1967. 22 R. BIANCHINI (†), A. MASPERO (⋆), AND S. PASQUALI (∗)
work page 1967
Show all 41 references
-
[9]
Berti, L
M. Berti, L. Corsi, A. Maspero, and P. Ventura. Infinitely many isolas of modulational instability for Stokes waves, 2024
2024
-
[10]
Berti, A
M. Berti, A. Maspero, and P. Ventura. Full description of Benjamin-Feir instability of Stokes waves in deep water. Inventiones mathematicae, 230(2):651–711, 2022
2022
-
[11]
Berti, A
M. Berti, A. Maspero, and P. Ventura. Benjamin–Feir instability of Stokes waves in finite depth.Archive for Rational Mechanics and Analysis, 247(5):91, 2023
2023
-
[12]
Berti, A
M. Berti, A. Maspero, and P. Ventura. Stokes waves at the critical depth are modulationally unstable.Communications in Mathematical Physics, 405(3):56, 2024
2024
-
[13]
Berti, A
M. Berti, A. Maspero, and P. Ventura. First isola of modulational instability of Stokes waves in deep water.EMS Surv. Math. Sci., 2025
2025
-
[14]
Bianchini, A.-L
R. Bianchini, A.-L. Dalibard, and L. Saint-Raymond. Near-critical reflection of internal waves.Analysis & PDE, 14(1):205–249, 2021
2021
-
[15]
Bianchini and T
R. Bianchini and T. Paul. The Triadic Resonant Instability of internal waves in stably stratified fluids.hal-04135944, 2023
2023
-
[16]
Bianchini and T
R. Bianchini and T. Paul. Reflection of internal gravity waves in the form of quasi-axisymmetric beams.Journal of Functional Analysis, 286(1):110189, 2024
2024
-
[17]
Bourget, T
B. Bourget, T. Dauxois, S. Joubaud, and P. Odier. Experimental study of parametric subharmonic instability for internal plane waves.Journal of Fluid Mechanics, 723:1–20, 2013
2013
-
[18]
Bourget, H
B. Bourget, H. Scolan, T. Dauxois, M. Le Bars, P. Odier, and S. Joubaud. Finite-size effects in parametric subhar- monic´ ainstability.Journal of Fluid Mechanics, 759:739–750, 2014
2014
-
[19]
T. J. Bridges and A. Mielke. A proof of the Benjamin-Feir instability.Archive for rational mechanics and analysis, 133:145–198, 1995
1995
-
[20]
R. P. Creedon, B. Deconinck, and O. Trichtchenko. High-frequency instabilities of stokes waves.Journal of Fluid Mechanics, 937:A24, 2022
2022
-
[21]
R. P. Creedon, H. Q. Nguyen, and W. A. Strauss. Transverse instability of stokes waves at finite depth.arXiv preprint arXiv:2408.07169, 2024
2024
-
[22]
R. P. Creedon, H. Q. Nguyen, and W. A. Strauss. Proof of the transverse instability of stokes waves.Annals of PDE, 11(1):4, 2025
2025
-
[23]
Dauxois, S
T. Dauxois, S. Joubaud, P. Odier, and A. Venaille. Instabilities of internal gravity wave beams.Annual review of fluid mechanics, 50(1):131–156, 2018
2018
-
[24]
R. E. Davis and A. Acrivos. The stability of oscillatory internal waves.Journal of Fluid Mechanics, 30(4):723–736, 1967
1967
-
[25]
Y. C. De Verdi` ere and L. Saint-Raymond. Attractors for two-dimensional waves with homogeneous hamiltonians of degree 0.Communications on Pure and Applied Mathematics, 73(2):421–462, 2020
2020
-
[26]
Desjardins, D
B. Desjardins, D. Lannes, and J.-C. Saut. Normal mode decomposition and dispersive and nonlinear mixing in stratified fluids.Water Waves, 3:153–192, 2021
2021
-
[27]
T. M. Elgindi and F. Pasqualotto. From Instability to Singularity Formation in Incompressible Fluids.arXiv: 2310.19780v1, 2023
2023 arXiv
-
[28]
Hibiya, M
T. Hibiya, M. Nagasawa, and Y. Niwa. Nonlinear energy transfer within the oceanic internal wave spectrum at mid and high latitudes.Journal of Geophysical Research: Oceans, 107(C11):28–1, 2002
2002
-
[29]
V. M. Hur and Z. Yang. Unstable stokes waves.Archive for Rational Mechanics and Analysis, 247(4):62, 2023
2023
-
[30]
Z. Jiao, L. M. Rodrigues, C. Sun, and Z. Yang. Small-amplitude finite-depth stokes waves are transversally unstable. arXiv preprint arXiv:2409.01663, 2024
2024 arXiv
-
[31]
Joubaud, J
S. Joubaud, J. Munroe, P. Odier, and T. Dauxois. Experimental parametric subharmonic instability in stratified fluids. Physics of Fluids, 24(4), 2012
2012
-
[32]
Jurja and K
C. Jurja and K. Widmayer. Long-time stability of a stably stratified rest state in the inviscid 2D Boussinesq equation. arXiv:2408.15154, 2024
2024 arXiv
-
[33]
Kato.Perturbation theory for linear operators
T. Kato.Perturbation theory for linear operators. Springer, 1966
1966
-
[34]
Kiselev, J
A. Kiselev, J. Park, and Y. Yao. Small scale formation for the 2-dimensional Boussinesq equation.Anal. PDE, 18(1):171– 198, 2025
2025
-
[35]
Lighthill
M. Lighthill. Contributions to the theory of waves in non-linear dispersive systems.IMA Journal of Applied Mathe- matics, 1(3):269–306, 1965
1965
-
[36]
R. R. Long. On the Boussinesq approximation and its role in the theory of internal waves.Tellus, 17(1):46–52, 1965
1965
-
[37]
MacKinnon, M
J. MacKinnon, M. H. Alford, O. Sun, R. Pinkel, Z. Zhao, and J. Klymak. Parametric subharmonic instability of the internal tide at 29 N.Journal of Physical Oceanography, 43(1):17–28, 2013
2013
-
[38]
J. W. McLean. Instabilities of finite-amplitude water waves.Journal of Fluid Mechanics, 114:315–330, 1982
1982
-
[39]
H. Q. Nguyen and W. A. Strauss. Proof of modulational instability of Stokes waves in deep water.Communications on Pure and Applied Mathematics, 76(5):1035–1084, 2023
2023
-
[40]
Varma, M
D. Varma, M. Mathur, and T. Dauxois. Instabilities in internal gravity waves.Math. Eng., 5(1):Paper No. 016, 34, 2023
2023
-
[41]
Zakharov
V. Zakharov. The instability of waves in nonlinear dispersive media.Sov. Phys. JETP, 24(4):740–744, 1967. (†)IAC, Consiglio Nazionale delle Ricerche (CNR), via dei Taurini 19, 00185, Rome, Italy Email address,†:roberta.bianchini@cnr.it (⋆)International School for Advanced Stud...
1967
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.