Quality monitoring of real-time GNSS precise positioning service system

ABSTRACT The Real-Time Global Navigation Satellite System (GNSS) Precise Positioning Service (RTPPS) is recognized as the most promising system by providing precise satellite orbit and clock corrections for users to achieve centimeter-level positioning with a stand-alone receiver in real-time. Although the products are available with high accuracy almost all the time, they may occasionally suffer from unexpected significant biases, which consequently degrades the positioning performance. Therefore, quality monitoring at the system-level has become more and more crucial for providing a reliable GNSS service. In this paper, we propose a method for the monitoring of real-time satellite orbit and clock products using a monitoring station network based on the Quality Control (QC) theory. The satellites with possible biases are first detected based on the outliers identified by Precise Point Positioning (PPP) in the monitoring station network. Then, the corresponding orbit and clock parameters with temporal constraints are introduced and estimated through the sequential Least Square (LS) estimator and the corresponding Instantaneous User Range Errors (IUREs) can be determined. A quality indicator is calculated based on the IUREs in the monitoring network and compared with a pre-defined threshold. The quality monitoring method is experimentally evaluated by monitoring the real-time orbit and clock products generated by GeoForschungsZentrum (GFZ), Potsdam. The results confirm that the problematic satellites can be detected accurately and effectively with missed detection rate and false alarm rate . Considering the quality alarms, the PPP results in terms of RMS of positioning differences with respect to the International GNSS Service (IGS) weekly solution in the north, east and up directions can be improved by 12%, 10% and 27%, respectively.


