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.     

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.     

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.     

Tightly combined GPS/Galileo RTK for short and long baselines: Model and performance analysis

To ensure the compatibility and interoperability with modernized GPS, Galileo satellites are capable of broadcasting navigation signals on carrier phase frequencies that overlap with GPS, i.e., GPS/Galileo L1-E1/L5-E5a. Moreover, the GPS/Galileo L2-E5b signals have different frequencies with wavelength differences smaller than 4.2?mm. Such overlapping and narrowly spaced signals between GPS and Galileo bring the opportunity to use the tightly combined double-differenced (DD) model for precise real-time kinematic (RTK) positioning, resulting in improved performance of ambiguity resolution and positioning with respect to the classical standard or loosely combined DD model. In this paper, we focus on the model and performance assessment of tightly combined GPS/Galileo L1-E1/L2-E5b/L5-E5a RTK for short and long baselines. We first investigate the tightly combined GPS/Galileo DD observational model for both short and long baselines with simultaneously considering the GPS/Galileo overlapping and non-overlapping frequencies. Particularly, we introduce a reparameterization approach to solve the rank deficiency that caused by the correlation between the DISB parameters and the DD ionospheric parameters for both overlapping and non-overlapping frequencies. Then we present performance assessment for the tightly combined GPS/Galileo RTK model with real-time estimation of the differential inter-system bias (DISB) parameters for short and long baselines in terms of ratio value, ambiguity dilution of precision (ADOP), ambiguity conditional number, decorrelation number, search count, empirical success rate, time-to-first-fix (TTFF), and positioning accuracy. Results from both static and kinematic experiments demonstrated that compared to the loosely combined model, the tightly combined model can deliver improved performance of ambiguity resolution and precise positioning with different satellite visibility. For the car-driven short baseline experiment with 10° elevation cut-off angle, the tightly combined model can not only significantly increase the ratio value by approximately 27.5% (from 16.0 to 20.4), but also reduce the ambiguity ADOP, the conditional number, and the search count in LAMBDA by approximately 22.2% (from 0.027 to 0.021 cycles), 14.9% (from 199.2 to 169.6), and 25.4% (from 150.1 to 112.0), respectively. Comparable decorrelation number, empirical success rate, and positioning accuracy are also obtained. For the car-driven long baseline experiment, it is also observed that the ambiguity resolution performance in terms of the ratio value, the decorrelation number, the condition number, and the search count are significantly improved by approximately 18.5% (from 2.7 to 3.2), 22.0% (from 0.186 to 0.227), 55.9% (from 937.6 to 413.7), and 10.3% (from 43.8 to 39.3), respectively. Moreover, comparable ADOP, empirical success rate, and positioning accuracy are obtained as well. Additionally, the TTFF can be reduced (from 54.1 to 51.8 epochs with 10° elevation cut-off angle) as well from the results of static experiments.     

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).     

Analysis of the code and phase between-receiver inter-system biases of the overlapping frequencies for GPS/Galileo/BDS

