arxiv: 2605.04284 · v1 · submitted 2026-05-05 · 🌌 astro-ph.SR

Recognition: unknown

Ensemble modeling of Coronal Mass Ejection dynamics and forecasts at 1 AU with a semi-analytic flux-rope model

S. Stamkos , S. Patsourakos , A. Vourlidas , E. Paouris

Authors on Pith no claims yet

Pith reviewed 2026-05-08 17:06 UTC · model grok-4.3

classification 🌌 astro-ph.SR

keywords coronal mass ejectionsflux rope modelensemble modelingspace weather forecastinguncertainty propagationsolar wind interaction1 AU propagationMonte Carlo simulation

0 comments

The pith

Uncertainties in erupting flux-rope inputs produce 2-8 hour spreads in CME arrival times at 1 AU

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper runs Monte Carlo ensembles of a semi-analytic erupting flux rope model to trace how uncertainties in CME eruption parameters and solar wind conditions affect forecasts of arrival time, speed, magnetic field strength, and impact duration at 1 AU. Across six events with plus or minus 20 percent input sampling, it measures the resulting spreads and identifies the dominant controlling parameters such as poloidal-flux injection history for timing and background flow for speed. A sympathetic reader would care because the work shows which physical quantities most limit forecast precision, indicating where improved measurements would tighten predictions of when and how a CME will hit Earth. The model includes updated sheath mass and drag to better capture propagation forces.

Core claim

The ensembles reveal event-dependent dispersions where, for 20 percent input variations, time-of-arrival spreads reach 2.4 to 7.7 hours controlled mainly by poloidal-flux injection history, upstream wind speed, and drag coefficient. Leading-edge speed spreads are 28 to 53 km/s driven by background flow. Sheath magnetic fields spread 1 to 3.5 nT while internal flux-rope fields spread 1 to 7.6 nT, and impact durations spread 2.4 to 6.3 hours governed by geometric size and expansion scaling.

What carries the argument

The semi-analytic erupting flux rope (EFR) model with updated sheath mass and drag term, embedded in a Monte Carlo framework with truncated-normal sampling of key eruption and background solar-wind inputs

If this is right

Poloidal-flux injection history, upstream wind speed, and drag coefficient are the main controls on time-of-arrival uncertainty at 1 AU.
Background solar-wind flow properties dominate uncertainty in leading-edge speed forecasts.
Sheath magnetic field predictions remain relatively tight compared with internal flux-rope field predictions.
Geometric size and expansion scaling most strongly affect the duration of CME impact at Earth.
Eruption-driving and flux-content parameters limit the precision of internal magnetic field forecasts.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Tighter observational constraints on poloidal flux evolution near the Sun would narrow arrival-time prediction windows.
The identified sensitivities point to priority targets for future solar missions measuring near-Sun CME properties.
The same ensemble method applied to a larger event sample could test whether the uncertainty patterns generalize across different CME types.
Separate handling of sheath versus internal field uncertainties suggests distinct observation strategies for each diagnostic.

Load-bearing premise

The updated semi-analytic erupting flux rope model accurately represents the dominant forces and geometry from eruption through 1 AU for the sampled events.

What would settle it

In-situ observations at 1 AU for the six events that show arrival-time or magnetic-field spreads outside the model's 1-sigma ensemble ranges would falsify the claimed quantitative links between inputs and forecast spreads.

Figures

Figures reproduced from arXiv: 2605.04284 by A. Vourlidas, E. Paouris, S. Patsourakos, S. Stamkos.

**Figure 1.** Figure 1: Leading-edge evolution of the CME flux rope as a func view at source ↗

**Figure 2.** Figure 2: Endpoint spreads at 1 AU for the comprehensive ensemble of the 23 February 1997 event, plotted relative to the standard view at source ↗

**Figure 3.** Figure 3: Top five input parameters ranked by the Pearson correlation–based F score for each 1 AU diagnostic in the comprehensive view at source ↗

