@essrt/physics 22.1.0 → 23.2.0

package/types/constants/coefficients.d.ts

@@ -2,7 +2,7 @@
 // The fitted coefficients' shapes for TypeScript consumers (§2g); values live
 // in coefficients.js, emitted VERBATIM from the fit artifacts under data/.
 
-export declare const COEFFICIENTS_HASH: "91e2d041c7b6e038";
+export declare const COEFFICIENTS_HASH: "ff77fbc9baa370e4";
 
 export declare const FITTED_COEFFICIENTS: {
   readonly BALANCE_PRESETS: Array<Array<number | string>>;
@@ -275,70 +275,70 @@ export declare const FITTED_COEFFICIENTS: {
             "ratio": number;
             "label": string;
           };
-          "27.2691": {
+          "27.2689": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "28.8426": {
+          "28.8424": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "29.0050": {
+          "29.0048": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "29.6174": {
+          "29.6172": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "29.8094": {
+          "29.8092": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "39.5564": {
+          "39.5560": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "40.2297": {
+          "40.2293": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "40.9244": {
+          "40.9239": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "41.6450": {
+          "41.6446": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;
             "ratio": number;
             "label": string;
           };
-          "53.6182": {
+          "53.6175": {
             "period_kyr": number;
             "lr04_post_mpt_amp": number;
             "epica_amp": number;

package/types/constants/generated.d.ts

@@ -2,14 +2,14 @@
 // Gives the TypeScript website full type safety at the boundary (§2g) while
 // packages/physics stays JavaScript.
 
-export declare const CONSTANTS_HASH: "f4f43386df2781b8";
+export declare const CONSTANTS_HASH: "896558c2d878571d";
 
 export declare const MODEL_VERSION: string;
 
 export declare const PREPRINT_DOI: string;
 
 export declare const DEFAULT_CONSTANTS: {
-  readonly hash: "f4f43386df2781b8";
+  readonly hash: "896558c2d878571d";
   readonly additionalBodiesReference: {
     "pluto": {
       "solarYearInput": number;
@@ -83,24 +83,11 @@ export declare const DEFAULT_CONSTANTS: {
   };
   readonly deepTime: {
     "alpha1PerMa": number;
-    "alpha3PerMa3": number;
-    "alpha4PerMa4": number;
     "alphaGiaRateJ2000PerYr": number;
     "alphaGiaRelaxationKyr": number;
     "j2RateJ2000PerYr": number;
     "dtStackTaperFullHalfwidthYr": number;
     "dtStackTaperTotalHalfwidthYr": number;
-    "recessionRegime": {
-      "jointMa": number;
-      "knotAgesMa": number[];
-      "knotDistancesKm": number[];
-      "genesisMa": number;
-      "rocheLimitKm": number;
-      "solarOceanLeakBeta0": number;
-      "thermalPumpStartMa": number;
-      "thermalPumpEndMa": number;
-      "thermalPumpFactor": number;
-    };
   };
   readonly earthOrbital: {
     "j2000EpochYear": number;
@@ -207,6 +194,8 @@ export declare const DEFAULT_CONSTANTS: {
     "earthJ2": number;
     "earthEquatorialRadiusKm": number;
     "earthParallaxRadiusKm": number;
+    "earthMeanRadiusKm": number;
+    "solarConstantWm2": number;
     "solarParallaxArcsec": number;
     "arcsecDisplacementKm": number;
   };
@@ -566,7 +555,6 @@ export declare const REFERENCE_DATA: {
     "wellsRecessionCmYr": number;
     "modernLLRRecessionCmYr": number;
     "tidalLockApproachGyr": number;
-    "essrtFormulaHorizonGyr": number;
     "iceAlbedoShareLongPct": number;
     "testAObliquityLagPercentile": number;
     "testAObliquityPeakMatch": string;

package/types/constants/units.d.cts

@@ -0,0 +1,8 @@
+/** The SI day, seconds. */
+export const SECONDS_PER_DAY: 86400;
+/** The Julian year, days (IAU). */
+export const DAYS_PER_JULIAN_YEAR: 365.25;
+/** The Julian century, days. */
+export const JULIAN_CENTURY_DAYS: 36525;
+/** The Julian year, SI seconds. */
+export const SECONDS_PER_JULIAN_YEAR: number;

package/types/deltat/deep-time-chain.d.ts

@@ -0,0 +1,266 @@
+/**
+ * The ΔT stack's switches. `resonator` gates the swing's δLOD terms; the
+ * swing's ΔT term is gated by `resonatorDeltaT` when given (the Node twin
+ * gates it on BOTH env vars, DT_CORRECTIONS and DT_RESONATOR), else by
+ * `resonator`. The four cycles' δLOD terms follow their flags only when
+ * `cycleLodFollowsFlags` is true (the browser's research toggles switch a
+ * cycle's ΔT and δLOD together); the Node twin sums them unconditionally.
+ * Their ΔT terms always follow their flags. Identical with every flag ON.
+ * @typedef {{bond: boolean, hallstatt: boolean, jose5: boolean, jose4: boolean, resonator: boolean, resonatorDeltaT?: boolean, cycleLodFollowsFlags?: boolean}} DeltaTFlags
+ */
+/**
+ * @param {{
+ *   constants: any,
+ *   fitted: any,
+ *   deepModes: any,
+ *   moonDistanceArtifact: any,
+ *   ageAnchorYear?: number,
+ *   flagsAt?: () => DeltaTFlags,
+ *   planetAuRatios?: Record<string, number>,
+ *   stellarDayRaProjection?: number,
+ * }} deps — `constants` the generated DEFAULT_CONSTANTS, `fitted` the
+ *   generated FITTED_COEFFICIENTS, `deepModes` the DEEP_MODES_ARTIFACT,
+ *   `moonDistanceArtifact` the MOON_DISTANCE_ARTIFACT — the cited distance
+ *   AND mean-obliquity columns, plan 12 option A / §8 (a) (injected so a
+ *   counterfactual context can replace any of them); `ageAnchorYear` the
+ *   scene epoch (default constants.timeReference.sceneEpochYear);
+ *   `flagsAt` the ΔT stack's switches (default all ON); `planetAuRatios`
+ *   the planets' J2000 semi-major axes in AU (for the per-planet family);
+ *   `stellarDayRaProjection` cos ε₀-class factor for the projected stellar day.
+ */
+export function createDeepTimeChain(deps: {
+    constants: any;
+    fitted: any;
+    deepModes: any;
+    moonDistanceArtifact: any;
+    ageAnchorYear?: number;
+    flagsAt?: () => DeltaTFlags;
+    planetAuRatios?: Record<string, number>;
+    stellarDayRaProjection?: number;
+}): Readonly<{
+    meanSolarYearDays: number;
+    meanSiderealYearDays: any;
+    meanSiderealYearSeconds: number;
+    meanLengthOfDay: any;
+    meanTropicalYearJ2000Seconds: number;
+    siTropicalYearDays: number;
+    axialPrecessionJ2000Years: number;
+    axialPrecessionJ2000Days: number;
+    G_CONSTANT: any;
+    MASS_RATIO_EARTH_MOON: any;
+    currentAUDistance: any;
+    earthMoiFactorJ2000: any;
+    moonDistanceKm: any;
+    moonSiderealMonthInput: any;
+    moonOrbitalShift: number;
+    moonDistanceCorrected: number;
+    GM_EARTH_MOON_SYSTEM: number;
+    GM_EARTH_ALONE: number;
+    GM_MOON_ALONE: number;
+    GM_SUN_PLUS_EARTH: number;
+    GM_SUN: number;
+    M_SUN: number;
+    moonApsidalJ2000Seconds: number;
+    moonNodalJ2000Seconds: number;
+    moonSiderealMonthJ2000Seconds: number;
+    moonPeriods: {
+        moonApsidalPrecessionDaysEarth: number;
+        moonNodalPrecessionDaysEarth: number;
+        moonApsidalJ2000Seconds: number;
+        moonNodalJ2000Seconds: number;
+        moonSiderealMonthJ2000Seconds: number;
+        moonAnomalisticMonthDays: number;
+        moonNodalMonthDays: number;
+        moonSynodicMonthDays: number;
+        moonTropicalMonthDays: number;
+        moonFullMoonCycleEarthDays: number;
+        moonFullMoonCycleICRFDays: number;
+        moonApsidalMeetsNodalDays: number;
+        moonLunarLevelingCycleDays: number;
+        moonDraconicYearICRFDays: number;
+        moonDraconicYearEarthDays: number;
+    };
+    certifiedAxialPrecessionJ2000Years: () => number;
+    precessionRateJ2000ArcsecPerYr: () => number;
+    EPOCH_PARAMS: import("../layer0/index.js").EpochParams;
+    ageAnchorYear: any;
+    alphaGia: {
+        alphaAt: (year: number) => number;
+        j2RatioAt: (year: number) => number;
+        laggedL1PermilAt: (year: number) => number;
+        kPerPermille: number;
+        relaxationKyr: number;
+        laggedL1Terms: Array<{
+            period_kyr: number;
+            a: number;
+            b: number;
+        }>;
+    };
+    earthMoiFactorAtAge: (tMa: number) => number;
+    withLatticeAlpha: (build: () => void) => void;
+    j2RatioAtYear: (year: number) => number;
+    iEarthAtAge: (tMa: number) => number;
+    recession: {
+        distanceMetresAtAge: (tMa: number) => number | null;
+        slopeMetresPerMaAtAge: (tMa: number) => number | null;
+        ratioAtAge: (tMa: number) => number | null;
+        spanMa: number;
+    };
+    meanObliquity: {
+        meanObliquityDegAtAge: (tMa: number) => number | null;
+        cosRatioAtAge: (tMa: number) => number | null;
+        meanObliquityJ2000Deg: number;
+        spanMa: number;
+    };
+    dtCycles: {
+        taperAt: (year: number) => number;
+        taperDerivativeAt: (year: number) => number;
+        cycleDeltaTSecondsAt: (cycleKey: string, year: number) => number;
+        cycleLodSecondsAt: (cycleKey: string, year: number) => number;
+        swingDeltaTSecondsAt: (year: number) => number;
+        swingLodSecondsAt: (year: number) => number;
+        swingLodRateAt: (year: number) => number;
+    };
+    s3ArcsecPerYr: number;
+    deepLod: {
+        lodSecondsAtAge: (t_Ma: number) => number | null;
+        lodSecondsAtAgeWithAlpha: (t_Ma: number, alpha: number) => number | null;
+        lodHoursAtAge: (t_Ma: number) => number | null;
+        lunisolarPrecessionRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
+        lunisolarPrecessionPeriodYearsAtAge: (t_Ma: number) => number | null;
+        lunisolarTorqueRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
+        lunisolarTorquePeriodYearsAtAge: (t_Ma: number) => number | null;
+        precessionObliquityFactorAtAge: (t_Ma: number) => number | null;
+        lunarTorqueFactorAtAge: (tMa: number) => (number | null);
+        precessionTorqueTermAtAge: (tMa: number) => (number | null);
+        siderealYearSecondsAtAge: (t_Ma: number) => number;
+        tropicalYearSecondsAtAge: (t_Ma: number) => number;
+        tropicalYearDaysAtAge: (t_Ma: number) => number | null;
+        tropicalYearSecondsFromLod: (t_Ma: number, LOD_s: number | null) => number;
+        eclipticLodCorrectionSecondsAtAge: (t_Ma: number) => number | null;
+        yearInDaysAtAge: (t_Ma: number) => number | null;
+        deltaTRawSecondsAtAge: (t_Ma: number) => number;
+        lodSecondsWithCorrectionsAtAge: (t_Ma: number) => number | null;
+        dLodDtDecompositionAtAge: (t_Ma: number) => {
+            tidal: number | null;
+            gia: number | null;
+            stack: number | null;
+            resonator: number | null;
+            net_L2: number | null;
+            net_L3: number | null;
+            net_L4: number | null;
+        };
+    };
+    l0: import("../layer0/index.js").EpochPrimitives;
+    deepEcc: {
+        eccAt: (tYr: number) => number;
+        eccRateAt: (tYr: number) => number;
+        modulation: (tMa: number, s: number) => number;
+        channelIntegral: (T: number, s: number) => number;
+        eFactorAt: (tYr: number) => number;
+        e0: number;
+        g0: number;
+    };
+    eCompModulation: (tMa: number, s: number) => number;
+    moonChain: {
+        distanceMetresAtAge: (tMa: number) => number | null;
+        distanceKmAtAge: (tMa: number) => number | null;
+        solarDeltaAKmAtAge: (tMa: number, aApparentKm: number) => number | null;
+        distanceCorrectedKmAtAge: (tMa: number) => number | null;
+        siderealMonthSecondsAtAge: (tMa: number) => number | null;
+        synodicMonthSecondsAtAge: (tMa: number) => number | null;
+        tropicalMonthSecondsAtAge: (tMa: number) => number | null;
+        perigeePrecessionSecondsAtAge: (tMa: number) => number | null;
+        nodePrecessionSecondsAtAge: (tMa: number) => number | null;
+        anomalisticMonthSecondsAtAge: (tMa: number) => number | null;
+        nodalMonthSecondsAtAge: (tMa: number) => number | null;
+        apsidalMeetsNodalSecondsAtAge: (tMa: number) => number | null;
+        lunarLevelingSecondsAtAge: (tMa: number) => number | null;
+    };
+    meanMoonDistanceMetresAtAge: (tMa: number) => number | null;
+    meanMoonDistanceAtAge: (tMa: number) => number | null;
+    meanLodSecondsAtAge: (tMa: number) => number | null;
+    meanLodHoursAtAge: (t_Ma: number) => number | null;
+    meanLunisolarPrecessionRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
+    meanLunisolarPrecessionPeriodYearsAtAge: (tMa: number) => number | null;
+    meanLunisolarTorqueRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
+    meanLunisolarTorquePeriodYearsAtAge: (t_Ma: number) => number | null;
+    precessionObliquityFactorAtAge: (t_Ma: number) => number | null;
+    meanObliquityCitedDegAtAge: (tMa: number) => number | null;
+    meanAuAtAge: (tMa: number) => number;
+    meanSiderealYearSecondsAtAge: (tMa: number) => number;
+    meanTropicalYearSecondsAtAge: (t_Ma: number) => number;
+    meanTropicalYearDaysAtAge: (t_Ma: number) => number | null;
+    meanYearInDaysAtAge: (t_Ma: number) => number | null;
+    meanDeltaTSecondsAtAge: (tMa: number) => number;
+    bondCycleDeltaTCorrection: (year: number) => number;
+    hallstattCycleDeltaTCorrection: (year: number) => number;
+    jose5CycleDeltaTCorrection: (year: number) => number;
+    jose4CycleDeltaTCorrection: (year: number) => number;
+    bondCycleLodCorrection: (year: number) => number;
+    hallstattCycleLodCorrection: (year: number) => number;
+    jose5CycleLodCorrection: (year: number) => number;
+    jose4CycleLodCorrection: (year: number) => number;
+    dtCycleLodCorrectionSum: (year: number) => number;
+    resonatorSwingDeltaTCorrection: (year: number) => number;
+    resonatorSwingLodCorrection: (year: number) => number;
+    resonatorSwingLodRate: (year: number) => number;
+    meanLodSecondsWithCorrectionsAtAge: (t_Ma: number) => number | null;
+    dLodDtDecompositionAtAge: (t_Ma: number) => {
+        tidal: number | null;
+        gia: number | null;
+        stack: number | null;
+        resonator: number | null;
+        net_L2: number | null;
+        net_L3: number | null;
+        net_L4: number | null;
+    };
+    h5Correction: (year: number) => number;
+    computeLodKinematicSecondsAtEpoch: (year: number) => number;
+    computeLodRealSecondsAtEpoch: (year: number) => number;
+    meanSolarDeltaAAtAge: (tMa: number, aApparentKm: number) => number | null;
+    meanMoonDistanceCorrectedAtAge: (tMa: number) => number | null;
+    meanMoonSiderealMonthAtAge: (tMa: number) => number | null;
+    meanSynodicMonthAtAge: (tMa: number) => number | null;
+    meanTropicalMonthAtAge: (tMa: number) => number | null;
+    meanLunarPerigeePrecessionAtAge: (tMa: number) => number | null;
+    meanLunarNodePrecessionAtAge: (tMa: number) => number | null;
+    meanAnomalisticMonthAtAge: (tMa: number) => number | null;
+    meanNodalMonthAtAge: (tMa: number) => number | null;
+    meanApsidalMeetsNodalAtAge: (tMa: number) => number | null;
+    meanLunarLevelingCycleAtAge: (tMa: number) => number | null;
+    meanStellarDayAtAge: (tMa: number) => number | null;
+    meanStellarDayProjectedAtAge: (tMa: number) => number | null;
+    meanSiderealDayAtAge: (tMa: number) => number | null;
+    meanPlanetSemiMajorAxisAtAge: (planetName: string, tMa: number) => number | null;
+    /** @param {number} tMa @param {number} periodJ2000Seconds @returns {number} */
+    meanPlanetOrbitalPeriodAtAge: (tMa: number, periodJ2000Seconds: number) => number;
+    /** @param {number} tMa @param {number} periodJ2000Seconds @returns {number} */
+    meanEarthPlanetSynodicPeriodAtAge: (tMa: number, periodJ2000Seconds: number) => number;
+}>;
+/**
+ * Moon-channel eccentricity sensitivities (perigee/node), the [g/g₀]^s
+ * exponents of the factored deep-time law — Meeus-effective (S_N moved
+ * 1.0 → 1.018 with the v4 frame-attribution batch). ONE home since B2 (they
+ * were literals in model.js, tools/lib/deep-time.js and script.js `_FW_MOON`).
+ */
+export const MOON_ECC_SENSITIVITY_PERIGEE: 2.407;
+export const MOON_ECC_SENSITIVITY_NODE: 1.018;
+/**
+ * The ΔT stack's switches. `resonator` gates the swing's δLOD terms; the
+ * swing's ΔT term is gated by `resonatorDeltaT` when given (the Node twin
+ * gates it on BOTH env vars, DT_CORRECTIONS and DT_RESONATOR), else by
+ * `resonator`. The four cycles' δLOD terms follow their flags only when
+ * `cycleLodFollowsFlags` is true (the browser's research toggles switch a
+ * cycle's ΔT and δLOD together); the Node twin sums them unconditionally.
+ * Their ΔT terms always follow their flags. Identical with every flag ON.
+ */
+export type DeltaTFlags = {
+    bond: boolean;
+    hallstatt: boolean;
+    jose5: boolean;
+    jose4: boolean;
+    resonator: boolean;
+    resonatorDeltaT?: boolean;
+    cycleLodFollowsFlags?: boolean;
+};

package/types/deltat/deep-time.d.cts

@@ -22,9 +22,10 @@ export type DeepTimeLodConstants = {
      * - Farhat polynomial anchor
      */
     aMoonNowMetres: number;
+    /**
+     * - the LLR recession rate (the J2000 slope of the cited history)
+     */
     alpha1PerMa: number;
-    alpha3PerMa3: number;
-    alpha4PerMa4: number;
     /**
      * - ω₀, the J2000 day length (the SI day since plan 09; identifier kept)
      */
@@ -39,7 +40,10 @@ export type DeepTimeLodConstants = {
 };
 export type DeepTimeLodDeps = {
     constants: DeepTimeLodConstants;
-    moonDistanceMetresAtAge: (tMa: number) => number;
+    /**
+     * - null beyond the cited history's span
+     */
+    moonDistanceMetresAtAge: (tMa: number) => number | null;
     /**
      * - GIA α(t)
      */
@@ -77,11 +81,26 @@ export type DeepTimeLodDeps = {
     nodalPeriodYearsFn: () => number;
     /**
      * - OPTIONAL time-dependent
-     * Earth-Moon angular momentum (the Driver-1½ solar channels,
-     * recession-history.cjs). Absent → the J2000 constant, which the budget
-     * module equals exactly for t ≤ jointMa — the pure-twin convention.
+     * Earth-Moon angular momentum. Absent → the J2000 constant: conservation
+     * (plan 12 option A — the Driver-1½ solar-channel budget is deleted;
+     * measured: conservation on the cited history reproduces Farhat's own LOD
+     * to 0.01 h at 100 Ma, 0.14 h at 2.46 Ga).
      */
     lEmAtAgeKgm2S?: ((tMa: number) => number) | undefined;
+    /**
+     * -
+     * OPTIONAL analytic slope of the distance history (metres per Ma, positive
+     * into the past means a growing a... sign as da/dt_Ma); absent → a central
+     * difference on moonDistanceMetresAtAge (the LOD-rate decomposition's tidal term).
+     */
+    moonDistanceSlopeMetresPerMaAtAge?: ((tMa: number) => number | null) | undefined;
+    /**
+     * -
+     * OPTIONAL cos ε̄(t)/cos ε̄(0) of the cited mean obliquity (earth/mean-
+     * obliquity-history) — the factor the precession RATE of date carries on the
+     * constant (plan 12 §8 option (a)); absent → 1, the constant alone, bit-identical.
+     */
+    meanObliquityCosRatioAtAge?: ((tMa: number) => number | null) | undefined;
 };
 /**
  * @typedef {Object} DeepTimeLodConstants
@@ -93,8 +112,7 @@ export type DeepTimeLodDeps = {
  * @property {number} eFactorMoon - √(1−e²) factor of the lunar orbit
  * @property {number} aLockMetres - tidal-lock distance (LOD denominator → 0)
  * @property {number} aMoonNowMetres - Farhat polynomial anchor
- * @property {number} alpha1PerMa @property {number} alpha3PerMa3
- * @property {number} alpha4PerMa4
+ * @property {number} alpha1PerMa - the LLR recession rate (the J2000 slope of the cited history)
  * @property {number} lodNowH13Seconds - ω₀, the J2000 day length (the SI day since plan 09; identifier kept)
  * @property {number} meanSiderealYearJ2000Seconds
  * @property {number} solarMassLossFracPerYear
@@ -104,7 +122,7 @@ export type DeepTimeLodDeps = {
 /**
  * @typedef {Object} DeepTimeLodDeps
  * @property {DeepTimeLodConstants} constants
- * @property {(tMa: number) => number} moonDistanceMetresAtAge
+ * @property {(tMa: number) => number|null} moonDistanceMetresAtAge - null beyond the cited history's span
  * @property {(tMa: number) => number} moiFactorAtAge - GIA α(t)
  * @property {(year: number) => number} cycleLodSumAt - gated δLOD sum (incl. swing)
  * @property {(year: number) => number} swingLodAt - gated swing δLOD alone
@@ -123,9 +141,18 @@ export type DeepTimeLodDeps = {
  *   formerly the spin unit's H/5, 67,063 yr, 2.5 % off and scaling with the
  *   spin at deep time). Lazy, like the precession anchor.
  * @property {(tMa: number) => number} [lEmAtAgeKgm2S] - OPTIONAL time-dependent
- *   Earth-Moon angular momentum (the Driver-1½ solar channels,
- *   recession-history.cjs). Absent → the J2000 constant, which the budget
- *   module equals exactly for t ≤ jointMa — the pure-twin convention.
+ *   Earth-Moon angular momentum. Absent → the J2000 constant: conservation
+ *   (plan 12 option A — the Driver-1½ solar-channel budget is deleted;
+ *   measured: conservation on the cited history reproduces Farhat's own LOD
+ *   to 0.01 h at 100 Ma, 0.14 h at 2.46 Ga).
+ * @property {(tMa: number) => number|null} [moonDistanceSlopeMetresPerMaAtAge] -
+ *   OPTIONAL analytic slope of the distance history (metres per Ma, positive
+ *   into the past means a growing a... sign as da/dt_Ma); absent → a central
+ *   difference on moonDistanceMetresAtAge (the LOD-rate decomposition's tidal term).
+ * @property {(tMa: number) => number|null} [meanObliquityCosRatioAtAge] -
+ *   OPTIONAL cos ε̄(t)/cos ε̄(0) of the cited mean obliquity (earth/mean-
+ *   obliquity-history) — the factor the precession RATE of date carries on the
+ *   constant (plan 12 §8 option (a)); absent → 1, the constant alone, bit-identical.
  */
 /** @param {DeepTimeLodDeps} deps */
 export function createDeepTimeLod(deps: DeepTimeLodDeps): {
@@ -134,6 +161,9 @@ export function createDeepTimeLod(deps: DeepTimeLodDeps): {
     lodHoursAtAge: (t_Ma: number) => number | null;
     lunisolarPrecessionRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
     lunisolarPrecessionPeriodYearsAtAge: (t_Ma: number) => number | null;
+    lunisolarTorqueRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
+    lunisolarTorquePeriodYearsAtAge: (t_Ma: number) => number | null;
+    precessionObliquityFactorAtAge: (t_Ma: number) => number | null;
     lunarTorqueFactorAtAge: (tMa: number) => (number | null);
     precessionTorqueTermAtAge: (tMa: number) => (number | null);
     siderealYearSecondsAtAge: (t_Ma: number) => number;

package/types/deltat/moon-distance-artifact.d.cts

@@ -0,0 +1,22 @@
+export const MOON_DISTANCE_ARTIFACT: Readonly<{
+    source: {
+        citation: string;
+        doi: string;
+        distributor: string;
+        fileSha256: {
+            data: string;
+            readme: string;
+        };
+    };
+    earthRadiusKm: 6378.1366;
+    ageStepMa: 1;
+    spanMa: 4429;
+    ageMa: number[];
+    aEarthRadii: number[];
+    aMinEarthRadii: number[];
+    aMaxEarthRadii: number[];
+    obliquityDeg: number[];
+    obliquityMinDeg: number[];
+    obliquityMaxDeg: number[];
+}>;
+export const MOON_DISTANCE_ARTIFACT_HASH: "850cd0b266d46a71";

package/types/deltat/moon-distance-history.d.cts

@@ -0,0 +1,36 @@
+export type MoonDistanceTable = {
+    /**
+     * - ascending, 0 first, uniform 1-Myr step
+     */
+    ageMa: ReadonlyArray<number>;
+    /**
+     * - the distance column (any consistent unit; only ratios are read)
+     */
+    aEarthRadii: ReadonlyArray<number>;
+};
+/**
+ * @typedef {Object} MoonDistanceTable
+ * @property {ReadonlyArray<number>} ageMa - ascending, 0 first, uniform 1-Myr step
+ * @property {ReadonlyArray<number>} aEarthRadii - the distance column (any consistent unit; only ratios are read)
+ */
+/**
+ * @param {{ table: MoonDistanceTable, aMoonNowMetres: number, alpha1PerMa: number }} deps
+ *   aMoonNowMetres — the model's J2000 anchor (the observed mean distance);
+ *   alpha1PerMa — the LLR rate as a fraction of a₀ per Ma (negative into the past: a shrinks).
+ * @returns {{
+ *   distanceMetresAtAge: (tMa: number) => number|null,
+ *   slopeMetresPerMaAtAge: (tMa: number) => number|null,
+ *   ratioAtAge: (tMa: number) => number|null,
+ *   spanMa: number,
+ * }}
+ */
+export function createMoonDistanceHistory(deps: {
+    table: MoonDistanceTable;
+    aMoonNowMetres: number;
+    alpha1PerMa: number;
+}): {
+    distanceMetresAtAge: (tMa: number) => number | null;
+    slopeMetresPerMaAtAge: (tMa: number) => number | null;
+    ratioAtAge: (tMa: number) => number | null;
+    spanMa: number;
+};

package/types/earth/deep-orbital-history.d.cts

@@ -25,9 +25,11 @@
  *   measured at the 500-yr artifact cadence: 0.16″/0.18″ rms vs IAU-2006
  *   over 1900–2100/1600–2400 vs the era tier's 0.31″/0.63″, and 0.0396°
  *   vs La2004 over −200 kyr vs the deep tier's 0.069° — the C-1 verdict,
- *   plan 02) and zetaModes serves only as the TAIL beyond the span (the
- *   seam at the span edge is the extraction residual, far outside every
- *   certified window). solarTorqueShareJ2000 and
+ *   plan 02) and zetaModes serves only as the TAIL beyond the span, joined
+ *   to the series at each edge so n̂ is continuous there (plan 12 §8.3 — the
+ *   J2000-point anchor left a 1.6° step at −10 Myr that kicked the hybrid's
+ *   mean obliquity +0.9° for the whole tail; far outside every certified window).
+ *   solarTorqueShareJ2000 and
  *   dynamicalEllipticityRatioAtYearFn: THE PRECESSION OF DATE — the two
  *   factors on the precession constant (see alphaAtGeneral below); absent,
  *   the constant is the secular law alone, bit-identical.