REVIEW 4 major objections 4 minor 56 references
XRISM observations of solar flare X-ray emission reflected in the Earth's atmosphere
T0 review · 4 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read One year of XRISM day-Earth occultation data measures solar flare metal abundances and separates Fe-K fluorescence lines.
desk verdict The paper's core point—XRISM can do solar flare science from reflected atmospheric X-rays—holds up, but the monotonic abundance trends and Fe-K anti-correlation are stated more strongly than the paper's own tables support. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing method is an equivalent-width-to-abundance conversion: line-to-continuum ratios are measured from stacked day-Earth spectra, then converted to elemental abundances assuming a plasma whose emission measure follows one power law in temperature, with the power-law slope pinned by the observed Si XIV/XIII flux ratio. This is the same route used by the Suzaku study the paper extends. For Fe-K, the key object is the atmospheric reflection model: a Monte Carlo calculation of Rayleigh and Compton scattering by N, O, and Ar that predicts the scattered Fe XXIV/XXV line shapes and the Rayleigh/Compton ratio (≈0.543), leaving a 6.4–6.5 keV excess that is modeled as low-ionized Fe XXI Kα
What would settle it
A higher-statistics Resolve observation of one bright X-class flare that resolves the 6.4–6.5 keV excess into the predicted Rayleigh/Compton scattered-line structure rather than a discrete Fe XXI Kα line would undercut the fluorescence-line identification; alternatively, a direct multi-temperature DEM reconstruction from the same spectra yielding abundances that differ by more than the quoted systematic would falsify the abundance trends.
Extended reading notes
Core claim
Central claim: reflected solar flare X-rays recorded during XRISM's day-Earth occultations form a usable solar dataset. Stacked Xtend spectra give abundances of Mg, Si, S, Ar, Ca, and Fe for M1–X10 flares; the pattern shows the inverse-FIP effect (low first-ionization-potential elements depleted), with Si, S, and Ar decreasing as flare magnitude rises, matching ponderomotive-model predictions while Ca rises. Resolve's Fe-K spectra resolve Rayleigh- and Compton-scattered Fe XXIV/XXV lines from neutral or low-ionized Fe Kα, whose equivalent width anti-correlates with 7.11–9.20 keV flux (slope -0.14 ± 0.09), favoring hard X-ray photoionization as the fluorescence driver while the electron-colli
Load-bearing premise
The abundance results rest on assuming the flare's temperature distribution is a single power law with slope fixed by one line ratio; if the true distribution is shaped differently, the abundances and their flare-class trends could shift beyond the quoted errors, and the roughly 20% systematic is borrowed from a previous study rather than measured here.
Editorial extensions
If this is right
- Day-Earth occultation data become a free, long-running solar flare monitor for the full XRISM mission.
- Abundance changes can be followed at roughly 100-second resolution, revealing when low-FIP elements enter the flaring loop relative to the flare peak.
- The Fe-K band provides a geometric diagnostic that separates atmospheric scattering from solar fluorescence, usable to probe hard X-ray irradiation of the lower solar atmosphere.
- The flare-class abundance trends for Si, S, Ar, and Ca give quantitative targets for chromospheric evaporation and dredge-up models.
- The Fe Kα equivalent width can serve as an X-ray flux proxy for flares when direct hard X-ray measurements are unavailable.
Reading between the lines
- If the reflection model is right, equivalent-width measurements are insensitive to atmospheric column fluctuations, so other pointed X-ray missions with day-Earth data could apply the same method without solar-dedicated hardware.
- The solar Fe Kα slope agrees with the stellar-flare slope, suggesting a common photoionization fluorescence mechanism across roughly five orders of magnitude in X-ray flux; a broader stellar sample would test whether the scaling holds.
- The pre-peak low-FIP enrichment seen in three X-class flares could be tied to the onset of chromospheric evaporation; correlating it with microwave or hard X-ray burst onset would test whether the enrichment is a cause or a consequence of flare triggering.
- A longer baseline through the declining solar cycle would extend the abundance-flare-class relation down to C-class flares and test whether the trends continue, saturate, or reverse.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper analyzes XRISM day-Earth occultation data to study solar flare X-rays reflected by the Earth's atmosphere. The authors stack one year of Xtend and Resolve spectra by GOES flare class, identify emission lines from Mg, Si, S, Ar, Ca, and Fe, and convert line equivalent widths to elemental abundances using a single power-law differential emission measure (DEM) model with the slope fixed by the Si XIV/XIII ratio. They report an inverse-FIP abundance pattern and a flare-magnitude dependence for Si, S, and Ar, trace short-term abundance variations around several X-class flares, and decompose the Fe-K region into Rayleigh- and Compton-scattered Fe XXIV/XXV lines plus a neutral/low-ionized Fe K-alpha component. They find an anti-correlation between the Fe K-alpha equivalent width and hard X-ray flux, and argue that XRISM can serve as a solar flare observatory by reusing occultation data.
Significance. If the results hold, the paper demonstrates a new, high-capability window into solar flare spectroscopy from XRISM, extending to the Fe-K band with both a large-grasp CCD (Xtend) and a high-resolution microcalorimeter (Resolve). The data reduction is transparent, the reflected origin of the emission is convincingly validated through vignetting profiles (Figures 6-7), and the Fe-K line decomposition is an interesting technical achievement. The paper also provides a large stacked flare dataset that can be used for further studies. However, the headline claims about abundance trends and the Fe K-alpha anti-correlation are not supported by the data with the stated significances, and the abundance analysis rests on an unvalidated single-power-law DEM assumption. These issues must be addressed before the conclusions can be accepted.
major comments (4)
- [§3.3, Table 3] The abstract and §4.2 claim that the abundances of Si, S, and Ar decrease with increasing flare magnitude. Table 3 does not show this monotonic behavior: Ar rises from 2.75±0.60 (M1–5) to 3.59±0.45 (M5–10) before falling to 1.54±0.43 (X5–10), and Si rises from 0.82±0.04 (X1–5) to 1.00±0.08 (X5–10). Mg is also non-monotonic (1.24, 0.77, 0.67, 1.86). Since the X5–10 bin uses only 0.9 ks of exposure (Table 1), the trend is fragile. This claim should be weakened to a possible trend or re-derived with a more appropriate fitting/binning scheme.
- [§3.3, 'Abundance measurement of flare loops'] The central abundance analysis uses a single power-law DEM model (EM∝(kT_e)^α) with α pinned by the Si XIV/XIII flux ratio. The paper itself notes that the 1.3–5.0 keV and 5–9 keV bands require different power-law continua, 'likely indicat[ing] that the observed spectra can be approximated with multiple DEM components.' This internal inconsistency shows that the single power-law DEM is not validated over the full bandpass. The associated systematic uncertainty (<20%) is inherited from Katsuda et al. (2020) and not recomputed for the XRISM line set and reflection geometry. A multi-thermal DEM could bias the derived abundances and the flare-class trends, especially for high-temperature elements S, Ar, Ca, and Fe. The authors should implement a multi-thermal or full DEM reconstruction or explicitly quantify the bias.
- [§3.3, Figure 13] For the M1–5 class, the pre-flare background exceeds the flare spectrum (Figure 13 alt text). Background subtraction in this regime leaves a low-S/N net spectrum and can bias line equivalent widths. The M1–5 abundances in Table 3 are therefore not reliable as a baseline for the flare-class trend. The authors should assess this bias (e.g., via simulations) or exclude the M1–5 bin from the trend analysis.
- [§4.1, Figure 12] The anti-correlation between the neutral/low-ionized Fe Kα equivalent width and 7.11–9.20 keV flux is based on a best-fit power-law slope of −0.14±0.09. This is consistent with zero at <2σ, so the data do not establish an anti-correlation. The comparison with Inoue et al. (2025) (−0.27±0.10) still leaves both measurements consistent with no trend. The conclusion that hard X-rays stimulate Fe Kα fluorescence should be presented as tentative and not as a detection.
minor comments (4)
- [Abstract, §2] The date 'November 31, 2024' appears in the abstract and Section 2; November has 30 days. Please correct to November 30, 2024.
- [§2] The word 'Chimera' is capitalized in the text ('a Chimera of multiple observation modes'); it should likely be lowercase 'chimera' unless intended as a proper noun.
- [References] The reference to Inoue et al. (2025) is incomplete; it lacks the journal, volume, and page/article number. Please provide full citation details.
- [Figure 3 caption] The caption text is garbled: 'Ar K , XVII, K' and 'Fe K , XXV' are not readable. Please fix the notation and clarify the line labels.
Circularity Check
Fe-K anti-correlation is partly built into the EW definition; DEM systematics are inherited from a self-cited prior paper.
-
self definitional
[Sec. 4.1, Figure 12 (and Sec. 3.3 spectral model)]
"Figure 12 shows the Fe-Kα equivalent width vs. 7.11–9.20 keV flux for four flare magnitudes. The decreasing Fe-Kα equivalent width as a function of flare magnitude, while the values of FeXXIV+XXV are almost constant, suggests that the neutral/low-ionized Fe Kα is stimulated by photoionization."
The EW is measured as the flux of the FeXXI line at 6.537 keV divided by the fitted power-law continuum, and the 7.11–9.20 keV flux is the integral of that same fitted continuum. For a power-law continuum C(E)=N E^{-Γ}, EW = F_line/[N (6.537)^{-Γ}] and H = N I(Γ), so EW = F_line * const(Γ) / H. Thus plotting EW against H is plotting a quantity inversely proportional to the continuum against the continuum itself. An anti-correlation is mathematically expected unless F_line grows faster than linearly with H. The paper does not use an independent hard X-ray measurement, so this 'anti-correlation' is not an independent observable and partially reduces to the definition of equivalent width.
-
self citation load bearing
[Sec. 3.3, paragraph on time variation and DEM model]
"According to Katsuda et al. (2020), the systematic uncertainties in the resultant abundances due to the uncertain temperature distribution would be <20%, which is similar to the statistical uncertainties as one can see below, and thus would not affect our conclusion."
The conclusion that the DEM-model assumption does not threaten the abundance results rests on a <20% systematic uncertainty estimate imported from Katsuda et al. (2020), a paper with overlapping authors that used the same power-law DEM method. The XRISM data themselves are not used to validate the DEM shape against, e.g., a two-temperature model in this bandpass. The self-citation is load-bearing for the error bars on the central abundance claims, though it is an external, published estimate rather than a derivation.
full rationale
The abundance measurement is a forward-model inversion: equivalent widths are converted to abundances using a power-law DEM whose slope α is pinned by the observed Si XIV/XIII ratio. This is not circular because the Si ratio constrains the DEM independently of the abundances, and the XRISM data are new. The main concern is the Fe-K anti-correlation: the EW and the 7.11–9.20 keV flux are both derived from the same fitted power-law continuum, so the anti-correlation is partly built into the ratio definition. A secondary issue is that the <20% DEM systematic is inherited from a self-cited prior paper rather than derived from XRISM data; this is load-bearing for the error bars but the cited work is an external, published estimate. Overall the central claims (inverse-FIP pattern, flare-class trends, Fe-K decomposition) have independent content, so the circularity is mild.
Assumptions & free parameters
free parameters (4)
- DEM slope alpha (power-law DEM model) =
-0.85 (M1-5), -0.398 (M5-10), 0.0352 (X1-5), -0.48 (X5-10); Table 2
- Fe-K reflection model temperature =
1.5 keV (fixed)
- Fe-K reflection scattering angle =
90 degrees (observed mean about 95 plus or minus 20 degrees)
- Power-law slope of Fe K alpha equivalent width vs 7.11-9.20 keV flux =
-0.14 plus or minus 0.09
assumptions (5)
- domain assumption The day-Earth occultation X-ray signal is dominated by solar flare X-rays reflected from a spatially uniform (flat-field) Earth atmosphere.
- domain assumption Line equivalent widths are not changed by reflection in Earth's atmosphere, so they can be converted to abundances using an unreflected plasma model.
- ad hoc to paper A single power-law DEM model (EM proportional to (kT_e)^alpha) with alpha from the Si XIV/XIII ratio is an adequate representation of the flaring plasma.
- domain assumption Pre-flare spectra (up to 10^4 seconds before each flare) are a clean background representing non-flare emission.
- standard math The atomic and plasma codes (AtomDB v3.0.9, apec, cevmkl, XSPEC) produce correct line emissivities and continua.
Cite this review
Pith. "Pith review of XRISM observations of solar flare X-ray emission reflected in the Earth's atmosphere." pith.science (2026). https://pith.science/paper/QRV7NH6M
@misc{pith2026250905029,
author = {Pith},
title = {Pith review of: XRISM observations of solar flare X-ray emission reflected in the Earth's atmosphere},
year = {2026},
howpublished = {\url{https://pith.science/paper/QRV7NH6M}},
note = {Machine review of arXiv:2509.05029}
}
abstract
The X-ray Imaging and Spectroscopy Mission (XRISM), launched into low-Earth orbit in 2023, observes the reflection of solar flare X-rays in the Earth's atmosphere as a by-product of celestial observations. Using a $\sim$one-year data set covering from October 2023 to November 2024, we report on our first results of the measurement of the metal abundance pattern and high-resolution Fe-K spectroscopy. The abundances of Mg, Si, S, Ar, Ca, and Fe measured with the CCD detector Xtend during M- and X-class flares show the inverse-first-ionization-potential (inverse-FIP) effect, which is consistent with the results of Katsuda et al., ApJ, 2020 using the Suzaku satellite. The abundances of Si, S, and Ar are found to decrease with increasing flare magnitude, which is consistent with the theoretical model by Laming (Laming, ApJ, 2021), whereas Ca exhibits an opposite trend. The large effective area and field of view of Xtend allow us to trace the evolution of the abundances in several X-class flare loops on a timescale of a few 100 s, finding an enrichment of low-FIP elements before flare peaks. The high-resolution Fe-K spectrum obtained with the microcalorimeter Resolve successfully separates the Rayleigh- and Compton-scattered Fe XXIV/XXV lines and neutral or low-ionized Fe-K$\alpha$ lines. The neutral/low-ionized Fe-K$\alpha$ equivalent width shows an anti-correlation with hard X-ray flux with the best-fit power-law slope of $-0.14 \pm 0.09$, suggesting that hard X-rays from flare loops are stimulating the Fe K$\alpha$ fluorescence. This work demonstrates that XRISM can be a powerful tool in the field of solar physics, offering valuable high-statistic CCD data and high-resolution microcalorimeter spectra in the energy range extending to the Fe-K band.
Figures
Figures from the paper (11 more)
Reference graph
Works this paper leans on
-
[1]
Arnaud, K.A., 1996, in Astronomical Data Analysis Software and Systems V , eds. G.H. Jacoby & J. Barnes, volume 101 ofAstronomical Society of the Pacific Conference Series, 17
work page 1996
- [2]
-
[3]
Basko, M.M., 1978, ApJ, 223, 268
work page 1978
-
[4]
Brooks, D.H., Ugarte-Urra, I., & Warren, H.P., 2015, Nature Communications, 6, 5947
work page 2015
-
[5]
Caspi, A., McTiernan, J.M., & Warren, H.P., 2014, ApJL, 788, 2, L31
work page 2014
-
[6]
Churazov, E., Sazonov, S., Sunyaev, R., & Revnivtsev, M., 2008, MNRAS, 385, 2, 719
work page 2008
-
[7]
Nittler, L.R., 2015, ApJ, 803, 2, 67
work page 2015
-
[8]
Doschek, G.A., Meekins, J.F., Kreplin, R.W., Chubb, T.A., & Friedman, H., 1971, ApJ, 170, 573
work page 1971
Show all 56 references
-
[9]
& Warren, H.P., 2019, ApJ, 884, 2, 158
Doschek, G.A. & Warren, H.P., 2019, ApJ, 884, 2, 158
2019
-
[10]
Doschek, G.A., Warren, H.P., & Feldman, U., 2015, ApJL, 808, 1, L7
2015
-
[11]
Ercolano, B., Drake, J.J., Reale, F., Testa, P., & Miller, J.M., 2008, ApJ, 688, 2, 1315
2008
-
[12]
Scr., 46, 3, 202
Feldman, U., 1992, Phys. Scr., 46, 3, 202
1992
-
[13]
Feldman, U., Doschek, G.A., Behring, W.E., & Phillips, K.J.H., 1996, ApJ, 460, 1034
1996
-
[14]
Feldman, U., Doschek, G.A., & Kreplin, R.W., 1980, ApJ, 238, 365
1980
-
[15]
Foster, A.R., Ji, L., Smith, R.K., & Brickhouse, N.S., 2012, ApJ, 756, 2, 128 HEASARC, 2014, HEAsoft: Unified Release of FTOOLS and XANADU
2012
-
[16]
Gendreau, K., 2025 Publications of the Astronomical Society of Japan(0000), Vol. 00, No. 011
2025
-
[17]
et al., 2018, Journal of Low Temperature Physics, 193, 5, 991
Ishisaki, Y . et al., 2018, Journal of Low Temperature Physics, 193, 5, 991
2018
-
[18]
& Mewe, R., 1993, A&AS, 97, 2, 443
Kaastra, J.S. & Mewe, R., 1993, A&AS, 97, 2, 443
1993
-
[19]
Katsuda, S., Ohno, M., Mori, K., et al., 2020, ApJ, 891, 2, 126
2020
-
[20]
Katsuda, S., Shinagawa, H., Fujiwara, H., et al., 2024, Geophys. Res. Lett., 51, 20, e2024GL112025
2024
-
[21]
Kepa, A., Sylwester, B., Sylwester, J., Gryciuk, M., & Siarkowski, M., 2018, Journal of Atmospheric and Solar-Terrestrial Physics, 179, 545
2018
-
[22]
Laming, J.M., 2004, ApJ, 614, 2, 1063
2004
-
[23]
Laming, J.M., 2012, ApJ, 744, 2, 115
2012
-
[24]
Laming, J.M., 2015, Living Reviews in Solar Physics, 12, 1, 2
2015
-
[25]
Laming, J.M., 2021, ApJ, 909, 1, 17
2021
-
[26]
Lodders, K., 2003, ApJ, 591, 2, 1220
2003
-
[27]
Mithun, N.P.S., Vadawale, S.V ., Del Zanna, G., et al., 2022, ApJ, 939, 2, 112
2022
-
[28]
& Shinagawa, H., 2023, Earth, Planets and Space, 75, 1, 68
Miyoshi, Y . & Shinagawa, H., 2023, Earth, Planets and Space, 75, 1, 68
2023
-
[29]
Mondal, B., Sarkar, A., Vadawale, S.V ., et al., 2021, ApJ, 920, 1, 4 Mori et al., 2024, in Space Telescopes and Instrumentation 2024: Ultraviolet to Gamma Ray, eds. J.W.A. den Herder, S. Nikzad, & K. Nakazawa, International Society for Optics and Photonics (SPIE), volume 1309...
2021
-
[30]
et al., 2022, International Society for Optics and Photonics (SPIE), volume 12181, 121811T
Mori, K. et al., 2022, International Society for Optics and Photonics (SPIE), volume 12181, 121811T
2022
-
[31]
et al., 2018, Publications of the Astronomical Society of Japan, 70, 2, 21
Nakajima, H. et al., 2018, Publications of the Astronomical Society of Japan, 70, 2, 21
2018
-
[32]
Phys., 298, 4, 55
Nama, L., Mondal, B., Narendranath, S., & Paul, K.T., 2023, Sol. Phys., 298, 4, 55
2023
-
[33]
Phys., 18, 3, 474
Neupert, W.M., 1971, Sol. Phys., 18, 3, 474
1971
-
[34]
Neupert, W.M., Gates, W., Swartz, M., & Young, R., 1967, ApJL, 149, L79
1967
-
[35]
Noda, H., Mori, K., Tomida, H., et al., 2025, PASJ, accepted
2025
-
[36]
Osten, R.A., Godet, O., Drake, S., et al., 2010, ApJ, 721, 1, 785
2010
-
[37]
Phillips, K.J.H., 2004, ApJ, 605, 2, 921
2004
-
[38]
& Miyoshi, Y ., 2024, Earth, Planets and Space, 76, 1, 15
Shinagawa, H. & Miyoshi, Y ., 2024, Earth, Planets and Space, 76, 1, 15
2024
-
[39]
& Churazov, E.M., 1996, Astronomy Letters, 22, 5, 648
Sunyaev, R.A. & Churazov, E.M., 1996, Astronomy Letters, 22, 5, 648
1996
-
[40]
Suzuki, H., Tsuji, N., Kanemaru, Y ., et al., 2025, ApJL, 978, 2, L20
2025
-
[41]
Suzuki, H., Yoneyama, T., Kobayashi, S.B., et al., 2024, in Space Telescopes and Instrumentation 2024: Ultraviolet to Gamma Ray, eds. J.W.A. den
2024
-
[42]
Sylwester, J., Lemen, J.R., Bentley, R.D., Fludra, A., & Zolcinski, M.C., 1998, ApJ, 501, 1, 397
1998
-
[43]
Sylwester, J., Lemen, J.R., & Mewe, R., 1984, Nature, 310, 5979, 665
1984
-
[44]
Sylwester, J., Sylwester, B., Phillips, K.J.H., & K˛ epa, A., 2022, ApJ, 930, 1, 77
2022
-
[45]
Tanaka, K., Watanabe, T., & Nitta, N., 1984, ApJ, 282, 793
1984
-
[46]
& Zirin, H., 1985, ApJ, 299, 1036
Tanaka, K. & Zirin, H., 1985, ApJ, 299, 1036
1985
-
[47]
Tashiro, M., Maejima, H., Toda, K., et al., 2018, in Space Telescopes and Instrumentation 2018: Ultraviolet to Gamma Ray, eds. J.W.A. den Herder, S. Nikzad, & K. Nakazawa, volume 10699 ofSociety of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 1069922 Tashi...
2018
-
[48]
et al., 2020, International Society for Optics and Photonics (SPIE), volume 11444, 1144422
Tashiro, M. et al., 2020, International Society for Optics and Photonics (SPIE), volume 11444, 1144422
2020
-
[49]
Micela, G., & Garcia-Alvarez, D., 2008, ApJL, 675, 2, L97
2008
-
[50]
Uchida, H., Mori, K., Tomida, H., et al., 2025, PASJ, accepted
2025
-
[51]
Vadawale, S.V ., Mondal, B., Mithun, N.P.S., et al., 2021, ApJL, 912, 1, L12
2021
-
[52]
Warren, H.P., 2014, ApJL, 786, 1, L2
2014
-
[53]
Woods, T.N., Caspi, A., Chamberlin, P.C., et al., 2017, ApJ, 835, 2, 122
2017
-
[54]
Woods, T.N., Schwab, B., Sewell, R., et al., 2023, ApJ, 956, 2, 94
2023
-
[55]
Yi, F., Yu, C., Zhang, S., Yue, X., He, Y ., Huang, C., Zhang, Y ., & Huang, K., 2009, Journal of Geophysical Research (Atmospheres), 114, D1, D01301
2009
-
[56]
13.Pre-flare background spectra (dashed lines) and background-subtracted flare spectra (crosses) for four different flare magnitudes
Zarro, D.M., Dennis, B.R., & Slater, G.L., 1992, ApJ, 391, 865 0.01 0.1 1 10 2 3 4 5 6 7 8 9 Counts s -1 keV -1 Energy (keV) Xtend, all CCDs M1–5 M5–10 X1–5 X5–10 +: Background-subtracted spectra - -: Background (pre flare) Fig. 13.Pre-flare background spectra (dashed lines) an...
1992
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.