read the original abstract

This study quantifies how uncertainty in physically meaningful coronal mass ejection (CME) and solar-wind inputs propagates into forecast-relevant diagnostics from eruption to 1 AU. We use a semi-analytic erupting flux rope (EFR) model to simulate CME initiation and Sun-to-1 AU propagation under Lorentz, gravitational, and drag forces, driven by a prescribed time-dependent poloidal-flux injection. Relative to the original EFR formulation, we include sheath and pile-up effects through an effective mass and update the drag term for CME solar-wind coupling. The model is embedded in a Monte Carlo framework with truncated-normal sampling of key eruption and background solar-wind inputs. Across six CME events, the ensembles show event-dependent dispersion in the 1 AU diagnostics. For +/- 20% input sampling, all spreads are 1-sigma ensemble standard deviations. The time-of-arrival spread is 2.4-7.7 h and is mainly controlled by the poloidal-flux injection history, upstream wind speed, and drag coefficient. The leading-edge speed spread is 28-53 km/s and is primarily controlled by background-flow properties. Magnetic-field diagnostics show two regimes: the sheath field is relatively tightly distributed, with a spread of 1-3.5 nT and sensitivity to upstream wind, size, and expansion scaling, whereas the internal flux-rope field has a larger spread of 1-7.6 nT and is governed mainly by eruption-driving and flux-content parameters. The impact-duration spread is 2.4-6.3 h and is controlled mostly by geometric size and expansion scaling, with additional sensitivity to the driving timescale. These results establish a quantitative link between EFR input uncertainties and the resulting spread in CME arrival and impact diagnostics, identifying the physical parameters that most strongly limit forecast precision at 1 AU.

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.

Desk Editor's Note private letter to a colleague

This paper extends a semi-analytic flux-rope model with sheath mass and drag updates, then runs Monte Carlo ensembles on six events to rank input sensitivities for 1 AU spreads, but provides no comparison of its predictions to actual observations.

read the letter

The core result is a sensitivity study: across six CMEs, the model shows time-of-arrival spreads of 2.4-7.7 hours driven mainly by poloidal-flux history, upstream wind speed, and drag coefficient, with speed spreads of 28-53 km/s tied to background flow and magnetic-field spreads varying by sheath versus internal rope. The work adds effective mass for pile-up and revises the drag term relative to the base EFR setup, then samples inputs with truncated normals at +/-20% to produce these numbers. That gives a concrete mapping from input uncertainty to forecast diagnostics, which is useful for thinking about where better observations would help most. The presentation is straightforward and the parameter rankings are reported clearly for each diagnostic. It is incremental rather than foundational, but the exercise is done transparently on real events. The main limitation is the missing validation step. The abstract and results focus on internal ensemble spreads without showing how the median or mean predictions line up against measured arrival times, speeds, or in-situ B-fields for the same six events. That leaves open whether the dominant controls identified (flux for arrival, wind for speed) reflect real forecast error budgets or just the model's internal dynamics. The sampling range itself is stated without much justification against observed input uncertainties. This is the kind of paper that fits a space-weather modeling group that already uses or tests semi-analytic codes. It is worth sending to peer review because the method is reproducible, the question is operationally relevant, and the gaps are fixable with added comparisons or model-fidelity checks rather than fatal. Reviewers will likely press on the validation point, but the underlying sensitivity analysis is worth the time.

Referee Report

2 major / 2 minor

Summary. The manuscript presents an ensemble modeling study of six CME events using an updated semi-analytic erupting flux rope (EFR) model that incorporates sheath/pile-up effects via effective mass and a revised drag term. Uncertainties in eruption parameters (e.g., poloidal-flux injection history) and background solar-wind properties are sampled via truncated-normal distributions in a Monte Carlo framework; the resulting spreads in 1 AU diagnostics (ToA: 2.4-7.7 h, leading-edge speed: 28-53 km/s, sheath B: 1-3.5 nT, internal B: 1-7.6 nT, impact duration: 2.4-6.3 h) are reported and attributed to specific controlling parameters such as poloidal flux for arrival time and background flow for speed.

