Evaluation on real-time dynamic performance of BDS in PPP,RTK, and INS tightly aided modes

Since China’s BeiDou satellite navigation system (BDS) began to provide regional navigation service for Asia-Pacific region after 2012, more new generation BDS satellites have been launched to further expand BDS’s coverage to be global. In this contribution, precise positioning models based on BDS and the corresponding mathematical algorithms are presented in detail. Then, an evaluation on BDS’s real-time dynamic positioning and navigation performance is presented in Precise Point Positioning (PPP), Real-time Kinematic (RTK), Inertial Navigation System (INS) tightly aided PPP and RTK modes by processing a set of land-borne vehicle experiment data. Results indicate that BDS positioning Root Mean Square (RMS) in north, east, and vertical components are 2.0, 2.7, and 7.6?cm in RTK mode and 7.8, 14.7, and 24.8?cm in PPP mode, which are close to GPS positioning accuracy. Meanwhile, with the help of INS, about 38.8%, 67.5%, and 66.5% improvements can be obtained by using PPP/INS tight-integration mode. Such enhancements in RTK/INS tight-integration mode are 14.1%, 34.0%, and 41.9%. Moreover, the accuracy of velocimetry and attitude determination can be improved to be better than 1?cm/s and 0.1°, respectively. Besides, the continuity and reliability of BDS in both PPP and RTK modes can also be ameliorated significantly by INS during satellite signal missing periods.     

Ambiguity resolution in precise point positioning with hourly data for global single receiver

Integer ambiguity resolution (IAR) can improve precise point positioning (PPP) performance significantly. IAR for PPP became a highlight topic in global positioning system (GPS) community in recent years. More and more researchers focus on this issue. Progress has been made in the latest years. In this paper, we aim at investigating and demonstrating the performance of a global zero-differenced (ZD) PPP IAR service for GPS users by providing routine ZD uncalibrated fractional offsets (UFOs) for wide-lane and narrow-lane. Data sets from all IGS stations collected on DOY 1, 100, 200 and 300 of 2010 are used to validate and demonstrate this global service. Static experiment results show that an accuracy better than 1 cm in horizontal and 1–2 cm in vertical could be achieved in ambiguity-fixed PPP solution with only hourly data. Compared with PPP float solution, an average improvement reaches 58.2% in east, 28.3% in north and 23.8% in vertical for all tested stations. Results of kinematic experiments show that the RMS of kinematic PPP solutions can be improved from 21.6, 16.6 and 37.7 mm to 12.2, 13.3 and 34.3 mm for the fixed solutions in the east, north and vertical components, respectively. Both static and kinematic experiments show that wide-lane and narrow-lane UFO products of all satellites can be generated and provided in a routine way accompanying satellite orbit and clock products for the PPP user anywhere around the world, to obtain accurate and reliable ambiguity-fixed PPP solutions.     

地磁暴期间北半球高纬度地区电离层变化特征及对精密定位的影响

基于加拿大地区高纬度电离层观测网的电离层闪烁观测数据,分析了2018年8月26日地磁暴事件引发的北半球高纬度地区电离层总电子含量(TEC)异常变化、TEC变化率指数(ROTI)及电离层相位闪烁的变化特征.结果表明:加拿大地区最大异常值约6 TECU,磁暴引发全球电离层TEC异常峰值高达20 TECU;加拿大地区电离层相...     

观测资料和卫星钟差对PPP参数收敛影响分析

在GPS 单点定位中, 参数解算的收敛时间和收敛稳定性是重要的研究内容之一, 影响收敛时间和收敛稳定性的因素很多, 本文主要就观测资料的不同采样间隔、卫星钟差资料的不同采样间隔、不同的定位精度要求对精密单点定位中参数收敛时间的影响进行了深入的分析探讨, 以中国上海GPS综合应用网中的12个测站资料为例, 分析了采样间隔、定位精度要求与收敛时间的关系, 并对不同采样间隔的收敛时间进行了统计分析, 得出一些初步结论.      

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.     

Influence of different data fusion methods on the accuracy of three-dimensional displacements field