The overlapping-frequency signals from different GNSS constellations are interoperable and can be integrated as one constellation in multi-GNSS positioning when inter-system bias (ISB) is properly disposed. The look-up table method for ISB calibration can enhance the model strength, maximize the number of integer-estimable ambiguities, and thus is preferred. However, the characteristics and magnitudes of the receiver code ISB and phase fractional ISB (F-ISB) are not well known and the wrong values of the biases can seriously degrade the positioning results. In this contribution, we first estimate the between-receiver code ISB and phase F-ISB of hundreds of the baselines up to around 25km in the European Permanent GNSS Network (EPN) and the Multi-GNSS Experiment (MGEX) for the overlapping frequencies L1-E1 (L1), L5-E5a (L5) and E5b-B2b (L7). The data collected from 1st January 2016 to 1st January 2019. Second, the receiver-type and firmware-version combinations for the receivers of Trimble, Leica, Javad, Septentrio and NovAtel are carefully classified. Results show that the Septentrio receivers have consistent code and phase ISB values for the three overlapping frequencies i.e. only one value for each frequency and no receivers are different. The Leica, Trimble and Javad receivers have two or more ISB values for at least one of the three frequencies. A few receivers with biases to the groups are also found and listed. Third, the code ISB and phase F-ISB of the groups are adjusted by the least-squares method. The root mean square errors (RMSE) of the least square adjustment are 0.240 m, 0.250 m and 0.200 m for code of L1, L5 and L7 frequencies, respectively, and are 0.0009 m, 0.0015 m and 0.0031 m for phase of L1, L5 and L7 frequencies, respectively. Finally, the effects of code ISB errors on code positing are investigated with the zero-baseline MAT1_MATZ. The distance root mean square error (DRMS) of L1-E1 code positioning can be reduced by 48.2% with 5 GPS and Galileo satellites and the DRMS degrades quickly when the code ISB error is larger.     

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.     

Total electron content monitoring using triple frequency GNSS: Results with Giove-A/-B data

Triple frequency GNSS will be fully operational within the next decade, opening opportunities for new applications. Dual frequency GNSS already allow to study the ionosphere through the estimation of Total Electron Content (TEC). However, the precision is limited by the ambiguity resolution process. This paper studies a triple frequency TEC monitoring technique in which the use of Geometry-Free and Iono-Free linear combinations improves the ambiguity resolution process and therefore the precision of TEC. We have tested it on a set of triple frequency Giove-A/-B data from January and December 2008. The conclusions achieved are (1) TEC values are affected by an error of about 2–2.5 TECU produced through the ambiguity resolution process; (2) the error caused by the Geometric Free phase combination delays (hardware, multipath, noise, antenna phase center) on TEC is about 0.2 TECU; (3) the total error on TEC approximately reach 2–3 TECU.     

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).     

Tight integration of BDS-3/BDS-2/GPS/Galileo observations considering the new overlapping DISBs and its application in obstructed environments

Tight integration can enhance the model strength and positioning performance by considering the characteristic of differential inter-system bias (DISB), especially in obstructed environments. However, limited work emphasizes the comprehensive analysis of five-frequency DISBs between BDS-3 and other systems considering the receiver type, receiver configuration, and antenna type. In addition, the overlapping DISBs between BDS-3 and BDS-2 are also in great demand for further investigation since they are often regarded as one system. In this study, one DISB-float model is introduced to estimate the DISBs, and one DISB-fixed model and one DISB-free model are formulated to enhance the model strength of tight integration. Four dedicated datasets were collected to estimate the DISBs, which are also comprehensively analyzed considering the receiver type, receiver configuration, and antenna type. The results show that the DISBs between BDS-3 and other systems are rather stable over a certain period and are related to the receiver type and receiver configuration, whereas are not related to the antenna type. More interestingly, the B1I code DISB between BDS-3 and BDS-2 exhibits significant magnitude with a mean value of ?1.44 m for the baseline composed of two different receivers. In this case, the B1I code DISB must be considered and the tight integration between BDS-3 and BDS-2 considering its calibration can improve the positioning performance. Besides, the tight integration of the DISB-fixed model can significantly improve the positioning accuracy between multiple GNSS. Compared to the loose integration, the improvement of 60.6 %, 56.6 %, and 61.2 % can be obtained in the E, N, and U directions, when only two satellites are available for each system. In real obstructed environments, the tight integration of the DISB-free model can also improve the positioning performance in terms of positioning availability and accuracy, as well as the ambiguity resolution performance.     

基于GNSS的高轨卫星定位技术研究

