Lunar Eclipse Validation — α(t) GIA Correction vs 270 Primary-Source Observations

This page is the canonical home for the α(t) viscoelastic GIA physics and the lunar-timing validation. For the closed-form ΔT formula itself, see Timekeeping & Delta-T; for the Moon polynomial mechanism, see Sun, Moon & Planets; for the solar visibility test, see Historical Solar Eclipse Validation; for the broader validation tradition (Wells 1963, Cheng 2016, etc.), see Supporting Evidence.

1. Thesis

The lunar-timing record requires a non-tidal Earth-rotation contribution whose dominant component matches the GIA magnitude measured independently by satellite gravimetry, plus a smaller fractional non-tidal secular rate (~0.5 ms/century) beyond that. The full Munk-MacDonald (~5-6 ms/century) postulate is rejected by the historical record; the fractional non-tidal channel is quantitatively acknowledged.

The conventional treatment of long-term Earth rotation attributes the non-tidal component to glacial isostatic adjustment (GIA) plus core-mantle coupling, conventionally at ~5-6 ms/century (the Munk-MacDonald estimate). Stephenson’s empirical polynomial reproduces the historical eclipse record and is consistent with such a component being present.

The model’s claim is more precise. The framework’s α(t) captures the dominant GIA-scale portion (Cox & Chao satellite value, ~0.23 ms/cy) with viscoelastic time-dependence from Peltier ICE-5G(VM2). Follow-up analysis of the residual identifies a smaller fractional non-tidal secular rate of ~0.5 ms/century beyond that — approximately 2× the Cox-Chao GIA-only baseline, and about 10% of the full Munk-MacDonald postulate. The full postulate is rejected; a fractional non-tidal channel is present. Candidate mechanisms for the fractional contribution include time-varying mantle-core coupling (a constant-rate version was tested and rejected — see §4) and continental hydrology on centennial scale. The framework’s α(t) is derived from independent satellite gravimetry, with zero parameters fitted to eclipse data in the live model.

The empirical test is the higher-resolution lunar-timing comparison below. Lunar eclipses are visible across Earth’s entire night-side hemisphere, so geographic localization is irrelevant — the constraint reduces to timing of opposition at the observation site, resolving ΔT to within minutes per event rather than the ~50-100 s resolution of the solar visibility test (historical-eclipse-validation).

2. The 270-observation test

The dataset is 270 timed lunar observations from Stephenson, Morrison & Hohenkerk (2016) — Babylonian, Greek, Chinese, and Arab observers spanning -720 BCE to 1280 CE, drawn from supplementary tables S01 + S02 + S04 + S05 + S07 + S09. For each observation, three ΔT predictions are compared:

• the value Stephenson derives from the original observation,
• NASA’s empirical Espenak/Meeus polynomial value (averaged over the observation year),
• the model’s pure-tidal Farhat + α(t) GIA viscoelastic prediction.

The headline metric is computed over the 267 of those 270 events for which all three ΔT predictions are defined — the remaining three fall outside the NASA polynomial’s published validity range or have missing per-event ΔT in Stephenson 2016, so they cannot participate in the three-way comparison.

Headline (267 events with all three ΔT defined)

                                       Mean |residual| (s)    Mean |residual| (min)
NASA Espenak/Meeus ΔT:                            1199                20.0
Model pure-tidal + L1-orbital α(t):               1604                26.7

Events where model closer to obs than NASA: 78/267 (29.2%)
NASA closer to obs by: 25.2% on average

NASA’s polynomial is FIT to (essentially) this exact observation dataset; ours PREDICTS it from the named independent physical constants in §4. The 6.7-minute gap between model and NASA is the model’s distance from the empirically-fitted polynomial on a ~20-minute observation noise floor under the L1-orbital-coupled α(t) refinement (see doc 99 §“Deep-time refinement” for the α(t) formulation).

Adding the α(t) GIA correction to the pure-tidal model reduces the linear-in-time R² of the residual from 0.53 to 0.36 — the organised structure that pure-tidal physics leaves behind is largely absorbed once the non-tidal GIA contribution is included. What remains in the residual is a mix of observation noise, a linear secular trend (~1.9 s/yr going deeper in the past), and quasi-periodic content at solar-activity timescales (see §6 for the periodic-forcing analysis).

