Effects of density stratification on Rossby waves in deep atmospheres
Pith reviewed 2026-05-18 00:37 UTC · model grok-4.3
The pith
Rossby waves trapped in the Sun's radiative interior maintain nearly constant vorticity through the entire convection zone.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In a deep, generally stratified solar model the vertical wave equation for Rossby waves possesses two propagation cavities. The radial vorticity eigenfunctions belonging to the radiative-interior cavity remain nearly constant with radius throughout the overlying convection zone.
What carries the argument
The vertical structure equation expressed in the Lagrangian pressure fluctuation δP, which yields singularity-free eigenfunctions and exposes the two-cavity structure under the β-plane approximation.
If this is right
- Rossby-wave propagation occurs in separate radial cavities rather than throughout the whole interior.
- Modes trapped in the radiative interior can still produce observable surface vorticity signals because their eigenfunctions do not decay in the convection zone.
- The Lagrangian-pressure formulation removes internal singularities that appeared in earlier formulations.
- General density stratification produces wave behavior that differs from both purely radiative and purely convective models.
Where Pith is reading between the lines
- Surface observations of Rossby waves could be used to infer properties of the radiative interior without needing to resolve the convection zone directly.
- The same vertical-equation approach might be applied to other stars that possess both a radiative core and a convective envelope.
- If the near-constant vorticity holds, it suggests a dynamical decoupling between the radiative interior and the convection zone for these particular waves.
Load-bearing premise
The beta-plane approximation remains valid over the full radial extent of the solar interior when the vertical structure equation is derived for a generally stratified fluid.
What would settle it
A direct measurement, from helioseismic or surface-flow data, showing that the radial vorticity profile of a low-frequency Rossby wave varies significantly with depth inside the convection zone.
Figures
read the original abstract
Though Rossby waves have been observed on the Sun, their radial eigenfunctions remain a mystery. The prior theoretical work either considers quasi-2D systems, which do not apply to the solar interior, or only considers fully radiative or fully convective atmospheres. This project calculates the radial eigenfunctions for Rossby waves in a deep atmosphere for a general stratification. Here, we use the $\beta$-plane approximation to derive a vertical equation in terms of the Lagrangian pressure fluctuation $\delta P$, and we then calculate radial eigenfunctions for Rossby waves in a standard solar model, Model S. We find that working in the Lagrangian pressure fluctuation results in cleaner wave equations that lack internal singularities that have been encountered in prior work. The resulting radial wave equation makes it abundantly clear that there are two wave cavities in the solar interior, one in the radiative interior and another in the convection zone. Surprisingly, our calculated radial vorticity eigenfunctions for the radiative interior modes are nearly constant throughout the convection zone, raising the possibility that they may be observable at the solar surface.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives a vertical structure equation for Rossby waves in a generally stratified deep atmosphere by adopting the β-plane approximation and working with the Lagrangian pressure fluctuation δP. Applied to the standard solar model (Model S), the resulting eigenproblem reveals two distinct wave cavities (radiative interior and convection zone). The central result is that radial vorticity eigenfunctions for modes trapped in the radiative interior remain nearly constant throughout the convection zone, raising the prospect of surface observability.
Significance. If the derivation and eigenfunctions hold, the work supplies the first explicit radial structure for Rossby waves across a realistic solar stratification, moving beyond quasi-2D or uniform-atmosphere models. The Lagrangian-pressure formulation is credited for producing singularity-free equations, and the reported near-constant vorticity in the convection zone constitutes a falsifiable prediction that could link interior modes to surface observations.
major comments (2)
- [derivation of the vertical equation (abstract and §2)] The central claim that radiative-interior modes have nearly constant vorticity eigenfunctions across the convection zone rests on the vertical structure equation derived under the β-plane approximation extended over the full radial domain (radiative interior through convection zone). This approximation replaces spherical geometry with a local Cartesian frame and assumes linear latitudinal variation of the Coriolis parameter; the manuscript does not quantify the size of neglected curvature or radial-vertical reorientation terms in a strongly stratified medium, which directly affects whether the reported eigenfunction constancy is robust.
- [vertical wave equation for δP] The statement that the Lagrangian-pressure formulation removes internal singularities is load-bearing for the cleanliness of the wave equation and the subsequent eigenfunction calculation. Without an explicit side-by-side comparison (e.g., the singular terms that appear in the Eulerian formulation versus their absence here) or a demonstration that the new equation remains regular at all turning points in Model S, the improvement cannot be verified from the given derivation.
minor comments (2)
- [abstract and §3] The reference to 'Model S' should include the full citation (Christensen-Dalsgaard et al.) at first mention for reproducibility.
- [figures] Figure captions for the eigenfunction plots should explicitly state the normalization and the radial range shown (e.g., 0.2–1.0 R⊙) to allow direct comparison with the text claims.
Simulated Author's Rebuttal
We thank the referee for their thoughtful and constructive report. We address the two major comments below, clarifying the scope of the β-plane approximation and the advantages of the Lagrangian formulation while agreeing to strengthen the manuscript with additional discussion and comparisons where feasible.
read point-by-point responses
-
Referee: [derivation of the vertical equation (abstract and §2)] The central claim that radiative-interior modes have nearly constant vorticity eigenfunctions across the convection zone rests on the vertical structure equation derived under the β-plane approximation extended over the full radial domain (radiative interior through convection zone). This approximation replaces spherical geometry with a local Cartesian frame and assumes linear latitudinal variation of the Coriolis parameter; the manuscript does not quantify the size of neglected curvature or radial-vertical reorientation terms in a strongly stratified medium, which directly affects whether the reported eigenfunction constancy is robust.
Authors: We agree that a quantitative assessment of the neglected curvature and radial-vertical terms would strengthen the justification for extending the β-plane over the full domain. Our choice follows standard practice in deep-atmosphere Rossby-wave studies to isolate stratification effects while retaining the essential latitudinal variation of the Coriolis parameter. The resulting eigenfunctions are robust within the adopted framework, but we will add an explicit order-of-magnitude estimate of the omitted terms (using Model S profiles) in a new subsection of §2 in the revised manuscript. revision: yes
-
Referee: [vertical wave equation for δP] The statement that the Lagrangian-pressure formulation removes internal singularities is load-bearing for the cleanliness of the wave equation and the subsequent eigenfunction calculation. Without an explicit side-by-side comparison (e.g., the singular terms that appear in the Eulerian formulation versus their absence here) or a demonstration that the new equation remains regular at all turning points in Model S, the improvement cannot be verified from the given derivation.
Authors: The Lagrangian pressure δP was chosen precisely because it yields a second-order equation free of the density-gradient singularities that appear when the Eulerian pressure perturbation is used. We will insert a short appendix or expanded paragraph in §2 that writes the corresponding Eulerian equation, highlights the singular coefficients proportional to d ln ρ/dr, and confirms that those terms are absent in our δP formulation. Regularity at the turning points in Model S follows directly from the absence of those coefficients and will be noted explicitly. revision: yes
Circularity Check
No significant circularity; derivation uses external Model S and standard beta-plane equations
full rationale
The paper derives the vertical structure equation for Lagrangian pressure fluctuation δP by applying the standard β-plane approximation to the linearized equations of motion in a stratified atmosphere, then numerically integrates the resulting eigenvalue problem against the external standard solar model (Model S). No parameters are fitted to the reported eigenfunctions or their near-constancy property; the eigenfunctions are computed outputs. No self-citations, ansatzes smuggled via prior work, or definitional loops appear in the derivation chain. The choice of variable (δP) and the identification of two wave cavities follow directly from the approximated equations without reducing to the target result by construction. The beta-plane validity across the full domain is an assumption whose correctness can be assessed externally, but it does not create circularity within the paper's own logic.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The beta-plane approximation is applicable for deriving the vertical Rossby-wave equation in a deep, radially stratified atmosphere.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We use the β-plane approximation to derive a vertical equation in terms of the Lagrangian pressure fluctuation δP... d²δP/dz² + (1/Hρ) dδP/dz + N²/(g h) δP = 0
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The resulting radial wave equation makes it abundantly clear that there are two wave cavities... radiative interior and convection zone
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Albekioni, M., Zaqarashvili, T. V ., & Kukhianidze, V . 2023, A&A, 671, A91, doi: 10.1051/0004-6361/202243985
-
[2]
Alshehhi, R., Hanson, C. S., Gizon, L., & Hanasoge, S. 2019, A&A, 622, A124, doi: 10.1051/0004-6361/201834237
-
[3]
Bekki, Y ., Cameron, R. H., & Gizon, L. 2022a, A&A, 662, A16, doi: 10.1051/0004-6361/202243164 —. 2022b, A&A, 666, A135, doi: 10.1051/0004-6361/202244150
-
[4]
2002, Reviews of Modern Physics, 74, 1073, doi: 10.1103/RevModPhys.74.1073
Christensen-Dalsgaard, J. 2002, Reviews of Modern Physics, 74, 1073, doi: 10.1103/RevModPhys.74.1073
-
[5]
Christensen-Dalsgaard, J., Dappen, W., Ajukov, S. V ., et al. 1996a, Science, 272, 1286, doi: 10.1126/science.272.5266.1286 —. 1996b, Science, 272, 1286, doi: 10.1126/science.272.5266.1286
-
[6]
Cox, J. P. 1983, Theory of stellar pulsations. (Princeton University Press)
work page 1983
-
[7]
Damiani, C., Cameron, R. H., Birch, A. C., & Gizon, L. 2020, A&A, 637, A65, doi: 10.1051/0004-6361/201936251
-
[8]
Gizon, L., Cameron, R. H., Bekki, Y ., et al. 2021, A&A, 652, L6, doi: 10.1051/0004-6361/202141462
-
[9]
Gough, D. O. 1993, in Astrophysical Fluid Dynamics - Les Houches 1987, ed. J. P. Zahn & J. Zinn-Justin, 399–560
work page 1993
-
[10]
Hanasoge, S., & Mandal, K. 2019, ApJL, 871, L32, doi: 10.3847/2041-8213/aaff60
-
[11]
Hanson, C. S., Gizon, L., & Liang, Z.-C. 2020, A&A, 635, A109, doi: 10.1051/0004-6361/201937321
-
[12]
Hanson, C. S., & Hanasoge, S. 2024, Physics of Fluids, 36, 086626, doi: 10.1063/5.0216403
-
[13]
Hathaway, D. H., & Upton, L. A. 2021, ApJ, 908, 160, doi: 10.3847/1538-4357/abcbfa
-
[14]
1940, Transactions, American Geophysical Union, 21, 262, doi: 10.1029/TR021i002p00262
Haurwitz, B. 1940, Transactions, American Geophysical Union, 21, 262, doi: 10.1029/TR021i002p00262
-
[15]
Hindman, B. W., & Jain, R. 2022, ApJ, 932, 68, doi: 10.3847/1538-4357/ac6d64
-
[16]
Hindman, B. W., & Julien, K. 2024, ApJ, 960, 64, doi: 10.3847/1538-4357/ad0967
-
[17]
1956, Ordinary Differential Equations (USA: Dover
Ince, E. 1956, Ordinary Differential Equations (USA: Dover
work page 1956
-
[18]
1949, Memoires of the Societe Royale des Sciences de Liege, 9, 3 9
Ledoux, P. 1949, Memoires of the Societe Royale des Sciences de Liege, 9, 3 9
work page 1949
-
[19]
Liang, Z.-C., Gizon, L., Birch, A. C., & Duvall, T. L. 2019, A&A, 626, A3, doi: 10.1051/0004-6361/201834849
-
[20]
Longuet-Higgins, M. S. 1968, Philosophical Transactions of the Royal Society of London Series A, 262, 511, doi: 10.1098/rsta.1968.0003 L¨optien, B., Gizon, L., Birch, A. C., et al. 2018, Nature Astronomy, 2, 568, doi: 10.1038/s41550-018-0460-x
-
[21]
Mandal, K., & Hanasoge, S. M. 2024, ApJ, 967, 46, doi: 10.3847/1538-4357/ad391b
-
[22]
Papaloizou, J., & Pringle, J. E. 1978, MNRAS, 182, 423, doi: 10.1093/mnras/182.3.423
-
[23]
1987, Geophysical Fluid Dynamics (Springer-Verlag New York Inc)
Pedlosky, J. 1987, Geophysical Fluid Dynamics (Springer-Verlag New York Inc)
work page 1987
- [24]
-
[25]
Proxauf, B., Gizon, L., L¨optien, B., et al. 2020, A&A, 634, A44, doi: 10.1051/0004-6361/201937007
-
[26]
1939, Journal of Marine Research, 2, 38
Rossby, C.-G. 1939, Journal of Marine Research, 2, 38
work page 1939
-
[27]
Vallis, G. K. 2017, Atmospheric and Oceanic Fluid Dynamics: Fundamentals and Large-Scale Circulation (Cambridge University Press), doi: 10.1017/9781107588417
-
[28]
Waidele, M., & Zhao, J. 2023, ApJL, 954, L26, doi: 10.3847/2041-8213/acefd0
-
[29]
Wolff, C. L., & Blizard, J. B. 1986, SoPh, 105, 1, doi: 10.1007/BF00156371
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.