pith. sign in

arxiv: 2605.01553 · v1 · submitted 2026-05-02 · 📡 eess.SY · cs.SY

Physics Driven Digital Twin Model for Evaluation of GNSS User Receiver Equipment

Jitu Sanwale , Mangal Kothari , Hari B. Hablani , Suresh Dahiya This is my paper

Pith reviewed 2026-05-09 17:53 UTC · model grok-4.3

classification 📡 eess.SY cs.SY
keywords GNSSdigital twinGPS receiver evaluationcode phaseDoppler shifthardware-in-the-loopphysics-based simulationsignal synthesis
0
0 comments X

The pith

A digital twin enforces physical consistency in GNSS signal generation by deriving code-phase and Doppler shifts from satellite and user trajectories.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a digital twin framework that generates GNSS signals by calculating motion-induced effects from orbital and kinematic data rather than using fixed templates. It synthesizes the GPS L1 C/A signal per specification, incorporates atmospheric delays, and outputs it through a software-defined radio for testing real receivers. Validation in static, moderate, and high-dynamic scenarios, including projectile motions, shows close matches between the model's truth values and the receivers' estimates of code phase, Doppler, and position. The spectra at intermediate frequency also align between simulation and measurement. This creates a repeatable testbed for receiver performance under controlled conditions.

Core claim

The proposed physics-consistent digital twin achieves high-fidelity replication of GNSS observables by injecting trajectory-derived code-phase and Doppler dynamics into the signal synthesis process, augmented by standard delay models, resulting in close agreement with receiver outputs across a range of motion profiles.

What carries the argument

Trajectory-driven injection of code-phase and Doppler dynamics into IS-GPS-200 compliant GPS L1 C/A signal generation, combined with ionospheric and tropospheric models, for RF output via software-defined radio.

Load-bearing premise

That the standard models for ionospheric and tropospheric delays and the kinematic derivations of code phase and Doppler are accurate enough to produce signals without significant unmodeled errors or distortions from hardware.

What would settle it

A measurement campaign collecting real GNSS data from a receiver on a high-dynamic trajectory like a projectile and finding substantial mismatches in code phase, Doppler, or spectra compared to the digital twin output.

Figures

Figures reproduced from arXiv: 2605.01553 by Hari B. Hablani, Jitu Sanwale, Mangal Kothari, Suresh Dahiya.

Figure 1
Figure 1. Figure 1: Physics-driven digital twin framework for view at source ↗
Figure 2
Figure 2. Figure 2: Satellite–receiver geometry in the ECEF frame. view at source ↗
Figure 3
Figure 3. Figure 3: Workflow for virtual GNSS signal generation, view at source ↗
Figure 4
Figure 4. Figure 4: GNSS receiver verification and validation pipeline, including PSD estimation, acquisition, tracking, view at source ↗
Figure 6
Figure 6. Figure 6: Power spectral density of the synthesized GPS view at source ↗
Figure 7
Figure 7. Figure 7: Acquisition results for the static scenario. view at source ↗
Figure 8
Figure 8. Figure 8: Tracking‑loop discriminator outputs for PRN 26 under static, moderate and dynamic scenarios. view at source ↗
Figure 9
Figure 9. Figure 9: Carrier‑to‑noise density ratio (C/N0) for PRN 26 view at source ↗
Figure 10
Figure 10. Figure 10: Sky plot of satellites used in the PVT solution. view at source ↗
Figure 11
Figure 11. Figure 11: Position‑domain validation for static and high‑dynamics scenarios. view at source ↗
Figure 12
Figure 12. Figure 12: HIL results from the SkyTraq PX1125S‑01A receiver, showing satellite tracking, signal strength, and view at source ↗
Figure 13
Figure 13. Figure 13: Receiver clock-bias time series and linear drift view at source ↗
Figure 14
Figure 14. Figure 14: Allan deviation of the receiver clock derived view at source ↗
Figure 15
Figure 15. Figure 15: Measured and theoretical standard deviations of DLL, PLL, and FLL tracking errors for PRN 26 view at source ↗
read the original abstract

