{"id":"92a0c62e-067a-49a8-baa6-d28fe7f9b2cd","arxiv_id":"2502.03630","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Using a hydrostatic Lagrangian transformation and H∞-calculus for the associated hydrostatic Lamé and Stokes operators, the authors prove local and global strong well-posedness for compressible primitive equations under several pressure laws, including with gravity.","lead":"This paper introduces a hydrostatic Lagrangian coordinate system, where the flow follows only the vertically averaged horizontal velocity, and uses it to prove new well-posedness results for the compressible primitive equations of atmospheric dynamics. The main outcomes are local strong solutions for large data and global strong solutions for small data, including the isothermal case with gravity.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The γ=2 with-gravity analysis is built on the wrong hydrostatic profile: ∂z(ρ²)=−ρ forces ρ=ξ−z/2, not ξ+z/2, so (2.6) and Theorem 2.2(b) do not cover the stated problem (1.1).","rationale":"The reader's weakest_assumption identifies precisely the most load-bearing point. The γ=2 hydrostatic balance is derived in the paper, not imported from an external convention, and the sign is wrong. The claim in Remark 2.1 that (2.6) is equivalent to (1.1) for γ=2 is therefore false as written. Since Theorem 2.2(b) is proved only by a sketch that uses the wrong coefficient c=1/(ξ0+1/2z), the local well-posedness claim for p=ρ², g=1 is not established. The correction has a non-trivial consequence: the physical density ξ−z/2 requires ξ>1/2 for positivity on the full layer, whereas assumption (A) permits M1 arbitrarily small. This is exactly the kind of hidden assumption that can break the maximal-regularity argument. Because the γ=1 global theorem (Theorem 2.3) and the no-gravity theorem (Theorem 2.4, though only sketched) rely on separate correct derivations, the overall verdict remains CONDITIONAL: the paper should be accepted only after the γ=2 section is corrected and the extra lower bound is stated. No stronger objection to the γ=1 argument was found in this pass.","tokens_in":25976,"tokens_out":11386,"duration_ms":107236,"concrete_test":"Verify by direct computation: integrate ∂z(ρ²)=−ρ with ρ(0)=ξ and replace (2.6) by the resulting system with ρ=ξ−z/2, pressure term (2ξ−z)∇Hξ, and continuity term −(1/2)z divH v. Then re-run the γ=2 fixed-point step (proof of Theorem 2.2(b)) with c=1/(ξ0−z/2). If the H∞-calculus/parameter-ellipticity argument for this operator requires the additional hypothesis ξ0>1/2, or if the nonlinear pressure estimate fails near 2ξ0=z, the stated theorem is false as written.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Confirmed: in §2.2 the paper states 'For γ=2 we obtain ρ = ξ + 1/2 z.' This is inconsistent with the hydrostatic equation (1.1)_3: ∂z p = −gρ with p=ρ² and g=1 gives 2ρ∂zρ = −ρ, hence ∂zρ = −1/2 and ρ = ξ − z/2. The subsequent system (2.6) — coefficients ξ+1/2z, pressure term (2ξ+z)∇Hξ, and the vertical-velocity reconstruction with +1/2 z divH v — is therefore not the γ=2, g=1 case of (1.1), contrary to Remark 2.1. The operator B in the proof of Theorem 2.2(b) uses c=1/(ξ0+1/2z), so the maximal-regularity argument is for the wrong operator. With the correct sign, the coefficient ξ0−z/2 is only uniformly positive on [0,1] if ξ0>1/2; assumption (A) only requires M1>0. Thus either Theorem 2.2(b) needs the corrected sign and an additional lower bound, or it does not apply to the original equations. The γ=1 results and the no-gravity results are not affected.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":26315,"tokens_out":12463,"duration_ms":108220,"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":[{"comment":"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).","section":"§2.2, Eq. (2.6) and Remark 2.1"},{"comment":"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.","section":"§6, Proof of Theorem 2.4"}],"minor_comments":[{"comment":"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.","section":"§2.2, Eq. (2.6)_1"},{"comment":"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.","section":"§2.1 and §2.2, Eqs. (2.4) and (2.6)"},{"comment":"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.","section":"§6, Eq. (6.15)"},{"comment":"In the introduction, the name 'Cao and Titi' appears with a typo as 'Cao ant Titi'.","section":"§1, references"}],"recommendation":"major_revision","confidential_remarks":"The sign error in the gamma=2 with-gravity part is localized and appears fixable by correcting rho = xi - z/2, adjusting (2.6), adding the condition xi0 > 1/2, and rerunning the perturbative argument for the corrected operator. The gamma=1 and no-gravity results are not affected. The sketch for Theorem 2.4 needs to be expanded into a real proof. I recommend major revision rather than rejection because the core strategy and the gamma=1 results seem defensible."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The hydrostatic Lagrangian transformation, which follows the vertically averaged horizontal velocity rather than the full three-dimensional flow, is a genuinely new device, and it earns its keep in the gamma=1 analysis. As far as I can tell, the global small-data strong well-posedness result for the isothermal compressible primitive equations with gravity (Theorem 2.3) is the first of its kind. The operator-theoretic core is real work: parameter ellipticity of the hydrostatic Lame operator, the H-infinity-calculus, and the fixed-point argument in Section 5 are laid out in enough detail to be checked. I did not find a circular dependency. The self-citations are to general tools, not to the target theorems, so that is not a problem.\n\nThe soft spot is the gamma=2 case, and it is not minor. From (1.1) with p = rho^2 and g = 1, the hydrostatic balance gives 2rho*partial_z rho = -rho, so rho = xi - z/2. The paper instead uses rho = xi + z/2 in (2.6) and throughout the proof of Theorem 2.2(b), including the operator B with c = 1/(xi0 + (1/2)z). That means Theorem 2.2(b), as stated, does not cover the original system (1.1) for gamma=2; it covers a different system with a sign flip. With the corrected sign, the coefficient xi0 - z/2 is positive only when xi0 > 1/2, which assumption (A) does not require. So either the sign is fixed and a lower bound is added, or the claim must be withdrawn. The gamma=1 and no-gravity results are unaffected.\n\nA smaller issue: Theorem 2.4 is explicitly not proved (\"We will not provide a detailed proof\"). The sketch is plausible, but the global part is asserted rather than demonstrated. For a paper of this size, that is acceptable if it is labeled as a short discussion, but the abstract sells it as a theorem.\n\nThe paper deserves a serious referee: the gamma=1 contribution is substantial and the method is reusable. The referee report should ask for the gamma=2 correction before acceptance. I would cite the gamma=1 result with the usual caveats and would bring the paper to reading group because the transformation itself is worth discussing.","headline":"Strong new Lagrangian tool and a solid gamma=1 global theorem, but the gamma=2 gravity case is invalidated by a sign error.","tokens_in":26823,"tokens_out":2953,"would_cite":true,"duration_ms":23917,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q35","35Q86","76N10","47D06"],"pacs":[],"model":"deepseek-v4-flash","headline":"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…","keywords":["compressible primitive equations","hydrostatic Lagrangian coordinates","strong well-posedness","hydrostatic Lame operator","hydrostatic Stokes operator","maximal regularity","H-infinity calculus","isothermal atmosphere"],"falsifier":"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.","tokens_in":25789,"feed_emoji":"🌍","tokens_out":12749,"duration_ms":103109,"temperature":0.7,"pith_summary":"Atmospheric and oceanic dynamics are modelled by the compressible primitive equations, but strong well-posedness for these equations has been much harder to obtain than for their incompressible counterpart. This paper introduces a hydrostatic Lagrangian coordinate change: the flow is generated by the vertical average of the horizontal velocity rather than by the full velocity, and the density and velocity are pulled back along this flow. After the change of variables, the hyperbolic continuity equation becomes a parabolic equation for the surface density, so the system can be studied with maximal-regularity theory for the linearised problem. The authors establish local strong well-posedness for large data under the pressure laws $p=\\rho$ and $p=\\rho^2$ with gravity, and for general monotone pressure laws without gravity, as well as global strong well-posedness for small data in the isothermal case $p=\\rho$ with gravity. A sympathetic reader would care because previous global strong results were only obtained without gravity, and this is a unified functional-analytic treatment that includes gravity.","feed_headline":"Hydrostatic Lagrangian method solves compressible primitive equations","feed_subtitle":"A change of variables along the averaged horizontal flow turns the system into a parabolic problem with exponential decay.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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)."],"supporting_citations":[{"why":"Supplies the Lp-theory of cylindrical boundary value problems, including parameter-ellipticity and bounded H-infinity calculus, used for the horizontal and vertical components of the hydrostatic Lame operator.","marker":"[31]"},{"why":"Provides the perturbation result for products of non-commuting operators that transfers the H-infinity calculus from the model operator to the hydrostatic Lame operator with variable coefficient.","marker":"[18]"},{"why":"Supplies the inf-sup condition and the two-dimensional Stokes regularity used to prove invertibility and higher regularity of the compressible hydrostatic Stokes operator.","marker":"[19]"},{"why":"Provides the product estimates and Sobolev embeddings for the horizontal-velocity nonlinearities in the local well-posedness proof.","marker":"[17]"},{"why":"Gives the R-boundedness and H-infinity calculus results for elliptic boundary value problems used for the vertical operator and for bootstrap arguments.","marker":"[6]"},{"why":"Supplies the classical H-infinity calculus result invoked to conclude the bounded H-infinity calculus for the shifted model operator in Lemma 4.1.","marker":"[5]"}],"fun_headline_variants":["Hydrostatic Lagrangian tames compressible primitive equations","Lagrangian hydrostatics crack primitive equations","Compressible primitive equations fall to Lagrangian hydrostatics","Global well-posedness via Lagrangian hydrostatics","Primitive equations solved using hydrostatic Lagrangian method"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Hydrostatic Lagrangian tames compressible primitive equations","Lagrangian hydrostatics crack primitive equations","Compressible primitive equations fall to Lagrangian hydrostatics","Global well-posedness via Lagrangian hydrostatics","Primitive equations solved using hydrostatic Lagrangian method"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000686,"raw_usage":{"total_tokens":3074,"prompt_tokens":874,"completion_tokens":2200,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":490,"completion_tokens_details":{"reasoning_tokens":2127}},"tokens_in":490,"tokens_out":2200,"duration_ms":15060,"temperature":1.0,"reasoning_tokens":2127,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T04:18:28.452836+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Nau, Lp-Theory of Cylindrical Boundary Value Problems","cited_arxiv_id":null,"evidence_quote":"Supplies the Lp-theory of cylindrical boundary value problems, including parameter-ellipticity and bounded H-infinity calculus, used for the horizontal and vertical components of the hydrostatic Lame operator."},{"cited_title":"Haller-Dintelmann, M","cited_arxiv_id":null,"evidence_quote":"Provides the perturbation result for products of non-commuting operators that transfers the H-infinity calculus from the model operator to the hydrostatic Lame operator with variable coefficient."},{"cited_title":"Hieber, T","cited_arxiv_id":null,"evidence_quote":"Supplies the inf-sup condition and the two-dimensional Stokes regularity used to prove invertibility and higher regularity of the compressible hydrostatic Stokes operator."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the product estimates and Sobolev embeddings for the horizontal-velocity nonlinearities in the local well-posedness proof."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the R-boundedness and H-infinity calculus results for elliptic boundary value problems used for the vertical operator and for bootstrap arguments."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the classical H-infinity calculus result invoked to conclude the bounded H-infinity calculus for the shifted model operator in Lemma 4.1."}],"review_version":1}