Introduction
The Real-Time Global Navigation Satellite System (GNSS) Precise Positioning Service (RTPPS) can offer services with scalable accuracy for a stand-alone receiver by providing various products including real-time precise orbits and clocks, Uncalibrated Phase Delay (UPD) and regional augmentation information (Ge et al. 2012;Li et al. 2013). The real-time precise orbit and clock products are the prerequisite for preforming the centimeter-level real-time Precise Point Positioning (PPP) (Zumberge et al. 1997). In addition, PPP ambiguity fixing with UPD corrections can significantly shorten its convergence time and improve its accuracy (Ge et al. 2008;Laurichesse et al. 2009;Collins et al. 2010). The convergence time can be further reduced by using precise atmospheric delay corrections retrieved from a regional reference network (Li, Zhang, and Ge 2011;Ge et al. 2012), which is also called PPP-RTK (Real Time Kinematic) service. PPP-RTK has comparable performance to Network Real-Time Kinematic (NRTK) and has become a more attractive and efficient tool (Li et al. 2020;Hexagon 2021). Based on the above advantages, RTPPS is recognized as the most promising service mode in the future.
The core feature of RTPPS is to provide real-time precise orbit and clock products, which are generated by processing data from 50 to 100 globally distributed tracking stations (Laurichesse et al. 2013;Fu et al. 2019;Li et al. 2019). These stations are equipped with highperformance GNSS receivers and transfer data stream via the Internet. The products are sent to users using Networked Transport of Radio Technical Commission for Maritime Service (RTCM) via Internet Protocol (NTRIP) (Weber et al. 2007). At present, the real-time satellite products are available with high accuracy. For example, the assessment of the CLK93 products shows that in November 2019 the Signal-in-Space Ranging Error (SISRE) values for GPS, GLONASS, Galileo, BeiDou MEO and BeiDou IGSO satellites are 2.3, 5.2, 1.6, 5.5 and 3.9 cm, respectively (Kazmierski, Zajdel, and Sośnica 2020). Therefore, based on these real-time orbit and clock products, a centimeter-level accuracy can be achieved by real-time PPP, which meets the requirements of many applications (Blewitt et al. 2006;Defraigne, Baire, and Guyennon 2007;Liu et al. 2018).
However, it is inevitable that there are outliers and/ or accuracy degradations in the real-time precise satellite products. On the one hand, the real-time orbit products are generated by orbit prediction based on the estimated orbits. The accuracy of real-time orbits is usually comparable to that of the estimated orbits, whereas the degradation in accuracy may occur due to unannounced orbit maneuvers and improper orbit modeling, such as an imperfect Solar Radiation Pressure (SRP) model for BeiDou satellites (Duan et al. 2019;Du et al. 2021). On the other hand, the satellite clock products are vulnerable to anomalies since they should be estimated precisely and updated in an interval of 5 s or less with well-distributed stations. The quality of the clock products can be affected by the availability and quality of tracking data, such as undetected cycle slips and bad station-satellite geometry (Du et al. 2021) or insufficient continuous observation due to internet problems. For example, the empirically derived probability of failure for the Trimble CenterPoint RTX satellite orbit and clock products is 1 � 10 À 6 for GPS and Galileo and 1 � 10 À 5 for GLONASS (Rodriguez-Solano et al. 2019). Thus, it is important to perform quality monitoring for the realtime satellite products.
The quality monitoring of the real-time precise satellite products at the system-level has not yet been widely discussed in the literature (Du et al. 2021), although its importance is well-known. Currently, none of the real-time products from International GNSS Service (IGS) real-time service or its Real-Time Analysis Center (RTAC) is provided along with quality information. Cheng et al. (2018) performed a preliminary analysis of the User Ranging Accuracy (URA), which is an important indicator for product quality, based on the real-time orbit and clock products provided by Center National d'Études Spatiales (CNES). Besides, Trimble CenterPoint RTX service constructed their own integrity monitoring system that uses the carrier phase observation residuals from 20 to 25 monitoring stations to check product quality both on newly generated corrections before broadcasting and on the already broadcast corrections (Weinbach et al. 2018). Wang and Shen (2020) proposed a real-time integrity monitoring method for the Wide Area Precise Positioning System (WAPPS) based on the satellite corrected residuals using ionosphere-free pseudorange and carrier-phase observations, which can monitor GPS and BeiDou orbit and clock corrections in real time. With more and more research interests in real-time PPP, it is necessary to provide comprehensive quality information of the real-time satellite orbit and clock products for guaranteeing a good real-time GNSS positioning services.
In this study, a method is proposed for quality monitoring of the real-time orbit and clock products in the state domain to fully consider different impacts of satellite orbit errors on ranges at different stations. This method and its processing procedure are developed based on PPP algorithm and the Quality Control (QC) theory and applied to a monitoring station network. The possible problematic satellites are detected according to the number of "outliers" votes based on the residuals from the parallel PPP processing lines. The product biases are represented by clock biases and orbit biases in the radial, along and cross directions, where the satellite clock bias can be assimilated by the radial component for better estimability. These parameters are estimated in the extended PPP model using data of the whole monitoring station network and the corresponding IUREs can be determined. A quality indicator can be calculated based on the IUREs in the monitoring station network and compared with the pre-defined threshold δ. An alarm for the satellite products can be issued to users if the quality indicator is larger than the threshold δ. This study is organized as follows: after this introduction, the mathematic model for identifying orbit and clock products is presented and a data processing procedure of the quality monitoring is developed and depicted in detail. Then, the monitoring network and processing strategy for experiments are described. The statistic results for the monitoring of the real-time products generated by the GFZ IGS RTAC are analyzed and a typical case is investigated carefully by comparing with the GeoForschungsZentrum (GFZ) final multi-GNSS product GBM (Deng et al. 2016) and comparing PPP with and without identified problematic satellites. Finally, the conclusions are summarized.