Significance. If the underlying EFR model is shown to be faithful, the work supplies a useful quantitative mapping from input uncertainties to forecast-relevant output spreads and highlights which physical quantities most constrain precision at 1 AU. The forward Monte Carlo sampling of physically motivated parameters and the separation of sensitivity regimes for sheath versus flux-rope fields are clear strengths that could inform targeted observational campaigns and model development.

major comments (2)

[Abstract and Results] The central claim that the ensembles identify parameters that 'most strongly limit forecast precision at 1 AU' rests on the premise that the updated EFR model accurately reproduces the dominant dynamics for the sampled events. However, the manuscript provides no comparison of ensemble means or medians against in-situ 1 AU arrival times, speeds, or magnetic-field measurements for the six events, nor any residual analysis. Without this anchor, the reported sensitivity rankings remain intra-model results whose relevance to actual forecast error budgets is untested (Abstract; Results section describing the six events).
[Methods] The 1-sigma ensemble standard deviations are presented as the primary uncertainty measures, yet no uncertainties on these statistics themselves (e.g., via bootstrap or jackknife estimates) are supplied, and the truncated-normal sampling with fixed ±20% bounds is adopted without demonstration that the resulting input distributions are consistent with observed variability in eruption parameters or solar-wind properties (Methods section on Monte Carlo framework).

minor comments (2)

[Abstract] Clarify in the text whether the reported spreads are computed identically for all diagnostics or whether any event-specific weighting or filtering is applied before quoting the 2.4-7.7 h ToA range.
[Model description] The description of the drag-term update and effective-mass formulation would benefit from an explicit equation reference or brief derivation to allow readers to assess how the changes differ from the original EFR implementation.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed review of our manuscript. We address each major comment below and have revised the manuscript accordingly to improve clarity and statistical rigor while preserving the core focus of the study on uncertainty propagation within the EFR model framework.

read point-by-point responses

Referee: [Abstract and Results] The central claim that the ensembles identify parameters that 'most strongly limit forecast precision at 1 AU' rests on the premise that the updated EFR model accurately reproduces the dominant dynamics for the sampled events. However, the manuscript provides no comparison of ensemble means or medians against in-situ 1 AU arrival times, speeds, or magnetic-field measurements for the six events, nor any residual analysis. Without this anchor, the reported sensitivity rankings remain intra-model results whose relevance to actual forecast error budgets is untested (Abstract; Results section describing the six events).

Authors: The primary objective of this work is a sensitivity analysis that quantifies how uncertainties in physically motivated input parameters propagate through the updated EFR model to produce spreads in 1 AU diagnostics. We agree that the manuscript does not include direct comparisons of ensemble statistics to in-situ observations for these six events, as such validation lies outside the stated scope. The sensitivity rankings are therefore intra-model results, but they remain useful for identifying which parameters most influence forecast-relevant outputs and for guiding targeted observations and model refinements. Prior publications have validated the base EFR model against observations; we will revise the Abstract and Results to explicitly state the intra-model nature of the analysis, clarify that the reported spreads assume model fidelity, and add references to those validation studies. This addresses the concern without expanding the manuscript into a full validation study. revision: partial
Referee: [Methods] The 1-sigma ensemble standard deviations are presented as the primary uncertainty measures, yet no uncertainties on these statistics themselves (e.g., via bootstrap or jackknife estimates) are supplied, and the truncated-normal sampling with fixed ±20% bounds is adopted without demonstration that the resulting input distributions are consistent with observed variability in eruption parameters or solar-wind properties (Methods section on Monte Carlo framework).

