We tested whether annual nitrate plus nitrite load became more sensitive to flow after 2000 in the Mississippi and Atchafalaya rivers. We used USGS version 3.0 load estimates and complete annual flow records for 1980-2025. We fitted changes in log combined load against changes in log combined flow, with separate slopes for transitions ending in 1981-2000 and 2001-2025. Slopes were 1.021 and 1.037. The late minus early slope was 0.016, with conditional 95% three-transition block limits of -0.291 to +0.502. All nine planned checks also had limits that included zero. The data give no clear increase, but do not establish equal slopes. Load estimates use discharge as an input, so the relation is partly built into the source method. Annual averaging, two connected branches, short periods, and absent joint source-error draws limit physical inference. Code, small aggregate inputs, all fitted checks, bootstrap draws, a figure, source hashes, retrieval steps and an independent arithmetic audit are included.
Filters: arxiv: physics.geo-ph (Geophysics) ×
Question: USGS operational aftershock forecasts outside California start from generic Reasenberg-Jones parameters for 15 tectonic regimes (Page et al. 2016; OpenSHA OAF file dated 2016-03-21). Are these forecasts calibrated for later mainshocks, and for earlier ones? Approach: we used the USGS ComCat catalog (M 4.5 or more, 1985 to November 2026). We selected shallow (depth 70 km or less) M 6.0 or more mainshocks, and counted M 4.6 or more events within the OAF Wells-Coppersmith radius in days 1-7, 7-30 and 30-365. We copied the OAF forecast model (Poisson counts mixed over a Gaussian a-value). We wrote the plan before scoring. Results: for 779 mainshocks in 2016-2025 (2,337 window forecasts), 2.9% of counts were below the 2.5% forecast quantile and 2.8% above the 97.5% quantile; the PIT histogram did not differ clearly from uniform (p = 0.18); the best-fit shift of the mean a-value was +0.02 (95% CI -0.04 to +0.08). For 2,413 mainshocks in 1990-2015 the forecasts were too high: shift -0.16 (-0.18 to -0.12), observed counts 75% of the forecast mean. The difference, +0.18 (0.12 to 0.24), remained at M 5.0 or more, with a double radius, and without mainshocks that a larger event followed. In a post-hoc check, the shift by period followed the ComCat ratio of small to large background events (r = 0.95, slope 1.08, 10 periods); the yearly number of M 4.6 or more events rose by about 50% while M 5.5 or more did not change. A regional set fitted to 1990-2015 did not score better than the published set in 2016-2025. Limitations: our mainshock and counting rules are not those of Page et al.; calibration depended on the radius (too low with a double radius) and failed the PIT test at M 5.0; the catalog explanation is post hoc and does not separate detection from magnitude practice; we tested only the first generic forecast, not the updated forecasts that USGS issues.
We tested whether a fixed annual cycle improves later earthquake month probabilities. We used a shared USGS catalog with preferred magnitude at least 5.5. The main sample has depth 0-70 km and absolute latitude at least 15 degrees. We removed events within 100 km and 30 days after an equal or larger event. A fit to 2,243 events in 2001-2012 was scored on 2,316 events in 2013-2025. Southern months were shifted by six months. The model conditioned on year and hemisphere totals and used month lengths. The log score gain over a constant rate per day was -0.00377 bits per event. The conditional 95% three-year block interval was -0.00659 to -0.00112. The separate-year interval was -0.00741 to +0.00005. Wider event-removal windows and a magnitude 6 cut also gave negative point gains, with three-year intervals that included zero. This simple local-season fit gave no held-out gain. The result does not test local water loads or rule out regional seasonal effects. A revised catalog, a short test period, fixed fitted coefficients and approximate removal of related events limit inference. Code, source hashes, aggregate scores, bootstrap draws, a separate count audit and a figure are included. External event records are not bundled.
Schorlemmer et al. (2005) found that the Gutenberg-Richter b-value is highest for normal faults, lowest for thrust faults and between for strike-slip faults. We tested this relation out of sample in the Global CMT catalog. We fitted b as a linear function of the folded rake for shallow (centroid depth 50 km or less) Mw 5.7+ earthquakes in 1977-2003 (4,726 events) and tested the frozen model on 2004 to mid-2026 (4,838 events). The plan was written before the b-rake statistics. The training slope was -0.138 per 90 degrees of rake (95% tile-bootstrap CI -0.203 to -0.072). The test slope was -0.073 (-0.129 to -0.022), half as large. The frozen rake model gave no test log-likelihood gain over a constant b (+0.0002 nats per event, CI -0.0037 to +0.0029), so the pre-set transfer rule failed. With one b intercept per Bird (2003) plate-boundary setting, the test slope was -0.005 (-0.048 to +0.034); in training it was -0.098. Spreading ridges and oceanic transforms have high b and few thrust events; subduction zones have b near 1.0 and most thrust events. Inside subduction zones there was no rake slope in either period. In an exploratory comparison, a frozen setting model gave +0.0078 nats per event (CI +0.0032 to +0.0125) in the test period, and rake added nothing. The within-setting slope was negative in all three pre-2004 windows and near zero in all three later windows. This step coincides with GCMT procedure changes in 2004, which we did not test as the cause. Limitations: one catalog, Mw 5.7+ only, a coarse nearest-boundary setting label, strong rake-setting correlation, and a low-power test of the corner-magnitude effect. Code, plan, aggregate tables and figures are included; raw catalogs are retrieved by script.
We tested whether shallow and deep earthquakes occur close in space and time more often than expected under annual relative time shifts. We used a fixed global USGS ComCat snapshot for 2000-2025. Shallow events had depth 0-70 km. Deep events had depth greater than 300 km. Both had preferred magnitude M5.5 or more. The primary count was 26 cross-depth pairs within 500 km and one day in the same year. In 9,999 shifts of the full deep series within each year, the mean count was 24.32 and the central 95% range was 13-44. The two-sided Monte Carlo p-value was 0.7292. Ten shallow events followed deep events and 16 preceded them, counting pairs. Radius, time-window and split-period checks gave no clear excess. This is a conditional null result for large-event coincidence, not proof that deep and shallow systems are independent. The shift model changes relative seasonal phase and cuts clusters at year ends. Sparse counts, mixed magnitudes, location errors and catalog revision limit inference. Code, aggregate results, shift draws, an independent count audit and source hashes are included. External event records remain outside the package.
We tested whether recent monthly level variation improves six-month 90% intervals for Superior and Michigan-Huron. Both rules used the same local linear forecasts. Annual calibration used only earlier completed forecast errors. Each lake had 426 origins in 1990-2025. Splitting errors at the median of recent 12-month variation gave an equal-lake interval score reduction of -0.95%. The conditional 95% five-year block limits were -2.87% to +1.18%. Thus the primary test gave no clear gain. Total coverage stayed near 91%, while split-rule low-state coverage fell below 90% and high-state coverage rose above 93%. All sensitivity checks had negative point gains. Two lakes, revised data, dependent errors and one split rule limit inference. The test does not identify physical causes or evaluate operational forecasts. Code, numeric inputs, all forecast intervals, resampling draws, an audit, source hashes and a figure are included.
We tested whether four-season rainfall terms remove a post-drought runoff deficit in CAMELS-AUS v2.03. We selected 70 Victorian gauges by fixed data coverage rules. We fitted annual and seasonal rainfall models to log-transformed annual runoff before 1997. In 2010-2019, the equal-gauge mean residual was -1.10 baseline standard deviations for annual rainfall and -1.02 for seasonal rainfall. The change was +0.08, with a conditional 95% spatial and time block interval of -0.09 to +0.26. The remaining seasonal residual had an interval of -1.61 to -0.46. Seasonal terms raised earlier held-out prediction error by 3.00% in the primary model. Checks with SILO rainfall, a linear annual response and stricter coverage retained negative mean seasonal residuals. This result limits one statistical explanation. It does not exclude other effects of rainfall timing or show a physical cause. Short records, model form, shared weather, data error and fixed fitted parameters limit inference.
We tested the claim that the Gutenberg-Richter b-value depends on faulting style (Schorlemmer et al. 2005: normal highest, thrust lowest). We used shallow (centroid depth 50 km or less) Global CMT earthquakes of Mw 5.8 and larger. We fitted on 1977-2004 (n = 3,687) and tested on 2005 to June 2026 (n = 3,320). We used an exponential magnitude likelihood with ln b = alpha + gamma sin(rake), 5-degree cell fixed effects and a cluster bootstrap. Globally, gamma stayed negative but became half as large: -0.117 (95% CI -0.176 to -0.064), then -0.061 (-0.110 to -0.006). After 2005, normal b was not clearly higher than thrust b. Within 5-degree cells, gamma was -0.070 (-0.157 to +0.006), then +0.022 (-0.038 to +0.092). For 5.8 <= Mw < 6.5, the rake slope was zero or positive in both periods. For Mw 6.5 and larger (post hoc), gamma was -0.16 in both periods, globally and within cells (intervals exclude 0). Thrust b was 0.97 and 0.87; strike-slip b was 1.23 and 1.12. Class b-values fitted on 1977-2004 did not improve the 2005-2026 magnitude likelihood (+0.8 millinats per event, -1.8 to +3.4). Cell b-values did improve it (+7.0, +3.2 to +11.2), and the faulting class added nothing beyond the cell b-values. Thus the faulting-style effect is in the large-event tail and is not a general change of slope. This agrees with rupture-width saturation more than with a stress control that acts at all magnitudes, but we did not test the mechanism. Limitations: the GCMT method changed in 2004, the Mw 6.5 split is post hoc, and the within-cell intervals are wide.
We tested whether large deep earthquakes have greater April-September odds than smaller deep earthquakes after matching year and broad location. A fixed USGS catalog for 2016-2025 contains 202 events deeper than 500 km with preferred magnitude at least 5.5. Nine have magnitude at least 7.0, and five of these occurred in April-September. The matched null expects 5.10. An exact conditional test gives an upper-tail probability of 0.6851. The common odds ratio is 0.94, with an exact 95% interval of 0.15 to 6.19. Six descriptive checks also give intervals that include one. This extends a prior campaign test of the same nine large events by using smaller-event times as a background control. It is not independent event evidence. The result gives no clear seasonal size contrast and permits a substantial effect. Small counts, exchangeable size labels, broad location bins, mixed magnitude scales and catalog revisions limit inference. Code, aggregate results, source hashes, an independent exact arithmetic audit and a figure are included. External event data are not bundled.
We tested a 24-hour local mean solar cycle in a fixed USGS ComCat sample of 10,473 earthquakes with preferred magnitude at least 5.5 and depth at most 70 km in 2000-2025. The fitted first harmonic amplitude was 3.85%. An independent-event test gave p = 0.0208. A test that rotated each complete month as one group gave p = 0.19575. Its null root mean square amplitude was 53.3% larger than the independent-event value. Weekly and yearly rotations also gave p above 0.1. The noon harmonic coefficient was +0.17%, with a conditional 95% year resampling interval of -2.52% to +2.80%. These tests do not resolve a common daily cycle. They do not exclude local effects or identify a physical cause. Group phase assumptions, mixed magnitude types, catalog changes, and only 26 years limit inference. The package includes code, source hashes, aggregate results, simulation draws, a figure and an independent arithmetic audit. Raw external records remain outside the package.
Many studies report that the Gutenberg-Richter b-value is highest for normal faults and lowest for thrust faults, and use b as a stress meter. We tested this in 10,194 shallow (depth 40 km or less) Global CMT earthquakes of Mw 5.6 or more from 1977 to mid-2026, with a plan written before the analysis. For all years, b was 1.14 for normal, 1.06 for strike-slip and 0.98 for thrust events (thrust/normal ratio 0.86, 95% cluster-bootstrap CI 0.80-0.92). The full order was present in 1977-2004, but not in 2005-2017, and in the new 2018-2026 data the ratio (0.88, 0.74-1.02) did not exclude 1. When we compared events only inside the same Flinn-Engdahl region (exact exponential likelihood with region fixed effects), the ratio was 0.98 (0.91-1.06): only 13% of the raw log contrast remained. Mc 5.8 or 6.0, depth to 70 km and plunge-based classes gave the same result; 10-degree cells and a declustered catalog kept a small part (ratios 0.93 and 0.91). The difference is not a corner-magnitude effect: it stayed in Mw 5.6-6.5 (ratio 0.87, 0.78-0.96) and in tapered fits with class-specific corners. An exploratory plate-boundary split (PB2002) showed the cause: normal events at oceanic spreading ridges had b = 2.06 and at oceanic transforms 1.50, while subduction normal and thrust events both had b near 1. Inside subduction-zone regions, thrust b was 15% higher than normal b, the opposite of the stress-meter order. Limitations: only Mw 5.6+ global data, so small-earthquake regional results are not tested; Flinn-Engdahl regions are large and mix stress environments; completeness and the 2004 GCMT method change can differ by region and period; the setting split is post hoc.
We tested a fixed spatial window in a global ComCat snapshot. Centers had preferred magnitude M7 or more, depth below 70 km, and origin years 2001-2024. We removed centers with an equal or larger event within 600 km and 365 days on either time side. This left 167 centers. We counted M5.5+ neighbors at gaps of 1 to below 30 days and subtracted before counts from after counts. For 20 M8+ centers, the net count was 198 within 200 km and 283 within 600 km. The signed share outside 200 km was 30.0%, with a conditional 95% three-year block interval of 16.1% to 40.0%. It was 21.5% after the two M9+ centers were removed. The M7 to below 7.5 group had a -9.1% share. The large-minus-small group contrast was 39.1 percentage points, with an interval of 12.9 to 74.2 points. The late-days contrast was uncertain. This is a window-size benchmark, not a new physical law. The shares are signed count differences, not identified causal aftershocks or all events beyond 600 km. Regional differences, catalog errors, mixed magnitude types and overlapping windows limit inference. Code, aggregate counts, source hashes, a scalar audit and a figure are included. External event records are excluded.
We tested whether late M5.5+ event excess is larger near the depth of a deep earthquake than far from it. We used a fixed global ComCat snapshot and 82 centers with depth greater than 300 km and M6.5 to below 8 in 2001-2024. Centers had no equal or larger event within 200 km and 365 days on either side. At horizontal distance up to 200 km, days 7-30 gave 2 before and 12 after counts within 50 km of the center depth. Depth differences above 150 km gave 0 before and 7 after counts. The near-minus-far excess was 0.0366 events per center, with a conditional 95% three-year block interval of -0.1081 to 0.1482. All seven parameter settings gave contrast intervals that included zero. The data do not resolve a larger near-depth excess. They do not show equal band rates or a shallow crustal response. Sparse counts, regional dependence, selection, location errors and mixed magnitude types limit inference. Code, aggregate counts, bootstrap weights, independent checks and exact source retrieval instructions are included. External catalog records are excluded.
We tested whether a known rise in earthquake spatial spread was limited to the most distant pairs. We used the GEOAZUR relocated catalog for the 2021 Haiti earthquake. For magnitude at least 3 within 100 km, median epicentral pair distance rose from 28.03 to 40.30 km between days [7,30) and [30,90). The late-to-early ratio was 1.438. A conditional 95% interval from three-day block weights was 1.171-1.938. The seven-day interval was 1.156-1.814. Thus the spatial change reached the median pair. The 90th-percentile-to-median ratio had an unresolved change. Within 50 km, the median distance ratio was 1.065 and its interval included one. These results permit changes in cluster weights and do not establish local fault growth, fluid motion, triggering, or stable detection. Pair distances are dependent. The weight intervals are sensitivity ranges for this catalog, not validated process confidence bounds. The package gives aggregate evidence, code, source hashes, a figure, and an independent arithmetic check. External raw records remain outside the package.
Question: after a long earthquake swarm, is a large earthquake more probable than the number of swarm events alone predicts? We used the JMA unified hypocenter catalog for 1998-2024 (depth 25 km or less, M 2.0 or more, far-offshore regions removed). An objective rule found sequences with high monthly rates in 0.2-degree windows that lasted 6 months or more. We kept swarm-like sequences, where the largest event did not start the sequence. The outcome was the first M 6.0 or larger event within 10 km of the sequence area, during the sequence or within one year after it. Of 26 swarm-like sequences, 2 had this outcome (7.7%, 95% interval 0.9-25%): Niijima 1998-1999 (followed by the M 6.5 Izu Islands event of July 2000) and Noto 2021-2023 (followed by the M 6.5 event of May 2023 and the M 7.6 event of January 2024). This count was 4.3 times the count expected from the mean Japanese rate of independent crustal M 6+ events for the same area and time (expected 0.47, p = 0.08). But it was 0.44 times the count expected if each M 2+ swarm event had the catalog-wide chance of being M 6+ (expected 4.5, 95% interval of the ratio 0.05-1.6). It was also 0.8 times the count expected from the long-term activity of the same areas. In 18 design variants, the ratio to the event-count expectation was 0.26-0.89 and never significantly above 1. The ratio to the area expectation was 2.4-8.1 (p < 0.05 in 8 variants). Thus long swarms raised the local chance of an M 6+ event, but by an amount that their event counts explain. We found no extra swarm-specific precursor effect. Exploratory comparisons of duration, event count, largest magnitude, b-value, rate trend and depth did not clearly separate the two cases from the 24 other sequences. Limitations: only 2 positive cases, so the upper bound is wide; one case was a volcanic dike intrusion; the detection rule and the trigger flag are simple; the JMA catalog is not homogeneous in time.
We tested whether broad geographic mix explains a larger large-event fraction at greater depth in a global ComCat snapshot. In nine common 20-degree cells in 2000-2025, the fixed-weight M6.5+ share among preferred M5.5+ events was 17.36% at depths of 300 km or more and 9.02% at 70 to below 300 km. The difference was 8.34 percentage points, with a conditional 95% three-year block interval of 3.54 to 13.52 points. Common support retained 1,140 intermediate and 687 deep events. Moment-family magnitudes and two other cell widths gave positive contrasts. The M7+ share among M6+ events gave a less certain result: 5.15 points, with an interval of -0.26 to 10.80. Thus the tested geographic adjustment did not remove the main contrast. This extends a prior campaign size-forecast test with direct standardized shares; the data overlap. Mixed magnitude methods, within-cell geography, depth errors, event sequences, data-selected support and conditional intervals prevent a causal depth or hazard claim. Code, aggregate evidence, uncertainty draws, an arithmetic audit, a figure and exact external-source provenance are included. External event records remain outside the package.
We tested whether relatively large nearby events were more common after selected shallow earthquakes than before them. A fixed global ComCat snapshot gave 171 centers of preferred magnitude 7 to below 8 in 2001-2024 after equal isolation on both time sides. Each center counted at most one success per side. A success had at least one event within one magnitude unit of the center, within a 300 km hypocentral chord distance, at least one day and less than 30 days away. There were 30 after successes and 14 before successes. The paired fraction difference was 9.36 percentage points, with a conditional 95% three-year block interval of 2.78 to 16.38 points. Radius, isolation and spatial resampling checks retained positive primary intervals. A half-unit gap and days 7-30 alone gave intervals that touched or crossed zero. Most centers had neither event. The primary excess thus remains when rich sequences cannot contribute many successes per center. It does not identify causal aftershocks, test the mean Bath magnitude gap or give a hazard forecast. Selection, sparse discordant centers, catalog revisions and dependence limit inference. The package gives code, aggregate evidence, source hashes, a figure, independent count checks and replay instructions.
We tested how group boundaries affect a seasonal earthquake test in a fixed USGS ComCat sample. We learned an annual phase from 2,568 shallow extratropical M5.5+ events in 2000-2012 and tested it on 2,558 events in 2013-2025. Connected event groups with pair links of at most 200 km and 30 days gave an exact rotation-null standard deviation of 7.65 percentage points. Fixed 20-degree cell-year groups gave 6.22 points. This is a 22.9% increase in standard deviation and a 51.1% increase in variance. Splitting the same event groups at cell-year boundaries reduced their variance by 24.5%. All nine distance-time checks gave larger spread than cell-years. The seasonal coefficient stayed at +1.25%, with two-sided conditional p = 0.874 in the primary event-group model. The earlier finding of no clear phase transfer remains. The result measures sensitivity to grouping assumptions. It does not show true uncertainty coverage or a physical load mechanism. Exact variance identities, aggregate results, source hashes, retrieval and analysis code, an independent audit and a figure are included. Raw and record-level external data remain outside the bundle.
Goswami et al. (2006) reported more heavy monsoon rain events over central India from 1951 to 2000. Lin and Huybers (2019) said that a station-network change after 1975 can cause part of this rise. We did an out-of-sample test with IMD 0.25-degree daily rainfall for 2001-2025. We used the Goswami definitions (June-September, 74.5-86.5E and 16.5-26.5N, 100 and 150 mm/day). We fitted a stationary mean, a log-linear trend and a step at 1976 to 1951-2000 only. For 1-degree box means, the 2001-2025 heavy-event mean was 0.97 of the trend forecast (95% season-bootstrap interval 0.80-1.17) and 1.39 of the step forecast (1.15-1.67). Without the three wet seasons 2005-2007, the ratios were 0.82 (0.72-0.92) and 1.18 (1.04-1.33). For 0.25-degree cells, the counts were near the step forecast and below the trend forecast. Very heavy cell events were 0.80 of the trend forecast (0.65-0.97), and 0.67 without 2005-2007. Inside 2001-2025, heavy counts had no clear trend (-8% per decade, -32% to +23%). Moderate events stopped their decline, and seasonal rain did not change. Box-mean extremes rose faster than cell extremes. A station-network proxy (local extrema in the interpolated fields) shows the 1971 station loss and a thinner network in 2005-2015. It explains about +3% to +12% of the +39% rise in 1-degree heavy events from 1976-2000 to 2001-2025. PERSIANN-CDR agrees with IMD seasonal totals (r = 0.92), but its daily extremes fall after 1997, so it cannot check the trend. Limitations: one gridded product, no public station counts, an indirect network proxy, and results that depend on a few seasons.
Question: Do aftershock sequences decay faster (larger Omori-Utsu p) where surface heat flow is higher? Kisslinger and Jones (1991) reported this for southern California. An 18-sequence Japanese study found no link. Approach: We selected 450 isolated shallow (depth 30 km or less) JMA mainshocks of magnitude 4.5 or more from 1998-2023. We fitted the Omori-Utsu rate with a fixed background by maximum likelihood over 0.02-365 days. We related p to the median IHFC GHFDB R2024 heat flow within 50 km. Intervals come from a one-degree-cell block bootstrap. Result: 191 sequences passed quality control and had heat-flow data. p did not change clearly with heat flow. Spearman rho = 0.005 (95% interval -0.16 to 0.19). The random-effects slope was 0.10 p per decade of heat flow (-0.11 to 0.34), and 0.16 (-0.05 to 0.40) with magnitude, depth, onshore and Mc controls. Median p was 1.11, 1.09 and 1.11 in the low, middle and high heat-flow terciles. All 15 Spearman intervals (primary and 14 pre-specified variants) include zero. Only the 30 km radius variant gave slope intervals above zero. log10 c, normalized productivity and distance to Holocene volcanoes also gave no clear link. In synthetic catalogs, the test detected an injected slope of 0.25 p per decade in 19 of 20 runs. Limitations: Surface heat flow is a coarse proxy for temperature at depth, and the data cover Japan unevenly. The windows are simple, without ETAS declustering. We amended the selection rules once, before any heat-flow statistic. A weak slope (up to about 0.4 p per decade) is not excluded.