GSFC DORIS contribution to ITRF2008

The NASA GSFC DORIS analysis center has provided weekly DORIS solutions from November 1992 to January 2009 (839 SINEX files) of station positions and Earth Orientation Parameters for inclusion in the DORIS contribution to ITRF2008. The NASA GSFC GEODYN orbit determination software was used to process the orbits and produce the normal equations. The weekly SINEX gscwd10 submissions included DORIS data from Envisat, TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5. The orbits were mostly seven days in length (except for weeks with data gaps or maneuvers). The processing used the GRACE-derived EIGEN-GL04S1 gravity model, updated modeling for time-variable gravity, the GOT4.7 ocean tide model and tuned satellite-specific macromodels for SPOT-2, SPOT-3, SPOT-4, SPOT-5 and TOPEX/Poseidon. The University College London (UCL) radiation pressure model for Envisat improves nonconservative force modeling for this satellite, reducing the median residual empirical daily along-track accelerations from 3.75 × 10⁻⁹ m/s² with the a priori macromodel to 0.99 × 10⁻⁹ m/s² with the UCL model. For the SPOT and Envisat DORIS satellite orbits from 2003 to 2008, we obtain average RMS overlaps of 0.8–0.9 cm in the radial direction, 2.1–3.4 cm cross-track, and 1.7–2.3 cm along-track. The RMS orbit differences between Envisat DORIS-only and SLR & DORIS orbits are 1.1 cm radially, 6.4 cm along-track and 3.7 cm cross-track and are characterized by systematic along-track mean offsets due to the Envisat DORIS system time bias of ±5–10 μs. We obtain a good agreement between the geometrically-determined geocenter parameters and geocenter parameters determined dynamically from analysis of the degree one terms of the geopotential. The intrinsic RMS weekly position repeatability with respect to the IDS-3 combination ranges from 2.5 to 3.0 cm in 1993–1994 to 1.5 cm in 2007–2008.     

Impact of GPS antenna phase center variations on precise orbits of the GOCE satellite

The first European Space Agency Earth explorer core mission GOCE (Gravity field and steady-state Ocean Circulation Explorer) has been launched on March 17, 2009. The 12-channel dual-frequency Global Positioning System receiver delivers 1 Hz data and provides the basis for precise orbit determination (POD) on the few cm-level for such a very low orbiting satellite (254.9 km). As a member of the European GOCE Gravity Consortium, which is responsible for the GOCE High-level Processing Facility (HPF), the Astronomical Institute of the University of Bern (AIUB) provides the Precise Science Orbit (PSO) product for the GOCE satellite. The mission requirement for 1-dimensional POD accuracy is 2 cm. The use of in-flight determined antenna phase center variations (PCVs) is necessary to meet this requirement. The PCVs are determined from 154 days of data and the magnitude is up to 3-4 cm. The impact of the PCVs on the orbit determination is significant. The cross-track direction benefits most of the PCVs. The improvement is clearly seen in the orbit overlap analysis and in the validation with independent Satellite Laser Ranging (SLR) measurements. It is the first time that SLR could validate the cross-track component of a LEO orbit.     

Generating precise and homogeneous orbits for Jason-1 and Jason-2

Driven by the GMES (Global Monitoring for Environment and Security) and GGOS (Global Geodetic Observing System) initiatives the user community has a strong demand for high-quality altimetry products. In order to derive such high-quality altimetry products, precise orbits for the altimetry satellites are a necessity. With the launch of the TOPEX/Poseidon mission in 1992 a still on-going time series of high-accuracy altimetry measurements of ocean topography started, continued by the altimetry missions Jason-1 in 2001 and Jason-2/OSTM in 2008. This paper contributes to the on-going orbit reprocessing carried out by several groups and presents the efforts of the Navigation Support Office at ESA/ESOC using its NAPEOS software for the generation of precise and homogeneous orbits referring to the same reference frame for the altimetry satellites Jason-1 and Jason-2. Data of all three tracking instruments on-board the satellites (beside the altimeter), i.e. GPS, DORIS, and SLR measurements, were used in a combined data analysis. About 7 years of Jason-1 data and more than 1 year of Jason-2 data were processed. Our processing strategy is close to the GDR-C standards. However, we estimated slightly different scaling factors for the solar radiation pressure model of 0.96 and 0.98 for Jason-1 and Jason-2, respectively. We used 30 s sampled GPS data and introduced 30 s satellite clocks stemming from ESOC’s reprocessing of the combined GPS/GLONASS IGS solution. We present the orbit determination results, focusing on the benefits of adding GPS data to the solution. The fully combined solution was found to give the best orbit results. We reach a post-fit RMS of the GPS phase observation residuals of 6 mm for Jason-1 and 7 mm for Jason-2. The DORIS post-fit residuals clearly benefit from using GPS data in addition, as the DORIS data editing improves. The DORIS observation RMS for the fully combined solution is with 3.5 mm and 3.4 mm, respectively, 0.3 mm better than for the DORIS-SLR solution. Our orbit solution agrees well with external solutions from other analysis centers, as CNES, LCA, and JPL. The orbit differences between our fully combined orbits and the CNES GDR-C orbits are of about 0.8 cm for Jason-1 and at 0.9 cm for Jason-2 in the radial direction. In the cross-track component we observe a clear improvement when adding GPS data to the POD process. The 3D-RMS of the orbit differences reveals a good orbit consistency at 2.7 cm and 2.9 cm for Jason-1 and Jason-2. Our resulting orbit series for both Jason satellites refer to the ITRF2005 reference frame and are provided in sp3 file format on our ftp server.     

