Recognition: no theorem link
Polarimetric Tomography Applied to Synthetic Multi-Spacecraft White-Light Images: Observing Coronal Mass Ejections in 3D
Pith reviewed 2026-05-15 15:07 UTC · model grok-4.3
The pith
Polarimetric tomography from multiple spacecraft reconstructs coronal mass ejection density in 3D with improving accuracy.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
A discrete tomography method reconstructs three-dimensional coronal mass ejection density from synthetic multi-spacecraft white-light images. The mean relative absolute error decreases consistently as the number of observing spacecraft increases and is lower for polarimetric reconstructions than for non-polarimetric ones. The CME front can be identified with 52-72 percent accuracy and its radial position constrained to 0.003-0.005 au using three spacecraft at L1, L4, and L5 with polarimetric data. At least four spacecraft are required for accurate 3D structure, while out-of-ecliptic observers increase the invertible volume.
What carries the argument
Discrete tomography algorithm that inverts multi-viewpoint polarimetric and non-polarimetric coronagraph images into three-dimensional electron density.
Load-bearing premise
The synthetic coronagraph images generated by the CORHEL model accurately capture the true scattering and density properties of real coronal mass ejections.
What would settle it
Applying the reconstruction to actual multi-spacecraft observations and comparing the derived densities and front positions against in-situ measurements from a probe that encounters the CME.
read the original abstract
A discrete tomography method has been developed that is able to reconstruct three-dimensional coronal mass ejection (CME) density structure. We test the method by producing synthetic coronagraph imagery for three events using the CORona--HELiosphere (CORHEL) model. We combine images from different numbers of observing spacecraft and we perform the method separately using polarimetric and non-polarimetric reconstructions, as a means to test their relative effectiveness. We show that increasing the number of observing spacecraft consistently reduces the mean relative absolute error (MRAE) between the simulated and reconstructed density. Furthermore, the MRAE is generally lower when using polarimetric reconstructions compared to non-polarimetric reconstructions. Methods applied to localise the CME front work well for all spacecraft configurations, and are improved when using polarimetric, over non-polarimetric, reconstructions. The presence of a CME front can be identified with an accuracy of $(72\pm9)\%$, $(70\pm8)\%$ and $(52\pm12)\%$ for CME1, CME2 and CME3 via polarimetric reconstructions using only three spacecraft at L1, L4 and L5. The radial position of the CME front can be constrained to a high level of precision when using polarimetric reconstructions using the same three spacecraft; $0.003\pm0.004$\,au, $0.004\pm0.005$\,au and $0.005\pm0.004$\,au for CME1, CME2 and CME3, respectively. We expect that at least four spacecraft are required in order to derive accurate information about 3D CME structure. We find no strong evidence of improvement when including out-of-ecliptic observers, but that their inclusion increases the volume of space within which the inversion can be performed.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a discrete tomography method for reconstructing three-dimensional coronal mass ejection (CME) density structures from synthetic white-light coronagraph images generated by the CORHEL model. The approach is tested across varying numbers of spacecraft observers, with separate polarimetric and non-polarimetric reconstructions, showing that MRAE between simulated and reconstructed densities decreases with more spacecraft and is generally lower for polarimetric cases. CME front localization achieves 72±9%, 70±8%, and 52±12% identification accuracy for three events with three spacecraft (L1/L4/L5) using polarimetry, and radial front positions are constrained to 0.003–0.005 au precision; the authors conclude that at least four spacecraft are needed for accurate 3D structure and note limited benefit from out-of-ecliptic observers beyond increased inversion volume.
Significance. If the reported MRAE reductions and front-localization precisions hold under broader conditions, the method offers a concrete path toward 3D CME density mapping with planned multi-spacecraft missions (e.g., L4/L5 constellations), providing quantitative benchmarks absent from single-spacecraft or non-polarimetric approaches. The use of independent CORHEL forward-model data for validation and the systematic comparison of spacecraft counts and polarimetry are clear strengths that support reproducibility within the synthetic setting.
major comments (3)
- [Abstract / Results] Abstract and results: the assertion that 'at least four spacecraft are required in order to derive accurate information about 3D CME structure' is not directly supported by the quantitative results, which demonstrate low MRAE and 0.003–0.005 au front precision already with three spacecraft at L1/L4/L5; the threshold should be justified by additional tests or explicit criteria for 'accurate' 3D structure.
- [Methods] Methods: the discrete tomography inversion procedure, including the precise formulation of the forward operator, any regularization, and how polarimetric versus total-brightness constraints are implemented, is not described in sufficient detail to allow independent reproduction or assessment of conditioning for the reported MRAE values.
- [Results] Results: the reported front-position uncertainties (e.g., 0.003±0.004 au) lack an explicit derivation or error-propagation analysis; it is unclear whether these reflect reconstruction variance across events, noise realizations, or another source.
minor comments (2)
- [Abstract] Abstract: specify the exact spacecraft configurations tested (beyond the L1/L4/L5 triplet) and whether out-of-ecliptic positions were included in the MRAE and front-localization statistics.
- [Discussion] Discussion: the statement that out-of-ecliptic observers increase the inversion volume should be accompanied by quantitative volume metrics or a supporting figure.
Simulated Author's Rebuttal
We thank the referee for their constructive review and positive assessment of the manuscript's strengths. We address each major comment below with point-by-point responses and indicate the planned revisions.
read point-by-point responses
-
Referee: [Abstract / Results] Abstract and results: the assertion that 'at least four spacecraft are required in order to derive accurate information about 3D CME structure' is not directly supported by the quantitative results, which demonstrate low MRAE and 0.003–0.005 au front precision already with three spacecraft at L1/L4/L5; the threshold should be justified by additional tests or explicit criteria for 'accurate' 3D structure.
Authors: We agree that the phrasing in the abstract requires clarification to avoid overstatement. The expectation of at least four spacecraft stems from our assessment that three observers yield accurate front localization but insufficient fidelity for the full 3D density morphology (higher MRAE and incomplete core reconstruction in the synthetic tests). In the revision we will replace the statement with an explicit criterion (MRAE below 0.15 combined with visual agreement on both front and internal structure) and add a short justification in the discussion section based on the existing results, without performing new tests. revision: yes
-
Referee: [Methods] Methods: the discrete tomography inversion procedure, including the precise formulation of the forward operator, any regularization, and how polarimetric versus total-brightness constraints are implemented, is not described in sufficient detail to allow independent reproduction or assessment of conditioning for the reported MRAE values.
Authors: We acknowledge the need for greater methodological transparency. The revised Methods section will include the mathematical statement of the discrete tomography problem, the explicit form of the forward operator (line-of-sight integration weighted by Thomson scattering geometry), the regularization term (L2 gradient penalty), and the distinct constraint matrices applied to polarized-brightness versus total-brightness data. These additions will enable independent reproduction and conditioning assessment. revision: yes
-
Referee: [Results] Results: the reported front-position uncertainties (e.g., 0.003±0.004 au) lack an explicit derivation or error-propagation analysis; it is unclear whether these reflect reconstruction variance across events, noise realizations, or another source.
Authors: The quoted uncertainties are the standard deviations of the front-position estimates obtained across the three simulated CME events for the three-spacecraft (L1/L4/L5) configuration. We will insert a clarifying sentence in the Results section stating this source and noting that the study uses noise-free synthetic images, so no noise-realization variance is included. A formal error-propagation derivation is not present and will be flagged as future work rather than added at this stage. revision: partial
Circularity Check
Synthetic validation uses independent forward model with no fitted inputs or self-referential definitions
full rationale
The paper generates synthetic coronagraph images via the external CORHEL model for three specific CME events, then applies its discrete tomography inversion separately to polarimetric and non-polarimetric versions of those images. All reported metrics (MRAE between simulated and reconstructed density, CME front identification accuracy, and radial position precision of 0.003-0.005 au) are obtained by direct comparison against the known ground-truth density field supplied by CORHEL. No parameters of the inversion are tuned to the reconstruction targets, no self-citation supplies a uniqueness theorem or ansatz that forces the result, and the forward operator is not derived from the same data being inverted. The derivation chain therefore remains independent of the quantitative claims.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption CORHEL model images accurately represent real CME density and Thomson scattering
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.