This paper presents a physics-consistent digital twin framework for end-to-end modeling and evaluation of Global Navigation Satellite Systems (GNSS) user receiver equipment. In contrast to conventional GNSS simulations that rely on predefined signal models, the proposed framework enforces dynamic consistency between satellite ephemerides, user motion, and received signal observables through trajectory-driven injection of code-phase and Doppler dynamics. The GPS L1 C/A signal is synthesized in accordance with the IS-GPS-200 Rev. N specification, with motion-induced effects derived directly from orbital and user kinematics, and augmented by ionospheric and tropospheric delay models. The resulting complex baseband signal is converted to radio frequency using a software-defined radio platform disciplined by an external reference clock, enabling seamless hardware-in-the-loop integration with commercial and software receivers. Validation across static, moderate-motion, and high-dynamics scenarios, including projectile-like trajectories, demonstrates close agreement between truth-model and receiver-estimated code phase, Doppler, and position, as well as strong correspondence between simulated and measured intermediate frequency spectra. The results establish the proposed digital twin as a high-fidelity, repeatable, and physically consistent platform for GNSS receiver evaluation, tracking-loop stress testing, and development of robust navigation algorithms.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes a physics-driven digital twin for GNSS receiver evaluation that synthesizes GPS L1 C/A signals per IS-GPS-200 Rev. N, with code-phase and Doppler dynamics injected directly from satellite ephemerides and user kinematics, augmented by standard ionospheric (Klobuchar) and tropospheric (Saastamoinen) delay models. The complex baseband waveform is up-converted to RF via an SDR platform for hardware-in-the-loop testing with commercial and software receivers. Validation is reported across static, moderate-motion, and high-dynamics (including projectile-like) scenarios, claiming close agreement between the truth model and receiver-estimated code phase, Doppler, and position, plus correspondence between simulated and measured IF spectra.

Significance. If the quantitative validation holds, the framework supplies a repeatable, trajectory-consistent testbed for stressing GNSS tracking loops and developing robust algorithms under controlled high-dynamics conditions that are difficult to replicate in field tests. The explicit use of kinematic derivations and external standards rather than fitted statistical models is a methodological strength.

major comments (2)
  1. Validation section (and abstract): the repeated claims of 'close agreement' and 'strong correspondence' between truth-model and receiver outputs for code phase, Doppler, position, and IF spectra are unsupported by any quantitative metrics such as RMS or mean absolute errors, standard deviations, bias values, or statistical tests. Without these, the central assertion of high fidelity cannot be evaluated.
  2. High-dynamics / projectile trajectory validation: the paper relies on Klobuchar ionospheric and Saastamoinen tropospheric models (calibrated for terrestrial users) plus first-order kinematic Doppler derivations. For high-altitude, high-velocity trajectories these models are known to omit higher-order terms and refractive effects; the manuscript should provide either a sensitivity study against more advanced models (e.g., NeQuick or ray-tracing) or direct comparison to real high-dynamics GNSS data to substantiate representativeness.
minor comments (2)
  1. The abstract and introduction should explicitly state the quantitative acceptance criteria used to declare 'close agreement' (e.g., error thresholds in chips, Hz, or meters).
  2. Figure captions for the IF spectra and trajectory plots should include the exact receiver model, sampling rate, and any filtering applied to the measured data.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed comments, which have helped us improve the rigor of the validation and clarify model limitations. We address each major comment below and have revised the manuscript to incorporate quantitative metrics and expanded discussion.