Impact of DORIS ground antennas environment on their radio signal quality

Among the factors which may disrupt the DORIS measurements quality, the ground antennas environment is of high importance. For a set of 15 selected DORIS beacon, the differences between the effective and theoretical power received on-board the satellites (SPOT-5 and Envisat) have been analyzed in terms of spatial direction around the antenna. Such antenna maps have also been established regarding the Doppler residuals of the least-square precise orbit adjustment. Thanks to 360° views from the antennas and aerial views of the sites, the impact of the signal obstructions (trees, roofs, antennas …) on power attenuation and Doppler residuals is discussed. Depending on the nature of the obstructed object, the attenuation level can reach more than 5 dB, and the residual RMS of the orbit adjustment may be doubled from the nominal value, reaching 1 mm/s locally. The nature of the ground at the foot of the antennas has been correlated to DORIS signal quality at high elevation: reflections on flat surfaces (e.g. roofs) affect the signal more significantly than reflections on natural ground (e.g. soil). In particular, a modeling of the multipath phenomenon affecting Fairbanks site has been established and fits remarkably with the observations. Finally, an evaluation of the direct impact of obstructing objects on the orbit has also been performed. The example of a scaffolding at Kauai site displays a few millimeters error in the along-track position of the satellite.     

Hardware-in-the-loop simulations of GPS-based navigation and control for satellite formation flying

A relative navigation and formation control algorithm for satellite formation flying was developed, and a hardware-in-the-loop (HIL) simulation testbed was established and configured to evaluate this algorithm. The algorithm presented is a relative navigation estimation algorithm using double-difference carrier-phase and single-difference code measurements based on the extended Kalman filter (EKF). In addition, a state-dependent Riccati equation (SDRE) technique is utilized as a nonlinear controller for the formation control problem. The state-dependent coefficient (SDC) form is formulated to include nonlinearities in the relative dynamics. To evaluate the relative navigation and control algorithms developed, a closed-loop HIL testbed is configured. To demonstrate the performance of the testbed, a test formation flying scenario comprising formation acquisition and keeping in a low earth orbit (LEO) has been established. The relative navigation results from the closed-loop simulations show that a 3D RMS of 0.07 m can be achieved for position accuracy. The targeted leader–follower formation flying in the along-track separation of 100 m was maintained with a mean position error of approximately 0.2 m and a standard deviation of 0.9 m. The simulation results show that the HIL testbed is capable of successful demonstration of the GPS-based satellite autonomous formation flying mission.     

Precise orbit determination using the batch filter based on particle filtering with genetic resampling approach

In this study, genetic resampling (GRS) approach is utilized for precise orbit determination (POD) using the batch filter based on particle filtering (PF). Two genetic operations, which are arithmetic crossover and residual mutation, are used for GRS of the batch filter based on PF (PF batch filter). For POD, Laser-ranging Precise Orbit Determination System (LPODS) and satellite laser ranging (SLR) observations of the CHAMP satellite are used. Monte Carlo trials for POD are performed by one hundred times. The characteristics of the POD results by PF batch filter with GRS are compared with those of a PF batch filter with minimum residual resampling (MRRS). The post-fit residual, 3D error by external orbit comparison, and POD repeatability are analyzed for orbit quality assessments. The POD results are externally checked by NASA JPL’s orbits using totally different software, measurements, and techniques. For post-fit residuals and 3D errors, both MRRS and GRS give accurate estimation results whose mean root mean square (RMS) values are at a level of 5 cm and 10–13 cm, respectively. The mean radial orbit errors of both methods are at a level of 5 cm. For POD repeatability represented as the standard deviations of post-fit residuals and 3D errors by repetitive PODs, however, GRS yields 25% and 13% more robust estimation results than MRRS for post-fit residual and 3D error, respectively. This study shows that PF batch filter with GRS approach using genetic operations is superior to PF batch filter with MRRS in terms of robustness in POD with SLR observations.     

