{"id":"a2301b8a-7922-44c0-add7-fc943bb3a24b","arxiv_id":"2505.04450","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Strong interplanetary shocks make downstream solar wind turbulence more balanced and more magnetically dominated, while quasi-parallel shocks have the most Alfvénic upstream regions.","lead":"This study analyzes 378 interplanetary shock waves measured by the Wind and Solar Orbiter spacecraft to see how they change the turbulence of the surrounding solar wind. Strong shocks make downstream fluctuations more balanced between sunward and antisunward waves and more dominated by magnetic energy, which matters for how particles are accelerated at shocks.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"AF/SFR occurrence trends in Figs. 8–9 are derivative of the σ_c/σ_r PDF shifts and rest on unvalidated hand-picked thresholds; a threshold sweep or independent structure validation is needed to confirm the central claim.","rationale":"The reader correctly identified the hand-chosen AF/SFR thresholds as the weakest assumption. I agree, and I would sharpen the point: the AF/SFR trends are not independent corroboration because they use the same parameters whose PDFs are claimed to shift. The AF criteria target the high-|σ_c|, near-equipartition region (which shrinks downstream), while the SFR criteria target the low-|σ_c|, magnetically dominated region (which grows downstream). Therefore the occurrence statistics are, to first order, a relabeling of the σ_c/σ_r PDF changes. What would make the AF/SFR results physically meaningful is a demonstration that the selected intervals correspond to genuine flux ropes and Alfvénic fluctuations; the paper provides no such validation. This is a testable concern, not a demonstration of error: a threshold sweep and an independent structure-identification check would settle it. The underlying σ_c/σ_r PDF comparisons remain plausible and are supported by a large public dataset, so the concern does not justify rejecting the paper. It does justify keeping the CONDITIONAL verdict pending the requested validation, so I do not change the reader's verdict.","tokens_in":25968,"tokens_out":10397,"duration_ms":109363,"concrete_test":"Run a threshold-sensitivity analysis on the Fig. 8–9 pipeline: vary the |σ_c| threshold over 0.7–0.95, the |σ_r| threshold over 0.2–0.5, and the |σ_m| threshold over 0.5–0.8, then recompute the upstream/downstream AF and SFR occurrence differences and their dependence on rg, βu, ΔV, and θBn. On a random subset of about 50 shocks, additionally compare with an independent structure detector (e.g., minimum-variance or rotation-angle analysis for flux ropes and a Walén test for Alfvén waves). If the sign and parameter scaling of the trends are stable across the threshold grid and reproduced by the independent detector, the concern is resolved; if not, the AF/SFR claims should be downgraded.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The AF/SFR results (Sect. 3.3.3, Figs. 8–9) are presented as independent physical-structure evidence, but they are not independent of the σ_c/σ_r PDF changes that form the headline claim. The AF criteria (|σ_m|>0.7 or >0.3, |σ_c|>0.9, |σ_r|<0.3) select exactly the parameter-space region that shrinks downstream, while the SFR criteria (|σ_m|>0.7, |σ_c|<0.4, |σ_r|<−0.5) select the region that grows downstream. Thus a decrease in AFs and an increase in SFRs will occur for almost any fixed thresholds once the reported PDF shifts are present, so the occurrence trends do not by themselves demonstrate that the selected intervals correspond to real Alfvénic fluctuations and flux ropes. The thresholds are adopted from prior work without validation in shock-modified turbulence, where compressive fluctuations and current sheets may produce high |σ_m| and low |σ_r| without being flux ropes. In addition, the criterion text in Sect. 3.3.3 is internally ambiguous for linearly polarized waves, which have near-zero σ_m, yet appears to require |σ_m|>0.3 for their identification. If the thresholds are mis-calibrated, the downstream increase in SFRs, one of the paper's headline statements, could be an artifact of the PDF shifts rather than a real increase in flux-rope structures.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript analyzes normalized cross-helicity, residual energy, and magnetic helicity in 1-h intervals upstream and downstream of 371 fast forward interplanetary shocks detected by Wind over 1995-2023 and 7 shocks observed by Solar Orbiter below 0.5 au. For each shock it computes the change across the shock and compares PDFs and occurrence statistics in subsets defined by the 20th and 80th percentiles of gas compression ratio, upstream plasma beta, velocity jump, and shock angle. The headline results are that strong shocks (large velocity jump and gas compression ratio) produce more balanced cross-helicity and more magnetically dominated fluctuations downstream, that the occurrence of Alfvénic-fluctuation periods decreases while small-scale-flux-rope periods increase across the shock, and that these trends depend on the shock parameters; the near-Sun sample is reported to show similar behavior.","tokens_in":26289,"tokens_out":5562,"duration_ms":53446,"significance":"If the claimed trends hold, this would be the first large-sample demonstration that shock strength and obliquity systematically control the reprocessing of solar-wind turbulence at interplanetary shocks, with implications for turbulence-driven particle acceleration at collisionless shocks. The paper's strengths are its large, public-data-based sample spanning 2.5 solar cycles, its use of publicly catalogued shocks and the SerPyShocks package, and its transparent presentation of median trends and PDFs that can be reproduced and extended. The main claims are physically plausible and consistent with prior case and small-sample studies; however, the absence of significance testing and the dependence of the Alfvénic-fluctuation/flux-rope component on unvalidated fixed thresholds currently prevent the paper from fully supporting the stronger structural interpretation in the abstract.","major_comments":[{"comment":"The conclusions that PDF differences are 'significant' and that trends with r_g, beta_u, Delta V, and theta_Bn are 'clear' are based on visual inspection of PDFs and 40-event running medians; no significance tests are reported. With 371 events split into many subsets and with the large event-to-event scatter visible in Fig. 4, sampling fluctuations could produce apparent PDF differences of the size shown. I ask for a quantitative assessment, for example two-sample Kolmogorov-Smirnov or Anderson-Darling tests between the upstream and downstream PDFs, bootstrap confidence intervals for the running medians, or permutation tests for the PDF differences, together with a statement of how many events drive each trend.","section":"§3.3.2, Figs. 5-7"},{"comment":"The AF and SFR occurrence trends are presented as independent physical-structure evidence, but the fixed criteria (|sigma_m|>0.7 or >0.3, |sigma_c|>0.9, |sigma_r|<0.3 for AFs; |sigma_m|>0.7, |sigma_c|<0.4, |sigma_r|<-0.5 for SFRs) select exactly the parameter-space regions that shift between the upstream and downstream PDFs in Figs. 6-7. Consequently, a decrease in AF occurrence and an increase in SFR occurrence would be expected for almost any fixed thresholds once the reported sigma_c and sigma_r shifts are present, so the occurrence trends do not by themselves validate that the selected intervals correspond to real Alfvénic fluctuations and small-scale flux ropes. Because the thresholds are adopted from Zhao et al. (2021) and Ruohotie et al. (2022) without validation in shock-modified turbulence, I request either a threshold sweep showing that the trends are robust, or a comparison with an independent structure-identification method, to support the abstract's claim that AF populations decrease and SFR populations increase across the shock.","section":"§3.3.3, Figs. 8-9 and Table 1"},{"comment":"The subset definitions are chosen post hoc as the 20th and 80th percentiles of each shock parameter, explicitly 'to give sufficiently distinct populations' (§3.2). All parameter-dependence claims in the abstract and in Figs. 6-9 are therefore conditioned on this arbitrary split. The paper should demonstrate that the conclusions are robust to the choice of percentile - for example by testing quartiles or terciles, or by performing a continuous regression of the upstream-to-downstream changes in sigma_c and sigma_r on the shock parameters - and should report the number of events in each extreme subset.","section":"§3.2, Figs. 3, 6-7"}],"minor_comments":[{"comment":"The sentence 'For AFs, we require that |sigma_m|>0.7 or |sigma_m|>0.3 for circular or linearly polarized waves' is internally ambiguous, since linearly polarized Alfvén waves have sigma_m approximately zero; please clarify which modes satisfy which threshold and whether the |sigma_m|>0.3 condition is really intended for linear polarization.","section":"§3.3.3"},{"comment":"The entry '1 .36±1.22' contains an errant space and should read '1.36±1.22'.","section":"Table 2"},{"comment":"The x-axis ranges in the main figure are truncated to 'exclude a few outlier points'; please state explicitly how many events are excluded in each panel, since the running medians near the plot boundaries may depend on the truncation.","section":"§3.3.1, Fig. 4"},{"comment":"With only seven Solar Orbiter shocks, the statement that 'the results show an overall similar behaviour as at 1 au' should be explicitly framed as exploratory and accompanied by the individual event values listed in Table 2 rather than only by combined PDFs.","section":"§3.4, Fig. 10"},{"comment":"The citation 'Pitna et al. (2023)' appears with inconsistent diacritics ('Pitna' versus 'Pitˇna') in the text and reference list.","section":"References"},{"comment":"The caption says '80th and 20th quartiles' but the text and figure markers refer to percentiles; please correct the wording.","section":"Fig. 3 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is a competent and useful statistical study for the journal's scope, and I do not see grounds for rejection. My main concern is that the Alfvénic-fluctuation and small-scale-flux-rope section is being presented as independent structural evidence when it is, in effect, a relabeling of the same parameter-space regions whose PDFs shift across the shock; the requested threshold-sensitivity or validation analysis should be feasible within a revision. I would also emphasize the need for significance testing before the headline parameter-dependence claims can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know. The central result — that strong shocks with large compression ratios and velocity jumps make downstream turbulence more balanced in cross-helicity and more magnetically dominated in residual energy — looks solid and is a genuine extension of the Borovsky/Soljento/Good work. The AF/SFR occurrence trends in Figs. 8–9, however, are not an independent line of evidence. They are the same sigma_c/sigma_r/sigma_m PDF shifts re-encoded through thresholds inherited from earlier papers, and those thresholds are not validated for shock-modified turbulence.\n\nWhat's new: a systematic stratification of inertial-range turbulence parameters by shock gas compression ratio, upstream beta, velocity jump, and obliquity, on 371 Wind shocks. That's a meaningful step beyond the earlier aggregate studies. The methods are careful: public data, public shock catalogues and fitting tools, clear wavelet parameters, rectified cross-helicity to the Parker spiral, and honest presentation of scatter via running medians. The paper writes plainly and acknowledges prior work properly.\n\nSoft spots, in proportion. The biggest is the absence of significance tests: the PDF comparisons in Figs. 6–7 have no error bars, KS tests, or bootstrap intervals, and the 20th/80th percentile subset definition is post hoc. The main trends are visually strong, but a referee should ask for a robustness check on the percentile choice and some formal test on the PDF differences. The AF/SFR analysis inherits thresholds from Zhao et al. and Ruohotie et al. and applies them to bin counts in the very same parameter space whose PDFs are shifting. A decrease in AF-like bins and increase in SFR-like bins is what any fixed thresholds would do once the PDFs move as shown. That doesn't prove the intervals correspond to real Alfvénic fluctuations and flux ropes downstream. A threshold sweep or independent structural validation would be enough to close this. Minor: the criterion text for linearly polarized AFs requiring |sigma_m|>0.3 is ambiguous since those waves have sigma_m near zero. And the Solar Orbiter part rests on seven events; the authors frame it as qualitative, which is fair, but the wording in the abstract (\"similar distributions\") is a bit strong.\n\nNone of these are load-bearing. The direct PDF evidence for the shock-parameter dependence stands on its own, and the paper is a clear step forward for shock-turbulence interaction studies.\n\nWho this is for: anyone working on solar wind turbulence at shocks, particle-scattering models at interplanetary shocks, or statistical heliospheric data analysis. I'd use it in a reading group and would cite it for the shock-parameter stratification, with a caveat about the structure counts.\n\nRecommendation: send it to peer review. It deserves refereeing; expect revisions on significance testing and threshold sensitivity, not rejection.","headline":"A large, careful statistical study whose core claim about shock-parameter control of sigma_c and sigma_r holds up; the AF/SFR structure counts are derivative of the same PDF shifts and need threshold validation.","tokens_in":26878,"tokens_out":3278,"would_cite":true,"duration_ms":32723,"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":"Strong shocks rebalance solar wind turbulence, turning Alfvénic upstream fluctuations into balanced, magnetically dominated downstream ones.","keywords":["normalized cross-helicity","residual energy","magnetic helicity","interplanetary shocks","solar wind turbulence","Alfvénic fluctuations","small-scale flux ropes","Solar Orbiter"],"falsifier":"Take the same 371 Wind shocks and label intervals of known content by an independent method, such as magnetic flux-rope reconstruction or particle-in-cell simulations of shock crossing with embedded coherent structures, and compare whether the $\\sigma_c$, $\\sigma_r$, $\\sigma_m$ thresholds recover the known structures; alternatively, recalculate the AF and SFR occurrence trends with thresholds varied over a plausible range and see whether the upstream-to-downstream decrease and increase survive.","tokens_in":25783,"feed_emoji":"🛰️","tokens_out":7458,"duration_ms":66166,"temperature":0.7,"pith_summary":"The paper aims to show that interplanetary shock waves do not merely compress the solar wind but actively rewrite the turbulence that passes through them, and that the rewriting is controlled by the shock's strength and orientation. Using 371 shocks measured by Wind at 1 au and seven measured by Solar Orbiter closer to the Sun, it argues that strong shocks with large velocity jumps and density compressions turn imbalanced, Alfvénic upstream fluctuations into balanced, magnetically dominated downstream fluctuations, while weak shocks leave the turbulence essentially unchanged. It further claims that shock passage suppresses periods identified as Alfvénic fluctuations and increases periods identified as small-scale flux ropes. These changes matter because the balance between sunward and antisunward waves, and between magnetic and kinetic fluctuation energy, sets the scattering environment that governs how efficiently shocks accelerate charged particles.","feed_headline":"371 shocks: strong shocks rebalance solar wind turbulence","feed_subtitle":"Alfvénic periods drop and flux-rope periods rise after strong shocks, across 1 au and closer to the Sun.","key_machinery":"The analysis is carried by three normalized turbulence parameters computed from wavelet spectrograms in the 1.67–16.7 mHz inertial-range band: cross-helicity $\\sigma_c$ (the normalized difference in power between fluctuations $z^\\pm = v \\pm b$ propagating along and against the mean magnetic field), residual energy $\\sigma_r$ (normalized difference between kinetic and magnetic fluctuation energy), and magnetic helicity $\\sigma_m$. These are compared upstream and downstream of each shock and binned by four shock parameters: gas compression ratio, upstream plasma $\\beta$, velocity jump, and shock angle. The structural interpretation rests on threshold-based identification of Alfvénic fluctuations (requiring $|\\sigma_m| > 0.7$ or $> 0.3$, $|\\sigma_c| > 0.9$, $|\\sigma_r| < 0.3$) and small-scale flux ropes (requiring $|\\sigma_m| > 0.7$, $|\\sigma_c| < 0.4$, $|\\sigma_r| < -0.5$), applied to the same wavelet spectrograms.","core_discovery":"The central claim, drawn from wavelet-derived inertial-range values of the normalized cross-helicity $\\sigma_c$, residual energy $\\sigma_r$, and magnetic helicity $\\sigma_m$ in 1-h windows around each shock, is that the shock transition systematically reprocesses solar wind turbulence, with the effect controlled by shock parameters. Shocks with large gas compression ratios and large velocity jumps produce downstream fluctuations that are considerably more balanced in cross-helicity ($\\sigma_c$ near 0) and more magnetically dominated (more negative $\\sigma_r$) than upstream, and the rectified cross-helicity shows this balance comes primarily from a strong reduction of antisunward-propagating fluctuation power. Across the shock, the occurrence of time periods meeting the adopted criteria for Alfvénic fluctuations decreases on average, while periods meeting the criteria for small-scale flux ropes increase, with the Alfvénic periods peaking for quasi-parallel shocks with large velocity jumps and small upstream $\\beta$, and flux-rope periods increasing with gas compression ratio and upstream $\\beta$. Magnetic helicity averages near zero and shows little change across the shock. The same upstream-to-downstream trends appear in the seven Solar Orbiter shocks below 0.5 au, indicating the effect is not confined to 1 au.","pith_inferences":["A testable extension the paper does not perform: if the downstream flux-rope increase comes from shock-triggered reconnection of current sheets, the increase should correlate with the density of current sheets in the upstream plasma, which can be counted from the same Wind data.","The AF decrease downstream may be partly a selection effect of the threshold criteria, because the downstream shift in $\\sigma_r$ toward more negative values automatically pushes intervals out of the AF box; a control analysis with scale-dependent or locally normalized thresholds could separate selection from a true structural change.","The rectified cross-helicity result—that the antisunward population is strongly reduced downstream of strong shocks—suggests either preferential damping of antisunward waves or generation of sunward waves; a direct examination of the sunward and antisunward fluctuation power spectra in the downstream would discriminate between these two mechanisms.","If the shock-angle dependence is confirmed with more quasi-parallel events, the observed trends imply that the foreshock region itself, not just the downstream, should show turbulence parameters organized by shock geometry, which could be tested with multi-spacecraft crossings of the same shock."],"forward_implications":["If strong shocks systematically balance cross-helicity and steepen magnetic dominance downstream, particle-acceleration models at shocks must use shock-dependent turbulence inputs rather than a universal upstream spectrum.","Because weak shocks leave the turbulence parameters nearly unchanged, their efficiency at accelerating particles should be governed by the pre-existing upstream turbulence rather than by shock-generated fluctuations.","The downstream increase in small-scale flux ropes implies that shock passage creates or amplifies coherent magnetic structures, which can act as scattering centres or reservoirs for energetic particles.","The similarity of trends at 0.3–0.5 au and at 1 au suggests these shock-processing effects are generic across the inner heliosphere, not an artefact of solar-wind evolution.","Quasi-parallel shocks with large velocity jumps and low upstream beta are the configuration that most enriches upstream turbulence with Alfvénic fluctuations, which is exactly the condition under which shock self-generated waves should be strongest."],"supporting_citations":[{"why":"Establishes the prior statistical result that Alfvénicity decreases from upstream to downstream at shocks, which this work extends to the shock-parameter dependence.","marker":"Borovsky (2020)"},{"why":"Provides the upstream/downstream cross-helicity and residual energy distributions for ICME-driven shocks that this study's PDFs are compared against.","marker":"Soljento et al. (2023)"},{"why":"Superposed epoch analysis of cross-helicity in CME sheaths that supports the trend toward more balanced downstream turbulence.","marker":"Good et al. (2022)"},{"why":"Case study showing balanced cross-helicity and enhanced magnetic helicity at an ICME shock; source of the AF/SFR criteria approach.","marker":"Zhao et al. (2021)"},{"why":"Supplies the statistical finding of more small-scale flux ropes in CME sheaths and the specific AF/SFR threshold criteria adopted here.","marker":"Ruohotie et al. (2022)"},{"why":"Theoretical basis for how cross-helicity influences scattering-centre compression and first-order Fermi acceleration at shocks.","marker":"Vainio and Schlickeiser (1998)"},{"why":"Provides the heliospheric shock database from which the Wind events are drawn, defining the analysed population.","marker":"Kilpua et al. (2015)"},{"why":"Supplies the shock-parameter estimation routine (mixed-mode normal determination) used to compute shock angle and speed jumps.","marker":"Trotta et al. (2022)"}],"fun_headline_variants":["Shock waves rebalance solar wind turbulence","Strong shocks rebalance turbulence, shift Alfvénic to flux ropes","Shock strength controls turbulence rebalancing","371 shocks reveal turbulence shift at shocks","Shocks flip Alfvénic to flux rope periods"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire Alfvénic-fluctuation and flux-rope trend analysis depends on the hand-chosen numerical thresholds that decide whether a time interval counts as one structure or the other, so if those thresholds do not cleanly separate the intended physical structures in shock-modified turbulence, the reported occurrence trends would not reflect real structure populations.","fun_headline_variants_meta":{"raw":{"variants":["Shock waves rebalance solar wind turbulence","Strong shocks rebalance turbulence, shift Alfvénic to flux ropes","Shock strength controls turbulence rebalancing","371 shocks reveal turbulence shift at shocks","Shocks flip Alfvénic to flux rope periods"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001077,"raw_usage":{"total_tokens":4607,"prompt_tokens":1142,"completion_tokens":3465,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":758,"completion_tokens_details":{"reasoning_tokens":3391}},"tokens_in":758,"tokens_out":3465,"duration_ms":24083,"temperature":1.0,"reasoning_tokens":3391,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:28:14.753706+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the same 371 Wind shocks and label intervals of known content by an independent method, such as magnetic flux-rope reconstruction or particle-in-cell simulations of shock crossing with embedded coherent structures, and compare whether the $\\sigma_c$, $\\sigma_r$, $\\sigma_m$ thresholds recover the known structures; alternatively, recalculate the AF and SFR occurrence trends with thresholds varied over a plausible range and see whether the upstream-to-downstream decrease and increase survive.","supporting_citations":[],"review_version":1}