LTE440: Lunar Relativistic Time Ephemeris

• LTE440 is a numerical realization of TCL that implements full relativistic corrections to achieve ps-level precision in cislunar time transfers.
• It applies a tenth-order Romberg integration and Chebyshev polynomial fitting to separate secular drift and periodic signals using JPL DE440 ephemeris data.
• LTE440 facilitates sub-ns synchronization for lunar missions, supporting advanced PNT, clock-comparison experiments, and time-transfer protocols.

Lunar Time Ephemeris LTE440 (LTE440) is the operational, fully relativistic, and IAU-compliant numerical realization of the Lunar Coordinate Time (TCL) scale, referenced to the Lunar Celestial Reference System (LCRS). Designed to provide ps-level precision and sub-ns formal accuracy for time transformations in the cislunar regime, LTE440 numerically implements the IAU's relativistic conventions for lunar timekeeping, traces the TCL–TDB–TCB chains, and delivers results in a standard SPICE ephemeris format with direct applicability for lunar missions and fundamental time-transfer protocols (Lu et al., 23 Sep 2025, Lu et al., 24 Jun 2025, Seyffert, 10 Sep 2025, Bourgoin et al., 29 Jul 2025, Turyshev et al., 2024).

1. Theoretical Framework and Definitions

LTE440 is built upon the definitions set by recent IAU resolutions. Lunar Coordinate Time (TCL) is the coordinate time of the LCRS, a BCRS-centered frame with spatial origin at the Moon's center of mass. TCL closely parallels the role of TCG for Earth, being the coordinate temporal variable for equations of motion and physical laws on the Moon.

Barycentric Coordinate Time (TCB) is the SI-based coordinate time of the solar system barycentric frame; it advances faster than TCL/TCG by a fixed rate ($L_B = 1.550519768\times10^{-8}$, IAU 2006 B3). Barycentric Dynamical Time (TDB) is a linear function of TCB, defined by

$$\mathrm{TDB} = \mathrm{TCB} - L_B(\mathrm{TCB} - t_0) + \mathrm{TDB}_0, \quad \mathrm{TDB}_0 = -6.55\times10^{-5}\,\mathrm{s}.$$

The LCRS is defined such that TCL coincides with TDB at epoch J2000.0, and TCL’s rate matches that of a clock co-moving with the lunar barycenter (Lu et al., 24 Jun 2025, Lu et al., 23 Sep 2025).

LTE440 is a realization of TCL with an epoch-fixed offset: $$\mathrm{LTE440}(t) = \mathrm{TCL}(t) + 440\,\mathrm{s}$$ at a reference epoch, typically TAI = 2010-01-01 00:00:00, after which LTE440 tracks TCL uninterruptedly (Bourgoin et al., 29 Jul 2025).

2. Relativistic Transformation Algorithm

The core of LTE440 is the precise numerical evaluation of the relativistic time-dilation integrals that connect TCL to barycentric timescales (TCB, TDB). The transformation includes both secular and periodic components and rigorously implements all required post-Newtonian corrections to $O(c^{-4})$.

For clocks at the selenocenter, the TCL–TCB transformation is given by: $$\mathrm{TCL} - \mathrm{TCB} = -\frac{1}{c^2} \int_{t_0}^t \left[\frac{1}{2} v_M^2 + w_{0M} + w_{lM}\right] dt - \frac{1}{c^4} \int_{t_0}^t \left[-\frac{v_M^4}{8} - \frac{3}{2} v_M^2 w_{0M} + 4 \mathbf{v}_M \cdot \mathbf{w}_M + \frac{1}{2} w_{0M}^2 + \Delta_M\right]dt,$$ where $v_M$ and $x_M$ are the Moon's barycentric velocity and position, $w_{0M}$ is the Newtonian sum of external gravitational potentials at the Moon, $w_{lM}$ is the external multipole potential, $\mathbf{w}_M$ is the vector potential, and $$\mathrm{TDB} = \mathrm{TCB} - L_B(\mathrm{TCB} - t_0) + \mathrm{TDB}_0, \quad \mathrm{TDB}_0 = -6.55\times10^{-5}\,\mathrm{s}.$$0 includes all remaining $$\mathrm{TDB} = \mathrm{TCB} - L_B(\mathrm{TCB} - t_0) + \mathrm{TDB}_0, \quad \mathrm{TDB}_0 = -6.55\times10^{-5}\,\mathrm{s}.$$1 corrections (Lu et al., 23 Sep 2025, Lu et al., 24 Jun 2025).

