GNSS | |
GNSS Constellation Specific Monthly Analysis Summary: April 2025
The analysis performed in this report is solely his work and own opinion. State Program: U.S.A (G); EU (E); China (C) “Only MEO- SECM satellites”; Russia (R); Japan (J); India (I) |
![]() |
|
Introduction
The article is a continuation of monthly performance analysis of the GNSS constellation. Please refer to previous issues for past analysis. As an addition, there is a section that explores mathematical concepts of variational equations used in the estimation of satellite clocks and orbits. This is in turn linked to the challenges of estimating orbits of inclined geosynchronous and geo-stationary satellites in QZSS and Beidou constellations.
Analyzed Parameters for April 2025
(Dhital et. al, 2024) provides a brief overview of the necessity and applicability of monitoring the satellite clock and orbit parameters.
a..Satellite Broadcast Accuracy, measured in terms of Signal-In Space Range Error (SISRE) (Montenbruck et. al, 2010).
b. SISRE-Orbit (only orbit impact on the range error), SISRE (both orbit and clock impact), and SISRE-PPP (as seen by the users of carrier phase signals, where the ambiguities absorb the unmodelled biases related to satellite clock and orbit estimations. Satellite specific clock bias is removed) (Hauschlid et.al, 2020)
c. Clock Discontinuity: The jump in the satellite clock offset between two consecutive batches of data uploads from the ground mission segment. It is indicative of the quality of the satellite atomic clock and associated clock model.
d. URA: User Range Accuracy as an indicator of the confidence on the accuracy of satellite ephemeris. It is mostly used in the integrity computation of RAIM.
e. GNSS-UTC offset: It shows stability of the timekeeping of each constellation w.r.t the UTC
f. Variational Equations and Satellite Orbit Estimation: The variational equations are a key set of differential equations to reliably capture the dynamics of the satellite motion in the orbit determination process. The solution to these equations is indicative of quality of the initial orbit determination and force models
Note:– for India’s IRNSS there are no precise satellite clocks and orbits as they broadcast only 1 frequency which does not allow the dual frequency combination required in precise clock and orbit estimation; as such, only URA and Clock Discontinuity is analyzed.
Satellite orbit and clock parameters form the backbone of broadcast navigation messages, providing essential data for users worldwide. Since January 2024, monthly performance reports in mycoordinates have analyzed their characteristics and quality. This series of articles now shifts focus to the computational aspects of these parameters—not to rehash existing mathematical methods, but to distill key concepts in a simplified, yet insightful manner. By doing so, this analysis aims to enhance comprehension of the performance variations observed in different GNSS constellations. In this issue, a brief linkage to the performance characteristics of IGOS and GEO satellites is provided.
Variational equations are fundamental to modern orbit determination and adjustment algorithms. They enable efficient and accurate refinement of estimated parameters by providing sensitivity information, which is crucial for adjusting the satellite’s trajectory to match observations.
Linearization and the Role of the Jacobian
In the first step of orbit determination, the highly nonlinear GNSS observation equations are linearized around an initial state. This yields a Jacobian matrix that contains the partial derivatives of the measurements with respect to the estimated parameters—primarily the satellite’s position and velocity.
However, the satellite’s state itself evolves according to nonlinear equations of motion. Therefore, to understand how changes in the initial state affect the satellite’s position at later times, we must also linearize the dynamics. This is where the state transition matrix and sensitivity matrix come into play.
Dynamics and Variational Equations
The satellite’s trajectory is governed by a set of differential equations where acceleration is a function of position, velocity, and dynamical parameters (e.g., drag coefficient, solar radiation pressure). A particular orbit solution is uniquely determined by its initial state and dynamical parameters.
Numerical integration of the equations of motion yields a predicted orbit. However, this orbit is only accurate if the initial conditions and force models are close to the truth. To refine this orbit using observations, we need to understand how small changes in these inputs affect the output trajectory.
Why Variational Equations?
Using finite differences to compute these sensitivities of the orbit trajectory to small changes in initial values is computationally expensive and numerically unstable. Instead, variational equations provide a continuous and efficient way to compute:
Quantifies how uncertainties in dynamical parameters affect the trajectory.
These matrices are integrated alongside the equations of motion using the same force models. At each observation epoch, they provide the necessary partial derivatives for the Jacobian matrix used in batch least squares estimation, for example.
Constructing the Jacobian Matrix
The full Jacobian matrix (J) used in orbit determination is built using the chain rule:
This formulation ensures that the measurement sensitivity accounts for how errors in the initial state propagate to the observation epoch.
What If We Ignore Variational Equations?
• Loss of sensitivity to initial conditions
• Poor convergence in least squares estimation
• Inaccurate orbit corrections and higher residuals
Practical Implications
• Large values in Φ indicate strong dependence on initial conditions—small errors grow significantly over time.
• Stable orbit solutions are achieved when the sensitivity to initial conditions and dynamical parameters becomes small.
• Variational equations reduce the number of iterations needed and improve convergence speed and accuracy.
Taking examples of IGSO and GEO satellites of Beidou and QZSS, the changes in observation geometry for GEO and IGSO satellites are
much smaller compared to MEO satellites. As a result, strong correlations often arise between:
• State parameters (e.g., orbital elements),
• Dynamical parameters (e.g., solar radiation pressure coefficients),
• And measurement model parameters (e.g., ambiguities, differential code biases).
The reduced orbit quality of IGSO satellites can also be attributed to:
• Limited and regional ground station coverage,
• Radiation pressure modeling issues, especially due to large communication antennas not accounted for in standard box-wing models.
For GEO, additional SRP modeling difficulties from large communication antennas with undisclosed shape and material properties impact the orbit quality.
Moreover, delta-V maneuvers (for orbit maintenance) and attitude control maneuvers (for yaw and orbit normal attitudes) are hard to model precisely. Both GEO and IGSO satellites need frequent maneuvers in comparison to MEO satellites. All above mentioned characteristics introduce discontinuities or biases that are not easily captured by variational equations.
In fact, analysis of the sensitivity matrix often reveals large values associated with dynamical parameters—an indicator of modeling deficiencies. Ideally, a well converged solution should show:
• Diagonal elements close to 1 (indicating stable propagation),
• Off-diagonal elements that are small (indicating low parameter correlation).
To dig deeper into these challenges, future work will involve using advanced tools for detailed analysis of the sensitivity and covariance matrices in reduced dynamics orbit determination. In this regard, this section serves as an introductory overview, laying the foundation for deeper investigations into these modeling and estimation challenges.
Monthly Performance Remarks:
1. Satellite Clock and Orbit Accuracy:
• The performance of all constellations is relatively stable with some minor changes from previous month.
• The satellite clock jumps identified a couple of issues in Galileo satellites. E11 and E12 had large jumps on DOY 91 and DOY 105, respectively. Satellite E19 continues to have the noisiest clock (95 percentile is 0.36 ns vs 95 percentile of 0.2 ns of the whole constellation).
• The improvement in the QZSS satellite clock and orbit accuracy is visible. It is improved by 5 cm. Further investigation (once the QZSS operational history is available) is needed to correlate to the Sun angle and the switches between the yaw attitude and orbit normal modes of the satellites.
• The URA for I02 showed a little more scatter in comparison to previous months. It suggests a degraded confidence in its satellite orbit.
2. UTC Prediction (GNSS-UTC):
• GPS showed some variations as in previous months. Glonass UTC prediction started to deviate significantly in the second half of the month.
References
Alonso M, Sanz J, Juan J, Garcia, A, Casado G (2020) Galileo Broadcast Ephemeris and Clock Errors Analysis: 1 January 2017 to 31 July 2020, MDPI
Alonso M (2022) Galileo Broadcast Ephemeris and Clock Errors, and Observed Fault Probabilities for ARAIM, Ph.D Thesis, UPC
Bento, M (2013) Development and Validation of an IMU/GPS/Galileo Integration Navigation System for UAV, PhD Thesis, UniBW.
BIMP (2024 a) https://e-learning.bipm. org/pluginfile.php/6722/mod_label/ intro/User_manual_cggtts_analyser. pdf?time=1709905608656
BIMP (2024 b) https://e-learning. bipm.org/mod/folder/view. php?id=1156&forceview=1
BIMP (2024 c) https://cggtts analyser.streamlit.app
Bruggemann, Troy & Greer, Duncan & Walker, R.. (2011). GPS fault detection with IMU and aircraft dynamics. IEEE Transactions on Aerospace and Electronic Systems – IEEE TRANS AEROSP ELECTRON SY. 47. 305-316. 10.1109/TAES.2011.5705677.
Cao X, Zhang S, Kuang K, Liu T (2018) The impact of eclipsing GNSS satellites on the precise point positioning, Remote Sensing 10(1):94
Chen, K., Chang, G. & Chen, C (2021) GINav: a MATLAB-based software for the data processing and analysis of a GNSS/IMU integrated navigation system. GPS Solut 25, 108. https:// doi.org/10.1007/s10291-021-01144-9
Curran, James T. & Broumendan, Ali. (2017). On the use of Low-Cost IMUs for GNSS Spoofing Detection in Vehicular Applications.
Dhital N (2024) GNSS constellation specific monthly analysis summary, Coordinates, Vol XX, Issue 1, 2, 3, 4
Dhital N (2025) GNSS constellation specific monthly analysis summary, Coordinates, Vol XXI, Issue 1
GINAv (2025). https://geodesy.noaa. gov/gps-toolbox/GINav.shtml
Goercke, L (2017) GNSS-denied navigation of fixed-wing aircraft using low-cost sensors and aerodynamic motion models, PhD Thesis, TUM.
Guo F, Zhang X, Wang J (2015) Timing group delay and differential code bias corrections for BeiDou positioning, J Geod,
Hauschlid A, Montenbruck O (2020) Precise real-time navigation of LEO satellites using GNSS broadcast ephemerides, ION
IERS C04 (2024) https://hpiers.obspm.fr/ iers/eop/eopc04/eopc04.1962-now
IGS (2021) RINEX Version 4.00 https://files.igs.org/pub/data/ format/rinex_4.00.pdf
Kj, Nirmal & Sreejith, A. & Mathew, Joice & Sarpotdar, Mayuresh & Suresh, Ambily & Prakash, Ajin & Safonova, Margarita & Murthy, Jayant. (2016). Noise modeling and analysis of an IMU-based attitude sensor: improvement of performance by filtering and sensor fusion. 99126W. 10.1117/12.2234255.
Li M, Wang Y, Li W (2023) performance evaluation of real-time orbit determination for LUTAN-01B satellite using broadcast earth orientation parameters and multi GNSS combination, GPS Solutions, Vol 28, article number 52
Li W, Chen G (2023) Evaluation of GPS and BDS-3 broadcast earth rotation parameters: a contribution to the ephemeris rotation error Montenbruck
Liu, Yue & Liu, Fei & Gao, Yang & Zhao, Lin. (2018). Implementation and Analysis of Tightly Coupled Global Navigation Satellite System Precise Point Positioning/Inertial Navigation System (GNSS PPP/IMU) with IMUufficient Satellites for Land Vehicle Navigation. Sensors. 18. 4305. 10.3390/s18124305.
O, Steigenberger P, Hauschlid A (2014) Broadcast versus precise ephemerides: a multi-GNSS perspective, GPS Solutions
Liu T, Chen H, Jiang Weiping (2022) Assessing the exchanging satellite attitude quaternions from CNES/CLS and their application in the deep eclipse season, GPS Solutions 26(1)
Montenbruck O, Steigenberger P, Hauschlid A (2014) Broadcast versus precise ephemerides: a multi-GNSS perspective, GPS Solutions
Montenbruck O, Hauschlid A (2014 a) Differential Code Bias Estimation using Multi-GNSS Observations and Global Ionosphere Maps, ION
Niu, Z.; Li, G.; Guo, F.; Shuai, Q.; Zhu, B (2022) An Algorithm to Assist the Robust Filter for Tightly Coupled RTK/IMU Navigation System. Remote Sens. 2022, 14, 2449. https://doi.org/10.3390/rs14102449
Schmidt, G, Phillips, R (2010) IMU/ GPS Integration Architecture Performance Comparisons. NATO.
Steigenberger P, Montenbruck O, Bradke M, Ramatschi M (2022) Evaluation of earth rotation parameters from modernized GNSS navigation messages, GPS Solutions 26(2)
Suvorkin, V., Garcia-Fernandez, M., González-Casado, G., Li, M., & Rovira-Garcia, A. (2024). Assessment of Noise of MEMS IMU Sensors of Different Grades for GNSS/IMU Navigation. Sensors, 24(6), 1953. https:// doi.org/10.3390/s24061953
Sylvain L, Banville S, Geng J, Strasser S (2021) Exchanging satellite attitude quaternions for improved GNSS data processing consistency, Vol 68, Issue 6, pages 2441-2452
Tanil, Cagatay & Khanafseh, Samer & Pervan, Boris. (2016). An IMU Monitor agaIMUt GNSS Spoofing Attacks during GBAS and SBAS-assisted Aircraft Landing Approaches. 10.33012/2016.14779.
Walter T, Blanch J, Gunning K (2019) Standards for ARAIM ISM Data Analysis, ION
Wang, C & Jan, S (2025). Performance Analysis of MADOCA-Enhanced Tightly Coupled PPP/ IMU. NAVIGATION: Journal of the IMUtitute of Navigation March 2025, 72 (1) navi.678; DOI: https://doi.org/10.33012/navi.678
Wang N, Li Z, Montenbruck O, Tang C (2019) Quality assessment of GPS, Galileo and BeiDou-2/3 satellite broadcast group delays, Advances in Space Research
Wang J, Huang S, Lia C (2014) Time and Frequency Transfer System Using GNSS Receiver, Asia Pacific Radio Science, Vol 49, Issue 12
https://cggtts-analyser.streamlit.app
Yang N, Xu A, Xu Z, Xu Y, Tang L, Li J, Zhu H (2025) Effect of WHU/GFZ/CODE satellite attitude quaternion products on the GNSS kinematic PPP during the eclipse season, Advances in Space Research, Volume 75, Issue 1,
Note: References in this list might also include references provided to previous issues.
Data sources and Tools:
https://cddis.nasa.gov (Daily BRDC); http://ftp.aiub. unibe.ch/CODE_MGEX/CODE/ (Precise Products); BKG “SSRC00BKG” stream; IERS C04 ERP files
(The monitoring is based on following signals- GPS: LNAV, GAL: FNAV, BDS: CNAV-1, QZSS:LNAV IRNSS:LNAV GLO:LNAV (FDMA))
Time Transfer Through GNSS Pseudorange Measurements: https://e-learning.bipm.org/login/index.php
Allan Tools, https://pypi.org/project/AllanTools
gLAB GNSS, https://gage.upc.edu/en/learning materials/software-tools/glab-tool-suite.
Leave your response!