11,505 Publications · 65,152 Citations

Stefani et al. (2025-2026) state that ground-level enhancements (GLEs) and large X-class flares occur in the positive phase of a phase-stable 629.29-day beat of planetary tides. They give p-values of 3e-5 to 7e-6, and they report peaks at nine predicted tidal periods. Those p-values treat each event as independent, but these events occur in bursts. Also, the period was first found in the same GLE data. We used 72 GLEs (1956-2025), 37 flares of class X10 or more, and NOAA GOES flare lists for 1975-2026. A 0.69 flux-scale ratio between GOES-15 and GOES-16 shows that the X1 threshold changed in 2017; we corrected this. We then repeated the tests with nulls that give a random phase to each 27-day or 81-day burst or to each solar cycle. We reproduce the published phases (2.04, 2.96 and 2.58 rad against 2.01, 2.96 and 2.58). With clustered nulls, the fixed-period p-values become 0.024-0.051 for GLEs and 0.07-0.31 for X1+ flares. After a search of the 450-900 day band, no set has p < 0.05. The 629-day component of the parameter-free planetary function has a phase within 3 days of the GLE phase, a chance coincidence of about 1 percent. The GLE V-test against that phase gives p = 0.004 for 27-day bursts. In solar cycle 25, which was not used to choose the period, X1+ flares lean to the planetary phase (mean cosine 0.35-0.38, p about 0.05). All 20 out-of-sample combinations are positive, but with clustered nulls p is 0.03-0.34. Five cycle-25 GLEs give p = 0.22. The seven new predicted periods do no better than random periods (p = 0.08-0.24). Limitations: few independent bursts per cycle, a time-gap proxy for active regions, and a post-hoc phase-coincidence estimate. Cycle 26 can give a decisive test.

Question: Marchitelli et al. (2020) reported that global M5.6+ earthquake counts increase on the first day after the daily SOHO proton density falls below about 15.5 per cm3 (1996-2016, p < 1e-5). Does this effect survive a different spacecraft, removal of single aftershock days, and new data from 2017-2026? Approach: We used the paper's relative-rate statistic R, its interval-shuffle null and a circular-shift null that keeps clustering. Data: SOHO CELIAS/MTOF PM hourly density, OMNI2 (Wind/ACE) density, USGS ComCat M5.6+ (11,851 events) and ISC-GEM magnitudes from ComCat (7,459 events, 1996-2016). Results: The in-sample effect replicates (SOHO, T = 15.5: R = 1.96 ComCat and 1.97 ISC-GEM, p < 1e-5 under the paper's null). One day causes it: 2011-03-11 (Tohoku), with 105 ComCat events, is a condition day because the SOHO daily mean on 2011-03-10 was 17.98. Without that day, R = 1.03 (p = 0.40). For days with any event, R = 1.00. OMNI for the same 21 years gives R = 0.89. The paper's null is too narrow (SD of R 0.108 against 0.171 for the circular null). In the 2017-2026 out-of-sample data, OMNI at 15.5 gives R = 1.13 (95% CI 0.85-1.52, p = 0.19). SOHO with a quantile-matched threshold gives R = 1.25 (0.84-1.77, p = 0.10). SOHO with the paper's min-max scaling gives R = 0.94 (p = 0.77). Declustered OMNI gives R = 1.01. The SOHO density scale fell relative to OMNI after 2015 (median ratio 1.00-1.23 before, 0.76-0.91 after). Thus SOHO reached 15.5 on only 4 days in 2017-2026. Limitations: We did not use ISC-GEM v7.0 itself. We tested only density and the first-day-below condition. The out-of-sample data reject R near 2 but not R of about 1.2-1.4, which is the R size in the 1996-2005 SOHO subset.