Orbit determination and prediction of GEO satellite of BeiDou during repositioning maneuver

In order to establish a continuous GEO satellite orbit during repositioning maneuvers, a suitable maneuver force model has been established associated with an optimal orbit determination method and strategy. A continuous increasing acceleration is established by constructing a constant force that is equivalent to the pulse force, with the mass of the satellite decreasing throughout maneuver. This acceleration can be added to other accelerations, such as solar radiation, to obtain the continuous acceleration of the satellite. The orbit determination method and strategy are illuminated, with subsequent assessment of the orbit being determined and predicted accordingly. The orbit of the GEO satellite during repositioning maneuver can be determined and predicted by using C-Band pseudo-range observations of the BeiDou GEO satellite with COSPAR ID 2010-001A in 2011 and 2012. The results indicate that observations before maneuver do affect orbit determination and prediction, and should therefore be selected appropriately. A more precise orbit and prediction can be obtained compared to common short arc methods when observations starting 1 day prior the maneuver and 2 h after the maneuver are adopted in POD (Precise Orbit Determination). The achieved URE (User Range Error) under non-consideration of satellite clock errors is better than 2 m within the first 2 h after maneuver, and less than 3 m for further 2 h of orbit prediction.     

CHAMP and GRACE accelerometer calibration by GPS-based orbit determination

Current and planned Earth observation missions are equipped with highly sensitive accelerometers. Before using the data, the instrument has to be calibrated by determining scale and bias parameters for each axis. Here, the accelerometer measurements are used in a GPS-based reduced-dynamic orbit determination approach, replacing the non-gravitational force models, and nominally daily calibration parameters are estimated. Additional empirical accelerations are estimated to account for deficiencies in the applied force models. This method is applied to 5 years of CHAMP and GRACE data, resulting in an orbit precision at the level of a few centimeters. In along-track direction the calibration parameters can be estimated freely, scale factors of 0.96 ± 0.014 and 0.95 ± 0.015 are obtained for GRACE A and B, and 0.85 ± 0.024 for CHAMP. A constant scale factor results in the smoothest bias series, with clear trends and occasional jumps. In radial and cross-track direction tight constraints to a priori biases have to be applied. Furthermore, the determined orbits are analyzed with respect to reference trajectories, and SLR, phase and KBR residuals are presented.     

圆轨道航天器八面体编队构型轨道动力学分析

提出一个全新的八面体航天器编队构型,该构型体现出当前编队飞行多种轨道构型的特征,同时也适于作为空间演示试验的编队飞行模式。八面体编队构型的设计思路是基于C-W方程,轨道平面内沿航向编队构型可利用轨道动力学自然保持,正上方或正下方编队构型则需要依靠平面内控制来实现,垂直轨道平面的编队构型需要施加法向控制来实现。对基于C-W方程的悬停动力学模型进行了精度分析,最后以低轨道航天器的八面体编队构型为例进行了数学仿真验证。     

BeiDou-3 orbit and clock quality of the IGS Multi-GNSS Pilot Project

Within the Multi-GNSS Pilot Project (MGEX) of the International GNSS Service (IGS), precise orbit and clock products for the BeiDou-3 global navigation satellite system (BDS-3) are routinely generated by a total of five analysis centers. The processing standards and specific properties of the individual products are reviewed and the BDS-3 orbit and clock product performance is assessed through direct inter-comparison, satellite laser ranging (SLR) residuals, clock stability analysis, and precise point positioning solutions. The orbit consistency evaluated by the signal-in-space range error is on the level of 4–8 cm for the medium Earth orbit satellites whereas SLR residuals have RMS values between 3 and 9 cm. The clock analysis reveals sytematic effects related to the elevation of the Sun above the orbital plane for all ACs pointing to deficiencies in solar radiation pressure modeling. Nevertheless, precise point positioning with the BDS-3 MGEX orbit and clock products results in 3D RMS values between 7 and 8 mm.     