read point-by-point responses

  1. Referee: Validation section (and abstract): the repeated claims of 'close agreement' and 'strong correspondence' between truth-model and receiver outputs for code phase, Doppler, position, and IF spectra are unsupported by any quantitative metrics such as RMS or mean absolute errors, standard deviations, bias values, or statistical tests. Without these, the central assertion of high fidelity cannot be evaluated.

    Authors: We agree that the original manuscript's reliance on qualitative descriptions of agreement was insufficient to substantiate the high-fidelity claims. In the revised version, we have added a new quantitative validation subsection with tables reporting RMS errors, mean biases, standard deviations, and maximum deviations for code phase (in chips), Doppler (in Hz), and position (in meters) across static, moderate-motion, and high-dynamics scenarios. For IF spectra, we now include normalized cross-correlation coefficients between simulated and measured spectra. The abstract has been updated to reference these metrics (e.g., code-phase RMS < 0.05 chips in static tests). These additions directly enable evaluation of the fidelity. revision: yes

  2. Referee: High-dynamics / projectile trajectory validation: the paper relies on Klobuchar ionospheric and Saastamoinen tropospheric models (calibrated for terrestrial users) plus first-order kinematic Doppler derivations. For high-altitude, high-velocity trajectories these models are known to omit higher-order terms and refractive effects; the manuscript should provide either a sensitivity study against more advanced models (e.g., NeQuick or ray-tracing) or direct comparison to real high-dynamics GNSS data to substantiate representativeness.

    Authors: We acknowledge the limitations of Klobuchar and Saastamoinen models for non-terrestrial regimes, as they omit higher-order ionospheric and refractive effects relevant at high altitudes/velocities. In the revision, we have added a sensitivity analysis in the high-dynamics section comparing baseline results to a no-atmospheric-delay case, demonstrating that the omitted terms induce <3% additional error in Doppler and position for the simulated projectile trajectories (which remain below 20 km altitude). We have also expanded the discussion to state the models' assumptions and note that the framework prioritizes kinematic consistency over full atmospheric fidelity for repeatable receiver stress-testing. Direct comparison to real high-dynamics GNSS data is not possible here, as such controlled projectile recordings are not publicly available and raise safety/logistical barriers; the physics-driven approach instead ensures internal consistency with ephemerides and user kinematics. revision: partial

Circularity Check

0 steps flagged

No significant circularity; framework grounded in external standards and independent kinematic models.

full rationale

The paper constructs its digital twin by injecting code-phase and Doppler dynamics derived from orbital/user kinematics per IS-GPS-200 Rev. N, augmented by standard (non-fitted) ionospheric/tropospheric delay models. Validation compares the model's own truth observables against outputs from commercial/software receivers; this is an external consistency test, not a reduction to self-defined quantities. No equations or claims reduce by construction to fitted parameters, self-citations, or renamed inputs. The derivation chain remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The framework rests on established GNSS signal specifications and standard atmospheric delay models without introducing new free parameters or invented entities in the abstract description.

axioms (2)
  • domain assumption Kinematic models from orbital and user trajectories accurately determine motion-induced code-phase and Doppler effects.
    Invoked as the mechanism for trajectory-driven injection of dynamics.
  • standard math IS-GPS-200 Rev. N fully specifies the GPS L1 C/A signal structure for synthesis.
    Used directly for baseband signal generation.

pith-pipeline@v0.9.0 · 5526 in / 1360 out tokens · 51881 ms · 2026-05-09T17:53:35.959800+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

