Pith. sign in

REVIEW 4 major objections 6 minor 8 references

Forecasting Shock-associated Energetic Particle Intensities in the Inner Heliosphere: A Proof-of-Concept Capability for the PUNCH Mission

T0 review · 4 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Shock speed jump can forecast solar particle intensities

desk verdict A useful proof-of-concept that consolidates a known ΔV–peak flux correlation and lays out a PUNCH forecasting workflow, but the forecasting pipeline depends on an unvalidated identification between imaging-derived and in-situ speed jumps. read the letter →

arxiv 2411.16510 v1 pith:R4RIIN2V submitted 2024-11-25 astro-ph.SR astro-ph.HE

Maher A. Dayeh , Michael J. Starkey , Heather A. Elliott , Raphael Attie , Craig E. DeForest , Radoslav Bučik , Mihir I. Desai This is my paper · ORCID
classification astro-ph.SRastro-ph.HE
keywords solarenergeticparticlesshockspeedjumpCME-drivenshockscorotatinginteractionregionsstormheliosphericimagingspaceweatherforecastingPUNCHmission
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper argues that the peak flux of shock-accelerated particles reaching Earth can be estimated from a single shock property: the jump in solar wind speed across the shock, $\Delta V$. Analyzing 59 CME-driven shocks and 74 corotating interaction region shocks observed by Wind between 1997 and 2023, the authors find a strong, statistically stable correlation between $\Delta V$ and the peak flux of associated energetic particles. They then propose that the upcoming PUNCH mission, which will track solar wind density structures in polarized white light, can measure the same $\Delta V$ remotely while the shock is still close to the Sun. That would turn the correlation into a probabilistic forecast window for particle intensities, with lead time gained from imaging rather than waiting for in situ arrival.

What carries the argument

The load-bearing quantity is $\Delta V$, the difference between upstream and downstream solar wind speed at the shock, obtained in situ from Wind plasma measurements. The argument's forward half is the empirical correlation of $\Delta V$ with peak flux; its forecasting half is Magnetic Balltracking, a technique that tracks small density structures in coronagraph images, adapted to estimate radial velocities and speed jumps of CMEs and CIRs from PUNCH polarized-light imagery. The paper treats the imaging-derived $\Delta V$ as a proxy for the in situ quantity, with an estimated velocity accuracy of 5% to 15% in 10-degree azimuthal bins based on STEREO/COR2 test data.

What would settle it

Compare imaging-derived shock speed jumps from PUNCH with in situ Wind or ACE measurements for the same interplanetary shocks; a systematic offset beyond the estimated 5% to 15% accuracy would falsify the forecasting claim, as would a failure of the $\Delta V$--peak-flux correlation to hold on an independent set of new events.

Watch

Extended reading notes

Core claim

The central discovery is an empirical power-law relation between the shock speed jump $\Delta V$ and the peak intensity of shock-accelerated energetic particles observed at 1 AU, built from 59 energetic storm particle events at CME-driven shocks and 74 particle-associated CIR shocks. The correlation is stable: leave-one-out cross-validation changes the correlation coefficient by at most 4% for ESPs and 8% for CIRs, and the fitted exponent is small enough that a given error in measuring $\Delta V$ produces a proportionally smaller error in predicted peak flux. On this basis the paper claims that the speed jump is a reliable probabilistic predictor of particle intensity, and that PUNCH's three-dimensional tracking of density structures can supply the needed $\Delta V$ early in the shock's propagation.

Load-bearing premise

The entire forecasting pipeline depends on the assumption that the shock speed jump PUNCH reconstructs from white-light images of density structures is the same physical quantity Wind measures in situ; if imaging-based $\Delta V$ is systematically biased relative to in situ $\Delta V$, the forecast fails even though the historical correlation is real.

Editorial extensions

