{"id":"4c7c0da6-6a94-4bfd-a8b2-696a889fe268","arxiv_id":"2502.00175","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A constant-turn flux rope MHD simulation reproduces the rotating magnetic field of a stealth ICME observed at radially aligned Solar Orbiter and Earth, with +/-5 hour arrival time errors.","lead":"This paper simulates a stealth coronal mass ejection with a data-constrained flux rope model and compares the magnetic field and arrival time at Solar Orbiter and Earth. The model reproduces the rotating magnetic field signature at both spacecraft, with arrival time errors of about five hours in opposite directions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Assumed poloidal flux (10^22 Mx) sets B0 and thereby the magnetic field magnitude; the 'very good accuracy' claim is not an independent test because the flux is unconstrained for this stealth CME.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the unconstrained poloidal flux sets the magnetic field amplitude, so the agreement in field magnitude is not an independent validation. My reading of the paper confirms this is the central soft spot. The paper is honest about the assumption, but the abstract and conclusion nevertheless assert 'very good accuracy' and the ability to capture radial evolution, which overstates the strength of the evidence. The model has other weaknesses—the selection of WSA realization R006, the ad-hoc energy definition, the 5-hour time shift at SolO, and the acknowledged cancellation of insertion-time and solar-wind-speed errors—but these primarily affect arrival timing and plasma properties, not the magnetic field magnitude claim. The magnetic field rotation shape is a genuine prediction from the GCS geometry and the assumed helicity, and the rare radial alignment is a valuable testbed. However, the field amplitude agreement is essentially a single-parameter fit to an unconstrained quantity, so the paper's strongest claim should be tempered. The conditional verdict is appropriate; no further downgrade is warranted because the paper explicitly discloses the assumption and the qualitative rotation pattern is still a meaningful success. The proposed test—fitting the flux at SolO and predicting Earth—would directly settle whether the radial-evolution claim has predictive content beyond the chosen flux value.","tokens_in":14940,"tokens_out":5509,"duration_ms":57971,"concrete_test":"Re-run the simulation with the poloidal flux inferred by fitting B0 to the observed SolO in situ magnetic field data (e.g., minimizing the RMS difference between simulated and observed BT/BN inside the ejecta at SolO). Then compare the predicted magnetic field profile at Earth to the observations without further tuning. If the Earth field amplitude and rotation are reproduced, the radial-evolution claim is robust. If not, the original agreement was contingent on the arbitrary flux assumption. Additionally, a sensitivity scan over poloidal flux values in the range 0.3–3×10^22 Mx, reporting the peak B and RMS error at Earth, would quantify how strongly the 'very good accuracy' depends on the assumed flux.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that the CTFR model reproduces the rotating magnetic field 'with very good accuracy' at both SolO and Earth—depends critically on the assumed poloidal flux of 10×10^21 Mx stated in Section 3.2. This value sets B0 in the Vandas & Romashets (2017) field solution, so the simulated BT and BN amplitudes scale linearly with it. The paper explicitly acknowledges that 'observations needed to constrain the magnetic flux of the flux rope are not available' and that the value was assumed. Given that the magnetic field profile is produced by inserting a fully formed flux rope with this B0, matching the observed field magnitude at both spacecraft shows that a plausible flux was adopted, not that the model predictively determines the field strength. The smooth rotation itself is largely a geometric consequence of a uniformly twisted flux rope with the observed GCS orientation and the assumed helicity sign (via the hemispheric rule, which holds only 60–75% of the time). Thus the load-bearing inference that the model 'can correctly capture the radial evolution of the magnetic field' is not independently established: the radial-evolution test is entangled with the choice of the flux normalization. The paper would need to demonstrate that the same assumed flux is required by other evidence, or quantify how the Earth-spacecraft comparison degrades when the flux is varied over its plausible range.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper simulates the 2020 April 14 stealth CME using the MS-FLUKSS inner-heliosphere MHD model with an inserted constant-turn flux rope (CTFR), and compares the results with in situ magnetic field and plasma data at Solar Orbiter and Earth, which were in a near-radial alignment. The authors also compare synthetic and observed STEREO HI J-maps to assess the CME's heliospheric kinematics. The central claim is that the CTFR model reproduces the smoothly rotating magnetic field signature, especially BT and BN, at both spacecraft with very good accuracy, and that it can correctly capture the radial evolution of the magnetic field of this ICME. The paper explicitly acknowledges that the poloidal flux and helicity sign of the stealth CME could not be observed and were assumed, and that the simulated ambient solar wind is too fast, leading to compensating arrival-time errors.","tokens_in":15223,"tokens_out":3628,"duration_ms":36281,"significance":"If the central claim were established, the study would strengthen the case for using flux-rope-based MHD models rather than hydrodynamic cone models for space-weather forecasting of ICME magnetic fields. The use of two radially aligned spacecraft is a valuable and relatively rare test bed, and the synthetic J-map comparison is a useful diagnostic that goes beyond single-point in situ comparisons. The authors are also commendably explicit about the limitations of their setup, including the assumed poloidal flux, the ad-hoc total energy definition, the high simulated temperatures, and the cancellation of kinematic errors. However, because the magnetic field amplitude is set by an assumed flux that is not independently constrained for this event, the abstract's 'very good accuracy' claim and the Section 5 conclusion about correctly capturing radial magnetic-field evolution are currently stronger than the evidence supports.","major_comments":[{"comment":"The assumed poloidal magnetic flux of 10 x 10^21 Mx, explicitly stated in Section 3.2 to be unconstrained by observations, directly determines B0 = 0.0031 G in the Vandas & Romashets (2017) field solution and therefore the amplitude of the simulated BT and BN components. Matching the observed field magnitude at SolO and Earth is thus a consistency check on a chosen normalization, not an independent validation of the model's predictive capability. To support the Section 5 claim that the CTFR model 'can correctly capture the radial evolution of the magnetic field associated with this ICME,' the authors should vary the assumed poloidal flux (and hence B0) over a plausible range and show how the Earth comparison degrades, or provide an independent observational constraint on the flux with quantified uncertainty.","section":"Section 3.2, Section 4.2, Section 5"},{"comment":"The two-spacecraft radial-evolution test is weakened by the 5-hour time shift applied to the SolO simulation data and by the compensating kinematic errors documented in Section 4.3: the CME is inserted 12 hours late, the simulated ambient solar wind is too fast, and the synthetic J-map slope is steeper than observed. The magnetic field comparison at SolO is therefore only performed after shifting the time axis, and the agreement at Earth benefits from an accidental cancellation of insertion-time and ambient-wind errors. The paper should show the unshifted SolO comparison, and should explicitly separate the magnetic-field-shape agreement from the arrival-time/radial-localization agreement when claiming that radial evolution is captured.","section":"Section 4.2, Figure 5, Section 5"},{"comment":"The helicity sign is assumed from the hemispheric rule, which the paper itself notes holds for only 60-75% of flux ropes, and the toroidal flux is derived from an empirical relation to the assumed poloidal flux. Since the sense of the magnetic field rotation in the simulation depends directly on this assumed sign, the successful reproduction of the rotation direction is not an independent confirmation of the model's helicity treatment. A sensitivity test in which the opposite helicity sign is used would help quantify how much of the claimed agreement depends on this assumption.","section":"Section 4.2, Section 4.3"}],"minor_comments":[{"comment":"The text contains a typo: 'The the conclusions are presented in Section 5' should read 'The conclusions are presented in Section 5.'","section":"Section 1"},{"comment":"The sentence 'we assumed the total mass in the flux rope to be 10 10 kg' has a missing superscript; it should read 10^10 kg.","section":"Section 4.2"},{"comment":"The text refers twice to the 'top-left panel' of Figure 6 when describing the synthetic J-map; the second occurrence should refer to the top-right panel.","section":"Section 4.3, Figure 6"},{"comment":"The ad-hoc total energy density formula includes the solar wind pressure and magnetic energy but no explicit thermal pressure term for the flux rope plasma itself; if this is intentional, a brief explanation of how the flux rope temperature is effectively set would help the reader.","section":"Section 4.2, equation for etotal"},{"comment":"The phrase 'the simulated ICME arrival was underestimated by 5 hours at SolO and overestimated by 5 hours at Earth' is confusing because the simulation arrives late at SolO and early at Earth; consider rephrasing to 'the simulated arrival was 5 hours late at SolO and 5 hours early at Earth.'","section":"Section 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and addresses a timely topic in heliospheric modeling and space weather. The main weakness is not a computational error but an overstatement of the validation strength given the assumed magnetic flux, which is acknowledged in the text. I would encourage the editor to request the sensitivity analysis described in Major Comment 1 before reconsidering the paper; without it, the central claim is not independently supported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a careful single-event application of the CTFR model to a stealth ICME seen at SolO and Earth in a rare near-radial alignment. The rotating field signature at both spacecraft is reproduced, and the synthetic J-map error-cancellation analysis is genuinely new. The soft spot is the one flagged in the stress test: the poloidal flux is assumed, not measured, and that sets the field amplitude. So the good magnitude agreement is not an independent test of the model. The authors say this plainly in Section 3.2, which I credit. The helicity sign is also assumed from the hemispheric rule (60-75% reliability). Still, given those assumptions, the match on sign and rotation of BT and BN is significant.\n\nOther limitations: the ambient SW is too fast, density and speed jumps are not reproduced, temperature is too high, and the 5-hour shifts at SolO and the error cancellation at Earth are acknowledged. None of this is hidden. The paper is honest about what the model can and cannot do.\n\nMy judgment: as a validation of CTFR for forecasting, it is incremental, not decisive, because the flux assumption weakens the amplitude claim. But as a modeling study of a stealth CME with a well-chosen multi-spacecraft event, it is solid and worth publishing after minor revision. I would send it to review, with a request that the authors either test sensitivity to the assumed flux over a plausible range or soften the 'very good accuracy' wording. I would also bring it to my group's reading session as an example of how model-data comparison can be done honestly with acknowledged assumptions.","headline":"A careful, honest single-event application of the CTFR model to a stealth ICME in a rare SolO–Earth radial alignment; the field rotation is reproduced well, but the assumed poloidal flux means the amplitude match is not an independent test.","tokens_in":15746,"tokens_out":2341,"would_cite":true,"duration_ms":23703,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A flux-rope MHD model reproduces the rotating magnetic field of a stealth CME at Solar Orbiter and Earth.","keywords":["coronal mass ejection","ICME","flux rope","constant-turn flux rope","magnetohydrodynamics","Solar Orbiter","space weather forecasting","stealth CME"],"falsifier":"Compute the axial magnetic flux through the simulated ejecta at 1 au and compare it with the flux obtained by integrating the observed $B_T$ profile across the in-situ ICME at Earth; if the simulated flux disagrees with the observed flux by more than the model's claimed accuracy, the field magnitude is inherited from the assumed $10\\times10^{21}$ Mx poloidal flux rather than being a dynamical prediction. A second check is to run the same CTFR setup for a radially aligned ICME pair whose magnetic profiles are not well correlated (as in Regnault et al. 2023); if the model still produces a smooth, matching rotation, the current success is specific to this event's simple structure rather than a general property of the model.","tokens_in":14732,"feed_emoji":"☀️","tokens_out":13479,"duration_ms":113274,"temperature":0.7,"pith_summary":"The paper tests whether the constant-turn flux rope (CTFR) model, embedded in a time-dependent inner-heliosphere magnetohydrodynamic simulation, can reproduce the magnetic field of an interplanetary coronal mass ejection (ICME) at two spacecraft lying nearly along the same radial line from the Sun. The event is a stealth CME, invisible in EUV, so the flux rope's magnetic flux cannot be measured: the authors insert a fully formed flux rope with an assumed poloidal flux of $10\\times10^{21}$ Mx and compare the result with Solar Orbiter and Earth in situ data plus STEREO-A heliospheric images. They find that the simulated magnetic field reproduces the sign, magnitude, and smooth rotation of the transverse and normal components, with the smaller radial component also matching, at both spacecraft. If correct, this supports using flux-rope MHD models to forecast the southward magnetic field component that controls whether an ICME drives a geomagnetic storm; the arrival-time errors of +5 hours at Solar Orbiter and -5 hours at Earth illustrate that kinetic agreement can be partly accidental, because a late insertion and an over-fast simulated solar wind cancel.","feed_headline":"MHD model reproduces stealth CME's magnetic field at SolO and Earth","feed_subtitle":"It reproduces the rotating magnetic field at Solar Orbiter and Earth, supporting flux-rope MHD forecasts of storm-driving fields.","key_machinery":"The load-bearing machinery is the constant-turn flux rope (CTFR) model: a croissant-shaped flux rope with circular cross-section and twin legs (the FRiED geometry) whose internal magnetic field is the uniform-turn analytic solution of Vandas & Romashets, configured with poloidal and toroidal fluxes, helicity sign, tilt, half-angle, and aspect ratio. The flux rope supplies the initial magnetic field and low-plasma-$\\beta$ energy density of a fully formed ejecta inserted at 0.1 au into the ambient solar wind; the paper's test is whether that initial field, advected and expanded by the MHD evolution, reproduces the observed rotating field at Solar Orbiter and Earth. The assumed poloidal flux of $10\\times10^{21}$ Mx, converted to toroidal flux via the empirical relation of Qiu et al. (2007), sets the field amplitude $B_0$; the helicity sign follows from the hemispheric helicity rule.","core_discovery":"On its own terms, the paper's discovery is that the CTFR model can correctly capture the radial evolution of the magnetic field associated with this ICME: when simulation output is extracted along the trajectories of Solar Orbiter and Earth, the signs and magnitudes of $B_T$ and $B_N$ within the ejecta match in situ observations, and the $B_R$ component stays small as observed, so the smooth rotation of the magnetic field across the ICME is preserved at both heliocentric distances. The paper presents this as evidence that the CTFR model is a good candidate for forecasting the magnetic structure, and therefore the geo-effectiveness, of ICMEs with simple magnetic structure. It also reports that the simulated ICME arrives 5 hours late at Solar Orbiter and 5 hours ahead at Earth, and that the J-map comparison shows the CME entered the model 12 hours late but then moved too fast in an over-speedy simulated solar wind, with the two errors canceling near Earth.","pith_inferences":["Beyond the paper, the close match of $B_T$ and $B_N$ at two radial distances suggests the interior magnetic profile of an ICME is largely set by the initial flux-rope parameters and only weakly modified by interaction with the solar wind; this could be tested by repeating the simulation with the same flux rope in different ambient-wind realizations and checking whether the in-ejecta field profile ","A testable extension is to compare the simulated and observed total axial flux at 1 au for this event; because the poloidal flux was assumed, such a comparison would separate the model's structural skill from the amplitude set by the assumption.","The authors leave implicit that their success is for an ICME with a simple, well-correlated magnetic structure; the harder forecasting case is events where radial alignment still shows poor correlation, and their own discussion points to applying CTFR to those events as the next test.","Another implication, not pursued in the paper, is that the cancellation of errors in arrival time is itself a warning for operational forecasting: matching the observed shock time at Earth does not constrain the simulated CME's speed history, so ensemble forecasts should sample both the ambient solar wind and the insertion time."],"forward_implications":["If the central claim is right, a data-constrained flux-rope MHD model can predict the sign and magnitude of the magnetic field components that determine an ICME's geo-effectiveness, rather than only its arrival time.","The same model should be applied to other radially aligned ICME observations, especially events where the field profiles at the two spacecraft are poorly correlated, to see whether the present match reflects the model or the event's unusually simple flux-rope structure.","Forecasters should treat an agreement in arrival time as a compound result: here a 12-hour late insertion and a faster-than-real simulated solar wind canceled, so a single arrival-time score does not validate the ambient-wind or the CME initialization.","For stealth CMEs, where magnetic flux cannot be measured, the assumed poloidal flux and the hemispheric helicity rule are enough to reproduce the observed rotation pattern, which is the structure needed for a storm forecast; the amplitude remains conditional on the assumed flux."],"supporting_citations":[{"why":"Introduces the CTFR model and its self-similar expansion velocity profile, the method this paper applies.","marker":"Singh et al. (2022)"},{"why":"Supplies the uniform-turn analytic solution used for the flux rope's internal magnetic field.","marker":"Vandas & Romashets (2017)"},{"why":"Provides the FRiED croissant-like global flux rope geometry.","marker":"Isavnin (2016)"},{"why":"Documented the closely correlated magnetic field profiles of this ICME at Solar Orbiter and Earth, defining the target the simulation must match.","marker":"Davies et al. (2021)"},{"why":"Gives the empirical relation used to set toroidal flux from the assumed poloidal flux.","marker":"Qiu et al. (2007)"},{"why":"Supplies the hemispheric helicity rule used to choose the flux rope's twist sign.","marker":"Pevtsov & Balasubramaniam (2003)"},{"why":"The GCS model used to fit coronagraph observations and set direction, tilt, half-angle, and aspect ratio.","marker":"Thernisien et al. (2009)"},{"why":"The drag-based model used to estimate the CME insertion time and speed at $70\\,R_\\odot$.","marker":"Vršnak & Žic (2007)"}],"fun_headline_variants":["CTFR model reproduces stealth CME field rotation at SolO and Earth","MHD simulation matches CME field rotation at SolO and Earth","Stealth CME field rotation reproduced by MHD model at two points","Data-constrained model reproduces CME magnetic field at two spacecraft","Simulation captures rotating magnetic field of stealth CME at SolO and Earth"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the flux rope's magnetic flux is $10\\times10^{21}$ Mx, an assumed value chosen because the stealth CME offers no EUV or magnetogram signature with which to measure the flux; if the true flux differs, the simulated field magnitude, and hence the claimed accuracy, changes even if the field geometry is right.","fun_headline_variants_meta":{"raw":{"variants":["CTFR model reproduces stealth CME field rotation at SolO and Earth","MHD simulation matches CME field rotation at SolO and Earth","Stealth CME field rotation reproduced by MHD model at two points","Data-constrained model reproduces CME magnetic field at two spacecraft","Simulation captures rotating magnetic field of stealth CME at SolO and Earth"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000474,"raw_usage":{"total_tokens":2413,"prompt_tokens":1064,"completion_tokens":1349,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":680,"completion_tokens_details":{"reasoning_tokens":1253}},"tokens_in":680,"tokens_out":1349,"duration_ms":11106,"temperature":1.0,"reasoning_tokens":1253,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T19:55:41.679494+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the axial magnetic flux through the simulated ejecta at 1 au and compare it with the flux obtained by integrating the observed $B_T$ profile across the in-situ ICME at Earth; if the simulated flux disagrees with the observed flux by more than the model's claimed accuracy, the field magnitude is inherited from the assumed $10\\times10^{21}$ Mx poloidal flux rather than being a dynamical prediction. A second check is to run the same CTFR setup for a radially aligned ICME pair whose magnetic profiles are not well correlated (as in Regnault et al. 2023); if the model still produces a smooth, matching rotation, the current success is specific to this event's simple structure rather than a general property of the model.","supporting_citations":[{"cited_title":"A., & Balasubramaniam, K","cited_arxiv_id":null,"evidence_quote":"Supplies the hemispheric helicity rule used to choose the flux rope's twist sign."}],"review_version":1}