Technical Note
(Letter Of Method)

Overview

Purpose: provide a predictive, testable tide model driven by Ӕ saturation (field geometry), not “pulling forces”.

Method: define an Ӕ driver Σ at Earth’s boundary and a compliant ocean response η; fit and forecast at tide gauges.

Driver: Ӕ Saturation at the Boundary

Let x=(φ,λ) and t be time; unit surface normal \hat{r}(x).

Let n_☉(t), n_◐(t) be Sun/Moon directions (topocentric).

Angle of Encounter factors: Ψ_j(x,t) = \hat{r}(x)·n_j(t),   j∈{☉, ◐}.

Quadrupole kernel: P₂(u) = ½(3u²−1).

Ӕ driver (dimensionless):  Σ(x,t) = A_☉ f_☉(t) P₂(Ψ_☉) + A_◐ f_◐(t) P₂(Ψ_◐).

A_j = saturation weights (fit) → expect A_◐/A_☉ ≈ 2:1.  f_j(t) = slow framing (distance/zenith normalization, EMF transparency).

Response
Local Compliance of the Boundary

Station-wise, use a compact linear response capturing phase and memory:

η(t) = α Σ(t) + β dΣ/dt + Σ_k γ_k Σ(t−τ_k).   (2–5 terms typically suffice)

Interpretation: α = in-phase compliance; β = inertial/phase lead; γ_k with delays τ_k = local basin memory.

Harmonic Content and Precision

Because Σ uses P₂ on a rotating sphere, η inherits classical constituents (M2, S2, K1, O1, …).

Thus we retain mariner-grade precision while replacing the cause with Ӕ saturation geometry.

Framing Factors f_j(t)

f_j(t) = c_{j0} + c_{j1} D_j(t) + c_{j2} M(t).

D_j(t): normalized distance/zenith framing (slow).

M(t): EMF transparency (e.g., Kp, Dst, dip-latitude).

Effect: allows testing Ӕ-unique predictions (geomagnetic storm modulation, latitude patterns).

Fitting Procedure (per station)

Inputs: gauge series η_obs(t); ephemerides (n_☉, n_◐, distances); geomagnetic indices (Kp, Dst).

1. Build Σ(t) from P₂(Ψ_j) and f_j(t).

2. Regress η_obs on Σ, dΣ/dt, and 2–4 delayed Σ(t−τ_k) terms → solve {A_j, α, β, γ_k, c_{jm}} by least squares.

3. Validate on hold‑out months; compare RMS to official harmonic predictions.

4. Examine A_◐/A_☉; expect ~2:1.

Test EMF terms: amplitude/phase vs Kp/Dst (small, systematic).

Validation Metrics

• RMS error vs official harmonics.

• Amplitude/phase residuals.

• Stability of coefficients month→month.

• Ӕ‑unique effects: storm‑time modulation; eclipse‑track residual patterns; dip‑latitude amplitude structure.

Falsifiability

If Σ‑driven fits cannot match gauge precision comparable to harmonic baselines → model fails.

If Ӕ‑unique signatures (Kp/Dst modulation, eclipse residuals) do not appear at predicted levels → extras fail.

If fitted A_◐/A_☉ does not cluster near ~2 across stations → geometry claim fails.

Implementation Outline (pseudo‑code)

Given station coords (φ,λ) and time vector t:
• Compute n_☉(t), n_◐(t) (ephemerides). 

• Ψ_j(t)=\hat{r}·n_j(t). 

• Σ(t)=A_☉ f_☉ P₂(Ψ_☉)+A_◐ f_◐ P₂(Ψ_◐).

• Build design matrix X=[Σ, dΣ/dt, Σ(t−τ₁), …, framing terms].

• Solve θ=(XᵀX)⁻¹Xᵀ η_obs.

• Forecast by rolling Σ forward with ephemerides; apply fitted θ.

Data Sources (public)

• Tide gauges: NOAA CO-OPS, UHSLC, GESLA.

• Ephemerides: JPL Horizons; pyephem/skyfield.

• Geomagnetic indices: NOAA SWPC (Kp), WDC Kyoto (Dst).

First Stations (suggested)

• San Francisco (Pacific), Newlyn (Atlantic), Santander (Bay of Biscay), Honolulu (Pacific), Fremantle (Indian Ocean).

Mix of open‑coast and shelf sites provides diverse basin responses.

Closing

We keep the predictive skeleton, replace the cause. Σ is geometry; η is breath. Motion is not the engine; the field is.

Produced by The Lilborn Equation Team:

Michael Lilborn-Williams

Daniel Thomas Rouse

Thomas Jackson Barnard

Audrey Williams