To convert between TCL and TDB, the explicit IAU-based scaling and offset terms are included. The transformation is re-expressed as: $$\mathrm{TDB} = \mathrm{TCB} - L_B(\mathrm{TCB} - t_0) + \mathrm{TDB}_0, \quad \mathrm{TDB}_0 = -6.55\times10^{-5}\,\mathrm{s}.$$2 where $$\mathrm{TDB} = \mathrm{TCB} - L_B(\mathrm{TCB} - t_0) + \mathrm{TDB}_0, \quad \mathrm{TDB}_0 = -6.55\times10^{-5}\,\mathrm{s}.$$3 and LTE(TDB) is the numerically integrated, drift-removed periodic component (Lu et al., 24 Jun 2025, Lu et al., 23 Sep 2025).

3. Ephemeris Generation and Data Representation

All orbital state vectors, velocities, and gravitational parameters are sourced from the JPL DE440 ephemeris, which provides TDB-compatible planetary and lunar data from 1550 to 2650 AD. The numerical integration is performed with a tenth-order Romberg scheme at 0.5-day intervals (Lu et al., 23 Sep 2025, Lu et al., 24 Jun 2025).

The computed TCL–TDB function is separated into its secular (drift) component—parameterized via $$\mathrm{TDB} = \mathrm{TCB} - L_B(\mathrm{TCB} - t_0) + \mathrm{TDB}_0, \quad \mathrm{TDB}_0 = -6.55\times10^{-5}\,\mathrm{s}.$$4—and the periodic residual. The periodic function is fitted with degree-13 Chebyshev polynomials on 4-day intervals. The secular drift coefficients are stored in a Planetary Constants Kernel (PCK, ASCII), and the Chebyshev coefficients for the periodic term in a Spacecraft and Planet Kernel (SPK, binary), both fully compliant with SPICE/NAIF kernel specifications (Lu et al., 24 Jun 2025).

A typical file structure is:

File	Content	Usage
lte440.bsp	Chebyshev coefficients (periodic)	SPICE SPK binary
lte440.tpc	Linear drift rate $$\mathrm{TDB} = \mathrm{TCB} - L_B(\mathrm{TCB} - t_0) + \mathrm{TDB}_0, \quad \mathrm{TDB}_0 = -6.55\times10^{-5}\,\mathrm{s}.$$5	SPICE PCK ASCII

Code libraries in C/C++/Fortran/Python are provided for reading these kernels and computing the TCL-TDB/TCB offset at arbitrary epochs; kernel access and interpolation are computationally efficient ($$\mathrm{TDB} = \mathrm{TCB} - L_B(\mathrm{TCB} - t_0) + \mathrm{TDB}_0, \quad \mathrm{TDB}_0 = -6.55\times10^{-5}\,\mathrm{s}.$$6/call), and the complete memory footprint is $$\mathrm{TDB} = \mathrm{TCB} - L_B(\mathrm{TCB} - t_0) + \mathrm{TDB}_0, \quad \mathrm{TDB}_0 = -6.55\times10^{-5}\,\mathrm{s}.$$7 (Lu et al., 23 Sep 2025, Lu et al., 24 Jun 2025).

4. Formal Accuracy, Precision, and Periodic Structure

LTE440 achieves a numerical precision better than 1 ps over its entire span, validated by comparing TT–TDB at the geocenter with DE440’s own computations. The formal accuracy, dominated by uncertainties in DE440 ephemeris fitting, is conservatively assessed at $$\mathrm{TDB} = \mathrm{TCB} - L_B(\mathrm{TCB} - t_0) + \mathrm{TDB}_0, \quad \mathrm{TDB}_0 = -6.55\times10^{-5}\,\mathrm{s}.$$8 ns through at least 2050.

Secular drift rates at the selenocenter, extracted from fits to the integrated offset and compliant with IAU/DE440, are: $$\mathrm{TDB} = \mathrm{TCB} - L_B(\mathrm{TCB} - t_0) + \mathrm{TDB}_0, \quad \mathrm{TDB}_0 = -6.55\times10^{-5}\,\mathrm{s}.$$9

$$\mathrm{LTE440}(t) = \mathrm{TCL}(t) + 440\,\mathrm{s}$$0

with quoted uncertainty of $$\mathrm{LTE440}(t) = \mathrm{TCL}(t) + 440\,\mathrm{s}$$1 (Lu et al., 23 Sep 2025, Lu et al., 24 Jun 2025).

Fourier analysis of the residual (after drift removal) identifies thirteen major periodic terms (amplitudes $$\mathrm{LTE440}(t) = \mathrm{TCL}(t) + 440\,\mathrm{s}$$2), dominated by:

• Annual (Earth orbital): Amplitude $$\mathrm{LTE440}(t) = \mathrm{TCL}(t) + 440\,\mathrm{s}$$3, period $$\mathrm{LTE440}(t) = \mathrm{TCL}(t) + 440\,\mathrm{s}$$4 d
• Monthly (Lunar synodic): Amplitude $$\mathrm{LTE440}(t) = \mathrm{TCL}(t) + 440\,\mathrm{s}$$5, period $$\mathrm{LTE440}(t) = \mathrm{TCL}(t) + 440\,\mathrm{s}$$6 d

