pith. sign in

arxiv: 2604.14063 · v1 · submitted 2026-04-15 · 🌌 astro-ph.SR

Prominence Plasma Parameters Maps Inferred From Lyman β and Lyman γ Observations and Non-LTE Modelling

Y. Zhang , N. Labrosse , T. A. Kucera , S. Parenti This is my paper

Pith reviewed 2026-05-10 12:00 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords solar prominencesLyman betaLyman gammanon-LTE modelingparameter mapsSPICE instrumentoff-limb observationsplasma parameters
0
0 comments X

The pith

A new method uses contribution functions to optimize non-LTE models and create temperature, pressure, and column mass maps for solar prominences from Lyman line 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 technique to infer plasma parameters across a solar prominence by analyzing Lyman beta and Lyman gamma line intensities observed by the SPICE instrument. It generates many random 1D non-LTE models, constrains altitude and velocity, and selects those that match the observed intensities in both lines. The key innovation is using contribution functions to guide the optimization of temperature and pressure profiles, ensuring a consistent physical match. This allows pixel-by-pixel modeling to produce spatial maps of temperature, pressure, and column mass.

Core claim

By constraining altitude and radial velocity and employing a 1D non-LTE radiative transfer code to produce 1000 random models, the authors compute line profiles and compare integrated intensities with SPICE observations. They then select models that simultaneously match both Lyman β and Lyman γ lines, using contribution functions to guide the optimization of temperature and pressure profiles. This produces a set of physically consistent solutions from which parameter maps are generated for the prominence region.

What carries the argument

Contribution-function-guided selection and optimization of temperature and pressure profiles within ensembles of 1D non-LTE radiative transfer models.

Load-bearing premise

That randomly sampled 1D non-LTE models, after fixing altitude and radial velocity, can represent the actual three-dimensional prominence without major biases from the modeling assumptions or the contribution function guidance.

What would settle it

Observing significant mismatches between the derived parameter maps and independent temperature or density measurements from other instruments or techniques at the same locations.

Figures

Figures reproduced from arXiv: 2604.14063 by N. Labrosse, S. Parenti, T. A. Kucera, Y. Zhang.

