REVIEW 2 major objections 6 minor 33 references
Earth-baseline VLBI restores absolute observability of a lunar surface station when the satellite constellation is too sparse to carry the datum alone.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · grok-4.5
2026-07-12 10:28 UTC pith:M3USOQWZ
load-bearing objection Clean algebraic observability result for lunar joint OD+clock: VLBI restores the absolute station datum when the constellation is sparse and only sharpens it when rich; the three-station floor is real inside the snapshot model. the 2 major comments →
Earth-baseline VLBI restores the observability of a lunar surface station in joint orbit-and-clock determination
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In a snapshot joint orbit-and-clock fit, internal lunar observables leave a purely positional six-dimensional rigid-body datum defect. An indirect Earth-to-satellite tie reaches the surface station only for rich constellation geometry; a direct Earth-baseline VLBI delay restores absolute observability when the constellation cannot, and merely sharpens the Cramér–Rao bound when it can. The absolute position becomes fully observable only once three non-collinear Earth stations supply plane-of-sky baselines.
What carries the argument
Theorem 1 (internal datum defect): the six rigid-body generators of the lunar cluster lie in the null space of the Fisher information built from station-to-satellite and inter-satellite ranges, so the residual defect after a clock-sync is purely positional and must be closed by an Earth-frame observable.
Load-bearing premise
The numerical design law and three-station threshold rest on a single-epoch closed-loop simulation that generates truth from the identical measurement model used for recovery, without force-model dynamics or a delay-rate arc.
What would settle it
Real Earth-baseline VLBI delays to an emplaced lunar surface beacon that leave the absolute-position Fisher information rank-deficient even with three non-collinear Earth stations, or that produce station errors systematically far above the predicted Cramér–Rao bound under comparable geometry and noise.
If this is right
- Sparse or early lunar constellations need a direct surface VLBI beacon to make absolute station position observable at all.
- Once satellite geometry is rich enough for the indirect Earth-to-satellite anchor to reach the station, VLBI becomes an optional accuracy refinement rather than a requirement.
- Single-epoch absolute-datum designs must budget at least three non-collinear Earth stations; fewer leave a residual null space.
- Because a standard Gauss–Newton batch estimator attains the bound, geometry, not solver choice, sets the achievable accuracy.
- Worst-axis (E-optimal) design is the matching criterion, since plane-of-sky baselines starve one absolute axis.
Where Pith is reading between the lines
- Adding a multi-epoch VLBI delay-rate over an Earth-rotation arc would supply a second rank-one term per baseline and could soften the three-station threshold.
- The same free-network rigid-body analysis applies to surface beacons on any airless body whose local ranging cannot see the inertial frame.
- The first decisive empirical test once a surface beacon exists is whether the observed null-space dimension follows the predicted ladder of three, then one, then zero residual defects.
- Any multi-provider lunar navigation architecture that shares only local ranging and an absolute frame deferred to later documents inherits the same conditional VLBI necessity.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies absolute observability of a lunar surface station in a joint orbit-and-clock snapshot batch fit. It proves that station-to-satellite and inter-satellite ranging plus a clock-sync leave a six-dimensional rigid-body datum defect (Theorem 1), so an Earth/inertial-frame tie is required. Two ties are distinguished: an indirect Earth-to-satellite ranging path that reaches the station only when constellation geometry is rich, and a direct Earth-baseline VLBI delay to the station beacon that fixes it regardless. The resulting conditional design law is that VLBI restores absolute observability for a sparse three-satellite constellation (station in ker M until VLBI is added; CRLB 20.1 m) and only sharpens the bound for a rich six-satellite constellation (23.2 m to 9.7 m). A single-epoch baseline informs at most two axes, so the datum closes at three non-collinear Earth stations. The Gauss–Newton estimator attains the CRLB (efficiency 1.02); a 91× median station-error improvement is reported for the sparse ensemble. The FIM/CRLB engine is externally validated; the lunar network application is labelled Modelled.
Significance. If the geometric claim holds inside the stated model—and the algebra of Theorem 1 plus the validated FIM engine make that case strong—the paper supplies a usable conditional design rule for lunar PNT infrastructure: when an Earth-baseline VLBI tie (e.g., NovaMoon-class) is necessary versus merely helpful. That is a concrete architecture input for Moonlight/LCNS, LunaNet, and related programmes, and it fills a gap left by prior work that either assumed a known reference station or estimated only satellite states. Strengths that should be credited explicitly include: the machine-checked rigid-body null-space argument; external validation of the FIM/CRLB engine against NumPy and Kay closed forms; the information-not-count control; honest Validated/Modelled labelling; open reproducible tooling with fixed seeds; and a clear complementary framing relative to Pöhlmann et al. and Iiyama & Gao (Table 1). The result is geometric and falsifiable within the snapshot model.
major comments (2)
- §8 (Discussion) and (2)/§4.4: The operational sizing guidance (three-station floor, E-optimality, diminishing-returns CRLB curve) is derived under the single-epoch snapshot model. Limitations §7 correctly notes that a delay-rate observable over an Earth-rotation arc would add a second rank-one term per baseline and could soften the three-station threshold. As written, §8 presents the floor and the design curve as mission-sizing rules without restricting them to single-epoch architectures or quantifying how a multi-epoch/delay-rate extension would change them. Please either (a) explicitly scope the design law to snapshot/single-epoch ties in §8 and the abstract, or (b) add a short multi-epoch sensitivity argument so a mission team cannot take the three-station floor as absolute.
- §3.1 and §4.2: The station clock-sync pseudo-observable is introduced so that the residual datum defect is purely positional and the range-only configuration is not trivially rank-deficient. That is methodologically clean for the rank analysis, but the manuscript should state more clearly what real-system counterpart is assumed (e.g., Earth two-way time transfer, a prior absolute clock tie) and whether residual clock–position coupling could re-enter the null space if that absolute clock pin is imperfect. Without that sentence, readers may over-read Theorem 1 as applying unchanged to a pure range-only lunar network.
minor comments (6)
- §2, “Frames, time and tracking”: the heading is broken as “F rames” in the source; fix the typography.
- Table 2 and §5.1: when the station is unobservable without VLBI, the seed-42 “2183 m → 3.55 m (615×)” contrast is start-point residual versus estimate. The ensemble median 91× and the null-space language are the honest figures; consider demoting the 615× seed-42 number further in the text so it cannot be quoted as the headline.
- Figure 1 caption and §5.1: state explicitly that the range-only “error” when the state is in ker M is not a CRLB but a solver start-point residual, to match the theory section’s language.
- §5.6 / Table 3: the ELFO re-run is valuable; add one sentence on how many planes and what Earth-station geometry were used so the 0.4 m / 0.48 m bounds can be regenerated without reading the scripts.
- References: several arXiv preprints are cited as 2025–2026; ensure final DOIs/versions are updated at camera-ready if available, and that the kshana Zenodo DOI and engine tag match the committed v0.23.0 claim in §10.
- Notation: H_int, M, and ker M are clear in §4; introduce the same symbols once in the abstract or introduction so the “null space of the Fisher information” claim is not purely verbal for non-FIM readers.
Circularity Check
No significant circularity: Theorem 1 and the restore/sharpen law follow from range invariance under SE(3) and FIM rank, not from fitted inputs or load-bearing self-citation.
specific steps
-
self citation load bearing
[§6 / Availability; citations [3],[11]]
"The engine is open source (kshana, AGPL-3.0-only... Version 0.23.0 (the release carrying the Fisher-information/CRLB observability layer)... every figure is deterministic and rebuilds from a committed scenario, seed list, and engine version"
The figures and FIM computations are produced by the author’s own toolkit. This is ordinary self-citation for reproducibility and is not load-bearing: the algebraic content of Theorem 1 and the rank statements do not depend on the engine’s uniqueness, and the engine itself is cross-checked against independent NumPy/Kay oracles. Mild elevation of score only.
full rationale
The load-bearing claim (six-dimensional rigid-body datum defect of internal ranges, closed only by an Earth-frame tie, with VLBI restoring observability for sparse geometry and merely sharpening for rich geometry, and the three-station threshold) is exact algebra of the measurement model. Pairwise ranges are invariant under rigid motions by construction (Theorem 1); the null-space dimension of M = HᵀWH is then read off directly, and the engine that computes it is externally validated against NumPy and Kay’s closed forms. The numerical CRLBs and 91× factor are labelled Modelled closed-loop simulations of that geometry, not predictions fitted to data. Self-citations are to the author’s open kshana engine (for reproducibility of the same algebra) and to standard external references; none supply a uniqueness theorem or ansatz that forces the central result. The single-epoch modelling boundary is openly stated and does not create a circular reduction. Score 1 only for the minor presence of author-engine citations that are not load-bearing.
Axiom & Free-Parameter Ledger
free parameters (4)
- VLBI delay sigma =
1e-11 s
- lunar range / ISL sigma =
0.1 m
- constellation size and orbit shape =
3 or 6 satellites
- Earth-station count and geometry =
1–10 stations
axioms (5)
- standard math Pairwise ranges are invariant under rigid translations and infinitesimal rotations of the entire lunar cluster (station + satellites).
- domain assumption Observations are Gaussian with diagonal covariance; Fisher information is therefore H^T W H.
- domain assumption A single station clock-sync pseudo-observable renders all clocks observable, leaving a purely positional datum defect.
- domain assumption Earth stations subtend a tiny angle from the Moon, so each single-epoch VLBI baseline gradient lies in the plane of the sky and informs at most two absolute axes.
- ad hoc to paper Snapshot batch least-squares with no force-model propagation, no unmodelled forces, and closed-loop truth generation is sufficient to expose the geometric observability structure.
read the original abstract
Lunar positioning, navigation, and timing (PNT) is moving from concept to hardware, ESA's Moonlight/LCNS, NovaMoon reference stations, LunaNet, and Coordinated Lunar Time, all reducing to one estimation core: fix the orbits and clocks of the lunar infrastructure and tie them to an Earth/inertial frame. We ask which measurements make a surface station's absolute position observable, and prove the answer. In a snapshot batch fit, the internal observables (station-to-satellite and inter-satellite ranging plus clock-sync) constrain only relative geometry and leave a six-dimensional rigid-body datum defect: three translations and three rotations of the cluster. The clocks are fully observable, so the defect is purely positional, and closing it needs a tie to the Earth frame. Two such ties exist and are not interchangeable. An indirect tie (Earth-to-satellite ranging through the constellation) reaches the station only when the satellite geometry is rich; a direct tie (an Earth-baseline VLBI delay to the station beacon) fixes it regardless. This gives a conditional design law, not a single number: VLBI restores absolute observability when the constellation cannot supply it, and merely sharpens the bound when it can. For a sparse three-satellite constellation the station lies in the null space of the Fisher information until VLBI is added, reaching a Cramer-Rao bound of 20.1 m; for a rich six-satellite constellation VLBI tightens the bound from 23.2 m to 9.7 m. A single-epoch baseline informs at most two of three axes, so the datum closes at three non-collinear Earth stations. The Gauss-Newton estimator attains the bound (efficiency 1.02), with a 91x median station-error improvement in the sparse regime. The FIM/CRLB engine is validated against NumPy and published closed forms; the lunar application stays modelled, every figure deterministic and reproducible.
Figures
Reference graph
Works this paper leans on
-
[1]
Brent A. Archinal, Charles H. Acton, Michael F. A’Hearn, et al. Report of the IAU working group on cartographic coordinates and rotational elements: 2015.Celestial Mechanics and Dynamical Astronomy, 130(3):22, 2018. doi: 10.1007/s10569-017-9805-5
-
[2]
Neil Ashby and Bijunath R. Patla. A relativistic framework to establish coordinate time on the moon and beyond.The Astronomical Journal, 168(3):112, 2024. doi: 10.3847/ 1538-3881/ad643a. Canonical surface-to-surface secular Earth–Moon rate≈56.02µs/day; arXiv:2402.11150
Pith/arXiv arXiv 2024
-
[3]
Kshana: a PNT-resilience simulator with quantum-sensor performance models
Chakshu Baweja. Kshana: a PNT-resilience simulator with quantum-sensor performance models. Zenodo, 2026. URLhttps://github.com/ashfordeOU/kshana. Version 0.23.0 (the release carrying the Fisher-information/CRLB observability layer); concept DOIhttps: //doi.org/10.5281/zenodo.20528627resolves to the latest release
-
[4]
Stefano Bertone, Daniel Arnold, Vincent Girardin, Adrian Jäggi, et al. Assessing reduced- dynamic parametrizations for GRAIL orbit determination and the recovery of independent lu- nar gravity field solutions.Earth and Space Science, 8(6), 2021. doi: 10.1029/2020EA001454. Establishes that empirical accelerations absorb force/frame/ephemeris mismodelling; ...
-
[5]
Sriramya Bhamidipati, Tara Mina, Adrian Sanchez, and Grace Gao. Satellite constellation design for a lunar navigation and communication system.NA VIGATION: Journal of the Institute of Navigation, 70(4):navi.613, 2023. doi: 10.33012/navi.613. 18
-
[6]
Adrien Bourgoin, Pascale Defraigne, and Frédéric Meynadier. Lunar reference timescale. Metrologia, 63(1):015003, 2026. doi: 10.1088/1681-7575/ae2c03. BIPM Time Department and Royal Observatory of Belgium; arXiv:2511.02709
-
[7]
Orekit: an accurate and efficient core layer for space flight dynamics applications.https://www.orekit.org, 2024
CS GROUP, Luc Maisonobe, et al. Orekit: an accurate and efficient core layer for space flight dynamics applications.https://www.orekit.org, 2024
2024
-
[8]
Curkendall and James S
David W. Curkendall and James S. Border. Delta-DOR: The One-Nanoradian Navigation Measurement System of the Deep Space Network — History, Architecture, and Componentry. Interplanetary Network Progress Report 42-193, Jet Propulsion Laboratory, 2013
2013
-
[9]
Moonlight: connecting Earth with the Moon.https://www.esa
European Space Agency. Moonlight: connecting Earth with the Moon.https://www.esa. int/Applications/Connectivity_and_Secure_Communications/Moonlight, 2024. Lunar Communications and Navigation Service (LCNS): five lunar satellites (one communications, four navigation) in elliptical lunar frozen orbits, south-pole priority
2024
-
[10]
Celestial Time Standardization in Cislunar Space
Executive Office of the President, Office of Science and Technology Policy. Celestial Time Standardization in Cislunar Space. Policy memorandum, 2 April 2024, 2024. Directs NASA to deliver a Coordinated Lunar Time (LTC) standard by end of 2026
2024
-
[11]
Continuous reproducibility in GNSS signal processing.IEEE Access, 6:20451–20463, 2018
Carles Fernández-Prades, Jordi Vila-Valls, Javier Arribas, and Antonio Ramos. Continuous reproducibility in GNSS signal processing.IEEE Access, 6:20451–20463, 2018. doi: 10.1109/ ACCESS.2018.2822835
arXiv 2018
-
[12]
Lunar Reference Systems, Frames and Time-scales in the context of the ESA Programme Moonlight, 2024
Agnès Fienga, Nicolas Rambaux, and Krzysztof Sośnica. Lunar Reference Systems, Frames and Time-scales in the context of the ESA Programme Moonlight, 2024. URLhttps://arxiv. org/abs/2409.10043
Pith/arXiv arXiv 2024
-
[13]
IAU SOFA Software Collection.http://www.iausofa.org, 2023
IAU SOFA Board. IAU SOFA Software Collection.http://www.iausofa.org, 2023
2023
-
[14]
GNSS-basedLunarOrbitandClockEstimationWithaStochas- tic Cloning UD Filter, 2026
KeidaiIiyamaandGraceGao. GNSS-basedLunarOrbitandClockEstimationWithaStochas- tic Cloning UD Filter, 2026. URLhttps://arxiv.org/abs/2601.16393. Submitted to AIAA J. Guidance, Control, and Dynamics
arXiv 2026
-
[15]
LuPNT: An open-source simulator for lunar positioning, navigation, and timing
Keidai Iiyama, Guillem Casadesús Vila, and Grace Gao. LuPNT: An open-source simulator for lunar positioning, navigation, and timing. InProc. ION GNSS+ 2023. Institute of Navigation,
2023
-
[16]
Stanford NavLab open-source lunar PNT simulator
-
[17]
Resolutions of the XXXII General Assembly: Lunar Ce- lestial Reference System (LCRS) and Lunar Coordinate Time (TCL)
International Astronomical Union. Resolutions of the XXXII General Assembly: Lunar Ce- lestial Reference System (LCRS) and Lunar Coordinate Time (TCL). XXXII IAU General Assembly, Cape Town, August 2024, 2024. Establishes a lunar celestial reference system and a lunar coordinate time scale (TCL) consistent with the IAU relativistic framework
2024
-
[18]
Kay.Fundamentals of Statistical Signal Processing: Estimation Theory
Steven M. Kay.Fundamentals of Statistical Signal Processing: Estimation Theory. Prentice Hall, Upper Saddle River, NJ, 1993. ISBN 978-0133457117
1993
-
[19]
Kanhere, Daniel Neamati, and Grace Gao
Derek Knowles, Ashwin V. Kanhere, Daniel Neamati, and Grace Gao. gnss_lib_py: Analyzing gnss data with python.SoftwareX, 27:101811, 2024. doi: 10.1016/j.softx.2024.101811
-
[20]
Lemoine, Sander Goossens, Terence J
Frank G. Lemoine, Sander Goossens, Terence J. Sabaka, Joseph B. Nicholas, Erwan Mazarico, David D. Rowlands, Bryant D. Loomis, Douglas S. Chinn, Douglas S. Caprette, Gregory A. Neumann, David E. Smith, and Maria T. Zuber. High-degree gravity models from GRAIL 19 primary mission data.Journal of Geophysical Research: Planets, 118(8):1676–1698, 2013. doi: 10...
-
[21]
Precise orbit determination of LEO satellites: A systematic review.GPS Solutions, 27, 2023
Xingxing Li et al. Precise orbit determination of LEO satellites: A systematic review.GPS Solutions, 27, 2023. doi: 10.1007/s10291-023-01441-5
-
[22]
Serena Molli, Agnès Fienga, Pascale Defraigne, Krzysztof Sośnica, Luigi Cacciapuoti, Javier Ventura-Traveset, et al. NovaMoon: A Strategic Lunar Reference Station for Positioning, Timing, and Largely Enhanced Science in the Earth–Moon System, 2026. URLhttps: //arxiv.org/abs/2602.08432. ESA-adopted Argonaut PNT payload (CM25, Nov. 2025); design combines an...
arXiv 2026
-
[23]
Springer-Verlag, Berlin, Heidelberg, 2000
Oliver Montenbruck and Eberhard Gill.Satellite Orbits: Models, Methods and Applications. Springer-Verlag, Berlin, Heidelberg, 2000. doi: 10.1007/978-3-642-58351-3
-
[24]
General Mission Analysis Tool (GMAT).https:// gmat.gsfc.nasa.gov, 2020
NASA Goddard Space Flight Center. General Mission Analysis Tool (GMAT).https:// gmat.gsfc.nasa.gov, 2020
2020
-
[25]
LunaNet Signal-in- Space Recommended Standard — Augmented Forward Signal (AFS), Volume A
NASA Space Communications and Navigation (SCaN) Program. LunaNet Signal-in- Space Recommended Standard — Augmented Forward Signal (AFS), Volume A. LNIS Applicable Document, 29 January 2025, 2025. URLhttps://www.nasa.gov/general/ lunanet-interoperability-specification/. S-band one-way ranging navigation signal; Lunar Augmented Navigation Service (LANS); ta...
2025
-
[26]
LunaNet Interoperabil- ity Specification (LNIS), Version 5
NASA Space Communications and Navigation (SCaN) Program. LunaNet Interoperabil- ity Specification (LNIS), Version 5. Baseline effective 29 January 2025, 2025. URL https://www.nasa.gov/general/lunanet-interoperability-specification/. Cross- agency (NASA/ESA/JAXA) lunar communications and navigation interoperability frame- work
2025
-
[27]
Ryan S. Park, William M. Folkner, James G. Williams, and Dale H. Boggs. The JPL planetary and lunar ephemerides DE440 and DE441.The Astronomical Journal, 161(3):105, 2021. doi: 10.3847/1538-3881/abd414
-
[28]
Joel J. K. Parker, Fabio Dovis, Lauren Konitzer, Nathan Esantsi, Benjamin Ashman, Alex Minetto, Andrea Nardin, Oliviero Vouch, Simone Zocca, Fabio Bernardi, Matilde Boschiero, Samuele Fantinato, Efer Miotti, Claudia Facchinetti, Mario Musmeci, and Giancarlo Varacalli. GNSS reception at the moon: First results of the lunar GNSS receiver experiment (LuGRE)....
2026
-
[29]
Towards Hybrid Lunar PNT: Error Models, Lower Bounds and Algorithms, 2025
Robert Pöhlmann, Emanuel Staudinger, and Gonzalo Seco-Granados. Towards Hybrid Lunar PNT: Error Models, Lower Bounds and Algorithms, 2025. URLhttps://arxiv.org/abs/ 2508.10699. Cramér–Rao lower bounds for lunar PNT; sub-metre accuracy with minimal satellite visibility when a reference station is available
Pith/arXiv arXiv 2025
-
[30]
Nyx: a high-fidelity astrodynamics and orbit-determination toolkit (with ANISE and Hifitime).https://nyxspace.com, 2025
Christopher Rabotin and Nyx Space contributors. Nyx: a high-fidelity astrodynamics and orbit-determination toolkit (with ANISE and Hifitime).https://nyxspace.com, 2025. Open- source Rust astrodynamics/OD; flown on Firefly Blue Ghost 1 and used for CAPSTONE. 20
2025
-
[31]
Development of the low-cost rtk-gps receiver with an open source program package rtklib.International Symposium on GPS/GNSS, 2009
Tomoji Takasu and Akio Yasuda. Development of the low-cost rtk-gps receiver with an open source program package rtklib.International Symposium on GPS/GNSS, 2009
2009
-
[32]
Telespazio signs a contract worth 123 million euro with ESA for the Moonlight pro- gramme
Telespazio. Telespazio signs a contract worth 123 million euro with ESA for the Moonlight pro- gramme. Press release, 15 October 2024 (IAC Milan), 2024. URLhttps://www.telespazio. com
2024
-
[33]
Vallado, Paul Crawford, Richard Hujsak, and T
David A. Vallado, Paul Crawford, Richard Hujsak, and T. S. Kelso. Revisiting spacetrack report #3. InAIAA/AAS Astrodynamics Specialist Conference and Exhibit, 2006. doi: 10. 2514/6.2006-6753. AIAA 2006-6753. 21
2006
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.