Estimation of differential code biases with Jason-2/3 onboard GPS observations

With the continuous deployment of Low Earth Orbit (LEO) satellites, the estimation of differential code biases (DCBs) based on GNSS observations from LEO has gained increasing attention. Previous studies on LEO-based DCB estimation are usually using the spherical symmetry ionosphere assumption (SSIA), in which a uniform electron density is assumed in a thick shell. In this study, we propose an approach (named the SHLEO method) to simultaneously estimate the satellite and LEO onboard receiver DCBs by modeling the distribution of the global plasmaspheric total electron content (PTEC) above the satellite orbit with a spherical harmonic (SH) function. Compared to the commonly used SSIA method, the SHLEO model improves the GPS satellite DCB estimation accuracy by 13.46% and the stability by 22.34%, respectively. Compared to the GPS satellite DCBs estimated based on the Jason-3-only observations, the accuracy and monthly stability of the satellite DCBs can be improved by 14.42% and 26.8% when both Jason-2 and Jason-3 onboard observations are jointly processed. Compared with the Jason-2 solutions, the GPS satellite DCB estimates based on the fusion of Jason-2 and Jason-3 observations have an improved consistency of better than 18.26% and 9.71% with the products provided by the Center for Orbit Determination in Europe (CODE) and Chinese Academy of Sciences (CAS). Taking the DCB products provided by the German Aerospace Center (DLR) as references, there is no improvement in accuracy of the GPS satellite DCB estimates based on the fusion of Jason-2 and Jason-3 observations than the Jason-2 solutions alone. A periodic variation is found in the time series of both the Jason-3 and Jason-2 onboard receiver DCB estimates. Preliminary analysis of the PTEC distribution based on the estimated SH coefficients are also presented.     

一种地基GNSS接收机差分码偏差估算方法

GNSS不同频点间的码伪距作差会引入信号的差分码偏差(DCB),包括GNSS卫星及地面接收机的DCB.本文提出一种地基GNSS接收机差分码偏差参数估算方法,首先由电离层文件参数作线性插值,计算出电离层延迟误差.之后对IGS站观测文件进行加权最小二乘法估计,得到GPS卫星和地面GNSS接收机的L1C频点和L2P频点间码偏...     

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.     

电离层周日变化对解算GPS硬件延迟稳定性的影响

针对电离层周日变化特征分析了其可能对SCORE方法估算的硬件延迟稳定性的影响. 利用BJFS以及XIAM台站的GPS观测数据, 解算了位于太阳活动高年(2001年)和太阳活动低年(2009年)的卫星硬件延迟并分析了估算的硬件延迟的稳定性. 研究发现, 电离层周日变化对估算的硬件延迟稳定性具有一定影响, 但是利用不同台站所得到的卫星硬件延迟稳定性在昼夜不同时间上的解算结果存在一定差异. 电离层周日变化对利用 BJFS台站数据解算的硬件延迟稳定性日夜差异较为明显, 在太阳活动高年利用XIAM 台站数据解算的硬件延迟日夜稳定性差异不很明显, 由于XIAM台站处于电离层赤道异常峰附近, 夜间电离层变化很大, 因此对比中纬度地区, 电离层周日变化对赤道异常峰附近地区硬件延迟稳定性解算结果的影响相对较小, 但在太阳活动低年, 其影响仍较为显著.      

基于三频数据的北斗卫星导航系统DCB参数精度评估方法

差分码偏差（Differential Code Biases,DCB）参数作为导航电文中重要的一项,是影响用户PNT服务的主要误差源之一。北斗卫星导航系统（后文简称“北斗系统”）发射三个频点的导航信号,在导航电文中需要发播卫星的2个TGD（Timing Group Delay）参数。文章首先介绍了北斗系统卫星DCB参数最小二乘解算与形式误差评估;其次根据北斗系统三频特点,提出了不同频点组合计算垂直方向电离层电子总含量（VTEC）互差的DCB精度定量评估方法,并与IGS(International GNSS Service)提供的GPS卫星DCB精度进行比较;最后,详细分析了DCB参数精度对用户等效距离误差（UERE）计算和定位计算的影响,分别采用卫星出场标定DCB参数和经过解算DCB参数进行评估。实测数据分析结果表明,北斗系统卫星DCB参数解算形式误差与IGS解算GPS卫星DCB参数形式误差相当,但受卫星类型和解算测站的几何分布限制,北斗系统卫星DCB参数解算不确定度相比IGS略差,估计精度优于0.5ns,不同频率组合计算VTEC互差绝对值均值优于0.6TECU。相比采用卫星出场标定值,采用系统解算DCB参数后,双频用户三维位置误差改善13.80%～47.42%。     