Extended PPP model for quality monitoring
The raw observations of pseudorange and carrier phase can be expressed as: where the superscript i represents the satellite number, and the subscripts j and f represent the station and frequency, respectively; ρ i j denotes the geometric distance between satellite i and station j; c denotes the speed of light in vacuum; Δt i and Δt i correspond to the receiver clock offset and satellite clock offset, respectively; m i j and T i j denote the tropospheric mapping function and Zenith Tropospheric Delay (ZTD), respectively; I i j;f denotes the ionospheric delay with frequency f ; b j;f and b s f denote the receiver-side and satellite-side pseudorange hardware delays, respectively; B j;f and B i f denote the receiver-side and satellite-side carrier-phase hardware delays, respectively; N i j;f denotes the integer carrier phase ambiguity; ε i P;j;f and ε i L;j;f represent the noise of pseudorange and carrier-phase observations, respectively. Other error items, such as phase wind-up (Wu et al. 1993), Phase Center Offset (PCO), Phase Center Variation (PCV) and relativistic effect, should be precisely corrected to obtain high-accuracy solution (Petit and Luzum 2010).
In this study, the ionosphere-free observations are used and satellite orbits and clocks are fixed to the precise products provided by the service system, therefore, the linearized ionosphere-free observation equations read as: where ΔL i j and ΔP i j denote the observed minus computed measurements for the ionosphere-free pseudorange and carrier-phase, respectively; u i j is the unit vector from the satellite i to the receiver j; Δx j and cΔt j denote the increments with respect to a priori station position vector and receiver clock offset, respectively; � N i j is the new ambiguity parameters in unit of cycle contaminated by both the pseudorange and carrier-phase hardware delays and λ the wavelength of the ionosphere-free carrier-phase observations; ε i P;j and ε i L;j represent the noise of ionosphere-free pseudorange and carrierphase observations, respectively.
To deal with the orbit and/or clock products with possible biases for satellite i, the functional model of PPP is extended with additional orbit and clock parameters and can be written in the following matrix form:  where e i is a unit vector for the i th component, i.e. all the elements are zero except that the i th component equals to 1; R rsw is the rotation matrix that converts coordinate in the radial, along and cross directions to the Earth-Centered, Earth-Fixed (ECEF) system; Δr i , Δa i and Δc i are the orbital biases in the radial, along and cross directions, respectively. c � Δt i is the bias in the satellite clock product in the unit of meter.
Because of the high correlation between the orbital radial bias and satellite clock bias, Equation (3) can be further simplified as: where Δr i contains orbital radial biases and clock biases. A monitoring station network with well-distributed stations is employed for the estimation of the product biases as well as the station parameters for quality monitoring purpose. In principle, the estimability of the parameters depends on the coverage of the monitoring network, in other words, the larger the geographical coverage of the monitoring network, the better the estimability. Although the radial orbital error dominates the line-of-sight ranging error, the effects of the other two directions are also not negligible according to the contribution factors of radial, along-track and crosstrack errors to observations which are 0.98, 0.14 and 0.14, respectively (Heng et al. 2011;Montenbruck, Steigenberger, and Hauschild 2015). Furthermore, the along and cross biases can hardly be compensated by the estimated satellite clock parameter in a large or global network (Lou et al. 2014), as they are mapped differently into observations at different stations. That is why along and cross orbit biases must be monitored for RTPPS to reduced their impact on the PPP solution.
For the discussion of the estimation at a single epoch, the observation equations of the receiver j expressed as Equation (4) is, for brevity, written as where l j is the observation vector; A j is the design matrix of the station parameters; x j is the vector of the station parameters i.e. station coordinates, receiver clock offset, tropospheric wet delay and ambiguity parameters; C i j is the design matrix of the additional orbit and clock parameters for satellite i; b i is the vector of the orbit and clock parameters for satellite i.
Involving all the stations for the estimation, the observation model with additional product bias parameters for satellite i can be written as: . . . where the subscript (1, 2, . . ., n) indicates the index of the station. v is residual vector. Assuming uncorrelated observations, the individual weight matrices P nn are diagonal matrices and consequently the entire weight matrix P is also a diagonal matrix. It is worthwhile to point out that only the parameters representing orbital biases are common for different tracking ground stations. From Equation (6), the observation contribution to the Least Square (LS) adjustment can be easily derived and the contribution of the state vector can be easily considered by accumulating the weight of the station parameters P x j to current epoch: The normal equation for estimating b i can be derived by eliminating the station parameters and written as: with Equations (8) and (9) indicate that the components of normal equation for estimating b i can be calculated individually for each station and then accumulated together for the final estimation. In addition, the satellite product biases can be estimated more precisely using observations over a certain time period. We assume that the product biases can be modeled by a random walk process as: where t denotes the epoch number; w t is the process noise; q is the power spectral density of the process noise which should be selected fine-tuned according to the estimated time series. The state equation is introduced into the LS estimation as the following pseudo observation equation: In case that the residual v is larger than 3 times of the ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi D v ð Þ p , which means that an abrupt orbit or clock jump may occur, the D v ð Þ can be set to a larger value, e.g. 3600 and the product bias can be re-estimated.
When the product bias of satellite i is estimated, the calculation of the corresponding IURE for station j can be written as: where IURE est represents the IURE calculated using the estimated product bias. In addition, the IURE calculated using the product bias bi , which is generated from the comparison between the final product and the real-time product, can be written as: IURE ref can be calculated when the final product is available and is used for the validation of the real-time quality monitoring results. The IURE values differ for different tracking stations since the projection coefficients, i.e. À u i j , differ. A quality indicator QI i is defined to represent the maximum absolute IURE est for satellite i among its tracking stations and its expression is as follow: where S i is a set which contains the indexes of all tracking stations for satellite i. To make sure that the QI i is larger than the magnitude of actual IURE, additional 0.03 m, which is the SISRE value of this realtime GPS satellite product, is added. In case the QI i exceeds the pre-defined threshold δ, an alarm will be triggered and the corresponding satellite observations would be removed from the data processing at the user-end. The four possible situations of the quality monitoring for satellite i are summarized in Table 1. On one hand, the first line of Table 1 indicates that when the maximum absolute IURE ref is less than or equal to the δ, the product of satellite i is normal. Meanwhile, if the QI i is less than or equal to the δ, the quality monitoring result is correct. But if the QI i is larger than the δ, it means a false alarm occurs. On the other hand, the second line indicates that when the maximum IURE ref is larger than the δ, the product of satellite i is biased. If the QI i is larger than the δ, an alarm is triggered correctly. But if the QI i is less than or equal to the δ, it means a missed detection occurs. The determination of the threshold δ will be discussed in Section 3.

