Performance Analysis of Multi-GNSS Real-Time PPP-AR Positioning Considering SSR Delay

: For real-time (RT) precise point positioning (PPP), the state space representation (SSR) information is often delayed due to possible communication delays and specific broadcast intervals. In this case, the positioning results will diverge and re-converge due to the increase of SSR products extrapolation errors. In addition, RT orbit and clock offset accuracy, as well as their extrapolation errors, will vary in different systems and satellites. We propose a PPP with ambiguity resolution (PPP-AR) method that combines a time-differenced carrier phase (TDCP) model, in which the characteristics of the orbit and clock are considered. Under normal communication, the PPP-AR solution is obtained by fixing satellites with small SSR product errors. When the communication is abnormal, the TDCP model is utilized to extrapolate user coordinates by considering different extrapolation error characteristics of satellites. The experimental results show that GPS and Galileo SSR products have better accuracy than BDS, with signal-in-space user ranger errors (SISREs) of 2.7, 2.2, and 8.6 cm, respectively. Optimizing the PPP stochastic model based on SISREs can effectively reduce the convergence time. Under 5 min SSR delay, SISREs caused by clock and orbit extrapolation for GPS/Galileo/BDS are 3.5, 1.4, and 2.6 cm, respectively. After optimizing the TDCP stochastic model based on extrapolation errors, the horizontal and vertical positioning accuracies can be maintained at 0.7 cm and 5.0 cm. For multi-GNSS, the combination of the TDCP and PPP-AR can overcome the influence of short delay. After optimizing the stochastic model, the GPS/Galileo/BDS positioning accuracy can be maintained at about 2.4 cm under 3 min delay, showing an accuracy improvement rate of 59.3% compared with the traditional method using only PPP. Additionally, the rapid PPP convergence results can be obtained by inheriting previous filter state information when the communication recovers normally.