Multi-GNSS contributions to differential code biases determination and regional ionospheric modeling in China

Due to the limited number and uneven distribution globally of Beidou Satellite System (BDS) stations, the contributions of BDS to global ionosphere modeling is still not significant. In order to give a more realistic evaluation of the ability for BDS in ionosphere monitoring and multi-GNSS contributions to the performance of Differential Code Biases (DCBs) determination and ionosphere modeling, we select 22 stations from Crustal Movement Observation Network of China (CMONOC) to assess the result of regional ionospheric model and DCBs estimates over China where the visible satellites and monitoring stations for BDS are comparable to those of GPS/GLONASS. Note that all the 22 stations can track the dual- and triple-frequency GPS, GLONASS, and BDS observations. In this study, seven solutions, i.e., GPS-only (G), GLONASS-only (R), BDS-only (C), GPS + BDS (GC), GPS + GLONASS (GR), GLONASS + BDS (RC), GPS + GLONASS + BDS (GRC), are used to test the regional ionosphere modeling over the experimental area. Moreover, the performances of them using single-frequency precise point positioning (SF-PPP) method are presented. The experimental results indicate that BDS has the same ionospheric monitoring capability as GPS and GLONASS. Meanwhile, multi-GNSS observations can significantly improve the accuracy of the regional ionospheric models compared with that of GPS-only or GLONASS-only or BDS-only, especially over the edge of the tested region which the accuracy of the model is improved by reducing the RMS of the maximum differences from 5–15 to 2–3 TECu. For satellite DCBs estimates of different systems, the accuracy of them can be improved significantly after combining different system observations, which is improved by reducing the STD of GPS satellite DCB from 0.243 to 0.213, 0.172, and 0.165 ns after adding R, C, and RC observations respectively, with an increment of about 12.3%, 29.4%, and 32.2%. The STD of GLONASS satellite DCB improved from 0.353 to 0.304, 0.271, and 0.243 ns after adding G, C, and GC observations, respectively. The STD of BDS satellite DCB reduced from 0.265 to 0.237, 0.237 and 0.229 ns with the addition of G, R and GR systems respectively, and increased by 10.6%, 10.4%, and 13.6%. From the experimental positioning result, it can be seen that the regional ionospheric models with multi-GNSS observations are better than that with a single satellite system model.     

Estimation and analysis of GPS receiver differential code biases using KGN in Korean Peninsula

The total electron content (TEC) estimation by the Global Positioning System (GPS) can be seriously affected by the differential code biases (DCB), referred to as inter-frequency biases (IFB), of the satellite and receiver so that an accuracy of GPS–TEC value is dependent on the error of DCBs estimation. In this paper, we proposed the singular value decomposition (SVD) method to estimate the DCB of GPS satellites and receivers using the Korean GPS network (KGN) in South Korea. The receiver DCBs of about 49 GPS reference stations in KGN were determined for the accurate estimation of the regional ionospheric TEC. They obtained from the daily solution have large biases ranging from +5 to +27 ns for geomagnetic quiet days. The receiver DCB of SUWN reference station was compared with the estimates of IGS and JPL global ionosphere map (GIM). The results have shown comparatively good agreement at the level within 0.2 ns. After correction of receiver DCBs and knowing the satellite DCBs, the comparison between the behavior of the estimated TEC and that of GIMs was performed for consecutive three days. We showed that there is a good agreement between KASI model and GIMs.     

LEO enhanced Global Navigation Satellite System (LeGNSS) for real-time precise positioning services