If this is right

  • If the relation holds, PUNCH images could yield probabilistic estimates of ESP and CIR-associated particle peak fluxes while the disturbance is still in the inner heliosphere.
  • Forecasts can be updated as new images arrive, refining the speed jump estimate as the shock approaches Earth.
  • The small power-law exponent means that even rough speed measurements give usable flux estimates, making the method forgiving of imaging noise.
  • Because the same correlation holds for both CME-driven and CIR shocks, a single forecasting pipeline could cover two major particle-event classes.
  • The correlation provides a benchmark analytic form for future operational use, even without full physical modeling of particle acceleration.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If PUNCH imaging validates the $\Delta V$ proxy, the same pipeline could be extended to forecast heavy-ion composition or higher-energy channels, since the mechanism is shock-speed-driven.
  • Events whose fluxes deviate strongly from the $\Delta V$ relation may flag shocks where geometry (quasi-parallel versus quasi-perpendicular) or seed population matters more, offering a selection tool for case studies.
  • The claimed 5% to 15% velocity accuracy from COR2 may not transfer directly to PUNCH's different field of view and cadence; comparing image-derived $\Delta V$ to in situ Wind measurements for the same shocks would be a straightforward early test.
  • A probabilistic forecast window could be tuned into a threshold-based alert system for space weather operations, rather than a single deterministic prediction.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. The paper analyzes 59 CME-driven shocks (54 in Figure 1) and 74 CIR shocks observed by Wind/STEP between 1997 and 2023, correlating the observed helium peak flux in the 0.1–0.5 MeV range with the solar wind speed jump ΔV across the associated shock. The authors report a strong, stable power-law correlation for both ESP and CIR events, supported by a leave-one-out cross-validation check, and they propose that PUNCH imaging of leading-edge speeds and density ratios can provide an estimate of ΔV close to the Sun, enabling probabilistic forecasts of shock-associated particle intensities. The manuscript presents the empirical relation as the basis for a proof-of-concept forecasting pipeline.

Significance. If the empirical ΔV–peak-flux relation holds, it provides a useful, operationally motivated benchmark for space-weather forecasting, and the use of 25 years of Wind/STEP data with a LOOCV stability check is a strength. The paper is also honestly framed as a proof-of-concept rather than a validated forecast system, and the suggestion that PUNCH can cross-check shock speed jumps via two independent imaging methods (tracking and photometric density ratios) is a reasonable scientific direction. However, the forecasting claim currently rests on an uncalibrated identification of two different quantities—the imaging-derived leading-edge speed difference and the in-situ ΔV—and the statistical reporting omits the numerical correlation coefficients and fit uncertainties needed to assess predictive power.

major comments (4)
  1. [Section 3 and Figure 3] The forecasting pipeline equates the imaging-derived 'speed difference between the leading edge/shock and the background wind' with the in-situ ΔV used in the Figure 1 fits. For a shock with compression ratio r, the in-situ downstream-minus-upstream speed jump is ΔV = (1 - 1/r)(V_shock - V_upstream), whereas the imaging proxy is (V_shock - V_upstream); the two differ by a factor that varies with shock strength. The manuscript nowhere calibrates this mapping, nor does it demonstrate that a speed jump measured at 6–180 R_s predicts the 1-AU speed jump used to build the correlation. This is a load-bearing gap for the central forecasting claim. Please add a quantitative conversion (e.g., using an assumed or independently measured density ratio from the photometric method) and validate the near-Sun-to-1-AU propagation, or explicitly reframe the forecast as a testable hypothesis that will be validated with PUNCH data.
  2. [Abstract, Section 2, and Figure 1 caption] The event count for CME-driven shocks is inconsistent: the abstract and Section 2 state 59 ESP events, while the Figure 1 caption and the fitted scatter plots use 54 ESP events. The LOOCV description also refers to removing events from a sample of unspecified size. This inconsistency undermines the reproducibility of the analysis. Please reconcile the number and ensure the figure, text, and abstract agree.
  3. [Section 2 and Figure 1e–1h] The paper claims a 'strong' and 'robust' correlation that is 'statistically significant,' but it reports no numeric Pearson r, no p-values, no fit exponent, and no uncertainties on the power-law amplitude or exponent. The LOOCV variation is given only as a percentage range, which does not quantify the scatter or the confidence interval of the fitted relation. Without these numbers, the reader cannot judge the predictive value of the relation or compare it with prior studies. Please report the fitted parameters, their uncertainties, the correlation coefficients, and the p-values for each energy channel and for the combined samples.
  4. [Section 2, selection criteria] The coefficient-of-variation thresholds (σ_n<0.15, σ_V<0.38, σ_B<0.4) are explicitly described as 'chosen to optimize the number of events used and maintain the intended quiet conditions.' This threshold optimization can introduce selection bias that is not quantified. Please perform a sensitivity analysis (for example, varying each threshold by ±20% and reporting the resulting change in the fitted power-law parameters and in the event counts) so that the robustness of the correlation to the selection choices can be assessed.
minor comments (6)
  1. [Abstract] The phrase 'equipped with photometric that enables 3D tracking' appears to be missing a noun; presumably 'photometric imaging' or 'polarimetric imaging.'
  2. [Section 1] 'Whitman et al. 2022 and references there in' should be 'references therein.'
  3. [Section 2] 'to insure the ICME is not in a fast solar wind stream' should be 'to ensure'; also the subscript in avg(V_sw) is garbled in the text.
  4. [Section 2] In the definition of the flux enhancement for CIR events, the notation F_i and the energy pass band ΔE_j are not fully defined; please provide a clear equation with all symbols defined.
  5. [Figure 1 caption] 'Scattered plots' should be 'Scatter plots', and the caption refers to 'bottom panels' without labeling them (g,h); please add panel labels for clarity.
  6. [Section 3] The reference to 'Figure 1b,c' in the discussion of the peak-intensity relation should be 'Figure 1e,f,' matching the scatter plots.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the ΔV–peak-flux relation is recomputed directly from Wind/STEP data, and the PUNCH forecast step is an unvalidated lookup application rather than a derivation that relies on its own output.

