{"id":"218ade07-95d8-4696-b3de-6f7c74b9a880","arxiv_id":"2607.27477","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"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.","lead":"This paper proves that, when the Mach number, aspect ratio, and Froude number are all small and equal, solutions of the isothermal compressible Navier-Stokes equations exist uniformly in the small parameter and converge, up to a subsequence, to the incompressible primitive equations used in ocean and atmosphere modeling. The result holds without smallness restrictions on the initial waves, but it is proved for an eddy-viscosity version of the equations rather than the origina","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Main theorem proves convergence only for the eddy-viscosity system (1.1), not the physically derived anisotropic-viscosity system (2.11); Remark 3.4 explicitly leaves that model untreated.","rationale":"The reader's weakest assumption correctly identifies the central scope limitation: the rigorous theorems are stated for the eddy-viscosity system (1.1), not for the anisotropic-viscosity system (2.11) that arises from the nondimensionalization of the physical Navier-Stokes equations. This is not a hidden defect—Remark 3.4 states it explicitly—but it is decisive for the central claim as advertised. The abstract says solutions of 'the compressible Navier-Stokes system' converge to the primitive equations, yet the proof only covers the system with isotropic eddy viscosity. Since the limit system is the same for both viscosity models, the asymptotic limit may well hold for (2.11), but the uniform existence and acoustic estimates are not established there. The paper is otherwise a serious and technically plausible contribution: the linear three-wave analysis is explicit, the energy and acoustic structures are spelled out, and the main estimates appear internally consistent for the eddy-viscosity model. No fatal contradiction was found. The appropriate verdict remains CONDITIONAL: the authors should either prove the result for the physical viscosity (2.11) or clearly frame the paper's contribution as a theorem for the eddy-viscosity regularization.","tokens_in":29949,"tokens_out":33227,"duration_ms":292280,"concrete_test":"Re-derive the acoustic energy estimate in Section 3.4 with the viscosity terms of (2.11) instead of the eddy viscosity. In (3.31c), replace ε²ΔW/ρ with ρ^{-1}[div_h(ε²∇_h W + ∂_z V) + 2∂_zz W + ∂_z(div_h V + ∂_z W)] and track the term ε∂_t W through (3.44)-(3.56). If the resulting dissipation in (3.47)-(3.48) is non-coercive in the horizontal direction—i.e., if the ε²Δ_h W term vanishes at leading order and cannot control ∥ε∂_t W∥—then Theorem 3.1 does not extend to (2.11), confirming that the advertised physical model is not covered.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The advertised result—justification of the hydrostatic-incompressible approximation for the compressible Navier-Stokes system—is proved only for the modified eddy-viscosity system (1.1), where the viscous stress is replaced by the isotropic operators Δv and Δw. The system obtained by nondimensionalizing the physical compressible Navier-Stokes equations is (2.11), whose viscosity tensor is anisotropic and, in the vertical momentum equation, contains the degenerate horizontal term div_h(ε²∇_h w) (see (2.11c)). Remark 3.4 states that this degeneracy obstructs the acoustic-wave estimates, and the paper does not treat (2.11). This gap is load-bearing because the uniform-in-ε estimates in Section 4, especially the acoustic energy estimate Lemma 3.4, rely on the isotropic ε²Δw term in (3.2c); without it the vertical acoustic wave lacks horizontal dissipation and the control of ε∂_t w / high vertical derivatives used in (4.9) and Proposition 4.4 does not close. The abstract and Theorem 1.1 nonetheless state the result for 'the compressible Navier-Stokes system,' which a reader would identify with (2.11) from Section 2.1, not with the modified model (1.1). The convergence theorem is internally coherent for (1.1), but the physically motivated system (2.11) is outside the theorem's scope.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the singular limit, as σ=Ma=Fr=ε→0, of a dimensionless isothermal compressible Navier-Stokes system in a periodic cylinder, with the goal of rigorously deriving the incompressible hydrostatic primitive equations. After reformulating the density through q=g+(1/ε)logρ, the authors identify a three-scale linear wave structure: a slow divergence-free mode, a fast horizontal acoustic mode, and a very fast vertical acoustic mode. The main results, Theorem 3.1 and Theorem 3.2, state uniform-in-ε existence of strong solutions on an ε-independent time interval and convergence, up to subsequence, to a solution of the primitive equations (1.2). The proof combines a weighted energy estimate, a second-order acoustic-wave estimate, and a compactness argument using non-orthogonal ε-dependent projections. The paper also contains formal asymptotics (Section 2) and two linear instability analyses (Section 6) indicating why the inviscid and general-gravity cases are excluded.","tokens_in":30376,"tokens_out":12039,"duration_ms":119149,"significance":"If correct, the result is a substantive contribution to the hydrostatic/incompressible limit for weakly stratified compressible flows: it goes beyond well-prepared data, handles large initial waves, and identifies the three-scale projection structure caused by the anisotropic aspect ratio. The paper is unusually explicit about its limitations: Remark 3.4 states that the physical anisotropic viscosity is not treated, and Section 6 gives concrete linear mechanisms for non-convergence in the inviscid and constant-gravity-gradient regimes. These admissions are a strength, but they also delimit the actual scope of the proof. The main theorems are internally coherent for the eddy-viscosity system (1.1), but the advertised claim for 'the compressible Navier-Stokes system' is broader than what is proved.","major_comments":[{"comment":"The central theorems are proved for the eddy-viscosity system (1.1), not for the physical anisotropic system (2.11) obtained by nondimensionalization. This is not a harmless simplification: in the vertical momentum equation (2.11c) the horizontal viscosity is the degenerate operator div_h(ε²∇h w), whereas the analysis uses the isotropic ε²Δw in (3.2c). The acoustic energy estimate Lemma 3.4 and the elliptic bound (4.9) rely on the isotropic Laplacian in (3.2c), and Remark 3.4 explicitly says the anisotropic case is not treated. The abstract and Theorem 1.1 should therefore be qualified: the result is for the eddy-viscosity model (1.1), or the anisotropic case must be handled.","section":"§3.1 (Remark 3.4), §2.1 (2.11), Theorems 1.1/3.1"},{"comment":"The assumption ∂z g, ∂zz g, ∂zzz g = O(ε) is load-bearing, not a harmless normalization. In the dimensionless physical setting of Section 2.1, ∂z g is O(1) for a standard gravity potential. Section 6.2 shows that for constant ∂z g=N≠0 the linearized vertical acoustic modes satisfy Re λ ∼ |N|/(2ε)>0, so the claimed convergence can fail outside the assumed weak-gravity regime. Thus the phrase 'weakly stratified' in the title and the abstract must be understood as this very specific O(ε) gravity-gradient condition, and the paper should state clearly that the result does not apply to the usual constant-gravity case.","section":"Theorem 3.1 (3.4), §6.2"},{"comment":"The final step from the differential inequality dE_total/dt + ½D_total ≤ H(E_total)(1+D_total^{1/2}) to the uniform bound (4.32) is not justified as written. Since H(E_total) is not uniformly bounded in ε or t, the D_total^{1/2} term cannot simply be absorbed into the left-hand side. A bootstrap or nonlinear Gronwall argument, with a time T depending on E0, is needed to control E_total before absorption. This is likely fixable, but it is a necessary step for the claimed uniform existence interval.","section":"§4.3, Eqs. (4.31)–(4.32)"}],"minor_comments":[{"comment":"The statement that ∂z q(0)=O(ε) 'is not a restriction on the initial data' is misleading. It is a compatibility/uniform-boundedness condition for the vertical acoustic energy; without it the data are singular in the ε→0 limit. The wording should be adjusted.","section":"Remark 3.5"},{"comment":"Several key estimates (e.g., I_1,...,I_8, Eqs. (4.22)–(4.30)) are summarized as 'tedious but straightforward.' For a rigorous journal, at least the nontrivial terms should be displayed with the needed Sobolev embeddings.","section":"§4.2, Proposition 4.4"},{"comment":"There are small presentational issues: 'consistent of' should be 'consisting of'; 'H¨older' and 'Cauchy-Schwartz' should be 'Hölder' and 'Cauchy-Schwarz'; and the abstract's phrase 'without any restriction on the size of the initial waves' should be reconciled with the boundedness assumption E0≤M in (3.4).","section":"Abstract and text"},{"comment":"The proof of L²-orthogonality of Q_va^{ker} and Q_va^{⊥} is compressed. A short display of the integration by parts in z would clarify the claim.","section":"Section 5, Lemma 5.2"}],"recommendation":"major_revision","confidential_remarks":"The mathematical core for the eddy-viscosity model is promising and the paper is honest about its limitations, but the main theorems overclaim relative to the physical system (2.11) derived from the Navier-Stokes viscosity tensor. The weak-gravity assumption is also more restrictive than the title suggests. These issues require rephrasing or extension before acceptance. I recommend major revision rather than rejection, since the eddy-viscosity result appears defensible and the scope can be made accurate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe paper does something genuinely hard: it proves uniform-in-ε existence for a three-scale compressible system and a subsequential limit to the primitive equations, without smallness on the initial waves, by tracking a non-orthogonal three-wave decomposition with ε-dependent projectors. That part is new and is a real step beyond the regular-projection analyses in [5,34,60] and the separate hydrostatic or incompressible limits in [15,47,48,56]. The linear wave analysis in Sections 3.2–3.4 is coherent, the weighted energy with e^{−εg} is clever, and the closing estimate (4.31) is plausible. Section 6, which exhibits instabilities in the inviscid case and for constant gravity, is honest and useful.\n\nThe soft spots are real, and the biggest one is exactly the stress-test note. The system obtained from the nondimensionalized compressible Navier–Stokes equations, (2.11), has anisotropic viscosity with a degenerate horizontal term ε²∇_h w in the vertical momentum equation. Remark 3.4 states that this degeneracy obstructs the acoustic estimates, so the entire proof is done for the modified eddy-viscosity system (1.1). The theorem statements themselves are for (1.1), which is the eddy-viscosity model; the abstract and title, however, say “compressible Navier–Stokes system” without qualification, and a reader naturally identifies that with (2.11). Since the acoustic-wave estimates in Lemma 3.4 and Proposition 4.4 rely on the isotropic ε²Δw, the advertised physical result is not proved. It is a theorem about a model with eddy viscosity, not about the physically derived viscosity tensor. The authors flag this in Remark 3.4, which is to their credit, but it remains a load-bearing gap between claim and theorem.\n\nTwo smaller issues. First, the weak-gravity condition ∂_z g = O(ε) is artificial; Section 6.2 shows that for constant N the linearized problem is unstable as ε→0, so it is needed within their framework, but it restricts the physical regime. Second, the abstract says solutions converge; Theorem 3.2 only gives subsequential convergence in (3.7), a mismatch in emphasis though not a flaw in the theorem.\n\nNet: the technical core is serious and the proof strategy is worth understanding. But the central advertised claim — justification of the hydrostatic-incompressible approximation for the compressible Navier–Stokes system with the natural anisotropic viscosity — is not established. A referee should be sent, with the expectation that either the anisotropic case is treated or the title and abstract are rewritten to claim the eddy-viscosity result. I would cite it for the eddy-viscosity theorem, with the caveat, and I would bring it to reading group if the group works on singular limits.\n\nBest regards","headline":"A genuinely hard three-scale singular-limit proof, but the advertised convergence is for the eddy-viscosity model, not the anisotropic system the nondimensionalization actually yields.","tokens_in":30769,"tokens_out":3938,"would_cite":true,"duration_ms":33529,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q30","35Q86","76D05","76N10","76N30"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that, for weakly stratified isothermal compressible flows in the small-Mach, small-aspect-ratio regime, solutions exist uniformly in the small parameter and converge to the incompressible primitive equations, with no restr","keywords":["Singular limit","Incompressible limit","Hydrostatic approximation","Primitive equations","Isothermal flow","Acoustic waves","Three-scale wave separation","Weak stratification"],"falsifier":"Simulate (1.1) with data satisfying (3.4) for ε = 10^{-2}, 10^{-3}, 10^{-4}; if the H^2 error between the horizontal velocity and a reference solution of (1.2) does not decay to zero (or the uniform bound (3.6) fails for some ε), the convergence theorem is false. As a complementary check, impose a small vertical acoustic perturbation and verify the dispersion relation (3.12): the vertical acoustic frequency must scale as |k_3|/ε², not ε^{-1}; a different scaling would falsify the three-wave separation.","tokens_in":29849,"feed_emoji":"🌊","tokens_out":9564,"duration_ms":87063,"temperature":0.7,"pith_summary":"This paper sets out to justify the hydrostatic-incompressible approximation—the passage from the isothermal compressible Navier–Stokes equations to the incompressible primitive equations—in the physical regime where the Mach number, Froude number, and vertical-to-horizontal aspect ratio are all of the same small order ε and stratification is weak. The authors claim that, for any fixed small ε, the system admits a unique strong solution on a time interval that does not shrink as ε→0, with no restriction on the size of the initial wave components. They further claim that, as ε→0, the solutions converge (along a subsequence) to the unique solution of the primitive equations (1.2). The reason to care: the primitive equations are the standard model for large-scale ocean and atmosphere dynamics, and this gives a proof that the model is the actual limit of the compressible system in this regime.","feed_headline":"Weak-stratified flow provably reduces to primitive equations","feed_subtitle":"Uniform existence on a fixed time interval and convergence to the ocean–atmosphere model, for general initial data.","key_machinery":"The machinery is the wave-decomposition framework of Definition 3.1: the horizontal acoustic projection P_ha, built on ∇_h Δ_h^{-1} div_h; the vertical acoustic projection P^ε_va, built on the ε-dependent elliptic operator (ε² Δ_h + ∂_zz)^{-1}; and the slow-wave projection P^ε_σ. These projections split the linear system (3.8) into the three separated time scales. The nonlinear analysis then runs on two lemmas: the weighted energy estimate (Lemma 3.3), where the factor e^{-εg} removes the singular gravity coupling, and the acoustic-wave energy estimate (Lemma 3.4), which provides the missing control of ∂_z q/ε, ε∂_t q, and w via equation (3.3).","core_discovery":"The central claim is Theorem 3.1 and Theorem 3.2: for the eddy-viscosity system (1.1), under the symmetry (SYM) and the weak-stratification hypothesis ∂_z g, ∂_zz g, ∂_zzz g = O(ε), initial data with bounded energy functional (3.4) yield unique strong solutions on an ε-independent time interval with the uniform bound (3.6); and as ε→0 the slow part of the solution converges to a strong solution of the incompressible primitive equations (1.2). The mechanism is a three-wave decomposition of the linearized dynamics—a slow/mean wave, a fast horizontal acoustic wave at frequency O(1/ε), and a very fast vertical acoustic wave at frequency O(1/ε²)—whose projection operators are non-orthogonal and d","pith_inferences":["Editorial: The main theorem does not apply to the anisotropic viscosity system (2.11) derived from the physical scaling; because Remark 3.4 identifies the degenerate horizontal viscosity in the vertical momentum equation as the obstruction, a natural test is whether a modified weighted estimate can restore the acoustic bound and extend the result to the original viscosity.","Editorial: The ε-dependent, non-orthogonal wave projections introduced here may serve as a template for other anisotropic singular limits, where the natural slow/fast/very-fast split cannot be made ε-independent.","Editorial: The instability calculations in Section 6 suggest a sharp borderline: one could try to prove that if ∂_z g decays more slowly than O(ε), the primitive-equations limit fails for the same class of data, making the paper's gravity condition necessary as well as sufficient.","Editorial: A concrete, cheap test of the three-scale separation would be a spectral simulation of the linear system (3.8) measuring the vertical acoustic frequency at O(ε^{-2}) for a fixed mode; agreement with (3.12) would independently confirm the central mechanism."],"forward_implications":["The incompressible primitive equations (1.2) are the rigorous asymptotic limit of the isothermal compressible Navier–Stokes system (1.1) in the weakly stratified, small-Mach, small-aspect-ratio regime, for general initial data without restrictions on wave size.","For any fixed small ε the system is uniformly well-posed on an ε-independent time interval, with the quantitative bound (3.6); this makes the singular limit amenable to numerical schemes that filter fast waves.","The weak-stratification condition ∂_z g = O(ε) is more than technical: Section 6.2 shows that with ∂_z g = N ≠ 0 the vertical acoustic mode grows at rate ≈ |N|/(2ε), so the leading-order limit is not the primitive equations for such data.","Viscosity is indispensable: the inviscid version (Section 6.1) contains growing modes with λ = 1/ε, so no convergence can be expected without dissipation.","The isentropic (γ > 1) case is not covered: Remark 4.1 shows the acoustic estimate fails to close there, leaving the extension to non-isothermal or isentropic flows open."],"fun_headline_variants":["Rigorous proof: weak stratification yields primitive equations","Primitive equations proven for weakly stratified flow","Weak-stratified limit: Navier-Stokes to ocean-atmosphere model","Hydrostatic-incompressible limit validated for weak stratification"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The central theorem applies to the modified isotropic eddy-viscosity system (1.1), not to the physically derived anisotropic viscosity (2.11); if that artificial viscosity is the reason the acoustic estimates close, the literal 'compressible Navier–Stokes' justification remains incomplete.","fun_headline_variants_meta":{"raw":{"variants":["Rigorous proof: weak stratification yields primitive equations","Primitive equations proven for weakly stratified flow","Weak-stratified limit: Navier-Stokes to ocean-atmosphere model","Hydrostatic-incompressible limit validated for weak stratification"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000307,"raw_usage":{"total_tokens":1590,"prompt_tokens":736,"completion_tokens":854,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":480,"completion_tokens_details":{"reasoning_tokens":788}},"tokens_in":480,"tokens_out":854,"duration_ms":8018,"temperature":1.0,"reasoning_tokens":788,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T07:12:46.330990+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate (1.1) with data satisfying (3.4) for ε = 10^{-2}, 10^{-3}, 10^{-4}; if the H^2 error between the horizontal velocity and a reference solution of (1.2) does not decay to zero (or the uniform bound (3.6) fails for some ε), the convergence theorem is false. As a complementary check, impose a small vertical acoustic perturbation and verify the dispersion relation (3.12): the vertical acoustic frequency must scale as |k_3|/ε², not ε^{-1}; a different scaling would falsify the three-wave separation.","supporting_citations":[],"review_version":1}