Authors: We agree that reporting uncertainties on the ensemble standard deviations would strengthen the statistical presentation. In the revised Methods section we will add bootstrap resampling (with 1000 resamples) to provide 1-sigma uncertainties on each reported spread. For the choice of truncated-normal distributions with ±20% bounds, this range was selected to reflect typical observational uncertainties cited in the literature for CME eruption parameters (e.g., poloidal flux, size) and solar-wind properties. We will expand the Methods text to include this justification together with supporting references to observational studies that document the variability ranges used. revision: yes

Circularity Check

0 steps flagged

No circularity: ensemble spreads derived from forward sampling of inputs through model equations

full rationale

The paper's central results are obtained by Monte Carlo forward propagation of physically motivated input uncertainties (e.g., poloidal-flux injection, upstream wind speed, drag coefficient) through the updated semi-analytic EFR model equations, yielding 1 AU diagnostic spreads. No output is defined in terms of itself, no parameter is fitted to a subset of results and then relabeled as a prediction, and no load-bearing premise reduces to a self-citation chain or ansatz. The model modifications (effective mass for sheath/pile-up, updated drag) are presented as explicit updates to a prior formulation and applied uniformly; the sensitivity rankings follow directly from the ensemble statistics without circular reduction. This is a standard self-contained sensitivity study.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The model rests on standard MHD force balance plus two ad-hoc updates (effective sheath mass and revised drag) whose functional forms are not derived from first principles in the abstract; input distributions are chosen rather than measured.

free parameters (3)

poloidal-flux injection history
Time-dependent driving term sampled with +/-20% uncertainty; directly controls arrival-time spread.
drag coefficient
Updated coupling term between CME and solar wind; listed as a primary controller of arrival-time spread.
background solar-wind speed and density
Upstream flow properties sampled; control leading-edge speed and sheath-field spreads.

axioms (2)

domain assumption Lorentz, gravitational, and aerodynamic drag forces dominate CME propagation from Sun to 1 AU
Invoked to justify the semi-analytic force balance in the EFR model.
ad hoc to paper Truncated-normal distributions adequately represent real uncertainties in eruption and solar-wind parameters
Sampling method chosen without reference to measured error distributions for the six events.

pith-pipeline@v0.9.0 · 5665 in / 1528 out tokens · 38589 ms · 2026-05-08T17:06:28.719996+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

55 extracted references

[1]

A., et al

Amerstorfer, T., Hinterreiter, J., Reiss, M. A., et al. 2021, Space Weather, 19, e02553

2021
[2]

2018, Space Weather, 16, 784

Amerstorfer, T., Möstl, C., Hess, P., et al. 2018, Space Weather, 16, 784

2018
[3]

K., DeV ore, C

Antiochos, S. K., DeV ore, C. R., & Klimchuk, J. A. 1999, ApJ, 510, 485

1999
[4]

N., Li, X., Pulkkinen, A., et al

Baker, D. N., Li, X., Pulkkinen, A., et al. 2013, Space Weather, 11, 585

2013
[5]

Bobra, M. G. & Couvidat, S. 2015, ApJ, 798, 135

2015
[6]

Cargill, P. J. 2004, Sol. Phys., 221, 135

2004
[7]

D., Biesecker, D

Cash, M. D., Biesecker, D. A., Pizzo, V ., et al. 2015, Space Weather, 13, 611

2015
[8]

Chen, J. 1996, J. Geophys. Res., 101, 27499

1996
[9]

& Garren, D

Chen, J. & Garren, D. A. 1993, Geophys. Res. Lett., 20, 2319

1993
[10]

& Krall, J

Chen, J. & Krall, J. 2003, Journal of Geophysical Research (Space Physics), 108, 1410

2003
[11]

& Kunkel, V

Chen, J. & Kunkel, V . 2010, ApJ, 717, 1105

2010
[12]

Chen, J., Marqué, C., V ourlidas, A., Krall, J., & Schuck, P. W. 2006, ApJ, 649, 452

2006
[13]

B., Del Moro, D., & Erdélyi, R