full rationale

The paper's central empirical result is a fit between measured in-situ solar-wind speed jumps (ΔV) and Wind/STEP peak particle fluxes for 59 ESP and 74 CIR events. The inputs (shock speed jumps from the Helsinki shock database and Wind solar-wind data) and outputs (peak fluxes from STEP) are independent measurements; neither quantity is defined in terms of the other, so the correlation is not self-definitional. The forecasting proposal uses this fitted relation as a lookup step: a future PUNCH imaging-derived speed jump would be inserted into the empirical power law to obtain an estimated peak-flux range. That is an in-sample application of a previously fit relation to a new proxy input, not a fitted parameter renamed as a prediction. The paper does not validate that the imaging leading-edge-minus-background speed is quantitatively equal to the in-situ ΔV used in the fit, and the missing (1−1/r) compression-ratio conversion is a real validation gap, but this is an extrapolation/correctness concern rather than circular reasoning. The self-citations (Dayeh et al. 2018, Bučík et al. 2009, Attie & Innes 2015, DeForest et al. 2018) supply prior context and algorithm provenance; the central correlation is recomputed here with LOOCV stability checks and therefore does not reduce to the authors' earlier claims. No derivation step equates the predicted quantity with an input by construction.

Assumptions & free parameters 5 free parameters · 5 assumptions · 0 invented entities

No new physical entities are postulated. The analysis relies on fitted power-law parameters, hand-chosen selection thresholds, and the assumption that remote imaging will recover the same speed jumps measured in situ. The free parameters are numerous because the paper reports no coefficients and uses several ad hoc data-selection choices to build the event sample.

free parameters (5)
  • Power-law amplitude and exponent for peak flux vs Delta V fits = not reported
    The weighted fits in Figure 1 are fitted to the 1997-2023 Wind/STEP data; the coefficients are not given in the text.
  • Quiet-time coefficient-of-variation thresholds (sigma_n < 0.15, sigma_V < 0.38, sigma_B < 0.4) = 0.15, 0.38, 0.4
    Chosen to optimize the number of ESP events while maintaining quiet pre-shock conditions; different thresholds would change the sample and likely the correlation.
  • CME association windows (12-hour after shock, 24-hour before shock) = 12 h and 24 h
    Ad hoc association window choices that define which shocks are classified as CME-driven.
  • CIR flux enhancement window (3-hour centered on reverse boundary) = 3 h
    Window used to define a particle enhancement for CIR events; affects which events enter the sample.
  • Background wind speed cut (avg(V_sw) < 500 km/s) = 500 km/s
    Used to exclude ICMEs in fast wind streams; boundary chosen by hand.
assumptions (5)
  • standard math Pearson correlation and leave-one-out cross-validation are appropriate for measuring and validating the peak-flux vs Delta V relation.
    These are standard statistical tools used without proof in Section 2.
  • domain assumption Wind/STEP particle measurements, the Helsinki shock database, the Richardson-Cane ICME list, and the Broiles CIR list are reliable inputs.
    The analysis depends on these external catalogs and instruments being accurate.
  • domain assumption The energetic particle enhancements in the selected events are causally associated with the identified CME and CIR shocks.
    Selection criteria try to enforce this, but contamination by other particle sources remains possible.
  • domain assumption The solar wind speed jump measured in situ equals the shock speed jump that would be observed remotely.
    This equivalence is required for the PUNCH forecasting concept but is not demonstrated.
  • ad hoc to paper The Magnetic Balltracking velocity accuracy of 5-15% measured on STEREO/COR2 will be similar for PUNCH.
    The paper states an expectation of similar accuracy for PUNCH without direct evidence.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Forecasting Shock-associated Energetic Particle Intensities in the Inner Heliosphere: A Proof-of-Concept Capability for the PUNCH Mission." pith.science (2026). https://pith.science/paper/R4RIIN2V

@misc{pith2026241116510,
  author       = {Pith},
  title        = {Pith review of: Forecasting Shock-associated Energetic Particle Intensities in the Inner Heliosphere: A Proof-of-Concept Capability for the PUNCH Mission},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/R4RIIN2V}},
  note         = {Machine review of arXiv:2411.16510}
}
read the original abstract