Combining multi-source data can improve the accuracy and the spatial resolution of the three-dimensional (3-D) displacements field. How to effectively integrate multi-source data to obtain high-precision and high spatial resolution 3-D displacements field is worthy of further study. The stochastic model and fusion model of integrating multi-source data affect the accuracy of data fusion. In this paper, based on the least squares method, the effects of different stochastic models and data fusion models on the 3-D displacements field’s accuracy are studied. The optimal method for estimating large-scale 3-D displacements field from integrated InSAR, leveling and GPS measurements is obtained. Then we realize the integrating InSAR, leveling and GPS measurements to obtain the high-precision 3-D displacements velocity field in Tianjin (China) from 2016 to 2018. The results are validated with GPS measurements at 6 independent stations, with the root mean squares (RMS) residuals of the discrepancies being 2.39 mm/yr, 2.54 mm/yr and 2.83 mm/yr in eastern, northern and vertical directions, respectively. By comparing different stochastic models, the 3-D displacements field obtained from multi-source data is optimized by the variance component estimation-least squares method, which is better than weighted least squares (WLS) method. By comparing different data fusion models, the accuracy of the horizontal displacements velocity is better than that of interpolated GPS results. The horizontal displacements component has a great influence on the vertical displacements velocity accuracy in the process of acquiring the 3-D displacements velocity by integrating InSAR, GPS and leveling measurements. This study provides a reference method for integrating multi-source data to obtain 3-D displacements field. This method effectively utilizes the advantages of GPS, InSAR and leveling measurements, and extends the limitations of single technical in describing surface-time scale applications. The 3-D displacements information with a large spatial scale and high spatial resolution provide a reliable data basis for studying the crustal movement and its dynamic mechanism.     

Preliminary analysis of the quality and positioning performance of BDS-3 global interoperable signal B1C&B2a

BeiDou-3 Navigation Satellite System (BDS-3) satellites are equipped with the new generation GNSS signals B1C and B2a, which support the interoperability with GPS and Galileo systems. In this study, the pseudo-range multipath error and carrier phase observation noise of the BDS-3 B1C and B2a signals were evaluated based on zero baseline measurements from the day of year (DOY) 113 to 116, 2020. Further, the precision and performance of the single point positioning (SPP) and precise point positioning (PPP) are assessed at 9 stations. This assessment manifests that the standard deviations (STDs) of the pseudo-range multipath error are about 0.09 ~ 0.22 m, while STDs of the carrier phase observation noise are about 0.075 mm. For the single-frequency SPP, its positioning precision is about 2.03 ~ 4.85 m and 3.29 ~ 10.73 m at the 99.99% confidence level in horizontal and vertical directions, respectively, while the dual-frequency SPP precision is about 1.92 ~ 8.02 m and 4.81 ~ 12.77 m in horizontal and vertical directions, respectively. For the daily static PPP, the convergence time is about 20 ~ 30 min, while the daily positioning precision can reach 1.38 ~ 4.42 cm and -1.31 ~ 4.34 cm in horizontal and vertical directions, respectively. In general, the quality and the SPP and PPP performance of the BDS-3 B1C&B2a signals are comparable to the GPS and Galileo.     

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.     

Contribution analysis of QZSS to single-frequency PPP of GPS/BDS/GLONASS/Galileo

The Quasi-Zenith Satellite System (QZSS) established by the Japan Aerospace Exploration Agency mainly serves the Asia-Pacific region and its surrounding areas. Currently, four in-orbit satellites provide services. Most users of GNSS in the mass market use single-frequency (SF) receivers owing to the low cost. Therefore, it is meaningful to analyze and evaluate the contribution of the QZSS to SF precise point positioning (PPP) of GPS/BDS/GLONASS/Galileo systems with the emergence of GNSS and QZSS. This study compares the performances of three SF PPP models, namely the GRoup and PHase Ionospheric Correction (GRAPHIC) model, GRAPHIC with code observation model, and an ionosphere-constrained model, and evaluated the contribution of the QZSS to the SF PPP of GPS/BDS/GLONASS/Galileo systems. Moreover, the influence of code bias on the SF PPP of the BDS system is also analyzed. A two-week dataset (DOY 013–026, 2019) from 10 stations of the MGEX network is selected for validation, and the results show that: (1) For cut-off elevation angles of 15, 20, and 25°, the convergence times for the static SF PPP of GLONASS + QZSS are reduced by 4.3, 30.8, and 12.7%, respectively, and the positioning accuracy is similar compared with that of the GLONASS system. Compared with the BDS single system, the convergence times for the static SF PPP of BDS + QZSS under 15 and 25° are reduced by 37.6 and 39.2%, the horizontal positioning accuracies are improved by 18.6 and 14.1%, and the vertical components are improved by 13.9 and 21.4%, respectively. At cut-off elevation angles of 15, 20, and 25°, the positioning accuracy and precision of GPS/BDS/GLONASS/Galileo + QZSS is similar to that of GPS/BDS/GLONASS/Galileo. And the convergence times are reduced by 7.4 and 4.3% at cut-off elevation angles of 20 and 25°, respectively. In imitating dynamic PPP, the QZSS significantly improves the positioning accuracy of BDS and GLONASS. However, QZSS has little effect on the GPS-only, Galileo-only and GPS/BDS/GLONASS/Galileo. (2) The code bias of BDS IGSO and MEO cannot be ignored in SF PPP. In static SF PPP, taking the frequency band of B1I whose multipath combination is the largest among the frequency bands as an example, the vertical component has a systematic bias of approximately 0.4–1.0 m. After correcting the code bias, the positioning error in the vertical component is lower than 0.2 m, and the positioning accuracy in the horizontal component are improved accordingly. (3) The SF PPP model with ionosphere constraints has a better convergence speed, while the positioning accuracy of the three models is nearly equal. Therefore the GRAPHIC model can be used to get good positioning accuracy in the absence of external ionosphere products, but its convergence speed is slower.     