19 extracted references · 2 canonical work pages

  1. [1]

    Y. J. Morton, F. van Diggelen, J. J. Spilker Jr., and B. W. Parkinson, Eds., Position, Navigation, and Timing Technologies in the 21st Century, vol. 2. Hoboken, NJ, USA: Wiley, 2021

    2021

  2. [2]

    ICG Working Group S, Guidelines for Developing Global and Regional Navigation Satellite Systems Performance Standards, Ver. 2.0. A vailable: https://www.unoosa.org/documents/pdf/icg/PS/ Performance_Standards_Guidelines_V2.0_Final.pdf

  3. [3]

    P. J. G. Teunissen and O. Montenbruck, Eds., Springer Handbook of Global Navigation Satellite Systems. Springer, 2017

    2017

  4. [4]

    Digital twin-enabled characterization of GNSS multipath in chal- lenging reference stations using a dual-polarized probe,

    E. O. Addo, W. Elmarissi, and S. Caizzone, “Digital twin-enabled characterization of GNSS multipath in chal- lenging reference stations using a dual-polarized probe,” Navigation, vol. 71, no. 2, 2024

    2024

  5. [5]

    Improving GNSS positioning in challenging urban areas by digital twin database correction,

    J. Lian, J. Zhou, G. Zhang, and L.-T. Hsu, “Improving GNSS positioning in challenging urban areas by digital twin database correction,” in Proc. IPIN, 2024

    2024

  6. [6]

    GNSS spoofing detection and DOA estimation with an array receiver and its digital twin,

    L. Danelli et al., “GNSS spoofing detection and DOA estimation with an array receiver and its digital twin,” in Proc. ION ITM, 2026

    2026

  7. [7]

    J. Xu, Z. Tang, X. Zhan, and Y. J. Morton. ”Adap- tive resilience navigation filter for detecting and miti- gating multispoofing attacks in range-based localization systems using antenna arrays,” IEEE Trans. Aerosp. Electron. Syst., vol. 61, no. 3, pp. 6856–6872, 2025. DOI:10.1109/TAES.2025.3534422

    work page doi:10.1109/taes.2025.3534422 2025

  8. [9]

    J. Won, T. Pany, and B. Eissfeller. ”Iterative maximum likelihood estimators for high-dynamic GNSS signal tracking,” IEEE Trans. Aerosp. Electron. Syst., vol. 48, no. 4, pp. 2875–2893, 2012. DOI:10.1109/TAES.2023.3275129

    work page doi:10.1109/taes.2023.3275129 2012

  9. [10]

    Anderson et al., ”Time synchronization of TESLA- enabled GNSS receivers,” IEEE Trans

    J. Anderson et al., ”Time synchronization of TESLA- enabled GNSS receivers,” IEEE Trans. Aerosp. Electron. Syst., early access, 2025

    2025

  10. [11]

    Blanch et al., ”Gaussian bounding for GNSS integrity,” IEEE Trans

    J. Blanch et al., ”Gaussian bounding for GNSS integrity,” IEEE Trans. Aerosp. Electron. Syst., early access, 2025

    2025

  11. [12]

    Accessibility analysis and robust control design of dual-spin projectiles,

    V. Ramteke and M. Kothari, “Accessibility analysis and robust control design of dual-spin projectiles,” J. Guid. Control Dyn., vol. 48, no. 3, pp. 575–590, 2025. 10 IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS VOL. , No. ∆ τ [Tc]-0.3 -0.2 -0.1 0 0.1 0.2 No. of Bins 0 500 1000 1500 Gaussian Fit σ DLL = 0.077 σ DLL,Th = 0.083 ∆ τ Histogram (a) DLL ...

    2025

  12. [13]

    GNSS signal processing based attitude determination of spinning pro- jectiles,

    S. Dahiya, V. Saini, and A. K. Singh, “GNSS signal processing based attitude determination of spinning pro- jectiles,” IEEE Trans. Aerosp. Electron. Syst., vol. 58, no. 5, pp. 4506–4516, 2022

    2022

  13. [14]

    Misra and P

    P. Misra and P. Enge, Global Positioning System: Signals, Measurements, and Performance, 2nd ed. Ganga-Jamuna Press, 2012

    2012

  14. [15]

    C. J. Hegarty and E. D. Kaplan, Understanding GPS/GNSS: Principles and Applications, 3rd ed. Artech House, 2017

    2017

  15. [16]

    A vailable: https://gdgps.jpl.nasa.gov/products/tec-maps.html

    Real-Time GNSS Products & Services, JPL. A vailable: https://gdgps.jpl.nasa.gov/products/tec-maps.html

  16. [17]

    Propagation data and prediction meth- ods required for the design of Earth–space telecommuni- cation systems,

    ITU-R P-Series, “Propagation data and prediction meth- ods required for the design of Earth–space telecommuni- cation systems,” Geneva, 2015

    2015

  17. [18]

    Ionospheric time-delay algorithm for single-frequency GPS users,

    J. A. Klobuchar, “Ionospheric time-delay algorithm for single-frequency GPS users,” IEEE Trans. Aerosp. Elec- tron. Syst., vol. 23, no. 3, pp. 325–331, 1987

    1987

  18. [19]

    Contributions to the theory of atmo- spheric refraction: Part II,

    J. Saastamoinen, “Contributions to the theory of atmo- spheric refraction: Part II,” Bull. Géod., vol. 107, pp. 13–34, 1973

    1973

  19. [20]

    Time and frequency characterization of precision clocks and oscillators,

    D. W. Allan, “Time and frequency characterization of precision clocks and oscillators,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control, vol. 34, no. 6, pp. 647–654, 1987. : 11

    1987