利用全球卫星导航系统(GNSS)进行导航定位具有全球、全天候、实时和高精度的优点,应用于高地球轨道(HEO)卫星的定位,能够提供精确的轨道和姿态确定,并且可以克服目前主要利用地面测控系统对HEO卫星进行定位的设备复杂、投资高等缺点,使得自主导航成为可能.本文对利用GNSS的高轨卫星定位相关技术进行了研究,分析了单一GNSS系统和多个GNSS组合系统的卫星可见性、动态性和几何精度因子(GDOP).通过仿真分析表明,利用组合GNSS系统并通过提高GNSS接收机灵敏度的方法,可以解决GNSS进行HEO卫星定位的相关问题,并能保证HEO卫星定位精度的要求.      

Low-cost precise measurement of oscillator frequency instability based on GNSS carrier observation

Global navigation satellite systems (GNSS) receivers can be used in time and frequency metrology by exploiting stable GNSS time scales. This paper proposes a low-cost method for precise measurement of oscillator frequency instability using a single-frequency software GNSS receiver. The only required hardware is a common radio frequency (RF) data collection device driven by the oscillator under test (OUT). The receiver solves the oscillator frequency error in high time resolution using the carrier Doppler observation and the broadcast ephemeris from one of the available satellites employing the onboard reference atomic frequency standard that is more stable than the OUT. Considering the non-stable and non-Gaussian properties of the frequency error measurement, an unbiased finite impulse response (FIR) filter is employed to obtain robust estimation and filter out measurement noise. The effects of different filter orders and convolution lengths are further discussed. The frequency error of an oven controlled oscillator (OCXO) is measured using live Beidou-2/Compass signals. The results are compared with the synchronous measurement using a specialized phase comparator with the standard coordinated universal time (UTC) signal from the master clock H226 in the national time service center (NTSC) of China as its reference. The Allan deviation (ADEV) estimates using the two methods have a 99.9% correlation coefficient and a 0.6% mean relative difference over 1–1000 s intervals. The experiment demonstrates the effectiveness and high precision of the software receiver method.     

Testing a new multivariate GNSS carrier phase attitude determination method for remote sensing platforms

GNSS (Global Navigation Satellite Systems)-based attitude determination is an important field of study, since it is a valuable technique for the orientation estimation of remote sensing platforms. To achieve highly accurate angular estimates, the precise GNSS carrier phase observables must be employed. However, in order to take full advantage of the high precision, the unknown integer ambiguities of the carrier phase observables need to be resolved. This contribution presents a GNSS carrier phase-based attitude determination method that determines the integer ambiguities and attitude in an integral manner, thereby fully exploiting the known body geometry of the multi-antennae configuration. It is shown that this integral approach aids the ambiguity resolution process tremendously and strongly improves the capacity of fixing the correct set of integer ambiguities. In this contribution, the challenging scenario of single-epoch, single-frequency attitude determination is addressed. This guarantees a total independence from carrier phase slips and losses of lock, and it also does not require any a priori motion model for the platform. The method presented is a multivariate constrained version of the popular LAMBDA method and it is tested on data collected during an airborne remote sensing campaign.     

Estimating inter-system biases for tightly combined Galileo/BDS/GPS RTK

The overlapping carrier frequencies L1/E1, L5/E5a and B2/E5b from GPS/Galileo/BDS allow inter-system double-differencing of observations, which shows a clear advantage over differencing of the observations of each constellation independently. However, the inter-system biases destroy the integer nature of the inter-system double-differencing ambiguities. Two methods of direct rounding and parameter estimation are used to determine the ISB value. By analyzing data collected from Curtin University from 2015 to 2018, the phase and code inter-system bias (ISB) are related to the receiver type, firmware version and the selected overlapping frequency. Upgrade of receiver firmware version results in changes of ISB values. For example, the upgrade of Javad firmware in Dec, 15, 2017 causes the difference of 0.5 cycles ISB between BDS GEO and non-GEO satellites. By comparing the three dynamic models which include white noise process, random walk process, and random constant in the parameter estimation method, the ISB determined by the random constant model is consistent with the value obtained by the direct rounding method. After the calibration of ISBs, the performances of tightly combined positioning are assessed. The success rate of ambiguity resolution and accuracy of positioning for the tight combination (TC) are significantly improved in comparison with that for the loose combination (LC) over short baselines. For L5/E5a, on which only few satellites can be observed, the maximum increase in success rates of ambiguity resolution can reach 31.7%, i.e., from 54.9% of LC to larger than 86.6% of TC, and the positioning accuracies can even be increased by 0.13 m, i.e., from 0.208 of LC to 0.074 m of TC in East direction for the mix-receiver TRIMBLE NETR9-SEPT POLARX4 in 2018.     