Chierichini, S., Liu, J., Korsós, M. B., Del Moro, D., & Erdélyi, R. 2024, ApJ, 963, 121

2024
[14]

Crowe, C. T. 2011, Multiphase Flows with Droplets and Particles, Second Edi- tion, 2nd edn. (CRC, Boca Raton FL, USA) Dumbovi´c, M., ˇCalogovi´c, J., Vršnak, B., et al. 2018, ApJ, 854, 180

2011
[15]

P., Biffis, E., Hapgood, M

Eastwood, J. P., Biffis, E., Hapgood, M. A., et al. 2017, Risk Analysis, 37, 206

2017
[16]

& Gibson, S

Fan, Y . & Gibson, S. E. 2003, ApJ, 589, L105

2003
[17]

Forbes, T. G. & Isenberg, P. A. 1991, ApJ, 373, 294

1991
[18]

Illing, R. M. E. & Hundhausen, A. J. 1986, J. Geophys. Res., 91, 10951

1986
[19]

J., & Palmerio, E

Kay, C., Nieves-Chinchilla, T., Hofmeister, S. J., & Palmerio, E. 2022, Space Weather, 20, e2022SW003165

2022
[20]

2024, Space Weather, 22, e2024SW003951

Kay, C., Palmerio, E., Riley, P., et al. 2024, Space Weather, 22, e2024SW003951

2024
[21]

D., Lynch, B

Kazachenko, M. D., Lynch, B. J., Welsch, B. T., & Sun, X. 2017, ApJ, 845, 49

2017
[22]

Kilpua, E., Koskinen, H. E. J., & Pulkkinen, T. I. 2017, Living Reviews in Solar Physics, 14, 5

2017
[23]

& Török, T

Kliem, B. & Török, T. 2006, Phys. Rev. Lett., 96, 255002

2006
[24]

T., Howard, R

Krall, J., Chen, J., Duffin, R. T., Howard, R. A., & Thompson, B. J. 2001, ApJ, 562, 1045

2001
[25]

2000, ApJ, 539, 964

Krall, J., Chen, J., & Santoro, R. 2000, ApJ, 539, 964

2000
[26]

2012, PhD thesis, George Mason University, Virginia

Kunkel, V . 2012, PhD thesis, George Mason University, Virginia

2012
[27]

& Chen, J

Kunkel, V . & Chen, J. 2010, ApJ, 715, L80

2010
[28]

& Chen, J

Lin, C.-H. & Chen, J. 2022, Journal of Geophysical Research (Space Physics), 127, e28744

2022
[29]

2018, ApJ, 855, 109

Liu, J., Ye, Y ., Shen, C., Wang, Y ., & Erdélyi, R. 2018, ApJ, 855, 109

2018
[30]

D., Luhmann, J

Liu, Y . D., Luhmann, J. G., Kajdiˇc, P., et al. 2014, Nature Communications, 5, 3481

2014
[31]

Manchester, W., Kilpua, E. K. J., Liu, Y . D., et al. 2017, Space Sci. Rev., 212, 1159

2017
[32]

B., Gombosi, T

Manchester, W. B., Gombosi, T. I., Roussev, I., et al. 2004, Journal of Geophys- ical Research (Space Physics), 109, A02107

2004
[33]

L., Taktakishvili, A., Pulkkinen, A., et al

Mays, M. L., Taktakishvili, A., Pulkkinen, A., et al. 2015, Sol. Phys., 290, 1775 Möstl, C., Isavnin, A., Boakes, P. D., et al. 2017, Space Weather, 15, 955

2015
[34]

2022, Space Weather, 20, e2021SW002925

Napoletano, G., Foldes, R., Camporeale, E., et al. 2022, Space Weather, 20, e2021SW002925

2022
[35]

2018, Journal of Space Weather and Space Climate, 8, A11

Napoletano, G., Forte, R., Del Moro, D., et al. 2018, Journal of Space Weather and Space Climate, 8, A11

