REVIEW 2 major objections 4 minor 1 cited by
The Lagrangian approach to the compressible primitive equations
T0 review · 2 major / 4 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read The paper proves local strong well-posedness for large data and global strong well-posedness for small data of the compressible primitive equations, using a hydrostatic Lagrangian coordinate change that converts the system into a…
desk verdict Strong new Lagrangian tool and a solid gamma=1 global theorem, but the gamma=2 gravity case is invalidated by a sign error. 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 object that carries the argument is the hydrostatic Lagrangian coordinate system: the flow $X(t,y_H)$ solves $\partial_t X = \bar v(t,X)$, where $\bar v$ is the vertical average of the horizontal velocity, and all functions are pulled back by this two-dimensional flow. This makes the continuity equation linear in the surface density $\zeta$ up to quadratic error terms, while the momentum equation becomes a parabolic equation whose leading part is the hydrostatic Lame operator $A_{\mathrm{HL}}$, a degenerate elliptic operator with coefficients depending on the transformed vertical coordinate $z' = (1-e^{-z})/\delta$. The linear theory rests on proving that $-A_{\mathrm{HL}} + \omega$ and the companion compressible hydrostatic Stokes operator $-A_{\mathrm{CHS}} + \omega$ have bounded $H^\infty$-calculi of angle $<\pi/2$, via cylindrical parameter-ellipticity and a perturbation result for products of non-commuting operators. On the mean-zero subspace, $A_{\mathrm{CHS}}$ is invertible and generates an exponentially stable semigroup, which is what the global small-data theorem needs.
What would settle it
Integrate $\partial_z p = -g\rho$ with $p=\rho^2$, $g=1$, and $\rho(x,y,0)=\xi(x,y)$: since $2\rho\,\partial_z\rho = -\rho$, one obtains $\rho = \xi - z/2$. Substituting the paper's profile $\rho = \xi + z/2$ into the balance gives $2(\xi+z/2)(1/2)=\xi+z/2$ on the left and $-(\xi+z/2)$ on the right, which are equal only when $\xi+z/2=0$. Checking the derivation of (2.6) against this one-line integration settles whether the $\gamma=2$ local well-posedness theorem applies to the stated $p=\rho^2$ system.
Extended reading notes
Core claim
The central claim is that the hydrostatic Lagrangian transform reduces the compressible primitive equations to a quasilinear parabolic system whose linearisation is the compressible hydrostatic Stokes operator $A_{\mathrm{CHS}}$ on $X_0 = H^{1,q}_{\mathrm{per}}(G) \times L^q(\Omega)$. The authors prove that $-A_{\mathrm{CHS}}$ admits a bounded $H^\infty$-calculus of angle strictly less than $\pi/2$, and that on the mean-zero subspace it is invertible and generates an exponentially stable semigroup. These spectral and functional-calculus facts yield maximal $L^p$-$L^q$ regularity on $\mathbb{R}_+$. Consequently, for data satisfying assumption (A), Theorem 2.2 gives a unique local strong solution to the transformed systems (2.4) and (2.6), and Theorem 2.3 gives a unique global strong solution in the isothermal case with gravity for data sufficiently close to a constant state $(\bar\xi,0)$, with exponential decay of the perturbations and the lower bound $\xi(t,x,y)\ge \bar\xi/2$. Theorem 2.4 states the analogous local and global results for general monotone pressure laws in the absence of gravity. The paper states in Remarks 2.1 and 2.5 that the transformed systems are equivalent to the original ones, so the well-posedness carries over.
Load-bearing premise
The $\gamma=2$ part of Theorem 2.2 depends on the identity $\rho = \xi + z/2$ obtained from the hydrostatic balance $\partial_z p = -g\rho$ with $p=\rho^2$ and $g=1$; direct integration gives $\rho = \xi - z/2$, so the transformed system (2.6) may not be equivalent to the original equations (1.1) as stated in Remark 2.1.
Editorial extensions
If this is right
- Local strong well-posedness for large data now holds for the compressible primitive equations in the isothermal case $\gamma=1$ and the quadratic case $\gamma=2$ with gravity, and for general increasing pressure laws without gravity, under periodic lateral boundary conditions.
- Small initial data near a constant state in the isothermal case with gravity yield a unique global strong solution that decays exponentially in time, with no vacuum formation since the surface density stays at least $\bar\xi/2$.
- The local and global theorems transfer to the flat torus $\mathbb{T}^2\times(0,1)$ and, for local well-posedness, to the infinite layer $\mathbb{R}^2\times(0,1)$, as stated in Remark 2.8.
- The authors state that the local result extends to the full power law $p=\rho^\gamma$ for $\gamma\ge 1$ by a modification of the $\gamma=2$ argument (Remark 2.6).
- Total mass is conserved and an explicit energy identity holds in the isothermal-with-gravity case and in the no-gravity general-pressure case (Remark 2.7).
Reading between the lines
- The approach avoids the pressure-coordinate assumption of earlier coupled atmosphere-ocean models, so the same machinery could in principle be applied to the interface problem where the surface pressure is not constant; this is an extension the authors do not carry out.
- The exponential decay rate in Theorem 2.3 is tied to the strictly negative spectral bound of the compressible hydrostatic Stokes operator; estimating that bound for realistic viscosity and density profiles would turn the qualitative theorem into quantitative atmospheric adjustment times.
- Because the hydrostatic Lame operator is already anisotropic, the same Lagrangian transform is a plausible tool for the anisotropic primitive equations with horizontal viscosity only, a regime relevant to ocean models.
- The local well-posedness claim for $p=\rho^2$ with gravity rests on a specific vertical profile for the density; if that profile is corrected, the same Banach fixed-point scheme may still work but for a slightly different transformed system, so the theorem as written should be checked against the corrected identity.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a hydrostatic Lagrangian transformation for the compressible primitive equations, in which the flow is defined along the vertical average of the horizontal velocity. The transformed equations are then treated as a quasilinear parabolic system, with linear theory based on an H-infinity calculus for the hydrostatic Lamé operator and the compressible hydrostatic Stokes operator. The main results are local strong well-posedness for large data for the isothermal (gamma=1) and quadratic (gamma=2) pressure laws with gravity, local and global strong well-posedness for small data near a constant state for gamma=1 with gravity, and local and global well-posedness for a general monotone pressure law without gravity. The proof strategy combines maximal Lp-Lq regularity, fixed-point contraction arguments, and the inversion of the Lagrangian transformation. The author's main claim is that this yields a new, systematic approach to strong well-posedness for the compressible primitive equations, including the first global small-data result in the isothermal case with gravity.
Significance. If the results were fully correct, the paper would make a substantial contribution. The hydrostatic Lagrangian viewpoint is novel and provides a unified operator-theoretic framework for both local and global well-posedness, and the global small-data theorem for gamma=1 with gravity appears to go beyond previous energy-method results. The detailed treatment of the hydrostatic Lamé operator and the compressible hydrostatic Stokes operator is a useful structural contribution, and the formal derivation is coherent in the gamma=1 and no-gravity cases. I find no circularity: the cited operator-theoretic results from [16,17,18,19] are prior general tools concerning H-infinity calculus and primitive equations, not the specific compressible well-posedness theorems being proved. However, the gamma=2 with-gravity reduction contains a sign error that invalidates Theorem 2.2(b) as stated, and the proof of Theorem 2.4 is only a sketch despite that theorem being one of the main results. These two points are load-bearing and need to be addressed before the paper can be accepted.
major comments (2)
- [§2.2, Eq. (2.6) and Remark 2.1] The asserted density profile rho = xi + (1/2)z is inconsistent with the hydrostatic balance in (1.1). With p = rho^2 and g = 1, equation (1.1)_3 gives 2 rho d_z rho = -rho, hence d_z rho = -1/2 and rho = xi - z/2. Consequently the averaged system (2.6), which uses coefficients xi + (1/2)z, pressure gradient (2xi+z) grad_H xi, and vertical-velocity reconstruction with +(1/2)z div_H v, is not equivalent to the original gamma=2, g=1 system, contrary to what is claimed in Remark 2.1. In the proof of Theorem 2.2(b), the operator B is defined with c = 1/(xi0 + (1/2)z), so the maximal-regularity argument is performed for the wrong operator. If the sign is corrected, the coefficient becomes xi0 - z/2, which is only uniformly positive on [0,1] when xi0 > 1/2; this condition is not implied by Assumption (A), which only requires M1 > 0. Thus Theorem 2.2(b) either needs to be corrected with an additional lower bound xi0 > 1/2, or it does not apply to the stated problem (1.1).
- [§6, Proof of Theorem 2.4] Theorem 2.4 is stated as a main result, covering both local strong well-posedness and global small-data well-posedness for a general pressure law without gravity, but its proof is not supplied. The text says 'We will not provide a detailed proof of Theorem 2.4' and then gives only a short discussion, including a change of variables zeta-tilde = P(zeta) for which the trace-space compatibility, the nonlinear estimates for the term P'(zeta+xi)(Z^T - I) grad_H zeta, and the global contraction argument are not verified. Since this theorem is one of the three principal results in the paper, the missing details are load-bearing and should be provided rather than deferred to a sketch.
minor comments (4)
- [§2.2, Eq. (2.6)_1] The first equation of (2.6) is stated on G x (0,T) but contains the vertical variable z along with div_H v; the averaged equation used later in the proof of Theorem 2.2(b) instead contains the term (1/2)∫_0^1 z div_H V dz. Please make the notation and the domain of the equation consistent.
- [§2.1 and §2.2, Eqs. (2.4) and (2.6)] In equations (2.4)_1 and (2.6)_1, the vertical-average bars on v appear to be missing in some places; for example, the terms v·grad_H xi and xi div_H v should presumably involve the vertical average of v. Please ensure that the averaged quantities are marked consistently throughout.
- [§6, Eq. (6.15)] The right-hand side of (6.15) reads 'f1 - (1/lambda) grad_H f2', but f1 is scalar-valued while f2 is vector-valued; from (6.8) the correct expression should be f2 - (1/lambda) grad_H f1. Please check and correct this and the related formulas in Lemma 6.3.
- [§1, references] In the introduction, the name 'Cao and Titi' appears with a typo as 'Cao ant Titi'.
Circularity Check
No load-bearing circularity: the derivation is a contraction argument around linear operators built from the transformed PDE, with only routine non-circular use of prior operator-theoretic results; the γ=2 with-gravity section contains a separate sign error that is a correctness issue, not a circularity.
full rationale
I find no circular step that reduces a claimed prediction to an input by construction. The main theorems are proved by transforming (1.1) via the hydrostatic Lagrangian flow, linearizing the transformed systems (2.4), (2.6), and (2.9), building the hydrostatic Lamé operator AHL and the compressible hydrostatic Stokes operator ACHS, proving maximal regularity, and then closing with contraction maps. The linearized operators are derived from the PDE rather than defined to force the conclusion: for example, AHL is the Lq-realization of A(x,y,D)v = µa∆Hv + ∂z(µb∂zv) + µ'a∇HdivHv with a = 1/((1−δz)ξ0) and b = (1−δz)/(δ²ξ0), and the maximal-regularity/resolvent estimates are established by parameter ellipticity and Babuška–Brezzi/Lax–Milgram arguments. The proof uses [18, Theorem 3.1], [31], [5], [6], and [16,17,19]. Some of these involve a coauthor of the present paper ([16,17,18,19] include Hieber), but they are general H∞-calculus, cylindrical parameter-ellipticity, function-space, and incompressible-primitive-equation tools; none of them assumes the compressible well-posedness theorem being proved. Hence the self-citations are not load-bearing circular premises, and the central claims retain independent content. There are also no fitted parameters renamed as predictions: δ is the smallness of the initial perturbation, and all constants come from resolvent, embedding, and product estimates. Separately, and independently of circularity, §2.2 contains a correctness error: the text states "For γ=2 we obtain ρ=ξ+1/2z," but integrating the hydrostatic balance ∂z(ρ²) = −ρ with g=1 gives ∂zρ = −1/2, hence ρ = ξ−z/2. Consequently system (2.6) and the operator B in the proof of Theorem 2.2(b) are not equivalent to (1.1) with γ=2 as claimed in Remark 2.1; with the corrected sign one would need ξ0>1/2 to keep ξ0−z/2 positive under Assumption (A). This affects the validity of Theorem 2.2(b), but it is an inconsistency rather than a circular reduction, so it does not raise the circularity score.
Assumptions & free parameters
assumptions (5)
- standard math Bounded H∞-calculus and maximal regularity for parameter-elliptic cylindrical boundary value problems.
- standard math Babuška-Brezzi and Lax-Milgram theory for the resolvent of the Stokes-type operator.
- domain assumption Initial density is bounded away from zero: M1 ≤ ξ0 ≤ M2 and ξ0 ≥ c > 0.
- domain assumption Pressure law conditions: γ=1, γ=2, or P with c1 ≤ P' ≤ c2.
- ad hoc to paper For γ=2 with gravity, ρ = ξ + z/2 follows from hydrostatic balance.
Cite this review
Pith. "Pith review of The Lagrangian approach to the compressible primitive equations." pith.science (2026). https://pith.science/paper/LTSK56KD
@misc{pith2026250203630,
author = {Pith},
title = {Pith review of: The Lagrangian approach to the compressible primitive equations},
year = {2026},
howpublished = {\url{https://pith.science/paper/LTSK56KD}},
note = {Machine review of arXiv:2502.03630}
}
read the original abstract
This article develops the hydrostatic Lagrangian approach to the compressible primitive equations. A fundamental aspect in the analysis is the investigation of the compressible hydrostatic Lam\'{e} and Stokes operators. Local strong well-posedness for large data and global strong well-posedness for small data are established under various assumptions on the pressure law, both in the presence and absence of gravity.
Forward citations
Cited by 1 Pith paper
-
Rigorous justification of the hydrostatic-incompressible approximation for weakly stratified isothermal flow
Weakly stratified isothermal compressible Navier-Stokes flows with small Mach number and aspect ratio converge rigorously, for general initial data, to the incompressible primitive equations.
Reference graph
Works this paper leans on
-
[1]
W. Arendt, Ch. Batty, M. Hieber, F. Neubrander, Vector-Valued Laplace Transforms and Cauchy Problems . Mono- graphs in Mathematics, vol. 96, 2nd-edition, Birkhäuser, Basel, 2011
work page 2011
-
[2]
Amann, Linear and Quasilinear Parabolic Problems
H. Amann, Linear and Quasilinear Parabolic Problems . Monographs in Mathematics, vol. 89, Birkhäuser, 1995
work page 1995
-
[3]
Amann, Linear and Quasilinear Parabolic Problems
H. Amann, Linear and Quasilinear Parabolic Problems. Vol. II. Functi on Spaces. Monographs in Mathematics, vol. 106, Birkhäuser/Springer, Cham, 2019
work page 2019
- [4]
-
[5]
R. Denk, G. Dore, M. Hieber, J. Prüss, A. Venni, New though ts on old results of R. T. Seeley. Math. Ann. 328 (2004), 545–583
work page 2004
-
[6]
R. Denk, M. Hieber, J. Prüss, R-Boundedness, Fourier Multipliers and Problems of Ellipti c and Parabolic Type . Mem. Amer. Math. Soc. 788 2003
work page 2003
-
[7]
R. Denk, M. Hieber, J. Prüss, Optimal Lp-Lq-estimates for parabolic boundary value problems with inho mogeneous data. Math. Z. 257 (2007), 193–224
work page 2007
-
[8]
Danchin, A Lagrangian approach for the compressible N avier-Stokes equations
R. Danchin, A Lagrangian approach for the compressible N avier-Stokes equations. Ann. de l’Institut Fourier 64 (2014), 753-791
work page 2014
Show all 33 references
-
[9]
Danchin, Fourier analysis methods for the compressib le Navier-Stokes equations
R. Danchin, Fourier analysis methods for the compressib le Navier-Stokes equations. Handbook of Math. Anal. in Mechanics of Viscous Fluids , Y. Giga, A. Novotny (eds.), Springer, 2018, 1843–1903. COMPRESSIBLE PRIMITIVE EQUATIONS 21
2018
-
[10]
Danchin, P
R. Danchin, P. Mucha, Compressible Navier-Stokes equa tions with ripped density. Comm. Pure Appl. Math. 76 (2023), 3437–3492
2023
-
[11]
R. Denk, T. Nau, Discrete Fourier multipliers and cylin drical boundary-value problems. Proc. Roy. Soc. Edinburgh 143 (2013), 1163–1183
2013
-
[12]
Enomoto, Y
Y. Enomoto, Y. Shibata, Global Existence of Classical S olutions and Optimal Decay Rate for Compressible Flows via the Theory of Semigroups. Handbook of Math. Anal. in Mechanics of Viscous Fluids , Y. Giga, A. Novotny (eds.), Springer, 2018, 2085–2181
2018
-
[13]
Ersoy, T
M. Ersoy, T. Ngom, Existence of a global weak solution to compressible primitive equations. C.R. Acad Sci. 350, (2012), 379–382
2012
-
[14]
Ersoy, T
M. Ersoy, T. Ngom, M. Sy, Compressible primitive equati ons: formal derivation and stability of weak solutions. Nonlinearity 24, (2011), 79–96
2011
-
[15]
Feireisl, A
E. Feireisl, A. Novotný and H. Petzeltová, On the existe nce of globally defined weak solutions to the Navier-Stokes equations, J. Math. Fluid Mech. 3 (2001), no. 4, 358–392
2001
-
[16]
Y. Giga, M. Gries, M. Hieber, A. Hussein, T. Kashiwabara , Bounded H∞ − calculus for the hydrostatic Stokes operator on Lp− spaces and applications. Proc. Amer. Math. Soc. 145 (2017), 3865–3876
2017
-
[17]
Y. Giga, M. Gries, M. Hieber, A. Hussein, T. Kashiwabara , Analyticity of solutions to the primitive equations. Math. Nachr. 293 (2020), 284–304
2020
-
[18]
Haller-Dintelmann, M
R. Haller-Dintelmann, M. Hieber, H∞ -calculus for products of non-commuting operators. Math. Z. 251 (2005), 85–100
2005
-
[19]
Hieber, T
M. Hieber, T. Kashiwabara, Global strong well-posedne ss of the three dimensional primitive equations in Lp-spaces. Arch. Rational Mech. Anal. 221 (2016), 1077–1115
2016
-
[20]
Hoff, Dynamics of singularity surfaces for compressi ble, viscous flows in two space dimensions
D. Hoff, Dynamics of singularity surfaces for compressi ble, viscous flows in two space dimensions. Comm. Pure Appl. Math. 55 (2002), 1365–1407
2002
-
[21]
Hytönen, J
T. Hytönen, J. van Neerven, M. Veraar, L. W eis, Analysis in Banach spaces, vol. II , Springer, Cham, 2016
2016
-
[22]
Kukavica, M
I. Kukavica, M. Ziane, On the regularity of the primitiv e equations of the ocean. Nonlinearity 20 (2007), 2739–2753
2007
-
[23]
Lions, R
J.-L. Lions, R. Temam, S.H. W ang, New formulations of th e primitive equations of atmosphere and applications. Nonlinearity 5 (1992), 237–288
1992
-
[24]
Lions, R
J.-L. Lions, R. Temam, S.H. W ang, On the equations of the large-scale ocean. Nonlinearity 5 (1992), 1007–1053
1992
-
[25]
Lions, R
J.-L. Lions, R. Temam, S.H. W ang, Mathematical theory f or the coupled atmosphere-ocean models. (CAO III). J. Math. Pures Appl. (9) 74 (1995), 105-163
1995
-
[26]
Lions, Mathematical topics in fluid mechanics
P.-L. Lions, Mathematical topics in fluid mechanics. Vo l. 2, Oxford Lecture Series in Mathematics and its Applications, Oxford Science Publications , 10 , Oxford Univ. Press, New York, 1998
1998
-
[27]
X. Liu, E. Titi, Global existence of weak solutions to th e compressible primitive equations of atmospheric dynamic s with degenerate viscosities. SIAM J. Math. Anal. 51 (2019), 1913–1964
2019
-
[28]
X. Liu, E. Titi, Zero Mach number limit of the compressib le primitive equations part I: W ell-prepared initial data. Arch. Rational Mech. Anal. 238 (2020), 705–747
2020
-
[29]
X. Liu, E. Titi, Local well-posedness of strong solutio ns to the three-dimensional compressible primitive equati ons, Arch. Rational Mech. Anal. 241 (2021), 729–764
2021
-
[30]
Matsumura, T
A. Matsumura, T. Nishida, Initial boundary value probl ems for the equations of motion of compressible viscous and heat-conductive fluids, Comm. Math. Phys. 89 (1983), 445–464
1983
-
[31]
Nau, Lp-Theory of Cylindrical Boundary Value Problems
T. Nau, Lp-Theory of Cylindrical Boundary Value Problems . Springer Spektrum, Wiesbaden, 2012, xvi+188 pp
2012
-
[32]
Prüss, G
J. Prüss, G. Simonett, Moving Interfaces and Quasilinear Parabolic Evolution Equ ations. Monographs in Mathematics, vol. 105, Birkhäuser, 2016
2016
-
[33]
Triebel, Interpolation Theory, Function Spaces, Differential Operat ors
H. Triebel, Interpolation Theory, Function Spaces, Differential Operat ors. North-Holland, 1978. Technische Universität Darmstadt, F achbereich Mathematik, Schlossgartenstr. 7, 64289 Darmstadt, Ger- many Email address : hieber@mathematik.tu-darmstadt.de Department of Pure and A...
1978
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.