Near real-time PPP-based monitoring of the ionosphere using dual-frequency GPS/BDS/Galileo data

Ionosphere delay is very important to GNSS observations, since it is one of the main error sources which have to be mitigated even eliminated in order to determine reliable and precise positions. The ionosphere is a dispersive medium to radio signal, so the value of the group delay or phase advance of GNSS radio signal depends on the signal frequency. Ground-based GNSS stations have been used for ionosphere monitoring and modeling for a long time. In this paper we will introduce a novel approach suitable for single-receiver operation based on the precise point positioning (PPP) technique. One of the main characteristic is that only carrier-phase observations are used to avoid particular effects of pseudorange observations. The technique consists of introducing ionosphere ambiguity parameters obtained from PPP filter into the geometry-free combination of observations to estimate ionospheric delays. Observational data from stations that are capable of tracking the GPS/BDS/GALILEO from the International GNSS Service (IGS) Multi-GNSS Experiments (MGEX) network are processed. For the purpose of performance validation, ionospheric delays series derived from the novel approach are compared with the global ionospheric map (GIM) from Ionospheric Associate Analysis Centers (IAACs). The results are encouraging and offer potential solutions to the near real-time ionosphere monitoring.     

Mitigation of receiver biases in ionospheric observables from PPP with ambiguity resolution