Near real-time BDS GEO satellite orbit determination and maneuver analysis with reversed point positioning

The Geostationary Earth Orbit (GEO) satellite is a crucial part of the BeiDou Navigation Satellite System (BDS) constellation. However, due to various perturbation forces acting on the GEO satellite, it drifts gradually over time. Thus, frequent orbit maneuvers are required to maintain the satellite at its designed position. During the orbit maneuver and recovery periods, the orbit quality of the maneuvered satellite computed with broadcast navigation ephemeris will be significantly degraded. Furthermore, the conventional dynamic Precise Orbit Determination (POD) approach may not work well, because of a lack of publicly available satellite information for modeling the thrust forces. In this paper, a near real-time approach free of thrust forces modeling is proposed for BDS GEO satellite orbit determination and maneuver analysis based on the Reversed Point Positioning (RPP). First, the station coordinates and receiver clock offsets are estimated by GPS/BDS combined Single Point Positioning (SPP) with single-frequency phase-smoothed pseudorange observations. Then, with the fixed station coordinates and receiver clock offsets, the RPP method can be conducted to determine the GEO satellite orbits. When no orbit maneuvers occur, the proposed method can obtain orbit accuracies of 0.92, 2.74, and 8.30?m in the radial, along-track, and cross-track directions, respectively. The average orbit-only Signal-In-Space Range Error (SISRE) is 1.23?m, which is slightly poorer than that of the broadcast navigation ephemeris. Using four days of GEO maneuvered datasets, it is further demonstrated that the derived orbits can be employed to characterize the behaviors of GEO satellite maneuvers, such as the time span of the maneuver as well as the satellite thrusting accelerations. These results prove the efficiency of the proposed method for near real-time GEO satellite orbit determination during maneuvers.     

Precise orbit determination and baseline consistency assessment for Swarm constellation

Precise orbit determination (POD) and precise baseline determination (PBD) of Swarm satellites with 4 years of data are investigated. Ambiguity resolution (AR) plays a crucial role in achieving the best orbit accuracy. Swarm POD and PBD based on single difference (SD) AR and traditional double difference (DD) AR methods are explored separately. Swarm antenna phase center variation (PCV) corrections are developed to further improve the orbit determination accuracy. The code multipath of C1C, C1W and C2W observations is first evaluated and clear variations in code noise related to different receiver settings are observed. Carrier phase residuals of different time periods and different loop tracking settings of receiver are studied to explore the effect of ionospheric scintillation on POD. The reduction of residuals in the polar and geomagnetic equator regions confirms the positive impact of the updated carrier tracking loops (TLs) on POD performance. The SD AR orbits and orbits with float ambiguity (FA) are compared with the Swarm precise science orbits (PSOs). An average improvement of 27 %, 4 % and 16 % is gained in along-track, cross-track and radial directions by fixing the ambiguity to integer. For Swarm-A/B and Swarm-B/C formations, specific days are selected to perform the DD AR-based POD during which the average distance of the formation satellites is less than 5000 km. Satellite laser ranging (SLR) observations are employed to validate the performance of FA, SD AR and DD AR orbits. The consistency between the SD AR orbits and SLR data is at a level of 10 mm which shows an improvement of 25 % when comparing with the FA results. An SLR residuals reduction of 15 % is also achieved by the DD AR solution for the selected days. Precise relative navigation is also an essential aspect for spacecraft formation flying missions. The closure error method is proposed to evaluate the baseline precision in three dimensions. A baseline precision of 1–3 mm for Swarm-A/C formation and 3–5 mm for Swarm-A/B and Swarm-B/C satellite pairs is verified by both the consistency check and closure error method.     

近圆轨道卫星编队捕获技术研究

基于近圆参考轨道的假设,研究处于同一入轨点多颗卫星的编队捕获方法.首先由高斯型拉格朗日轨道摄动运动方程得到轨道坐标系中控制冲量与轨道根数偏差的关系,基于近圆轨道的条件简化并带入相对运动方程,得到控制冲量与相对运动的关系表达式;通过深入分析各个方向(径向、沿迹向与轨道面法向)的控制冲量对相对运动的影响,给出了分别用径向与轨道面法向控制冲量组合和沿迹向与轨道面法向控制冲量组合实现编队捕获的两种控制策略;最后给出了一个空间圆编队捕获实例,并从燃料消耗、施加冲量次数及捕获时间等角度对比研究了两种控制策略的特点.仿真结果表明,这两种控制策略简单、实用,能够较好地解决近圆轨道卫星编队的捕获问题.      