3. Cross-cultural agreement — the strongest evidence

More striking than the headline residual is the agreement across observation traditions. After detrending the small remaining linear slope (−0.77 s/yr), the four independent traditions in the dataset — Babylonian, Greek, Chinese, Arab — agree on the magnitude of the model’s residual to within ±200 s. A regional observational bias would show up as a strong source-specific mean; a real model–data residual would show consistent magnitude across sources. The data shows the latter.

Per-table cross-cultural consistency

The three highest-precision, largest-sample sources — S01 Babylonian (n=125), S05 Chinese (n=69), and S09 Arab (n=32) — all detrend to within ±200 s. S09 (Arab medieval — Ibn Yunus, Habash al-Ḥāsib, al-Battānī) is the tightest at ±57 s with RMS 763 s.

Four independent observation traditions, separated by thousands of years and tens of thousands of kilometres, agree on the magnitude of the model’s residual to within the noise floor. Cross-cultural agreement at this level cannot be accidental — it confirms the residual is a real property of the model–data fit, not a regional observational bias.

Three thousand years deep — the Babylonian convergence

The deepest, hardest-to-fit observations — cuneiform tablets from Babylon spanning -800 to -300 BCE — converge to within ~2 minutes per century:

Three thousand years deep, observations from clay tablets, reproduced by a model whose every physical constant comes from independent literature — zero fitting parameters. That this works is the headline.

4. How it works — the α(t) GIA correction

Long-term Earth-rotation evolution has two distinct channels, acting in opposite directions. The tidal channel — the Moon raises tides on Earth, the tides drag Earth’s bulge slightly ahead of the Moon, the Moon recedes outward and Earth’s spin slows — is the dominant secular contribution and is already captured by the Farhat 2022 lunar-distance evolution model used elsewhere in the framework. The non-tidal channel is what α(t) handles — α being Earth’s polar moment coefficient, the dimensionless number that captures how much of Earth’s mass sits along the rotation axis. This channel acts the other way: it speeds rotation up slightly. The two channels partly cancel; what’s left over after subtracting the dominant tidal slowing is the small residual signal recorded in the historical eclipse data.

The physical mechanism is glacial isostatic adjustment (GIA): 12,000 years after the last glacial maximum, continents that once carried massive ice sheets — Scandinavia, Hudson Bay, the Antarctic margins — are still rebounding upward. As mass redistributes from former-equatorial ocean basins back toward polar continents, Earth’s moment of inertia decreases. By conservation of angular momentum, the rotation rate ticks up slightly, and the day shortens. This shows up in the historical eclipse record as a measurable shift in ΔT.

Critically, GIA is a purely Earth-internal mass redistribution. It does not transfer any angular momentum to the Moon, so the Moon-distance evolution (Farhat 2022) and Kepler’s third law for the Moon’s orbit are completely untouched. Only the partition of angular momentum within Earth shifts; the total Earth-Moon system angular momentum is conserved exactly.

The relaxation is not instantaneous. Earth’s mantle behaves viscoelastically — each layer (upper mantle, transition zone, lower mantle) relaxes back to equilibrium on its own timescale, from about 1,500 years for the upper mantle out to about 14,000 years for the lower mantle. The model uses the standard three-mode decomposition from Peltier 2004 (ICE-5G(VM2)), with the transition-zone mode (≈ 5,000 yr) dominating the response in the historical-record window.

All the inputs come from independent measurements: the modern value of α from IERS gravity-field observations, its rate of change from satellite laser ranging (Cox & Chao 2002), and the three mantle-layer timescales from the standard Peltier GIA model. Nothing here is fitted to the eclipse data — those numbers are fixed by independent geodesy and rheology, and the eclipse record is the test.

α(t_age) = α_J2000 + Σᵢ Δαᵢ · (1 − exp(−t_age / τᵢ))

The deep-time refinement — α(t) as a climate signal

The three-mode viscoelastic form above passes the ±3 kyr historical-eclipse test cleanly. Extending the same form into deep-time work — hundreds of thousands to millions of years — exposes a limitation. The modes saturate symmetrically around J2000: α evolves the same way going past or future, which produces a small slope discontinuity in α at year 2000. The discontinuity shows up as a visible kink in the length-of-day curve across the year-2000 line in the ±12 kyr Formula Verification view.