QC in PPP
As is well known, QC is a key component in GNSS data processing to detect observation outliers, especially for real-time precise positioning. In this study, the output of PPP QC is used as basic information for detecting possible problematic satellites according to their frequentness marked with outliers by the stations. There are several approaches of QC, in this study the well-known Baarda method is employed and will be introduced in this subsection for completeness. The standard Gauss-Markov model can be written as (Yang et al. 2013): where σ 0 is the theoretical reference standard deviation and set to 1 in this study; Q and P are the cofactor and weight matrix, respectively. If any misspecification occurs, such as outliers, the model must be adapted to consider such misspecifications properly. If there is an outlier in the k th observation, the extended functional model takes the form (Teunissen 2000 where c k is a m-dimensional unit vector with the k-th element equals to 1 and b k represents the bias. The significance of b k can be determined by the wellknown w-test (Baarda 1968) where Q vv is the cofactor matrix of the residuals. A robust estimation method is employed for PPP QC. A robust equivalent weight matrix � P k is used to adapt the weight of observations whose w-test value exceeds a threshold and its expression can be written as (Yang and Xu 2016): with k 0 ¼ 1:5 The satellite product biases can also be identified as outliers by PPP QC. For a single station, it is impossible to distinguish whether an identified outlier is from the observation itself or from the satellite orbit and clock products. Therefore, a voting method is used for the preliminary detection of the satellite product biases. Specifically, it is very likely that the outliers come from its orbit and clock product if the observations from one satellite are identified as outliers by several stations. It should be noted that once the w-test of the satellite exceeds k 0 it means this satellite is identified as outlier. This voting method will also be described in Section 2.3. Figure 1 shows the flow chart of the processing procedure for the system-level quality monitoring. First, PPP data processing is performed in parallel for all the stations by using a monitoring network with station coordinates fixed to IGS weekly solutions. The aforesaid QC approach is implemented in PPP to identify observation outliers and the detected observations will be down-weighted to eliminate their impacts on the PPP solutions. Then, the problematic satellite will be voted out according to the observation residuals from all the tracking stations. Specifically, if the observations of satellite i are marked as outliers by most of its tracking stations, it is reasonable to assume that its orbit and clock products may contain biases. In this study, when one satellite is detected by PPP QC in more than 30% of its tracking stations, it will be marked as a problematic satellite. All the marked satellites will be sorted according to this percentage and the satellite with the highest percentage is the most likely to have problems and the estimation of its orbit and clock biases will be given very high priority. Before calculating the matrix P kk for estimating the product biases, the original observation weights of the problematic satellites need to be restored since they are likely to be down-weighted by PPP QC. Thereafter, with the C matrix for specified problematic satellite(s), the normal equation for solving the product biases is constructed with P kk and l k from all the stations. Beware that the matrix P kk is accumulated over time which is denoted by using covariance matrix of the station parameters of the last epoch. The estimation of the satellite product biases can be achieved via Equation (8). Meanwhile, the QC should be employed again after the estimation of the satellite product biases in order to detect and exclude outliers. The quality indicator can be generated using the estimated product biases. If the quality indicator is smaller than the threshold δ, the satellite will be excluded from the problematic satellites and the voting step will be restarted to find out another possible problematic satellite. Otherwise, the estimate product biases are used to update the observation residuals and this satellite would be kept in the estimator and estimated together with the next possible problematic satellite. This procedure will run iteratively until no satellite is selected by the voting step.

Procedure of the quality monitoring
Finally, the quality alarm information for the problematic satellites will be sent to users. The same procedure will be restarted at the next epoch.

Processing strategy and statistics of 20-day solution
The proposed quality monitoring method is combined with the GFZ PPP software according to the processing procedure in Section 2.3 for the monitoring of real-time orbit and clock products. The real-time products generated at the GFZ IGS RTAC are selected as monitoring targets in this experimental evaluation. As one of the RTACs of the IGS real-time service, GFZ has been providing real-time products for years using its homemade software package (Ge et al. 2012). Recently, the real-time products of all four global GNSS systems including BDS-3 are generated and contributed. The orbits are processed in a batch mode and updated every 2 h using about 100 IGS stations, while the clocks are updated every 5 s using about 85 IGS real-time stations (Zuo et al. 2021). The monitoring network is shown in Figure 2, where data streams of all the stations are available in real-time via the internet from IGS Real-Time Service (RTS) caster. We intentionally select the monitoring stations that are not involved in the precise clock estimation. However, in Africa, South America and some marine areas, some overlapping stations are included for the quality monitoring as there are almost no redundant real-time stations available.
Some details about the data and network, as well as the parameters of the observation model and data processing strategies used in our experiments are summarized in Table 2.
In order to determine the threshold δ beyond which an alarm should be triggered, the impacts of different user range errors caused by product biases on the positioning are investigated for each satellite at each station by adding a simulated error to the observed minus computed measurement Ñl, and the expression is shown below: Where Ñl i k is the measurement vector with all members zero except for the members corresponding to the i-th satellite which is set to be the simulated error, i.e. the value in the horizontal axis of Figure 3, Ñx i k is the positioning bias caused by the product bias of satellite i, the subscript, superscript and other symbols have the same meaning with Equation (6) and (7). The statistical results in Figure 3 show that the positioning bias increases with the magnitude of the user range errors. When the magnitude of the user range error reaches 0.2 m, the median values of both the horizontal and vertical positioning biases are larger than 0.05 m and the maximum values are about 0.15 and 0.25 m, respectively. When the magnitude of the product bias is 0.15 m, most of the horizontal and vertical positioning biases are within 0.10 and 0.15 m, respectively. Considering the impacts of the different product biases and the sensitivity of the quality monitoring method, δ is set to be 0.2 m, which means a satellite product is seen as normal when the its product bias is less than 0.2 m.
Before discussing the statistic results of the quality monitoring system, we take a look at an example of IURE est and IURE ref values as well as the corresponding quality indicators shown in Figure 4. When the satellite product can pass the detection of PPP QC for at least 70% of its tracking stations, the product bias of this satellite will not be estimated and their values are all zeros. Therefore, the quality indicator of this satellite is set to be 0.03 m according to Equation (14). From Figure 4 we can find that before about 13:00 UTC the satellite product biases of G15 are not estimated and their IURE est values are all zeros and the corresponding quality indicators are 0.03 m. Meanwhile, the IURE ref values are smaller than 0.1 m, which means the satellite product is normal. After 13:00 UTC, the IURE ref values of different stations get larger, which means the satellite products are biased. Meanwhile, the IURE est values of different stations get larger and are close to the IURE ref , which means the product bias estimation is correct. In addition, the quality indicators based on the IURE est can cover all the IURE ref values when the satellite products are biased and thus the alarms can be triggered effectively. Figure 5 shows the 20-day quality monitoring results. One point corresponds to one GPS observation, whose abscissa value represents its absolute IURE ref value and the ordinate value represents the quality indicator generated by the quality monitoring system. About 98.5% points located in the square area    in the lower left corner, which means that both the satellite products and the monitoring results are normal. The points located in the square area in the upper right corner mean that the satellite products are biased and the alarms are triggered. In addition, because the quality indicators are larger than or equal to the absolute IURE ref values, the points mainly locate in the upper triangle area. The points located in the rectangular area in the lower right corner mean that the satellite products are biased but no alarms are triggered. These points are classified as missed detections and the missed detection rate is 4 � 10 À 6 for this experiment. is larger than the δ, these points are classified as correct alarms according to the definition in Table 1. Otherwise, these points are classified as false alarms and the false alarm rate is 1:2 � 10 À 5 for this experiment. Furthermore, about 91.7% of satellite products can pass the detection of PPP QC for at least 70% of their tracking stations and the ordinate values of the corresponding points are 0.03 m. The TTA of the quality monitoring system is shown in Figure 6. The TTA is defined as the time elapsed from a quality event occurring to the alarm received by the users. It consists of three parts: 1) observation data and satellite products collection, 2) PPP data processing and product bias estimation, 3) uploading the quality information and sending to users. The three parts can run in parallel and the TTA values are mainly within 5 s and the mean value is 4.47 s as shown in Figure 6.

