{"id":"4828e9bc-665b-4a54-973d-5129e958cc70","arxiv_id":"2607.11142","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Persistent seismic swarms show strong anomalous multiscaling of interevent distances, while the Landers sequence saturates at high moments, suggesting a diagnostic migration fingerprint.","lead":"Ten Southern California earthquake swarms share a multiscaling migration signature: short and long interevent distances evolve under different effective laws, unlike the Landers tectonic sequence. That contrast is offered as a statistical fingerprint of swarm-style diffusion.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.5","headline":"The fingerprint claim rests on one tectonic baseline (Landers); without more mainshock–aftershock sequences the nonlinear-vs-saturation contrast may be sequence-specific rather than a class diagnostic.","rationale":"The reader’s weakest_assumption is precisely the load-bearing concern: one tectonic baseline is insufficient to underwrite a diagnostic class boundary. The empirical pattern among the ten swarms is clear and the methods (consensus clustering, fixed temporal window across q, AICc model selection) are transparent and reproducible via the Zenodo code, so the observational claim that these swarms show strong anomalous diffusion is solid. What remains conditional is the stronger language that ε(q) “distinguishes swarm migration from conventional tectonic earthquake sequences.” That language requires the Landers saturation to be representative. Expanding the tectonic sample is the single concrete check that settles the issue; until it is done the CONDITIONAL verdict is appropriate and no adjustment is needed. No other internal inconsistency (e.g., in the multifractal formalism or moment definitions) rises to the same load-bearing level.","tokens_in":10721,"tokens_out":598,"duration_ms":6856,"concrete_test":"Apply the identical multiscaling pipeline (same q range, r_ij > 0.05 km cut, common temporal-interval criteria of ≥6 bins / ≥0.7 decades / ≥50 pairs, AICc polynomial selection) to at least two additional large SCEDC tectonic sequences (e.g., Hector Mine 1999 and Ridgecrest 2019 aftershock catalogs). If either sequence yields a nonlinear positive branch comparable to the swarms rather than high-order saturation, the class-boundary fingerprint claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that ε(q) is a quantitative fingerprint distinguishing swarm migration from conventional tectonic sequences. That claim requires the observed contrast—nonlinear positive branch for all ten swarms versus high-order saturation for Landers (Fig. 3; §3)—to mark a class boundary, not a catalog- or sequence-specific difference. The paper supplies only a single tectonic baseline (Landers), chosen because it is large and previously studied by Godano & Pingue (2005). No other mainshock–aftershock sequences (e.g., Hector Mine, Ridgecrest, or other SCEDC tectonic clusters of comparable size) are analyzed with the same moment-scaling pipeline, common temporal-interval selection, and AICc polynomial protocol. Consequently the diagnostic interpretation in the Abstract and Conclusions is under-supported: saturation could be a Landers-specific feature of a long, quasi-linear fault projection rather than the generic signature of tectonic sequences. The reader correctly flags this as the weakest assumption; it is load-bearing because every stronger claim about “distinguishing” and “fingerprint” collapses if additional tectonic sequences also produce nonlinear ε(q).","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript identifies ten persistent seismic swarms in the relocated Southern California catalogue via a consensus DBSCAN-like clustering procedure and analyzes their spatial migration with the multiscaling spectrum ε(q) of interevent-distance moments M_q(θ). All ten swarms show an approximately linear ε(q) for negative q and a systematically nonlinear positive-q branch, interpreted as strong anomalous diffusion in which small and large separations evolve with different effective exponents. By contrast, the Landers mainshock–aftershock sequence exhibits high-order saturation of ε(q), attributed to spatially bounded (caged) diffusion along the fault projection. The authors conclude that ε(q) supplies a quantitative fingerprint distinguishing swarm migration from conventional tectonic sequences.","tokens_in":11095,"tokens_out":1278,"duration_ms":18304,"significance":"If the reported contrast generalizes, the multiscaling spectrum would be a useful, catalog-based diagnostic for swarm versus tectonic migration, complementing moment-release shape criteria and classical single-exponent diffusion models. Strengths include: (i) a clearly stated moment definition and shared temporal fitting window across q (avoiding artificial curvature); (ii) AICc-based selection of the positive-branch polynomial; (iii) open Zenodo code for swarm identification and multiscaling; and (iv) a reproducible consensus clustering workflow that reduces dependence on a single (τ, r0, ε) choice. The qualitative split between negative- and positive-order branches is visible across the ten swarm panels and is of genuine interest for anomalous-diffusion studies of seismicity.","major_comments":[{"comment":"Abstract, §3 (Fig. 3), and Conclusions: the claim that ε(q) is a fingerprint “capable of distinguishing swarm migration from conventional tectonic earthquake sequences” rests on a single tectonic baseline (Landers). Saturation may reflect Landers’ long, quasi-linear fault geometry rather than a generic tectonic signature. At least one or two additional large SCEDC mainshock–aftershock sequences (e.g., Hector Mine, Ridgecrest, or other N≳200 tectonic clusters) should be processed with the same pipeline (shared temporal interval, AICc protocol, r_ij cut) so that the nonlinear-vs-saturation contrast can be tested as a class boundary rather than a sequence-specific difference.","section":null},{"comment":"§2 Methods: only 10 of 95 persistent swarms (those with N>200) enter the multiscaling analysis. The central “all analyzed swarms” claim is therefore conditioned on a size cut that may select the most spatially extended or longest-lived systems. The manuscript should either (a) show that smaller swarms with adequate pair counts in the admissible θ windows yield the same qualitative ε(q) shape, or (b) explicitly restrict the fingerprint claim to large, persistent swarms and discuss possible size bias.","section":null},{"comment":"§3 and Fig. 2, Swarm 043: this sequence is flagged as a “notable exception” yet is still counted among the ten that exhibit the common signature. Its negative-order slope (H_−≈0.73) and positive-branch shape differ markedly from the other panels. Either reclassify 043 (and restate “all” claims accordingly) or provide a quantitative criterion for when a spectrum is considered consistent with the swarm fingerprint versus Landers-like saturation.","section":null},{"comment":"§3, selection of the common temporal scaling interval: candidate windows must span ≥0.7 decades, ≥6 bins, and ≥50 pairs/bin, with the “best overall linear scaling” chosen for a representative set of q. Because the same interval is then imposed on all q, the procedure can still favor intervals that look linear at low |q| while leaving high-q curvature under-constrained. Report sensitivity of ε(q) (especially the positive branch and AICc degree) to alternative admissible windows, or show that the nonlinear-vs-linear split is stable across the admissible set.","section":null}],"minor_comments":[{"comment":"Eq. (1) and surrounding text: the CV formula is typeset as “Cy — bu”; restore standard notation c_v = σ/μ and define symbols cleanly.","section":null},{"comment":"Eq. (2)–(5): the space–time metric, membership probabilities p_i and p_ij, and Mo conversion are incompletely or garbled in the text (missing powers, incomplete formulas). Provide complete, self-contained expressions.","section":null},{"comment":"Figure 2: independent vertical scales make cross-swarm comparison of H_− and high-q curvature difficult; consider a common y-range or a summary table of H_−, AICc degree, and ε(q) at fixed q (e.g., q=2, 5, 10).","section":null},{"comment":"Plain Language Summary and Key Points: useful, but “classical tectonic earthquake sequence” should be qualified as “the Landers sequence” until more baselines are shown.","section":null},{"comment":"References: Passarelli et al. (2026) and Ross et al. (2025) appear with overlapping titles/DOIs in the list; verify bibliographic entries.","section":null},{"comment":"Open Research: Zenodo DOI is given with a space (“zenodo. 21297044”); fix the link for reproducibility.","section":null}],"recommendation":"major_revision","confidential_remarks":"The empirical swarm multiscaling result is interesting and the methods are largely transparent; the main risk is overclaiming a class diagnostic from one tectonic control. If the authors add 1–2 more tectonic sequences and clarify Swarm 043 / N>200 scope, this could be a solid contribution. Fit for a geophysics / seismology journal is appropriate; novelty is incremental relative to Godano & Pingue (2005) but the swarm-focused application and consensus clustering are worthwhile."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The useful takeaway is simple: ten large, consensus-selected Southern California swarms share a qualitative ε(q) shape—roughly linear for negative moments, systematically nonlinear for positive ones—while Landers saturates at high q. That is a clean observational contrast, not a re-derivation of multiscaling itself.\n\nWhat is new is the application, not the formalism. Godano & Pingue (2005) already put multiscaling of interevent distances on the map and already saw Landers-like caging. Here they add a careful consensus DBSCAN-style swarm catalog (Passarelli-style moment-release criteria plus membership/co-membership thresholds), restrict to N>200, fix a common temporal window across all q, and use AICc for the positive-branch polynomial. The figures make the shared swarm pattern easy to see, code is on Zenodo, and the moment definition is standard. That is real, reproducible work.\n\nThe soft spot is load-bearing for the language they chose. “Fingerprint capable of distinguishing swarm migration from conventional tectonic sequences” needs more than one tectonic baseline. Landers is large and familiar, but it is also a long, quasi-linear fault projection; saturation could be that geometry rather than the generic tectonic signature. Hector Mine, Ridgecrest, or other SCEDC mainshock–aftershock clusters of comparable size are not run through the same pipeline. Swarm 043 is already an admitted exception, only 10 of 95 swarms meet the size cut, and the shared-interval selection still has discretion. None of that kills the empirical pattern; it just means the class-boundary claim is ahead of the evidence.\n\nMath and citations look honest: ε(q) is estimated, not forced; self-citation is to prior multiscaling and completeness work that the paper actually uses. Free parameters are many but mostly declared.\n\nThis is for people who already care about fluid-driven vs stress-driven clustering and want a statistical handle beyond single diffusion exponents. I would bring it to reading group as a methods-plus-pattern paper, cite the swarm spectra if I am writing on anomalous diffusion or swarm diagnostics, and send it to peer review. Referees should demand more tectonic baselines before the fingerprint wording sticks; the core observation is still worth the time.","headline":"Solid empirical multiscaling contrast for ten consensus swarms vs Landers, but the “fingerprint” claim is one tectonic baseline short of a class diagnostic.","tokens_in":11684,"tokens_out":553,"would_cite":true,"duration_ms":7049,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Earthquake swarms share a multiscaling signature that separates their migration from ordinary tectonic sequences.","keywords":["earthquake swarms","anomalous diffusion","multiscaling spectrum","interevent distances","seismic migration","Southern California","tectonic sequences"],"falsifier":"Compute ε(q) for several other large, well-relocated tectonic mainshock–aftershock sequences of comparable size; if those sequences also show a nonlinear positive-q branch without high-order saturation, the claimed fingerprint fails as a diagnostic boundary between swarms and conventional tectonic sequences.","tokens_in":11603,"feed_emoji":"🌍","tokens_out":867,"duration_ms":24555,"temperature":0.7,"pith_summary":"Seismic swarms migrate in space, but ordinary single-exponent diffusion models have not captured how that migration works. This paper identifies ten persistent swarms in the relocated Southern California catalogue and tracks the statistical moments of distances between events as time lag grows. In every swarm the multiscaling spectrum is roughly linear for negative moments (short separations) yet systematically nonlinear for positive moments (long separations), meaning small and large interevent distances obey different effective scaling laws—strong anomalous diffusion. The Landers tectonic sequence instead saturates at high moments, as if epicenters remain spatially caged. The authors therefore present the multiscaling spectrum itself as a quantitative fingerprint that can distinguish swarm-style migration from conventional tectonic sequences, giving catalogue-based work a clearer statistical handle on fluid-driven or transient processes.","feed_headline":"Swarm quakes leave a multiscaling fingerprint tectonic ones lack","feed_subtitle":"Ten California swarms show short and long separations evolve under different laws; Landers saturates instead","key_machinery":"The multiscaling spectrum ε(q), defined as the scaling exponents of the temporal moments of interevent epicentral distances, M_q(θ) ~ θ^{ε(q)}. Linearity of ε(q) would imply a single diffusion exponent; the observed break between a linear negative-q branch and a nonlinear positive-q branch is the object that establishes strong anomalous diffusion and the swarm-versus-tectonic contrast.","core_discovery":"All ten persistent earthquake swarms display the same multiscaling spectrum ε(q): approximately linear for negative moments and systematically nonlinear for positive moments, so that short and long interevent distances evolve under different effective scaling laws and swarm migration is strong anomalous diffusion. By contrast, the Landers tectonic sequence shows high-order saturation of ε(q), consistent with spatially bounded diffusion. The multiscaling spectrum therefore supplies a quantitative fingerprint capable of distinguishing swarm migration from conventional tectonic earthquake sequences.","pith_inferences":["If the fingerprint is stable, automated pipelines could flag swarm-like migration from catalogue geometry alone without waiting for full moment-release statistics.","The same moment-scaling test might separate industrial induced-seismicity swarms from tectonic aftershocks in monitoring settings.","Repeating the analysis on three-dimensional hypocenters rather than epicenters would test whether Landers-style caging is only a fault-plane projection effect."],"forward_implications":["Multiscaling spectra can serve as a catalogue-based diagnostic to classify clustered seismicity as swarm-like versus mainshock–aftershock.","Swarm migration cannot be reduced to one diffusion exponent and must be treated as heterogeneous, multi-scale transport.","Negative-order moments mainly track compact pairs while positive-order moments track extremes, so the two branches separately constrain short-range and long-range migration.","Differences in high-order curvature among swarms may later be linked to distinct driving or geometric regimes once more sequences are measured."],"fun_headline_variants":["Swarms share multiscaling fingerprint tectonic quakes lack","Short and long swarm distances follow different scaling laws","Multiscaling spectrum flags anomalous diffusion in swarms","Earthquake swarms show nonlinear ε(q); Landers saturates","ε(q) fingerprint separates swarm migration from tectonic sequences"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The contrast that makes the spectrum a class fingerprint rests on treating the single Landers sequence as a fair stand-in for ordinary tectonic aftershock sequences in general.","fun_headline_variants_meta":{"raw":{"variants":["Swarms share multiscaling fingerprint tectonic quakes lack","Short and long swarm distances follow different scaling laws","Multiscaling spectrum flags anomalous diffusion in swarms","Earthquake swarms show nonlinear ε(q); Landers saturates","ε(q) fingerprint separates swarm migration from tectonic sequences"]},"model":"grok-4.5","effort":"low","cost_usd":0.003974,"raw_usage":{"total_tokens":1149,"prompt_tokens":730,"num_sources_used":0,"completion_tokens":64,"cost_in_usd_ticks":39740000,"prompt_tokens_details":{"text_tokens":730,"audio_tokens":0,"image_tokens":0,"cached_tokens":0},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":355,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":730,"tokens_out":64,"duration_ms":3930,"temperature":1.0,"reasoning_tokens":355,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T06:45:44.687155+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Compute ε(q) for several other large, well-relocated tectonic mainshock–aftershock sequences of comparable size; if those sequences also show a nonlinear positive-q branch without high-order saturation, the claimed fingerprint fails as a diagnostic boundary between swarms and conventional tectonic sequences.","supporting_citations":[],"review_version":1}