High-precision GPS orbit determination by integrating the measurements from regional ground stations and LEO onboard receivers

High-precision Global Navigation Satellite System (GNSS) orbit and clock products are crucial for precise applications. An evenly distributed global network enables continuous tracking for GNSS satellites, while a regional network may result in tracking gaps in the areas where monitoring stations are not deployed. This also means that the orbit determination accuracy based on a regional network is not comparable to that with a global network. Integrating the measurements from regional ground stations and Low Earth Orbit (LEO) satellites onboard receivers is a potential approach for generating GNSS orbit and clock products with centimeter-level accuracy, which is particularly important for BDS and the local commercial providers relying on a regional network. In the integrated Precise Orbit Determination (POD), LEO satellites are used to compensate for the drawback of regional ground stations in the precise orbit and clock determination of GNSS satellites. To validate the role of LEO satellites in the orbit determination with a regional network, 6 International GNSS Service stations around China and 13 LEO satellites from January 20 to 26, 2019, including GRACE-C/D, SWARM-A/B/C, Jason-3, Sentinel-3A/B, and SAT-A/B/C/D/E are selected in this study to perform the integrated POD. The orbit and clock accuracies of GPS and LEO satellites are evaluated by comparison with precise products. The average Root Mean Square (RMS)of GPS orbit errors in the radial (R), along-track (T) and cross-track (N) directions are 2.27 cm, 3.45 cm, and 3.08 cm, respectively, and the clock accuracy is better than 0.15 ns based on a comparison with the final products provided by Center for Orbit Determination in Europe (CODE). The LEO orbit accuracy is better than 2 cm in the R direction, and the position errors are mostly within 4 cm. The results indicate that the integrated POD can generate high-precision orbit and clock products for GPS and LEO satellites based on regional network stations. Finally, the integrated POD products are assessed for Precise Point Positioning (PPP). Simulated kinematic PPP has a comparable performance in terms of the convergence time and positioning accuracy. With more LEO satellites available, the orbit and clock determination accuracy and PPP positioning accuracy can be improved.


