{"id":"244a5374-bb0f-4ab9-af24-7719cc6a6663","arxiv_id":"2411.16510","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Peak fluxes of shock-accelerated particles at Earth correlate with shock speed jumps, and the PUNCH mission could turn that correlation into a forecasting tool.","lead":"The paper analyzes decades of solar wind and particle data and finds that the size of a shock's speed jump correlates with the peak intensity of the energetic particles it accelerates. It argues that the upcoming PUNCH mission can measure those speed jumps from images and use the correlation for earlier space weather forecasts.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The PUNCH forecast pipeline has not been validated: imaging-derived speed jumps are not the same quantity as the in-situ ΔV used in the correlation, and no test shows they can predict the 1 AU peak flux.","rationale":"The paper's empirical correlation between in-situ ΔV and peak particle flux may well be real, and the LOOCV check is a useful stability test. The reader correctly identified the PUNCH velocity-accuracy premise as a weak point. My concern is sharper: even if the 5-15% tracking accuracy transfers, the observable that PUNCH provides is not the in-situ ΔV used in the fit. The conversion from shock-speed-minus-background-wind to downstream-minus-upstream speed jump requires the compression ratio, and that conversion is neither derived nor validated. Moreover, the timing mismatch between near-Sun imaging and the 1 AU in-situ correlation is unaddressed. These are concrete, load-bearing gaps in the forecasting claim, but they do not invalidate the underlying correlation; they just show the proof-of-concept is incomplete. The reader's CONDITIONAL verdict already captures this, so no verdict change is needed. My proposed dress rehearsal would turn the condition into a concrete validation step.","tokens_in":5147,"tokens_out":4338,"duration_ms":46239,"concrete_test":"Perform a historical PUNCH dress rehearsal: for the 59 ESP events (or a subset with COR2/HI1 coverage), reconstruct PUNCH-like ΔV by tracking the leading edge and background wind in white-light images at radii 30-120 R_s, optionally corrected by photometric density ratios, and compare these values with Wind in-situ ΔV for the same shocks. Then apply the paper's fitted power law to predict peak He flux at 1 AU and compare with observations. If imaging-derived ΔV differs systematically from in-situ ΔV (notably by the missing (1-1/r) factor) or if the predicted-peak-flux rank correlation is not significantly better than a simple solar-wind-speed climatology, the PUNCH forecast proof-of-concept is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central forecasting claim rests on identifying the speed jump PUNCH will measure remotely with the ΔV that enters the empirical correlation. Section 3 states that 'the speed difference between the leading edge/shock and the background wind is proportional to the speed jump,' but this is not the in-situ ΔV used in Figure 1. For a shock with compression ratio r, the in-situ downstream-minus-upstream solar wind speed jump is ΔV = (1 - 1/r)(V_shock - V_upstream), whereas the leading-edge-minus-background proxy is V_shock - V_upstream; the two differ by the factor (1 - 1/r). This is not a 5-15% tracking error; it is a systematic conversion that depends on shock strength and is not calibrated anywhere in the paper. The photometric density-ratio method mentioned in Section 2 could supply r, but no quantitative relation is given between imaging-based ΔV and the Wind-based ΔV used to fit the power law. In addition, all events in the correlation were measured at 1 AU after shock arrival, while PUNCH would estimate ΔV at 6-180 R_s before arrival; the paper does not demonstrate that a near-Sun speed jump predicts the 1 AU speed jump used in the correlation. Without this two-step validation, the proof-of-concept pipeline is an unsupported extrapolation of a plausible empirical correlation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":5487,"tokens_out":2535,"duration_ms":26619,"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":[{"comment":"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.","section":"Section 3 and Figure 3"},{"comment":"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.","section":"Abstract, Section 2, and Figure 1 caption"},{"comment":"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.","section":"Section 2 and Figure 1e–1h"},{"comment":"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.","section":"Section 2, selection criteria"}],"minor_comments":[{"comment":"The phrase 'equipped with photometric that enables 3D tracking' appears to be missing a noun; presumably 'photometric imaging' or 'polarimetric imaging.'","section":"Abstract"},{"comment":"'Whitman et al. 2022 and references there in' should be 'references therein.'","section":"Section 1"},{"comment":"'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.","section":"Section 2"},{"comment":"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.","section":"Section 2"},{"comment":"'Scattered plots' should be 'Scatter plots', and the caption refers to 'bottom panels' without labeling them (g,h); please add panel labels for clarity.","section":"Figure 1 caption"},{"comment":"The reference to 'Figure 1b,c' in the discussion of the peak-intensity relation should be 'Figure 1e,f,' matching the scatter plots.","section":"Section 3"}],"recommendation":"major_revision","confidential_remarks":"The empirical correlation is a reasonable and potentially useful result, but the forecasting claim—the stated proof-of-concept for PUNCH—is currently supported only by an uncalibrated equivalence between imaging-derived and in-situ speed jumps. The authors should be encouraged to add a quantitative calibration or reframe the claim as a future validation exercise. The event-count inconsistency and missing statistical details are straightforward to fix and should be addressed in revision. The paper fits the scope of Solar Physics well."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is worth a read. It reproduces the established correlation between shock speed jump ΔV and energetic particle peak flux on a larger combined sample of ESPs and CIRs, and it adds a concrete workflow for how PUNCH could use imaging to estimate ΔV near the Sun. The leave-one-out cross-validation is a responsible check, and the authors are honest that the relation provides a probabilistic forecast window rather than precise per-event predictions. That is real value.\n\nThe soft spots are in the forecasting application, not the underlying correlation. First, the paper never reports the correlation coefficients, p-values, or fit uncertainties, even though it mentions Pearson and statistical significance. The reader also caught that the abstract says 59 ESP events while Figure 1e shows 54; that inconsistency needs fixing. The selection thresholds are tuned to maximize event numbers, which is fine but should be disclosed more explicitly.\n\nThe bigger issue, which the stress-test note gets right, is that the PUNCH pipeline identifies the speed difference between the CME leading edge and the background wind with the ΔV used in the correlation. For a shock, the in-situ ΔV is (1 − 1/r) times the leading-edge speed difference, where r is the compression ratio. The paper says this is “proportional” but never calibrates the factor or shows that a near-Sun imaging proxy predicts the 1 AU ΔV used to fit the power law. Without that validation, the forecast pipeline is an extrapolation of a plausible empirical relation. This does not sink the paper, because the authors frame it as a proof-of-concept, but it is the load-bearing assumption that determines whether the forecast actually works.\n\nWho should read this: space weather forecasters, the PUNCH team, and researchers working on SEP and CIR empirical correlations. The paper deserves a serious referee. It is coherent, useful, and the empirical core is solid, but the forecasting claim needs major revision: add the missing statistics, reconcile the event counts, and most importantly, address the conversion between imaging-based and in-situ ΔV, either with a calibration or an explicit validation plan.","headline":"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.","tokens_in":6052,"tokens_out":1929,"would_cite":true,"duration_ms":21662,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Shock speed jump can forecast solar particle intensities","keywords":["solar energetic particles","shock speed jump","CME-driven shocks","corotating interaction regions","energetic storm particles","heliospheric imaging","space weather forecasting","PUNCH mission"],"falsifier":"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.","tokens_in":73,"feed_emoji":"🌞","tokens_out":8588,"duration_ms":126553,"temperature":0.7,"pith_summary":"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.","feed_headline":"Shock speed jump predicts particle storm intensity","feed_subtitle":"Wind data tie peak particle flux to a shock's speed jump; PUNCH imaging could measure that jump days before arrival.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Showed that speed jump correlates with ESP peak fluxes, the relation this paper extends to more events.","marker":"(Dayeh et al. 2018)"},{"why":"Established speed jump as a correlative parameter for CIR-associated particle enhancements.","marker":"(Bučík et al. 2009)"},{"why":"Introduced Magnetic Balltracking, the tracking method adapted here for imaging-based speed estimates.","marker":"(Attie & Innes 2015)"},{"why":"Describes the PUNCH mission and its polarized-light capability that motivates the forecasting concept.","marker":"(DeForest et al., 2022)"},{"why":"Describes the Wind/STEP instrument that provides the energetic particle data for the correlation.","marker":"(von Rosenvinge et al. 1995)"},{"why":"Supplies the ICME list used to select the 59 ESP events.","marker":"(Richardson and Cane 2010)"},{"why":"Supplies the CIR list used to select the 74 CIR events.","marker":"(Broiles et al. 2013)"},{"why":"Provides the leave-one-out cross-validation method showing the correlation is not driven by single events.","marker":"(Hastie et al. 2009)"}],"fun_headline_variants":["Shock speed jump predicts solar particle storm intensity","PUNCH mission to forecast particle storms via shock jumps","Shock jump size foretells energetic particle peak flux","New forecast method: shock speed jump ties to particle flux","Speed jump at shocks predicts space weather particle storms"],"cache_read_input_tokens":8064,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Shock speed jump predicts solar particle storm intensity","PUNCH mission to forecast particle storms via shock jumps","Shock jump size foretells energetic particle peak flux","New forecast method: shock speed jump ties to particle flux","Speed jump at shocks predicts space weather particle storms"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000305,"raw_usage":{"total_tokens":1735,"prompt_tokens":912,"completion_tokens":823,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":528,"completion_tokens_details":{"reasoning_tokens":746}},"tokens_in":528,"tokens_out":823,"duration_ms":8349,"temperature":1.0,"reasoning_tokens":746,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:02:01.339362+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"W., Desai, M","cited_arxiv_id":null,"evidence_quote":"Supplies the CIR list used to select the 74 CIR events."},{"cited_title":"A., Desai, M","cited_arxiv_id":null,"evidence_quote":"Showed that speed jump correlates with ESP peak fluxes, the relation this paper extends to more events."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduced Magnetic Balltracking, the tracking method adapted here for imaging-based speed estimates."}],"review_version":1}