Orbit determination and time synchronization for BDS-3 satellites with raw inter-satellite link ranging observations

To provide competitive global positioning and timing services under the condition that monitoring stations are confined to Chinese territory, inter-satellite link (ISL) technology is used by the third-generation BeiDou Navigation Satellite System (BDS-3). The ISL, together with the dual one-way links between satellites and anchor stations, may enable autonomous navigation for BDS-3. In this paper, we propose a general observation model for orbit determination (OD) and time synchronization (TS) directly using non-simultaneous observations, such as raw ISL pseudoranges. With the proposed model, satellite orbits, clocks, and hardware delay biases of ISL equipment can be determined simultaneously by jointly processing inter-satellite one-way pseudorange data and observation data from ground monitoring stations. Moreover, autonomous OD and TS are also achievable with one-way pseudorange data from anchor stations and satellites. Data from eight BDS-3 satellites, two anchor stations, and seven monitoring stations located in China were collected to validate the proposed method. It is shown that by jointly processing data from the ISL and seven monitoring stations, the RMS of overlap orbit differences in radial direction is 0.019 m, the overlap clock difference (95%) is 0.185 ns, and the stability of the estimated hardware delay biases for each satellite is greater than 0.5 ns. Compared with the results obtained with the seven stations, the improvements of orbits in radial direction and clocks are 95.7% and 90.5%, respectively. When the hardware delay biases are fixed to predetermined values, the accuracies of orbits and clocks are further improved. By jointly processing pseudoranges from the satellites and the two anchor stations, the RMS of overlap orbit differences is 0.017 m in the radial direction, and the overlap clock difference (95%) is 0.037 ns. It has also been demonstrated that under the condition of one-way ranging links, the accuracies of orbits and clocks obtained by the above two modes are still significantly better than those obtained by using the data from the monitoring stations alone.


Introduction
To enhance orbit determination (OD) and time synchronization (TS) under the condition of regional monitoring stations, the third-generation BeiDou Navigation Satellite System (BDS-3) employs inter-satellite link (ISL) technology which also enable autonomous navigation. The feasibility and performance of the ISL in BDS-3 have been verified with the in-orbit validation satellites [1][2][3][4][5][6][7]. Recent findings show that with the ISL ranging data, the radial orbit accuracy of 10 cm is achievable [2].
Phased array antennas on BDS-3 satellites receive and transmit ISL signals in Ka-band and control the beam pointing of the signal. The Ka-band signal with frequencies more than 10 times the L-band signal would provide greater communication bandwidth and make the ISL antennas smaller in size, which make it preferred for ISL. The ISL of BDS-3 satellites operates in a time-division multiple access (TDMA) mode, i.e., at any time (slot),

