{"id":"365b5090-e7e1-4468-b1c6-00dcc10b3d7e","arxiv_id":"1908.04450","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"All seven fast CME-driven shocks analyzed here completed their major deceleration inside 1 au and then propagated with gradual deceleration beyond Earth's orbit.","lead":"This paper tracks seven fast coronal mass ejections from the Sun to beyond Earth's orbit using coronagraph images, radio bursts, and multi-spacecraft measurements. It concludes that all seven shocks finished their strong deceleration before 1 au and then slowed gradually on the way to Ulysses.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Beyond-1 au gradual-deceleration claim rests on unproven Earth-Ulysses shock identity; the paper itself admits a possible alternative CME for 1997 Nov 4 and a substituted ICME edge for 2005 May 13.","rationale":"The paper's Sun-to-Earth conclusion rests on measured endpoint speeds plus type II drift and is reasonably secure. The new, headline contribution is the beyond-1 au behavior, and that part depends on treating the Ulysses-observed shock as the same physical front tracked from Earth. The manuscript itself concedes the two weakest identifications: for 1997 November 4 a later CME from the same AR could explain the Ulysses shock, and for 2005 May 13 the Ulysses signature is not a shock at all. Because Table 1 lists separations up to 138 degrees in longitude, the Earth-versus-Ulysses speed difference could also be a front-geometry effect rather than a temporal deceleration of the same shock element. The reader's weakest-assumption statement captures exactly this concern, and I agree with it. I would not change the CONDITIONAL verdict: the central claim is plausible and the MHD model gives rough consistency, but the beyond-1 au gradual-deceleration statement should be explicitly conditioned on shock identification and, for at least one event, on excluding the alternative CME. A directed re-analysis of the 1997 November 6 CME propagation would settle the highest-risk case.","tokens_in":19864,"tokens_out":13643,"duration_ms":141581,"concrete_test":"Run a 3D heliospheric MHD model (e.g., ENLIL or EUHFORIA) initialized with both the 1997 November 4 and 1997 November 6 CMEs from AR 08100, using their LASCO-derived initial speeds and source positions, and predict the shock arrival time and speed at Ulysses. If the November 6 CME can produce a shock at Ulysses near 1997 November 24 15:07 UT with speed about 410 km/s, then event 1's Earth-Ulysses shock association is ambiguous and the beyond-1 au conclusion for that event is unsupported. If the November 6 contribution arrives at a clearly different time or speed, the paper's identification is upheld for the one event where the text explicitly concedes an alternative.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step for the novel beyond-1 au part of the claim is the identification of each Ulysses shock as the same CME-driven shock observed at Earth, and the further assumption that the Earth-to-Ulysses speed drop is temporal deceleration rather than a spatial (nose-to-flank) property of the front. This is not independently secured. Table 1 shows Ulysses separated by up to 138 degrees in longitude and as much as 73 degrees in latitude from Earth; with such separations, the same expanding shock can arrive at different local speeds for geometrical reasons. The paper's own text weakens the association in two places: Section 3.1 states that the 1997 November 6 CME from the same active region could have reached Ulysses and 'cannot be completely excluded', and Section 3.6 states for the 2005 May 13 event that 'the CME-driven shock missed Ulysses, and the ICME leading edge has been used as a substitute'. If the 1997 November 24 Ulysses shock is actually the November 6 CME shock, or if the 2005 Ulysses signature is an ICME leading edge rather than the shock, then the speeds listed in Table 2 for those events are not shock speeds of the same CME, and the statement that all 7 shocks 'thereafter propagated with a gradual deceleration' loses its support. The Sun-to-Earth part of the claim would survive such an error, but the beyond-1 au gradual-deceleration conclusion would not be established for the affected events.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a kinematic analysis of 7 fast CME-driven shocks observed in solar cycle 23, combining SOHO/LASCO coronagraph data via GCS forward modeling, Wind/WAVES type II radio bursts, in-situ solar wind measurements at L1 (Wind/ACE/Genesis) and Ulysses, a simple analytical propagation model with three input parameters, and a 1D MHD model. For each event the authors derive the near-Sun speed v0, the shock speed and transit time at Earth (vs, ts), and the shock speed at Ulysses, then use the analytical model to obtain the deceleration a, deceleration duration ta, and deceleration cessation distance ra within 1 au. They report that all 7 shocks decelerated strongly before reaching 1 au and then continued with gradual deceleration beyond 1 au, with cessation distances ranging from 0.28 to 0.60 au. They also report correlations among the kinematic parameters and discuss the geomagnetic consequences of several events, including a shock–ejecta interaction that enhanced southward fields.","tokens_in":20236,"tokens_out":7505,"duration_ms":75411,"significance":"If the claims are accepted, the paper adds a small multi-event sample supporting the picture that fast CME-driven shocks complete most of their Sun–Earth deceleration within roughly 0.3–0.6 au and continue to slow gradually beyond 1 au, and it demonstrates a simple analytical model as a useful forecasting tool. The Sun-to-Earth part of the evidence is strong: for every event, the near-Sun GCS speed (1200–2600 km/s) exceeds the in-situ speed at L1 (560–1040 km/s), and the type II drift curves broadly match the model after per-event adjustment of the 1-au density scale n0. The paper also makes a transparent comparison with a 1D MHD model and is honest about several limitations, notably the possible alternative CME for the 1997 November 4 event and the substituted ICME leading edge for 2005 May 13. The beyond-1 au conclusion and the statistical correlations, however, are weakened by the issues detailed below.","major_comments":[{"comment":"The correlations in Figure 22 are presented as empirical findings, but several are mathematical consequences of the analytical model's definitions. In Section 2.3, ta = 2(ds - vs*ts)/(v0 - vs), a = (vs - v0)/ta, and ra = (v0 + vs)*ta/2 are deterministic functions of the three inputs v0, vs, and ts. Consequently, the near-perfect anti-correlation between a and ta (r = -0.91) is essentially a restatement of ta = (vs - v0)/a, not an independent physical relationship. The abstract's claim that \"the faster the CME, the larger the deceleration, and the shorter the deceleration time period\" is therefore partly a property of the model, not an independent empirical trend. The authors should either compute the correlations using independently measured quantities (e.g., deceleration fits to type II drift without n0 tuning) or explicitly state that the correlations are algebraic consequences of the model. This issue is load-bearing for the statistical conclusions in Section 4 and the abstract.","section":"Section 4 and Section 2.3"},{"comment":"The identification of the Ulysses shock as the same CME-driven shock observed at Earth is not securely established for at least two of the seven events, which undermines the \"all 7\" claim of gradual deceleration beyond 1 au. In Section 3.1 the authors note for 1997 November 4 that a later CME from the same active region could have reached Ulysses and \"cannot be completely excluded\"; in Section 3.6 they state for 2005 May 13 that \"the CME-driven shock missed Ulysses, and the ICME leading edge has been used as a substitute.\" Including these events in the statement that each shock \"thereafter propagated with a gradual deceleration\" is not supported by the presented evidence. Moreover, Ulysses was separated from Earth by up to 138 degrees in longitude (Table 1), so a lower speed at Ulysses could reflect the flank of the shock rather than temporal deceleration. The authors should either exclude these two events from the beyond-1 au conclusion or provide stronger evidence for the association (e.g., compositional signatures, multi-spacecraft timing consistency, or modeling that accounts for the shock's non-radial structure).","section":"Section 3.1, Section 3.6, Table 1"},{"comment":"The analytical model is repeatedly described as having \"no free parameters,\" but the validation against the type II radio bursts uses a per-event adjusted value of n0, the nominal 1-au density in the Leblanc model. Because n0 is tuned to make the model's frequency-time profile match the observations, the type II drift does not provide an independent verification of the model; it merely shows that a density scaling can reconcile the two. The statement in Section 5 that the model is \"verified\" by agreement with the type II bursts is therefore overstated. The authors should report the n0 values used for each event and discuss how sensitive the conclusions are to the density-model scaling, or else clearly separate the model's three input parameters from the validation step's tuning parameter.","section":"Sections 2.2, 2.3, and 5"},{"comment":"The MHD model is presented as an independent check of the Earth–Ulysses connection, but its predicted shock arrival times at Ulysses differ from the observed times by 21–39 hours (e.g., 39 hr early for 1997 November 4, 37 hr early for 2000 June 6, and 26 hr early for 2001 April 2). For the 2001 April 2 event, 26 hours is more than half of the roughly 46-hour Earth-to-Ulysses transit time, which the paper itself classifies as \"a larger uncertainty compared with other cases.\" These discrepancies, attributed to large longitudinal/latitudinal separations, are large enough that the MHD comparison does not convincingly confirm that the same shock was observed at both locations. The claim in Section 2.4 that the MHD model \"can also be used to evaluate the analytical model beyond 1 au\" is thus not strongly supported by the presented results.","section":"Sections 3.1–3.4 and 2.4"}],"minor_comments":[{"comment":"There is a typographical error: \"analyitcal model\" should be \"analytical model.\"","section":"Section 3.7"},{"comment":"The note to Table 1 states that the Ulysses latitude and longitude are given relative to Earth when the CME erupted, but Section 3.2 gives the Ulysses position \"when the shock arrived there\" (e.g., S58.4). The authors should specify which epoch is used consistently throughout the paper and in the table.","section":"Table 1 and Section 3.2"},{"comment":"The correlation coefficients are reported without p-values or confidence intervals. With only 7 events, the statistical significance is low, and the reader cannot assess the strength of the relations shown in panels (c), (d), (e), and (f); the paper should either add significance measures or temper the language about these correlations.","section":"Figure 22"},{"comment":"The statement that Figure 22(f) shows \"no obvious correlation\" would be better phrased as \"not statistically significant at this sample size,\" since with N=7 the absence of a visible trend is expected even if a weak relationship exists.","section":"Section 4"},{"comment":"The type II frequency-to-distance conversion depends on the Leblanc density model and the chosen n0 values. For reproducibility, the authors should list the n0 values used for each event in a table or in the text, rather than only in the figures.","section":"Section 2.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is built on a relatively small sample (7 events) and relies heavily on the analytical model proposed by the same group (Liu et al. 2017), with three of the seven events taken from previous papers by the same authors. The novelty therefore lies mainly in the systematic comparison and the statistical correlations, but the correlations are partly definitional, and the beyond-1 au conclusion depends on shock associations that are acknowledged to be questionable for two events. The Sun-to-Earth deceleration claim itself is well supported and likely publishable, but the current manuscript overstates the strength of the beyond-1 au and statistical conclusions. A major revision that re-frames the claims around the robust Sun-to-Earth result and explicitly separates model-built-in relationships from independent measurements would make this a sound contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper gives you seven fast CMEs tracked from the Sun to Ulysses using coronagraph/GCS, type II radio drift, and in-situ data, with 1D MHD as a bridge. The four new event studies (1997 Nov 4, 2000 Jun 6, 2001 Apr 2, 2001 Nov 4) are careful, and the height-time profiles are genuinely useful as a reference set. The finding that each shock decelerates hard inside 1 au is well supported: near-Sun speeds of 1200–2600 km/s against 1 au speeds of 560–1040 km/s, and the type II drift is broadly consistent after a per-event n0 adjustment. I believe the Sun-to-1 au part of the claim.\n\nThe soft spots match where the stress-test landed. The beyond-1 au claim—“thereafter propagated with a gradual deceleration”—rests on identifying the Ulysses shock as the same structure observed at Earth, and the paper itself concedes two cracks: for 1997 Nov 4, a later CME from the same active region “cannot be completely excluded” as the Ulysses shock; for 2005 May 13, the CME-driven shock missed Ulysses and the ICME leading edge was substituted. With Ulysses up to 138 degrees in longitude away, a shock's local speed at Ulysses can differ from the nose speed for geometrical reasons, so speed drops beyond 1 au are not unambiguous evidence of temporal deceleration. The 1D MHD model also predicts the Ulysses arrival 20–40 hours early in several cases, which the authors attribute to the Earth-Ulysses separation; that's plausible, but it means the beyond-1 au agreement is only qualitative. That doesn't kill the Sun-to-Earth result, but it does mean the beyond-1 au part is not established for those events. Also note the “no free parameters” slogan is only half true: the analytical model is parameter-free, but the radio verification tunes n0 per event.\n\nThe correlation analysis in Figure 22 should be read with caution. Since a, ta, and ra are deterministic functions of v0, vs, and ts, some of the correlations—especially the -0.91 between a and ta—are near definitional. The authors do warn that the sample is small, but the framing “verifies the analytical model” overstates what the correlations add. And there are no error bars anywhere; GCS speeds and type II drift widths give some sense of uncertainty, but the quoted model parameters come without it.\n\nWho is this for? Someone assembling event catalogs of CME propagation or planning arrival-time forecasts for fast shocks will get real value from the event-by-event profiles and the tabulated parameters. The paper deserves a serious referee—the workload is justified by the new event analyses, even though the statistical claims need heavy revision. I'd recommend engaging with it, but expect that the referee will require qualification of the beyond-1 au and correlation statements.","headline":"Useful multi-event compilation of fast CME shock propagation Sun-to-Ulysses, but the beyond-1 au gradual-deceleration conclusion leans on unproven shock identities and correlations that are partly baked into the model.","tokens_in":20747,"tokens_out":3116,"would_cite":true,"duration_ms":30532,"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":"Seven fast CME-driven shocks complete their major deceleration before reaching Earth's orbit.","keywords":["coronal mass ejections","interplanetary shocks","CME deceleration","Ulysses","type II radio bursts","solar wind","space weather","MHD modeling"],"falsifier":"Take the next fast CME observed at both L1 and a well-separated outer-heliosphere spacecraft, and check the predicted arrival time from the analytical model against the observed arrival while confirming the two shocks are the same structure using energetic particle signatures or magnetic field orientation. A systematic discrepancy beyond the model's stated uncertainties, or identification of a different event at Ulysses in any of the seven cases, would overturn the claim.","tokens_in":19676,"feed_emoji":"☀️","tokens_out":5809,"duration_ms":52189,"temperature":0.7,"pith_summary":"Coronal mass ejections (CMEs) are large expulsions of magnetized plasma from the Sun, and the shocks they drive can disturb space weather all the way past Earth. This paper examines seven fast CMEs from solar cycle 23 that were observed by coronagraphs, tracked by interplanetary type II radio bursts, and measured in situ at both the Earth and Ulysses, and it argues that in every case the shock's major deceleration was completed before 1 au, followed by gradual deceleration beyond Earth. The evidence is assembled into a simple analytical model with no free parameters, and the resulting Sun-to-Ulysses height-time profile is roughly consistent with all available data. If the result holds, a fast shock's long-range propagation can be captured by three measured quantities: its speed near the Sun, its speed at 1 au, and its Sun-Earth transit time.","feed_headline":"Fast CME shocks finish major braking inside 1 au","feed_subtitle":"Seven Earth-to-Ulysses events show a simple three-input model tracks the slowdown all the way out.","key_machinery":"The load-bearing tool is an analytical, no-free-parameter model of shock propagation. It assumes the shock starts with initial speed $v_0$, decelerates at constant rate $a$ for time $t_a$, and then moves at constant speed $v_s$; combining $v_0$, $v_s$, and the transit time fixes $a$, $t_a$, and the cessation distance $r_a=(v_0+v_s)t_a/2$. The Leblanc density model converts the model's distances into radio frequencies to check against type II burst drift, a 1D MHD model propagates near-Earth solar wind data outward to test the same shock at Ulysses, and the graduated cylindrical shell model supplies the near-Sun speed and direction.","core_discovery":"The central discovery is a propagation pattern for fast CME-driven shocks. Each event is reduced to three measured quantities: the shock speed near the Sun from a graduated cylindrical shell fit, the shock speed at 1 au from in-situ data, and the Sun-Earth transit time. From these inputs the model derives a constant deceleration lasting a time $t_a$ and ending at a distance $r_a$, followed by motion at roughly constant speed or with gradual deceleration. Across the seven events the deceleration cessation distance lies between 0.28 and 0.60 au, so the rapid braking is over before the shock reaches Earth. The paper reads the agreement of the same profile with type II radio burst drift inside 1 au, with MHD model output, and with Ulysses arrival times beyond 1 au as verification that the simple model captures real behavior.","pith_inferences":["If the same three-input profile scales to the outer heliosphere, shocks observed near Voyager distances should arrive slightly later than a constant-speed extrapolation from 1 au would predict; checking this is a direct test the paper does not carry out.","The Earth-Ulysses pairing is the load-bearing link beyond 1 au; independent confirmation through energetic particle signatures or magnetic field orientation at both spacecraft would make the gradual-deceleration conclusion much harder to doubt.","Because the sample contains only fast CMEs, the claimed correlations between initial speed, deceleration, and deceleration time should be tested on slower events and on events observed by more than one outer-heliosphere spacecraft."],"forward_implications":["A fast shock's arrival farther out can be predicted from just its near-Sun speed, its 1 au speed, and its Sun-Earth transit time.","The rapid-deceleration phase ends between 0.28 and 0.60 au for the seven events studied, so major braking is a Sun-to-sub-1-au phenomenon.","Faster CMEs decelerate more strongly and for a shorter time inside 1 au, with linear correlation coefficients near 0.79 and -0.91 respectively.","Southward sheath fields, rather than the ejecta itself, drove several of the associated geomagnetic storms, and shock-ejecta interaction can intensify them, as in the severe 2001 November 4 storm.","A three-input kinematic profile can describe the Sun-to-Ulysses transit of fast shocks using data that are already available in real time at L1."],"supporting_citations":[{"why":"Supplies the no-free-parameter analytical model and the 2001 November 22 event parameters used in the statistical sample.","marker":"Liu et al. (2017)"},{"why":"Provides the solar wind density model that converts type II radio frequencies into shock distances.","marker":"Leblanc et al. (1998)"},{"why":"Introduces the graduated cylindrical shell model used to estimate near-Sun CME speed and direction.","marker":"Thernisien et al. (2006, 2009)"},{"why":"Provides the 1D MHD model that propagates near-Earth solar wind data outward to test shock arrival at Ulysses.","marker":"Wang et al. (2000)"},{"why":"Gives the earlier Sun-to-1 au kinematic analysis and deceleration cessation distances that the results are compared with.","marker":"Reiner et al. (2007)"},{"why":"Establishes the typical fast-CME profile of impulsive acceleration, rapid deceleration, then constant or gradual speed that the seven events are interpreted to follow.","marker":"Liu et al. (2013)"},{"why":"Provides the previous analysis and parameters for the 2005 May 13 CME included in the statistics.","marker":"Zhao et al. (2017)"},{"why":"Provides the previous analysis and shock parameters for the 2006 December 13 CME.","marker":"Liu et al. (2008)"}],"fun_headline_variants":["CME shocks slam the brakes before reaching Earth","Fast CMEs decelerate hard inside 1 AU, study finds","Three numbers predict CME shock slowdown to Ulysses","Sun-to-Ulysses CME tracking: braking zone mapped","Major CME deceleration ends by 0.6 AU, model says"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument rests on identifying the shock seen at Ulysses as the same structure as the shock seen at Earth, based mainly on similar transit speeds even though Ulysses was up to roughly 138 degrees of longitude away; if that pairing is wrong for any event, the inferred gradual deceleration beyond 1 au loses its observational support.","fun_headline_variants_meta":{"raw":{"variants":["CME shocks slam the brakes before reaching Earth","Fast CMEs decelerate hard inside 1 AU, study finds","Three numbers predict CME shock slowdown to Ulysses","Sun-to-Ulysses CME tracking: braking zone mapped","Major CME deceleration ends by 0.6 AU, model says"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000617,"raw_usage":{"total_tokens":2901,"prompt_tokens":1021,"completion_tokens":1880,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":637,"completion_tokens_details":{"reasoning_tokens":1794}},"tokens_in":637,"tokens_out":1880,"duration_ms":14440,"temperature":1.0,"reasoning_tokens":1794,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:42:23.001769+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the next fast CME observed at both L1 and a well-separated outer-heliosphere spacecraft, and check the predicted arrival time from the analytical model against the observed arrival while confirming the two shocks are the same structure using energetic particle signatures or magnetic field orientation. A systematic discrepancy beyond the model's stated uncertainties, or identification of a different event at Ulysses in any of the seven cases, would overturn the claim.","supporting_citations":[{"cited_title":"D., Zhao, X., & Zhu, B","cited_arxiv_id":null,"evidence_quote":"Supplies the no-free-parameter analytical model and the 2001 November 22 event parameters used in the statistical sample."},{"cited_title":"A., & Bougeret, J.-L","cited_arxiv_id":null,"evidence_quote":"Provides the solar wind density model that converts type II radio frequencies into shock distances."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the graduated cylindrical shell model used to estimate near-Sun CME speed and direction."},{"cited_title":"J., Kaiser, M","cited_arxiv_id":null,"evidence_quote":"Gives the earlier Sun-to-1 au kinematic analysis and deceleration cessation distances that the results are compared with."},{"cited_title":"D., Hu, H., & W ang, R","cited_arxiv_id":null,"evidence_quote":"Provides the previous analysis and parameters for the 2005 May 13 CME included in the statistics."},{"cited_title":"G., M¨ uller-Mellin, R., et al","cited_arxiv_id":null,"evidence_quote":"Provides the previous analysis and shock parameters for the 2006 December 13 CME."}],"review_version":1}