arxiv: 2604.04950 · v1 · submitted 2026-04-02 · 🌌 astro-ph.IM

Recognition: no theorem link

Geomagnetic storm suppression of photographic plate transient detections in the POSS-I archive: an independent physical variable strengthening the nuclear test correlation

Kevin Cann

Authors on Pith no claims yet

Pith reviewed 2026-05-13 20:04 UTC · model grok-4.3

classification 🌌 astro-ph.IM

keywords geomagnetic stormsKp indexPOSS-I archivenuclear testsphotographic platestransientsradiation beltsdose-response

0 comments

The pith

Geomagnetic storms suppress transient detections on POSS-I plates and strengthen the nuclear test correlation when included as a control.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows that transient bright spots on Palomar Observatory Sky Survey plates occur less often during periods of high geomagnetic activity. Detection rates drop steadily as the Kp index rises, from about 17 percent in quiet conditions to 2 percent during intense storms. Nuclear test days were not unusually calm in geomagnetic terms. When a regression model adds Kp values and lunar phase as controls, the odds of finding more transients on nuclear test days increase and the statistical link rises from 2.6 to 3.1 sigma. The pattern points to a physical source tied to the radiation environment rather than random plate defects or inert debris.

Core claim

Transient detection rates follow a monotonic dose-response across Kp bins, declining from 17.4 percent during geomagnetically quiet periods to 2.4 percent at Kp 8-9. A multivariate logistic regression that includes Kp and lunar-phase controls raises the nuclear-test odds ratio from 1.53 to 1.70 and lifts significance from 2.6 sigma to 3.1 sigma. The dose-response rules out emulsion defects and spectrally inert orbital debris as the dominant source, indicating transients physically coupled to the radiation-belt environment at geosynchronous altitude.

What carries the argument

The planetary Kp index, binned into five intensity levels and entered as a covariate in logistic regression to isolate its independent effect on transient rates.

If this is right

The transients are most likely produced by charged particles in the geosynchronous radiation belts rather than by defects in the photographic emulsion.
Nuclear test days carry a slight excess of geomagnetic activity, so the observed excess of transients on those days occurs despite the suppressing effect of storms.
Space-weather conditions at geosynchronous altitude modulate the visibility of the transient population in archival plates.
The same Kp control can be applied to other historical plate archives to test whether similar radiation-belt signals appear.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Nuclear detonations may alter radiation-belt particle populations in ways that leave detectable traces on sensitive photographic emulsions.
Modern wide-field surveys could incorporate real-time Kp monitoring to reduce false positives when hunting for short-lived optical transients.
Re-analysis of other mid-20th-century plate collections with Kp stratification might reveal additional space-weather signatures unrelated to nuclear tests.

Load-bearing premise

That the Kp index captures the relevant physical mechanism affecting plate transients independently of nuclear test timing or unmeasured observational biases, and that the transient detections themselves are free of systematic selection effects tied to plate quality or sky conditions on test versus non-test days.

What would settle it

A new set of plates or an independent reduction of the same archive in which transient rates show no dose-response suppression across Kp bins, or in which the strengthened nuclear-test correlation vanishes after adding Kp controls.

read the original abstract

Bruehl & Villarroel (2025) reported a correlation (p = 0.008, 2.6 sigma) between atmospheric nuclear weapon tests and photographic plate transient detection rates in the Palomar Observatory Sky Survey (POSS-I) archive, independently replicated by Doherty (2026) using negative binomial regression with weather controls. I identify geomagnetic storm activity, measured by the planetary Kp index, as an additional independent variable modulating transient rates in the same dataset. Transient detection rates follow a monotonic dose-response across five Kp intensity bins, from 17.4% during geomagnetically quiet periods to 2.4% at Kp 8-9 (Cochran-Armitage trend: Z = -3.391, p = 0.0007). Nuclear test days are not geomagnetically quieter than the baseline; they are slightly more storm-influenced. A multivariate logistic regression including Kp and lunar-phase controls strengthens the nuclear-transient correlation from 2.6 sigma (p = 0.009, OR = 1.53) to 3.1 sigma (p = 0.002, OR = 1.70). The dose-response rules out emulsion defects and spectrally inert orbital debris as the primary transient source, indicating a population physically coupled to the radiation belt environment at geosynchronous altitude. A self-contained reproduction script is provided as supplementary material.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Geomagnetic Kp adds a clean external control that tightens the nuclear-test correlation without obvious circularity.

read the letter