An example of successful detection
For better understanding of the algorithm of the proposed method, we picked up a typical identified event for detailed investigation, for which large biases occur for both orbit and clock products.
In the product comparison, since the most radial biases can be assimilated by the clock biases (Lou et al. 2014;Montenbruck, Steigenberger, and Hauschild 2015), the radial biases should be eliminated from the clocks before comparison. It should be noted that the good clock products should be stable over the time and the systematic bias can be assimilated by the carrier-phase ambiguity and thus hardly affect the positioning. From the product comparison result, G10 and G32 contain large biases, the differences are shown in Figure 7 with other three satellites with good agreement as example. The orbits of G10 and G32 have large jumps at 20:40 Coordinated Universal Time (UTC). From 20:40 to 23:40 UTC, the biases for the orbit of G10 range from −0.37 to 0.1 m, 1.1 to 1.6 m and −0.8 to 0.3 m for the radial, along and cross directions, respectively, and the clock bias has become a linear drift since 20:40 UTC and reaches to 0.47 m at 23:40 UTC. From 20:40 to 23:40 UTC, the biases for the orbit of G32 range from −0.38 to −0.18 m, 2.2 to 2.3 m and −1.3 to 1.1 m for the radial, along and cross directions, respectively, and the clock bias has become a linear drift since 20:40 UTC and reaches to 0.49 m at 23:35 UTC. After 23:40 UTC, the orbit biases of G10 and G32 become normal and are smaller than 5 cm, 10 and 10 cm for the radial, along and cross directions, respectively, the clock bias of G10 keeps stable with value of 0.47 m, while the clock bias of G32 has a jump at 23:40 UTC and reaches to 0.81 m.
The quality monitoring procedure has detected the abnormal of products of G10 and G32 and the estimated satellite product biases are shown in Figure 8. The estimated radial biases for G10 have become larger since 20:37 UTC and reach 0.48 m at 23:40 UTC and    then keep stable. The estimated along and cross biases for G10 have opposite large jumps with value of about 1 m at 20:37 UTC, and from 20:37 to 23:40 UTC the values of along and cross biases range from 1.08 to 1.75 m and from −1.01 to 0.24 m, respectively, and after 23:40 UTC the along and cross biases both become normal and are smaller than 10 cm. The estimated radial biases for G32 have become larger since 20:36 UTC and reach to 0.51 m at 23:35 and afterward increase to 0.90 m and keep stable. In addition, the estimated along and cross biases for G32 have opposite large jumps with values of about 2.2 and 1.09 m at 20:37 UTC, and afterward the biases range between 2.0 and 2.5 m for the along component and between −1.0 and −1.5 m for the cross component. The IURE est generated using the estimated product biases are shown in the top part of Figure 9. The red triangles represent the positive and negative quality indicators calculated using Equation (14). We can find that the quality indicators can cover all the IURE ref value in the bottom part of Figure 9, which means that when the IURE ref exceeds the δ, an alarm can be triggered instantly.

PPP validation
We further use PPP and some IGS reference stations to validate the contribution of the quality monitoring for the case shown in Section 4 by comparing the positioning performances of PPP solutions with and without the detected problematic satellites. The raw kinematic PPP results using GPS observations for station OUS2 and MAC1 are shown in Figure 10. We can see that the positioning results have become worse since 20:37 UTC. The maximum errors for the horizontal and vertical directions reach up to 0.26 m and 0.41 m, respectively. Such large positioning errors are not allowed for high accuracy applications.

Data availability statement
The real-time satellite orbit and clock products used in this study are available from caster: products.igs-ip.net:2101 with mountpoint SSRA00GFZ0.