An effective automatic processing engine for improving the multi‑GNSS constellation precise orbit prediction

Orbit prediction (OP) recently tends to be a very crucial step for supporting real-time GNSS orbit services due to the dynamic stability of navigation satellite orbits. The OP performance depends on the length of the predicted orbits and the accuracy of precise orbit determination (POD) as basis. Considering this, a new automatic processing engine is established for improving the multiple global navigation satellite systems (multi-GNSS) constellation OP performance. From the architecture-oriented high-performance parallel processing perspective, the multi-node and multi-core computer sources are fully exploited to implement the hourly update of the current multi-GNSS POD. For MEO-type satellites (e.g., Galileo satellites), the accuracy of predicted orbits is improved from 3.8 cm, 6.5 cm, and 12.3 cm to 3.5 cm, 4.3 cm, and 6.3 cm, in the radial, cross, and along directions, respectively, compared to the three-hour POD update. Despite the shortened OP length, the OP performance of regional navigation satellite system (RNSS) satellites is still limited due to their regional observability. The BDS-IGSO and QZSS-IGSO satellitesexhibit radial directional orbital errors of up to 36.9 cm and 28.9 cm, respectively. Therefore, an orbit fitting (OF) processing method with orbit reconstruction is implemented into the processing engine. By utilizing this method, the radial orbital errors for BDS-IGSO and QZSS-IGSO satellites can be reduced to 7.0 cm and 10.4 cm, respectively. The mean real-time positioning errors are thus reduced from 28.3 to 18.4 cm and from 24.4 to 18.2 cm in the horizontal and vertical components, respectively