Measuring atmospheric density using GPS–LEO tracking data

We present a method to estimate the total neutral atmospheric density from precise orbit determination of Low Earth Orbit (LEO) satellites. We derive the total atmospheric density by determining the drag force acting on the LEOs through centimeter-level reduced-dynamic precise orbit determination (POD) using onboard Global Positioning System (GPS) tracking data. The precision of the estimated drag accelerations is assessed using various metrics, including differences between estimated along-track accelerations from consecutive 30-h POD solutions which overlap by 6 h, comparison of the resulting accelerations with accelerometer measurements, and comparison against an existing atmospheric density model, DTM-2000. We apply the method to GPS tracking data from CHAMP, GRACE, SAC-C, Jason-2, TerraSAR-X and COSMIC satellites, spanning 12 years (2001–2012) and covering orbital heights from 400 km to 1300 km. Errors in the estimates, including those introduced by deficiencies in other modeled forces (such as solar radiation pressure and Earth radiation pressure), are evaluated and the signal and noise levels for each satellite are analyzed. The estimated density data from CHAMP, GRACE, SAC-C and TerraSAR-X are identified as having high signal and low noise levels. These data all have high correlations with anominal atmospheric density model and show common features in relative residuals with respect to the nominal model in related parameter space. On the contrary, the estimated density data from COSMIC and Jason-2 show errors larger than the actual signal at corresponding altitudes thus having little practical value for this study. The results demonstrate that this method is applicable to data from a variety of missions and can provide useful total neutral density measurements for atmospheric study up to altitude as high as 715 km, with precision and resolution between those derived from traditional special orbital perturbation analysis and those obtained from onboard accelerometers.     

The impact of orbital errors on the estimation of satellite clock errors and PPP

Precise satellite orbit and clocks are essential for providing high accuracy real-time PPP (Precise Point Positioning) service. However, by treating the predicted orbits as fixed, the orbital errors may be partially assimilated by the estimated satellite clock and hence impact the positioning solutions. This paper presents the impact analysis of errors in radial and tangential orbital components on the estimation of satellite clocks and PPP through theoretical study and experimental evaluation. The relationship between the compensation of the orbital errors by the satellite clocks and the satellite-station geometry is discussed in details. Based on the satellite clocks estimated with regional station networks of different sizes (∼100, ∼300, ∼500 and ∼700 km in radius), results indicated that the orbital errors compensated by the satellite clock estimates reduce as the size of the network increases. An interesting regional PPP mode based on the broadcast ephemeris and the corresponding estimated satellite clocks is proposed and evaluated through the numerical study. The impact of orbital errors in the broadcast ephemeris has shown to be negligible for PPP users in a regional network of a radius of ∼300 km, with positioning RMS of about 1.4, 1.4 and 3.7 cm for east, north and up component in the post-mission kinematic mode, comparable with 1.3, 1.3 and 3.6 cm using the precise orbits and the corresponding estimated clocks. Compared with the DGPS and RTK positioning, only the estimated satellite clocks are needed to be disseminated to PPP users for this approach. It can significantly alleviate the communication burdens and therefore can be beneficial to the real time applications.     

SPOT-2 and TOPEX/Poseidon precise orbit determination from DORIS doppler tracking

The French earth observation satellite SPOT-2 has served as a testbed for precise orbit determination from DORIS doppler tracking in anticipation of the TOPEX/Poseidon mission. Using the most up-to-data gravity field model, JGM-2, a radial orbit accuracy of about 2–9 cm was achieved, with an rms of fit of the tracking data of about 0.64 mm/s. Furthermore, it was found that the coordinates of the ground stations can be determined with an accuracy of the order of 2–5 cm after removal of common rotations, and translations.

Using a slightly different model for atmospheric drag, but the same gravity model, precise orbits of TOPEX/Poseidon from DORIS tracking data were determined with a radial orbit accuracy of the order of 4–5 cm, which is far within the 13 cm mission requirement. This conclusion is based on the analysis of 1-day overlap of successive 11-day orbits, and the comparisons with orbits computed from satellite laser tracking (SLR) and from the combination of SLR and DORIS tracking. Results indicate a consistency between the different orbits of 1–4 cm, 4–20 cm, and 6–13 cm in the radial, cross-track, and along-track directions, respectively. The residual rms is about 4–5 cm for SLR data and 0.56 mm/s for DORIS tracking. These numbers are roughly twice as large as the system noise levels, reflecting the fact that there are still some modeling errors left.     