Comparison of convergence time and positioning accuracy among BDS,GPS and BDS/GPS precise point positioning with ambiguity resolution

Precise point positioning (PPP) usually takes about 30?min to obtain centimetre-level accuracy, which greatly limits its application. To address the drawbacks of convergence speed and positioning accuracy, we develop a PPP model with integrated GPS and BDS observations. Based on the method, stations with global coverage are selected to estimate the fractional cycle bias (FCB) of GPS and BDS. The short-term and long-term time series of wide-lane (WL) FCB, and the single day change of narrow-lane (NL) FCB are analysed. It is found that the range of GPS and BDS non-GEO (IGSO and MEO) WL FCB is stable at up to a 30-day-time frame. At times frame of up to 60?days, the stability is reduced a lot. Whether for short-term or long-term, the changes in the BDS GEO WL FCB are large. Moreover, BDS FCB sometimes undergoes a sudden jump. Besides, 17 and 10 stations were used respectively to investigate the convergence speed and positioning errors with six strategies: BDS ambiguity-float PPP (Bfloat), GPS ambiguity-float PPP (Gfloat), BDS/GPS ambiguity-float PPP (BGfloat), BDS ambiguity-fixed PPP (Bfix), GPS ambiguity-fixed (Gfix), and BDS/GPS ambiguity-fixed (BGfix). The average convergence speed of the ambiguity-fixed solution is greatly improved compared with the ambiguity-float solution. In terms of the average convergence time, the Bfloat is the longest and the BGfix is the shortest among these six strategies. Whether for ambiguity-float PPP or ambiguity-fixed PPP, the convergence reduction time in three directions for the combined system is the largest compared with the single BDS. The average RMS value of the Bfix in three directions (easting (E), northing (N), and up (U)) are 2.0?cm, 1.5?cm, and 5.9?cm respectively, while those of the Gfix are 0.8?cm, 0.5?cm, and 1.7?cm. Compared with single system, the BDS/GPS combined ambiguity-fixed system (BGfix) has the fastest convergence speed and the highest accuracy, with average RMS as 0.7?cm, 0.5?cm, and 1.9?cm for the E, N, U components, respectively.     

The realization and convergence analysis of combined PPP based on raw observation

In order to speed up Precise Point Positioning (PPP)’s convergence, a combined PPP method with GPS and GLONASS which is based on using raw observations is proposed, and the positioning results and convergence time have been compared with that of single system. The ionospheric delays and receiver’s Differential Code Bias (DCB) corrections are estimated as unknown parameters in this method. The numerical results show that the combined PPP has not caused significant impacts on the final solutions, but it greatly improved Position Dilution of Precision (PDOP) and convergence speed and enhanced the reliability of the solution. Meanwhile, the convergence speed is greatly influenced by the receiver’s DCB, positioning results in horizontal which are better than 10 cm can be realized within 10 min. In addition, the ionosphere and DCB products can be provided with high precision.     

Landslide monitoring by combining of CR-InSAR and GPS techniques