Evaluation of SARAL/AltiKa performance using GNSS/IMU equipped buoy in Sajafi,Imam Hassan and Kangan Ports

Performance of SARAL/AltiKa mission has been evaluated within 2016 altimeter calibration/validation framework in Persian Gulf through three campaigns conducted in the offshore waters of Sajafi, Imam Hassan and Kangan Ports, while the altimeter overflew the passes 470, 111 and 25 on 13 Feb, 7 March and 17 June 2016, respectively. As the preparation, a lightweight buoy was equipped with a GNSS receiver/choke-ring antenna and a MEMS-based IMU to measure independent datasets in the field operations. To obtain accurate sea surface height (SSH) time series, the offset of the onboard antenna from the equilibrium sea level was predetermined through surveying operations as the buoy was deploying in the onshore waters of Kangan Port. Accordingly, the double-difference carrier phase observations have been processed via the Bernese GPS Software v. 5.0 so as to provide the GNSS-derived time series at the comparison points of the calibration campaigns, once the disturbing effects due to the platform tilt and heave have been eliminated. Owing to comparing of the SSH time series and the associating altimetry 1?Hz GDR-T datasets, the calibration/validation of the SARAL/AltiKa has been performed in the both cases of radiometer and ECMWF wet troposphere corrections so as to identify potential land contamination. An agreement of the present findings in comparison with those attained in other international calibrations sites confirms the promising feasibility of Persian Gulf as a new dedicated site for calibration/validation of ongoing and future altimetry missions.     

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

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

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.     

Demonstration of a high sensitivity GNSS software receiver for indoor positioning

Advances in signal processing techniques contributed to the significant improvements of GNSS receiver performance in dense multipath environments and created the opportunities for a new category of high-sensitivity GNSS (HS-GNSS) receivers that can provide GNSS location services in indoor environments. The difficulties in improving the availability, reliability, and accuracy of these indoor capable GNSS receivers exceed those of the receivers designed for the most hostile urban canyon environments. The authors of this paper identified the vector tracking schemes, signal propagation statistics, and parallel processing techniques that are critical to a robust HS-GNSS receiver for indoor environments and successfully incorporated them into a fully functional high-sensitivity software receiver. A flexible vector-based receiver architecture is introduced to combine these key indoor signal processing technologies into GSNRx-hs™ – the high sensitivity software navigation receiver developed at the University of Calgary. The resulting receiver can perform multi-mode vector tracking in indoor environment at various levels of location and timing uncertainties. In addition to the obvious improvements in time-to-first-fix (TTFF) and signal sensitivity, the field test results in indoor environments surrounded by wood, glass, and concrete showed that the new techniques effectively improved the performance of indoor GNSS positioning. With fine GNSS timing, the proposed receiver can consistently deliver indoor navigation solution with the horizontal accuracy of 2–15 m depending on the satellite geometry and the indoor environments. If only the coarse GNSS timing is available, the horizontal accuracy of the indoor navigation solution from the proposed receiver is around 30 m depending on the coarse timing accuracy, the satellite geometry, and the indoor environments. From the preliminary field test results, it has been observed that the signal processing sensitivity is the dominant factor on the availability of the indoor navigation solution, while the GNSS timing accuracy is the dominant factor on the accuracy of the indoor navigation solution.     

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.