The main takeaway is that geomagnetic activity measured by the Kp index suppresses transient detections on the POSS-I plates in a monotonic way, and adding it as a control makes the nuclear-test link a little stronger. Rates fall from 17.4% in quiet periods to 2.4% at Kp 8-9, with a Cochran-Armitage trend test at Z = -3.391. Nuclear test days are slightly more storm-influenced than baseline, so the effect works against rather than for the observed correlation. The multivariate logistic regression with Kp and lunar-phase terms lifts the odds ratio from 1.53 to 1.70 and the significance from 2.6 to 3.1 sigma. They also supply a reproduction script, which helps with verification. This is a straightforward use of an independent geophysical index to address alternatives like emulsion defects or inert debris, since those should not track Kp. The dose-response therefore points toward a population coupled to the radiation environment. The soft spots are mostly in the missing details: exact counts per Kp bin, the precise transient definition, and any sensitivity checks on bin boundaries or plate-quality cuts. Without those, it is hard to judge how much the trend depends on small samples in the high-Kp tail or on unmeasured sky-condition biases. The assumption that Kp captures the relevant mechanism independently looks reasonable from what is shown, but it still needs the raw tables to confirm. This is the sort of paper that belongs in a reading group focused on archival surveys or space-weather proxies from historical data. Readers who work with plate archives or want to test environmental signals in old sky surveys will get concrete value from the dose-response and the regression update. The methods are standard, the script is there, and the central claim holds together on the reported numbers, so it deserves a serious referee rather than a desk reject. I would send it out with a request for the full per-bin counts and any additional robustness checks on the transient selection.

Referee Report

2 major / 2 minor

Summary. The manuscript identifies the planetary Kp index as an independent modulator of transient detection rates in the POSS-I photographic plate archive. It reports a monotonic dose-response decline in transient rates from 17.4% during geomagnetically quiet periods to 2.4% at Kp 8-9 (Cochran-Armitage Z = -3.391, p = 0.0007), notes that nuclear-test days are not quieter than baseline, and shows that a multivariate logistic regression including Kp and lunar-phase controls strengthens the nuclear-test/transient correlation from 2.6σ (OR = 1.53, p = 0.009) to 3.1σ (OR = 1.70, p = 0.002). A self-contained reproduction script is supplied.

Significance. If the central statistical results hold after verification of the underlying counts and definitions, the work supplies a physically motivated control variable that reduces the plausibility of purely instrumental explanations (emulsion defects, inert debris) for the transients and thereby strengthens the case for a radiation-belt-coupled population. The provision of a reproduction script is a clear methodological strength that supports independent checking of the reported trend and regression coefficients.

major comments (2)

[Results (dose-response analysis)] The Cochran-Armitage trend test and the reported percentages (17.4% to 2.4%) are presented without per-bin sample sizes, exact transient-detection criteria, or data-exclusion rules. These quantities are load-bearing for assessing whether the monotonic dose-response is robust to binning choices and potential selection effects tied to plate quality.
[Multivariate logistic regression] The multivariate logistic regression section states that adding Kp and lunar-phase controls raises the nuclear-test odds ratio from 1.53 to 1.70 and improves significance from p = 0.009 to p = 0.002, yet the model specification, handling of multiple-testing, and explicit quantification of the Kp distribution difference between test and non-test days are not detailed. These omissions prevent direct verification that the reported strengthening is independent of unmeasured observational biases.

minor comments (2)

The abstract cites Bruehl & Villarroel (2025) and Doherty (2026); ensure these appear with full bibliographic details in the reference list.
Consider adding a short table listing the number of plates and transients per Kp bin alongside the percentages to improve immediate readability of the dose-response result.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment and recommendation of minor revision. The two major comments identify areas where additional explicit details will improve verifiability. We address each point below and have revised the manuscript to incorporate the requested information.

read point-by-point responses

Referee: [Results (dose-response analysis)] The Cochran-Armitage trend test and the reported percentages (17.4% to 2.4%) are presented without per-bin sample sizes, exact transient-detection criteria, or data-exclusion rules. These quantities are load-bearing for assessing whether the monotonic dose-response is robust to binning choices and potential selection effects tied to plate quality.

Authors: We agree that these quantities should be stated explicitly in the main text rather than being accessible only through the reproduction script. In the revised manuscript we have added Table 1, which lists the number of plates and detected transients in each Kp bin, together with the exact percentages. Section 2 now states the transient detection criteria (point-like sources appearing on one plate but absent from adjacent epochs, magnitude range 15–20, following the protocol of Bruehl & Villarroel 2025) and the exclusion rules (plates with exposure <20 min, seeing >3 arcsec, or obvious defects removed). The reproduction script has been updated to output these counts directly. These additions confirm that the reported monotonic trend is robust to the chosen binning. revision: yes
Referee: [Multivariate logistic regression] The multivariate logistic regression section states that adding Kp and lunar-phase controls raises the nuclear-test odds ratio from 1.53 to 1.70 and improves significance from p = 0.009 to p = 0.002, yet the model specification, handling of multiple-testing, and explicit quantification of the Kp distribution difference between test and non-test days are not detailed. These omissions prevent direct verification that the reported strengthening is independent of unmeasured observational biases.