Global Navigation Satellite System (GNSS) has been widely used in many geosciences areas with its Positioning, Navigation and Timing (PNT) service. However, GNSS still has its own bottleneck, such as the long initialization period of Precise Point Positioning (PPP) without dense reference network. Recently, the concept of PNTRC (Positioning, Navigation, Timing, Remote sensing and Communication) has been put forward, where Low Earth Orbit (LEO) satellite constellations are recruited to fulfill diverse missions. In navigation aspect, a number of selected LEO satellites can be equipped with a transmitter to transmit similar navigation signals to ground users, so that they can serve as GNSS satellites but with much faster geometric change to enhance GNSS capability, which is named as LEO constellation enhanced GNSS (LeGNSS). As a result, the initialization time of PPP is expected to be shortened to the level of a few minutes or even seconds depending on the number of the LEO satellites involved. In this article, we simulate all the relevant data from June 8th to 14th, 2014 and investigate the feasibility of LeGNSS with the concentration on the key issues in the whole data processing for providing real-time PPP service based on a system configuration with fourteen satellites of BeiDou Navigation Satellite System (BDS), twenty-four satellites of the Global Positioning System (GPS), and sixty-six satellites of the Iridium satellite constellations. At the server-end, Precise Orbit Determination (POD) and Precise Clock Estimation (PCE) with various operational modes are investigated using simulated observations. It is found out that GNSS POD with partial LEO satellites is the most practical mode of LeGNSS operation. At the user-end, the Geometry Dilution Of Precision (GDOP) and Signal-In-Space Ranging Error (SISRE) are calculated and assessed for different positioning schemes in order to demonstrate the performance of LeGNSS. Centimeter level SISRE can be achieved for LeGNSS.     

A method of estimating GPS instrumental biases with a convolution algorithm

This paper presents a method of deriving the instrumental differential code biases (DCBs) of GPS satellites and dual frequency receivers. Considering that the total electron content (TEC) varies smoothly over a small area, one ionospheric pierce point (IPP) and four more nearby IPPs were selected to build an equation with a convolution algorithm. In addition, unknown DCB parameters were arranged into a set of equations with GPS observations in a day unit by assuming that DCBs do not vary within a day. Then, the DCBs of satellites and receivers were determined by solving the equation set with the least-squares fitting technique. The performance of this method is examined by applying it to 361?days in 2014 using the observation data from 1311 GPS Earth Observation Network (GEONET) receivers. The result was crosswise-compared with the DCB estimated by the mesh method and the IONEX products from the Center for Orbit Determination in Europe (CODE). The DCB values derived by this method agree with those of the mesh method and the CODE products, with biases of 0.091?ns and 0.321?ns, respectively. The convolution method's accuracy and stability were quite good and showed improvements over the mesh method.     

Calibration errors in determining slant Total Electron Content (TEC) from multi-GNSS data

The global navigation satellite system (GNSS) is presently a powerful tool for sensing the Earth's ionosphere. For this purpose, the ionospheric measurements (IMs), which are by definition slant total electron content biased by satellite and receiver differential code biases (DCBs), need to be first extracted from GNSS data and then used as inputs for further ionospheric representations such as tomography. By using the customary phase-to-code leveling procedure, this research comparatively evaluates the calibration errors on experimental IMs obtained from three GNSS, namely the US Global Positioning System (GPS), the Chinese BeiDou Navigation Satellite System (BDS), and the European Galileo. On the basis of ten days of dual-frequency, triple-GNSS observations collected from eight co-located ground receivers that independently form short-baselines and zero-baselines, the IMs are determined for each receiver for all tracked satellites and then for each satellite differenced for each baseline to evaluate their calibration errors. As first derived from the short-baseline analysis, the effects of calibration errors on IMs range, in total electron content units, from 1.58 to 2.16, 0.70 to 1.87, and 1.13 to 1.56 for GPS, Galileo, and BDS, respectively. Additionally, for short-baseline experiment, it is shown that the code multipath effect accounts for their main budget. Sidereal periodicity is found in single-differenced (SD) IMs for GPS and BDS geostationary satellites, and the correlation of SD IMs over two consecutive days achieves the maximum value when the time tag is around 4?min. Moreover, as byproducts of zero-baseline analysis, daily between-receiver DCBs for GPS are subject to more significant intra-day variations than those for BDS and Galileo.     

A MARS-based method for estimating regional 2-D ionospheric VTEC and receiver differential code bias