Introduction
GNSS plays an extremely important role in Positioning, Navigation, and Timing (PNT), and brings new opportunities for high-precision orbit determination of LEO satellites.With the successful launch of TOPEX/ POSEIDON, the era of utilizing onboard GNSS data for LEO satellite orbit determination has officially begun (Tapley et al., 1994).In this POD approach, precise orbit, and clock products of GNSS satellites are derived from ground station measurements.Subsequently, these precise products are utilized to conduct LEO POD by using onboard GNSS data.It is evident that the orbit determination accuracy for LEO satellites directly depends on the accuracy of GNSS satellite products.At present, the IGS Analysis Centers (ACs) utilize the data from many global ground stations to generate the realtime, ultra-rapid, rapid, and final products (https:// igs.org/ produ cts/# about) for GPS, GLONASS, Galileo, and BDS.The POD accuracy of GNSS satellites is significantly influenced by the quantity and distribution of participating ground stations.Although well-distributed global ground stations are accessible to the scientific community, some commercial operators or BDS (Tang et al., 2018;Yang et al., 2020) choose to establish their own regional monitoring stations for the sake of system stability and reliability.Based on the regional network, it is impossible to generate GNSS satellite orbits and clock offsets with high accuracy, leading to the difficulty in providing high-precision services with GNSS.
In recent years, numerous LEO constellation projects have been rapidly developed, including Iridium (Fossa et al., 1998), StarLink (Khalife et al., 2021;Osoro & Oughton., 2021), OneWeb (Osoro & Oughton, 2021), and Boeing (Reid et al., 2016(Reid et al., , 2018) ) in the USA, as well as CentiSpace (Yang, 2019) in China.These LEO constellations can provide communication services through Inter-Satellite Links (ISLs) while also offering PNT services by broadcasting ranging signals, message data and augmentation information to various users to meet the demands of precision, availability, integrity, and reliability (Yang et al., 2024).LEO satellites are typically equipped with onboard receivers and can function as mobile monitoring stations for achieving global tracking of GNSS satellites.This brings about new opportunities for enhancing GNSS satellite orbit determination accuracy by processing the observations from both LEO satellites and ground stations.Many studies have been performed using this approach, also named as "one-step" method (Zhu et al., 2004).This method has been adopted to improve the accuracy of the estimated orbits and Earth Orientation Parameters (EOP) (Geng et al., 2008;Huang et al., 2020;Hugentobler et al., 2005;König et al., 2005;Li et al., 2018;Zhao et al., 2017).Furthermore, the one-step method has also been proven effective in improving the stability of parameter estimations (Boomkamp & Dow, 2005), the rate of ambiguity fixing (Zoulida et al., 2016), and the accuracy of geocentric coordinates (Männel & Rothacher, 2017).Additionally, it can enhance data processing efficiency in various scientific and practical applications.However, in the existing studies, the accuracy of satellite clock offsets estimated with the one-step method has not been considered, and its performance in providing high-precision navigation services has also not been mentioned.
Given this background, we conduct the integrated POD processing experiments using the data from regional ground stations and LEO satellites.It is crucial to investigate the achievable accuracy of orbit and clock products with the integrated processing method.Furthermore, we need to analyze whether the products obtained with this approach can meet the demands of high-precision services.Although the monitoring station deployment of BDS is regional, it seems more meaningful to conduct integrated POD using BDS and LEO satellites.However, due to the data limitations, we currently cannot obtain redundant onboard BDS observation data from enough LEO satellites.Therefore, we take GPS satellites as an example, the integrated POD is conducted along with an accuracy assessment.The integrated POD strategy, which combines a regional network and 13 LEO satellites from January 20 to 26, 2019, is first introduced in Section "Strategy of integrated POD".Then, the accuracy of the derived GPS orbits, clock offsets and LEO orbits are assessed in section "Results and analysis".Besides, the GPS orbit and clock products are also validated with PPP using ground stations worldwide.A discussion and the conclusions are given in Section "Discussion and Conclusions".

Strategy of integrated POD
This section is divided into three parts.The first part briefly introduces the methodology of the integrated POD.The second part gives a detailed description of the experimental data used in this study, mainly the data from LEO satellites and ground stations.Finally, the models adopted and the parameter estimation strategy employed in the integrated POD are listed in detail.

Methodology
In the integrated POD, the Ionosphere-Free (IF) combinations of pseudo-range and phase observations of GNSS satellites tracked by ground stations and LEO satellites are used.The observation equations for a ground station and LEO satellite can be expressed as follows: where the subscript G represents ground stations, and L represents LEO satellites, the superscript S represents a specific GPS satellite.P s G,IF , L s G,IF , P s L,IF and L s L,IF are the pseudo-range and carrier phase observations received by (1) Then, based on the principle of statistical orbit determination, the observation model at time t i can be expressed as follows (Montenbruck et al., 2002): where X G,i and X L,i are the orbital parameters of GNSS and LEO satellites, respectively, and X STA,i represent the parameters related to the ground station, such as the station coordinates and tropospheric delay.P G,i and P L,i represent other parameters, such as clock offsets, ambiguities, solar radiation coefficients, atmospheric drag coefficients, and empirical parameters.ε STA,i and ε L,i represent the noise of ground station data and LEO onboard data, respectively.
According to the equation of motion and variational equations of GNSS and LEO satellites, their state transition matrices ψ G (t i , t 0 ) and ψ L (t i , t 0 ) can be obtained by numerical integration.They meet the following conditions: where x G,i and x L,i are the state vectors of GNSS and LEO satellites, respectively, x G,0 and x L,0 are the initial epoch states.Based on Eqs. ( 2) and (3), orbital and other parameters of the GNSS and LEO satellites can be solved by using the least squares estimation.

Experimental setup
For BDS or some commercial service providers who cannot control a global network, it is a good means to generate precise orbits by integrating the measurements from regional ground stations and LEO satellites.LEO satellites are equipped with onboard GNSS receivers, which can be regarded as a supplement to regional ground stations to achieve full-arc tracking of (2) GNSS satellites.Therefore, in this section, we select 6 regional stations around China and the onboard data of 13 LEO satellites from January 20 to 26, 2019 to simulate the scenario of integrated POD.

LEO satellites
In recent years, an increasing number of LEO satellites have played important roles in the fields of gravity recovery, ocean altimetry, and atmospheric sounding.Considering the limitations of data access, the onboard GPS data of 13 LEO satellites are selected to conduct integrated POD processing.These LEO satellites are divided into five different missions.GRACE Follow-On is the continuation of GRACE's legacy tracking of Earth's water movement across the planet.Monitoring the changes in ice sheets and glaciers, underground water storage, the amount of water in large lakes and rivers, and the changes in sea level provides a unique view of Earth's climate and has far-reaching benefits for people (https:// grace fo.jpl.nasa.gov/).
Swarm is the ESA's first mission for Earth Observation (EO).It consists of three identical satellites named Alpha, Bravo, and Charlie (A, B and C).Swarm is dedicated to creating a detailed survey of the Earth's geomagnetic field and its temporal evolution as well as the electric field in the atmosphere by using a satellite constellation that carries sophisticated magnetometers and other instruments (https:// earth.esa.int/ eogat eway/ missi ons/ swarm).
The Jason-3 satellite is the follow-up altimetry mission of Jason-2/OSTM.Its objective is to provide the unique accuracy and coverage of the TOPEX/Poseidon, Jason-1, and OSTM/Jason-2 missions in support of operational applications related to extreme weather events, operational oceanography, and climate applications and forecasting.
Sentinel-3 is a multi-instrument mission that measures sea-surface topography, sea-and land-surface temperature, ocean color and land color with high accuracy and reliability.The mission supports ocean forecasting systems, as well as environmental and climate monitoring (https:// www.esa.int/ Appli catio ns/ Obser ving_ the_ Earth/ Coper nicus/ The_ Senti nel_ missi ons).
SAT-A/B/C/D/E are Chinese satellites using GNSS Radio Occultation (GNSS-RO) to collect atmospheric data for weather prediction and for ionosphere, climate, and gravity research.These satellites are capable of profiling vertical parameters such as temperature, humidity, and pressure in the middle and lower latitudes globally.It provides high vertical resolution and accurate observations for numerical weather prediction.The detailed information is shown as follows (Table 1).

Regional ground stations
Some providers, including BDS and other commercial companies, are unable to autonomously control and maintain well-distributed global ground stations.Regional stations are more likely to be available.In this study, 6 IGS stations (Johnston et al., 2017) around China, namely, SHAO, BJFS, WUH2, URUM, LHAZ, and ULAB, are selected to simulate regional network tracking.Their distribution is shown in Fig. 1.

POD strategy
We use the software platform called ShangHai ORbit DEtermination (SHORDE) (He et al., 1988), which was developed by Shanghai Astronomical Observatory (SHAO), Chinese Academy of Sciences (CAS), to perform the integrated POD.SHORDE is capable of processing the data from various sources such as GNSS, Satellite Laser Ranging (SLR), Doppler Orbitography and Radio Positioning Intergrated by Satellite (DORIS), and ISLs to estimate the orbit and clock offsets of GNSS and LEO satellites, as well as performing PPP solution and other geodetic applications.For the configuration of POD strategy, the arc length is 24 h, the data sampling interval is 30 s, and the elevation cutoff is 5°.L1 and L2 dualfrequency IF combinations of pseudo-range and phase observations are used in the data processing.The detailed strategy is given in Tables 2 and 3. To correct for antenna Phase Center Offsets (PCOs) and Variations (PCVs), we utilize the absolute phase centers (Schmid et al., 2007) to obtain the corrections for GPS satellites and ground stations.PCOs are applied for LEO onboard receivers, while PCVs are disregarded in this study.PCOs of GRCC and GRCD are obtained from the products (VGN1B) published by JPL.SWMA, SWMB, and SWMC utilize the satellite-specific ionosphere-free combined PCOs offered by the European Space Agency (ESA) (Siemes, 2019).For JAS3, the PCOs are obtained from the prior values provided by Centre National d'Etudes Spatiales (CNES) (Couderc, 2015).SN3A and SN3B use the ground-calibrated reference point and PCOs (Fernández Martìn, 2016), and SATA, SATB, SATC, SATD and SATE use the prior values provided by satellite operators.Regarding the satellite attitude, the measured quaternion data provided by operators are used for GRCC, GRCD, SWMA,  SWMB, SWMC, JAS3, SN3A and SN3B.Meanwhile, the remaining 5 LEO satellites adopt the geocentric stabilization attitude mode.To ensure data quality, the Turboedit algorithm (Blewitt, 1990) is employed for cycle slip detection.A weight ratio of 1:10,000 is set for pseudo-range to phase observations.Additionally, a threshold of three times the mean square error is utilized to edit the residuals of orbit determination.Both GPS and LEO satellites experience perturbation forces, including conservative and nonconservative forces.As for conservative forces, the Earth's gravitational force is modeled using EGM2008 (Pavlis et al., 2012).Considering the magnitude of perturbations, low-order gravitational models are employed for GPS satellites, while high-order gravitational models are used for LEO satellites.IERS Conventions 2010 is adopted to describe solid tide and pole tide, as well as for the computation of relativity effects.FES 2004 (Lyard et al., 2006) model is used to account for ocean tides.However, accurate modeling of nonconservative forces such as atmospheric drag and solar radiation pressure remains challenging due to various factors like atmospheric density variations, solar activity, and spacecraft characteristics.To compensate for the model errors associated with nonconservative forces, it is necessary to estimate additional force model parameters, such as the solar radiation coefficient, drag scale coefficient, and empirical parameters.For the computation of solar radiation pressure, we utilize the 5-parameter ECOM model for GPS satellites and the box-wing model for LEO satellites.Atmospheric drag is modeled using the box-wing model, which considers the atmospheric density computed with DTM94 (Berger et al., 1998) for LEO satellites.In addition, a drag scale coefficient is estimated every 1.5 h to represent the variability of atmospheric drag more accurately.The perturbation of atmospheric drag on GPS satellites is negligible.Ambiguities are solved but not fixed to integers.The satellite clock parameters are treated as white noise,  2 and 3.

Visibility analysis
For the integrated POD, both LEO satellites and ground stations track the signals from GPS satellites.The number of GPS satellites tracked by a LEO satellite will affect the accuracy of determined LEO satellite orbits.The number of receivers tracking GPS satellites, regardless of the LEO on-board receivers or the ground station receivers, will affect the orbit accuracy of GPS satellites.Therefore, whether it is for GPS or LEO satellites, conducting satellite visibility analysis is essential to obtain high-precision orbit and clock products.

Number of GPS satellites tracked by a LEO onboard receiver
When using the onboard GPS data to perform LEO POD, a larger number of visible GPS satellites provides more redundant observations, ensuring greater stability in parameter estimation.We take the selected 13 LEO satellites from 00:00 to 24:00 on January 20, 2019 as an example, a comprehensive analysis of the minimum, maximum, average, and valid numbers of tracked GPS satellites is illustrated in Fig. 2. The satellites with dualfrequency observations are considered valid.From Fig. 2, the number of GPS satellites tracked is very similar for the same LEO satellite mission, mainly due to two reasons: (1) Satellites of the same LEO mission carry the same type of onboard receivers and have the same number of satellite tracking channels; (2) Satellites of the same mission also have similar orbital heights, inclinations, and orbital types.Therefore, for the same LEO satellite mission, the number of GPS satellites tracked by each individual LEO satellite within the mission is similar.From the satellite visibility results on January 20, 2019, a minimum of 5 satellites and a maximum of 10 satellites are tracked by GRCC and GRCD.On average, they track 8.7 and 8.9 GPS satellites, respectively, and the data efficiency reaches 100%.For SWMA, SWMB and SWMC, the number of tracked GPS satellites ranges from 6 to 8, with an average of 7.9 and the data efficiency of 97.5%.For JAS3 which circles in a near-polar orbit, the number of tracked satellites varies with the maximum of 12, the minimum of 2, the average of 9.7, and the data efficiency of 97.9%.For SN3A and SN3B, the number of tracked satellites ranges from 5 to 8 with an average of 7.7, and the data efficiency of 97.9%.For the remaining 5 satellites of SATA, SATB, SATC, SATD and SATE, the maximum number is 12, the minimum number is not less than 5, the average number is above 10, and the data efficiency of SATC is the lowest and reaches 91.8%.The maximum and average number of GPS satellites tracked of SATA, SATB, SATC, SATD, and SATE are equivalent to that of JAS3, but the data efficiency is slightly lower than that of JAS3.The average number of GPS satellites tracked by SWMA, SWMB, and SWMC are slightly larger than that of JAS3, but the data efficiency is almost the same.
The selected 13 LEO satellites are all equipped with on-board receivers, which can track 8 to 10 GPS satellites on average.Therefore, the onboard observations can serve as an important data source in the GPS satellites orbit determination, especially under the conditions of regional network deployment.

Number of Receivers tracking GPS satellites
The number of receivers tracking GPS satellites directly affects the orbit and clock determination accuracy of GPS satellites.Figure 3 takes G01 and G04 as an example to analyze the number of receivers on January 20, 2019.In Fig. 3, the ground stations denoted by a green star is abbreviated as "GRD", the LEO satellite denoted by a blue cross is abbreviated as "LEO", and the combination of ground station and LEO satellite denoted by a red linecircle is abbreviated as "GRD + LEO". Figure 3 illustrates that the 6 regional ground stations exhibit brief tracking arcs for G01 and G04 with prolonged periods during which the number of ground receivers drops to zero.Consequently, relying solely on regional ground receivers makes it impossible to estimate high-precision orbits for GPS satellites.After involving 13 LEO satellites, the receivers tracking satellite G04 are scarce for both ground stations and LEO satellites.Therefore, subsequent accuracy assessments will not consider satellite G04.Meanwhile, for satellite G01, the maximum number of tracking receivers rises from 5 to 13, and the proportion of receivers being zero has significantly decreased from 55% to 2%.It indicates that the issue of tracking gaps resulting from the limited coverage of 6 regional ground stations can be effectively resolved with the inclusion of LEO satellites.Typically, there is at least one receiver tracking G01 most time on January 20, 2019.However, there are still some epochs at which no receivers can track G01.The lack of continuous tracking may result in the initialization of all ambiguity parameters in the orbit determination process.Figure 4 gives the average number of receivers tracking the GPS satellites.
According to Fig. 4, the average number of receivers for satellite G04 remains unchanged before and after the inclusion of LEO satellites.This consistency arises due to the limited capability of the selected 13 LEO satellites to receive observations from satellite G04.Conversely, for the remaining 31 satellites, the average number surges from 2.3 to 5.9, marking a substantial increase of 61.4%.This underscores the notable enhancement in receiver visibility to GPS satellites facilitated by the integration of LEO satellites, consequently benefiting the orbit determination process for GPS satellites.

Accuracy of determined orbit and clock
When the integrated POD is completed, the final products released by CODE are selected as a reference to PRN Fig. 4 The average number of receivers tracking GPS satellites evaluate the orbit and clock accuracy of the GPS satellites.
For LEO satellites, the corresponding organizations of the former 8 satellites have released Post Science Orbits (PSOs).Therefore, we can compare the LEO orbit derived from the integrated POD with PSOs to evaluate the orbit accuracy.However, for SAT-A/B/C/D/E, due to the absence of official precise orbit products, this study evaluates the orbit accuracy of these 5 LEO satellites using two approaches.The first one is comparing overlapping orbits and the second one compares the orbits derived from the integrated POD with the orbits solved by introducing fixed COD final GPS orbit and clock products provided to IGS.To avoid the errors resulting from different spatial datum used, this study considered the Helmert 7-parameter transformation when conducting orbit accuracy assessments.

GPS satellite orbit accuracy
To evaluate the role of LEO satellites in GPS satellite orbit determination, we first conducted orbit determination based on the measurements from the regional network only.The primary distinction in the POD strategy for the regional network compared to Tables 2 and 3 is that the arc length is set to 72 h. Figure 5 shows the orbit comparisons in the R, T, N and 3-dimensional (3D) direction between the COD orbits and those of the other 31 GPS satellites excluding G04, calculated by the regional network POD from January 20 to 26, 2019.The average RMS value in the R, T, N, and 3D direction is 0.32 m, 1.07 m, 0.52 m, and 1.23 m, and orbit accuracy is relatively stable from day to day, indicating the appropriate orbit determination strategy and reliable observation data.Based on this, it clearly indicates that with the regional network POD, only meter-level precision orbit products can be obtained, which is obviously insufficient to meet the demand for high-precision navigation services.
Figure 6 shows the orbit comparisons in the R, T, N and 3D direction between the COD orbits and those orbits of the other 31 GPS satellites excluding G04 calculated by the integrated POD on January 20, 2019.The horizontal axis represents the PRN of GPS satellites, and the vertical axis represents the RMS of the orbit differences for each satellite.
From Fig. 6, the derived GPS satellite orbits achieve centimeter-level accuracy compared to the COD final products.Among the three components, the R direction exhibits the smallest RMS, averaging 2.32 cm.It's followed by the T direction with an average of 3.65 cm and the N direction with an average of 3.14 cm.The average RMS of 31 GPS satellites excluding G04 during January 20-26, 2019, are shown in Fig. 7.
In Fig. 7, we performed the integrated POD for 7 days, the average RMS in the R, T, and N direction is 2.27 cm, 3.45 cm, 3.08 cm, and the average RMS in the 3D direction is 5.15 cm.Compared with the orbits based on the regional network only in Fig. 5, the orbit accuracy in R, T, N, and 3D direction has been improved by 92.9%, 96.8%, 94.0%, and 95.8% respectively.The results demonstrate that the involution of LEO satellites can effectively compensate for the deficiency of insufficient tracking arcs, and the orbit errors are greatly reduced, which agree to the existing findings in Konig et al. (2005), Hugentobler et al. (2005), Geng et al. (2008), Zhao et al. (2017), Li et al. (2018), Huang et al. (2020).Li et al. (2022a) conducted GNSS orbit determination by processing more than 150 globally distributed Multi-GNSS Experiment (MGEX) stations.Among them, GPS orbits have the highest accuracy, with an average 3D RMS of 2.11 cm.Obviously, the integrated GPS orbit accuracy is slightly worse than the orbit accuracy generated with a large global network from Li et al. (2022a), which can be attributed to that the regional network and insensitivity of the inter-satellite link measurements to the orbit orientation angles account for the larger normal orbit error of GPS satellites.The measurements from LEO onboard receivers are in fact some kind of satelliteto-satellite ranging measurements, which are insensitive to any motion of the constellation as a rigid body and the orientation parameters including inclination and the right ascension of ascending node.Therefore, adding the LEO satellites does not give constraints in the normal direction as strong as in the radial and cross-track directions.The error in the normal direction for GPS orbits of the integrated POD products is larger than the IGS products due to the regional stations used.In addition, the problem of insufficient observations still exists, though 13 LEO satellites are involved, because there are some epochs when no receivers track GPS satellites.Moreover, the ambiguities solved in this contribution are float solution.All these factors affect the accuracy of GPS orbit determination.However, the accuracy is indeed at the same level, approximately within a few centimeters.

LEO orbit accuracy
Since there is no official organization to release precise orbits for SAT-A/B/C/D/E, we first evaluate the orbit accuracy by comparing overlapping arcs.Assuming the current time is day0, the first orbit arc is set from 00:00 to 24:00 on day0, and the second orbit arc is shifted forward by 12 h, i.e., from 12:00 on day0 to 12:00 on day1, resulting in a 12-h overlapping orbits between the two arcs.To avoid the boundary effects, 1-h of orbit data is removed from each end of the 12-h overlapping arc, focusing on the 10-h overlapping orbit.Figure 8 shows the RMS values of overlapping arc differences in the R, T, N, and 3D direction.The horizontal axis represents the DOY in 2019, and the vertical axis represents the RMS of the overlapping orbit differences for each satellite.Table 4 shows the average RMS of the 5 LEO satellites.The accuracy of the overlapping arcs in the R direction is better than 1 cm and better than 2 cm in the 3D direction.After a thorough comparative analysis of overlapping arcs, the orbit accuracy of SAT-A/B/C/D/E can achieve an internal consistency accuracy of better than 2 cm. Figure 9 shows the orbit comparisons in the R, T, N and 3D direction between the reference orbits and the orbits derived from the integrated POD of 13 LEO satellites on January 20, 2019.To establish the reference orbit, we utilized a zero-difference orbit determination algorithm, introducing the fixed GPS satellite precise orbit and clock products provided by CODE to compute the orbits of SAT-A/B/C/D/E satellites.Following that, we chosen this orbit as the reference for accuracy assessment.The horizontal axis represents the abbreviations of the LEO satellites, and the vertical axis represents the RMS of the orbital difference for each satellite.
As shown in Fig. 9, when orbit determination is conducted for 13 LEO satellites on January 20, 2019, the R direction has the smallest RMS value, better than 2 cm, while the RMS in the T direction is maximal.The RMS in the 3D direction is better than 4 cm compared to the reference orbit.Figure 10 gives the detailed statistical results over 7 days.The horizontal axis represents the DOY from 20 to 26 in 2019, and the vertical axis represents the RMS values in the R, T, N, and 3D direction for each day.Table 5 shows the average RMS of 13 LEO satellites derived from the integrated POD.
From Fig. 10 and Table 5, the RMS values for a period of 7 days are better than 2 cm and 5 cm in the R and 3D direction.The orbit accuracy in the T direction is the worst, primarily due to the influence of atmospheric drag perturbation, which is challenging to be accurately modeled.LEO POD is usually done by fixing precise orbit and clock products of GNSS satellites generated with the measurements from global stations as known values.
The results indicate that with the integrated POD it is possible to generate precise LEO orbits with centimeterlevel in a regional network.Kang et al. (2020), van den Ijssel et al. (2015), Duan & Hugentobuler (2021) and Li et al., (2022b) reported that 1-3 cm orbit determination accuracy of LEO satellites is achievable by fixing the IGS GNSS satellite orbits generated with well-distributed global stations.Compared with the above results, the orbit accuracy of LEO satellites calculated by the integrated POD is slightly worse.The small number of regional receivers involved in the integrated orbit determination accounts for the larger orbit errors.Also, the refined processing strategy can help improve the LEO orbit determination, such as ambiguity resolution (Duan & Hugentobuler, 2021;Jiang et al., 2023;Zhang et al., 2021), which will be done soon.c) GPS satellite clock accuracy.
In the process of the integrated POD, the orbital parameters and clock offsets of the GPS satellites are estimated simultaneously.To evaluate the accuracy of clock offsets, 2-order difference method is introduced between  the results calculated by the integrated POD and COD final precise clock products, and the Standard Deviation (STD) of the 2-order differences is given as the final accuracy evaluation.The evaluation steps are expressed as: (1) Select one GPS satellite as reference, and the clock offsets of other satellites are compared with the ref-erence satellite to obtain 1-order differences at time t i .
(2) Calculate the difference of the 1-order difference series, so that the influence of the reference difference can be eliminated, and the 2-order differences of the clock offsets can be obtained, (4)  (3) Compute the STD value of the 2-order differences.
where the subscript pod and cod denote the inte- grated POD products and COD final products, respectively.The subscript i indicates the sampling epoch.and ∇ means to make 1-and 2-order dif- ference.The subscript n represents the total num- ber of epochs.∇ is the mean value of the 2-order differences. (7 Continuous tracking of GPS satellites is crucial for ensuring the accuracy of satellite clock offsets.In visibility analysis for a specific GPS satellite there are the epochs when the number of receivers is zero.In this case, the clock offsets cannot be solved.When the satellite is tracked again, the 2-order difference of the clock offsets at this epoch will experience a certain degree of jump.This is mainly because for the estimation of clock offsets a reference datum needs to be selected first.When data gaps occur in GPS tracking, it will lead to a switch in the selected reference, resulting in discontinuity in the estimated clock offsets.In this case, G10 is chosen as the reference datum for computing the 2-order differences in clock offsets.Figure 11 takes G01 as an example, the time series of 2-order differences between G01 and G10 is presented.
From Fig. 11, the 2-order difference series exhibit a certain degree of jump when the number of receivers is zero.This indicates a switch in the reference datum for clock offsets.Therefore, when evaluating the accuracy of clock offsets, it is necessary to perform segmented statistical analysis on the 2-order difference sequence.Figure 12 shows the RMS and STD values of the segmented 2-order differences between the clock offsets calculated by the integrated POD and the COD clock products.
Due to the limited observations of G04, G10 has been selected as the reference datum.Consequently, these two satellites are excluded from the statistical analysis.From Fig. 12, it is evident that among the remaining 30 satellites, although systematic biases exist between the clock offsets derived from the integrated POD and the COD final precise clock products, these biases can be entirely absorbed by the receiver clock offsets in the positioning model without affecting the positioning accuracy (Ye, 00:00 04:00 08:00 12:00 16:00 20:00 24:00 Time The relationship between 2-order differences and the number of receivers Fig. 12 The RMS and STD of 2-order clock differences of integrated POD compared with the COD final precise clock offsets 2002).Most of the STD values are better than 0.15 ns.Thus, with the integrated processing of the measurements from LEO satellites and regional ground stations, GPS satellite clock offsets with centimeter-level precision can be achieved.

PPP validation
After evaluating the accuracy of orbit and clock products, it is necessary to carry out positioning experiment using the products derived from the integrated POD.We selected the observation data of 18 IGS stations and conducted simulated kinematic PPP scenarios and discussed the impact of the discontinues in clock parameters estimation of GPS satellites.The PPP strategy is shown in Table 6.Although the stations are static, the station coordinates are treated as white noise and estimated epoch wise to assess the performance of the integrated POD products.The performance is assessed in terms of the convergence time and positioning accuracy.For comparison, simulated kinematic PPP validation is also carried out using the COD precise final products.The distribution of selected stations used for PPP validation is shown in the figure below (Fig. 13).

Convergence time
PPP is a convergence process that requires a long period of data accumulation to obtain an ideal positioning result.The convergence condition requires that the horizontal absolute positioning error remain below 10 cm and the vertical remain below 20 cm for at least 5 min in the corresponding direction.The average convergence time for the 18 IGS stations is given in Fig. 14.For the stations CHUR, CRO1, HOB2, and WHIT, the convergence time is significantly shorter when using POD products compared to using COD products.However, for most stations, the convergence time using POD products is longer than or equal to the convergence time using COD products.

Positioning accuracy
To verify the effects of the integrated POD products on kinematic PPP positioning accuracy, we compare the positioning results with the precise coordinates provided in the SINEX file released by IGS analysis center and give the differences in the north (N), east (E) and up (U) directions.The kinematic PPP with the COD final products is also calculated for comparison purpose.The station-specific positioning accuracy is given in Fig. 15.The average accuracy is summarized in Table 7.
In Fig. 15, the feasibility of the integrated products for PPP applications is verified and decimeter-level positioning accuracy can be obtained.From Table 7, positioning accuracy with the COD products are slightly better than those with the POD products but at the same level.It can be attributed to the larger orbit errors and clock offset errors of the integrated POD products compared to the COD products.Based on the visibility analysis, the GPS satellites are not always tracked by two or more receivers.There are the epochs when some GPS satellites are not tracked by any

Ionosphere delay
The first-order ionosphere delay is eliminated by the dual-frequency combinations

Receiver clock
One receiver clock offset is set and estimated as white noise.The initial value is given as 0 and prior covariance is given as 9*10 12 m 2

Receiver position Estimated as white noise
Phase ambiguities Estimated as float solutions and are not fixed as integers Estimation Forward Kalman filter receivers, resulting in unattainable clock offsets for the GPS satellites and the discontinues in clock parameters estimation.In the PPP processing with the integrated POD products, the number of available satellites is smaller than the number with COD products at a few epochs, the ambiguity parameters should be reset and estimated, resulting in larger positioning errors.With more LEO satellites participating in the integrated orbit determination, it is believed that the data gap of GPS satellite tracking will no longer exist.The estimated clock offsets will also become continuous, and the accuracy of both estimated orbit and clock will be improved.As a result, the performance of PPP validation will be more prominent.

Discussion and conclusions
This study initially validates the potential for generating high-precision GPS orbit and clock products by integrating the measurements from a regional network and LEO satellites.The orbit accuracy of GPS satellites in the R, T and N directions is 2.27 cm, 3.45 cm, and 3.08 cm, respectively, with the clock accuracy better than 0.15 ns.The LEO orbit accuracy is better than 2 cm in the R direction, and the position errors are mostly within 4 cm.Then based on the integrated POD products, kinematic PPP validation was carried out in terms of convergence time and positioning accuracy.Additionally, COD final products were utilized for comparative analysis.The results indicate that the PPP performance using POD products is generally inferior to that using COD products for most stations.This discrepancy is primarily attributed to larger orbit and clock errors in POD products compared to COD products.Moreover, the periods of missing observations can lead to the discontinuities in clock estimation.While it is evident that high precision GNSS orbits and clocks can be obtained by integrating the observations from regional stations and LEO satellites, it is important to note that the integrated solutions exhibit slightly lower accuracy in GNSS orbit, LEO orbit, and PPP performance compared to the solutions generated with the measurements of a well-distributed global network.Future efforts will focus on enhancing performance with ambiguity resolution and increased involvement of LEO satellites.

Fig. 2
Fig. 2 Number of GPS satellites tracked by 13 LEO satellites

Fig. 8
Fig. 8 RMS values in the R, T, N, and 3D direction of overlapping orbits derived from the integrated POD during DOY 20-26, 2019

Fig. 9 Fig. 10
Fig. 9 Orbit comparisons of 13 LEO satellites between the integrated POD products and COD final products on January 20, 2019 Fig. 11The relationship between 2-order differences and the number of receivers

Fig. 13 Fig. 14
Fig. 13 Distribution of 18 IGS stations used for PPP validation

Fig. 15
Fig. 15 Positioning accuracy of kinematic PPP based on COD final products and integrated POD products the ground station and LEO satellite.ρ s G,IF and ρ s L,IF represent the satellite-to-receiver distances for ground stations and LEO satellites relative to GNSS satellites.δt G and δt L denote the receiver clock offsets, respectively, whereas δt s is the satellite clock offsets.T s G represents the tropospheric delay experienced by the ground station.

Table 1
Detailed information on the 13 selected LEO satellites Fig. 1 Distribution of 6 IGS ground stations

Table 2
The models adopted in the integrated POD

Table 3
The parameter estimation strategy in the integrated POD

Table 4
Average RMS of overlapping orbits for SAT-A/B/C/D/E derived from the integrated POD during DOY20-26, 2019 (unit:  cm)

Table 5
Average RMS of 13 LEO satellites derived from the integrated POD during DOY20-26, 2019 (unit: cm)

Table 6
PPP strategy adopted in this research

Table 7
RMS values of kinematic PPP solutions based onCOD  and POD products (unit: cm)