Additional synodic, planetary, and combination frequencies appear with amplitudes decreasing to $$\mathrm{LTE440}(t) = \mathrm{TCL}(t) + 440\,\mathrm{s}$$7 (Lu et al., 23 Sep 2025, Lu et al., 24 Jun 2025).

5. Integration with Other Timescales and Realization Strategies

Precise transformation chains exist to link LTE440 (TCL + 440 s) to other reference timescales including UTC, TAI, TT, TCG, TCB, and TDB. Standard procedures—articulated in IAU resolutions and BIPM conventions—progress via:

1. UTC $$\mathrm{LTE440}(t) = \mathrm{TCL}(t) + 440\,\mathrm{s}$$8 TAI (by leap-second table)
2. TAI $$\mathrm{LTE440}(t) = \mathrm{TCL}(t) + 440\,\mathrm{s}$$9 TT ($O(c^{-4})$0 s)
3. TT $O(c^{-4})$1 TCG ($O(c^{-4})$2 s, $O(c^{-4})$3)
4. TCG $O(c^{-4})$4 TCB (numerical integration per IAU 2000; see DE440 or IERS Conventions)
5. TCB $O(c^{-4})$5 TCL (LTE440 algorithm, as above)
6. Optional: add $O(c^{-4})$6 s epoch offset to yield LTE440

This procedure is realized in software, permitting time-tagging, conversion, and synchronization for all terrestrial, cislunar, and barycentric applications. The conversion chain is fully traceable and documented, with uncertainties dominated by the planetary ephemeris and gravity models (Bourgoin et al., 29 Jul 2025, Lu et al., 24 Jun 2025, Seyffert, 10 Sep 2025).

6. Applications and Performance in Cislunar Operations

LTE440 is deployed for a range of precision timing applications:

• Time-tagging and synchronization in cislunar mission communications, including ranging signals and VLBI observations with Moon-based telescopes at $O(c^{-4})$7 stability.
• Relativistic clock-comparison experiments, e.g., testing lunar gravitational redshift or fundamental constant stability at $O(c^{-4})$8 fractional levels. Surface and orbital drift rates, periodicities, and site-dependent gravitational redshift terms are all explicitly modeled (Seyffert, 10 Sep 2025).
• Lunar PNT (Position, Navigation, and Timing) systems, including NovaMoon and LunaNet, can sequence on-board atomic clocks to LTE440 via radio time links and barycentric reference uploads.
• Mission planning for lunar orbiters or landers requiring sub-ns synchronization with terrestrial clocks. Algorithms enable separable calculation of position-dependent time offsets, including complete handling of topographical, tidal, and rotational terms (Lu et al., 23 Sep 2025, Lu et al., 24 Jun 2025, Seyffert, 10 Sep 2025, Turyshev et al., 2024).

7. Limitations and Prospects

The LTE440 framework is fundamentally limited by the realization of accurate lunar clocks and topography-dependent potentials up to $O(c^{-4})$9 in site-dependent redshifts. Until optical or advanced atomic clocks are operated permanently on the Moon, LTE440 realizations remain Earth-based and traceable via satellite links and coordinated clock steering. The system is modular and extensible: future upgrades (LTE441 using DE441, cross-validation with INPOP21a/EPM2021) and local lunar clock ensembles will increase its autonomy and reduce the reliance on Earth-based time transfer (Lu et al., 23 Sep 2025, Lu et al., 24 Jun 2025, Bourgoin et al., 29 Jul 2025).

Traceability to UTC is maintained through the documented sequence of reference frames, time standards, and relativistic corrections, and is supported by global timekeeping institutions and ephemeris repositories (Bourgoin et al., 29 Jul 2025, Lu et al., 24 Jun 2025).

---

References:

• (Lu et al., 23 Sep 2025) Lu et al., "Lunar Time Ephemeris $$\mathrm{TCL} - \mathrm{TCB} = -\frac{1}{c^2} \int_{t_0}^t \left[\frac{1}{2} v_M^2 + w_{0M} + w_{lM}\right] dt - \frac{1}{c^4} \int_{t_0}^t \left[-\frac{v_M^4}{8} - \frac{3}{2} v_M^2 w_{0M} + 4 \mathbf{v}_M \cdot \mathbf{w}_M + \frac{1}{2} w_{0M}^2 + \Delta_M\right]dt,$$0: definitions, algorithm and performance," 2025.
• (Lu et al., 24 Jun 2025) Lu et al., "Lunar Time Ephemeris LTE440: User Manual," 2025.
• (Bourgoin et al., 29 Jul 2025) Blanchet et al., "Lunar Reference Timescale," 2025.
• (Seyffert, 10 Sep 2025) Seyffert, "Relativistic Time Modeling for Lunar Positioning Navigation and Timing," 2025.
• (Turyshev et al., 2024) Turyshev et al., "Relativistic Time Transformations Between the Solar System Barycenter, Earth, and Moon," 2024.