The geometry-free linear combination of dual-frequency GNSS reference station ground observations are currently used to build the Vertical Total Electron Content (VTEC) model of the ionosphere. As it is known, besides ionospheric delays, there are differential code bias (DCB) of satellite (SDCB) and receiver (RDCB) in the geometry-free observation equation. The SDCB can be obtained using the International GNSS Service (IGS) analysis centers, but the RDCB for regional and local network receivers are not provided. Therefore, estimating the RDCB and VTEC model accurately and simultaneously is a critical factor investigated by researchers. This study uses Multivariate Adaptive Regression Splines (MARS) to estimate the VTEC approximate model and then substitutes this model in the observation equation to form the normal equation. The least squares method is used to solve the RDCB and VTEC model together. The research findings show that this method has good modeling effectiveness and the estimated RDCB has good reliability. The estimated VTEC model applied to GPS single-frequency precise point positioning has better positioning accuracy in comparison to the IGS global ionosphere map (GIM).     

星载GPS观测数据的模拟研究

重点分析了星载GPS观测值模拟的原理、数学模型等，并用模拟软件模拟了CHAMP卫星的星载GPS观测值情况．结果表明，根据本文提供模拟方法所模拟的CHAMP星载GPS观测值，与实测值相比，无论是在大小还是观测噪声水平上都很接近，因此能满足不同层次人员对星载GPS模拟观测值的需要．     

Regional 4-D modeling of the ionospheric electron density from satellite data and IRI

Accurate knowledge of the electron density is the key point in correcting ionospheric delays of electromagnetic measurements and in studying ionosphere physics. During the last decade Global Navigation Satellite Systems (GNSS) have become a promising tool for monitoring ionospheric parameters such as the total electron content (TEC). In this contribution we present a four-dimensional (4-D) model of the electron density consisting of a given reference part, i.e., the International Reference Ionosphere (IRI), and an unknown correction term expanded in terms of multi-dimensional base functions. The corresponding series coefficients are calculable from the satellite measurements by applying parameter estimation procedures. Since satellite data are usually sampled between GPS satellites and ground stations, finer structures of the electron density are modelable just in regions with a sufficient number of ground stations. The proposed method is applied to simulated geometry-free GPS phase measurements. The procedure can be used, for example, to study the equatorial anomaly.     

A real-time ionospheric model based on GNSS Precise Point Positioning

This paper proposes a method of real-time monitoring and modeling the ionospheric Total Electron Content (TEC) by Precise Point Positioning (PPP). Firstly, the ionospheric TEC and receiver’s Differential Code Biases (DCB) are estimated with the undifferenced raw observation in real-time, then the ionospheric TEC model is established based on the Single Layer Model (SLM) assumption and the recovered ionospheric TEC. In this study, phase observations with high precision are directly used instead of phase smoothed code observations. In addition, the DCB estimation is separated from the establishment of the ionospheric model which will limit the impacts of the SLM assumption impacts. The ionospheric model is established at every epoch for real time application. The method is validated with three different GNSS networks on a local, regional, and global basis. The results show that the method is feasible and effective, the real-time ionosphere and DCB results are very consistent with the IGS final products, with a bias of 1–2 TECU and 0.4 ns respectively.     

单接收机GNSS数据组合硬件延迟的联合求解方法

利用GNSS观测数据解算TEC的最大误差源是硬件延迟,包括卫星硬件延迟和接收机硬件延迟.在单接收机情况下,由于数据稀疏以及接收到的卫星信号时间不对齐等特点,已有的解算硬件延迟方法的求解结果往往不理想.在应用局域模式拟合方法和SCORE方法求解单接收机数据基础上,利用局域模型拟合法在电离层平静期拟合较准确的优点,提出一种联合改进方法,同时改正了SCORE方法解算过程中约束过强的缺点.通过利用GPStation-6接收机的GPS和BDS实际观测数据进行解算分析,验证了所提方法的有效性与准确性.      

NeQuick G model based scale factor determination for using SBAS ionosphere corrections at low earth orbit

A space-based augmentation system (SBAS) provides real-time correction data for global navigation satellite system (GNSS) users near ground. In order to use the SBAS ionosphere correction for low Earth orbit (LEO) satellites, the correction should be scaled down for the LEO altitude. This scale factor varies with ionosphere distribution and it is hard to determine the value at LEO in real time. We propose a real-time scale factor determination method by using Galileo GNSS’s NeQuick G model. A LEO satellite GPS data and SBAS data received on ground were used to evaluate the performance of the NeQuick G derived variable scale factor. The NeQuick G derived scale factor shows a significant accuracy improvement over NeQuick G model or pre-determined constant scale factor. It improves a vertical positioning accuracy of the LEO satellite. The error mean reductions of the vertical positioning over NeQuick G and the constant scale factor are 31.5% and 11.7%, respectively.     