The physical reason for refining the form is deeper than the kink. Glacial isostatic adjustment from the last glacial maximum is only part of α’s evolution. The larger, longer-timescale driver is the glacial cycle itself — ice sheets build up and melt on ~100 kyr timescales driven by orbital forcing, and that ice-mass redistribution changes α at the same period. The dα/dt rate Cox & Chao measured at J2000 is one snapshot of that cycle, not the whole story.

The refined α(t) binds to the same L1 orbital layer that already fits the LR04 δ¹⁸O record over the last 1 Myr in the Climate Formula. One mechanism, two observables: the L1 orbital signal drives both the ice-volume proxy in the sediment record and the α-driven length-of-day oscillation.

α(t) = α_J2000 − k · [ L1(year) − L1(2000) ]

Here L1(year) is the L1 orbital layer of the Climate Formula (the sum of cosines on H-lattice divisors of 8H ≈ 2.68 Myr) evaluated in δ¹⁸O ‰ units, and k = −5.24 × 10⁻⁷ per ‰ is calibrated so that dα/dt at J2000 matches Cox & Chao’s −1.8 × 10⁻¹¹ /yr exactly. That single scale factor is the only new parameter, and it is not fitted to eclipse data — it is fixed by the satellite gravimetry measurement.

The historical-eclipse validation is preserved to first order. The refined form preserves every property the three-mode Peltier form had within the ±3 kyr eclipse window:

• α(J2000) = 0.3306947 exact — IERS Conventions 2010 anchor preserved.
• dα/dt(J2000) = −1.8 × 10⁻¹¹ /yr exact — Cox & Chao anchor preserved via the k calibration.
• The integrand over the ±3 kyr eclipse window differs from the three-mode form by less than a percent per year, because the L1 signal near J2000 is dominated by the same low-frequency components that the three-mode timescales captured.

The 26.7-minute mean |residual| against 267 primary-source observations, the 6.7-minute gap to NASA’s polynomial, and the solar cross-validation results in §5 remain valid. This is a refinement of α’s shape at deep time and at the ±100 kyr Milankovitch band — not a replacement of the historical-eclipse physics. The small increase in mean residual under L1-orbital α(t) (vs 24.4 min under the |t|-symmetric form) reflects the slightly different α(t) trajectory at deep past — a systematic that emerges over 2,000+ years of ΔT integration — and is more than compensated for by (a) the removal of the derivative discontinuity at J2000 and (b) the new physically-motivated glacial-cycle-driven α trajectory suitable for deep-time work.

What the refinement adds: a testable ~100-kyr LOD oscillation

Under the refined form α oscillates with the ~100 kyr, ~41 kyr, and precession-band periodicities that drive climate. This translates to a peak-to-peak length-of-day variation of about 250 ms riding on top of the smooth tidal-recession trend. The extrema line up with well-dated Marine Isotope Stages:

The underlying physical mechanism — ice mass measurably shifting Earth’s J₂ oblateness on human-observable timescales — has been directly confirmed by Cheng, Tapley & Ries 2011 . Their LAGEOS satellite J₂ record transitions from a linear GIA-driven decrease to an accelerating trend around 1998, attributed to polar ice-sheet mass loss. That decadal observation validates the “ice → J₂ → α → LOD” link at satellite timescales. Extrapolating the same mechanism to the ~100-kyr Milankovitch cycle is a physics-based extrapolation, not a direct observation — modern satellite records span only ~50 years, too short for the glacial-cycle band. A definitive test would require either paleoclimate LOD indicators sensitive enough to resolve the ~250 ms band on top of the tidal-recession trend, or a much longer future satellite record.

The Fibonacci Universe simulation exposes this directly. In the Formula Verification → Solar Day Length view, an Export Cycles button renders the model’s LOD prediction over −248 kyr to +102 kyr with the MIS peak-age labels above aligned to the L1 orbital extrema.

Why the refinement matters — the unifying claim

The refined α(t) elevates a small implementation detail into a unifying claim about three coupled Earth systems:

Planetary gravity, Earth’s climate, and Earth’s rotation are not three independent systems. They are one system, connected by ice mass as the mediator.