Solar energetic particles (SEPs) associated with shocks driven by fast coronal mass ejections (CMEs) or shocks developed by corotating interaction regions (CIRs) often extend to high energies, and are thus key elements of space weather. The PUNCH mission, set to be launched in 2025, is equipped with photometric that enables 3D tracking of solar wind structures in the interplanetary space through polarized light. Tracking techniques are used to estimate speeds and speed gradients of solar structures, including speed jumps at fast shocks. We report on a strong and a robust relation between the shock speed jump magnitude at CME and CIR shocks and the peak fluxes of associated energetic particles from the analysis of 59 CME-driven shocks and 74 CIRs observed by Wind/STEP between 1997-2023. We demonstrate that this relation, along with PUNCH anticipated observations of solar structures can be used to forecast shock-associated particle events close to the Sun; thus, advancing and providing a crucial input to forecasting of SEP fluxes in the heliosphere.

Figures

Figures reproduced from arXiv: 2411.16510 by the authors.

Figure 1
Figure 1. (a) Temporal profiles of H intensities between 0.1-0.5 MeV in an ESP and a CIR (b) event, along with the solar wind profiles (c,d). (e,f) Scattered plots showing the correlation and the weighted fits between observed peak fluxes and the speed jump in 54 ESP events and 74 CIRs. Bottom panels demonstrate how the correlation varies if one event is removed, and colored boxes and arrows indicate the events affecting the … view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

8 extracted references · 8 canonical work pages

  1. [1]

    Accepted in Solar Physics Forecasting Shock-associated Energetic Particle Intensities in the Inner Heliosphere: A Proof-of-Concept Capability for the PUNCH Mission M. A. Dayeh1,2, M. J. Starkey1, H. A. Elliott1,2, R. Attie4,5, C. E. DeForest3, R. Bučik1, and M. I. Desai1,2 1 Southwest Research Institute, San Antonio, TX 78238 (maldayeh@swri.edu) 2 Univers...

    work page 2025

  2. [2]

    background corona

    Introduction Energetic particles in interplanetary space play a critical role in space weather, posing significant challenges for deep space exploration (Whitman et al. 2022 and references there in). These particles originate from various sources within and beyond our heliosphere, spanning solar and galactic sources (Desai and Giacalone, 2016). Galactic c...

    work page 2022

  3. [4]

    Magnetic Balltracking

    (a) Temporal profiles of H intensities between 0.1-0.5 MeV in an ESP and a CIR (b) event, along with the solar wind profiles (c,d). (e,f) Scattered plots showing the correlation and the weighted fits between observed peak fluxes and the speed jump in 54 ESP events and 74 CIRs. Bottom panels demonstrate how the correlation varies if one event is removed, a...

    work page 2015

  4. [5]

    References Attie, R., & Innes, D. E. 2015, Astron Astrophys, 574, A106 Benz, A. O., & Güdel, M. 2010, ARA&A, 48, 241 Bučík, R., Space Sci Rev 216, 24,

    work page 2015

  5. [1995]

    We utilize the Helsinki shock database (http://ipshocks.fi), an ICME list (Richardson and Cane 2010), and a CIR list (Broiles et al

    onboard the Wind spacecraft. We utilize the Helsinki shock database (http://ipshocks.fi), an ICME list (Richardson and Cane 2010), and a CIR list (Broiles et al

    work page 2010

  6. [2009]

    A., Desai, M

    Dayeh, M. A., Desai, M. I., Ebert, R. W., et al. 2018, JPhCS, 1100, 012008 Chancellor, J. C., Scott, G. B. I., & Sutton, J. P. 2014, Life, 4, 491 Decker, R. B. 1981, J Geophys Res: Space Phys, 86, 4537 DeForest, C. E., Howard, R. A., Velli, M., Viall, N., & Vourlidas, A. 2018, Astrophys J, 862, 18 DeForest, C., Killough, R., Gibson, S., et al. 2022, 2022 ...

    work page 2018

  7. [2012]

    Rev., 173(1–4), 247–281

    Space Sci. Rev., 173(1–4), 247–281. Mason, G. M., Leske, R. A., Desai, M. I., et al. 2008, Astrophysical J, 678, 1458 Maurer, R. H., Fretz, K., Angert, M. P., et al. 2017, ITNS, 64, 2782 Moreland, K., Dayeh, M. A., Li, G., et al. 2023, Astrophys J, 956, 107 Onorato, G., Di Schiavi, E., & Di Cunto, F. 2020, FrP, 8, 1 Richardson, I. G., & Cane, H. V . 2010,...

    work page 2008

  8. [2020]

    W., Desai, M

    Broiles, T. W., Desai, M. I., Lee, C. O., & MacNeice, P. J. 2013, J Geophys Res Space Phys, 118, 4776 Bučík, R., Mall, U., Korth, A., and Mason, G. M., Ann. Geophys., 27, 3677–3690,

    work page 2013

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.