Historical Solar Eclipse Validation — Pure-Tidal + α(t) ΔT vs the Documented Record

This page is the empirical confirmation that the model’s predictions actually hold against the solar eclipse record. For the parallel lunar eclipse validation (which has higher temporal resolution because lunar eclipses don’t depend on geographic localization), see Lunar Eclipse Validation. For the closed-form ΔT formula itself, see Timekeeping & Delta-T; for the Moon polynomial’s implementation and parallax limit at modern epochs, see Sun, Moon & Planets; for the broader validation tradition (Wells 1963, Cheng 2016, etc.), see Supporting Evidence.

1. Thesis

Pure-tidal Farhat-based Moon orbital evolution + a viscoelastic α(t) GIA correction from satellite gravimetry produces a ΔT formula that explains documented solar-eclipse visibility across at least 2,400 years of historical records, without phenomenological fitting to the eclipse data itself.

The conventional interpretation of the long-term ΔT record holds that pure tidal physics alone cannot account for the observed eclipse data — that an additional non-tidal Earth-rotation component is required. Stephenson’s empirical polynomial reproduces the observations and is consistent with this non-tidal component being present. The conventional Munk-MacDonald magnitude for this component is ~5-6 ms/century, attributed to glacial isostatic adjustment plus core-mantle coupling.

The model takes a more careful position. Pure tidal physics alone, derived from the Farhat 2022 deep-time evolution model applied via angular momentum conservation, fits the solar-eclipse visibility record at least as well as Stephenson. The model adds a single physically-derived non-tidal channel — a viscoelastic α(t) GIA correction from independent satellite gravimetry (Cox & Chao 2002 dα/dt + Peltier ICE-5G(VM2) multi-mode mantle rheology) — that does apply, but with much smaller magnitude (~0.6 ms/century) than the conventional Munk-MacDonald assumption.

The higher-resolution lunar timing test additionally identifies a smaller fractional non-tidal secular rate ~0.5 ms/century in the residual beyond α(t) (about 2× Cox-Chao’s satellite baseline; ~10% of full Munk-MacDonald), detected but not currently modelled. The full Munk-MacDonald-scale assumption is rejected; the dominant GIA-scale channel is included; the fractional channel is quantitatively acknowledged (see Lunar Eclipse Validation §1 for the quantitative statement and §6 for the three-component decomposition of the medieval residual).

2. Three independent validations

The validation runs in three layers, each independently informative.

Layer 1 — Moon polynomial vs NASA in TT space (foundation)

The simulation’s Moon position uses Meeus Ch. 47 (60 longitude + 60 latitude perturbation terms) on top of a 5-layer geometric precession hierarchy. For 11 canonical eclipses from NASA’s Five Millennium Catalog of Solar Eclipses spanning -524 to 985 CE, the model’s computed time of Moon-Sun conjunction was compared to NASA’s published Terrestrial Dynamical Time (TD) of greatest eclipse.

The comparison is ΔT-independent: in TT-space, the astronomical event is fixed regardless of which ΔT either side assumes. Any residual is purely a Moon polynomial accuracy question.

This is the expected Meeus Ch. 47 polynomial residual at these timescales (≈ 0.13° in Moon ecliptic longitude at year 980, the worst case). The polynomial is sound at every epoch tested.

Layer 2 — Per-event same-day conjunction check

For 19 documented historical solar eclipses with known calendar date and observation site, the model’s Moon-Sun ecliptic longitude separation Δλ was evaluated at noon UT on the documented date. The question: is there a conjunction within ±12 hours of that noon (same calendar day)?

Result: 19/19 same-day matches. Every documented eclipse date corresponds to a conjunction in the model. The visibility-at-site question — does the model’s eclipse path actually reach the observation site? — is the subject of Layer 3 below.

Layer 3 — Visibility window: pure-tidal vs Stephenson (the headline)

For each event, the question becomes: what range of ΔT values would put the model’s eclipse path within penumbra reach (< 7,500 km) of the observation site? And does either (a) the model’s pure-tidal ΔT or (b) Stephenson’s empirical ΔT fall inside that range?

The model wins both the visibility count and the mean-residual test. The penumbra count is the headline: if Stephenson’s empirical curve were “the truth”, it should clearly win the broad visibility test — Stephenson’s coefficients were calibrated to make eclipses visible at observed sites. Instead, the model explains more of the visibility, not less.

The two events where Stephenson loses penumbra are Ibn Yunus 979 May 28 and 1004 Jan 24 — both fail because Stephenson’s ΔT is too low at those medieval epochs, pushing sub-solar east of Cairo. The model’s higher ΔT recovers them.

The umbra-count tie at 6/13 each is the noise floor — totality strips are narrow (~270 km wide) and ancient localization is too coarse to discriminate at sub-100 s ΔT precision.

