{"id":"16fbc023-cc91-4e89-a42f-f5fb064b8ace","arxiv_id":"2505.22283","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"ICMEs show no radial trend in intermittency, while their sheaths become more intermittent with distance from the Sun.","lead":"Researchers analyzed magnetic field fluctuations in 49 coronal mass ejections (ICMEs) seen by Parker Solar Probe and Solar Orbiter. They found the burstiness of these fluctuations stays constant as ICMEs travel outward, while the shocked plasma ahead of them becomes more bursty with distance.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Radial invariance is an equivalence claim, but the paper only shows a non-significant Kendall correlation for ICMEs (τ=0.11) and provides no confidence intervals or equivalence test.","rationale":"The reader's stated weakest assumption concerned contamination by flux-rope rotation and boundary structures; that is a valid threat but relies on assumptions about what the residual contains. I find a more immediate, testable gap: the key empirical support for radial invariance is a null correlation reported without any uncertainty. The paper gives no p-values or confidence intervals for any entry in Table 1; for n=49, Kendall τ=0.11 has a 95% confidence interval that comfortably includes the sheath value of 0.30. Therefore the contrast that drives the headline is not established. This is not a claim of misconduct; it is a missing inferential step. The appropriate fix is straightforward and would either support or weaken the conclusion. I therefore keep the reader's CONDITIONAL verdict: the analysis is plausible but the central interpretation is currently under-constrained. If the equivalence test fails, the verdict should move toward REJECT or UNVERDICTED for the radial-invariance claim, even though the descriptive statistics remain useful.","tokens_in":15751,"tokens_out":7960,"duration_ms":92621,"concrete_test":"Reanalyze Table 1 with uncertainty and equivalence testing. For each interval type, compute bootstrap 95% confidence intervals for Kendall's τ (or the least-squares slope of ακ versus R) using the 49 events, stratified by spacecraft; then run a two one-sided test (TOST) for the ICME slope with an equivalence bound of ±0.05 in ακ per au (or using the sheath/upstream slope as a benchmark), and a permutation test comparing the ICME slope with the sheath slope. If the ICME confidence interval includes the sheath slope, or the equivalence test fails, the 'radially invariant' conclusion is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that ICME intermittency is radially invariant and that ICMEs are 'relatively static, well-developed turbulent environments' (Abstract, Sections 4–5) rests on Table 1, where the Kendall correlation between ακ and heliocentric distance for ICMEs is τ=0.11, while the sheath value is τ=0.30. No p-values, confidence intervals, or equivalence tests are reported anywhere. With 49 events, the approximate 95% confidence interval for the ICME τ is about −0.08 to 0.30, so a radial trend as strong as the sheath trend is not excluded. A non-significant correlation cannot establish invariance; the distinction between 'invariant ICME' and 'evolving sheath' is an under-determined null result, not a demonstrated difference. The problem is compounded by the bimodal radial sampling (PSP covers R<0.6 au, SolO covers R>0.6 au): a spacecraft-dependent offset in ακ could mask or create an apparent radial trend in the combined sample. The authors' statement in Section 3.2 that 'the values remain low in other kdi ranges' does not supply the missing uncertainty quantification. This is the load-bearing weakness because the physical interpretation ('fully developed turbulence') is only as strong as the evidence for invariance.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes scale-dependent magnetic-field intermittency in 49 ICMEs observed by Parker Solar Probe and Solar Orbiter between 0.25 and 1 au, together with their sheath, upstream, and downstream solar-wind intervals. It uses structure-function-based kurtosis κ(τ) and its scaling exponent ακ, with a wavelet-based local intermittency measure as a cross-check, and reports Kendall correlations between ακ and heliocentric distance R, speed, spectral index, residual energy, cross helicity, and proton beta. The central claims are that ICMEs show a radially invariant level of intermittency, interpreted as relatively static, fully developed turbulence, while sheaths show increasing intermittency with distance, interpreted as still-evolving turbulence.","tokens_in":16032,"tokens_out":5504,"duration_ms":56509,"significance":"If the central claim is substantiated, this would be a valuable multi-spacecraft statistical result: it would show that large-scale ICME expansion does not, by itself, modify inertial-range intermittency, and it would provide a concrete observational contrast between ICMEs and sheaths. The study has clear strengths: a 49-event sample from two modern spacecraft, a transparent event list in Table 2, use of standard structure-function and wavelet diagnostics, and explicit comparison with earlier case studies. The analysis is reproducible in principle because the events are listed and the data are public. However, the headline distinction between 'radially invariant ICME' and 'evolving sheath' currently rests on low Kendall coefficients without significance tests, confidence intervals, or equivalence tests, and the radial sampling is bimodal between the two spacecraft. These issues are load-bearing for the paper's conclusions.","major_comments":[{"comment":"The conclusion that ICMEs are 'radially invariant' rests on the single Kendall coefficient τ=0.11 between ακ and R for ICMEs, while the sheath value is τ=0.30. No p-values, confidence intervals, or equivalence tests are reported anywhere. With 49 events, the approximate 95% CI for τ=0.11 is about −0.08 to 0.30, so a radial trend as strong as the sheath trend is not excluded by the ICME data, and the claimed ICME-versus-sheath difference is not demonstrated. The statement in §3.2 that 'the values remain low in other kdi ranges' does not supply the missing uncertainty quantification. Please add bootstrap confidence intervals and/or equivalence tests (e.g., two one-sided tests against a pre-specified margin), and report p-values for all entries in Table 1. Without this, 'radial invariance' is an under-determined null result rather than a supported claim.","section":"§3.2, Table 1; §5"},{"comment":"The radial comparison is confounded by the bimodal spacecraft sampling. In Table 2, PSP events are concentrated at R ≤ 0.69 au and SolO events at R ≥ 0.61 au, with few events in the overlap region. If PSP and SolO yield systematically different ακ values due to different noise levels, cadences, or data processing, a spacecraft-dependent offset could either create or mask the apparent radial trend in ακ versus R. Please test for a spacecraft dependence by including spacecraft as a covariate, computing within-spacecraft Kendall coefficients, or overlaying the ακ distributions by spacecraft in Figure 4. This is particularly important for the sheath trend, where τ=0.30 could in principle be driven by a PSP/SolO offset rather than by a true radial evolution.","section":"§3.2, Figure 4; Table 2"},{"comment":"The choice of the fitting range 10^-2.5 < kdi < 10^-1.5 appears post hoc. The text states that power-law fitting was performed in seven one-decade bins, that the two outermost bins were excluded because they include parts of the injection and kinetic scales, and that the analysis then concentrates on one middle range. The paper says that correlation values remain low in other kdi ranges but does not show them. Because the central claim concerns ακ and its correlations, please report the Kendall coefficients for all seven bins (for example, in a supplementary table) and either justify the chosen range a priori or demonstrate that the ICME/sheath distinction is stable across the full inertial range. Small shifts in the fitting range could alter the apparent radial invariance.","section":"§3.2, Figure 3"},{"comment":"The interpretation of ακ as an intermittency measure of homogeneous turbulence assumes that the residual fluctuations, after removing running averages of 1 to 4 hours, are stationary and are not dominated by the large-scale flux-rope rotation or boundary layers. The paper itself notes in §3.1 that 'some systematic inhomogeneities in the distribution of enhanced fluctuation amplitudes across the subintervals may be enhancing the LIM', and it removes only 10% of the ICME interval at each boundary. In addition, the shortest ICMEs have durations near the 3-hour minimum, for which a 4-hour running average is ill-defined, and the stated robustness check is not shown quantitatively. Please provide a stationarity analysis of the residuals (e.g., splitting each interval into subintervals and comparing ακ) or a flux-rope model subtraction for a subset of events. Without such a check, the measured kurtosis may reflect non-turbulent large-scale structure rather than turbulence intermittency, and the radial invariance could be an artifact of averaging over heterogeneous intervals.","section":"§3.1–3.2"}],"minor_comments":[{"comment":"Please specify the fitting procedure used to obtain ακ and αPSD (log-log least squares, weighting, number of points per bin) and how the 'center of each bin' is defined.","section":"§3.2"},{"comment":"The sentence 'only events with average proton β of less than 0.5 are included' should state the time averaging window and data source used for β; Table 2 does not list β, so the reader cannot verify this selection criterion.","section":"§2"},{"comment":"For the example event, the ακ value of 0.40 is reported without an uncertainty, and the fitting range (10^-3 to 10^-1 Hz) differs from the range used in the statistical analysis; please give the fit uncertainty and clarify the relationship between the two ranges.","section":"§3.1"},{"comment":"The color scales for αPSD and |σc| are not described in the caption or in the text; please add color-bar labels and explain the five-point moving average used for the trend lines.","section":"Figure 4"},{"comment":"The phrase 'T able 2' should be corrected, and Table 2 should clarify that shock times are absent for some events rather than implying a complete shock-time column.","section":"Table 1 and Table 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and uses a sensible, standard methodology with a useful multi-spacecraft event sample. The main weakness is not the method itself but the statistical support for the headline invariance claim; this is fixable with additional inference and robustness tests. I therefore recommend major revision rather than rejection. The authors should also be encouraged to make the robustness checks quantitative, since the current text asserts insensitivity to the kdi range and running-average removal without showing the supporting numbers."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First statistical look at intermittency in ICMEs across 0.25–1 au, with 49 events from PSP and SolO. The kurtosis scaling exponent ακ is a standard measure, and the comparison to sheaths and ambient wind is a useful step forward. Credit where due: the authors use public data, list all events in the appendix, handle boundary regions by shaving 10% off intervals, and check running-average removal up to 4 hours. The observation that sheath ακ rises with distance while ICME ακ does not is genuinely interesting, even if not formally nailed down.\n\nThe soft spot is the interpretation. Table 1 gives Kendall τ = 0.11 between ακ and heliocentric distance for ICMEs, versus 0.30 for sheaths and upstream. No p-values, confidence intervals, or equivalence tests are reported. A non-significant correlation is not evidence of invariance. With 49 events, the approximate 95% CI on that 0.11 is –0.08 to 0.30, so a trend as strong as the sheath's is not excluded. The bimodal sampling (PSP mostly inside 0.6 au, SolO mostly outside) could also hide a spacecraft-dependent offset. The paper's conclusion that ICMEs are 'relatively static, well-developed turbulent environments' is therefore stronger than the statistics support. I would want bootstrap CIs on τ, an explicit equivalence test, or language that simply says 'no significant radial trend was detected.'\n\nThe kdi range choice (10^-2.5 to 10^-1.5) is post hoc, but they do report that other ranges give low correlations, so that is a minor concern. Also minor: the 'fully developed turbulence' interpretation is speculative; the null could just mean the ICME is isolated from the forcing that drives sheath evolution.\n\nWho does this matter for? People working on solar wind turbulence and CME propagation. The event list and method are reusable. It deserves peer review; a good referee should push on the uncertainty quantification and the wording of the invariance claim. My recommendation: send to review, expect a major revision or at least a careful rewrite of the abstract and conclusions.","headline":"A useful first statistical sample of ICME intermittency, but the radial-invariance claim needs significance testing or softer wording before it can stand.","tokens_in":16525,"tokens_out":2901,"would_cite":true,"duration_ms":30605,"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":"In coronal mass ejections, turbulence intermittency is radially invariant between 0.25 and 1 au, while sheath regions become more intermittent with distance.","keywords":["intermittency","interplanetary coronal mass ejections","magnetic flux rope","solar wind turbulence","kurtosis","structure functions","Parker Solar Probe","Solar Orbiter"],"falsifier":"Apply the same kurtosis analysis after subtracting a fitted smooth flux-rope model (for example a force-free cylindrical field) instead of a running average, and check whether ακ in ICMEs remains radially invariant; if the invariance vanishes or the exponent values shift systematically, the result is an artifact of the background-removal choice.","tokens_in":15591,"feed_emoji":"🌀","tokens_out":7306,"duration_ms":71290,"temperature":0.7,"pith_summary":"The paper asks whether the small-scale magnetic fluctuations inside interplanetary coronal mass ejections behave like developed turbulence and whether that turbulence changes as the ejections travel from 0.25 to 1 au. Using kurtosis and its scaling exponent in 49 ICMEs observed by Parker Solar Probe and Solar Orbiter, it finds that intermittency inside the ICMEs is radially invariant, while in the sheath regions ahead of the ICMEs intermittency increases with distance. The same scale-dependent kurtosis behavior appears in ICMEs, sheaths, and the surrounding solar wind, meaning ICME interiors are no less intermittent than the wind around them. The paper interprets the ICME result as evidence that these structures are relatively stable, fully developed turbulent environments, and the sheath result as evidence that sheath turbulence is still developing close to the Sun and may contain non-turbulent coherent structures that inflate kurtosis.","feed_headline":"ICME turbulence is already fully formed at 0.25 au","feed_subtitle":"49 ejections show flat kurtosis scaling to 1 au; sheath turbulence still grows with distance.","key_machinery":"The analysis is carried by scale-dependent kurtosis of magnetic-field increments, $\\kappa(\\tau) = S_4^B(\\tau)/(S_2^B(\\tau))^2$, computed from second- and fourth-order structure functions of the field magnitude; $\\kappa = 3$ marks a Gaussian distribution and higher values indicate intermittency. The intermittency level is quantified by $\\alpha_\\kappa$, the exponent of the power-law fit $\\kappa \\propto f^{\\alpha_\\kappa}$ (or $\\kappa \\propto k\\,d_i^{\\alpha_\\kappa}$) across bins in the inertial range, with steeper exponents meaning a more developed turbulent cascade. Spacecraft frequencies are converted to plasma-frame scales $k\\,d_i$ through Taylor's hypothesis, and a wavelet-based local intermittency measure (LIM) is used as a consistency check and to locate the fluctuations that drive the non-Gaussian tails. The comparison across radial distances is made by averaging event-wise $\\alpha_\\kappa$ in a fixed bin, $10^{-2.5} < k\\,d_i < 10^{-1.5}$, away from injection and kinetic scales.","core_discovery":"The central claim is that the level of intermittency in ICME interiors, measured by the scaling exponent $\\alpha_\\kappa$ of magnetic-field kurtosis in the MHD inertial range ($10^{-2.5} < k\\,d_i < 10^{-1.5}$), does not vary with heliocentric distance between 0.25 and 1 au. In the ICME sheath regions, by contrast, $\\alpha_\\kappa$ increases with distance. Kurtosis behaves similarly across all four interval types (upstream wind, sheath, ICME, downstream wind), with values well above the Gaussian level at small scales and a gradual return toward Gaussian statistics at larger scales; the average intermittency level is comparable in all intervals. Correlations between $\\alpha_\\kappa$ and distance, speed, spectral index, residual energy, cross helicity, and proton $\\beta$ are generally low, with the only clear radial trends appearing in upstream and sheath intervals. The authors conclude that ICMEs are relatively static, well-developed turbulent environments, while sheaths are younger structures whose turbulence is not yet fully developed at small heliocentric distances; unusually high absolute kurtosis with low scaling exponent in some sheaths may reflect non-turbulent structures rather than an evolved cascade.","pith_inferences":["Editorial inference: if ICME turbulence is already fully developed at 0.25 au, particle-acceleration and cosmic-ray transport models may be able to use a fixed fluctuation-scattering prescription inside ejecta rather than a distance-dependent one.","Editorial inference: the sheath trend predicts even lower intermittency at distances below 0.25 au, so close perihelion passes could test whether the increase continues or saturates.","Editorial inference: separating coherent structures from turbulence in sheaths (for example by thresholding the local intermittency measure or checking phase coherence) on the same events would test whether the high-κ, low-ακ population is truly non-turbulent.","Editorial inference: re-running the analysis separately for magnetic clouds and complex ejecta would show whether radial invariance is generic to ICMEs or specific to smooth flux ropes."],"forward_implications":["ICME interiors can be treated as statistically stationary turbulence over 0.25–1 au, so observations at different heliocentric distances can be merged when studying their small-scale fluctuations.","The absence of radial evolution means the turbulent state of an ICME is established early and is not reset by expansion or interaction with the surrounding solar wind.","Sheath turbulence is still developing near the Sun, so radial distance must be accounted for when comparing sheath intermittency between events or missions.","In sheaths, absolute kurtosis and scaling exponent can disagree, so studies that report only kurtosis values may overstate the development of the turbulent cascade.","Because intermittency levels are similar across ICMEs, sheaths, and ambient wind, intermittency alone cannot identify an ICME interval in the magnetic-field data."],"supporting_citations":[{"why":"Defines structure functions and kurtosis as standard tools for quantifying intermittency in turbulence.","marker":"Frisch (1995)"},{"why":"Earlier case study of kurtosis in ICMEs whose scaling exponents and behavior this paper reproduces and extends statistically.","marker":"Sorriso-Valvo et al. (2021)"},{"why":"Second case study of ICME intermittency providing comparison values for kurtosis and ακ.","marker":"Márquez Rodríguez et al. (2023)"},{"why":"Established that slow solar wind shows no radial evolution of intermittency, the analogue used to interpret the ICME invariance as fully developed turbulence.","marker":"Bruno et al. (2003)"},{"why":"Reported related indications of the static, non-evolving inertial-range properties of ICMEs.","marker":"Good et al. (2023)"},{"why":"Provided the approach of using ακ, the kurtosis scaling exponent, as the intermittency measure.","marker":"Sorriso-Valvo et al. (2018)"},{"why":"Applied the same ακ-based comparison of intermittency between solar wind intervals.","marker":"Telloni et al. (2021)"},{"why":"Supplies the ICME event catalog with boundaries and shock times used to build the 49-event sample.","marker":"Möstl et al. (2020)"},{"why":"Documents the FIELDS magnetic field data from Parker Solar Probe used for the fluctuation analysis.","marker":"Bale et al. (2016)"},{"why":"Documents the MAG magnetometer data from Solar Orbiter used for the fluctuation analysis.","marker":"Horbury et al. (2020)"}],"fun_headline_variants":["ICME interiors show mature turbulence from 0.25 au","ICMEs: no radial change in turbulence, sheaths still growing","Sheath turbulence evolves with distance; ICME interiors stay constant","Kurtosis scaling shows ICMEs are fully turbulent already at 0.25 au"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion rests on assuming that, after subtracting a 1–4 hour running average, the residual magnetic-field fluctuations inside each ICME are stationary, homogeneous turbulence; if flux-rope rotation, interval boundaries, or non-turbulent coherent structures contaminate the residual, the measured kurtosis and its radial invariance would not describe turbulent intermittency.","fun_headline_variants_meta":{"raw":{"variants":["ICME interiors show mature turbulence from 0.25 au","ICMEs: no radial change in turbulence, sheaths still growing","Sheath turbulence evolves with distance; ICME interiors stay constant","Kurtosis scaling shows ICMEs are fully turbulent already at 0.25 au"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000595,"raw_usage":{"total_tokens":2816,"prompt_tokens":1005,"completion_tokens":1811,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":621,"completion_tokens_details":{"reasoning_tokens":1733}},"tokens_in":621,"tokens_out":1811,"duration_ms":14040,"temperature":1.0,"reasoning_tokens":1733,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:10:30.012690+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Apply the same kurtosis analysis after subtracting a fitted smooth flux-rope model (for example a force-free cylindrical field) instead of a running average, and check whether ακ in ICMEs remains radially invariant; if the invariance vanishes or the exponent values shift systematically, the result is an artifact of the background-removal choice.","supporting_citations":[{"cited_title":"1995, Turbulence: The Legacy of A","cited_arxiv_id":null,"evidence_quote":"Defines structure functions and kurtosis as standard tools for quantifying intermittency in turbulence."}],"review_version":1}