Recognition: 1 theorem link
· Lean TheoremChromospheric and photospheric properties of sunspots as inferred from Stokes inversions under magneto-hydrostatic and non-local-thermodynamic equilibrium
Pith reviewed 2026-05-15 13:33 UTC · model grok-4.3
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.
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
- 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
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.
Referee Report
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)
- [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.
- [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)
- [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
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
-
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
-
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
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
axioms (2)
- domain assumption The sunspot atmosphere is in 3D magneto-hydrostatic equilibrium
- domain assumption Non-LTE effects must be solved for the selected spectral lines
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking echoes?
echoesECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.
three-dimensional magneto-hydrostatic (MHS) equilibrium to constrain the gas pressure and density
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]
Anstee, S. D. & O’Mara, B. J. 1995, MNRAS, 276, 859
work page 1995
-
[2]
Barklem, P. S. & O’Mara, B. J. 1997, MNRAS, 290, 102
work page 1997
- [3]
-
[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
work page 1969
-
[5]
Borrero, J. M. & Ichimoto, K. 2011, Living Reviews in Solar Physics, 8, 4
work page 2011
- [6]
-
[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
work page 2019
-
[8]
M., Pastor Yabar, A., & Ruiz Cobo, B
Borrero, J. M., Pastor Yabar, A., & Ruiz Cobo, B. 2021, A&A, 647, A190
work page 2021
-
[9]
M., Pastor Yabar, A., Schmassmann, M., et al
Borrero, J. M., Pastor Yabar, A., Schmassmann, M., et al. 2025, A&A, 699, A149
work page 2025
-
[10]
Borrero, J. M. & Solanki, S. K. 2008, ApJ, 687, 668
work page 2008
-
[11]
Borrero, J. M., Solanki, S. K., Bellot Rubio, L. R., Lagg, A., & Mathew, S. K. 2004, A&A, 422, 1093
work page 2004
-
[12]
Bose, S., Henriques, V . M. J., Rouppe van der V oort, L., & Pereira, T. M. D. 2019, A&A, 627, A46
work page 2019
-
[13]
Bruls, J. H. M. J., Rutten, R. J., & Shchukina, N. G. 1992, A&A, 265, 237
work page 1992
-
[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
work page 2018
-
[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
work page 2013
-
[16]
Deming, D., Boyle, R. J., Jennings, D. E., & Wiedemann, G. 1988, ApJ, 333, 978
work page 1988
- [17]
- [18]
-
[19]
Gary, G. A. & Hagyard, M. J. 1990, Sol. Phys., 126, 21
work page 1990
-
[20]
Georgoulis, M. K. 2005, ApJ, 629, L69
work page 2005
-
[21]
Henriques, V . M. J., Mathioudakis, M., Socas-Navarro, H., & de la Cruz Rodr´ıguez, J. 2017, ApJ, 845, 102
work page 2017
- [22]
-
[23]
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., &
work page 2015
- [24]
-
[25]
Kaithakkal, A. J., Borrero, J. M., Pastor Yabar, A., & de la Cruz Rodr ´ıguez, J. 2023, MNRAS, 521, 3882
work page 2023
-
[26]
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
work page 2008
-
[27]
J., Reardon, K., Carlsson, M., & Hansteen, V
Leenaarts, J., Rutten, R. J., Reardon, K., Carlsson, M., & Hansteen, V . 2010, ApJ, 709, 1362
work page 2010
- [28]
-
[29]
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...
work page 1993
-
[30]
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
work page 2003
-
[31]
Morosin, R., de la Cruz Rodr´ıguez, J., Vissers, G. J. M., & Yadav, R. 2020, A&A, 642, A210
work page 2020
-
[32]
N`emec, N.-E., Porqueras-Le´on, `O., Ribas, I., & Shapiro, A. I. 2026, A&A, 705, A111
work page 2026
-
[33]
Nye, A., Bruning, D., & Labonte, B. J. 1988, Sol. Phys., 115, 251
work page 1988
-
[34]
Osborne, C. M. J. & Mili´c, I. 2021, ApJ, 917, 14
work page 2021
-
[35]
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
work page 1979
-
[36]
Rempel, M. & Borrero, J. 2021, Oxford Research Encyclopedia of Physics S´anchez Cuberes, M., Puschmann, K. G., & Wiehr, E. 2005, A&A, 440, 345
work page 2021
-
[37]
B., Narayan, G., Hillberg, T., et al
Scharmer, G. B., Narayan, G., Hillberg, T., et al. 2008, ApJ, 689, L69
work page 2008
-
[38]
Shine, R. A. & Linsky, J. L. 1974, Sol. Phys., 39, 49
work page 1974
- [39]
-
[40]
Socas-Navarro, H., de la Cruz Rodr´ıguez, J., Asensio Ramos, A., Trujillo Bueno, J., & Ruiz Cobo, B. 2015, A&A, 577, A7
work page 2015
-
[41]
Socas-Navarro, H., Ruiz Cobo, B., & Trujillo Bueno, J. 1998, ApJ, 507, 470
work page 1998
-
[42]
Socas-Navarro, H., Trujillo Bueno, J., & Ruiz Cobo, B. 2000, Science, 288, 1396
work page 2000
-
[43]
Solanki, S. K. 2003, A&A Rev., 11, 153
work page 2003
-
[44]
Stanchfield, II, D. C. H., Thomas, J. H., & Lites, B. W. 1997, ApJ, 477, 485
work page 1997
-
[45]
Sumida, V . Y . D., Estrela, R., Swain, M., & Valio, A. 2026, A&A, 706, A281
work page 2026
-
[46]
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
work page 1993
-
[47]
Vaquero, J. M. & V ´azquez, M. 2009, The Sun Recorded Through History: Scientific Data Extracted from Historical Documents, V ol. 361 (Springer New
work page 2009
-
[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
work page 2007
-
[49]
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
work page 1981
-
[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...
work page 1969
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.