{"id":"8e80f6a1-b23f-420d-857e-a1080688ac98","arxiv_id":"2412.11852","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":16,"one_line_summary":"PARASOL couples SOLPACS-calibrated foreshock physics to PARADISE transport, and reproduces the 12 July 2012 ESP event intensity to within an order of magnitude for energies below about 5 MeV.","lead":"A new model called PARASOL combines a fast, semi-analytical description of particle acceleration at shock waves with a 3D particle transport simulation, using solar wind and shock parameters from MHD models. In a test on the 12 July 2012 solar energetic particle event, it reproduced the measured energetic storm particle intensities near the shock within about a factor of ten.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Agreement rests on three free parameters labeled 'illustrative' in §3.3; without a priori constraints or a second event, the single-event reproduction does not establish forecast skill.","rationale":"The reader's weakest assumption correctly identifies the free parameters as the most load-bearing point. The paper's own text supports this: Section 3.3 calls the simulation an 'illustrative example' and the three parameters 'free', Section 4 states that the high-energy overestimate is 'largely' due to the assumed Δt and that Δt is 'not necessarily constant', and Section 5 and 3.3 defer constraint of ϵ_inj and Δt to future multi-event studies. The reproduction within one order of magnitude is a genuine forward-modeling result for that parameter choice, and the authors are transparent about the model's limitations, including the inverted low-energy upstream spectra. However, the central implication for forecasting—that PARASOL can predict ESP intensities without in-situ SEP observations—requires that the parameters be set a priori. A single event with chosen parameters does not establish this. The proposed parameter scan would directly quantify how much of the agreement is due to parameter selection; if a large fraction of plausible parameter space yields agreement, the concern is mitigated, whereas a narrow tuned region would confirm that forecast skill is not yet shown. This concern does not invalidate the model-development contribution, but it does justify the conditional verdict: the model is promising and the reproduction is evidence of capability, yet forecasting skill remains unverified. The 1 au mean free path assumption (Eqs. 8, 21) is a secondary limitation that also favours a conditional rather than a full acceptance, but the free-parameter issue is the primary load-bearing concern. Overall, the reader's CONDITIONAL verdict with high confidence is appropriate and does not require revision.","tokens_in":22949,"tokens_out":6035,"duration_ms":57197,"concrete_test":"Perform a systematic parameter scan for the July 12, 2012 event over physically plausible ranges: ϵ_inj ∈ [1e-4, 1e-3], Δt ∈ [1 h, 100 h], x_M ∈ [5 R_sun, 15 R_sun]. For each combination, run the PARASOL chain and record the ratio of the peak simulated intensity at 1–5 MeV at Earth to the observed value. If the set of parameter combinations that reproduces the observation within a factor of 10 is a small or hand-picked subset of the scanned volume, then the claimed agreement is parameter-tuned rather than an a priori forecast. Ideally, also fix the parameters from a first event (or from SOLPACS physics) and run a second, unseen ESP event without re-tuning; a second event outside one order of magnitude would decisively show that forecast skill is not yet demonstrated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of reproducing the July 12, 2012 ESP event within one order of magnitude is contingent on the choices ϵ_inj = 5e-4, Δt = 10 h, and x_M = 10 R_sun, which Section 3.3 explicitly calls 'free parameters' and 'an illustrative example'. Section 4 states that Δt is 'not necessarily constant' and that the overestimated high-energy intensities are 'largely' due to the assumed Δt = 10 h, while Figure 11 shows that reducing ϵ_inj by a factor of five lowers intensities across all channels and weakens the ESP peak. Thus the agreement is not a robust model prediction but a result of selecting parameter values that are not fixed by any a priori rule. The paper itself defers constraining these parameters to future multi-event studies (§3.3, §5). Without such constraints, the headline claim that PARASOL can produce a 2-hour ESP forecast without using in-situ SEP observations is unsupported; the single-event reproduction is a demonstration of model capability, not of forecasting skill. A secondary but related issue is the use of the SOLPACS-calibrated 1 au mean free path λ^(S)_1au for shocks at other distances (Eqs. 8, 21), which is an approximation that further ties the result to the 1 au test geometry.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces PARASOL, a physics-based model for solar energetic particle (SEP) forecasting that couples a semi-analytical description of the inner foreshock, calibrated against SOLPACS self-consistent simulations of shock acceleration, with the PARADISE test-particle transport model. The model is driven by MHD parameters from EUHFORIA and is tested on the 12 July 2012 SEP event, with the central claim that the simulated energetic storm particle (ESP) component (E ≲ 5 MeV) near the shock is reproduced within one order of magnitude in intensity. The paper also discusses discrepancies in the event onset and in the upstream low-energy spectrum, and outlines future work to constrain the model's three free parameters (injection efficiency ε_inj, acceleration timescale Δt, and matching distance x_M).","tokens_in":23379,"tokens_out":3701,"duration_ms":36422,"significance":"If the claimed reproduction is robust, PARASOL represents a valuable step toward operational, physics-based ESP forecasting because it does not use in-situ SEP observations of the event to set parameters, and the full chain can run in roughly two hours on a supercomputer. The model's construction from SOLPACS simulations of self-consistent wave-particle interactions is physically well motivated, and the use of an independent EUHFORIA MHD run for the 2012 event gives the validation a forward-modeling character. However, the current single-event demonstration with hand-set free parameters does not yet establish forecast skill, and the paper's honest acknowledgment of this limitation does not fully remove the need to qualify the headline claim.","major_comments":[{"comment":"The three free parameters ε_inj = 5×10⁻⁴, Δt = 10 h, and x_M = 10 R_sun are not fixed by any a priori rule; Section 3.3 explicitly labels the simulation an illustrative example, and Section 4 states that Δt is 'not necessarily constant' and that the overestimated high-energy intensities are 'largely' due to the assumed Δt = 10 h. Figure 11 shows that reducing ε_inj to 10⁻⁴ lowers intensities across all channels and weakens the ESP peak. Because the agreement therefore depends on parameter choices that are not independently constrained, the statement that the model 'reproduced' the observed ESP event should be presented as a capability demonstration, not as validation of forecast skill, and the paper should either provide a sensitivity analysis or explicitly restrict the claim to the illustrative parameter set.","section":"Section 3.3 and Eq. (27)"},{"comment":"The simulated upstream low-energy spectrum is inverted (intensity increasing with energy), while the observed spectrum eight hours before shock arrival is nearly flat and the observed time-intensity profiles at low energies overlap. The paper itself calls this a 'striking discrepancy' and attributes it to too-efficient trapping of low-energy particles near the shock. This is a substantive failure to reproduce the pre-shock low-energy enhancement, which is part of the ESP event, and it should be reflected in the scope of the headline claim: the 'within one order of magnitude' agreement applies to the shock-peak intensities, not to the full ESP foreshock behavior.","section":"Section 4 and Fig. 10b"},{"comment":"The PARADISE parallel mean free path λ^(P) is taken from the SOLPACS-calibrated 1 au fit and used for shocks at other heliocentric distances along the shock surface. This is an extrapolation: the SOLPACS scaling properties in Appendix A show that the physical mean free path depends on the local ambient parameters (density, magnetic field, Mach number), so the 1 au calibration is not obviously valid for inner heliosphere shock portions, especially when particles are injected along the entire shock front. The manuscript should justify this approximation with a scaling argument or with a test at a different heliocentric distance, as it directly affects the demodulation in Eq. (7) and the source spectrum in Eq. (20).","section":"Eqs. (8) and (21)"},{"comment":"The semi-analytical model relies on a large set of fitted coefficients (K, β_λ, E_b, δ_λ, C_1, C_2, ΔE_1, ΔE_2, α_1, α_2, and the cutoff-energy fit) derived from only fourteen SOLPACS runs, with no uncertainty quantification or goodness-of-fit statistics reported. For example, C_2 is taken as an average over runs with ~20% scatter, and α_2 is set to a constant. While such fits are acceptable for a first model version, the absence of uncertainty information makes it difficult to assess how much of the 2012 agreement is attributable to the fitted parameter values rather than to the model physics; at minimum, the fit residuals should be shown for the key dependencies.","section":"Section 2.2.4 and Eqs. (31)-(36)"}],"minor_comments":[{"comment":"The phrase 'within one order magnitude' should read 'within one order of magnitude', and the sentence in Section 4 beginning 'A test PARASOL simulation... have reproduced' has a subject-verb agreement error ('simulation... have' should be 'simulation... has').","section":"Abstract and Section 4"},{"comment":"The title contains a stray space in 'P ARASOL' in the manuscript header; the paper body otherwise uses 'PARASOL' consistently.","section":"Section 1 and Fig. 1"},{"comment":"The notation x_M is introduced as the matching distance and later fixed to 10 R_sun in Section 3.3, but Section 2.1 discusses x_M as energy-dependent and states that x_M > 10 R_sun for the quasi-steady-state energies; this tension should be clarified, for instance by explaining how the illustrative fixed value is chosen relative to the energy-dependent criterion in Eq. (11).","section":"Section 2.1"},{"comment":"Figure 10b is discussed only briefly; the reader would benefit from an explicit statement of the time offset relative to shock arrival and from labeling the energy channels used to define 'flat' versus 'inverted' spectra, since these definitions are central to the discussion in Section 4.","section":"Section 3.1 and Fig. 10"},{"comment":"The downstream mean free path in Eq. (44) contains an exponential damping scale d_0 = 1 R_sun, but no physical justification or sensitivity test is given for this choice; a one-sentence rationale or a reference would improve reproducibility.","section":"Section 2.3"}],"recommendation":"major_revision","confidential_remarks":"The paper is transparent about its limitations and the self-citations provide the simulation tools rather than forcing the validation result, so there is no circularity problem. The main issue for the journal is whether the 'reproduced the observed ESP event' claim is appropriately qualified given the hand-set free parameters and the unvalidated 1 au mean-free-path extrapolation. If the authors rephrase the headline claim as a capability demonstration, show sensitivity to the free parameters, and clarify the scope of the pre-shock discrepancy, the paper would be suitable for publication. I would not reject it, but the current framing overstates what a single illustrative event can establish."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. The paper does what it claims: it presents the first version of PARASOL, coupling a SOLPACS-calibrated semi-analytical inner foreshock with PARADISE transport, and shows a forward simulation of the July 12, 2012 ESP event within an order of magnitude. The fitting functions (Eqs. 17-36) are new and grounded in the SOLPACS simulations; the authors show the fits and give parameter dependencies. The validation is genuinely forward in that no in-situ SEP data from that event are used to set the parameters. The paper is also unusually transparent: it calls the three parameters 'free' and labels the run 'an illustrative example', and it discusses the inverted low-energy spectrum at length. Credit where due.\n\nThe soft spots are real but not fatal. The three free parameters - injection efficiency, acceleration timescale, and matching distance - are not constrained a priori. Delta_t is the most problematic because it directly controls the cutoff energy, and Section 4 admits the overestimated high-energy intensities are 'largely' due to the assumed Delta_t = 10 h. The sensitivity run with reduced epsilon_inj shows how much the result moves. So the agreement with the 2012 event is not yet a model prediction in the forecasting sense; it is a demonstration that the machinery can produce realistic ESP intensities with reasonable parameter choices. The paper says this itself and defers the parameter study to future work. That is a legitimate way to publish a first model, provided the abstract doesn't oversell it. It mostly doesn't.\n\nI also note the inverted upstream spectrum is a genuine unresolved discrepancy, and the 1 au SOLPACS mean free path is being applied to shocks at other distances, which is an approximation. Neither sinks the paper; both are clearly flagged. The citation pattern is fine: the self-citations point to the actual tools and the earlier foreshock model, not to pad the reference list.\n\nVerdict: send it to referees. A good referee should ask for a sensitivity study of the three parameters and ideally a second event to begin justifying the forecasting claim. But this is a useful, honest building block, better than most first model papers.","headline":"A solid, transparent first-model paper: worth refereeing, but the single-event agreement is a capability demonstration, not forecast skill, given three unconstrained free parameters.","tokens_in":23898,"tokens_out":2292,"would_cite":true,"duration_ms":21478,"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":"The paper introduces PARASOL, a physics-based forecasting chain for energetic storm particle events, and reports that a test simulation of the 12 July 2012 event reproduced observed intensities within one order of magnitude.","keywords":["solar energetic particles","energetic storm particles","CME-driven shocks","diffusive shock acceleration","space weather forecasting","particle transport","pitch-angle scattering","self-generated Alfvén waves"],"falsifier":"Run PARASOL with the fixed parameter set $\\epsilon_{\\rm inj}=5\\times10^{-4}$, $\\Delta t=10$ h, $x_M=10\\,R_\\odot$ on several well-observed ESP events at different heliocentric distances; if any event's peak proton intensity at $\\lesssim 5$ MeV near the shock deviates from observations by more than one order of magnitude, or if matching required changing the fixed parameters, the central forecasting claim is refuted.","tokens_in":22747,"feed_emoji":"☀️","tokens_out":15338,"duration_ms":117602,"temperature":0.7,"pith_summary":"PARASOL is a new modeling chain aimed at forecasting energetic storm particle (ESP) events: the sharp spike in low-energy proton intensity when a CME-driven shock sweeps past a spacecraft. Its central idea is to split the problem. A semi-analytical description of the inner foreshock — the region just upstream of the shock where self-generated Alfvén waves dominate scattering — is fitted once to detailed kinetic simulations and then supplies the accelerated-particle spectrum to a fast test-particle transport code for the outer heliosphere. The chain requires only magnetohydrodynamic (MHD) solar wind and shock parameters as input, not in-situ SEP observations of the event being forecast. For the 12 July 2012 event, the paper reports that a test simulation reproduced the observed ESP intensity at energies $\\lesssim 5$ MeV in the close vicinity of the shock within one order of magnitude, and that the full chain runs in about two hours on a supercomputer. If that accuracy holds with fixed calibration parameters, physics-based near-real-time ESP forecasting is within reach.","feed_headline":"PARASOL reproduces the July 2012 proton storm within 10x intensity","feed_subtitle":"The model needs only MHD solar wind inputs and runs in about two hours — a step toward operational ESP forecasting.","key_machinery":"The load-bearing object is the semi-analytical inner-foreshock model. From fourteen SOLPACS simulation runs spanning shock Alfvénic Mach number, Alfvén speed, shock-normal angle, and injection efficiency, the paper fits the upstream proton mean free path as $\\lambda^{(S)}(x,E) = \\Lambda(E)\\,[1 + (x/\\Delta x(E))^{q(E)}]^{1/q(E)}$, where $\\Lambda(E)$ is a modified Bell-theory power law with an exponential rollover and $\\Delta x(E)$, $q(E)$ are themselves analytical functions of shock parameters and of the rollover energy $E_b$. This mean free path is inserted into the steady-state Parker equation to connect the SOLPACS shock spectrum to a matching distance $x_M$; the spectrum is then demodulated and remodulated with the 1 au-calibrated mean free path used in PARADISE, yielding the shock emission spectrum injected into the focused transport equation. The cutoff energy of the injected spectrum depends on the free acceleration timescale $\\Delta t$ through Eq. (27), and the full scheme has three free parameters: injection efficiency $\\epsilon_{\\rm inj}$, the timescale $\\Delta t$, and the matching distance $x_M$.","core_discovery":"The central claim is that a semi-analytical foreshock model, calibrated once on self-consistent simulations of coupled proton acceleration and Alfvén wave generation, can serve as the shock source for a transport code and reproduce a real ESP event without event-specific tuning to observed particle intensities. Concretely, the test simulation of the 12 July 2012 event, driven by EUHFORIA MHD output, matched the observed proton intensity at energies up to about 5 MeV near the shock within one order of magnitude, and produced an ESP peak at the modeled shock arrival as observed. The model also reproduces the observed low-energy spectral slope at the shock while overestimating higher-energy intensities, which the paper attributes to the assumed acceleration timescale $\\Delta t = 10$ h. Along the shock surface the model predicts a spatial pattern in which low-energy protons concentrate at the quasi-parallel eastern flank and high-energy protons at the quasi-perpendicular nose and western flank, controlled by the local shock obliquity and Alfvénic Mach number.","pith_inferences":["A direct consequence the authors leave for future work is a multi-event test: run PARASOL with $\\epsilon_{\\rm inj}=5\\times10^{-4}$, $\\Delta t=10$ h, and $x_M=10\\,R_\\odot$ on several well-observed ESP events and measure the scatter in peak intensity; the paper's forecasting claim stands or falls on that scatter.","The model's predicted inverted (intensity increasing with energy) upstream spectrum at low energies is a distinctive signature of self-generated turbulence; observed flat or overlapping upstream spectra would imply that the mean free path's energy dependence near the shock is not simply increasing with energy, suggesting missing physics such as Alfvén wave damping or distance-dependent resonance b","Basing the outer-foreshock mean free path on a SOLPACS calibration at 1 au, together with the $\\alpha^{1/2}\\epsilon_{\\rm inj}=\\mathrm{const.}$ scaling property, implies a definite prediction for how foreshock intensities and spectra change with heliocentric distance; radially separated spacecraft observations could test that scaling directly."],"forward_implications":["If the fixed parameter choices hold across events, ESP peak intensities near the shock can be forecast from MHD solar wind and shock inputs before the shock arrives, without using in-situ SEP observations of the event.","The complete chain consumes about 1500 CPU-hours, about two hours of wall-clock time on a supercomputer, which the paper regards as adequate for timely ESP prediction.","Because the injected particle population is derived from MHD shock properties rather than by accelerating particles inside a finite-resolution MHD shock, the model avoids the overly thick shock that softened high-energy spectra in the earlier PARADISE-only simulation of this event.","The model makes spatial predictions along the shock: the quasi-parallel eastern flank emits the highest low-energy intensities, while the quasi-perpendicular nose and west flank emit the highest high-energy intensities, reflecting the local shock geometry and Alfvénic Mach number.","Adding a coronal MHD model should extend the same chain from ESP events toward full SEP event forecasting, including the early onset phase that the present inner boundary at 21.5 solar radii cannot capture."],"supporting_citations":[{"why":"Supplies the SOLPACS model and the full-resonance pitch-angle scattering formalism behind the simulated mean free paths that the semi-analytical foreshock is fitted to.","marker":"Afanasiev et al. (2015)"},{"why":"Validates SOLPACS against an observed ESP event and provides the calibrated parameterizations and scaling property used to build the fitting grid.","marker":"Afanasiev et al. (2023)"},{"why":"Gives the earlier semi-analytical foreshock model and the Bell-theory parameterizations that PARASOL modifies for the injection and cutoff-energy formulas.","marker":"Vainio et al. (2014)"},{"why":"Provides the steady-state shock acceleration theory with self-generated Alfvén waves that underlies the analytical mean free path forms.","marker":"Bell (1978)"},{"why":"Supplies the steady-state Parker equation and acceleration-rate theory used in the demodulation from the SOLPACS spectrum to the matching distance.","marker":"Drury (1983)"},{"why":"Defines the EUHFORIA MHD model that produces the solar wind and shock parameters driving the test simulation.","marker":"Pomoell and Poedts (2018)"},{"why":"Provides the EUHFORIA simulation of the 12 July 2012 event and the earlier PARADISE-only ESP modeling baseline that PARASOL is compared against.","marker":"Wijsen et al. (2022)"},{"why":"Supplies the spheromak CME parameters used to set up the EUHFORIA simulation of the 12 July 2012 event.","marker":"Scolini et al. (2019)"},{"why":"Defines the PARADISE focused transport model that carries particles through the outer foreshock and ambient solar wind.","marker":"Wijsen (2020)"}],"fun_headline_variants":["PARASOL forecast nails July 2012 proton storm within 10x","New model predicts SEP events from MHD inputs alone","PARASOL: 2-hour SEP forecast matches real event","Solar storm forecasting gets a boost with PARASOL model","PARASOL reproduces 2012 SEP event, no tuning needed"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model's forecasting claim rests on the three free parameters — $\\epsilon_{\\rm inj}=5\\times10^{-4}$, $\\Delta t=10$ h, and $x_M=10\\,R_\\odot$ — being set in advance rather than tuned to the event being forecast, and on the SOLPACS-calibrated 1 au mean free path describing scattering near shocks at other distances.","fun_headline_variants_meta":{"raw":{"variants":["PARASOL forecast nails July 2012 proton storm within 10x","New model predicts SEP events from MHD inputs alone","PARASOL: 2-hour SEP forecast matches real event","Solar storm forecasting gets a boost with PARASOL model","PARASOL reproduces 2012 SEP event, no tuning needed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000463,"raw_usage":{"total_tokens":2364,"prompt_tokens":1045,"completion_tokens":1319,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":661,"completion_tokens_details":{"reasoning_tokens":1231}},"tokens_in":661,"tokens_out":1319,"duration_ms":12341,"temperature":1.0,"reasoning_tokens":1231,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:31:40.123120+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run PARASOL with the fixed parameter set $\\epsilon_{\\rm inj}=5\\times10^{-4}$, $\\Delta t=10$ h, $x_M=10\\,R_\\odot$ on several well-observed ESP events at different heliocentric distances; if any event's peak proton intensity at $\\lesssim 5$ MeV near the shock deviates from observations by more than one order of magnitude, or if matching required changing the fixed parameters, the central forecasting claim is refuted.","supporting_citations":[{"cited_title":"Self-consistent Monte Carlo simulations of proton acceleration in coronal shocks: Effect of anisotropic pitch-angle scattering of particles","cited_arxiv_id":"1603.08857","evidence_quote":"Supplies the SOLPACS model and the full-resonance pitch-angle scattering formalism behind the simulated mean free paths that the semi-analytical foreshock is fitted to."},{"cited_title":"Vainio , D","cited_arxiv_id":null,"evidence_quote":"Validates SOLPACS against an observed ESP event and provides the calibrated parameterizations and scaling property used to build the fitting grid."},{"cited_title":"P \\\"o nni , M","cited_arxiv_id":null,"evidence_quote":"Gives the earlier semi-analytical foreshock model and the Bell-theory parameterizations that PARASOL modifies for the injection and cutoff-energy formulas."},{"cited_title":"Observation-based modelling of the energetic storm particle event of 14 July 2012","cited_arxiv_id":"2201.06454","evidence_quote":"Provides the EUHFORIA simulation of the 12 July 2012 event and the earlier PARADISE-only ESP modeling baseline that PARASOL is compared against."},{"cited_title":"Observation-based modelling of magnetised Coronal Mass Ejections with EUHFORIA","cited_arxiv_id":"1904.07059","evidence_quote":"Supplies the spheromak CME parameters used to set up the EUHFORIA simulation of the 12 July 2012 event."}],"review_version":1}