Introduction
The positioning service is implemented by synchronously processing the distance from receiver to four or more satellites based on the space resection method when the satellite positions and clock offsets are known.Therefore, the errors of satellite orbits and clock offsets are crucial for the performance of the Positioning Navigation Timing (PNT) services.Currently, there are four global navigation satellite systems: the US GPS, Russian GLONASS, and the newly developed European Galileo system and the third generation Chinese BeiDou System (BDS-3).All the four systems are now in full operation and provide their strong space-based PNT services.This brings great opportunities and challenges for more precise and reliable GNSS applications (Prange et al. 2017).Besides, there are also regional navigation satellite systems (RNSSs), for example, the Indian NavIC and the Japanese Quasi-Zenith Satellite System (QZSS).In 2022, QZSS with a constellation including one geostationary earth orbit (GEO) and three inclined geosynchronous orbit (IGSO) satellites is available for enhancing PNT services in the Asia-Pacific region.The BDS is developed from a regional system to the current global system BDS-3, and a number of GEO and IGSO satellites are still involved in the constellation for regional augmentation.
The precise orbit determination (POD) accuracy for medium and high orbit satellites does not only depend on the observation quality, observation modeling, and satellite orbit modeling, but it is also strongly related to the distribution of ground stations, the length of observing time and so forth.The orbit prediction (OP) performance is determined by the length of the prediction and the accuracy and length of the estimated orbits.These aspects should be well investigated for the multi-GNSS OP processing.
A high accuracy of real-time orbits is usually guaranteed through the rapid orbit prediction algorithm with the highprecision dynamic force model of GNSS satellites.The IGS ultra-rapid (IGU) orbit products are presently generated by using the fast-update mode.The current IGU orbit products are updated every six hours (https:// igs.org/ produ cts/).The three-dimensional root-mean-square (RMS) values can reach up to about 3 cm and 5 cm for the estimated orbits and the predicted orbits (combined orbits for real-time services), respectively (Lou et al. 2022).The real-time satellite orbit accuracy can be further improved due to the shortened orbit prediction arc.The update interval of the IGU orbit products was accordingly reduced from 12 to 6 h.Accordingly, the orbit prediction arc has been further shortened when the update interval is decreased to 3 h or even 1 h by several research institutions through the improvement in ultra-rapid POD strategies (Liu 2016;Deng et al. 2016).Currently, the procedure of parameter estimation becomes time-consuming in the multi-GNSS POD when dealing with a massive number of unknown parameters within a single estimator.To address this challenge, researchers have focused on simplifying the algorithms and adapting the models to reduce the number of parameters and observations without significant loss of accuracy (Zhao et al. 2018;Li et al. 2019).However, these improvements in the processing efficiency are limited because a certain number of parameters and observations should be kept to ensure the accuracy of POD and OP.Hence, researchers have started to shift their focus toward using high-performance computing techniques for large GNSS network solutions (Jiang et al. 2019;Cui et al. 2021).Using the existing multi-core and multi-node computing platforms, a more efficient and accurate processing engine is needed to realize the fast update of the current multi-GNSS POD.
To guarantee the estimated orbit accuracy, various orbit fitting (OF) processing methods have been used to improve OP performance.A normal equation combination of recent short multiple arcs was employed to continuously and rapidly update the initial orbit state parameters associated with the latest sliding short arc (Lou 2008).Results show that the three-dimensional RMS value for predicted GPS orbits is about 5 cm when the update rate of the sliding orbit arcs is one hour.The impact of different fitting arc lengths of estimated orbits on the GPS orbit prediction was evaluated for IGU orbit products.Choi et al. (2013) found that the optimal arc length of estimated orbits is 40-45 h, which can yield the most stable and accurate orbit predictions during 2010.Using different fitting arc lengths of estimated orbits, the predicted orbits are generated for the satellite clock estimation and precise point positioning ambiguity resolution (PPP-AR).Results confirm that the 42-h arc length of estimated orbits is optimal for the satellite clock estimation and PPP-AR (Li et al. 2015).The global PNT services have already benefited from the improvements in MEO satellite OP.However, the global distribution of ground stations for GEO and IGSO satellites is not as uniform as that for medium Earth orbit (MEO) satellites because of their regional coverage, which degrades the OP performance for a constellation with GEO and IGSO satellites.To address this issue due to regional observability, the OF processing should be specially investigated to improve the OP performance for RNSS.
This study implements the hourly update of the five system GPS + GLONASS + Galileo + BDS + QZSS POD using a new system-specific parallel processing strategy.To further improve the RNSS OP performance, we develop a new OF processing method based on the hourly updated estimated orbits.First, the parallel processing engine is established based on multi-core and multi-node computer sources for realizing the fast update of multi-GNSS POD.Then, the OF processing method with orbit reconstruction is further implemented into the processing engine to enhance the RNSS OP performance.Third, the validations of this newly developed engine are carried out to demonstrate the feasibility of the fast update of multi-GNSS POD and the improvements in the RNSS OP performance due to the OF processing method.Finally, the conclusions and contributions are summarized for this study.

Parallel processing for multi-GNSS POD
We have introduced a new parallel processing capability to our multi-GNSS POD processing engine.This advancement is accompanied by a detailed description of the hardware and software infrastructure at GFZ.In addition, we provide an in-depth explanation of the underlying principles of our novel parallel processing method to enhance comprehension.

Hardware and software
For parallel processing of multi-GNSS POD, three exclusive servers 'srtg4,' 'sigs52,' and 'srtg3' are used as hardware devices from the real-time GNSS group at German Research Centre for Geosciences (GFZ).The CPU configuration information of these hardware resources is given in Table 1.The servers are used for realizing the multi-core parallelization in the following POD processing due to their multi-core architectures.
The Positioning And Navigation Data Analyst (PANDA) software was designed and developed at Wuhan University (Liu and Ge 2003;Shi et al. 2012;Shi et al. 2013).A branch based on the PANDA 2008 version has been further developed at GFZ, mainly by several PhD students from China for their PhD studies.The new software package is named as PANDA-G.The current PANDA-G software can support both the standard IGS data processing and the real-time multi-GNSS precise positioning services.Furthermore, it also supports the low earth orbit (LEO) POD, the GNSS/ LEO integrated POD and RNSS POD.The accuracy of the PANDA-G products is comparable to that of the IGS final orbit products and LEO precise orbit products published by the European Space Agency (ESA) and the NASA's Jet Propulsion Laboratory (JPL).Therefore, all the data processing and validation in this study can be carried out with the PANDA-G software package for enhancing orbit prediction.

Data and processing strategy
To complete the task, this study requires an automatic POD and OP processing engine based on the PANDA-G software package.The whole processing is implemented in a sliding window mode with a certain length of the window (session/ arc length) and a specified moving step (update interval).For example, using a 24-h window and a moving step of one hour means the single session processing is performed every hour using data of the latest 24 h.
Apart from the processing line, a line to collect observation files of the tracking network and a line to retrieve the navigation messages are needed as well.The implementation of the processing engine is shown in Fig. 1.For the data collection for multi-GNSS data processing, the hourly data files required for the processing are downloaded from all the multi-GNSS experiment (MGEX) ftp servers (Montenbruck et al. 2017;Johnston et al. 2017).The multi-thread technique is used for downloading, merging hourly FTP files, and data preprocessing, which significantly reduces the computational time.Table 1 shows that the hardware cannot allow more than one thread to run on each core due to the processor without hyper threading.Therefore, the multi-thread parallel processing can be categorized into the multi-core parallelization.
Although various multi-thread techniques are introduced into the processing system, the minimum update interval of three hours can be kept only through computer resources of a single server when the number of ground stations exceeds 100 in the POD processing for the current multi-GNSS constellation with about 120 satellites.To implement a faster update of POD processing and thus shorten the orbit prediction arc length, a new system-specific parallel processing is further developed for the multi-GNSS POD, as shown in Fig. 1.Taking the common GPS system as a baseline to keep the consistency of inter-system biases (ISBs), the number of other satellites can be properly adjusted among the three servers listed in Table 1 in order to evenly divide the computational burden for parallel computing.
For each session POD processing, the global distribution of IGS/MGEX stations is divided into multiple regions.Within each region, station priorities are initially determined based on their historical long-term data availability and quality.Then, the quality of stations is assessed through data preprocessing during the current POD processing.As a result, stations exhibiting good quality and higher priority are selected for inclusion in this study.The one-way float ambiguities estimated in the LSQ procedure can be used for establishing the wide-lane and narrow-lane double differenced (DD) ambiguities for ambiguity fixing.The information about fixed ambiguities will be recorded and imposed as constraints in the LSQ estimation, i.e., the fixed solution.
For the three session PODs, the normal equations (NEQs) derived from the fixed solutions are extracted for the NEQ combination.The NEQ combination is realized based on a Fig. 1 Implementation of the automatic processing engine to generate predicted orbits based on the PANDA-G software package.For the three session PODs, the same processing procedures are carried out in parallel and enclosed in a red frame.The arc length of estimated orbits and the sampling interval are set to 24 h and 300 s, respectively parallel algorithm of parameter elimination.The observation equations for the LSQ procedure are denoted as follows with the global parameters y and the time-dependent parameters x(i) at the epoch i.
where x j (i) is the state transition matrix of x from epoch i-1 to epoch i in the state equation.The design matrices C j (i) and D j (i) are employed to establish the equation for observation l j .The contribution of new observations and state equations at epoch i to the NEQ can be written as follows: where the NEQ at epoch i is expressed as follows after the above contribution is introduced.
Following the parameter elimination principle (Chen et al. 2022), the time-dependent parameters x can be removed from the NEQ, while their effects are corrected as follows to generate the NEQ where only the global parameters y are kept for the final NEQ combination.
where the parallel computing of (4) has been implemented based on multi-core computing sources.The NEQs from the three session PODs are processed through (4) and then are stacked together as a combined NEQ in (5).
To solve the combined NEQ, we involve the global parameters y, i.e., the earth rotation parameters (ERP), the initial orbit parameters including satellite positions, velocities, and force model parameters, ISBs, and ground station coordinates, respectively.After the orbit combination is finished, the satellite orbits and ERP will be updated again for the (1) x(i) (5) following OF and OP processing.Using this system-specific strategy, the computer resources in the three GFZ servers are fully exploited in parallel to guarantee the hourly update of the POD processing for the current multi-GNSS constellation.
In 2014, GFZ started contributing the GFZ multi-GNSS rapid product GBM to the IGS MGEX project based on the Earth Parameter and Orbit System (EPOS) software.After a few years of development, the GBM product includes all the GNSS, i.e., GPS, GLONASS, Galileo, BDS-2, BDS-3, and QZSS.The GBM is well known in the GNSS world for its high quality (Deng et al. 2016;Montenbruck et al. 2017;Li et al. 2019;Kawate et al. 2023).Besides the daily rapid GBM product, we have generated hourly ultra-rapid precise orbit and clock product GBU for real-time GNSS applications since 2015.More than 40 global research institutes and companies have registered to access the GBU products during the period from 2015 to 2019.The GBM/GBU orbit products are used as a reference for the feasibility test of the strategies for enhancing multi-GNSS OP in this study.For comparison, the observation models, force models, and estimated parameters in a single session of POD processing align with those for GBM/GBU products and are listed in Table 2.However, it should be noted that several processing items, including the arc length of estimated orbits, the Solar Radiation Pressure (SRP) model, or a more up-to-date tropospheric model etc., can be optimized in future studies.

Orbit fitting method for enhancing OP
Although the hourly update is completed, the orbit quality of RNSS satellites is still limited due to its regional coverage.To enhance the OP performance of RNSS satellites, we develop a new OF method with orbit reconstruction and introduce its principle in detail.

Multi-GNSS tracking network
Figure 2 shows the global distributions of IGS ground stations dedicated to the multi-GNSS constellations.For MEO-type satellites, good POD accuracy can be guaranteed through the IGS tracking network of uniformly distributed ground stations, as shown in Fig. 2.However, for RNSS satellites, the distribution of tracking ground stations is rather worse due to their regional coverage characteristics.From Fig. 2, it can be seen that the regional observability of GEO satellites is limited to the Asia-Pacific area.Although the IGSO tracking network is slightly extended from the Asia-Pacific area, its unobservable areas still include the North American east and south regions, most of South America, and the west and north parts of Africa, respectively.The POD accuracy for GEO and IGSO satellites is accordingly limited due to the regional observability, which has been confirmed in Section 'Validation.'Accordingly, for GEO and IGSO satellites with poor observability, the improvements in their OP performance are investigated by shortening the prediction orbit arc length and enhancing the accuracy of estimated orbits available for OP, respectively.The former can be implemented through the parallel processing strategy of multi-GNSS POD mentioned in the last section.The latter one is realized by using a new OF processing method, which will be described in detail within the following section.

OF processing method
To reduce the effects due to the regional observability, the estimated orbits with a certain arc length will be fully exploited in the OF processing procedure.As is well known, the middle part for a whole arc of estimated orbits has a higher accuracy than its boundary parts in onesession POD.Considering this fact, the Center for Orbit Determination in Europe (CODE) IGS analysis center extracts a one-day middle orbital arc in the three-day orbit solution to generate the daily rapid orbit products (COD) in a routine manner.The root-mean-square (RMS) values of orbit differences between the COD and GBM rapid orbit products are listed in Table 3 for different parts of estimated orbital arcs.As shown in Table 3, a better agreement between the COD and GBM orbits is achieved when the orbital arc is closer to its middle part, which indicates a reduction of mis-modeling errors.Taking the QZSS-IGSO satellites as a typical example, the RMS of orbit differences can be reduced from 9.2 cm, 10.4 cm, and 7.6 cm to 2.1 cm, 8.0 cm, and 2.3 cm, in the radial, cross, and along directions, respectively, for the 1-h middle part of orbital arc compared to that for the whole orbital arc.The middle parts of estimated orbits should be fully exploited in the following OF/OP processing procedures, especially for the RNSS satellites.
After the hourly update of multi-GNSS POD is realized by using the parallel processing strategy, a 1-h middle part of estimated orbits can be extracted from a one-session POD for the subsequent OF processing.Figure 3 shows the OF processing method for enhancing the OP performance of multi-GNSS satellites.The whole processing is implemented in a sliding window mode.The 1-h middle parts from estimated orbits in the previous windows are combined together for the reconstruction of the current window, which is called as orbit reconstruction.The arc length of the reconstructed window can be flexibly extended by increasing the number of previous windows involved for the orbit reconstruction.After the orbit reconstruction is finished with a certain  Fig. 3 Orbit fitting processing method for enhancing the orbit prediction of multi-GNSS satellites arc length, the OF processing procedure can be carried out based on the reconstructed orbits for re-estimation of initial orbital state parameters.Considering the orbit discrepancies of up to 5 mm at the boundary of connecting hourly sessions, a random-walk constraint is introduced into the orbital state parameters at the boundary in the OF processing and its pseudo-observation at the boundary (epoch k) is expressed as follows: where indices m represent the position of the parameter x in NEQ.p x m is the weight of constraint on the parameter with a priori unit weight variance 2 0 of 1.The time interval Δt and the priori power density of process noise p are 300 s and 1 × 10 -6 , respectively.The contribution of this constraint to the NEQ is given as follows: where the estimated initial orbital state parameters are then employed to generate the enhanced predicted orbits through obit integrator.

Validation
To demonstrate the effectiveness of the automatic processing engine, we mainly focus on two major issues related to the OP performance.On the one hand, the feasibility of the system-specific parallel processing strategy is confirmed at first and its impact on the OP performance will be evaluated as well.On the other hand, we will test the presented OF processing method in order to further improve the OP performance, especially for RNSS satellites.

Fast update of multi-GNSS POD
For one session, POD operated on a single computing server, introducing all the observations into the parameter estimation is still time-consuming due to limited computer sources, which cannot satisfy the hourly update of multi-GNSS POD.To release the computational burden, the processing engine is extended to multiple computing servers running in parallel by using a system-specific parallel processing strategy.This can guarantee the fast update of multi-GNSS POD.Unlike the time-consuming POD where all the observations are directly involved in a whole NEQ for parameter solution, the NEQs extracted from three session PODs in parallel are combined in the final solution for the new parallel processing strategy.In view of this, we should ensure that the combination of NEQ for the final parameter solution will not cause a significant degradation in orbit accuracy compared to the POD with a whole NEQ.
For comparison, the differences between the two POD processing strategies, referred to as S1 and S2, respectively, are listed in Table 4. Using the new parallel processing strategy, the hourly update of multi-GNSS POD can be realized in strategy S2 based on three Linux computer servers running in parallel.The mean and maximum computing time can be, respectively, reduced from 2.3 h to 45 min and from 2.6 h to 55 min, compared to S1.Among the sliding sessions of estimated orbits in S2, we can identify the session corresponding to that for S1 due to its faster update of POD.Then, the common sessions are compared, as shown in Fig. 4. Figure 4 shows the RMS values of differences between the multi-GNSS orbits estimated for strategy S1 and that for strategy S2, in the radial, cross, and along directions, respectively.For MEO-type satellite orbits, the averaged RMS values of orbit differences can reach up to 2 mm, 4 mm, and 5 mm in the radial, cross, and along directions, respectively.For IGSO satellite orbits, the orbit errors in the radial, cross, and along directions are about 3 mm, 5 mm, and 6 mm, respectively.Even for GEO satellites, a millimeter-level agreement is achieved after the combination NEQ is used for the final parameter solution.
Because of a faster update of POD, the orbit prediction arc derived from strategy S2 is accordingly reduced compared to that for strategy S1.To evaluate the improvements in real-time orbit accuracy due to the shortened orbit prediction arc, the predicted orbits generated using strategy S1 and strategy S2 are used for real-time services and compared with external GBM final satellite orbits.Figure 5 shows the RMS values of orbit differences with respect to GBM final orbits for real-time orbits derived from strategy S1 and S2 in the radial, cross, and along directions, respectively.Taking Galileo satellites as a typical example for MEO-type satellites, the orbit accuracy for S2 can be improved from 3.8 cm, 6.5 cm, and 12.3 cm to 3.5 cm, 4.3 cm, and 6.3 cm, in the radial, cross, and along directions, respectively, compared to that for S1.For IGSO and GEO satellites with poor observability, more significant improvements in orbit accuracy can be achieved by shortening the prediction arc.Despite the application of strategy S2, the orbit errors for RNSS satellites remain rather high, particularly in the radial component.For instance, the orbit errors in the radial component can reach up to 36.9 cm for BDS-IGSO satellites, which cannot serve well for a cm-level precise positioning.

Assessments of orbit fitting method
To further improve the OP performance of RNSS satellites, the OF processing method has been applied in the automatic processing engine.Figure 3 shows multiple hourly arcs over the middle of previous windows are used for orbit reconstruction.To confirm the positive contribution of this orbit reconstruction, we used the OF solution with a 24-h orbital arc length (called strategy S3) for comparison with strategy S2.The only difference between strategy S2 and S3 is the orbit reconstruction.Figure 8 shows the comparison of realtime orbit accuracy in the radial, cross, and along directions for strategy S2 and S3, respectively.For BDS-IGSO and BDS-GEO satellites, the averaged RMS values of real-time orbit differences in the radial, cross, and along directions are reduced by 2.2 cm, 1.7 cm, and 2.4 cm and by 2.1 cm, 5.7 cm, and 11.7 cm, respectively, by using the orbit reconstruction.The improvements due to the orbit reconstruction in the radial, cross, and along directions reach up to 1.8 cm, 1.2 cm, 0.8 cm and for QZSS-IGSO and 2.5 cm, 5.2 cm, and 0.7 cm for QZSS-GEO satellites, respectively.Figure 8 is included in 'Appendix.'Results have confirmed that the advantage of this orbit reconstruction for the OF processing with a 24-h orbital arc length in terms of the OP performance.In the OF processing method, the orbital arc length can be flexibly extended by increasing the number of pieces of previous windows for orbit reconstruction.With an increase in the orbital arc length, there are two possible scenarios that may occur during the OF&OP processing.In the first scenario, insufficient training data may be utilized due to a short orbital arc length.In the second scenario, the force model may not be complex enough to accurately capture the relationship between the training data and output parameters because the orbital arc length is too long.These two situations will cause underfitting phenomena for the OF&OP procedure.Therefore, a proper arc length is required in the OF processing in order to fully optimize the OP performance.
For comparison, the orbital arc length is set to 24 h (strategy S3), 48 h (strategy S4), and 72 h (strategy S5) in the OF processing, respectively.Figure 6 shows the accuracy of realtime orbits for strategies S3, S4, and S5, in the radial, cross, and along directions, respectively.As shown in Fig. 6, significant improvements can be achieved for RNSS satellites when the orbital arc length is increased from 24 h (S3) to 48 h (S4).For BDS-IGSO satellites, the averaged RMS values of real-time orbit differences are improved from 34.7 cm, 21.3 cm, and 19.9 cm to 7.0 cm, 9.6 cm, and 8.6 cm, in the radial, cross, and along directions, respectively.For QZSS-IGSO satellites, the orbit errors are reduced from 27.1 cm, 20.2 cm, and 14.6 cm to 10.4 cm, 15.2 cm, and 11.6 cm in the radial, cross, and along directions, respectively, due to the extension of orbital arc length.Results also confirm that underfitting problems occur when the orbital arc length is increased to 72 h (S5) in comparison with strategy S4.Using BDS-IGSO satellites as example, the real-time orbit accuracy is degraded by 1.3 cm, 1.4 cm, and 1.1 cm, in the radial, cross, and along directions, respectively, compared to that for S4.
Compared with S2 (without OF), the orbit accuracy for RNSS satellites has significantly improved due to S4 (with  satellite clock errors can be reduced from 1.10 to 0.22 ns due to S4 compared to S2, as shown in Table 5. Using the real-time products derived from S2 and S4, a real-time BDS-only precise point positioning (PPP) float solution filter is performed every two hours for the five real-time ground stations built for the HI-POS project, as shown in Fig. 7.There is a total of 251 sets of PPP results.To evaluate the PPP convergence performance at the 30-min filtering time, the absolute values of positioning errors for all 251 sets are calculated and arranged in ascending order.The values corresponding to the 99.7%, 95.5%, and 68.3% positions in this sorted list represent the positioning errors at the respective confidence levels.The horizontal positioning errors at the confidence levels of 99.7%, 95.5%, and 68.3% can be reduced from 327.6 cm, 92.0 cm, and 27.8 to 69.6 cm, 59.4 cm, and 23.7 cm due to S4 compared to S2, respectively.The mean positioning accuracy in the horizontal and vertical components is improved from 28.3 to 18.4 cm and from 24.4 to 18.2 cm, respectively.

Conclusions
In this contribution, an effective processing engine is established for improving the OP performance of the multi-GNSS constellation satellites.To improve the computational efficiency of GNSS POD and realize its hourly update, a new system-specific parallel processing strategy is developed based on the multi-node and multi-core computer sources.Using the parallel processing strategy, the update interval for multi-GNSS POD is reduced from 3 to 1 h, which can shorten the arc length for orbit prediction in each one-session POD.Because of the shortened orbit prediction arc length, the real-time orbit accuracy can be improved from 3.8 cm, 6.5 cm, and 12.3 cm to 3.5 cm, 4.3 cm, and 6.3 cm, in the radial, cross, and along directions, respectively, for MEO-type satellites (e.g., Galileo satellites).Although this parallel processing strategy reduces the update interval, the RNSS satellite still has rather high orbital errors due to its poor observability, especially for the radial component.For BDS-IGSO satellites, the orbital errors in the radial direction can reach up to 36.9 cm, which can degrade its performance for a cm-level precise positioning.To further improve the OP performance for the RNSS satellites, an OF processing method with orbit reconstruction is applied into the processing engine.
Using the OF processing method with an arc length of 48 h, the radial orbital errors can be reduced to 7.0 cm and 10.4 cm for BDS-IGSO satellites and QZSS-IGSO satellites, respectively.Due to this OF method, the mean positioning errors for BDS-only real-time PPP can be improved from 28.3 to 18.4 cm and from 24.4 to 18.2 cm in the horizontal and vertical components.

Fig. 2
Fig. 2 Global distributions of GNSS/RNSS tracking ground stations.The ground stations are shown as red dots for an individual constellation tracking network

Fig. 5
Fig. 5 RMS values of differences over April 2023 between real-time orbits and GBM final orbits for strategy S1 and S2, in the radial (top), cross (middle), and along (bottom) directions.Mean RMS values of orbit differences are shown in the legend box

Fig. 7 Fig. 8
Fig. 7 Distribution of real-time ground stations (top) and real-time BDS-only PPP positioning errors (bottom) in April 2023 for strategy S2 and S4

Table 2
Processing strategy, observation model, force model, and estimated parameters for the multi-GNSS POD

Table 4
Comparison of the two processing strategies for fast update of multi-GNSS POD

Table 5
Averaged STD values (unit: ns)of clock differences over April 2023 between GBM final clocks and real-time clocks for S2 and S4, respectively