Standard geophysics treats orbital dynamics, climate, and rotation as coupled only through the tidal channel — the well-understood tidal-friction chain that recedes the Moon and slows Earth’s rotation. The refined form of α(t) additionally invokes an orbital-eigenmode → ice → α → LOD channel, slower (100-kyr timescale), indirect (mediated by ice mass on Earth’s surface), and predicting an observable phase lock between length-of-day and the LR04 δ¹⁸O record. This coupling is included as a testable prediction in Predictions §13.

5. Solar cross-validation: same physics, different dataset

ΔT is a property of Earth rotation, not of the eclipse type — so the same α(t) physics that fits the lunar timing record should fit the solar timing record independently. The cross-validation runs the same three-way comparison pipeline against 89 primary-source solar observations from Stephenson 2016 supplementary tables S03 (Babylonian solar, 25 events), S06 (Chinese solar, 42 events), and S08 (Arab solar, 22 events), spanning -356 BCE to 1277 CE.

                                       Mean |residual| (s)    Mean |residual| (min)
NASA Espenak/Meeus ΔT:                             672                11.2
Model pure-tidal + L1-orbital α(t):                994                16.6

The absolute residuals are smaller than the lunar test (NASA 672 vs 1199 s; model 994 vs 1604 s) because solar observations have a tighter intrinsic precision — narrow totality paths give sharper timing. The model is 32.4% further from observations than NASA on average — a wider relative gap than the lunar 25.2% — because tighter observations expose the residual structure more visibly. In absolute terms, the model’s solar residual (994 s) is closer to observations than its lunar residual (1604 s), consistent with the per-observation solar timing being intrinsically sharper.

The per-century breakdown shows the same medieval overshoot direction in both datasets — confirming that ΔT is a property of Earth rotation, not an artifact of eclipse type. This is the type-independence requirement passing.

For the year 1000-1099 century specifically — where the medieval residual structure peaks — the model is closer to observations on 44% of solar events (vs 29.2% globally across L-7), indicating relative model strength in this era despite the medieval overshoot. This is direct independent evidence that the medieval-era residual structure is a real signal common to both eclipse types.

6. The medieval residual

A residual bump remains after the α(t) correction. Under the L1-orbital-coupled α(t) refinement, its shape is characterized as a ~1000 s peak in the 840–1020 CE window (the exact peak year is reference-polynomial-dependent — see below), FWHM ~660 yr, visible in both the lunar and solar records. This is the medieval residual.

Follow-up analysis decomposes the residual cleanly into three physical components:

1. A framework-native millennial-scale 8H lattice harmonic at n=1851 = 73 × Jupiter-Saturn synodic = 1449 yr. The period is a zero-fit structural prediction — it drops out of the framework’s 8H arithmetic with no calibration against eclipse data. An 8H integer-divisor scan across the sub-Milankovitch band identifies n=1851 as the top-ranked structural match at 0.001% error (73 × J-S synodic exact). Empirical fit against the eclipse residual gives R² = 0.975 in-sample, R²_test ≈ +0.97 cross-validated on CE→BCE prediction. Live integration into the model is deferred pending independent (non-eclipse) amplitude/phase calibration from a paleoclimate proxy such as Bond 1997 IRD or a SIM-driven solar-activity reconstruction. A research toggle in the sim demonstrates the fit but is OFF by default specifically to preserve the zero-eclipse-fitting claim.

2. A fractional non-tidal secular rate of ~0.5 ms/century. After the lattice-harmonic accounts for the “bump,” the remaining drift has a quadratic component consistent with a secular non-tidal rate about 2× the Cox-Chao satellite-measured GIA value, and about 10% of the full Munk-MacDonald postulate. The full Munk-MacDonald rate (~5 ms/cy) is decisively rejected by the ancient BCE data — the constant-Holme extrapolation over-corrects Babylonian ΔT by ~2,700 s. A fractional (~0.5 ms/cy) time-averaged non-tidal contribution IS present and is a candidate for time-varying mantle-core coupling or continental hydrology.

3. Observation noise and small artifacts. ~60 s RMS floor after both above corrections. This is essentially the Stephenson dataset’s per-event precision averaged into the sampled residual curve — the irreducible floor.

The reference-polynomial peak-year caveat: the exact peak year is reference-dependent — Stephenson 2016 places it near 1020 CE, a NASA-derived catalog average places it near 840 CE. The gross bump shape correlates at r = 0.82 between the two references and RMS magnitudes agree (890 s vs 907 s), so the “there is a broad medieval bump” claim is robust. Only the exact peak-year and peak-magnitude numbers are fragile.

