Lunar Celestial Reference System (LCRS)

• LCRS is a quasi-inertial, selenocentric reference system defined by the IAU, supporting high-precision lunar navigation, timing, and scientific research.
• It employs post-Newtonian corrections in its metric tensor and time transformations to achieve millimeter-level positioning and sub-nanosecond synchronization.
• LCRS underpins practical lunar applications by integrating data from LLR, VLBI, and GNSS for accurate cartography and operational coordination.

The Lunar Celestial Reference System (LCRS) is the quasi-inertial, selenocentric spacetime reference framework adopted by the IAU for all high-precision, relativistically consistent lunar and cislunar navigation, timing, cartography, and science. It provides the conceptual and operational infrastructure enabling millimeter-level positioning, sub-nanosecond timing synchronization, and rigorous interoperability between Earth, lunar, and planetary frames.

1. Definition, IAU Foundations, and Key Parameters

The LCRS is a kinematically non-rotating coordinate system whose spatial origin is at the Moon’s center of mass (CoM), and whose axes are maintained (to post-Newtonian order) parallel to those of the Solar System Barycentric Celestial Reference System (BCRS) (Turyshev et al., 15 Nov 2025, Turyshev, 29 Jul 2025, Sośnica et al., 17 Oct 2025, Fienga et al., 2024). The formal coordinates are $(\mathcal{T}, \mathcal{X})$, with time coordinate $\mathcal{T} \equiv$ TCL (Lunicentric Coordinate Time). This time coordinate is synchronized so that TCL equals TCB (Barycentric Coordinate Time) at a specified epoch (e.g., 1977-01-01 T00:00:32.184 at the lunar CoM).

The LCRS metric, at 1PN order, is:

$$ds^2_\text{LCRS} = \left[1 - \frac{2}{c^2}(U_M + U^*_\text{tid})\right] c^2 d\mathcal{T}^2 - \left[1 + \frac{2}{c^2}(U_M + U^*_\text{tid})\right] d\mathcal{X}^2 + \mathcal{O}(c^{-4})$$

where:

• $U_M(\mathcal{T},\mathcal{X}) = GM_M/R$ is the monopole lunar gravitational potential,
• $U^*_{\rm tid}$ is the Newtonian tidal potential from Sun, Earth, and planets, expanded about the lunar center.

Fundamental rate constants for time transformations (cf. (Turyshev et al., 2024), Table I in (Turyshev et al., 15 Nov 2025)):

• $L_L = 3.13905 \times 10^{-11}$ (TCL→TL rate at lunar surface),
• $L_H = 1.48254 \times 10^{-8}$ (secular BCRS→TCL rate),
• $L_M = 1.48568 \times 10^{-8}$ (BCRS→TL).

2. LCRS–Related Time Scales and Transformations

TCL (Lunicentric Coordinate Time) is defined analogously to TCG (Geocentric Coordinate Time), providing a coordinate time for the lunar reference system realized from TCB using:

$$\mathcal{T} = t - \frac{1}{c^2}\left\{ \int^t \left[ \frac{1}{2} v_M^2 + \sum_{B \neq M} \frac{G M_B}{r_{BM}} \right] dt' + v_M \cdot r_M \right\} + \mathcal{O}(10^{-16}(t-t_0))$$

Here $v_M$ and $\mathcal{T} \equiv$0 are the Moon’s barycentric velocity and position.

TL ("lunar proper time," or "lunar surface time") is related by a constant rate offset from TCL: $\mathcal{T} \equiv$1 The drift between TL and Earth TT is $\mathcal{T} \equiv$2, and periodic corrections $\mathcal{T} \equiv$3 (mean anomalistic period) must be included for sub-nanosecond precision (Turyshev et al., 2024, Turyshev, 29 Jul 2025).

The recommended LCRS operational timescale is to set TL = TCL with only the gravitational redshift difference, minimizing complications with coordinate scaling and parameters (Bourgoin et al., 29 Jul 2025).

3. Metric Tensor, Gravitational Potentials, and Reference Ephemerides

Full post-Newtonian (PN) treatment is required for mm-level and sub-nanosecond realization. The LCRS metric retains all contributions above fractional $\mathcal{T} \equiv$4 ($\mathcal{T} \equiv$50.1 ps) by expanding the lunar gravity field to spherical-harmonic degree $\mathcal{T} \equiv$6, with temporal Love number variations, and including external Earth/Sun tidal and inertial multipoles up to $\mathcal{T} \equiv$7 (Turyshev, 29 Jul 2025). For time and frame transformations, use Chebyshev polynomials representing $\mathcal{T} \equiv$8 (precession), $\mathcal{T} \equiv$9 (nutation), and $$ds^2_\text{LCRS} = \left[1 - \frac{2}{c^2}(U_M + U^*_\text{tid})\right] c^2 d\mathcal{T}^2 - \left[1 + \frac{2}{c^2}(U_M + U^*_\text{tid})\right] d\mathcal{X}^2 + \mathcal{O}(c^{-4})$$0 (proper rotation), as extracted from high-precision lunar ephemerides (INPOP21a, DE430, EPM2021).

Typical transformation chain (see (Sośnica et al., 17 Oct 2025), Eq. (13)): $$ds^2_\text{LCRS} = \left[1 - \frac{2}{c^2}(U_M + U^*_\text{tid})\right] c^2 d\mathcal{T}^2 - \left[1 + \frac{2}{c^2}(U_M + U^*_\text{tid})\right] d\mathcal{X}^2 + \mathcal{O}(c^{-4})$$1 where PA is the co-rotating principal-axis lunar frame, and LCRF is the (quasi-)inertial LCRS-aligned frame.

Transformations between BCRS, GCRS, and LCRS coordinate times, and their proper times, incorporate both secular and periodic terms, as detailed in (Turyshev et al., 2024, Turyshev, 29 Jul 2025, Bourgoin et al., 29 Jul 2025).

4. Materialization and Realization: Data, Frames, and NovaMoon

Material realization of the LCRS (and its body-fixed counterpart, the LRS) depends on accurate lunar ephemerides, Lunar Laser Ranging (LLR), Very Long Baseline Interferometry (VLBI), surface retroreflectors, and new multi-technique packages like NovaMoon (Molli et al., 9 Feb 2026). The International Lunar Reference Frame (ILRF) is the physical realization of the rotating body-fixed frame, using principal axes tied to the Moon’s inertia tensor, with its origin at the lunocenter (CoM) (Sośnica et al., 17 Oct 2025, Pavlov, 2019).

Key instrumentation at NovaMoon (ESA Argonaut):

• LLR retroreflector (mm-level ranging),
• VLBI transmitter (sub-mas angular tie to ICRF),
• LunaNet-compliant GNSS receiver (differential corrections in the local Lunar Reference Frame),
• atomic clock ensemble defining TL/TCL (frequency stability $$ds^2_\text{LCRS} = \left[1 - \frac{2}{c^2}(U_M + U^*_\text{tid})\right] c^2 d\mathcal{T}^2 - \left[1 + \frac{2}{c^2}(U_M + U^*_\text{tid})\right] d\mathcal{X}^2 + \mathcal{O}(c^{-4})$$2),
• DTE radio links for range/Doppler/time transfer.

Encapsulated covariance studies: NovaMoon co-located data can reduce LRF origin uncertainty to $$ds^2_\text{LCRS} = \left[1 - \frac{2}{c^2}(U_M + U^*_\text{tid})\right] c^2 d\mathcal{T}^2 - \left[1 + \frac{2}{c^2}(U_M + U^*_\text{tid})\right] d\mathcal{X}^2 + \mathcal{O}(c^{-4})$$35 cm and LRF orientation to $$ds^2_\text{LCRS} = \left[1 - \frac{2}{c^2}(U_M + U^*_\text{tid})\right] c^2 d\mathcal{T}^2 - \left[1 + \frac{2}{c^2}(U_M + U^*_\text{tid})\right] d\mathcal{X}^2 + \mathcal{O}(c^{-4})$$410 mas, with lunar clock–UTC time link $$ds^2_\text{LCRS} = \left[1 - \frac{2}{c^2}(U_M + U^*_\text{tid})\right] c^2 d\mathcal{T}^2 - \left[1 + \frac{2}{c^2}(U_M + U^*_\text{tid})\right] d\mathcal{X}^2 + \mathcal{O}(c^{-4})$$510 ns (Molli et al., 9 Feb 2026).

5. Coordinate and Time Transformations: Closed-Form and Operational Realizations

Screened forms for positions, velocities, and accelerations (in practice, used to support high-precision GNSS-class applications) are (Turyshev et al., 15 Nov 2025):

• Position ($$ds^2_\text{LCRS} = \left[1 - \frac{2}{c^2}(U_M + U^*_\text{tid})\right] c^2 d\mathcal{T}^2 - \left[1 + \frac{2}{c^2}(U_M + U^*_\text{tid})\right] d\mathcal{X}^2 + \mathcal{O}(c^{-4})$$6 in TDB): $$ds^2_\text{LCRS} = \left[1 - \frac{2}{c^2}(U_M + U^*_\text{tid})\right] c^2 d\mathcal{T}^2 - \left[1 + \frac{2}{c^2}(U_M + U^*_\text{tid})\right] d\mathcal{X}^2 + \mathcal{O}(c^{-4})$$7
• Velocity and acceleration carry similar post-Newtonian and scaling corrections.

The standard light-time model applies for one-way ranging: $$ds^2_\text{LCRS} = \left[1 - \frac{2}{c^2}(U_M + U^*_\text{tid})\right] c^2 d\mathcal{T}^2 - \left[1 + \frac{2}{c^2}(U_M + U^*_\text{tid})\right] d\mathcal{X}^2 + \mathcal{O}(c^{-4})$$8 with only Earth and Sun Shapiro terms needed for cm/ps accuracy (Turyshev et al., 15 Nov 2025).

6. Realization Uncertainties, Applications, and Performance

The best current realization of the ILRF (2010–2030 epoch) has combined origin–orientation error of 17.6 cm (origin: 15.3 cm, orientation: 8.6 cm). LLR-motivated principal axes determination achieves 2–3 cm RMS for the best stations, with mm-level local repeatability (Sośnica et al., 17 Oct 2025, Pavlov, 2019).

NovaMoon’s architecture enables differential GNSS positioning at 0.1–1 m accuracy in the South Polar region, with clock-link realization of TL monitored to $$ds^2_\text{LCRS} = \left[1 - \frac{2}{c^2}(U_M + U^*_\text{tid})\right] c^2 d\mathcal{T}^2 - \left[1 + \frac{2}{c^2}(U_M + U^*_\text{tid})\right] d\mathcal{X}^2 + \mathcal{O}(c^{-4})$$9 ns (Molli et al., 9 Feb 2026). Frame closure (BCRS$U_M(\mathcal{T},\mathcal{X}) = GM_M/R$0LCRS$U_M(\mathcal{T},\mathcal{X}) = GM_M/R$1BCRS) achieves $U_M(\mathcal{T},\mathcal{X}) = GM_M/R$2 m round-trip in operational simulations (Turyshev et al., 15 Nov 2025).

The LCRS and its associate frames support:

• Navigation and timing in cislunar and lunar surface operations (millimeter–centimeter level, sub-ns time transfer),
• Fundamental physics (tests of relativity, constraints on alternative gravity),
• Improved lunar ephemerides, orientation/EOP determination, and interior modeling,
• Interoperability for Earth-Moon system users (GNSS, LunaNet, Artemis, commercial missions).

7. Practical Implementation and Maintenance

Recommended practice is to use TDB-compatible ephemerides (DE440/INPOP21a) without rescaling, apply only those corrections needed above the mission’s required error budget (e.g., position $U_M(\mathcal{T},\mathcal{X}) = GM_M/R$3 m, velocity $U_M(\mathcal{T},\mathcal{X}) = GM_M/R$4 m/s, acceleration $U_M(\mathcal{T},\mathcal{X}) = GM_M/R$5 m/s²), and verify with fixed vector input–output closures (Turyshev et al., 15 Nov 2025, Fienga et al., 2024). Time-transfer chains UTC→TT→TCG→TCB→TCL→TL are constructed according to IAU and IERS conventions, incorporating redshift, periodic, and kinematic correction terms as required (Bourgoin et al., 29 Jul 2025, Turyshev, 29 Jul 2025).

Operationally, LLR/VLBI sessions, DTE radio tracking, and clock transfer are coordinated by the ground segment (ESA/ILRS/VGOS/LunaNet), with real-time differential corrections and time-tags distributed to user segments (landers, rovers, crew, orbiters). Product delivery includes updated lunar ephemerides, LCRS–ICRF rotation matrices, LRF realizations, and TL–UTC offsets. The system is fully SI-traceable through established Earth–Moon time links (Molli et al., 9 Feb 2026, Turyshev et al., 15 Nov 2025).

---

Relevant literature: (Turyshev et al., 15 Nov 2025, Molli et al., 9 Feb 2026, Sośnica et al., 17 Oct 2025, Bourgoin et al., 29 Jul 2025, Fienga et al., 2024, Turyshev et al., 2024, Turyshev, 29 Jul 2025, Pavlov, 2019)