2018
[36]

Odstrcil, D., Riley, P., & Zhao, X. P. 2004, Journal of Geophysical Research (Space Physics), 109, A02116

2004
[37]

& V ourlidas, A

Paouris, E. & V ourlidas, A. 2022, Space Weather, 20, e2022SW003070

2022
[38]

2015, Space Weather, 13, 734

Pulkkinen, A. 2015, Space Weather, 13, 734

2015
[39]

A., & Yurchyshyn, V

Qiu, J., Hu, Q., Howard, T. A., & Yurchyshyn, V . B. 2007, ApJ, 659, 758

2007
[40]

L., Andries, J., et al

Riley, P., Mays, M. L., Andries, J., et al. 2018, Space Weather, 16, 1245

2018
[41]

2024, A&A, 689, A187

Rodriguez, L., Shukhobodskaia, D., Niemela, A., et al. 2024, A&A, 689, A187

2024
[42]

P., Poirier, N., Lavarra, M., et al

Rouillard, A. P., Poirier, N., Lavarra, M., et al. 2020, ApJS, 246, 72

2020
[43]

I., Gombosi, T

Roussev, I. I., Gombosi, T. I., Sokolov, I. V ., et al. 2003, ApJ, 595, L57

2003
[44]

J., Kauristie, K., Aylward, A

Schrijver, C. J., Kauristie, K., Aylward, A. D., et al. 2015, Advances in Space Research, 55, 2745

2015
[45]

Schuck, P. W. 2010, ApJ, 714, 68

2010
[46]

Shafranov, V . D. 1966, Reviews of Plasma Physics, 2, 103

1966
[47]

2011, ApJ, 743, 101 Tóth, G., Sokolov, I

Temmer, M., Rollett, T., Möstl, C., et al. 2011, ApJ, 743, 101 Tóth, G., Sokolov, I. V ., Gombosi, T. I., et al. 2005, Journal of Geophysical Research (Space Physics), 110, A12226

2011
[48]

M., Thalmann, J

Tschernitz, J., Veronig, A. M., Thalmann, J. K., Hinterreiter, J., & Pötzi, W. 2018, ApJ, 853, 41 V ourlidas, A., Patsourakos, S., & Savani, N. P. 2019, Philosophical Transactions of the Royal Society of London Series A, 377, 20180096 Vršnak, B. 2021, J. Space Weather Space Clim., 11, 34 Vršnak, B., Žic, T., Vrbanec, D., et al. 2013, Sol. Phys., 285, 295

2018
[49]

Webb, D. F. & Howard, T. A. 2012, Living Reviews in Solar Physics, 9, 3

2012
[50]

2011, Fluid Mechanics (McGraw Hill, New York)

White, F. 2011, Fluid Mechanics (McGraw Hill, New York)

2011
[51]

& Sakurai, T

Wiegelmann, T. & Sakurai, T. 2021, Living Reviews in Solar Physics, 18, 1

2021
[52]

M., Mays, M

Wold, A. M., Mays, M. L., Taktakishvili, A., et al. 2018, Journal of Space Weather and Space Climate, 8, A17

2018
[53]

P., Howard, R

Zhang, J., Dere, K. P., Howard, R. A., Kundu, M. R., & White, S. M. 2001, ApJ, 559, 452

2001
[54]

G., Webb, D

Zhang, J., Richardson, I. G., Webb, D. F., et al. 2007, Journal of Geophysical Research (Space Physics), 112, A10102

2007
[55]

Zhu, C., Qiu, J., Liewer, P., et al. 2020, ApJ, 893, 141 Article number, page 12 of 20 S.Stamkos et al.: Ensemble Modeling of CME Dynamics Appendix A: 30 April 1997 event: comprehensive ensemble results This appendix presents the Monte Carlo ensemble results for the 30 April 1997 event. We show (i) the kinematic evolution of the flux-rope leading-edge spe...

2020