pith. sign in

arxiv: 2507.20488 · v4 · submitted 2025-07-28 · 🧮 math.AP · math-ph· math.MP· math.OC

Linear toroidal-inertial waves on a differentially rotating sphere with application to helioseismology: Modeling, forward and inverse problems

Tram Thi Ngoc Nguyen , Damien Fournier , Laurent Gizon , Thorsten Hohage This is my paper

Pith reviewed 2026-05-19 03:21 UTC · model grok-4.3

classification 🧮 math.AP math-phmath.MPmath.OC
keywords toroidal wavesinertial wavesdifferential rotationinverse problemhelioseismologyregularization methodswell-posednesssolar interior
0
0 comments X

The pith

A fourth-order scalar equation models linear toroidal waves on a differentially rotating sphere to reconstruct viscosity and rotation parameters from surface observations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper develops a framework for solar inertial waves under the idealized assumption that they are linear and purely toroidal. The stream function then satisfies a fourth-order scalar equation on the sphere. Well-posedness of the forward problem is shown for appropriate differential rotation. The inverse problem of jointly recovering viscosity and differential rotation from complete or partial surface data is solved using iterative regularization methods. Convergence is guaranteed by verifying the tangential cone condition, local uniqueness of parameters is proven, and numerical tests with the Nesterov-Landweber method show robustness to noise and incomplete data.

Core claim

The central claim is that linear toroidal-inertial waves on a differentially rotating sphere are described by a fourth-order scalar equation for the stream function. This forward model is well-posed under explicit conditions on the rotation. In the inverse setting, the tangential cone condition holds for the parameter-to-observation map, ensuring that iterative regularization methods converge and that viscosity together with differential rotation are locally uniquely determined from surface measurements.

What carries the argument

The fourth-order scalar equation for the stream function of the toroidal flow, which reduces the vector wave problem to a scalar one and underpins the analysis of both forward well-posedness and inverse identifiability.

If this is right

  • Explicit conditions on differential rotation guarantee existence and uniqueness of wave solutions.
  • Surface data suffice to reconstruct both viscosity and differential rotation simultaneously.
  • Iterative methods with regularization converge to the correct parameters even with noise.
  • Local uniqueness prevents multiple parameter sets from producing the same observations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This toroidal reduction might be relaxed in future work to include small poloidal flows for more complete solar models.
  • The identifiability proof could inspire similar analyses for other rotating fluid systems like planetary atmospheres.
  • Robustness to partial data implies that sparse observation networks may still enable interior parameter estimation in helioseismology.

Load-bearing premise

The waves must be purely toroidal and linear so that their dynamics reduce to the fourth-order scalar equation for the stream function.

What would settle it

Direct observation of significant poloidal velocity components in solar inertial waves or reconstruction algorithms failing to converge for synthetic data generated from the model itself would falsify the approach.

Figures

Figures reproduced from arXiv: 2507.20488 by Damien Fournier, Laurent Gizon, Thorsten Hohage, Tram Thi Ngoc Nguyen.