Introduction
The global satellite navigation system (GNSS) can provide high-quality positioning, navigation, and timing services for users around the world, with unique advantages in many fields of geoscientific research and applications [1,2].Precise point positioning (PPP) technology has gradually developed into a research front in the field of satellite navigation and positioning due to its own advantages and characteristics, demonstrating broad application prospects.In real-time (RT) positioning, PPP with ambiguity resolution (PPP-AR) highly relies on the support of external organizations, which provide precise products or broadcast correction information [3,4].Regardless of the communication mode, RT-PPP generally requires the user to receive correction information from the server continuously in real time, imposing rigorous requirements on the accuracy and timeliness of provided products [5,6].In order to address the high precision and real-time application needs of the global satellite navigation system, the international GNSS service (IGS) established the real-time working group in 2001.Subsequently, in 2007, the organization officially started the real-time service (RTS) to provide RT corrections based on state space representation (SSR).Indeed, corrections from the server are not broadcast all the times but are updated at a specific interval.SSR corrections are typically broadcast at intervals of 5 s or more.This implies that the correction information received by the client is not exactly "real-time" but rather possesses a certain data age.Moreover, there will be a certain delay in the corrections generated from the server sent to the user due to possible communication instability.The data age of SSR incorporates both the effects of delay and update frequency.In this case, the delay includes communication, computation, and broadcasting delay.The primary challenge in RT-PPP lies in maintaining the continuity and real-time performance of positioning, which is adversely affected by delays in receiving precise orbit and clock corrections.It is very common for synchronous precise orbit and clock offset corrections to be unattainable due to unstable communication links in practical applications [7,8].
The large SSR delay can directly lead to jumping and the re-convergence of PPP-AR results, failing to meet the demands of highly time-efficient and accurate applications.For this reason, in most studies, the effectiveness of RT-PPP results is preserved by forecasting predictions of precise orbit and clock offsets when communication delays occur.Yang et al. [9] proposed satellite clock prediction for different update intervals and applied it to RT-PPP.If a sudden communication break takes place, the interruption period for receiving orbit and clock corrections may extend from a few minutes to hours.El-Mowafy et al. [10,11] adopted prediction products with a specific method based on short-and long-term interruptions to address disruptions of precise orbit and clock corrections in RT-PPP applications.The three-dimensional (3D) accuracy of the method is maintained within a decimeter when such a break occurs.At the same time, positioning results in different scenarios also entirely proved the effectiveness of this method.Nie et al. [12] constructed a clock prediction model that eliminates the need to store received RTS clock offset corrections.This model facilitates clock offset prediction even with just one epoch of RTS data, addressing interruptions in the reception of RTS products.Different from the previous method of directly predicting orbits and clocks, Zhao et al. [13] proposed an alternative method that utilizes broadcast ephemeris to generate comprehensive corrections including orbit errors, satellite clock errors, and receiver correlation errors.Additionally, PPP-AR can also provide solutions by extrapolation using previous SSR information during communication interruptions.But, it is prone to the re-convergence phenomenon for data recovery, which affects the high precision and continuity of real-time positioning.
Therefore, it is important to overcome the problem of abnormal results caused by SSR delay in order to promote the application of RT-PPP.The time-differenced carrier phase (TDCP) method proposed by related scholars can overcome such problems to some extent.The TDCP is a method for accurate velocity estimation, while also obtaining accurate relative positional information during the calculation process [14][15][16].Similar to PPP, the time-differenced carrier phase (TDCP) is a method for accurate position estimation using only one receiver and can theoretically be combined with PPP.It can be used directly to estimate precise velocity without fixing the ambiguity, demonstrating nice positioning performance.Some studies have applied it to high-precision real-time earthquake monitoring [17,18].The PPP method, based on the mixed use of time-differenced and undifferenced carrier phase observations, can effectively improve positioning accuracy and computational efficiency, which has significant advantages in high-speed positioning [19].Using the TDCP model to extrapolate the receiver's position during SSR delay can protect the PPP filter from large orbit and clock errors.Consequently, PPP can converge rapidly by inheriting the previous state information when communication returns to normal.
Along with the development of a new generation of navigation satellite systems, such as BDS and Galileo, multi-GNSS has become the predominant trend in current high-precision positioning.Combined observations significantly increase the number of visible satellites, enhance space geometry configurations, and provide valuable advantages for high-precision positioning.Nonetheless, it is important to note that the accuracy of orbits and clocks for different satellites in multi-GNSS systems can vary significantly, which introduces complexities and challenges in error modeling and the reduction in positioning errors.For RT PPP-AR, the positioning performance highly relies on precise product quality as well as extrapolation errors under SSR delay [20].
As a result, this study proposes the PPP-AR method that takes into account differences in the performance of RT products to further enhance positioning reliability.In this study, the PPP-AR method considering the accuracy and extrapolation errors of precise products was introduced first.When SSR delay occurs, the method uses the TDCP model for positioning.Then, the quality of RT multi-GNSS precise products from the Centre National D'Etudes Spatiales (CNES) was evaluated, and the stochastic model was optimized based on signal-in-space user ranger errors (SISREs).Finally, by utilizing high-frequency observation data from the Multi-GNSS Experiment (MGEX), it has been verified that the proposed method is effective in avoiding abnormal positioning solutions and improving the positioning accuracy when the SSR delay occurs.

Methodology
In the data processing of real-time PPP-AR, orbit, and clock errors are important factors affecting positioning performance.Once the communication link is blocked, the accuracy of the calculated orbit and clock will degrade, and positioning results will diverge.In this study, we propose an RT multi-GNSS PPP-AR method that takes into account SSR delay, in which the characteristics of the orbit and clock themselves and their extrapolation errors are carefully considered.

General Idea on PPP-AR Considering SSR Delay
The general thought of this method is shown in Figure 1.When the communication is normal, the receiver position is solved by the ionosphere-free (IF) PPP-AR model, and the stochastic model is optimized by considering the accuracy of multi-GNSS orbit and clock products.When the SSR delay occurs, the TDCP model is used to acquire relative coordinates, and the stochastic model is optimized based on extrapolation errors of multi-GNSS products.Note that the SSR corrections are updated within specific intervals.Even if there is no communication delay, the orbit and clock will be frequently extrapolated with the fresh SSR corrections.As a result, the TDCP model will be started only when the time delay exceeds a certain threshold value.Also, it should be noted that the accuracy of coordinates derived by both the TDCP and PPP models will decrease with an increase in communication latency.However, using the TDCP model can protect the PPP model from large orbit and clock errors.Inheriting the previous filter state information, such as troposphere and ambiguity parameters, enables the PPP filter to recover and converge rapidly when communication becomes normal.