Considering the limitations related to the landslide monitoring by Interferometric Synthetic Aperture Radar (InSAR) technique, the method of integration of Globe Positioning System (GPS) with Corner Reflector Interferometric SAR (CR-InSAR) techniques is proposed in this paper. Firstly, deformation in radar line-of-slight (LOS) direction is optimized by introducing the GPS-measured height and atmospheric delay products into the CR-InSAR model. Then, GPS-measured horizontal deformation and CR-InSAR measured LOS deformation are combined to produce the more accurate vertical deformation. Finally, high precision three-dimensional deformation (N, E, U) is projected to the along-slope direction to monitor the actual movement of landslide. In order to test this method, four X-band stripmap-mode TerraSAR images, eight Trihedral Corner Reflectors (TCR) data and eight GPS observed data are collected to monitor the deformation of three potential landslide fields located at the north of Shaanxi province, China. The detailed analysis demonstrates that the estimated precision of along-slope direction is about two times better for proposed method (±1.1 mm) versus GPS (±2.1 mm) in this case. Meanwhile, our result indicates that almost all of the monitoring points present the trends of sliding down along the slope at the different levels from April 9 2011 to August 30 2011, showing the certain instability. Further investigation of the relationship between the magnitudes of displacement at CR points and the implementation of early control reflects the rationality of our result. Our proposed method could provide of the strong support in the high precision landslide deformation monitoring.     

An Integer Precise Point Positioning technique for sea surface observations using a GPS buoy

GPS data dedicated to sea surface observation are usually processed using differential techniques. Unfortunately, the precision of resulting kinematic positions is baseline-length dependent. So, high precision sea surface observations using differential GPS techniques are limited to coasts, lakes, and rivers. Recent improvements in GPS satellite products (orbits, clocks, and phase biases) make phase ambiguity fixing at the zero difference level achievable and opens up the observation of the sea surface without geographical constraints. This paper recalls the concept of the Integer Precise Point Positioning technique and discusses the precision of GPS buoy positioning. A sequential version of the GINS software has been implemented to achieve single epoch GPS positioning. We used 1 Hz data from a two week GPS campaign conducted in the Kerguelen Islands. A GPS buoy has been moored close to a radar gauge and 90 m away from a permanent GPS station. This infrastructure offers the opportunity to compare both kinematic Integer Precise Point Positioning and classical differential GPS positioning techniques to in situ radar gauge data. We found that Precise Point Positioning results are not significantly biased with respect to radar gauge data and that horizontal time series are consistent with differential processing at the sub-centimetre precision level. Nevertheless, standard deviations of height time series with respect to radar gauge data are typically [4–5] cm. The dominant driver for noise at this level is attributed to errors in tropospheric estimates which propagate into position solutions.     

Determining the precipitable water vapor thresholds under different rainfall strengths in Taiwan

Precipitable Water Vapor (PWV) plays an important role for weather forecasting. It is helpful in evaluating the changes of the weather system via observing the distribution of water vapor. The ability of calculating PWV from Global Positioning System (GPS) signals is useful to understand the special weather phenomenon. In this study, 95 ground-based GPS and rainfall stations in Taiwan were utilized from 2006 to 2012 to analyze the relationship between PWV and rainfall. The PWV data were classified into four classes (no, light, moderate and heavy rainfall), and the vertical gradients of the PWV were obtained and the variations of the PWV were analyzed. The results indicated that as the GPS elevation increased every 100?m, the PWV values decreased by 9.5?mm, 11.0?mm, 12.2?mm and 12.3?mm during the no, light, moderate and heavy rainfall conditions, respectively. After applying correction using the vertical gradients mentioned above, the average PWV thresholds were 41.8?mm, 52.9?mm, 62.5?mm and 64.4?mm under the no, light, moderate and heavy rainfall conditions, respectively. This study offers another type of empirical threshold to assist the rainfall prediction and can be used to distinguish the rainfall features between different areas in Taiwan.     

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.     

The impact of microwave absorber and radome geometries on GNSS measurements of station coordinates and atmospheric water vapour

We have used microwave absorbing material in different geometries around ground-based Global Navigation Satellite System (GNSS) antennas in order to mitigate multipath effects on the estimates of station coordinates and atmospheric water vapour. The influence of a hemispheric radome – of the same type as in the Swedish GPS network SWEPOS – was also investigated. Two GNSS stations at the Onsala Space Observatory were used forming a 12 m baseline. GPS data from October 2008 to November 2009 were analyzed by the GIPSY/OASIS II software using the Precise Point Positioning (PPP) processing strategy for five different elevation cutoff angles from 5° to 25°. We found that the use of the absorbing material decreases the offset in the estimated vertical component of the baseline from ∼27 mm to ∼4 mm when the elevation cutoff angle varies from 5° to 20°. The horizontal components are much less affected. The corresponding offset in the estimates of the atmospheric Integrated Water Vapour (IWV) decreases from ∼1.6 kg/m² to ∼0.3 kg/m². Changes less than 5 mm in the offsets in the vertical component of the baseline are seen for all five elevation cutoff angle solutions when the antenna was covered by a hemispheric radome. Using the radome affects the IWV estimates less than 0.4 kg/m² for all different solutions. IWV comparisons between a Water Vapour Radiometer (WVR) and the GPS data give consistent results.     

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.     