Eight statistical hypotheses were tested formally under L1-orbital α(t), plus two follow-up predictive tests of proposed physical mechanisms. Under a per-era stability check, all correlation-based hypotheses reduce to drift-tracking artifacts — their aggregate correlations exist but reflect ancient-BCE monotonic co-variation rather than causal per-observation links. The structural mechanisms that survive are all spectral or lattice-native:

Combining these tests with the three-component decomposition above: the medieval residual is structurally explained — a lattice-native 73×J-S synodic harmonic accounts for the bump; a small fractional non-tidal secular rate accounts for the residual drift beyond α(t); observation noise accounts for the rest. Correlation-based hypotheses that appeared strong under aggregate analysis do not survive per-era stability checking.

(Full statistical methodology — Bonferroni multiple-comparison correction, Lomb-Scargle FAP thresholds, jackknife robustness, white-noise null comparison, per-era stability check as a standard filter against drift-tracking artifacts, higher-order polynomial audit for fit-order artifacts — is in the simulation repository. See doc 102  for the complete analysis.)

7. What the validation establishes — and what it does not

Established

1. The non-tidal Earth-rotation contribution IS real and detectable in the historical lunar record. Lunar-timing resolution (sub-100 s ΔT per event) discriminates this signal where solar-eclipse visibility (50-100 s) cannot.
2. The non-tidal contribution decomposes into a dominant GIA-scale piece (captured by α(t)) plus a smaller fractional non-tidal secular rate detected in the residual; the full Munk-MacDonald postulate is rejected (see §1 for the quantitative statement).
3. All physical constants in the live model come from independent literature, zero eclipse-fitting parameters — IERS α + Cox & Chao dα/dt + Peltier ICE-5G(VM2) multi-mode decomposition + Climate Formula L1 orbital layer for the deep-time refinement. The live model predicts NASA’s empirical polynomial to within 6.7 min on a 20 min observation noise floor.
4. Earth-Moon angular momentum and Kepler’s 3rd law preserved exactly. α(t) is purely Earth-internal mass redistribution; the Moon orbit chain (Farhat 2022) is untouched.
5. Four independent observation traditions (Babylonian, Greek, Chinese, Arab) agree on the magnitude of the model’s residual to within ±200 s after detrending — the cross-cultural validation argument.
6. The medieval residual has structure that decomposes cleanly into three named physical components — a framework-native 8H/1851 lattice harmonic, a fractional non-tidal secular rate, and observation noise (see §6 for the full decomposition and spectral evidence at solar-activity periods).

Not claimed

• That NASA’s polynomial is “beaten.” NASA is closer to the observations by 25% on average. NASA’s polynomial is FIT to this dataset; ours PREDICTS it. The achievement is “predicting historical eclipse timing to within a few times the observation noise floor using only first-principles physical constants, with zero coefficients fitted to eclipse data in the live model” — not beating it.
• That the 26.7 min model residual is purely physical. The Stephenson 2016 dataset has a ~20 min irreducible per-observation scatter; the remaining 6.7-min gap to NASA includes both observation noise and small contributions from non-tidal channels not fully modelled (the millennial-scale lattice harmonic and the fractional non-tidal secular rate documented in §6).
• That α(t) GIA is the only non-tidal channel. A fractional secular non-tidal rate is present beyond α(t) (§6, component ii); its physical channel is unmodelled and is dedicated follow-up work (see §9(a)).
• That correlation-based hypotheses (H3 mass-balance, Path A solar-activity direct coupling, Test 5 planetary perihelion) have been shown to explain the medieval bump. Their aggregate correlations exist as statistical facts but are drift-tracking artifacts under per-era analysis — they do not causally explain the residual on a per-observation basis.

8. Limits

1. Stephenson 2016 per-observation noise is ~20 min RMS, dominant at the per-event level. Neither model can do better than this; the 6.7-min model-vs-NASA gap is the structural disagreement on top of the noise floor.
2. The medieval residual (years 800-1300, ~1000 s peak in the 840–1020 CE window under L1-orbital α(t); exact peak year is reference-conditional) has its structure fixed by the three-component decomposition of §6; the framework-native lattice component is not yet live in the model, and the fractional non-tidal channel is not yet modelled.
3. The Cox & Chao 2002 satellite measurement is a modern-era value (satellite era ~1979-present). The lunar-timing cross-validation tests whether this modern rate, viscoelastically extrapolated back over 2,000+ years and refined via the L1 orbital layer of the Climate Formula, matches the historical record — the dominant GIA component matches; the residual gap beyond it is component (ii) of §6.