Multi-GNSS PPP Aided with the TDCP Model
The dual-frequency IF combination is one of the commonly used PPP models [21].This model combines dual-frequency pseudo-ranges and carrier phases, respectively, to eliminate the impact of the ionosphere.The IF pseudo-range and carrier phase measurements can be expressed as follows: where  and  denote pseudo-range and carrier phase observations, respectively; superscripts  and  refer to the GNSS system and the satellite pseudo-random noise (PRN), respectively; subscript  refers to the receiver;  is the frequency,  and  are different frequency points;   , denotes the geometric distance between the satellite and receiver;  is the speed of light;    and  , are the satellite and receiver clock offset;   , is the orbit error;   , is the tropospheric delay;  ,  ,   , ,  ,  , and   , are the uncalibrated pseudo-range and phase delays for the receiver and satellite, respectively;   is the IF combination wavelength;  , , is the integer phase ambiguity; � , , � and � , , � are the sum of other unmodeled errors and noise for pseudo-range and carrier phase measurements.This observation equation is also applicable to multi-GNSS PPP.At this point, the number of receiver clock offset parameters is the number of participating systems [22].The satellite clock offset and pseudo-range hardware delays can be corrected with the IGS clock products.The receiver clock offset and ambiguity can be re-parameterized as follows: The slant troposphere wet delay for all satellites can be represented by one zenith troposphere parameter with mapping coefficients.The parameters to be estimated in the equation are as follows: where  is the receiver coordinates;  ̂  is the receiver clock offset;  , is the zenith tropospheric wet delay; and  , , * represents the ambiguity (assuming m satellites are observed).The IF ambiguity is commonly divided into wide-lane and narrow-lane ambiguities, each of which is individually fixed to restore the integer properties of the IF ambiguity.
When an SSR delay occurs, the receiver coordinates (  ,   ,   ) at the current epoch can be obtained by following equations: where �  0 ,   0 ,   0 � denotes the last absolute coordinates obtained by PPP at the moment before the SSR delay occurs; the term Δ  * represents the relative position derived from the TDCP model.In Equation ( 1), the carrier phase TDCP model between two epochs  0 and  can be expressed as follows: where Δ represents the differencing operation.For example, Δ  , =   , () −   , ( 0 ) represents the change in geometric distance between two epochs, and other terms in Equation ( 5) are defined accordingly.The ionospheric delay has been eliminated through the IF combination.In the absence of cycle jumps, ambiguities are also eliminated, along with the majority of common model errors that change gradually within a limited sampling rate [23].The unified expression of this model is as follows: where Δ primarily represents receiver noise and minimal uncanceled errors.The linearization of the equation above can be written as follows: where [      1] is the matrix of coefficients.Δ  , Δ  , Δ  , and Δ   represent the change components of the receiver position and clock between two epochs.