Figure 1
Figure 1. Figure 1: 1% data noise. Top to bottom: full and leaked data (filled area) at different [PITH_FULL_IMAGE:figures/full_fig_p022_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Leaked data. Top to bottom: 1%, 5%, 10%, 20% relative noise level. [PITH_FULL_IMAGE:figures/full_fig_p023_2.png] view at source ↗
read the original abstract

This paper develops a mathematical framework for interpreting observations of solar inertial waves in an idealized setting. Under the assumption of purely toroidal linear waves on the sphere, the stream function of the flow satisfies a fourth-order scalar equation. We prove well-posedness of wave solutions under explicit conditions on differential rotation. Moreover, we study the inverse problem of simultaneously reconstructing viscosity and differential rotation parameters from either complete or partial surface data. We establish convergence guarantee of iterative regularization methods by verifying the tangential cone condition, and prove local unique identifiability of the unknown parameters. Numerical experiments with Nesterov-Landweber iteration confirm reconstruction robustness across different observation strategies and noise levels.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper develops a mathematical framework for linear toroidal-inertial waves on a differentially rotating sphere under the assumption of purely toroidal flow. The stream function satisfies a fourth-order scalar PDE; well-posedness is proved under explicit conditions on the differential rotation profile. The inverse problem of recovering viscosity and differential rotation parameters from complete or partial surface observations of the stream function is studied, with convergence guarantees for Nesterov-Landweber iteration obtained by verifying the tangential cone condition and with a proof of local unique identifiability. Numerical experiments confirm reconstruction robustness across observation strategies and noise levels.

Significance. If the central claims hold, the work supplies a rigorous forward and inverse theory for an idealized model relevant to helioseismology, including parameter-free well-posedness conditions and a convergence proof for a nonlinear inverse problem. The combination of analytic identifiability results and reproducible numerical validation strengthens the case for using such models to interpret solar inertial-wave observations.

major comments (2)
  1. [Inverse-problem section (around the statement of the tangential cone condition and its verification)] The verification that the tangential cone condition holds for the forward map (defined via the fourth-order PDE) is stated to guarantee convergence of the Nesterov-Landweber iteration, yet the well-posedness result supplies existence only under explicit conditions on the rotation profile without deriving the necessary Lipschitz-type bounds on the Fréchet derivative that would keep the TCC constant finite and uniform as differential-rotation shear increases. This gap is load-bearing for the convergence claim.
  2. [Identifiability theorem and its proof] Local unique identifiability is asserted for the pair (viscosity, differential-rotation parameters) from surface data, but the proof sketch does not explicitly address whether the identifiability constant remains positive when the observation operator is restricted to partial surface data or when the background shear approaches the boundary of the well-posedness regime.
minor comments (2)
  1. [Model setup and notation] Notation for the differential-rotation profile and the viscosity coefficient should be introduced once and used consistently; several equations reuse symbols without redefinition.
  2. [Numerical experiments] Figure captions for the numerical reconstructions should state the exact noise level, the number of iterations, and whether the data are complete or partial surface observations.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address the two major comments point by point below and will revise the manuscript to strengthen the justification of the convergence result and the identifiability statement.

read point-by-point responses

  1. Referee: [Inverse-problem section (around the statement of the tangential cone condition and its verification)] The verification that the tangential cone condition holds for the forward map (defined via the fourth-order PDE) is stated to guarantee convergence of the Nesterov-Landweber iteration, yet the well-posedness result supplies existence only under explicit conditions on the rotation profile without deriving the necessary Lipschitz-type bounds on the Fréchet derivative that would keep the TCC constant finite and uniform as differential-rotation shear increases. This gap is load-bearing for the convergence claim.

    Authors: We acknowledge that the current presentation of the tangential cone condition verification would benefit from explicit Lipschitz bounds on the Fréchet derivative of the forward map. Under the well-posedness conditions already stated for the differential rotation profile, such bounds can be derived from the a priori estimates on the fourth-order operator; the resulting TCC constant remains finite and depends continuously on the shear. In the revised manuscript we will insert these estimates immediately after the well-posedness theorem, thereby making the convergence guarantee for Nesterov-Landweber iteration fully rigorous and uniform within the admissible parameter range. revision: yes

  2. Referee: [Identifiability theorem and its proof] Local unique identifiability is asserted for the pair (viscosity, differential-rotation parameters) from surface data, but the proof sketch does not explicitly address whether the identifiability constant remains positive when the observation operator is restricted to partial surface data or when the background shear approaches the boundary of the well-posedness regime.

    Authors: The local unique identifiability result holds for both complete and partial surface observations; the constant depends on the measure of the observed set and remains strictly positive whenever this set has positive surface measure. We will expand the proof sketch in the revision to make this dependence explicit. As the background shear approaches the boundary of the well-posedness regime the constant may deteriorate, and we will add a short remark together with supporting numerical diagnostics to clarify this limitation without altering the statement of the theorem. revision: partial

Circularity Check

0 steps flagged

No significant circularity; central claims rest on independent proofs of well-posedness and TCC verification

full rationale

The paper states an assumption of purely toroidal linear waves leading to a fourth-order scalar PDE, then proves well-posedness under explicit conditions on the differential rotation profile and separately verifies the tangential cone condition to guarantee convergence of Nesterov-Landweber iteration for the inverse problem. These steps are presented as new mathematical derivations and local identifiability results rather than reductions to fitted inputs, self-citations, or quantities defined by construction from the target parameters. No load-bearing step reduces by the paper's own equations to a prior self-citation or to a parameter fit renamed as a prediction; the derivation chain is self-contained against the stated PDE and observation operators.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The framework rests on the domain assumption of purely toroidal linear waves that reduces the vector problem to a scalar fourth-order equation; no free parameters are introduced in the forward model itself, and no new entities are postulated.

axioms (1)
  • domain assumption Purely toroidal linear waves on the sphere
    Invoked in the abstract to reduce the flow to a stream function satisfying a fourth-order scalar equation.

pith-pipeline@v0.9.0 · 5660 in / 1393 out tokens · 43930 ms · 2026-05-19T03:21:06.120125+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages · 2 internal anchors

  1. [1]

    Learning-informed parame- ter identification in nonlinear time-dependent PDEs

    C. Aarset, M. Holler, and T. T. N. Nguyen. “Learning-informed parame- ter identification in nonlinear time-dependent PDEs”. In: Applied Math- ematics and Optimization 88 (2023), 53 pp. doi: 10.1007/s00245-023- 10044-y

    work page doi:10.1007/s00245-023- 2023

  2. [2]

    R. A. Adams. Sobolev Spaces. New York: Academic Press, 1975

    work page 1975

  3. [3]

    Algebra properties for Sobolev spaces - applications to semilinear PDEs on manifolds

    N. Badr et al. “ Algebra properties for Sobolev spaces - applications to semilinear PDEs on manifolds”. In: JAMA 118 (2012), pp. 509–544. doi: 10.1007/s11854-012-0043-1

    work page doi:10.1007/s11854-012-0043-1 2012

  4. [4]

    Global Seismology of the Sun

    S. Basu. “Global seismology of the Sun”. In: Living Reviews in Solar Physics 13.1, 2 (Aug. 2016), p. 2. doi: 10.1007/s41116-016-0003-4 . arXiv: 1606.07071 [astro-ph.SR] . 29

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1007/s41116-016-0003-4 2016

  5. [6]

    Theory of solar oscillations in the inertial frequency range: Amplitudes of equatorial modes from a nonlinear rotating convection simulation

    Y. Bekki, R. H. Cameron, and L. Gizon. “Theory of solar oscillations in the inertial frequency range: Amplitudes of equatorial modes from a nonlinear rotating convection simulation”. In: Astronomy & Astrophysics 666, A135 (Oct. 2022), A135. doi: 10 . 1051 / 0004 - 6361 / 202244150. arXiv: 2208.11081 [astro-ph.SR]

    work page arXiv 2022

  6. [7]

    Nonlinear dynamics of Rossby waves in a western bound- ary current

    L. J. Campbell. “Nonlinear dynamics of Rossby waves in a western bound- ary current”. In: Advances in Fluid Mechanics VI 52 (2006), pp. 457–468. doi: 10.2495/AFM06045

    work page doi:10.2495/afm06045 2006

  7. [8]

    Science , year = 1996, month = may, volume = 272, pages =

    J. Christensen-Dalsgaard et al. “The Current State of Solar Modeling”. In: Science 272.5266 (1996). 10.1126/science.272.5266.1286, pp. 1286–

    work page doi:10.1126/science.272.5266.1286 1996

  8. [9]

    The non-radial oscillations of polytropic stars

    T. G. Cowling. “The non-radial oscillations of polytropic stars”. In: Monthly Notices of the Royal Astronomical Society 101 (Jan. 1941), p. 367. doi: 10.1093/mnras/101.8.367

    work page doi:10.1093/mnras/101.8.367 1941

  9. [10]

    Tangential Cone Condition for the Full Waveform Forward Operator in the Viscoelastic Regime: The Nonlocal Case

    M. Eller, R. Griesmaier, and A. Rieder. “Tangential Cone Condition for the Full Waveform Forward Operator in the Viscoelastic Regime: The Nonlocal Case”. In: SIAM Journal on Applied Mathematics 84.2 (2024), pp. 412–432. doi: 10.1137/23M1551845. eprint: https://doi.org/10. 1137/23M1551845. url: https://doi.org/10.1137/23M1551845

    work page doi:10.1137/23m1551845 2024

  10. [11]

    L. C. Evans. Partial Differential Equations . Graduate Studies in Mathe- matics 19. AMS, Providence, RI, 1998

    work page 1998

  11. [12]

    Viscous inertial modes on a dif- ferentially rotating sphere: Comparison with solar observations

    D. Fournier, L. Gizon, and L. Hyest. “Viscous inertial modes on a dif- ferentially rotating sphere: Comparison with solar observations”. In: As- tronomy & Astrophysics 664, A6 (2022), 16 pp. doi: 10 . 1051 / 0004 - 6361/202243473

    work page 2022

  12. [13]

    Meridional flow in the Sun’s convection zone is a single cell in each hemisphere

    L. Gizon et al. “Meridional flow in the Sun’s convection zone is a single cell in each hemisphere”. In: Science 368.6498 (June 2020), pp. 1469–

    work page 2020

  13. [14]

    doi: 10.1126/science.aaz7119

    work page doi:10.1126/science.aaz7119

  14. [15]

    Solar Inertial Modes

    L. Gizon et al. “Solar Inertial Modes”. In: Dynamics of Solar and Stellar Convection Zones and Atmospheres Proceedings IAU Symposium No. 365, A. V. Getling & L. L. Kitchatinov, eds . 2024

    work page 2024

  15. [16]

    , keywords =

    L. Gizon et al. “Solar inertial modes: Observations, identification, and diagnostic promise”. In: Astronomy & Astrophysics 652, L6 (2021), 22 pp. doi: 10.1051/0004-6361/202141462. 30

    work page doi:10.1051/0004-6361/202141462 2021

  16. [17]

    A convergence analysis of the Landweber iteration for nonlinear ill-posed problems

    M. Hanke, A. Neubauer, and O. Scherzer. “A convergence analysis of the Landweber iteration for nonlinear ill-posed problems”. In: Numer. Math. 72 (1995), pp. 21–37

    work page 1995

  17. [18]

    E. Hebey. Sobolev Spaces on Riemannian Manifolds. Springer Berlin, Hei- delberg, 1996. isbn: 978-3-540-69993-4. doi: 10.1007/BFb0092907

    work page doi:10.1007/bfb0092907 1996

  18. [20]

    Convergence analysis of a two-point gradi- ent method for nonlinear ill-posed problems

    S. Hubmer and R. Ramlau. “Convergence analysis of a two-point gradi- ent method for nonlinear ill-posed problems”. In: Inverse Problems 33.9 (2017), p. 095004. doi: 10.1088/1361-6420/aa7ac7

    work page doi:10.1088/1361-6420/aa7ac7 2017

  19. [21]

    Lam´ e Parameter Estimation from Static Displacement Field Measurements in the Framework of Nonlinear Inverse Problems

    S. Hubmer et al. “Lam´ e Parameter Estimation from Static Displacement Field Measurements in the Framework of Nonlinear Inverse Problems”. In: SIAM Journal on Imaging Sciences 11.2 (2018). doi: 10 . 1137 / 17M1154461

    work page 2018

  20. [22]

    Levenberg–Marquardt method for solving inverse problem of MRE based on the modified stationary Stokes system

    Y. Jiang, G. Nakamura, and K. Shirota. “Levenberg–Marquardt method for solving inverse problem of MRE based on the modified stationary Stokes system”. In: Inverse Problems 37.12, 125013 (2021)

    work page 2021

  21. [23]

    Kaltenbacher, A

    B. Kaltenbacher, A. Neubauer, and O. Scherzer. Iterative Regularization Methods for Nonlinear Problems . Radon Series on Computational and Applied Mathematics. Berlin, New York: de Gruyter, 2008

    work page 2008

  22. [24]

    Kaltenbacher, T

    B. Kaltenbacher, T. T. N. Nguyen, and O. Scherzer. The tangential cone condition for some coefficient identification model problems in parabolic PDEs. Springer, 2021, 121–163. doi: 10.1007/10.1007/978- 3- 030- 57784-1_5

    work page doi:10.1007/10.1007/978- 2021

  23. [25]

    IMA Journal of Numerical Analysis , volume =

    B. Kaltenbacher. “Convergence guarantees for coefficient reconstruction in PDEs from boundary measurements by variational and Newton-type methods via range invariance”. In: IMA Journal of Numerical Analysis drad044 (2023). doi: 10.1093/imanum/drad044

    work page doi:10.1093/imanum/drad044 2023

  24. [26]

    Convergence of the gradient method for ill-posed prob- lems

    S. Kindermann. “Convergence of the gradient method for ill-posed prob- lems”. In: Inverse Problems & Imaging 11.4 (2017), pp. 703–720

    work page 2017

  25. [27]

    On the tangential cone condition for electrical impedance tomography

    S. Kindermann. “On the tangential cone condition for electrical impedance tomography”. In: Electronic Transaction on Numerical Analysis57 (2022), pp. 17–34

    work page 2022

  26. [28]

    Quantitative passive imaging by iterative holography: The example of helioseismic holography

    B. M¨ uller et al. “Quantitative passive imaging by iterative holography: The example of helioseismic holography”. In: 40(4) (2024).doi: 10.1088/ 1361-6420/ad2b9a. 31

    work page 2024

  27. [29]

    J. C. N´ ed´ elec.Acoustic and Electromagnic Equations: Integral Respresen- tations for Harmonic Problems . New York: Springer, 2001

    work page 2001

  28. [30]

    On Nesterov acceleration for Landweber iteration of linear ill-posed problems

    A. Neubauer. “On Nesterov acceleration for Landweber iteration of linear ill-posed problems”. In: Journal of Inverse and Ill-posed Problems 25.3 (2017), pp. 381–390. doi: 10.1515/jiip-2016-0060

    work page doi:10.1515/jiip-2016-0060 2017

  29. [31]

    Sequential bi-level regularized inversion with applica- tion to hidden reaction law discovery

    T. T. N. Nguyen. “Sequential bi-level regularized inversion with applica- tion to hidden reaction law discovery”. In: Inverse Problems 41.6 (2025). doi: 10.1088/1361-6420/addf73

    work page doi:10.1088/1361-6420/addf73 2025

  30. [32]

    SpeedNet: dual- channel CNN-GRU network with adaptive noise fusion for robust vehicle dead reckoning in GNSS-denied environments,

    T. T. N. Nguyen. “Bi-level iterative regularization for inverse problems in nonlinear PDEs”. In: Inverse Problems 40.4 (2024), p. 045020. doi: 10.1088/1361- 6420/ad2905 . url: https://dx.doi.org/10.1088/ 1361-6420/ad2905

    work page doi:10.1088/1361- 2024

  31. [33]

    Astrophys.,670, A20

    J. Philidet and L. Gizon. “Interaction of solar inertial modes with turbu- lent convection. A 2D model for the excitation of linearly stable modes”. In: Astronomy & Astrophysics 673 (2023), 19 pp. doi: 10.1051/0004- 6361/202245666

    work page doi:10.1051/0004- 2023

  32. [34]

    Renardy and R

    M. Renardy and R. C. Rogers. An Introduction to Partial Differential Equations. Springer, 1992

    work page 1992

  33. [35]

    , keywords =

    P. H. Scherrer et al. “The Helioseismic and Magnetic Imager (HMI) Inves- tigation for the Solar Dynamics Observatory (SDO)”. In: Solar Physics 275.1-2 (Jan. 2012), pp. 207–227. doi: 10.1007/s11207-011-9834-2

    work page doi:10.1007/s11207-011-9834-2 2012

  34. [36]

    Gauss-Newton’s methods for solving linear inverse problems with neural network coders

    O. Scherzer, B. Hofmann, and Z. Nashed. “Gauss-Newton’s methods for solving linear inverse problems with neural network coders”. In: Sampl. Theory Signal Process. Data Anal. 21.25 (2023). doi: 10.1007/s43670- 023-00066-6

    work page doi:10.1007/s43670- 2023

  35. [37]

    M. E. Taylor. Partial Differential Equations I . Springer New York, NY,

    work page

  36. [38]

    doi: 10.1007/978-1-4419-7055-8

    isbn: 978-1-4614-2726-1. doi: 10.1007/978-1-4419-7055-8

    work page doi:10.1007/978-1-4419-7055-8

  37. [39]

    Disentangling the independently controllable factors of variation by interacting with the world

    M. J. Thompson et al. “The Internal Rota- tion of the Sun”. In: Annu. Rev. Astron. Astrophys. 41 (2003), pp. 599–643. doi: 10.1146/annurev. astro.41.011802.094848

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1146/annurev 2003

  38. [40]

    H. Triebel. Interpolation theory, function spaces, differential operators . Amsterdam ; New York : North-Holland Pub. Co., 1978

    work page 1978

  39. [41]

    Analysis of Finite Difference Discretization Schemes for Diffusion in Spheres with Variable Diffusivity

    A. Versypt and R. Braatz. “Analysis of Finite Difference Discretization Schemes for Diffusion in Spheres with Variable Diffusivity”. In: Comput Chem Eng. 71 (2014), pp. 241–252. doi: 10.1016/j.compchemeng.2014. 05.022

    work page doi:10.1016/j.compchemeng.2014 2014

  40. [42]

    Shear instability of differential rotation in stars

    M. Watson. “Shear instability of differential rotation in stars”. In: Geo- physical & Astrophysical Fluid Dynamics 16.1 (1981), pp. 285–298. doi: 10.1080/03091928008243663. 32

    work page doi:10.1080/03091928008243663 1981

  41. [43]

    E. Zeidler. Nonlinear Functional Analysis and Its Applications I: Fixed- point theorems. New York: Springer-Verlag, 1986, pp. xxi+897. doi: 10. 1007/978-1-4612-4838-5 . url: http://dx.doi.org/10.1007/978-1- 4612-4838-5. 33

    work page doi:10.1007/978-1- 1986