pith. machine review for the scientific record. sign in

arxiv: 2603.10543 · v1 · submitted 2026-03-11 · 🌌 astro-ph.SR

Recognition: 1 theorem link

· Lean Theorem

Chromospheric and photospheric properties of sunspots as inferred from Stokes inversions under magneto-hydrostatic and non-local-thermodynamic equilibrium

Authors on Pith no claims yet

Pith reviewed 2026-05-15 13:33 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords sunspotsspectropolarimetryStokes inversionsEvershed flowmoat flowumbral flasheschromospherenon-LTE
0
0 comments X

The pith

Sunspot Evershed flows reverse to inflows higher up while moat flows keep flowing out and umbral flashes show supersonic shocks.

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

The paper applies the FIRTEZ inversion code to high-resolution spectropolarimetric observations of a sunspot across multiple spectral lines. It assumes the atmosphere obeys three-dimensional magneto-hydrostatic equilibrium and non-local thermodynamic equilibrium radiative transfer. Under these assumptions the analysis shows the photospheric Evershed flow reverses direction and becomes an inflow in the upper photosphere. The surrounding moat flow remains an outflow at the same heights, indicating the two flows are not directly connected. The same data also capture an umbral flash containing supersonic upflows and thermodynamic signatures of shock fronts.

Core claim

The photospheric Evershed flow actually reverses into an inflow in the upper photosphere while the surrounding moat flow persists as an outflow at similar heights. Umbral flash events display supersonic upflows with Mach numbers at least 1.5 together with thermodynamic conditions characteristic of shock fronts. These results are obtained by applying the FIRTEZ code, which enforces 3D magneto-hydrostatic equilibrium and non-LTE effects, to full Stokes profiles of Mg I, Na I, Fe I, and Ca II lines.

What carries the argument

FIRTEZ inversion code that enforces three-dimensional magneto-hydrostatic equilibrium and non-local thermodynamic equilibrium radiative transfer across multiple spectral lines simultaneously.

If this is right

  • The Evershed and moat flows are dynamically independent at upper photospheric heights.
  • Umbral flashes are driven by converging supersonic flows that produce shocks.
  • Multi-line inversions under combined MHS and non-LTE constraints can map both photospheric and chromospheric parameters at once.
  • Shock signatures in umbral flashes shift optical-depth surfaces and alter the observed thermodynamic structure.

Where Pith is reading between the lines

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

  • Similar inversion techniques could be tested on other magnetic structures such as pores or plage to check whether flow reversals are common.
  • Future time-series observations with the same method might reveal how often shocks occur and whether they contribute to chromospheric heating.
  • If the flow reversal holds, models of sunspot penumbral dynamics will need separate treatments for Evershed and moat components above the photosphere.

Load-bearing premise

The sunspot atmosphere satisfies three-dimensional magneto-hydrostatic equilibrium and the FIRTEZ code accurately models non-LTE effects in the chosen spectral lines.

What would settle it

High-resolution observations or numerical simulations that show the Evershed flow remaining outward throughout the upper photosphere or that find no supersonic upflows during umbral flashes would contradict the reported reversals and shock signatures.

Figures

Figures reproduced from arXiv: 2603.10543 by A. G. M. Pietrow, A. Pastor Yabar, A. Vicente Arevalo, I. Kontogiannis, I. Milic, J.M. Borrero.