Methods for Evaluating Orbit and Clock Errors
Multi-GNSS positioning can further increase the number of visible satellites and optimize the geometric configuration of constellations [24,25].In the conventional stochastic model of the multi-GNSS PPP and TDCP, satellite elevation angles are utilized to characterize the random properties of observations.The altitude angle stochastic model represents a functional model that defines the weight-ratio relationship among observations, considering the variations in observation quality with the satellite altitude.The typical expression for the stochastic model based on satellite elevation angles is as follows [26]: where   2 represents receiver-side observation noise;  represents the satellite elevation angle; and  0 2 denotes the reference variance of observations.The PPP performance of the terminal is closely associated with the accuracy of the clock and orbit products broadcasted by the server.Due to factors such as satellite constellation characteristics, solar elevation angles, and orbit types, there are differences in the clock and orbit products and their extrapolation errors among different systems [27,28].By addressing errors in satellite clocks and orbits, optimizing the stochastic model can significantly reduce conventional errors and enhance convergence speed [29].The SISRE is not only the primary parameter describing spatial signal accuracy but also a critical indicator that affects GNSS positioning services' performance.This index reflects the combined impact of orbit and clock errors within the broadcast ephemeris on users ranging under specific data age conditions.It provides a more comprehensive assessment of system performance [30,31].We can optimize the stochastic models of the PPP and TDCP based on the calculated SISRE, which can be expressed as follows: where   2 represents receiver-side observation noise;   , 2 signifies the SISRE; and superscripts  and  refer to the GNSS system and the satellite PRN, respectively.During normal communication, the PPP model is used, and the stochastic model can be optimized based on the SISRE value of precise products themselves.The term   , 2 specifically represents the SISRE derived from orbit and clock offset errors.The TDCP performance is primarily influenced by the geometric distribution of satellites and the accuracy in extrapolating real-time precise ephemeris data [32,33].When SSR messages are delayed, the TDCP model is used, introducing additional extrapolation errors into the calculations of orbit and clock offset.The stochastic model can be optimized based on the SISRE value for extrapolation errors.The term   , 2 represents the SISRE derived from orbit and clock extrapolation during SSR delay.
The clock-only and orbit-only contributions to SISREs can be expressed as follows [34,35]: where   denotes the clock offset error and   ,   ,   represent orbit errors in the radial, along-track, and cross-track directions, respectively.The calculation coefficients  and  depend on the altitude of the satellite, with values for different systems provided in Table 1.Taking into account the influence of both clock offset and orbit errors, the SISRE can be expressed as follows: Table 1.SISRE calculation coefficients for the statistical contribution of radial, along-track, and cross-track errors to the line-of-sight ranging error [32].The root-mean-square (RMS) errors in Equations ( 10) and ( 11) can be derived as follows:

BDS (GEO/IGSO
where  = , , , ;   (  ) denotes the clock offset or orbit error at time   ; and  denotes the number of samples.
is calculated differently for the two cases of normal communication and the occurrence of SSR delay in Equations ( 10) and (11).Taking clock offset as an example, under normal communication, the clock error is calculated following the restoration of the realtime clock offset with SSR corrections as follows: where () denotes the satellite clock offset at time ; superscripts  and  indicate a satellite and a reference satellite, respectively; and subscripts  and  are for realtime and post-precise products, respectively.When SSR messages are delayed, additional extrapolation errors are introduced into the clock offset calculation.The term   can be determined by the following equation: where  and  0 denote the current time and the time stamp of the newest SSR messages, respectively;  −  0 means the time delay of SSR messages;  , 0  () is the extrapolated clock offset at time , which is calculated by the newest SSR messages;  , 0  ( 0 ) is the clock offset at time  0 computed through the broadcast ephemeris and SSR corrections; and    () and    ( 0 ) represent the clock offsets of post-precise product at times  and  0 , respectively.

Data Collection and Processing Strategies
As depicted in Figure 2, observation data from 10 MGEX stations during DOY 121-127, 2022, were utilized to assess the positioning performance under SSR delay.All the receivers were equipped to receive dual-frequency signals from GPS, Galileo, and BDS, with a data sample interval of 1 s.RT precise orbit, clock, and observable-specific signal bias (OSB) products for the same period were also acquired from the CNES.The sample intervals for orbit, clock, and OSB data were 5 min, 5 s, and 30 s, respectively.To evaluate the SISRE of precise products, post-precise products from Wuhan University were used as the reference.In order to assess the effectiveness of the proposed model under varying delays, SSR delay was simulated from 1 s to 300 s in this study.For the PPP-AR model, the processing strategies and data settings are outlined in Table 2.The TDCP also uses the IF combination to eliminate the impact of ionospheric delay during the positioning process.

Accuracy Evaluation of Real-Time Orbit and Clock Products
In this section, the quality of CNES precise orbit and clock products was evaluated, and the extrapolation errors for GPS, BDS, and Galileo systems were analyzed by simulating real-time SSR delay.The calculated SISREs were used to optimize the stochastic model in the PPP and TDCP models, thereby improving positioning reliability.
During normal communication, the precision of RT products plays a pivotal role in influencing PPP-AR positioning performance.To evaluate the accuracy of CNES products, we used the post-precise products provided by Wuhan University as the reference.Figure 3 presents the average SISREs caused by clock and orbit errors for each satellite in three systems during DOY 121-127, 2022.GPS and Galileo RT products exhibit superior quality in comparison to BDS; the average SISREs for GPS, Galileo, and BDS systems are 2.7 cm, 2.2 cm, and 8.6 cm, respectively.Due to significant clock and orbit errors associated with GEO satellites, BDS GEO satellites were excluded from the calculation of average SISREs.It is evident that, except for the individual satellites of BDS, the precision of both clock offset and orbit among satellites within the same system is nearly identical.This consistency makes it possible to directly optimize the PPP stochastic model based on SISREs of three systems in the experiment.Based on Equation ( 14) presented in Section 2.3, clock and orbit extrapolation errors for GPS, Galileo, and BDS were processed under 5 min delays.Figures 4 and 5 display the SISRE resulting from clock and orbit extrapolation for each satellite in three systems, considering delays ranging from 1 to 300 s.As can be observed in Figures 4 and 5, the extrapolation accuracy of clock and orbit products decreases as the SSR delay increases.It is noteworthy that the extrapolation error growth is of obvious variations among different systems, and the growth rate and trend for various satellites also differ significantly.Therefore, it is necessary to optimize the TDCP stochastic model based on the SISRE of each satellite in subsequent experiments.We also can see that clock extrapolation errors for Galileo and BDS are smaller than those for GPS, while the orbit extrapolation errors for all three systems exhibit similar characteristics.Additionally, the extent of variation remains minimal within the range of 1 to 10 s for extrapolation errors.Given that the CNES updates RT product information every 5 s, an SSR delay within 10 s can be negligible in terms of its impact.Consequently, this study sets 10 s as the threshold for switching to the TDCP positioning mode.Figure 6 illustrates the average SISREs resulting from clock and orbit extrapolation with a 5 min delay for each satellite in the GPS, Galileo, and BDS systems during DOY 121-127, 2022.Galileo demonstrates better accuracy than the others.SISREs attributed to extrapolation within 5 min delay remain below 3 cm, except for E14 and E19.Different from the results caused by precise products themselves (Figure 3), GPS exhibits notably lower extrapolation accuracy compared to Galileo and BDS, with some satellites showing errors exceeding 4 cm.However, Block IIF-type satellites, including G01, G03, G06, G09, G10, G25, G26, G27, G30, and G32, consistently display smaller clock extrapolation errors due to their rubidium atomic clocks.Additionally, BDS-3 exhibits slightly higher extrapolation accuracy than BDS-2, likely attributed to the implementation of domestically produced rubidium and hydrogen atomic clocks, which enhance stability [37].The average SISREs resulting from clock and orbit extrapolation under 5 min delay for GPS, Galileo, and BDS are 3.5 cm, 1.4 cm, and 2.6 cm, respectively.

Performance Analysis of GPS/Galileo/BDS PPP-AR Considering SSR Delay
The combination of multiple GNSS systems offers substantial enhancements in time to first fixed (TTFF), AR positioning accuracy, and reliability [38,39].Consequently, the fusion of GPS, Galileo, and BDS was adopted to validate the performance of the new method.By integrating SISREs into the stochastic model under normal communication and SSR delay for optimization, the positioning accuracy before and after optimization was compared.Subsequently, optimized models were utilized to verify the positioning accuracy of the new method.
In the normal case of communication, the PPP model is used, and the stochastic model is optimized based on SISREs calculated from the quality of precise products themselves in Section 3.2.Figure 7 illustrates GPS/BDS/Galileo PPP-AR positioning errors using the original and optimized stochastic model at the DYNG station.When the original stochastic model is utilized for positioning, floating point issues arise intermittently, with more pronounced effects observed in the up direction.After using the optimized stochastic model, stable fixed results can be obtained during these periods, effectively alleviating the negative impact on positioning accuracy caused by the low accuracy of individual satellite orbits and clocks.GPS/BDS/Galileo PPP-AR was conducted on 10 stations, with re-initialization occurring every two hours.Table 3 presents the mean TTFF and positioning accuracy of fixed solutions for these specific stations.The results unequivocally demonstrate that adopting the optimized stochastic model accelerates the convergence speed of solutions.Particularly, the reduction in TTFF from 21.5 min to 13.5 min for the BIK0 station is notable (8 min less).Moreover, the optimized stochastic model enhances the positioning accuracy of each station to varying degrees.The varying degrees of improvement among different stations are associated with their respective observational environments, observed satellites, and satellite geometrical configurations.When the SSR delay exceeds a certain threshold, the TDCP model is used, and the stochastic model is optimized based on the extrapolated SISRE results in Section 3.2.Figure 8 illustrates the disparities in the RMS values in the horizontal and vertical directions when using the original and optimized stochastic models for the TDCP under various delay conditions.It becomes evident that an increase in SSR delay leads to a gradual reduction in the real-time accuracy of multi-GNSS.In Figure 8, it is apparent that when using the optimized stochastic model, the accuracy remains at approximately 0.7 cm horizontally and 5.0 cm vertically for delays up to 10 min.Compared to the original stochastic model, the accuracy has improved by 20% and 40% in the horizontal and vertical directions, respectively.To highlight the advantages of multi-GNSS PPP aided by the TDCP method, we simulated multiple SSR delays of 30 s, 2 min, and 3 min for 24 h, respectively.We then compared the positioning performance of this approach with the traditional PPP-AR method, both utilizing the optimized stochastic model.
Figures 9-11 illustrate the positioning errors of the traditional PPP and PPP/TDCP combined method under 30 s, 2 min, and 3 min SSR delays at the DYNG station.Since the experiment involved 24 h high-frequency data, we enlarged the results for double SSR delay epochs to facilitate a more comprehensive analysis of positioning stability during SSR delays and data transmission recovery.For a 30 s delay, the positioning results for both methods are generally comparable in most periods, with accuracy consistently maintained at a stable level.However, there are still some periods where the results of the traditional method exhibited dispersion phenomena (e.g., as shown in Figure 9, right).Both methods quickly achieve convergent positioning results after communication recovery.With a 2 min delay shown in Figure 10, the traditional method displays a sudden change in results during the delay period, with noticeable discontinuity.Using the TDCP model ensures positioning reliability during this period.In Figure 11, as the delay extends to 3 min, the new method consistently maintains better positioning performance, with errors of 0.8 cm and 2.7 cm in the horizontal and vertical directions, respectively, while the accuracy of the traditional method is up to 2.7 cm and 9.2 cm.The discontinuity of the articulation at the 2 min delay (as shown in Figure 10, right) compared to the 3 min delay may be related to the state of its inherited previous epoch parameters.With an increase in SSR delay, the new method achieves convergent results in a shorter time by utilizing previous filter state information during communication recovery compared to traditional PPP.The solutions of the traditional method experience varying degrees of accuracy loss in all three directions, with the pronounced dispersion up to 30 cm in some periods, while errors of the PPP/TDCP combined model can be maintained within 10 cm.The above findings demonstrate that the new method effectively overcomes abnormal issues during SSR delay, ensuring that accuracy remains within a reliable range.Simultaneously, data from 10 stations were processed.Figure 12 presents the average RMS values of MGEX stations using the traditional PPP and PPP/TDCP combined method for all delay periods.It can be observed that the most significant improvement is evident with a 3 min delay.The average 3D accuracy for 10 stations using the traditional and new methods is 5.9 cm and 2.4 cm, respectively.The 3D positioning accuracy of the DYNG station increased from 5.7 cm to 1.7 cm, with a 70.2% enhancement.However, variations in positioning performance among different stations are observed, potentially influenced by factors such as the number of visible satellites, the constellation geometry, and the observation environment around the station.

Conclusions
In real-time positioning, PPP-AR heavily relies on correction information broadcast by the server; however, it frequently encounters SSR delays caused by communication anomalies.This study introduced the multi-GNSS PPP-AR aided with the TDCP method, which accounts for SSR delay.To maintain positioning stability and reliability during SSR delay, the TDCP was utilized for position extrapolation.Stochastic models for the PPP and TDCP were optimized based on the quality and extrapolation errors of RT products, and the performance of PPP-AR under SSR delay in the new method was analyzed.The experimental findings are summarized as follows: (1) Both the accuracy of precise products themselves and their extrapolation errors vary widely across different systems.GPS and Galileo products demonstrate better accuracy with SISREs of 2.7 cm and 2.2 cm, which are smaller than that of BDS (8.6 cm).However, when the SSR delay occurs, the extrapolation errors for BDS and Galileo are smaller.The average SISREs from clock and orbit extrapolation for GPS, Galileo, and BDS are 3.5 cm, 1.4 cm, and 2.6 cm, respectively.(2) Optimizing the PPP stochastic model based on the accuracy of the products themselves can effectively shorten the PPP-AR convergence time and improve the fixed solution accuracy.By optimizing the TDCP stochastic model based on extrapolation errors, the accuracy under a 10 min delay is maintained at 0.7 cm and 5.0 cm in the horizontal and vertical directions, resulting in improvements of 61.1% and 37.5%, respectively.The multi-GNSS PPP-AR method with stochastic model optimization has the capability to acquire precise position information by utilizing the TDCP during SSR delay.With a 3 min SSR delay, the average 3D positioning accuracy of 10 stations can be improved from 5.9 cm to 2.4 cm using the TDCP, representing a 59.3% improvement.Furthermore, positioning results can swiftly converge after communication recovery.The method presented here effectively addresses the issue of anomalous positioning results caused by SSR delay due to unstable communication links in practical applications, thus avoiding re-convergence.

Figure 2 .
Figure 2. Distribution of 10 selected MGEX stations in the PPP-AR experiments.

Figure 3 .
Figure 3.The average SISREs caused by clock and orbit errors for GPS, Galileo, and BDS satellites during DOY 121-127, 2022.

Figure 4 .
Figure 4.The SISREs for clock extrapolating under 1~300 s delay for each satellite in the GPS, Galileo, and BDS systems on DOY 121, 2022 (the right graphs show the variation in extrapolation errors from 1 s to 10 s for the highlighted sections shown in the left graphs, in the same units as the left graphs).

Figure 5 .
Figure 5.The SISREs for orbit extrapolating under 1~300 s delay for each satellite in the GPS, Galileo, and BDS systems on DOY 121, 2022 (the right graphs show the variation in extrapolation errors from 1 s to 10 s for the highlighted sections shown in the left graphs, in the same units as the left graphs).

Figure 6 .
Figure 6.The average SISREs for clock and orbit extrapolating under 5 min delay for each satellite in the GPS, Galileo, and BDS systems during DOY 121-127, 2022.

Figure 7 .
Figure 7.The positioning errors for GPS/BDS/Galileo PPP-AR with the original and optimized stochastic model at DYNG on DOY 121, 2022.

Figure 8 .
Figure 8. Positioning RMS errors in horizontal (h) and vertical (v) directions for the GPS/Galileo/BDS TDCP using the original (origin) and optimized (opt) stochastic model under different delays.

Figure 9 .
Figure 9. Positioning errors at DYNG under multiple 30 s delays for 24 h on DOY 121, 2022, with different PPP approaches (the right graphs are the results for double SSR delay epochs, corresponding to the highlighted region in the left graphs, in the same units as the left graphs).

Figure 10 .
Figure 10.Positioning errors at DYNG under multiple 2 min delays for 24 h on DOY 121, 2022, with different PPP approaches (the right graphs are the results for double SSR delay epochs, corresponding to the highlighted region in the left graphs, in the same units as the left graphs).

Figure 11 .
Figure 11.Positioning errors at DYNG under multiple 3 min delays for 24 h on DOY 121, 2022, with different PPP approaches (the right graphs are the results for double SSR delay epochs, corresponding to the highlighted region in the left graphs, in the same units as the left graphs).
Data processing flowchart of the real-time PPP-AR method considering SSR delay.

Table 2 .
Processing strategies and data settings for PPP-AR.

Table 3 .
TTFF and positioning accuracy for GPS/BDS/Galileo PPP-AR using the original and optimized stochastic model under no SSR delay.