SuNeRF-CME: Physics-Informed Neural Radiance Fields for Tomographic Reconstruction of Coronal Mass Ejections
Pith reviewed 2026-05-18 15:25 UTC · model grok-4.3
The pith
Physics-informed neural radiance fields recover the full 3D structure and motion of coronal mass ejections from only two viewpoints.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
SuNeRF-CME uses Neural Radiance Fields to estimate electron density in the heliosphere through a ray-tracing approach that accounts for Thomson scattering. The model is optimized by fitting time-dependent observational data while applying physical constraints on plasma continuity, propagation direction, and speed. On synthetic observations of a CME simulation the method estimates parameters from only two viewpoints with a mean velocity error of 3.01±1.94% and propagation direction errors of 3.39±1.94° in latitude and 1.76±0.79° in longitude. It also produces a full 3D reconstruction that correctly models the three-part structure, deformed CME front, and internal plasma variations, with extra
What carries the argument
SuNeRF-CME, a Neural Radiance Field framework that performs ray-tracing for Thomson scattering while enforcing continuity, propagation direction, and speed constraints on heliospheric plasma during time-dependent optimization.
If this is right
- The reconstruction correctly captures the three-part structure of the CME along with its deformed front and internal plasma variations from two viewpoints.
- Additional viewpoints integrate directly and improve the modeled plasma distribution without retraining the core framework.
- The same physical constraints enable reliable parameter recovery even for complex or time-varying plasma distributions in the simulation.
- The approach supports advanced space weather monitoring by providing dynamic 3D models of CMEs that can be updated as new observations arrive.
Where Pith is reading between the lines
- The method could be applied to real multi-spacecraft coronagraph data to test how well the constraints handle instrument noise and calibration differences absent from simulations.
- Similar continuity and speed constraints might extend the framework to reconstruct other large-scale heliospheric features such as solar wind streams from sparse viewpoints.
- Coupling the density reconstruction with separate magnetic field models could yield testable predictions of CME arrival times and geoeffectiveness.
Load-bearing premise
The physical constraints on continuity, propagation direction, and speed of the heliospheric plasma are sufficient to overcome limitations from the sparse number of viewpoints and allow accurate reconstruction even when the underlying plasma distribution is complex or time-varying.
What would settle it
A direct comparison of the reconstructed CME velocity, propagation angles, and internal density structure against independent multi-viewpoint observations or in-situ measurements that shows errors much larger than the reported 3 percent and few-degree levels would indicate the constraints are insufficient.
Figures
read the original abstract
Coronagraphic observations enable direct monitoring of coronal mass ejections (CMEs) through scattered light from free electrons, but determining the 3D plasma distribution from 2D imaging data is challenging due to the optically-thin plasma and the complex image formation processes. We introduce SuNeRF-CME, a framework for 3D tomographic reconstructions of the heliosphere using multi-viewpoint coronagraphic observations. The method leverages Neural Radiance Fields (NeRFs) to estimate the electron density in the heliosphere through a ray-tracing approach, while accounting for the underlying Thomson scattering of image formation. The model is optimized by iteratively fitting the time-dependent observational data. In addition, we apply physical constraints in terms of continuity, propagation direction, and speed of the heliospheric plasma to overcome limitations imposed by the sparse number of viewpoints. We utilize synthetic observations of a CME simulation to fully quantify the model's performance for different viewpoint configurations. The results demonstrate that our method can reliably estimate the CME parameters from only two viewpoints, with a mean velocity error of $3.01\pm1.94\%$ and propagation direction errors of $3.39\pm1.94^\circ$ in latitude and $1.76\pm0.79^\circ$ in longitude. We further show that our approach can achieve a full 3D reconstruction of the simulated CME from two viewpoints, where we correctly model the three-part structure, deformed CME front, and internal plasma variations. Additional viewpoints can be seamlessly integrated, directly enhancing the reconstruction of the plasma distribution in the heliosphere. This study underscores the value of physics-informed methods for reconstructing the heliospheric plasma distribution, paving the way for unraveling the dynamic 3D structure of CMEs and enabling advanced space weather monitoring.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces SuNeRF-CME, a physics-informed Neural Radiance Fields framework for 3D tomographic reconstruction of coronal mass ejections from coronagraphic observations. It models electron density via ray-tracing that accounts for Thomson scattering, optimizes the NeRF by fitting time-dependent multi-view data, and augments the loss with physical regularizers enforcing continuity, constant propagation direction, and constant speed of heliospheric plasma. Quantitative evaluation on synthetic CME simulation data shows that two viewpoints suffice to recover CME parameters with a reported mean velocity error of 3.01±1.94 % and direction errors of 3.39±1.94° (latitude) and 1.76±0.79° (longitude), while also reproducing the three-part structure, deformed front, and internal density variations.
Significance. If the central claims hold under broader testing, the work would be significant for heliospheric imaging and space-weather applications: it demonstrates that NeRFs augmented with domain-specific plasma constraints can mitigate the severe under-determination inherent in optically-thin line-of-sight integrals when only two viewpoints are available. The explicit quantitative error metrics on synthetic data and the seamless extensibility to additional viewpoints are clear strengths.
major comments (2)
- [Abstract and §4] Abstract and §4 (synthetic validation): the reported velocity and direction errors are obtained exclusively on synthetic data generated from a simulation whose plasma evolution obeys the same constant-speed, constant-direction, and continuity assumptions used as regularizers. This matching of forward model and constraint set risks understating the tomographic ambiguity that would arise for real CMEs exhibiting variable speed, strong deformation, or internal flows; a load-bearing test would require at least one mismatched simulation or real-event case.
- [§3.2] §3.2 (physical constraints): the continuity, direction, and speed terms are introduced as weighted regularizers, yet no explicit functional form, weighting schedule, or ablation study is provided to show how these terms interact with the NeRF photometric loss to resolve the inverse problem. Without these details it is impossible to judge whether the constraints are sufficient to overcome the sparsity of two viewpoints or merely reproduce the simulation priors.
minor comments (2)
- [Figure captions and §4.2] Figure captions and §4.2 should explicitly state the number of synthetic viewpoints, the exact simulation code or reference, and the precise definition of the three-part structure used for visual comparison.
- [§2] Notation for the electron-density field and the Thomson-scattering integral should be introduced once in §2 and used consistently thereafter; several symbols appear without prior definition in the results section.
Simulated Author's Rebuttal
We are grateful to the referee for their thorough review and valuable suggestions. We respond to each major comment in turn and indicate the changes we will make to the manuscript.
read point-by-point responses
-
Referee: [Abstract and §4] Abstract and §4 (synthetic validation): the reported velocity and direction errors are obtained exclusively on synthetic data generated from a simulation whose plasma evolution obeys the same constant-speed, constant-direction, and continuity assumptions used as regularizers. This matching of forward model and constraint set risks understating the tomographic ambiguity that would arise for real CMEs exhibiting variable speed, strong deformation, or internal flows; a load-bearing test would require at least one mismatched simulation or real-event case.
Authors: We acknowledge that the synthetic CME data follows the same physical assumptions encoded in the regularizers. This was a deliberate choice to enable quantitative error assessment against known ground truth while isolating the effect of viewpoint sparsity. The regularizers are implemented as soft penalties rather than hard constraints, allowing the photometric loss to drive recovery of deformed fronts and internal density variations. In the revised manuscript we will expand the discussion in §4 to explicitly address this limitation, including a qualitative assessment of expected performance under violated assumptions and a statement of planned future tests on simulations with variable speed. revision: partial
-
Referee: [§3.2] §3.2 (physical constraints): the continuity, direction, and speed terms are introduced as weighted regularizers, yet no explicit functional form, weighting schedule, or ablation study is provided to show how these terms interact with the NeRF photometric loss to resolve the inverse problem. Without these details it is impossible to judge whether the constraints are sufficient to overcome the sparsity of two viewpoints or merely reproduce the simulation priors.
Authors: We agree that the current description in §3.2 lacks sufficient implementation detail. In the revised version we will add the explicit mathematical expressions for each regularizer term, the numerical weighting coefficients, and the schedule by which their relative strength is increased during optimization. We will also include a new ablation experiment that quantifies reconstruction error when each term is removed individually, thereby demonstrating their contribution to resolving the under-determined inverse problem with two viewpoints. revision: yes
Circularity Check
No significant circularity; reconstruction validated on independent synthetic benchmarks with standard physical regularizers.
full rationale
The paper optimizes a NeRF model via iterative fitting to time-dependent synthetic coronagraphic observations while adding physical constraints on plasma continuity, propagation direction, and speed as regularizers to address sparse viewpoints. Performance metrics (velocity error 3.01±1.94%, direction errors ~2-3°) are computed against ground-truth quantities from the CME simulation, which serves as an external benchmark rather than a self-derived input. No equations or steps reduce by construction to tautological definitions, fitted parameters renamed as predictions, or self-citation chains; the Thomson scattering ray-tracing and optimization are standard inverse-problem machinery applied to independent simulation data. The derivation chain remains self-contained against these benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- constraint weights for continuity, direction, and speed
axioms (2)
- domain assumption Image formation in coronagraphs is governed by Thomson scattering from free electrons in optically thin plasma
- domain assumption Heliospheric plasma obeys continuity and propagates with consistent direction and speed
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We apply physical constraints in terms of continuity, propagation direction, and speed of the heliospheric plasma... minimizing the residuals of the continuity equation... Lradial = ∥v × r∥ / (∥v∥ ∥r∥) ... Lvelocity bounds
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The model is optimized by iteratively fitting the time-dependent observational data... Physics-Informed Neural Radiance Field
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]
Amerstorfer, T., Hinterreiter, J., Reiss, M. A., et al. 2021, Space Weather, 19, e02553, doi: 10.1029/2020SW002553
-
[2]
2004, Journal of Atmospheric and Solar-Terrestrial Physics, 66, 1295
Arge, C., Luhmann, J., Odstrcil, D., Schrijver, C., & Li, Y. 2004, Journal of Atmospheric and Solar-Terrestrial Physics, 66, 1295
work page 2004
-
[3]
Aschwanden, M. J. 2011, Living Reviews in Solar Physics, 8, 5, doi: 10.12942/lrsp-2011-5 Asensio Ramos, A. 2023, SoPh, 298, 135, doi: 10.1007/s11207-023-02226-2 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy Collaboration, Price-Whelan, A. M., Sip˝ ocz, B. M., et al. 2018, AJ...
-
[4]
Barnes, W. T., Bobra, M. G., Christe, S. D., et al. 2020, The Astrophysical Journal, 890, 68
work page 2020
-
[5]
2021, Space Weather, 19, e02873, doi: 10.1029/2021SW002873
Bauer, M., Amerstorfer, T., Hinterreiter, J., et al. 2021, Space Weather, 19, e02873, doi: 10.1029/2021SW002873
-
[6]
Billings, D. E. 1966, A guide to the solar corona
work page 1966
-
[7]
Brueckner, G. E., Howard, R. A., Koomen, M. J., et al. 1995, SoPh, 162, 357, doi: 10.1007/BF00733434
-
[8]
Cyr, O. C., & Seiden, J. A. 2004, Journal of Geophysical Research (Space Physics), 109, A03103, doi: 10.1029/2003JA010149
-
[9]
Cargill, P. J. 2004, SoPh, 221, 135, doi: 10.1023/B:SOLA.0000033366.10725.a2
-
[10]
Colaninno, R. C., & Vourlidas, A. 2009, ApJ, 698, 852, doi: 10.1088/0004-637X/698/1/852 de Koning, C. A., & Pizzo, V. J. 2011, Space Weather, 9, 03001, doi: 10.1029/2010SW000595 de Koning, C. A., Pizzo, V. J., & Biesecker, D. A. 2009, SoPh, 256, 167, doi: 10.1007/s11207-009-9344-7
-
[11]
CPEG: A Convex Predictor-corrector Entry Guidance Algorithm,
Deforest, C., Killough, R., Gibson, S., et al. 2022, in 2022 IEEE Aerospace Conference, 1–11, doi: 10.1109/AERO53065.2022.9843340
-
[12]
DeForest, C. E., de Koning, C. A., & Elliott, H. A. 2017, ApJ, 850, 130, doi: 10.3847/1538-4357/aa94ca
-
[13]
DeForest, C. E., Howard, T. A., & Tappin, S. J. 2013, ApJ, 765, 44, doi: 10.1088/0004-637X/765/1/44 D´ ıaz Baso, C. J., Asensio Ramos, A., de la Cruz Rodr´ ıguez, J., da Silva Santos, J. M., & Rouppe van der Voort, L. 2025, A&A, 693, A170, doi: 10.1051/0004-6361/202452172 16
-
[14]
Fletcher, L., Dennis, B. R., Hudson, H. S., et al. 2011, SSRv, 159, 19, doi: 10.1007/s11214-010-9701-8
-
[15]
Forbes, T. G., Linker, J. A., Chen, J., et al. 2006, SSRv, 123, 251, doi: 10.1007/s11214-006-9019-8
-
[16]
Frazin, R. A., & Kamalabadi, F. 2005, ApJ, 628, 1061, doi: 10.1086/430846
-
[17]
Gibson, S. E., & Low, B. C. 1998, ApJ, 493, 460, doi: 10.1086/305107
-
[18]
Howard, T. A. 2015, Journal of Space Weather and Space Climate, 5, A22, doi: 10.1051/swsc/2015024
-
[19]
Howard, T. A., & Tappin, S. J. 2009, SSRv, 147, 31, doi: 10.1007/s11214-009-9542-5
-
[20]
2006, Journal of Geophysical Research (Space Physics), 111, A04S91, doi: 10.1029/2004JA010942
Webb, D. 2006, Journal of Geophysical Research (Space Physics), 111, A04S91, doi: 10.1029/2004JA010942
-
[21]
Jackson, B. V., Hick, P. P., Buffington, A., et al. 2011, Journal of Atmospheric and Solar-Terrestrial Physics, 73, 1214, doi: 10.1016/j.jastp.2010.10.007
-
[22]
E., Tremblay , B., Centeno , R., and Rempel , M
Rempel, M. 2025, ApJL, 985, L7, doi: 10.48550/arXiv.2502.13924
-
[23]
2023, Nature Astronomy, doi: 10.1038/s41550-023-02030-9
Podladchikova, T. 2023, Nature Astronomy, doi: 10.1038/s41550-023-02030-9
-
[24]
2024a, ApJL, 963, L21, doi: 10.3847/2041-8213/ad2450
Jarolim, R., Tremblay, B., Rempel, M., et al. 2024a, ApJL, 963, L21, doi: 10.3847/2041-8213/ad2450
-
[25]
2024b, ApJL, 961, L31, doi: 10.3847/2041-8213/ad12d2
Jarolim, R., Tremblay, B., Mu˜ noz-Jaramillo, A., et al. 2024b, ApJL, 961, L31, doi: 10.3847/2041-8213/ad12d2
-
[26]
Kaiser, M. L., Kucera, T. A., Davila, J. M., et al. 2008, SSRv, 136, 5, doi: 10.1007/s11214-007-9277-0
-
[27]
Keller, C. U. 2025, ApJ, 989, 92, doi: 10.3847/1538-4357/add6a9
-
[28]
Kraft, S., Puschmann, K. G., & Luntama, J. P. 2017, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 10562, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 105620F, doi: 10.1117/12.2296100
-
[29]
Mildenhall, B., Srinivasan, P. P., Tancik, M., et al. 2021, Communications of the ACM, 65, 99 M¨ ostl, C., Rollett, T., Frahm, R. A., et al. 2015, Nature Communications, 6, 7135, doi: 10.1038/ncomms8135
-
[30]
J., Christe, S., Freij, N., et al
Mumford, S. J., Christe, S., Freij, N., et al. 2020, SunPy, v1.1.4, Zenodo, doi: 10.5281/zenodo.3871057
-
[31]
2003, Advances in Space Research, 32, 497, doi: 10.1016/S0273-1177(03)00332-6
Odstrcil, D. 2003, Advances in Space Research, 32, 497, doi: 10.1016/S0273-1177(03)00332-6
-
[32]
2019, in Advances in Neural Information Processing Systems 32, ed
Paszke, A., Gross, S., Massa, F., et al. 2019, in Advances in Neural Information Processing Systems 32, ed. H. Wallach, H. Larochelle, A. Beygelzimer, F. d'Alch´ e-Buc, E. Fox, & R. Garnett (Curran
work page 2019
-
[33]
2018, Journal of Space Weather and Space Climate, 8, A35, doi: 10.1051/swsc/2018020
Pomoell, J., & Poedts, S. 2018, Journal of Space Weather and Space Climate, 8, A35, doi: 10.1051/swsc/2018020
-
[34]
Provornikova, E., Merkin, V. G., Vourlidas, A., et al. 2024, ApJ, 977, 106, doi: 10.3847/1538-4357/ad83b1
-
[35]
Raissi, M., Perdikaris, P., & Karniadakis, G. E. 2019, Journal of Computational Physics, 378, 686, doi: 10.1016/j.jcp.2018.10.045
-
[36]
2020, Advances in Neural Information Processing Systems, 33, 7462
Wetzstein, G. 2020, Advances in Neural Information Processing Systems, 33, 7462
work page 2020
-
[37]
2016, ApJ, 830, 58, doi: 10.3847/0004-637X/830/2/58
Susino, R., & Bemporad, A. 2016, ApJ, 830, 58, doi: 10.3847/0004-637X/830/2/58
-
[38]
2021, Living Reviews in Solar Physics, 18, 4, doi: 10.1007/s41116-021-00030-3
Temmer, M. 2021, Living Reviews in Solar Physics, 18, 4, doi: 10.1007/s41116-021-00030-3
-
[39]
Temmer, M., Preiss, S., & Veronig, A. M. 2009, SoPh, 256, 183, doi: 10.1007/s11207-009-9336-7
-
[40]
Temmer, M., Scolini, C., Richardson, I. G., et al. 2023, arXiv e-prints, arXiv:2308.04851, doi: 10.48550/arXiv.2308.04851
-
[41]
2011, ApJS, 194, 33, doi: 10.1088/0067-0049/194/2/33
Thernisien, A. 2011, ApJS, 194, 33, doi: 10.1088/0067-0049/194/2/33
-
[42]
Thernisien, A. F. R., Howard, R. A., & Vourlidas, A. 2006, ApJ, 652, 763, doi: 10.1086/508254 V´ asquez, A. M., Frazin, R. A., Hayashi, K., et al. 2008, ApJ, 682, 1328, doi: 10.1086/589682
-
[43]
Verbeke, C., Mays, M. L., Kay, C., et al. 2023, Advances in Space Research, 72, 5243, doi: 10.1016/j.asr.2022.08.056
-
[44]
Vourlidas, A., & Howard, R. A. 2006, ApJ, 642, 1216, doi: 10.1086/501122 Vrˇ snak, B., Sudar, D., Ruˇ zdjak, D., &ˇZic, T. 2007, A&A, 469, 339, doi: 10.1051/0004-6361:20077175 Vrˇ snak, B.,ˇZic, T., Vrbanec, D., et al. 2013, SoPh, 285, 295, doi: 10.1007/s11207-012-0035-4
-
[45]
Webb, D. F., & Howard, T. A. 2012, Living Reviews in Solar Physics, 9, 1
work page 2012
-
[46]
Wiegelmann, T., Thalmann, J. K., & Solanki, S. K. 2014, A&A Rv, 22, 78, doi: 10.1007/s00159-014-0078-7
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.