The impact on tsunami detection from space using GNSS-reflectometry when combining GPS with GLONASS and Galileo

The devastating Sumatra tsunami in 2004 demonstrated the need for a tsunami early warning system in the Indian Ocean. Such a system has been installed within the German-Indonesian Tsunami Early Warning System (GITEWS) project. Tsunamis are a global phenomenon and for global observations satellites are predestined. Within the GITEWS project a feasibility study on a future tsunami detection system from space has therefore been carried out. The Global Navigation Satellite System Reflectometry (GNSS-R) is an innovative way of using GNSS signals for remote sensing. It uses ocean reflected GNSS signals for sea surface altimetry. With a dedicated Low Earth Orbit (LEO) constellation of satellites equipped with GNSS-R receivers, densely spaced sea surface height measurements could be established to detect tsunamis. Some general considerations on the geometry between LEO and GNSS are made in this simulation study. It exemplary analyzes the detection performance of a GNSS-R constellation at 900 km altitude and 60° inclination angle when applied to the Sumatra tsunami as it occurred in 2004. GPS is assumed as signal source and the combination with GLONASS and Galileo signals is investigated. It can be demonstrated, that the combination of GPS and Galileo is advantageous for constellations with few satellites while the combination with GLONASS is preferable for constellations with many satellites. If all three GNSS are combined, the best detection performance can be expected for all scenarios considered. In this case an 18 satellite constellation will detect the Sumatra tsunami within 17 min with certainty, while it takes 53 min if only GPS is considered.     

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.     

Rapid local ionosphere modeling based on Precise Point Positioning over Iran: A case study under 2014 solar maximum

Presently, the ionosphere effect is the main source of the error in the Global Positioning System (GPS) observations. This effect can largely be removed by using the two-frequency measurements, while to obtain the reasonable results in the single-frequency applications, an accurate ionosphere model is required. Since the global ionosphere models do not meet our needs everywhere, the local ionosphere models are developed. In this paper, a rapid local ionosphere model over Iran is presented. For this purpose, the GPS observations obtained from 40 GPS stations of the Iranian Permanent GPS Network (IPGN) and 16 other GPS stations around Iran have been used. The observations have been selected under 2014 solar maximum, from the days 058, 107, 188 and 271 of the year 2014 with different geomagnetic activities. Moreover, ionospheric observables based on the precise point positioning (PPP) have been applied to model the ionosphere. To represent our ionosphere model, the B-spline basis functions have been employed and the variance component estimation (VCE) method has been used to regularize the problem.To show the efficiency our PPP-derived local ionosphere model with respect to the International GNSS Service (IGS) global models, these models are applied on the single point positioning using single-frequency observations and their results are compared with the precise coordinates obtained from the double-differenced solution using dual-frequency observations. The results show that the 95th percentile of horizontal and vertical positioning errors of the single-frequency point positioning are about 3.1 and 13.6?m, respectively, when any ionosphere model are not applied. These values significantly improve when the ionosphere models are applied in the solutions. Applying CODE’s Rapid Global ionosphere map (CORG), improvements of 59% and 81% in horizontal and vertical components are observed. These values for the IGS Global ionosphere map (IGSG) are 70% and 82%, respectively. The best results are obtained from our local ionosphere model, where 84% and 87% improvements in horizontal and vertical components are observed. These results confirm the efficiency of our local ionosphere model over Iran with respect to the global models. As a by-product, the Differential Code Biases (DCBs) of the receivers are also estimated. In this line, we found that the intra-day variations of the receiver DCBs could be significant. Therefore, these variations must be taken into account for the precise ionosphere modeling.     

粗差拟准检定法在星载GPS低轨卫星定轨中的应用

利用自适应卡尔曼滤波进行星载GPS低轨卫星定轨时,必须解决量测方程中经常存在的粗差问题．在分析以往方法的优缺点后,用拟准检定法来探测和修正量测方程中存在的粗差．该法的优点是辨识粗差准确率高,能同时定位多个粗差．另外,为了克服星载GPS低轨卫星定轨的滤波器可能出现的数值不稳定性及发散现象,还采用了UD分解算法及Sage自适应滤波器．最后用一个CHAMP卫星的模拟算例验证本方法的可行性和有效性．