Open Access
Satellite Navigation https://satellite-navigation.springeropen.com/ *Correspondence: rrg2002me@163.com 1 Xi'an Research Institute of Surveying and Mapping, Xi'an 710054, China Full list of author information is available at the end of the article one satellite can only communicate and perform ranging measurement with another satellite according to the preset link planning. Pseudorange measurements observed at one satellite to others are generated in different time slots.
In existing methods of OD or TS for Global Navigation Satellite System (GNSS) satellites with ISL ranging data, the raw non-simultaneously observed pseudoranges are not directly used. Instead, in a typical approach originally proposed for autonomous OD and TS [8], observations are derived at specified times, and then the derived clock measurement (DCM) and the derived ephemeris measurement (DEM) are formed to separate geometry information and clock information [8][9][10][11]. The DEMs are used for OD alone [2,[5][6][7][8][9][12][13][14] or together with other data such as the monitoring station data [3] and the inter-satellite orientation observations [15], while the DCMs are used for TS [4]. Therefore, OD and TS are performed independently. In another approach, OD and TS are performed simultaneously by directly processing the one-way derived ISL observations together with the monitoring station data [16][17][18]. An advantage of the latter method is that the hardware delay biases of ISL equipment of each satellite can be estimated together with satellite orbits and clocks. However, the latter has only been used in simulation studies at present, and there have been no reports on the application for BDS-3.
A method to directly process raw one-way ISL pseudoranges for OD and TS of BDS-3 satellites was proposed by Ruan [19,20]. In this paper, the principal of this method will be presented in detail, and its feasibility will be demonstrated with two typical modes for OD and TS: (1) jointly processing the ISL pseudorange data and the L-band pseudorange and carrier phase data from monitoring stations, and (2) jointly processing the Ka-band ISL pseudorange data between satellites and the Ka-band pseudorange data between satellites and anchor stations. The anchor station refers to a station equipped with ISL equipment and can perform communication and ranging with BDS-3 satellites, just as the satellites do. Besides, it will be demonstrated that the proposed OD and TS method is achievable even under the condition of oneway ranging links between satellites (and between satellites and anchor stations).

Observation equation for ISL pseudoranges
The ISL pseudorange observations of satellite j observed by satellite i at time t i,j are expressed as (1) where τ i,j is the propagation time of ISL signal from satellite j to satellite i ; R i and R j are the instantaneous positions of satellite i and j at the receiving and transmitting times, respectively; δ i and δ j are the clock offsets of satellite i and j , respectively; j and β i are the transmitting and receiving hardware delays of satellite j and i , respectively, considered to be stable in time, and ω is the measurement noise [2,17]. Equation (1) omits terms such as relativistic delay of signal propagation and satellite antenna phase center offsets.
Equation (1) can also be used for the pseudorange observations between BDS-3 satellites and anchor stations, when the troposphere delay and ionosphere delay are taken into account. Since the BDS-3 ISL runs in Ka-band, pseudorange observations between satellites or between satellites and anchor stations are uniformly referred to as Ka-pseudoranges.

Observation equation for pseudoranges and phases at monitoring stations
Without losing generality, ionosphere-free carrier phase L i r and pseudorange P i r of satellite i observed by station (receiver) r at the time t k can be expressed as where k is the observation epoch; τ i r is the propagation time of signal from satellite i to station r ; R r is the instantaneous position of the station at the receiving time; θ r,k = δ r,k + b r and θ i k = δ i k − b i are the estimable clock parameters for the station and the satellite, i.e., superposition of the clock offset ( δ r,k and δ i k ) and the ionosphere free combination of hardware group delay ( b r and b i ); a i r is the ionosphere-free phase ambiguity parameter, including the hardware phase and group delays; T i r is the troposphere delay; ε and ξ are measuring errors of carrier phase and pseudorange measurements, respectively [21]. Equation (2) omits the correction terms of relativistic delay, phase wind-up effect and antenna phase centre offsets of the satellite and the station.

General observation equation for non-simultaneous observations
Pseudoranges and carrier phases [Eq. (2)] and Ka-pseudoranges [Eq. (1)] are radio measurements that contains clock offsets of the signal transmitting and receiving devices. In order to directly process the non-simultaneous observation data, such as the raw Ka-pseudorange observations of ISL, a general observation model is proposed in this section.
Continuous time is divided into discrete non-overlapping time slots with length of w . In any time slot k, the behavior of clocks on transmitting or receiving devices can be modeled with a polynomial, thus, clock offsets at any time t within time slot k are, according to Ruan [19], where t w,k is the selected reference time for the time slot k and is preferably at the middle point of the time slot. δ i k is the zero-order coefficient of the polynomial, i.e., clock bias at t w,k ; δ i k and δ i k are first and second order coefficients, i.e., drift and acceleration of the clock, respectively. Terms of higher order are omitted in Eq. (3) (usually terms higher than the second order can be safely ignored). In practice, accurate modeling of long-term changes of clocks using polynomials is challenging. In order to control the error of the clock polynomial models, the length of time slots should not be too long, depending on the frequency instability of an individual clock. For example, in a short time period, such as a few minutes, the clock of a GNSS satellite can usually be modeled using a linear polynomial, with adequate precision.
In this way, a radio measurement O r,e for a transmitting device e observed at t r,e by the receiving device r can be expressed as, according to Ruan [19], where t r,e is within the time slot k, t r,e k = t r,e − t w,k is the difference between observing time and the reference time t w,k , ρ r,e is the distance from a transmitting antenna to a receiving antenna, and ω 0 is the measuring error. The above formula does not include other factors that may affect the observation.
Specifically, if clock offsets are modeled with linear polynomials, the observation equation, Eq. (1), for Ka-pseudoranges observed at time t i,j , which is within the time slot k, can be rewritten as Thus, for any satellite-pair, all Ka-pseudoranges within the same time slot share the same clock parameters, although they are observed at different times. For example, a satellite's clock parameters are first and second order coefficients in Eq. (5). The non-simultaneously observed Ka-pseudoranges of satellites (and anchor stations) can be directly processed separately or jointly with other tracking data, estimation of orbits, clocks, and even hardware delay biases can be achieved simultaneously.
O r,e (t r,e ) = ρ r,e + δ r k +δ r k · (�t r,e )+ · · · − δ e k +δ e k · (�t r,e − τ r,e )+ · · · + ω o , In practice, the drift of satellite clocks can usually be predicted with a relatively high accuracy; therefore, the first order coefficients of clock models in Eq. (5) can be replaced by predicted values. For example, the uncertainty of predicted clock drift of a BDS satellite is less than 1 × 10 −13 s/s [2], and if the linear variation of the clock is corrected with the predicted clock drift, the resulting modeling error is less than 2 mm for time slots with length less than 1 min. In extreme cases without predicted clock drift, clock variation during short durations can even be directly ignored. Therefore, with predicted clock drift or by ignoring the clock variation, Eq. (5) is simplified to where ω is the sum of measuring error and additional errors caused by predicted clock drift or ignoring the clock drift.
The conventional observation model for pseudorange and phase observations observed at monitoring stations, such as in Eq. (2), can be regarded as a special case of Eq. (4) because each time slot contains only one epoch; thus, the model for the receiver clock is simplified with only the zero-order term. For any satellite, the time difference of signal transmitting corresponding to observations simultaneously observed by different stations at the same epoch is very small (i.e., much less than 1 s), so the variation of satellite clock can also be completely ignored.

Orbit determination and time synchronization of BDS-3 satellites
For OD and TS of the BDS-3 satellites, there are several types of tracking data that can be used, including pseudoranges and carrier phases observed at monitoring stations, mainly located within Chinese territory, Kapseudoranges observed at satellites and anchor stations, and two-way satellite time frequency transferring data between stations and satellites. For the proposed general observation model, the orbits and clocks of BDS-3 satellites can theoretically be determined by independently processing any type of tracking data or jointly processing several data types. This paper focuses on two modes for OD and TS: • Mode-1 Jointly processing ISL Ka-pseudoranges observed from satellites and observation data from monitoring stations; • Mode-2 Processing Ka-pseudoranges observed from both satellites and anchor stations.
In Mode-1, satellite orbits, clocks, and the hardware delay biases of each satellite are determined simultaneously [17]. This is a significant mode for BDS-3 and may be the best OD and TS mode, as it is an important and self-consistent method, which can resolve the inconsistency of hardware delays of the ISL equipment and the navigation signal generation equipment. In order to construct the joint observation model with both ISL data and the data from monitoring stations, the simplest method is to adopt the observation epoch time t k of the monitoring station as the reference time of specific time slots; i.e., t w,k = t k . Thus, the zero-order coefficient of the clock polynomial is exactly the clock offset θ i k at the time t k . Omitting detailed derivation, joint observation equation with ISL Ka-pseudoranges and ionosphere-free pseudoranges, and phases at time slot k centered on its reference time t k is expressed as follows: are the transmitting and receiving hardware delay biases for ISL equipment on satellite * ( = i, j ) [17], referred to as transmitting and receiving Ka-biases, respectively, in this context. The definition of Ka-biases and the estimation method are presented in Ruan et al. [16,17]. In Eq. (7), the observation time of the Ka-pseudoranges is usually different from the sampling time of the L-band pseudorange and phase data, but they share the same clock parameters. After modeling and linearization, unknown parameters such as satellite orbit, clock parameters, and Ka-biases are solved simultaneously, e.g., with the least square method.
Mode-2 uses Ka-pseudoranges observed from both satellites and anchor stations, and is usually referred to as (centralized) autonomous OD and TS [2,10]. The observation equations are where a represents anchor stations and T i a is the troposphere delay in Ka-pseudoranges between satellites and anchor stations. Ka-biases are not present in Eq. (8) because they are linearly correlated with the clock parameters when Ka-pseudoranges are used alone and must be calibrated with known values, e.g., estimates obtained by Mode-1, to achieve an unbiased clock solution for positioning and timing applications.

Data set and experimental scheme
To validate the proposed approach, the Ka-pseudorange data from eight BDS-3 satellites (B19-B22, B27-B30), two anchor stations (A01 and A02), and the pseudorange and carrier phase data from B1I and B3I frequency from seven monitoring stations of the International GNSS Monitoring and Assessment System (iGMAS) during days 152-163 of 2018 were collected. The distribution of the iGMAS stations used in the experiment is shown in Fig. 1. Table 1 shows the signal transmission and reception between satellites or between satellites and anchor stations during the data arc. It can be seen that each satellite has a bidirectional inter-satellite measurement link with 5-6 satellites, and each satellite has a bidirectional measurement link with the two anchor stations, except satellite B30.
With the above data, validation of OD and TS was performed by the following experiments: • EXP1: Only data from the iGMAS stations were used.
• EXP2: Data from the iGMAS stations and Ka-pseudoranges between satellites were jointly processed with Ka-biases of each satellite estimated simultaneously. • EXP3: Data from the iGMAS stations and Ka-pseudoranges between satellites were jointly processed with Ka-biases of each satellite calibrated with the mean estimated values obtained in EXP2. • EXP4: Ka-pseudoranges from satellites and anchor stations were processed with Ka-biases of each satel-lite and anchor station calibrated with the mean estimated values in EXP2 and the estimated values from literature [19], respectively.
Data were processed with the SPODS software [22,23] on a daily basis with 72-h data started on each day. Time slots of 1 min were assigned every 5 min and centered on the sampling epoch, i.e., for the data from monitoring stations, the sampling interval was 5 min. The coordinates of the iGMAS stations, anchor stations, and the earth rotation parameters were fixed to known values. The following orbital dynamics were carefully modeled: gravity of the earth computed with the EGM2008 up to degree and order 12, with linear variation of low-degree coefficients and periodical variation caused by the solid earth tides and the pole tide [24], gravity of the thirdbodies (the moon, sun, and planets) computed with JPL ephemeris DE405, and the general relativistic effects. For each satellite, without any a priori model, the ECOM model [25] was used to model the force acting on the satellites from solar radiation pressure with five parameters ( D 0 , Y 0 , B 0 , B c , and B s ), estimated together with the initial state vector. The Saastamoinen model [24] was used to calculated the hydrostatic zenith troposphere delay, with wet zenith delay parameter estimated every 2 h for the monitoring stations and every 24 h for anchor stations. The Global Mapping Function [26] was adopted to map the troposphere delay from the zenith to the line of sight. The satellite clock drift was replaced with the first order

Table 1 Measurement relationships between satellites or between satellites and anchor stations
Signal transmission B19 B20 B21 B22 B27 B28 B29 B30 A01 A02 coefficient of the broadcast clock correction model. In order to remove the rank defects caused by clock parameters, one clock parameter was chosed as the reference and constrained to 0 in each time slot. Since none of the monitoring stations, anchor stations, and satellites have observation data covering the whole arc, reference clock often changes in different time slots and can be a monitoring station, satellite, or an anchor station, at a certain time slot, depending on their observation number. In EXP2 and EXP5 (see "Orbit determination and time synchronization with single one-way ISL" section), one Ka-bias is chosen as the reference and constrained to 0 because not all Ka-biases are estimable, e.g., the receiving Ka-bias of B21. The pseudoranges and phases of the monitoring stations are weighted according to a prior precisions of 2 m and 0.02 m, respectively, and the Kapseudorange is weighted according to a prior precision of 0.1 m.
In order to evaluate the accuracy of the obtained orbits and clocks, differences of values in a 48-h overlap arc of adjacent solutions are calculated. The RMS differences for orbits are computed as an indicator of orbit accuracy. Mean value of differences over all satellites are firstly computed and deducted from the original overlapping differences at each epoch; then, the 95th percentile of the resultant differences is an indicator of clock accuracy. Figure 2 shows the RMSs in R, T, and N directions for each satellite obtained by the four experiments (note that the RMSs of EXP2, EXP3 and EXP4 have a tenfold magnification). It can be seen that in EXP1, the RMSs are larger than 0.3 m for R, larger than 1 m for T, and range from 1.2 to 2.3 m for N. In EXP2, the RMSs are smaller than 0.025 m in R, are not larger than 0.1 m for T, and are not larger than 0.15 m for N. With the Ka-biases of each satellite fixed to the multi-day average estimation in EXP2, the orbit accuracy in EXP3 is further improved. In EXP4, the orbit accuracy in R is very close to those in EXP2 and EXP3 for each individual satellite. However, the orbit accuracies in T and N are worse than those of EXP2 and EXP3, as there are only two anchor stations leading to looser constraints on the constellation orientation. Figure 3 shows the clock accuracy measured with overlap differences (95%) for each satellite in the four experiments. It can be seen that in EXP1 the clock accuracy for each satellite is between 1.5 and 2.5 ns. The accuracies for each satellite are better than 0.25 ns and 0.07 ns for EXP2 and EXP3, respectively, and the accuracy for each satellite in EXP4 is better than that in EXP3. Table 2 shows the accuracies of orbits and clocks obtained by different schemes. Compared with EXP1, improvements of orbit accuracies in EXP2 are 95.7%, 94.5%, 92.6%, and 93.4% in R, T, N, and 3D position, respectively. While, in EXP3, improvements of 95.7%, 94.5%, 92.6%, and 93.4% are achieved for R, T, N, and 3D position, respectively. In EXP4, the orbit accuracy in the R direction for each satellite is close to that in EXP3 but less than those of EXP2 and EXP3 in T and N directions. For satellite clocks, compared with EXP1, the accuracies are improved by 90.5% and 97.2% in EXP2 and EXP3, respectively. Since the overlap clock differences in EXP2 include the difference of estimated Ka-biases, the nominal accuracy of clocks is less than those in EXP3 and EXP4. The orbit accuracy in EXP4 is less than those in EXP3 and EXP2 but the clock accuracy is greater, implying model inaccuracies. For example, there may be biases of the antenna phase center for the Ka-band and the L-band, or there may be non-ignorable variation of the hardware delay for the Ka-band and the L-band signal. Figure 4 shows the variations of estimated receiving and transmitting Ka-biases of each satellite obtained in EXP2 with the bias removed for each individual satellite, so all the values are within the range of the y-axis. The Ka-biases of each satellite are stable in time. In Table 3, the standard deviation (SD) of the receiving and transmitting hardware biases for each satellite is smaller than 0.15 m (or equally 0.5 ns). Figure 5 shows the RMS of residuals for one-way Kapseudoranges observed at each satellite or anchor station in EXP2, EXP3, and EXP4. The RMSs of all satellites are smaller than 0.08 m; the RMSs for satellite B27-B30 are larger than those of other satellites, and satellite B28 has the largest RMS. The RMSs of individual satellites from different experiments are very close to each other. In EXP2, the average RMS is 0.054 m for all satellites, while it is 0.040 m for B19-B22 and 0.067 m for B27-B28. It can be seen that the RMSs of the two anchor stations are larger than those of the satellites. It is noticed that the RMSs for satellite B27-B30 in EXP2 and EXP3 are larger than those in EXP4, which may also be due to incorrect assumptions for the stability of hardware delays. Figure 6 shows the percentage of valid one-way Kapseudorange data in EXP3 and EXP4 for each satellite or anchor station. The percentage of valid data received or transmitted at satellites and anchor stations is more than 97% and the percentage of valid data received at B28 is the lowest. In EXP4, the percentage for each satellite is larger than that in EXP3. The average percentage for all satellites is 98.82% in EXP3 and 99.55% in EXP4. This also confirms that there are possible modeling errors, when the Ka-band data and the L-band monitoring data are jointly processed.

Orbit determination and time synchronization with single one-way ISL
Compared to the existing methods, only single one-way ranging data is necessary for OD and TS in the proposed method. In order to prove this point, EXP5, EXP6, and EXP7 were performed with the same data processing strategies as used in EXP2, EXP3, and EXP4, but only single one-way Ka-pseudorange data from measurement links were used (denoted by the red "√" in Table 1).
The accuracies of orbits for each satellite in EXP5, EXP6, and EXP7 are shown in Fig. 7. We can see that, the accuracies of orbits in each component in all experiments are all greater than 0.05 m. As expected, the accuracies in EXP6 are slightly higher than those in EXP5, and the accuracies are greater than 0.04 m, 0.1 m, and 0.2 m in R, T, and N, respectively, for both experiments. In EXP7, the orbit accuracy in R is similar to those in EXP5 and EXP6, but the accuracies in T and N are about 0.3 m and much worse than those in EXP5 and EXP6, for the same reason mentioned above. Figure 8 compares the clock accuracies obtained by different processing modes under the condition of single or double one-way Ka-pseudoranges. It can be seen that in EXP5, the clock accuracies of different satellites are between 0.2 and 0.35 ns, and the clock accuracies obtained by EXP6 and EXP7 are similar, and all of them are no more than 0.2 ns. Table 4 shows the accuracies for orbits and clocks obtained by EXP5, EXP6, and EXP7. It can be seen that the orbit accuracy in R is greater than 0.027 m in the three experiments, and the accuracies are decreased by up to 60% compared with those obtained with double one-way data. In other words, the radial orbital accuracies are improved by up to 40% using double one-way data. The clock accuracies obtained by the three experiments are   Table 3, the clock accuracies decreased by 42%, 183%, and 297% in EXP5, EXP6, and EXP7, respectively. This implies that the clocks benefit more from double one-way data than the orbits do. Table 5 provides the repeatability of the estimated Kabiases for each satellite obtained by EXP5 and the differences between the mean estimated Ka-biases obtained by EXP5 and those obtained by EXP2. The SDs of the obtained receiving biases in EXP5 do no exceed 0.18 m and the differences compared with those in EXP2 are less than 0.082 m. For the transmitting biases, the SDs are smaller than 0.14 m while the differences do not exceed 0.04 m.

Discussion
Since the ISL of BDS-3 runs in a TDMA mode, the Kapseudorange observations at a certain satellite are generated at different times, and each observation contains two unique clock offset parameters that remove redundant observations. Existing methods for OD and TS with ISL data firstly derive the raw Ka-pseudorange observations from different times to some specified time epoch, and then the DCMs and the DEMs are formed to separate the orbit and clock information so that OD and TS are separated [17]. This idea has been first proposed by Ananda et al. [8] and applied to OD and TS for BDS-3 Fig. 6 The percentage of valid one-way Ka-pseudorange data in EXP3 and EXP4 for each satellite or anchor station. R and E in the legend represent reception and transmission, respectively satellites in the operational control center [2][3][4]6]. It may have advantages in stability and computational efficiency for on-board distributed autonomous OD and TS, but for routine OD and TS at the operational control center or for the centralized autonomous OD and TS, its shortcomings are obvious. For example, data processing with this method is cumbersome, and it is difficult to calibrate the hardware delay of ISL equipment. To compute DEMs and DCMs, if either of the required double oneway pseudorange observations is not present (e.g., satellite observations do not arrive in time) or lacks quality, then the other will be discarded. If the method of polynomial interpolation or fitting is employed to compute derived pseudorange observations for a certain satellitepair, a sequence of one-way pseudorange observations with a certain length is needed for each derived observation, and the interval of the sequence must contain the "observing" time of the required derived observation. The method proposed in this paper directly uses the raw ISL pseudoranges for OD and TS simultaneously, and does not require either to generate derived observation or to match one-way pseudoranges observed at satellite-pairs. The procedure of data processing is very simple, almost the same as that for processing observation monitoring station data. Compared with the existing methods, the greatest advantage of this method is that it can simultaneously determine satellite orbits and clocks and the hardware delay biases of ISL equipment on satellites, if observations from monitoring stations are jointly processed. The obtained orbits and clocks are more selfconsistent, and the estimated ISL hardware biases can be used to correct the ISL observation data for autonomous OD and TS. Another advantage of directly processing raw one-way pseudorange data is that data quality or performance of each ISL transmitting or receiving device at each satellite can be analyzed independently.
The proposed method for OD and TS requires only single one-way pseudorange data, so there is no need to wait for the return of the data observed at those satellites that are beyond tracking by the ground stations, increasing the data availability. Only 20% of all the ISL data were used for validation (the length of time slots was 1 min and the sampling interval was 5 min). This implies that the frequency for ranging or transferring ISL data can be reduced to allow more time for inter-satellite communication or more channel resources for transferring other information.