Figure 1
Figure 1. Figure 1: SPICE observation of the integrated intensity of the Ly [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The process to construct initial set of prominence models that are close to observations of the Lyman [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The 9-pixel intensity map of the Lyman β and the Lyman γ line in the prominence region (the red box region in [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Comparison of best matched intensity map and observation for the red box region in Figure 1. The upper panels are obser [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: The schematic of 1D prominence slab model in the [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Temperature profiles of Model A and Model B and contribution functions of Model A for the Lyman [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: The lower-right panel of Figure 6 showing the details of [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: The comparison of calculated intensity by the best coe [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Parameter maps for the most important parameters in the Lyman [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: The electron density map and Hα intensity map as examples of other maps we can produce. Solar Orbiter is a space mission of international collaboration between ESA and NASA, operated by ESA. The EUI instru￾ment was built by CSL, IAS, MPS, MSSL/UCL, PMOD/WRC, ROB, LCF/IO with funding from the Belgian Federal Sci￾ence Policy Office (BELSPO/PRODEX PEA 4000112292 and 4000134088); the Centre National d’Etudes … view at source ↗
read the original abstract

The first dedicated observation of an off-limb prominence by the Spectral Imaging of the Coronal Environment (SPICE) instrument took place on April 15, 2023. We aim to create parameter maps on the prominence region, including temperature, pressure, and column mass, by studying the integrated intensity of the Lyman $\beta$ and Lyman $\gamma$ lines from SPICE data. After constraining the altitude and radial velocity in the prominence, we use a 1D non-LTE radiative transfer code to generate 1000 random models and compute the Lyman $\beta$ and Lyman $\gamma$ line profiles. The computed intensities are compared with observed integrated intensities from SPICE. Then, we create models which simultaneously give a reasonable match with the observed intensities in both lines. Unlike previous approaches, our method uses contribution functions to guide the optimisation of temperature and pressure profiles. Our approach enables a physically constrained and consistent match to both spectral lines. The method in this paper enables us to generate models from pixels on the prominence region and use this information to generate parameter maps. The results obtained have potential for future research.

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 manuscript presents a method to derive maps of temperature, pressure, and column mass across a solar prominence observed off-limb by SPICE on 15 April 2023. After fixing altitude and radial velocity, the authors generate 1000 random 1D non-LTE models, compute Lyman β and Lyman γ profiles, and retain those providing a reasonable simultaneous match to the observed integrated intensities; contribution functions are used to guide optimization of the temperature and pressure profiles. The resulting per-pixel models are then assembled into parameter maps for the prominence region.

Significance. If the central claim holds, the work would be moderately significant for solar physics: it offers a workflow that enforces consistency between two Lyman lines via contribution-function guidance, which could improve upon earlier non-LTE prominence analyses and be applicable to future SPICE datasets. However, the absence of any quantitative validation, error estimates, or uniqueness tests in the presented material limits the immediate impact.

major comments (2)
  1. [Abstract] Abstract: the claim that the method 'enables a physically constrained and consistent match' and produces usable parameter maps is unsupported by any reported quantitative results (goodness-of-fit values, uncertainties, or comparisons to independent diagnostics). Without these, it is impossible to judge whether the selected models are physically meaningful or merely consistent with the 1D forward model.
  2. [Modeling approach] Modeling section (implied by workflow description): the central assumption that random sampling of 1000 1D slab profiles, after fixing altitude and radial velocity, yields an unbiased ensemble whose contribution functions can be optimized to match both lines simultaneously is load-bearing. Prominences are known to be 3D with thread-like structure and possible transverse flows; no synthetic 3D tests, degeneracy analysis, or discussion of possible biases from the contribution-function guidance are provided.
minor comments (1)
  1. [Abstract] The abstract would be strengthened by including at least one example quantitative result (e.g., typical temperature or column-mass range recovered) or a statement of the achieved fit quality.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive feedback on our manuscript. We agree that strengthening the quantitative support for our claims and expanding the discussion of modeling limitations will improve the paper, and we outline targeted revisions below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim that the method 'enables a physically constrained and consistent match' and produces usable parameter maps is unsupported by any reported quantitative results (goodness-of-fit values, uncertainties, or comparisons to independent diagnostics). Without these, it is impossible to judge whether the selected models are physically meaningful or merely consistent with the 1D forward model.

    Authors: We accept this point. The current abstract and results section rely on qualitative descriptions of 'reasonable match.' In the revised manuscript we will add explicit quantitative metrics: the mean and standard deviation of the intensity ratio (modeled/observed) for the retained models across the map, the typical number of models retained per pixel from the 1000-member ensemble, and 1-sigma uncertainties on temperature, pressure and column mass derived from the spread of acceptable solutions. We will also reference typical values from prior prominence studies for context. These additions will allow readers to evaluate the physical consistency directly. revision: yes

  2. Referee: [Modeling approach] Modeling section (implied by workflow description): the central assumption that random sampling of 1000 1D slab profiles, after fixing altitude and radial velocity, yields an unbiased ensemble whose contribution functions can be optimized to match both lines simultaneously is load-bearing. Prominences are known to be 3D with thread-like structure and possible transverse flows; no synthetic 3D tests, degeneracy analysis, or discussion of possible biases from the contribution-function guidance are provided.

    Authors: The 1D slab geometry is a deliberate and standard simplification for off-limb integrated intensities, where the line of sight sums over multiple threads. Random sampling of 1000 models is used to broadly explore the temperature-pressure-column-mass space consistent with the observed altitude and velocity; the contribution-function weighting then selects only those profiles whose formation regions simultaneously reproduce both Lyman lines, thereby enforcing cross-line consistency that earlier single-line or ad-hoc methods lacked. We acknowledge that this does not capture full 3D effects or transverse flows. In revision we will insert a dedicated limitations subsection that (i) quantifies the range of acceptable models per pixel as a proxy for degeneracy, (ii) discusses how contribution-function guidance mitigates but does not eliminate biases from unresolved structure, and (iii) notes that full 3D forward modeling lies outside the present scope. This provides the requested discussion without claiming uniqueness. revision: partial

Circularity Check

0 steps flagged

No circularity: standard forward-model inversion from independent observations

full rationale

The derivation generates 1000 random 1D non-LTE models (after fixing altitude and radial velocity), computes Lyman β/γ intensities, and selects those matching observed SPICE integrated intensities via contribution-function-guided optimization of T and P profiles. The resulting parameter maps are outputs of this selection process against external data, not quantities defined by or equivalent to the inputs by construction. No self-definitional equations, fitted inputs renamed as predictions, load-bearing self-citations, uniqueness theorems, or smuggled ansatzes appear in the abstract or described chain. The method is a conventional radiative-transfer inversion and remains self-contained against the observations.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The approach rests on standard non-LTE radiative transfer assumptions and the premise that random sampling plus contribution-function guidance can yield physically consistent solutions.

free parameters (1)
  • temperature and pressure profiles
    Randomly generated and then optimized to match observed integrated intensities
axioms (1)
  • domain assumption Prominence can be adequately represented by 1D plane-parallel non-LTE models once altitude and radial velocity are constrained
    Explicitly stated as the starting point before generating the 1000 models

pith-pipeline@v0.9.0 · 5515 in / 1267 out tokens · 44203 ms · 2026-05-10T12:00:41.791548+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

24 extracted references · 24 canonical work pages

  1. [1]

    2020, A&A, 642, A14

    Anderson, M., Appourchaux, T., Auchère, F., et al. 2020, A&A, 642, A14

    work page 2020

  2. [2]

    & Heinzel, P

    Anzer, U. & Heinzel, P. 1999, A&A, 349, 974

    work page 1999

  3. [3]

    W., et al

    Barczynski, K., Schmieder, B., Peat, A. W., et al. 2021, A&A, 653, A94

    work page 2021

  4. [4]

    & Stein, R

    Carlsson, M. & Stein, R. F. 1997, The Astrophysical Journal, 481, 500

    work page 1997

  5. [5]

    & Heinzel, P

    Gouttebroze, P. & Heinzel, P. 2002, A&A, 385, 273

    work page 2002

  6. [6]

    1993, Astronomy and Astrophysics Sup- plement Series (ISSN 0365-0138), vol

    Gouttebroze, P., Heinzel, P., & Vial, J. 1993, Astronomy and Astrophysics Sup- plement Series (ISSN 0365-0138), vol. 99, no. 3, p. 513-543., 99, 513 Gunár, S., Heinzel, P., Anzer, U., & Schmieder, B. 2008, A&A, 490, 307 Gunár, S., Schwartz, P., Schmieder, B., Heinzel, P., & Anzer, U. 2010, Astron- omy & Astrophysics, 514, A43

    work page 1993

  7. [7]

    2015, in Astrophysics and Space Science Library, V ol

    Heinzel, P. 2015, in Astrophysics and Space Science Library, V ol. 415, Solar Prominences, ed. J.-C. Vial & O. Engvold, 103

    work page 2015

  8. [8]

    1994, Astronomy and Astrophysics (ISSN 0004-6361), vol

    Heinzel, P., Gouttebroze, P., & Vial, J.-C. 1994, Astronomy and Astrophysics (ISSN 0004-6361), vol. 292, no. 2, p. 656-668, 292, 656

    work page 1994

  9. [9]

    2023, The Astrophysical Journal Letters, 957, L10

    Heinzel, P., Jejˇciˇc, S., Št ˇepán, J., et al. 2023, The Astrophysical Journal Letters, 957, L10

    work page 2023

  10. [10]

    2001, A&A, 370, 281 Jejˇciˇc, S., Heinzel, P., Schmieder, B., et al

    Heinzel, P., Schmieder, B., Vial, J.-C., & Kotrˇc, P. 2001, A&A, 370, 281 Jejˇciˇc, S., Heinzel, P., Schmieder, B., et al. 2022, The Astrophysical Journal, 932, 3 Jejˇciˇc, S., Heinzel, P., Zapiór, M., et al. 2014, Solar Physics, 289, 2487

    work page 2001

  11. [11]

    2015, in Astrophysics and Space Science Library, V ol

    Labrosse, N. 2015, in Astrophysics and Space Science Library, V ol. 415, Solar Prominences, ed. J.-C. Vial & O. Engvold, 131

    work page 2015

  12. [12]

    2007, A&A, 463, 1171

    Labrosse, N., Gouttebroze, P., & Vial, J.-C. 2007, A&A, 463, 1171

    work page 2007

  13. [13]

    2010, Space Science Reviews, 151, 243

    Labrosse, N., Heinzel, P., Vial, J.-C., et al. 2010, Space Science Reviews, 151, 243

    work page 2010

  14. [14]

    & McGlinchey, K

    Labrosse, N. & McGlinchey, K. 2012, A&A, 537, A100

    work page 2012

  15. [15]

    S., Radziszewski, K., et al

    Labrosse, N., Rodger, A. S., Radziszewski, K., et al. 2022, Monthly Notices of the Royal Astronomical Society: Letters, 513, L30

    work page 2022

  16. [16]

    Levens, P. J. & Labrosse, N. 2019, A&A, 625, A30

    work page 2019

  17. [17]

    M., Del Zanna, G., & Mason, H

    Mulay, S. M., Del Zanna, G., & Mason, H. 2017, A&A, 606, A4 Müller, D., St. Cyr, O. C., Zouganelis, I., et al. 2020, A&A, 642, A1

    work page 2017

  18. [18]

    2014, Living Reviews in Solar Physics, 11, 1

    Parenti, S. 2014, Living Reviews in Solar Physics, 11, 1

    work page 2014

  19. [19]

    2015, Astronomy & Astrophysics, 577, A92 SPICE Consortium, Anderson, M., Appourchaux, T., et al

    Schwartz, P., Gunar, S., & Curdt, W. 2015, Astronomy & Astrophysics, 577, A92 SPICE Consortium, Anderson, M., Appourchaux, T., et al. 2020, A&A, 642, A14

    work page 2015

  20. [20]

    & Wiehr, E

    Stellmacher, G. & Wiehr, E. 2017, Solar Physics, 292, 83

    work page 2017

  21. [21]

    2019, IL NUOVO CIMENTO, 100, 42

    Susino, R., Bemporad, A., Heinzel, P., et al. 2019, IL NUOVO CIMENTO, 100, 42

    work page 2019

  22. [22]

    E., Avrett, E

    Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635

    work page 1981

  23. [23]

    2007, Solar Physics, 246, 327

    Vial, J.-C., Ebadi, H., & Ajabshirizadeh, A. 2007, Solar Physics, 246, 327

    work page 2007

  24. [24]

    2016, Central European Astrophysical Bulletin, 40, 79

    Wiehr, E., Stellmacher, G., & Bianda, M. 2016, Central European Astrophysical Bulletin, 40, 79

    work page 2016