Seismic deformation of the Mw 8.0 Wenchuan earthquake from high-rate GPS observations

High-rate GPS positioning has been recognized as a powerful tool in estimating epoch-wise station displacement which is particularly useful for seismology. In this study, station displacements during the 12 May 2008 M_w 8.0 Wenchuan earthquake are derived from the 1-Hz GPS data collected at a set of stations in China. The impacts of integer ambiguity resolution and station environment-dependent effects are investigated in order to yield more accurate results. The position accuracy of horizontal components of better than 1 cm suggests that GPS can sense the rapid position oscillation of about 2 cm in amplitude. Temporal and spatial analysis is applied to the surface displacement at station XANY and the characteristics of the movements due to Rayleigh and Love waves are detected and discussed. The comparison of GPS-derived displacement with relevant synthetic data computed based on a recently published rapture model shows a reasonable agreement in waveform. The various differences in amplitude need further investigation and also imply that rapture inversion might be improved if GPS-derived displacement is assimilated.     

Real-time clock comparison and monitoring with multi-GNSS precise point positioning: GPS,GLONASS and Galileo

The precise point positioning (PPP) technique is widely used in time and frequency applications. Because of the real-time service (RTS) project of the International GNSS Service, we can use the PPP technique for real-time clock comparison and monitoring. As a participant in the RTS, the Centre National d’Etudes Spatiales (CNES) implements the PPPWIZARD (Precise Point Positioning with Integer and Zero-difference Ambiguity Resolution Demonstrator) project to validate carrier phase ambiguity resolution. Unlike the Integer-PPP (IPPP) of the CNES, fixing ambiguities in the post-processing mode, the PPPWIZARD operates in the real-time mode, which is also called real-time IPPP (RT-IPPP). This paper focuses on applying the RT-IPPP for real-time clock comparison and monitoring. We review the principle of real-time clock comparison and monitoring, and introduce the methodology of the RT-IPPP technique. The observations of GPS, GLONASS and Galileo were processed for the experiments. Five processing modes were provided in the experiment to analyze the benefits of ambiguity resolution and multi-GNSS. In the clock comparison experiment, the average reduction ratios of standard deviations with respect to the G PPP mode range from 9.7% to 35.0%. In the clock monitoring experiment, G PPP mode can detect clock jumps whose magnitudes are larger than 0.9 ns. The RT-IPPP technique with GRE PPP AR (G) mode allows for the detection of any clock jumps larger than 0.6 ns. For frequency monitoring, G PPP mode allows detection of frequency changes larger than 1.1 × 10⁻¹⁴. When the RT-IPPP technique is applied, monitoring with GRE PPP AR (G) mode can detect frequency changes larger than 6.1 × 10⁻¹⁵.     

Real-time zenith tropospheric delays in support of numerical weather prediction applications

The Geodetic Observatory Pecný (GOP) routinely estimates near real-time zenith total delays (ZTD) from GPS permanent stations for assimilation in numerical weather prediction (NWP) models more than 12 years. Besides European regional, global and GPS and GLONASS solutions, we have recently developed real-time estimates aimed at supporting NWP nowcasting or severe weather event monitoring. While all previous solutions are based on data batch processing in a network mode, the real-time solution exploits real-time global orbits and clocks from the International GNSS Service (IGS) and Precise Point Positioning (PPP) processing strategy. New application G-Nut/Tefnut has been developed and real-time ZTDs have been continuously processed in the nine-month demonstration campaign (February–October, 2013) for selected 36 European and global stations. Resulting ZTDs can be characterized by mean standard deviations of 6–10 mm, but still remaining large biases up to 20 mm due to missing precise models in the software. These results fulfilled threshold requirements for the operational NWP nowcasting (i.e. 30 mm in ZTD). Since remaining ZTD biases can be effectively eliminated using the bias-reduction procedure prior to the assimilation, results are approaching the target requirements in terms of relative accuracy (i.e. 6 mm in ZTD). Real-time strategy and software are under the development and we foresee further improvements in reducing biases and in optimizing the accuracy within required timeliness. The real-time products from the International GNSS Service were found accurate and stable for supporting PPP-based tropospheric estimates for the NWP nowcasting.