9. Open question — the physical mechanism of the fractional non-tidal channel

Under the L1-orbital α(t) refinement and the follow-up drift-decomposition diagnostics, the structure of the medieval residual is well-characterized (see §6). Two open questions remain:

(a) Which physical channel produces the fractional non-tidal secular rate identified in §6 (component ii)? The residual demands its presence but the framework’s α(t) does not currently include it. Candidate mechanisms:

• Time-varying mantle-core electromagnetic coupling. The constant-Holme (1998) modern-era rate over-corrects Babylonian ΔT decisively, but a time-variable coupling whose 2000-year time-average is ~0.5 ms/cy is quantitatively consistent with the diagnostic. Requires derivation of the time variability from independent geomagnetic observations, not fitting to eclipses.
• Continental hydrology / groundwater / sea-level redistribution on centennial timescales. Not currently in the framework’s α(t) 3-mode decomposition. Measurable in modern satellite era; extrapolation over 2000+ years would require an independent hydrologic reconstruction.
• Regional GIA structure beyond the global three-mode average. A higher-resolution ICE-6G_C-type model with continental-resolution rebound profiles could absorb part of this rate. Calibrating against eclipse data would forfeit the first-principles independence the framework treats as load-bearing.

(b) Independent amplitude/phase calibration for the framework-native 8H/1851 lattice harmonic identified in §6 (component i). The period is a zero-fit prediction; the amplitude and phase would need to be calibrated against an independent (non-eclipse) source such as Bond 1997 IRD, Braun 2005 thermohaline reconstruction, or a SIM-driven solar-activity model calibrated on satellite-era data. Once independent physics gives amplitude/phase, the correction can be turned on live without compromising the zero-eclipse-fitting claim.

Spectral evidence at solar-activity periods (Gleissberg 88 yr, Jose 179 yr, de Vries 182 yr, partially-supported 14.2 yr peak — the coherent H6 forcing family; see §6 Table row H6) is consistent with SIM-mediated solar-activity forcing as the physical driver: 73 × J-S synodic dynamics modulate the Sun’s motion around the barycenter, driving solar-activity envelope modulation on multi-centennial timescales (Charvátová 1990–2007, Wilson 2013, Scafetta 2010). The direct-coupling variant (solar activity → ionospheric-thermospheric coupling → LOD, via the Holme–de Viron 2013 published coefficient) was tested predictively (Path A) and does NOT close the residual — sign inversion and per-era failure — so the mechanism is likely indirect (SIM-mediated climate response with mass redistribution) rather than direct atmospheric coupling.

Follow-up work: (i) independent amplitude/phase calibration of the lattice harmonic from a paleoclimate source; (ii) derivation of a time-variable mantle-core coupling from geomagnetic secular-variation observations; (iii) higher-resolution regional GIA modelling for the fractional non-tidal rate. The eclipse-record residual, taken seriously, becomes a measurement tool for whichever direction pays off first.

10. Reproducing the validation

The complete validation pipeline is available as developer-mode console tests in the simulation at Console Tests (F12) > Lunar Eclipses & Validation — 20 buttons covering foundation tests, predictive lunar and solar eclipse finders, NASA Canon cross-check (12,064 events), the primary observation tests (the 270-event lunar comparison + 89-event solar cross-validation), and residual-investigation diagnostics (regression, periodogram, mass-balance correlation, missing-signal shape characterization). The 14 documented modern lunar eclipses are also exposed in the tweakpane menu under Solar & Lunar Eclipses → Lunar Eclipses, where the user can step through each event with Prev/Next buttons and visually verify Moon-Sun alignment in the 3D scene.

Full reproducibility notes, methodology, and underlying numerical inputs are in the simulation repository at doc 102 .

Continue to Predictions for the model’s testable predictions across near-term, medium-term, and deep-time horizons. For the parallel solar-visibility test on 19 documented historical eclipses, see Historical Solar Eclipse Validation.