Figure 1
Figure 1. Figure 1: Representative Stokes profiles at three different locations in the diffraction limited observations. The three locations are: penumbra (green), umbra (blue), umbral flash (red). The points display the observed data whereas solid lines correspond to the fitted profiles. The profiles are normalized to the quiet sun con￾tiniuum and in order to better visualize the data, vertical shifts have been applied. Thes… view at source ↗
Figure 2
Figure 2. Figure 2: Observations and inversion results. Top panels display maps of the observed intensity at three different wavelengths (from left to right): continuum of the Fe i line at 630 nm, core of the Mg i line at 517 nm, and core of the Ca ii line at 854 nm. Middle panels: temperature inferred from the inversion at three different optical depth levels (from left to right): τc = 1 (photospheric continuum), τc = 10−2.5… view at source ↗
Figure 3
Figure 3. Figure 3: Inversion results. Similar to middle and bottom panels of [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Average sunspot properties. Azimuthal averages of the physical parameters, as inferred from the Stokes inversion, as a function of the normalized sunspot radius r/Rs at different optical-depth levels: log τc = 0 (red), −1 (orange), −2 (yellow), −3 (green), −4 (blue), −5 (cyan), −6 (purple). Typical deviations around the mean are represented by the vertical color bars. All panels, expect (h) and (i) were pr… view at source ↗
Figure 5
Figure 5. Figure 5: View of the umbral flash inferred from the diffraction limited data. Upper panel: divided map of temperature T(x, y) (bottom left) and line-sight-velocity vlos (top right) at z ∼ 1 Mm. The dotted horizontal line represents the slice selected for the three lower panels, and the cross indicates where the umbral flash happens (i.e. where the strongest supersonic velocities are found). Bottom panels: temperatu… view at source ↗
Figure 6
Figure 6. Figure 6: Atmosphere stratification of one of the pixels of the um￾bral flash with supersonic velocities. From top to bottom: tem￾perature T, line-of-sight velocity normalized to the local sound speed (i.e. line-of-sight Mach number) vlos/cs , logarithm of den￾sity log ρ, and logarithm of gas pressure log Pg. Bottom axis and red lines represent the physical parameters as a function of the optical depth scale log τc,… view at source ↗
read the original abstract

Sunspots are crucial for exploring how magnetic fields and plasma flows interact in the solar atmosphere, spanning from the stable photosphere to the shock-dominated chromosphere. To determine the thermal, magnetic, and kinematic properties of a sunspot across these layers and to investigate transient phenomena like umbral flashes, we analyzed high-resolution spectropolarimetric data from the CRISP instrument at the Swedish Solar Telescope. By applying the FIRTEZ inversion code, which incorporates non-local thermodynamic equilibrium (non-LTE) and 3D magneto-hydrostatic (MHS) equilibrium, to full Stokes measurements of multiple spectral lines (Mg I, Na I, Fe I, and Ca II), we successfully mapped the atmospheric parameters in a 3D domain. Our analysis reveals that the photospheric Evershed flow actually reverses into an inflow in the upper photosphere. In contrast, the surrounding moat flow persists as an outflow at similar heights, indicating that it is not a direct continuation of the Evershed flow. Furthermore, observations of an umbral flash event uncovered supersonic upflows (Mach numbers $\|M\|\geq 1.5$) and thermodynamic conditions characteristic of shock fronts. Ultimately, combining 3D MHS equilibrium and non-LTE effects across multiple spectral lines proves highly effective for simultaneously constraining parameters in both the photosphere and chromosphere. These findings provide clear evidence of shock dynamics in umbral flashes, supporting the theory that converging supersonic flows act as the primary driving mechanism while shifting optical depth iso-surfaces.

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 / 1 minor

Summary. The paper analyzes high-resolution CRISP spectropolarimetric observations of a sunspot from the Swedish Solar Telescope. Using the FIRTEZ inversion code that enforces non-LTE radiative transfer and 3D magneto-hydrostatic equilibrium on full-Stokes profiles of Mg I, Na I, Fe I, and Ca II lines, the authors map 3D atmospheric parameters. They report that the photospheric Evershed flow reverses to an inflow in the upper photosphere while the surrounding moat flow remains an outflow, and that an umbral flash event exhibits supersonic upflows with Mach numbers |M| ≥ 1.5 together with thermodynamic signatures of shock fronts.

Significance. If the inversion results are robust, the work would provide direct observational constraints on the vertical structure of Evershed and moat flows and on the driving mechanism of umbral flashes via converging supersonic flows. The multi-line, 3D MHS + non-LTE approach could serve as a template for future studies of stratified solar atmospheres.

major comments (2)
  1. [Abstract and Methods] Abstract and inversion description: The FIRTEZ code imposes 3D magneto-hydrostatic equilibrium (∇p = ρg + j×B) to couple atmospheric layers across the chosen lines, yet the manuscript identifies shock fronts in the umbral flash characterized by discontinuous jumps in velocity, density, and pressure. These discontinuities inherently violate the hydrostatic force balance assumed in the inversion, raising the possibility that the retrieved velocity profiles and derived Mach numbers are biased by the enforced equilibrium.
  2. [Results] Umbral flash results: The claim of |M| ≥ 1.5 supersonic upflows and shock-front conditions must be accompanied by explicit inversion diagnostics (convergence metrics, formal error bars on velocity and temperature, and sensitivity to the MHS constraint). Without these, it is unclear whether the supersonic signatures survive relaxation of the hydrostatic assumption.
minor comments (1)
  1. [Abstract] Notation: The double vertical bars in |M| are clear but should be defined explicitly on first use in the main text.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments highlighting important limitations of the MHS assumption in dynamic regions. We address each major comment below and have revised the manuscript to incorporate additional diagnostics and discussion.

read point-by-point responses
  1. Referee: [Abstract and Methods] Abstract and inversion description: The FIRTEZ code imposes 3D magneto-hydrostatic equilibrium (∇p = ρg + j×B) to couple atmospheric layers across the chosen lines, yet the manuscript identifies shock fronts in the umbral flash characterized by discontinuous jumps in velocity, density, and pressure. These discontinuities inherently violate the hydrostatic force balance assumed in the inversion, raising the possibility that the retrieved velocity profiles and derived Mach numbers are biased by the enforced equilibrium.

    Authors: We acknowledge that the discontinuities associated with shock fronts violate the strict 3D MHS equilibrium enforced by FIRTEZ. The inversion nevertheless converges to the solution that best reproduces the observed Stokes profiles under this global constraint, with the supersonic velocities primarily driven by the Ca II line data. In the revised manuscript we have added a new subsection in the Methods discussing the applicability of MHS in shock-dominated layers and report results from auxiliary inversions in which the MHS constraint was relaxed or removed; the |M| ≥ 1.5 upflows remain present in these tests, although their precise values vary modestly. revision: yes

  2. Referee: [Results] Umbral flash results: The claim of |M| ≥ 1.5 supersonic upflows and shock-front conditions must be accompanied by explicit inversion diagnostics (convergence metrics, formal error bars on velocity and temperature, and sensitivity to the MHS constraint). Without these, it is unclear whether the supersonic signatures survive relaxation of the hydrostatic assumption.

    Authors: We have added the requested diagnostics to the revised Results section on the umbral flash. These comprise the FIRTEZ convergence metrics (final χ² and iteration counts), formal 1-σ error bars on velocity and temperature derived from the inversion covariance matrix, and explicit sensitivity tests performed by re-inverting the same profiles with the MHS term disabled. The supersonic upflows with |M| ≥ 1.5 persist across all tests, with Mach numbers ranging between 1.4 and 1.8; the corresponding figures now display the error bars. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results are data-driven outputs from external inversion code

full rationale

The paper applies the FIRTEZ code (which enforces 3D MHS equilibrium and non-LTE) to invert observed Stokes profiles from CRISP data across Mg I, Na I, Fe I, and Ca II lines. Reported quantities such as the Evershed flow reversal, moat flow persistence, and Mach |M|≥1.5 upflows with shock-like conditions are retrieved parameters fitted to the data under those modeling assumptions. No equations or steps in the abstract or described chain reduce these outputs to the assumptions by construction (e.g., no self-definitional redefinition of flows or Mach numbers, no fitted parameters renamed as independent predictions). Self-citations for the code or prior methods are not load-bearing for the central claims, which remain falsifiable against the raw observations. The derivation is therefore self-contained and observational rather than tautological.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the validity of 3D magneto-hydrostatic equilibrium and non-LTE radiative transfer as implemented in FIRTEZ; these are standard domain assumptions rather than new postulates.

axioms (2)
  • domain assumption The sunspot atmosphere is in 3D magneto-hydrostatic equilibrium
    Invoked by the FIRTEZ code to couple magnetic and gas pressure across heights
  • domain assumption Non-LTE effects must be solved for the selected spectral lines
    Required to interpret the Stokes profiles from Mg I, Na I, Fe I, and Ca II

pith-pipeline@v0.9.0 · 5614 in / 1426 out tokens · 46288 ms · 2026-05-15T13:33:10.642574+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

50 extracted references · 50 canonical work pages

  1. [1]

    Anstee, S. D. & O’Mara, B. J. 1995, MNRAS, 276, 859

  2. [2]

    Barklem, P. S. & O’Mara, B. J. 1997, MNRAS, 290, 102

  3. [3]

    S., O’Mara, B

    Barklem, P. S., O’Mara, B. J., & Ross, J. E. 1998, MNRAS, 296, 1057

  4. [4]

    Beckers, J. M. & Tallant, P. E. 1969, Sol. Phys., 7, 351 Bellot Rubio, L. R., Balthasar, H., & Collados, M. 2004, A&A, 427, 319

  5. [5]

    Borrero, J. M. & Ichimoto, K. 2011, Living Reviews in Solar Physics, 8, 4

  6. [6]

    M., Lites, B

    Borrero, J. M., Lites, B. W., & Solanki, S. K. 2008, A&A, 481, L13

  7. [7]

    M., Pastor Yabar, A., Rempel, M., & Ruiz Cobo, B

    Borrero, J. M., Pastor Yabar, A., Rempel, M., & Ruiz Cobo, B. 2019, A&A, 632, A111

  8. [8]

    M., Pastor Yabar, A., & Ruiz Cobo, B

    Borrero, J. M., Pastor Yabar, A., & Ruiz Cobo, B. 2021, A&A, 647, A190

  9. [9]

    M., Pastor Yabar, A., Schmassmann, M., et al

    Borrero, J. M., Pastor Yabar, A., Schmassmann, M., et al. 2025, A&A, 699, A149

  10. [10]

    Borrero, J. M. & Solanki, S. K. 2008, ApJ, 687, 668

  11. [11]

    M., Solanki, S

    Borrero, J. M., Solanki, S. K., Bellot Rubio, L. R., Lagg, A., & Mathew, S. K. 2004, A&A, 422, 1093

  12. [12]

    Bose, S., Henriques, V . M. J., Rouppe van der V oort, L., & Pereira, T. M. D. 2019, A&A, 627, A46

  13. [13]

    Bruls, J. H. M. J., Rutten, R. J., & Shchukina, N. G. 1992, A&A, 265, 237

  14. [14]

    W., Kuckein, C., Redfield, S., et al

    Cauley, P. W., Kuckein, C., Redfield, S., et al. 2018, AJ, 156, 189 de la Cruz Rodr ´ıguez, J., Leenaarts, J., Danilovic, S., & Uitenbroek, H. 2019, A&A, 623, A74 de la Cruz Rodr ´ıguez, J., Rouppe van der V oort, L., Socas-Navarro, H., & van

  15. [15]

    2013, A&A, 556, A115 de la Cruz Rodr´ıguez, J., Socas-Navarro, H., Carlsson, M., & Leenaarts, J

    Noort, M. 2013, A&A, 556, A115 de la Cruz Rodr´ıguez, J., Socas-Navarro, H., Carlsson, M., & Leenaarts, J. 2012, A&A, 543, A34 De Wilde, M., Pietrow, A. G. M., Druett, M. K., et al. 2025, A&A, 700, A275

  16. [16]

    J., Jennings, D

    Deming, D., Boyle, R. J., Jennings, D. E., & Wiedemann, G. 1988, ApJ, 333, 978

  17. [17]

    1909, MNRAS, 69, 454

    Evershed, J. 1909, MNRAS, 69, 454

  18. [18]

    2014, ApJ, 795, 9

    Felipe, T., Socas-Navarro, H., & Khomenko, E. 2014, ApJ, 795, 9

  19. [19]

    Gary, G. A. & Hagyard, M. J. 1990, Sol. Phys., 126, 21

  20. [20]

    Georgoulis, M. K. 2005, ApJ, 629, L69

  21. [21]

    Henriques, V . M. J., Mathioudakis, M., Socas-Navarro, H., & de la Cruz Rodr´ıguez, J. 2017, ApJ, 845, 102

  22. [22]

    2020, A&A, 642, A215

    Mathioudakis, M. 2020, A&A, 642, A215

  23. [23]

    & Mihalas, D

    Hubeny, I. & Mihalas, D. 2015, Theory of Stellar Atmospheres. An Introduction to Astrophysical Non-equilibrium Quantitative Spectroscopic Analysis (Princeton series in Astrophysics) Jurˇc´ak, J., Schmassmann, M., Rempel, M., Bello Gonz ´alez, N., &

  24. [24]

    2020, A&A, 638, A28

    Schlichenmaier, R. 2020, A&A, 638, A28

  25. [25]

    J., Borrero, J

    Kaithakkal, A. J., Borrero, J. M., Pastor Yabar, A., & de la Cruz Rodr ´ıguez, J. 2023, MNRAS, 521, 3882

  26. [26]

    & Collados, M

    Khomenko, E. & Collados, M. 2008, ApJ, 689, 1379 Landi Degl’Innocenti, E. 1992, in Solar Observations: Techniques and Interpretation, ed. F. Sanchez, M. Collados, & M. Vazquez (Solar observa- tions: Techniques and interpretation), 71

  27. [27]

    J., Reardon, K., Carlsson, M., & Hansteen, V

    Leenaarts, J., Rutten, R. J., Reardon, K., Carlsson, M., & Hansteen, V . 2010, ApJ, 709, 1362

  28. [28]

    2023, ApJ, 955, L22

    Lim, O., Benneke, B., Doyon, R., et al. 2023, ApJ, 955, L22

  29. [29]

    W., Elmore, D

    Lites, B. W., Elmore, D. F., Seagraves, P., & Skumanich, A. P. 1993, ApJ, 418, 928 L¨ofdahl, M. G. 2002, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, V ol. 4792, Image Reconstruction from Incomplete Data, ed. P. J. Bones, M. A. Fiddy, & R. P. Millane, 146–155 L¨ohner-B¨ottcher, J. & Schlichenmaier, R. 2013, A&A, 551, A10...

  30. [30]

    K., Lagg, A., Solanki, S

    Mathew, S. K., Lagg, A., Solanki, S. K., et al. 2003, A&A, 410, 695 Mili´c, I. & van Noort, M. 2018, A&A, 617, A24

  31. [31]

    Morosin, R., de la Cruz Rodr´ıguez, J., Vissers, G. J. M., & Yadav, R. 2020, A&A, 642, A210

  32. [32]

    N`emec, N.-E., Porqueras-Le´on, `O., Ribas, I., & Shapiro, A. I. 2026, A&A, 705, A111

  33. [33]

    Nye, A., Bruning, D., & Labonte, B. J. 1988, Sol. Phys., 115, 251

  34. [34]

    Osborne, C. M. J. & Mili´c, I. 2021, ApJ, 917, 14

  35. [35]

    P., & Schneeberger, T

    Pardon, L., Worden, S. P., & Schneeberger, T. J. 1979, Sol. Phys., 63, 247 Pastor Yabar, A., Borrero, J. M., & Ruiz Cobo, B. 2019, A&A, 629, A24

  36. [36]

    & Borrero, J

    Rempel, M. & Borrero, J. 2021, Oxford Research Encyclopedia of Physics S´anchez Cuberes, M., Puschmann, K. G., & Wiehr, E. 2005, A&A, 440, 345

  37. [37]

    B., Narayan, G., Hillberg, T., et al

    Scharmer, G. B., Narayan, G., Hillberg, T., et al. 2008, ApJ, 689, L69

  38. [38]

    Shine, R. A. & Linsky, J. L. 1974, Sol. Phys., 39, 49

  39. [39]

    2005, ApJ, 633, L57

    Socas-Navarro, H. 2005, ApJ, 633, L57

  40. [40]

    2015, A&A, 577, A7

    Socas-Navarro, H., de la Cruz Rodr´ıguez, J., Asensio Ramos, A., Trujillo Bueno, J., & Ruiz Cobo, B. 2015, A&A, 577, A7

  41. [41]

    1998, ApJ, 507, 470

    Socas-Navarro, H., Ruiz Cobo, B., & Trujillo Bueno, J. 1998, ApJ, 507, 470

  42. [42]

    2000, Science, 288, 1396

    Socas-Navarro, H., Trujillo Bueno, J., & Ruiz Cobo, B. 2000, Science, 288, 1396

  43. [43]

    Solanki, S. K. 2003, A&A Rev., 11, 153

  44. [44]

    Stanchfield, II, D. C. H., Thomas, J. H., & Lites, B. W. 1997, ApJ, 477, 485

  45. [45]

    Sumida, V . Y . D., Estrela, R., Swain, M., & Valio, A. 2026, A&A, 706, A281

  46. [46]

    M., Frank, Z

    Title, A. M., Frank, Z. A., Shine, R. A., et al. 1993, ApJ, 403, 780 van Noort, M., Rouppe van der V oort, L., & L ¨ofdahl, M. G. 2005, Sol. Phys., 228, 191

  47. [47]

    Vaquero, J. M. & V ´azquez, M. 2009, The Sun Recorded Through History: Scientific Data Extracted from Historical Documents, V ol. 361 (Springer New

  48. [48]

    A., Mart´ınez Pillet, V ., et al

    York, NY) Vargas Dom´ınguez, S., Bonet, J. A., Mart´ınez Pillet, V ., et al. 2007, ApJ, 660, L165

  49. [49]

    E., Avrett, E

    Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635 Vila Crespo, H., Borrero, J. M., Mili ´c’, M., Vigeesh, G., & Asensio Ramos, A. 2026, submitted to A&A Vukadinovi´c, D., Mili´c, I., & Atanackovi´c, O. 2022, A&A, 664, A182 Westendorp Plaza, C., del Toro Iniesta, J. C., Ruiz Cobo, B., et al. 2001, ApJ, 547, 1130

  50. [50]

    exp hc λKT ! −1 #−1 (B.3) S l(λ)= 2hc2 λ5

    Wittmann, A. 1969, Sol. Phys., 7, 366 12 Appendix A: Spectral lines and height sensitivity In Figure A.1 we show spatially averaged response functions to several physical parameters for the observed spectral lines. They were calculated by averaging the absolute value of the response functions obtained from the atmospheric models as inferred from the Stoke...