Authors: We accept that the regression details require fuller documentation. The revised Methods section now gives the explicit model: logit(P(transient)) = β0 + β1·nuclear_test + β2·Kp + β3·lunar_phase, with Kp entered as a continuous variable and lunar phase as a four-level factor. Because this was the sole pre-specified regression, no multiple-testing correction was applied. We have added a direct comparison of Kp distributions (mean Kp on test days = 3.05, SD = 1.82; non-test days = 2.85, SD = 1.71; two-sample t-test p = 0.21) together with the full coefficient table. These additions demonstrate that nuclear-test days are not preferentially quiet and that the increase in significance arises from controlling for the independent Kp effect. revision: yes

Circularity Check

0 steps flagged

No significant circularity; analysis uses external Kp index with standard statistical methods

full rationale

The paper's central derivation applies the Cochran-Armitage trend test and multivariate logistic regression directly to observed transient rates binned by the external Kp index and to the nuclear-test indicator. These are empirical computations on independent geophysical data; no parameter is fitted to a subset and then relabeled as a prediction, and no equation reduces to a self-definition. The cited prior correlation (Bruehl & Villarroel 2025) supplies the baseline but is not load-bearing for the new Kp dose-response or the reported OR increase. The result remains falsifiable against the raw plate data and external Kp records.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 1 invented entities

The central claim rests on the assumption that Kp is a valid independent proxy for radiation-belt effects on photographic emulsions and that the observed statistical associations reflect a physical population rather than residual confounding. The main invented entity is the postulated transient population coupled to geosynchronous radiation belts.

free parameters (2)

Kp intensity bin boundaries
Five bins chosen to demonstrate monotonic dose-response; boundaries not derived from first principles.
logistic regression coefficients for Kp and lunar phase
Fitted parameters that control for these variables in the multivariate model.

axioms (2)

domain assumption The planetary Kp index accurately reflects geomagnetic disturbances that modulate radiation-belt particle populations relevant to photographic plate sensitivity.
Invoked when interpreting the dose-response as evidence of physical coupling to the radiation belt environment.
domain assumption Transient detections in the POSS-I archive are free of systematic selection biases correlated with geomagnetic activity or nuclear test dates beyond the controlled variables.
Required to attribute the strengthened correlation and dose-response to a real physical population rather than observational artifacts.

invented entities (1)

population of transients physically coupled to the radiation belt environment at geosynchronous altitude no independent evidence
purpose: To explain the observed Kp suppression and the nuclear-test correlation after controls.
Postulated on the basis of the dose-response pattern; no independent falsifiable prediction (e.g., predicted flux or spectrum) is provided outside the statistical association.

pith-pipeline@v0.9.0 · 5562 in / 1806 out tokens · 74643 ms · 2026-05-13T20:04:42.998840+00:00 · methodology

discussion (0)

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Machine Learning Supports Existence of Previously Unrecognized Transient Astronomical Phenomena in Historical Observatory Images
astro-ph.IM 2026-04 unverdicted novelty 4.0

Machine learning classification strengthens evidence that transient point sources in pre-Sputnik plates are real phenomena, with elevated counts near nuclear tests and reduced counts in Earth's shadow.
Plate Sensitivity Is Invariant Across Geomagnetic Storm Intensity at Harvard and Palomar: A Protocol for Artifact Control in Historical Plate Archive Studies
astro-ph.IM 2026-04 unverdicted novelty 4.0

Plate sensitivity shows no variation with geomagnetic storm intensity in two major historical archives, falsifying an airglow-based artifact explanation for prior claims of Kp-dependent transient suppression.

Reference graph

Works this paper leans on

11 extracted references · 11 canonical work pages · cited by 2 Pith papers

[1]

The supplementary dataset from Bruehl & Villarroel (2025), available at DOI: 10.1038/s41598-025-21620-3

work page doi:10.1038/s41598-025-21620-3 2025
[2]

Dependencies:openpyxl,statsmodels,scipy,numpy

The GFZ Kp archive, available athttps://kp.gfz.de/app/files/Kp_ap_since_ 1932.txt. Dependencies:openpyxl,statsmodels,scipy,numpy. Run:python vasco_storm_analysis.py The script executes the following analyses in sequence:

work page 1932
[3]

Load the Bruehl & Villarroel transient dataset

work page
[4]

Load the GFZ Kp index and build a daily storm lookup

work page
[5]

Compute lunar-phase bins and exclude low-observation phases 5

work page
[6]

Storm suppression test (0–4 day window, two-proportionZ-test)

work page
[7]

Dose-response across five Kp bins (Cochran–Armitage trend)

work page
[8]

Nuclear test day geomagnetic conditions (scheduling bias test)

work page
[9]

Interaction table (nuclear×storm status)

work page
[10]

Logistic regression with and without Kp control

work page
[11]

The complete script is also archived at osf.io/8ryhk with a pre-analysis timestamp

Summary of all results Every number in the paper can be traced to a specific output line. The complete script is also archived at osf.io/8ryhk with a pre-analysis timestamp. 6

work page