PPP (Precise Point Positioning) is a GNSS (Global Navigation Satellite Systems) positioning method that requires SSR (State Space Representation) corrections in order to provide solutions with an accuracy of centimetric level. The so-called RT-PPP (Real-time PPP) is possible thanks to real-time precise SSR products, for orbits and clocks, provided by IGS (International GNSS Service) and its associate analysis centers such as CNES (Centre National d'Etudes Spatiales). CNES SSR products also enable RT-PPP with integer ambiguity resolution. In GNSS related literature, PPP with ambiguity resolution (PPP-AR) in real-time is often referred as PPP-RTK (PPP – Real Time Kinematic). PPP-WIZARD (PPP - With Integer and Zero-difference Ambiguity Resolution Demonstrator) is a software that is made available by CNES. This software is capable of performing PPP-RTK. It estimates slant ionospheric delays and other GNSS positioning parameters. Since ionospheric effects are spatially correlated by GNSS data from active networks, it is possible to model and provide ionospheric delays for any position in the network coverage area. The prior knowledge ionospheric delays can reduce positioning convergence for PPP-RTK users. Real-time ionospheric models could benefit from highly precise ionospheric delays estimated in PPP-AR. In this study, we demonstrate that ionospheric delays obtained throughout PPP-AR estimation are actu ally ionospheric observables. Ionospheric observables are biased by an order of few meters caused by the receiver hardware biases. These biases prohibit the use of PPP-WIZARD ionospheric delays to produce ionospheric models. Receiver biases correction is essential to provide ionospheric delays while using PPP-AR based ionospheric observables. In this contribution, a method was implemented to estimate and mitigate receiver hardware biases influence on slant ionospheric observables from PPP-AR. In order to assess the proposed approach, PPP-AR data from 12 GNSS stations were processed over a two-month period (March and April 2018). A comparison between IGS ionospheric products and PPP-AR based ionospheric observables corrected for receiver biases, resulted in a mean of differences of −39 cm and 51 cm standard deviation. The results are consistent with the accuracy of the IGS ionospheric products, 2–8 TECU, considering that 1 TECU is ~16 cm in L1. In another analysis, a comparison of ionospheric delays from 5 pairs of short baselines GNSS stations found an agreement of 0.001 m in mean differences with 22 cm standard deviation after receiver biases were corrected. Therefore, the proposed solution is promising and could produce high quality (1–2 TECU) slant ionospheric delays. This product can be used in a large variety of modeling approaches, since ionospheric delays after correction are unbiased. These results indicate that the proposed strategy is promising, and could benefit applications that require accuracy of 1–2 TECU (~16–32 cm in L1).     

Performance improvement of real-time PPP ambiguity resolution using a regional integer clock

Obtaining reliable GNSS uncalibrated phase delay (UPD) or integer clock products is the key to achieving ambiguity-fixed solutions for real-time (RT) precise point positioning (PPP) users. However, due to the influence of RT orbit errors, the quality of UPD/integer clock products estimated with a globally distributed GNSS network is difficult to ensure, thereby affecting the ambiguity resolution (AR) performance of RT-PPP. In this contribution, by fully utilising the consistency of orbital errors in regional GNSS network coverage areas, a method is proposed for estimating regional integer clock products to compensate for RT orbit errors. Based on Centre National d’Études Spatiales (CNES) RT precise products, the regional GPS/BDS integer clock was estimated with a CORS network in the west of China. Results showed that the difference between the estimated regional and CNES global integer clock/bias products was generally less than 5 cm for GPS, whereas clock differences of greater than 10 cm were observed for BDS due to the large RT orbit error. Compared with PPP using global integer clock/bias products, the AR performance of PPP using the regional integer clock was obviously improved for four rover stations. For single GPS, the horizontal and vertical accuracies of ambiguity-fixed PPP solutions were improved by 56.2% and 45.3% on average, respectively, whereas improvements of 67.5% and 50.5% in the horizontal and vertical directions, respectively, were observed for the combined GPS/BDS situation. Based on a regional integer clock, the RMS error of a kinematic ambiguity-fixed PPP solution in the horizontal direction could reach 0.5 cm. In terms of initialisation time, the average time to first fix (TTFF) in combined GPS/BDS PPP was shortened from 18.2 min to 12.7 min. In view of the high AR performance realised with the use of regional integer clocks, this method can be applied to scenarios that require high positioning accuracy, such as deformation monitoring.     

Higher order ionospheric error correction in BDS precise orbit determination

The ionospheric error affects the accuracy of the Global Navigation Satellite Systems observation and precise orbit determination. Usually, only the first order ionospheric error is considered, which can be eliminated by the ionospheric-free linear combination observation. But the remaining higher order ionospheric error will affect the accuracy of observations and their applications. In this paper, the influence of the higher order ionospheric error have been studied by using the International Geomagnetic Reference Field 13 and the Global Ionosphere Maps model produced by the Center for Orbit Determination in Europe. Focus on ionospheric error, the experiment of paper at doy 302 in 2019, which show that the second order ionospheric error impacting BeiDou Navigation Satellite System (BDS) B1I and B3I observation is 6.3569 mm and 11.8484 mm, respectively. Whereas, the third order ionospheric error impacting BDS B1I and B3I observation is 0.1734 mm and 0.3977 mm, respectively. Due to the current measurement accuracy of BDS carrier-phase observation can reach 2 mm, the influence of high order ionospheric error on observation should be considered. For BDS precise orbit determination, the orbit overlapping results are indicated that its orbit accuracy can be improved approximately 5 mm with the higher order ionospheric error correction, which is also in agreement with the results of Satellite Laser Ranging in this work.     

Towards high accuracy GNSS real-time positioning with smartphones

The possibility to access undifferenced and uncombined Global Navigation Satellite System (GNSS) measurements on smart devices with an Android operating system allows us to manage pseudorange and carrier-phase measurements to increase the accuracy of real-time positioning. The goal is to perform real-time kinematic network positioning with smartphones, evaluating the positioning accuracy regarding an external mass-market device. The positioning of Samsung Galaxy S8+ and Huawei P10 plus smartphones was performed using a dedicated tool developed by the authors, considering a continuous operating reference station (CORS) network with a mean inter-station distance of about 50?km. The same positioning technique was also applied to an external GNSS low-cost single-frequency receiver (u-blox EVK-M8T) to compare performance between the receiver and antenna embedded in the previous smartphones and this low-cost receiver coupled with a mass-market antenna (Garmin GA38). Attention was also focused on the phase ambiguity resolution, that it is still a challenging aspect for mass-market devices: even if the two smartphones provide slightly different results, the accuracy obtainable today is greater than 60?cm with a precision of few centimetres in real-time, if a CORS network is available. For real-time applications using portable devices, decimetre-level accuracy is sufficient for many applications, such as rapid mapping and search and rescue activities: these results will open new frontiers in terms of real-time positioning with portable low-cost devices.     

Carrier phase prediction method for GNSS precise positioning in challenging environment

Since the signals of global navigation satellite system (GNSS) are blocked frequently in challenging environments, the discontinuous carrier phases seriously affect the application of GNSS precise positioning. To improve the carrier phase continuity, this paper proposes a carrier phase prediction method based on carrier open-loop tracking. In the open-loop tracking mode, the carrier numerically controlled oscillator (NCO) is controlled by the predicted Doppler, but not by the loop filter output. To improve the phase prediction effective time, accurate receiver clock drift estimation is studied in the prediction method. The phase prediction performance is tested on GNSS software receiver. In the phase prediction effective time tests, open-loop processes were set for the tested channel. The test results show that, when some satellite signals are blocked in 15?s, the probability of carrier phase error less than quarter cycles is more than 94%. In the real time kinematic (RTK) positioning tests, some satellite signals are blocked in 10–15?s repeatedly. The test results show that, the carrier phase continuity is basically not affected by the signal interruption, and the RTK can almost keep continuous centimeter-level positioning accuracy without re-fixing the integer ambiguity.     

附加Doppler预报坐标约束的模糊度 动态解算新方法

利用GPS(Global Positioning System)相位观测进行动态定位的主要困难就是法方程不适定问题(秩亏或病态问题),导致模糊度浮动解及其协方差阵不准,以此构造的模糊度搜索范围比较大,使得模糊度的搜索非常困难.提出利用Doppler高精度测速预报近似坐标并对其附加约束,解决法方程求逆中的不适定问题,提高模糊度浮动解及其方差阵的准确性,缩小模糊度的搜索范围,提高模糊度搜索的成功率.结果表明,短基线情况,新方法在模糊度动态解算中取得了明显效果.      

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.     

A real-time combined quality control method for GNSS precise positioning in harsh environments

Global Navigation Satellite System (GNSS) precise positioning can be significantly affected by severe multipath effects and outliers in harsh environments, and highly relies on quality control strategies. Previous studies mainly focus on the posterior residuals to check and exclude the outliers in GNSS observations, limited work emphasizes the combined quality control method considering both the prior and posterior knowledge simultaneously. This paper proposed a real-time combined quality control method to process the multipath effects and outliers in harsh environments simultaneously. Specifically, in the prior stage, a modified multipath processing strategy is proposed for both phase and code observations, then a modified detection, identification, and adaptation (DIA) method considering the maximum times of data snooping is studied in the posterior stage. Two dedicated experiments in real harsh environments were carried out to evaluate the performance of the proposed combined quality control method. For the static experiment, the proposed method exhibits smaller positioning errors, the best positioning accuracy, and the highest availability in this study. Specifically, the proposed method exhibits an improved percentage of 55.4 %, 56.3 %, and 59.7 % for positioning accuracy compared to those without the quality control method in the E, N, and U directions, respectively. Besides, the proposed method can further improve the performance of ambiguity resolution with an improved percentage of 32.2 %. For the kinematic experiment, the three-dimensional positioning accuracy of the proposed method is 0.577 m, which exhibits a 40.0 % improvement compared to those without the quality control method. Also, the proposed method exhibits better performance under relatively strong multipath effects. In this sense, the proposed real-time combined quality control method is highly appreciated in terms of positioning availability, accuracy, and ambiguity resolution for GNSS precise positioning, especially in harsh environments.     

北斗D2导航电文解析和完好性性能分析

为评估北斗D2导航电文对系统空间信号完好性影响，对D2导航电文结构和内容进行研究。首先，通过基本导航信息的解析，给出了确定卫星位置和用户位置过程所需的信息和时间，并对受卫星星历和时钟影响的用户测距误差做了说明，介绍了电离层误差、对流层误差改正模型。然后，通过增强服务信息的解析，论述了完好性信息（RURA、UDRE）及格网电离层延迟的计算方法。最后给出了使用完好性信息、用户局部误差获取用户定位误差保护级的算法。     

A novel predictive algorithm for double difference observations of obstructed BeiDou geostationary earth orbit (GEO) satellites

Transmission link disturbances and device failure cause global navigation satellite system (GNSS) receivers to miss observations, leading to poor accuracy in real-time kinematic (RTK) positioning. Previously described solutions for this problem are influenced by the length of the prediction period, or are unable to account for changes in receiver state because they use information from previous epochs to make predictions. We propose an algorithm for predicting double difference (DD) observations of obstructed BeiDou navigation system (BDS) GEO satellites. Our approach adopts the first-degree polynomial function for predicting missing observations. We introduce a Douglas-Peucker algorithm to judge the state of the rover receiver to reduce the impact of predictive biases. Static and kinematic experiments were carried out on BDS observations to evaluate the proposed algorithm. The results of our navigation experiment demonstrate that RTK positioning accuracy is improved from meter to decimeter level with fixed ambiguity (horizontal?<?2?cm, vertical?<?18?cm). Horizontal accuracy is improved by over 50%, and the vertical accuracies of the results of the static and kinematic experiments are increased by 47% and 27% respectively, compared with the results produced by the classical approach. Though as the baseline becomes longer, the accuracy is weakened, our predictive algorithm is an improvement over existing approaches to overcome the issue of missing data.     

北斗三频相位观测值线性组合研究

以北斗三频数据为基础,给出了北斗三频组合观测值的数学模型并系统分析了组合后的波长、电离层误差以及观测噪声误差;在保持整周模糊度整数特性的前提下,以整数线性变换法为基础,通过Matlab编程实现组合系数的优化选取。通过将不同组合系数应用于伪距相位法三频数据周跳探测中,发现组合系数之和不为0的组合周跳探测值均在0.5周以上,且组合系数之和越大,波动越大;通过将不同组合系数应用于无几何CIR法三频数据模糊度解算中,发现系数之和不为0,但波长较长的组合求得的模糊度残差值大致在0.6周,通过合适的模糊度搜索方法可以得到正确的模糊度固定解,在短基线模糊度解算中可以考虑使用。     

Evaluation of a regional real-time precise positioning system based on GPS/BeiDou observations in Australia

The performance of real-time (RT) precise positioning can be improved by utilizing observations from multiple Global Navigation Satellite Systems (GNSS) instead of one particular system. Since the end of 2012, BeiDou, independently established by China, began to provide operational services for users in the Asia-Pacific regions. In this study, a regional RT precise positioning system is developed to evaluate the performance of GPS/BeiDou observations in Australia in providing high precision positioning services for users. Fixing three hourly updated satellite orbits, RT correction messages are generated and broadcasted by processing RT observation/navigation data streams from the national network of GNSS Continuously Operating Reference Stations in Australia (AUSCORS) at the server side. At the user side, RT PPP is realized by processing RT data streams and the RT correction messages received. RT clock offsets, for which the accuracy reached 0.07 and 0.28?ns for GPS and BeiDou, respectively, can be determined. Based on these corrections, an accuracy of 12.2, 30.0 and 45.6?cm in the North, East and Up directions was achieved for the BeiDou-only solution after 30 min while the GPS-only solution reached 5.1, 15.3 and 15.5?cm for the same components at the same time. A further improvement of 43.7, 36.9 and 45.0 percent in the three directions, respectively, was achieved for the combined GPS/BeiDou solution. After the initialization process, the North, East and Up positioning accuracies were 5.2, 8.1 and 17.8?cm, respectively, for the BeiDou-only solution, while 1.5, 3.0, and 4.7?cm for the GPS-only solution. However, we only noticed a 20.9% improvement in the East direction was obtained for the GPS/BeiDou solution, while no improvements in the other directions were detected. It is expected that such improvements may become bigger with the increasing accuracy of the BeiDou-only solution.     

载波相位/捷联惯导组合对高轨飞行器的导航

高轨飞行器可用卫星数目较少,信号空间链路损耗大,使用伪距进行测量的精度较低.提出基于GNSS(Global Navigation Satellite System)卫星载波相位与捷联惯导紧组合方法对高轨机动飞行器进行自主导航.该方法将连续跟踪的卫星初次可用时的整周模糊度的浮点解作为状态变量,通过平方根UKF建立了组合导航非线性滤波模型,提出了基于整周浮点解交集的滤波器故障检测方法.研究表明,提出的组合导航方法充分利用了载波相位高测量精度和系统性误差缓变的特点,提高了系统的可靠性和精度.     

Minimum of PDOP and its applications in inter-satellite links (ISL) establishment of Walker-δ constellation

Within the next decade, there will be a number of GNSS (Global Navigation Satellite System) available, i.e. modernized GPS, Galileo, restored GLONASS, BeiDou and many other regional GNSS augmentation systems. Thus, measurement redundancies and geometry of the satellites can be improved. GDOP (Geometric Dilution of Precision) and PDOP (Position Dilution of Precision) are associated with the constellation geometry of satellites, and they are the geometrically determined factors that describe the effect of geometry on the relationship between measurement error and position error. GDOP and PDOP are often used as standards for selecting good satellites to meet the desired positioning precision. In this paper, the related conclusions of minimum of GDOP which was discussed are given, and it is used to study the minimum of PDOP for two cases that the receiver is on the earth’s surface and the receiver is on satellite. The corresponding theorem and constructive solutions of minimum of PDOP are given. Then, the rationality of the ISL (inter-satellite link) establishment criteria in Walker-δ constellation is discussed by using the theory of minimum of PDOP. Finally, the minimum of PDOP is calculated when the number of satellites is 4–10, and these results are verified by using Monte Carlo method.     

基于MLAMBDA的BDS/GPS单频RTK性能分析

GNSS RTK技术以其高精度、高效率、实时性的优点,被广泛应用于航空航天等领域.目前双频RTK技术已非常成熟并且应用较广.相比于双频,单频GNSS RTK在数据质量控制、定位误差处理等方面存在难点.因此单频RTK服务精度可能会受到限制,其定位性能有待研究.本文基于扩展卡尔曼滤波模型,通过MLAMBDA模糊度搜索方法和Ratio检验法,结合实测数据,对比分析BDS,GPS,BDS/GPS三种模式下的单频RTK定位性能.实验证明在静态场景下,三种模式的单频RTK定位精度都在厘米级,可满足高精度定位需求;动态场景下三种模式的模糊度固定率都在70%以上,可满足日常定位需求.在静态及动态应用场景下,北斗的模糊度固定率最高,模糊度解算所用时间短,能实现快速RTK定位.      

基于NTCM-BC模型的全球卫星导航系统单频电离层延迟修正

选择NTCM-BC模型作为单频电离层延迟修正模型,通过非线性最小二乘拟合的方法,利用提前一天预测的电离层图（COPG文件）,计算得到NTCM-BC模型修正系数;利用Klobuchar模型和IGS发布的GIM数据对NTCM-BC模型进行比较和分析.对太阳活动高、中、低年实测数据的分析结果表明:全球平均水平上,NTCM-BC模型的电离层延迟修正性能明显优于Klobuchar模型,NTCM-BC模型的TEC平均误差和均方根误差比Klobuchar模型分别下降了41%和30%;模型的TEC计算误差与太阳活动剧烈程度成正相关,即太阳活动高年模型误差较大,太阳活动低年误差相对较低.相较于磁静日,磁扰日期间Klobuchar模型和NCTM模型的误差均有一定程度的增加.此外,模型的电离层修正误差同时存在明显的纬度、季节和地方时差异.      

Estimation of QZSS differential code biases using QZSS/GPS combined observations from MGEX

As an important error source in Global Navigation Satellite System (GNSS) positioning and ionospheric modeling, the differential code biases (DCB) need to be estimated accurately, e.g., the regional Quasi-Zenith satellite system (QZSS). In this paper, the DCB of QZSS is estimated by adopting the global ionospheric modeling method based on QZSS/GPS combined observations from Multi-GNSS experiment (MGEX). The performance of QZSS satellite and receiver DCB is analyzed with observations from day of year (DOY) 275–364, 2018. Good agreement between our estimated QZSS satellite DCB and the products from DLR and CAS is obtained. The bias and root mean square (RMS) of DCB are mostly within ±0.3 ns. The day-to-day fluctuation of the DCB time series is less than 0.5 ns with about 96% of the cases for all satellites. However, the receiver DCB is a little less stable than satellite DCB, and their standard deviations (STDs) are within 1.9 ns. The result shows that the stability of the receiver DCBs is not significantly related to the types of receiver or antenna.     

Initial accuracy and reliability of current BDS-3 precise positioning,velocity estimation,and time transfer (PVT)

The primary system of Chinese global BeiDou satellite system (BDS-3) was completed to provide global services on December 27, 2018; this was a key milestone in the development process for BDS in terms of its provision of global services. Therefore, this study analyzed the current performance of BDS-3, including its precise positioning, velocity estimation, and time transfer (PVT). The datasets were derived from international GNSS monitoring and assessment system (iGMAS) tracking networks and the two international time laboratories in collaboration with the International Bureau of Weights and Measures (BIPM). With respect to the positioning, the focus is on the real-time kinematic (RTK) positioning and precise point positioning (PPP) modes with static and kinematic scenarios. The results show that the mean available satellite number is 4.8 for current BDS-3 system at short baseline XIA1–XIA3. The RTK accuracy for three components is generally within cm level; the 3D mean accuracy is 8.9 mm for BDS-3 solutions. For the PPP scenarios, the convergence time is about 4 h for TP01 and BRCH stations in two scenarios. After the convergence, the horizontal positioning accuracy is better than cm level and the vertical accuracy nearly reaches the 1 dm level. With respect to kinematic scenarios, the accuracy stays at the cm level for horizontal components and dm level for the vertical component at two stations. In terms of velocity estimation, the horizontal accuracy stays at a sub-mm level, and the vertical accuracy is better than 2 mm/s in the BDS-3 scenario, even in the Arctic. In terms of time and frequency transfer, the noise level of BDS-3 time links can reach 0.096 ns for long-distances link NT01–TP02 and 0.016 ns for short-distance links TP01–TP02. Frequency stability reaches 5E–14 accuracy when the averaging time is within 10,000 s for NT01–TP02 and 1E–15 for TP01–TP02.     

Fast calibration of between-receiver GPS/Galileo/BDS differential phase bias with millimeter accuracy based on short baseline mode

The differential code and phase biases induced by the receiver hardware (including receiver, antenna, firmware, etc.) of the Global Navigation Satellite System (GNSS) have significant effects on precise timing and ionosphere sensing, thus deserve careful treatment. In this contribution, we propose an approach to fast fix the single-difference ambiguity to finally obtain the unbiased estimates of between-receiver differential phase bias (BR-DPB) and between-receiver differential code-phase bias (BR-DCPB) based on the short baseline mode. The key to this method is that the error sources can be significantly eliminated due to the length of the baseline is very short. At the same time, the empirical constraints and random characteristics of BR-DPB/BR-DCPB were considered, which is conducive to the resolution of single-difference ambiguity. Several sets of GNSS data (GPS L1/L2, Galileo E1/E5b, and BDS B1/B3), recorded by the short baselines in an interval of 30 s and covered a broad range of receiver/antenna types (JAVA, SEPT, LEIC, and TRIM), were used to verify the effectiveness of the proposed method. The numerical tests show that the proposed method is capable of fast fixing the single-difference ambiguity successfully within a few epochs and then providing the unbiased estimates of BR-DPB and BR-DCPB in an epoch-by-epoch manner. Experiments show that the estimated BR-DPB is in millimeter accuracy, which is of great significance for the millimeter-accuracy phase time transfer and ionospheric delay estimation. Furthermore, the calibrated BR-DPB/BR-DCPB can be treated as the known products for long-distance precise timing and ionosphere sensing based on the inter-station single-difference model.