Conclusions
For OD and TS of BDS-3 satellites with the ISL operated in a TDMA mode, a general model using simultaneous or non-simultaneous observation data is proposed. The proposed approach was validated with experiments for two representative and important modes of OD and TS:   (1) jointly processing ISL Ka-pseudoranges between satellites and ionosphere-free pseudorange and phase data from monitoring stations, and (2) jointly processing Kapseudoranges between satellites and between satellites and anchor stations. The results show that by jointly processing ISL data and the data of seven iGMAS stations located in China, the average RMS overlap orbit difference of the eight BDS-3 satellites is 0.019 m for radial direction and 0.156 m for 3D position; the overlap clock differences (95%) are 0.185 ns, and the standard deviation of the estimated Ka-biases is greater than 0.5 ns. Compared with the results obtained with only the seven monitoring stations, the improvements for orbits in radial and clocks are 95.7% and 90.5%, respectively. If the Ka-biases are fixed to known (pre-estimated) values, the accuracies of orbits and clocks are further improved; e.g., the overlap clock difference (95%) reduces to 0.054 ns. By jointly processing Ka-pseudoranges observed by satellites and two anchor stations, the average RMSs of overlap orbit differences are 0.017 m in the radial direction and 0.235 m in the 3D position, while the overlap clock difference (95%) is 0.037 ns. It is also shown that under the condition of single one-way ranging link, the accuracy of orbits and clocks obtained by the two modes is still significantly better than that obtained by only using data from monitoring stations. The RMS orbit difference is no more than 0.03 m in R, the overlap clock difference (95%) is no more than 0.3 ns, and the SDs of estimated Ka-biases are less than 0.6 ns. As shown in the paper, the accuracy of clocks obtained by processing Ka-pseudoranges observed at anchor stations and satellites is greater than that obtained by jointly processing ISL data and monitoring station data. This may imply that not all error is modeled correctly; this could be due to the phase center offsets for the L-band and Ka-band signal or the time variation of hardware delay in the L-band and Ka-band equipment.