3. The ΔT gap — pre- and post-α(t) framing

Pre-α(t) pattern (the original finding)

When this validation was first performed (with pure-tidal-only physics, no GIA correction), the model showed a constant linear ΔT excess over Stephenson’s empirical fit, scaling at ~2 s/yr into the past:

The linear-in-time slope of ~1.96 s/yr corresponded geometrically to a constant Length-of-Day difference of ~5–6 ms between pure-tidal-only and Stephenson — the order of magnitude of canonical non-tidal Earth-rotation estimates (the Munk-MacDonald mechanism).

At the time, two readings of this gap were possible: either a real non-tidal Earth-rotation component (Munk-MacDonald-scale), or a phenomenological feature of Stephenson’s fit. The solar-eclipse visibility test couldn’t distinguish them.

Post-α(t) pattern (the resolved finding)

The higher-resolution lunar-timing test (doc 102 , see also Lunar Eclipse Validation) resolves the tension. The model now includes a viscoelastic α(t) GIA correction from independent satellite gravimetry, and the constant linear excess is absorbed in the ancient era:

The constant-LOD interpretation no longer applies: ancient-era ΔTs now agree with Stephenson to within ~50 s. A bump-shaped residual in the medieval window (~1,000 s peak in the 840–1020 CE window; exact peak year is reference-conditional) is the remaining structured signal — characterized in detail in the lunar validation page and decomposed there into a framework-native 8H/1851 lattice harmonic plus a fractional non-tidal secular drift plus observation noise.

The reading that survives: the non-tidal contribution IS real, and decomposes into a dominant GIA-scale channel (~0.6 ms/century, captured by α(t)) plus a smaller fractional non-tidal secular rate (~0.5 ms/century, detected in the lunar-timing residual but not modelled) — not at the full Munk-MacDonald magnitude (~5-6 ms/century) the pre-α(t) framing originally suggested. The lunar timing test, which has ΔT resolution of minutes rather than hours, distinguishes them.

4. What the validation establishes

What the model can claim with confidence:

1. The Moon polynomial used in the simulation is validated against NASA’s JPL reference at ±15 min back to 2,500 years before J2000 (n = 11 events, -524 to 985 CE).
2. The model’s pure-tidal + α(t) GIA ΔT formula explains documented solar-eclipse visibility for 19/19 events spanning -762 to 1654 CE.
3. The full Munk-MacDonald-scale (~5-6 ms/cy) non-tidal Earth-rotation postulate is rejected by the historical eclipse record. A dominant GIA-scale (~0.6 ms/cy) channel is included via the α(t) correction, and a smaller fractional non-tidal secular rate ~0.5 ms/century is detected in the residual by the higher-resolution lunar timing test (see §6 for the three-component decomposition).
4. The model’s deep-time grounding is independent of the eclipse record — anchored to Wells 1963 (Devonian coral growth bands), Wu et al. 2024 (650-Myr cyclostratigraphy), modern Lunar Laser Ranging, and Cox & Chao 2002 satellite gravimetry. None of this evidence is circular with the historical eclipses.

What the model is not claiming:

• That Stephenson’s curve is wrong. It fits the eclipses well too.
• That all non-tidal mechanisms are absent. A dominant GIA-scale (~0.6 ms/cy) channel is included via the α(t) GIA correction, and a smaller fractional non-tidal (~0.5 ms/cy) is detected in the lunar timing test residual but not currently modelled. The rejected claim is specifically the full Munk-MacDonald magnitude (~5-6 ms/cy).
• That solar-eclipse data alone settles every question. Lunar-eclipse timing — which doesn’t depend on geographic localization and has minute-scale ΔT resolution — is the stronger constraint and is covered in the companion Lunar Eclipse Validation page.

5. Limits

Caveats the reader should keep in mind:

1. n = 19 historical solar eclipses is small (originally 20 — one prior test entry labelled “Cambyses solar eclipse” was removed when Stephenson 1997  confirmed the Babylonian astronomical-diary references from Cambyses II’s reign are lunar, not solar). Statistical power for discriminating models with sub-100 s ΔT differences is limited.
2. Geographic localization of ancient eclipse paths has irreducible uncertainty. The 4,500 km umbra reach and 7,500 km penumbra reach are approximations based on the model’s sub-solar-point distance.
3. Ancient observation sites often have a latitude error larger than the umbra reach (e.g., 977 Dec 13: Cairo is 3,000 km north of the umbra path regardless of ΔT). For these events the test only tells us about penumbra visibility, not totality.
4. The visibility-window methodology cannot detect ΔT differences smaller than ~50 s because the penumbra window is wide (~30,000 s typical). For finer ΔT discrimination, lunar-eclipse timing-at-site is the stronger constraint.

6. Umbra-centerline tightening — ★ TOTAL matches at conventional documented dates