Assessing attitude error of FORMOSAT-3/COSMIC satellites and its impact on orbit determination

An attitude determination and control system (ADCS) is critical to satellite attitude maneuvers and to the coordinate transformation from the inertial frame to the spacecraft frame. This paper shows specific sensors in the ADCS of the satellite mission FORMOSAT-3/COSMIC (F3/C) and the impact of the ADCS quality on orbit accuracy. The selection of main POD antenna depends on the beta angles of the different F3/C satellites (for FM2 and FM4) during the inflight phase. In particular, under the eclipse, alternative attitude sensors are activated to replace the Sun sensors, and such a sensor change leads to anomalous GPS phase residuals and a degraded orbit accuracy. Since the nominal attitude serves as a reference for ADCS, the 3-dimensional attitude-induced errors in reduced dynamic orbits over selected days in 2010 show 9.35, 10.78, 4.97, 5.48, 7.18, and 6.89 cm for FM1–FM6. Besides, the 3-dimensional velocity errors induced by the attitude effect are 0.10, 0.10, 0.07, 0.08, 0.09, and 0.10 for FM1–FM6. We analyze the quality of the observed attitude transformation matrix of F3/C and its impact on kinematic orbit determination. With 249 days of GPS in 2008, the analysis leads to the following averaged 3-dimensional attitude-induced orbit errors: 2.72, 2.62, 2.37, 1.90, 1.70, and 1.99 cm for satellites FM1–FM6. Critical suggestions of geodetic payloads for the follow-on mission of F3/C are presented based on the current result.     

基于双频GPS观测的简化动力学最小二乘批处理精密定轨

星载双频GPS载波相位和伪距观测量已成为低轨卫星获取精确三维位置和速度信息的主要方式. 本文以非差消电离载波相位和伪距组合作为观测量,应用简化动力学最小二乘批处理方法进行地球低轨卫星的精密定轨,并给出完整定轨流程. 采用逐段常量的经验加速度对动力学模型误差进行补偿,描述了经验加速度敏感矩阵及稀疏带状矩阵求逆的有效计算方法. 利用GRACE-A卫星GPS观测数据对定轨位置精度进行分析,结果显示,三维位置定轨精度优于5cm,经验加速度在径向、切向和法向上的补偿水平不超过40nm·s-2,大气阻力系数和辐射光压系数的估计值符合物理实际,星载接收机钟差大致呈线性并具有短周期小波动.      

Time synchronization of new-generation BDS satellites using inter-satellite link measurements

Autonomous satellite navigation is based on the ability of a Global Navigation Satellite System (GNSS), such as Beidou, to estimate orbits and clock parameters onboard satellites using Inter-Satellite Link (ISL) measurements instead of tracking data from a ground monitoring network. This paper focuses on the time synchronization of new-generation Beidou Navigation Satellite System (BDS) satellites equipped with an ISL payload. Two modes of Ka-band ISL measurements, Time Division Multiple Access (TDMA) mode and the continuous link mode, were used onboard these BDS satellites. Using a mathematical formulation for each measurement mode along with a derivation of the satellite clock offsets, geometric ranges from the dual one-way measurements were introduced. Then, pseudoranges and clock offsets were evaluated for the new-generation BDS satellites. The evaluation shows that the ranging accuracies of TDMA ISL and the continuous link are approximately 4?cm and 1?cm (root mean square, RMS), respectively. Both lead to ISL clock offset residuals of less than 0.3?ns (RMS). For further validation, time synchronization between these satellites to a ground control station keeping the systematic time in BDT was conducted using L-band Two-way Satellite Time Frequency Transfer (TWSTFT). System errors in the ISL measurements were calibrated by comparing the derived clock offsets with the TWSTFT. The standard deviations of the estimated ISL system errors are less than 0.3?ns, and the calibrated ISL clock parameters are consistent with that of the L-band TWSTFT. For the regional BDS network, the addition of ISL measurements for medium orbit (MEO) BDS satellites increased the clock tracking coverage by more than 40% for each orbital revolution. As a result, the clock predicting error for the satellite M1S was improved from 3.59 to 0.86?ns (RMS), and the predicting error of the satellite M2S was improved from 1.94 to 0.57?ns (RMS), which is a significant improvement by a factor of 3–4.