Global Time Echoes: Empirical Synthesis
Abstract
Analysis of 25.3 years of GNSS timing data (2000–2025) reveals a persistent, distance-structured correlation in global atomic clock networks that tests whether corrected distributed GNSS timing residuals contain a coherent spatial covariance beyond the standard relativistic and geodetic measurement model. Examination of 165.2 million station-pair observations from 474 unique receivers (814 station codes) over a 25.3-year baseline demonstrates a spatial correlation signal decaying with distance over a finite scale (Gaussian-family kernels preferred in the 25-year dataset, exponential ΔAIC = 12.8; power law excluded at ΔAIC > 30; azimuth-averaged exponential scale λ = 3,210 ± 389 km, with cross-sector directional dispersion λ T = 4,201 ± 1,967 km across eight azimuthal sectors, a measurement of order thousands of kilometres rather than a single fixed scale), following cross-centre validation yielding R² = 0.92–0.97 across three distinct analysis-centre software architectures on a 62.7 million observation baseline. These findings emerge from a systematic multi-paper research program: theoretical framework development with pre-specified expectations and theoretical search ranges (Paper 0), multi-center consistency across software-diverse product families within the shared IGS tracking ecosystem (Paper 1), 25-year longitudinal analysis enabling long-period geophysical detection (Paper 2), raw-data consistency test strongly constraining precise-product processing artifacts (Paper 3), and cosmological extension testing whether temporal covariance or gradient effects contribute to dark-sector phenomenology (Paper 4). The evidence programme combines timing-network signatures, environmental tests and an optical sensitivity bound: distance-structured spatial decay over a finite scale; East-West/North-South anisotropy (ratio 2.16 ± 0.23, aggregate across the 25.3-year dataset; nominal pair-count p < 10 −15 ); orbital velocity coupling (r = −0.888, 3.65σ autocorrelation-corrected; surrogate p < 2 × 10⁻⁷ as a resolution bound); a recovered annual-phase direction near Earth's aphelion velocity tangent (RA = 6°, Dec = +4°; directional-rank permutation p < 10⁻⁵; the directed CMB-apex template is anti-phase and the Solar Apex is geometrically incompatible, predicting N-S anisotropy against the observed E-W dominance; no global full-sky-search surrogate p-value has been computed; on the MGEX product the scan is projection-degenerate); an optical-domain sensitivity bound (SLR sampling-matched controls reproduce the apparent spectral concentration; inter-station candidates remain exploratory); a planetary event catalogue (56/156 significant at ≥ 2σ under the free-centre pulse statistic — repeatable estimator behaviour whose excess is attributable to calendar-position sampling in the season-preserving RINEX controls; the 40-event opposition subset is null-validated at 5.5 × and a continuous planetary-potential correlation survives at r = +0.116, p = 0.012); 18.6-year lunar nutation coupling (R² = 0.641 ± 0.10 over 1.36 cycles); pulsar spin-down excess (exploratory channel); and host-potential H 0 correlation. On the held-out MGEX product the distance-structured coherence replicates (λ ≈ 1,862 km), while the anisotropy and primary orbital-velocity metrics do not replicate robustly. Raw RINEX consistency test using Single Point Positioning with broadcast ephemerides achieves consistent signal detection across all 72 metric combinations (t-statistics up to 112, Cohen's d up to 0.304), strongly constraining precise-product processing artifacts as the sole origin. Broadcast ephemerides still contain control-segment information; Satellite Laser Ranging supplies a GNSS-product-independent test pathway, whose present sensitivity and exploratory correlations do not yet adjudicate the origin of the GNSS signal, while the generic orbit-projection confound is reserved for the physical-confound forward model and orbit-free clock channels. The network's selectivity profile—sensitive to velocity-dependent dynamics while blind to GM/r² scaling and solar rotation—characterizes it as an inertial interferometer measuring correlation geometry rather than a gravimeter measuring Newtonian force. Table 5 scores these observations against the pre-specified expectations of the Temporal Equivalence Principle: the orbital-velocity sign, EW > NS ordering, band localization, and GM/r² null are met, while the 2,000–3,000 km screening-parameterized forecast, the < 5% cross-centre criterion, and—in the 25-year data—the exponential-kernel preference are not met. TEP is a bi-metric scalar-tensor framework in which proper time is a dynamical field governed by a conformal factor A(ϕ) = exp(β A ϕ/M Pl ). Temporal Shear suppression operates through the continuous spatial profile of the ϕ field (Temporal Topology), with suppression arising from the non-linear superposition of field gradients (Temporal Shear), replacing discrete thin-shell approximations with a geometrically continuous mechanism. The observed correlation length λ T represents the characteristic scale of Temporal Topology. Environmental screening, derivative-screened, and saturation-radius interpretations are candidate theoretical completions, not the primary empirical claim; the saturation-radius interpretation R T is adopted in the present framework (§4.4.2). The framework preserves local Lorentz invariance while predicting spatially correlated clock-network residuals in the conformal sector and closed-loop synchronization holonomy in the disformal/non-exact sector. The measured GNSS covariance is interpreted primarily as a conformal-sector A(ϕ) observable; it should not be identified with the non-zero closed-loop holonomy predicted only by the non-exact/disformal sector. The empirical evidence presented here therefore supports the conformal sector alone; the disformal-sector holonomy — the framework's definitive discriminating test — remains unmeasured at the canonical amplitude (the J0437 phase-closure detection attests non-exact transport but sits roughly two orders of magnitude above the canonical connection's ceiling on that sightline). Evidence Tiers: A (Paper 0/EXP): Foundational theory with pre-specified search ranges; B (GNSS I–III): Primary timing-network evidence; C (SLR): Optical-domain bound and exploratory correlation test; D (H₀, COS, WB): Environmental response tests; E (JWST, GL, UCD, RBH): Conditional astrophysical/cosmological extensions. These evidence channels are not statistically independent: the exponential decay, anisotropy, orbital coupling, preferred-axis recovery, nutation, and planetary-catalogue signatures share the same GNSS data and phase-alignment method, while the SLR check, pulsar spin-down, and host-potential H 0 correlation use distinct data — the effective number of independent tests is smaller than the number of listed results, and smaller still at channel level: the SLR check reads the same conformal clock sector through different hardware, the CMB axis is product-unstable, and the pulsar channel is exploratory (its standalone evidence uninformative — mildly GR-favouring — in the J0337 adjudication, where the kinetic impedance solved in Paper 0 suppresses the differential acceleration ratio to η ≈ 2.5 × 10⁻⁷, well within the experimental bound 2.6 × 10⁻⁶). Critically, the conformal sector responsible for clock-rate modulation is not directly constrained by photon–graviton differential-propagation tests such as GW170817, although local conformal-gradient/source-charge sectors remain constrained by PPN, clock, and equivalence-principle tests. If validated through independent replication, TEP predicts that the dynamical phenomenology attributed to dark matter is the observational projection of Temporal Topology and Temporal Shear, not evidence for a new particulate matter component. In this interpretation, GNSS-scale covariance, dynamical phantom mass — with lensing tracking the gravitating matter (Papers 4, 19) — and galactic mass discrepancies are different observational projections of the same time-field topology, with Papers 4 and 13 testing distinct optical and dynamical projections of the environmental-response architecture (Paper 4: the dynamical–lensing mass difference as a phantom-mass test; Paper 13: the SPARC-anchored wide-binary transition), not a derived Earth-to-galaxy covariance-to-shear transfer kernel. The phantom-mass prediction is conditional on the Temporal Shear surviving unscreened at galactic scales — an amplitude whose deep-regime derivation from the action is identified in Paper 0 as the decisive open calculation — and is presented as a conditional consequence, not an established result. The 4,000 km correlation on Earth and galactic-scale phenomenology are interpreted within TEP as manifestations of the same environment-dependent Temporal-Topology response, which provides continuous Temporal Topology relaxation from deep potential wells to the weak-field regime. The saturation-density scale ρ T ≈ 20 g/cm³, and the diagnostic radius R T = (3M/4πρ T ) 1/3 it defines for each body, remain testable through interplanetary missions. Explicit falsification criteria include: failure of independent groups to replicate the raw carrier-phase signal; Temporal Topology correlation length falling outside the 500–20,000 km range; adjudication of the ephemeris-artifact question — GNSS-product provenance via the Satellite Laser Ranging check, and the wider orbit-projection confound via the physical-confound forward model or an orbit-free clock channel; and null synchronization holonomy in closed-loop triangular time-transfer experiments.