The diagnostic effort identified an ad-hoc visualization-layer overlay in the simulation’s scene-graph Earth rotation that had been adding a ΔT × 2π/86400 rotation correction on top of standard GMST(UT1) Earth rotation. The overlay made the visible umbra disc render at Stephenson 2016’s documented locations but introduced a systematic offset relative to the model’s underlying pure-tidal physics. Once the overlay is disabled, the pure-tidal physics itself produces umbra centerlines passing through the documented observation sites at ★ TOTAL for the headline events above. The 19/19 penumbra result from §2 is unaffected because that test uses a Meeus-based sub-solar function independent of the scene-graph Earth orientation.

The one persistent residual: -135 Babylonian (1,159 km)

A single deep-time event refuses to match at the umbra-centerline criterion: the 15 April 136 BCE (= -135 astronomical) Babylonian eclipse, recorded in the Babylonian astronomical diaries BM 45745 and LBAT 1285. This is one of the most scholarly-secure attributions in the historical eclipse corpus (four-planet astronomical fingerprint, double-dated Arsacid/Seleucid eras, re-confirmed by Stephenson & Steele 2006). The model’s predicted umbra centerline lands ~1,159 km south-east of Babylon, over central Saudi Arabia (~24.8°N, 52.3°E).

Diagnostic decomposition of the 1,159 km gap traces it to three sources:

• ~270 km from the model’s ΔT being +669 s over NASA’s value (closable to ~80 km by tuning α(t) within Peltier ICE-5G vs ICE-6G uncertainty, but the empirical sensitivity is only ~3.3 km per 100 s of ΔT change — α(t) tuning alone is mathematically incapable of closing the full gap)
• ~440 km from Meeus Ch. 47 polynomial residual in Moon ecliptic latitude β at this specific JD (β_Meeus = 0.728° vs the more accurate β ≈ 0.66° from ELP-2000/82, a 0.07° polynomial residual amplified by the d_M / R_E ≈ 60 leverage on Earth’s surface)
• ~450 km from other Meeus polynomial terms (Lp residual, Sun constant residual, geometric ray-trace approximation)

The dominant component (~890 km) is Meeus Ch. 47 polynomial residual at this specific JD, characteristic of the truncated 60-term Meeus series summing constructively at a specific JD even when the full 37,000-term ELP-2000/82 series would cancel. NASA Five Millennium Canon (which uses ELP-2000/82) places greatest at (46.8°N, 58.9°E) — itself ~2,100 km from Babylon, with the path crossing Babylon at a non-greatest moment via ELP-2000/82’s more accurate Moon polynomial.

The interpretation: the model’s prediction (centerline over Saudi Arabia, deep partial at Babylon at ~95–99% magnitude) is the correct prediction given its Meeus polynomial. Babylon at 1,159 km from centerline would have observed a dramatically darkened sky consistent with the diary’s record of Venus, Mercury, and “Normal Stars” visible. Closing the gap to NASA’s path-through-Mesopotamia would require replacing Meeus Ch. 47 with a higher-precision lunar theory — documented as a forward-path proposal in the simulation repo, but not currently prioritized. This is consistent with IMCCE’s own published acknowledgement  that polynomial-precision limits at deep historical past are real and finite. Full diagnostic decomposition + per-scale α(t) sweep data + external references on the eclipse path are in the simulation repo’s doc 103 .

7. Reproducing the validation

The complete validation suite is available in the simulation as developer-mode console tests at Console Tests (F12) > Historical Eclipses & ΔT (10 buttons): the NASA catalog cross-check (Layer 1 above), the per-event same-day check (Layer 2), the visibility-window comparison (Layer 3), plus 7 supporting diagnostics. The 19 documented eclipses plus 8 modern reference events are also exposed as planetStats nav-buttons under Moon → Historical Solar Eclipses (validation), where the user can step through each event and visually verify Moon-Sun alignment in the 3D scene.

A separate group of buttons under Console Tests (F12) > Lunar Eclipses & Validation runs the higher-resolution lunar-timing track — 20 buttons in 5 subgroups (foundation, predictive finders, NASA Canon cross-check, primary observation tests, residual investigation) covering 270 primary-source lunar observations + 89 primary-source solar observations from Stephenson 2016. See the Lunar Eclipse Validation page for the results.

Full reproducibility notes, methodology, and underlying numerical inputs are in the simulation repo’s doc 101  (solar) and doc 102  (lunar + α(t) GIA physics). A companion analysis using a complementary residual-RMS methodology (Moon–Sun geocentric separation in degrees, evaluated across a broader 35-eclipse set) is in doc 100 .

Continue to Lunar Eclipse Validation for the higher-resolution test on 270 primary-source lunar observations, or Predictions for the model’s testable predictions across near-term, medium-term, and deep-time horizons.