Deformation monitoring using GNSS-R technology

GNSS reflectometry (GNSS-R) has been widely studied in recent years for various applications, such as soil moisture monitoring, biomass analysis, and sea state monitoring. This paper presents the concept of a novel application of using GNSS-R technology for deformation monitoring. Instead of installing GNSS on the deformation body to sense the movement, GNSS-R deformation monitoring system estimates the deformation from receiving GNSS signal reflected by the deformation body remotely. A prototype of GNSS-R deformation monitoring system has been developed based on GNSS software receiver technology. A 3D geometrical model of GNSS signal reflection has been used to reveal the relationship between the change of carrier phase difference and deformation. After compensating the propagation path delay changes caused by satellite movement, the changes in the remaining carrier phase difference are linked to the deformation. Field tests have been carried using the GNSS-R system developed and the results show sub-centimeter level deformation can be observed with the new technology. Unlike other GNSS deformation monitoring methods, GNSS-R receivers are not installed on the slope which makes this new technology more attractive.     

Noise analysis of continuous GPS coordinate time series for CMONOC

The Crustal Movement Observation Network of China (CMONOC) is one of the major scientific infrastructures, mainly using Global Positioning System (GPS) measurements, to monitor crustal deformation in the Chinese mainland. In this paper, decade-long coordinate time series of 26 continuous GPS sites of CMONOC are analyzed for their noise content using maximum likelihood estimation (MLE). We study the noise properties of continuous GPS time series of CMONOC for the unfiltered, filtered solutions and also the common mode signals in terms of power law plus white noise model. In the spatial filtering, we remove for every time series a common mode error that was estimated from a modified stacking of position residuals from other sites within ∼1000 km of the selected site. We find that the common mode signal in our network has a combination of spatially correlated flicker noise and a common white noise with large spatial extent. We demonstrate that for the unfiltered solutions of CMONOC continuous GPS sites the main colored noise is a flicker process, with a mean spectral index of ∼1. For the filtered solutions, the main colored noise is a general power law process, indicating that a major number of the filtered regional solutions have a combination of noise sources or local effects. The velocity uncertainties from CMONOC continuous GPS coordinate time series may be underestimated by factors of 8–16 if a pure white noise model is assumed. In addition, using a white plus flicker noise model, the median values of velocity errors for the unfiltered solutions are 0.16 (north), 0.17 (east) and 0.58 (vertical) mm/yr, and the median values for the filtered solutions are 0.09 (north), 0.10 (east) and 0.40 (vertical) mm/yr.     

Crustal deformation characteristics of Sichuan-Yunnan region in China on the constraint of multi-periods of GPS velocity fields

In order to obtain deformation parameters in each block of Sichuan-Yunnan Region (SYG) in China by stages and establish a dynamic model about the variation of the strain rate fields and the surface expansion in this area, we taken the Global Positioning System (GPS) sites velocity in the region as constrained condition and taken advantage of the block strain calculation model based on spherical surface. We also analyzed the deformation of the active blocks in the whole SYG before and after the Wenchuan earthquake, and analyzed the deformation of active blocks near the epicenter of the Wenchuan earthquake in detail. The results show that, (1) Under the effects of the carving from India plate and the crimping from the potential energy of Tibetan Plateau for a long time, there is a certain periodicity in crustal deformation in SYG. And the period change and the earthquake occurrence have a good agreement. (2) The differences in GPS velocity fields relative Eurasian reference frame shows that the Wenchuan earthquake and the Ya'an earthquake mainly affect the crustal movement in the central and southern part of SYG, and the average velocity difference is about 4–8?mm/a for the Wenchuan earthquake and 2–4?mm/a for the Ya'an earthquake. (3) For the Wenchuan earthquake, the average strain changed from 10 to 20 nanostrian/a before earthquake to 40–50 nanostrian/a after the earthquake, but before and after the Ya’an earthquake, the strain value increased from about 15 nanostrian/a to about 30 nanostrian/a. (4) The Wenchuan earthquake has changed the strain parameter of each active block more or less. Especially, the Longmen block and Chengdu block near the epicenter. The research provides fundamental material for the study of the dynamic mechanism of the push extrusion from the north-east of the India plate and the crimp from Qinghai Tibet Plateau, and it also provides support for the study of crustal stress variation and earthquake prediction in Sichuan Yunnan region.     

High temporal resolution global PWV dataset of 2005–2016 by using a neural network approach to determine the mean temperature of the atmosphere

Atmospheric water vapour plays an important role in phenomena related to the global hydrologic cycle and climate change. However, the rapid temporal–spatial variation in global tropospheric water vapour has not been well investigated due to a lack of long-term, high-temporal-resolution precipitable water vapour (PWV). Accordingly, this study generates an hourly PWV dataset for 272 ground-based International Global Navigation Satellite System (GNSS) Service (IGS) stations over the period of 2005–2016 using the zenith troposphere delay (ZTD) derived from global-scale GNSS observation. The root mean square (RMS) of the hourly ZTD obtained from the IGS tropospheric product is approximately 4 mm. A fifth-generation reanalysis dataset of the European Centre for Medium-range Weather Forecasting (ECMWF ERA5) is used to obtain hourly surface temperature (T) and pressure (P), which are first validated with GNSS synoptic station data and radiosonde data, respectively. Then, T and P are used to calculate the water vapour-weighted atmospheric mean temperature (Tm) and zenith hydrostatic delay (ZHD), respectively. T and P at the GNSS stations are obtained via an interpolation in the horizontal and vertical directions using the grid-based ERA5 reanalysis dataset. Here, Tm is calculated using a neural network model, whereas ZHD is obtained using an empirical Saastamoinen model. The RMS values of T and P at the collocated 693 radiosonde stations are 1.6 K and 3.1 hPa, respectively. Therefore, the theoretical error of PWV caused by the errors in ZTD, T and P is on the order of approximately 2.1 mm. A practical comparison experiment is performed using 97 collocated radiosonde stations and 23 GNSS stations equipped with meteorological sensors. The RMS and bias of the hourly PWV dataset are 2.87/?0.16 and 2.45/0.55 mm, respectively, when compared with radiosonde and GNSS stations equipped with meteorological sensors. Additionally, preliminary analysis of the hourly PWV dataset during the EI Niño event of 2014–2016 further indicates the capability of monitoring the daily changes in atmospheric water vapour. This finding is interesting and significant for further climate research.     

An improved cyclic multi model-eXtreme gradient boosting (CMM-XGBoost) forecasting algorithm on the GNSS vertical time series

The study of GNSS vertical coordinate time series forecasting is helpful for monitoring the crustal plate movement, dam or bridge deformation monitoring, and global or regional coordinate system maintenance. The eXtreme Gradient Boosting (XGBoost) algorithm is a machine learning algorithm that can evaluate features, and it has a great potential and stability for long-span time series forecasting. This study proposes a multi-model combined forecasting method based on the XGBoost algorithm. The method constitutes a new time series as features through the fitting and forecasting results of the forecasting model. The XGBoost model is then used for forecasting. In addition, this method can obtain higher precision forecasting results through circulation. To verify the performance of the forecasting method, 1095 epochs of data in the Up coordinate of 16 GNSS stations are selected for the forecasting test. Compared with the CNN-LSTM model, the experimental results of our forecasting method show that the mean absolute error (MAE) values are reduced by 30.23 %～52.50 % and the root mean square error (RMSE) values are reduced by 31.92 %～54.33 %. The forecasting results have higher accuracy and are highly correlated to the original time series, which can better forecast the vertical movement of the GNSS stations. Therefore, the forecasting method can be applied to the up component of the GNSS coordinate time series.     

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.     

Evaluation and improvement of coastal GNSS reflectometry sea level variations from existing GNSS stations in Taiwan

Global sea level rise due to an increasingly warmer climate has begun to induce hazards, adversely affecting the lives and properties of people residing in low-lying coastal regions and islands. Therefore, it is important to monitor and understand variations in coastal sea level covering offshore regions. Signal-to-noise ratio (SNR) data of Global Navigation Satellite System (GNSS) have been successfully used to robustly derive sea level heights (SLHs). In Taiwan, there are a number of continuously operating GNSS stations, not originally installed for sea level monitoring. They were established in harbors or near coastal regions for monitoring land motion. This study utilizes existing SNR data from three GNSS stations (Kaohsiung, Suao, and TaiCOAST) in Taiwan to compute SLHs with two methods, namely, Lomb–Scargle Periodogram (LSP)-only, and LSP aided with tidal harmonic analysis developed in this study. The results of both methods are compared with co-located or nearby tide gauge records. Due to the poor quality of SNR data, the worst accuracy of SLHs derived from traditional LSP-only method exceeds 1?m at the TaiCOAST station. With our procedure, the standard deviations (STDs) of difference between GNSS-derived SLHs and tide gauge records in Kaohsiung and Suao stations decreased to 10?cm and the results show excellent agreement with tide gauge derived relative sea level records, with STD of differences of 7?cm and correlation coefficient of 0.96. In addition, the absolute GNSS-R sea level trend in Kaohsiung during 2006–2011 agrees well with that derived from satellite altimetry. We conclude that the coastal GNSS stations in Taiwan have the potential of monitoring absolute coastal sea level change accurately when our proposed methodology is used.     

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.     

SHAKING: Adjusted spherical harmonics adding KrigING method for near real-time ionospheric modeling with multi-GNSS observations

Precise positioning based on Global Navigation Satellite System (GNSS) technique requires high accuracy ionospheric total electron content (TEC) correction models to account for the ionospheric path delay errors. We present an adjusted Spherical Harmonics Adding KrigING method (SHAKING) approach for regional ionospheric vertical TEC (VTEC) modeling in real time. In the proposed SHAKING method, the VTEC information over the sparse observation data area is extrapolated by the Adjusted Spherical Harmonic (ASH) function, and the boundary distortion in regional VTEC modeling is corrected by the stochastic VTEC estimated using Kriging interpolation. Using real-time GPS, GLONASS and BDS-2/3 data streams of the Crust Movement Observation Network of China (CMONOC), the SHAKING-based regional ionospheric VTEC maps are re-constructed over China and its boundary regions. Compared to GNSS VTECs derived from the independent stations, the quality of SHAKING solution improves by 13–31% and 6–33% with respect to the ASH-only solution during high and low geomagnetic periods, respectively. Compared to the inverse distance weighting (IDW) generated result, significant quality improved of SHAKING-based VTEC maps is also observed, especially over the edge areas with an improvement of 60–80%. Overall, the proposed SHAKING method exhibits notable advantage over the existing regional VTEC modeling techniques, which can be used for regional TEC modeling and associated high-precision positioning applications.     

Exploration and analysis of the factors influencing GNSS PWV for nowcasting applications

Precipitable water vapor (PWV) can be assimilated into a numerical weather model (NWM) to improve the prediction accuracy of numerical weather prediction. In this study, taking GNSS data for the Beijing Fangshan station (BJFS) as an example, based on the method of Pearson correlation coefficient combined with quantitative analysis, GNSS datasets are used to study the relationships between GNSS-derived PWV (GNSS PWV__Met) and its influencing factors, including the internal influencing factors zenith troposphere delay (ZTD), zenith hydrostatic delay (ZHD), zenith wet delay (ZWD), and surface temperature (T_s), and the external influencing factor haze (mainly PM_2.5). Firstly, based on the strong correlation between PWV__Met and ZTD hourly sequences from the International GNSS Service Network’s BJFS station for DOYS 182–212, 2015, the results of experiment prove that the reliability of GNSS ZTD is used to forecast PWV_Met in short-term forecasting. Secondly, based on hourly data of BJFS in 2016, the correlation between PWV__Met and ZTD, ZWD, ZHD, pressure (P) and T_s is analyzed, and then, with the rate of ZTD variation as the main factor, ZTD variation as auxiliary factor, the prediction success rate is 88.24% from hourly data of precipitation event for DOYs 183–213 in Beijing. The experiment indicates that ZTD can help forecast short-term precipitation. Thirdly, based on data from three hazy periods with relatively stable weather conditions, no heavy rainfall, and relatively continuous data in the past three years, the correlation between GNSS PWV_Met/ZTD and PM_2.5 hourly series is analyzed. The results of the experiments suggests that GNSS ZTD should be considered to assist in haze monitoring. So in the absence of radiosonde stations and meteorological elements, ZTDs on retrieval of GNSS stations have more application value in short-term forecast.     

Total column water vapor estimation over land using radiometer data from SAC-D/Aquarius

The aim of this study is retrieving atmospheric total column water vapor (CWV) over land surfaces using a microwave radiometer (MWR) onboard the Scientific Argentine Satellite (SAC-D/Aquarius). To research this goal, a statistical algorithm is used for the purpose of filtering the study region according to the climate type.A log-linear relationship between the brightness temperatures of the MWR and CWV obtained from Global Navigation Satellite System (GNSS) measurements was used. In this statistical algorithm, the retrieved CWV is derived from the Argentinian radiometer’s brightness temperature which works at 23.8?GHz and 36.5?GHz, and taking into account CWVs observed from GNSS stations belonging to a region sharing the same climate type. We support this idea, having found a systematic effect when applying the algorithm; it was generated for one region using the previously mentioned criteria, however, it should be applied to additional regions, especially those with other climate types.The region we analyzed is in the Southeastern United States of America, where the climate type is Cfa (Köppen - Geiger classification); this climate type includes moist subtropical mid-latitude climates, with hot, muggy summers and frequent thunderstorms. However, MWR only contains measurements taken from over ocean surfaces; therefore the determination of water vapor over land is an important contribution to extend the use of the SAC-D/Aquarius radiometer measurements beyond the ocean surface. The CWVs computed by our algorithm are compared against radiosonde CWV observations and show a bias of about ?0.6?mm, a root mean square (rms) of about 6?mm and a correlation of 0.89.     

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.     

Study of seasonal and long-term vertical deformation in Nepal based on GPS and GRACE observations

Lithospheric deformation signal can be detected by combining data from continuous global positioning system (CGPS) and satellite observations from the Gravity Recovery and Climate Experiment (GRACE). In this paper, we use 2.5- to 19-year-long time series from 35 CGPS stations to estimate vertical deformation rates in Nepal, which is located in the southern side of the Himalaya. GPS results were compared with GRACE observations. Principal component analysis was conducted to decompose the time series into three-dimensional principal components (PCs) and spatial eigenvectors. The top three high-order PCs were calculated to correct common mode errors. Both GPS and GRACE observations showed significant seasonal variations. The observed seasonal GPS vertical variations are in good agreement with those from the GRACE-derived results, particularly for changes in surface pressure, non-tidal oceanic mass loading, and hydrologic loading. The GPS-observed rates of vertical deformation obtained for the region suggest both tectonic impact and mass decrease. The rates of vertical crustal deformation were estimated by removing the GRACE-derived hydrological vertical rates from the GPS measurements. Most of the sites located in the southern part of the Main Himalayan Thrust subsided, whereas the northern part mostly showed an uplift. These results may contribute to the understanding of secular vertical crustal deformation in Nepal.     

Refinement of global ionospheric coefficients for GNSS applications: Methodology and results

We developed the methodology for the optimal estimation of global ionospheric coefficients of the current Global Navigation Satellite Systems (GNSSs), including the eight- and ten-parameter Klobuchar-like as well as NeQuick models. The ionospheric coefficients of those correction models are calculated from two sets of globally distributed tracking stations of the International GNSS Services (IGS). Performance of the re-estimated Klobuchar-like and NeQuick coefficients are validated during 2002–2014 over the continental and oceanic areas, respectively. Over the continental areas, GPS TECs derived from 40 ground GPS receivers are selected as reference. The eight-, ten-parameter Klobuchar-like and NeQuick models can mitigate the ionospheric delay by 65.8, 67.3 and 75.0%, respectively. Over the global oceans, the independent TECs derived from Jason-1&2 altimeters are used as reference. The re-estimated ionospheric correction models can mitigate 56.1–66.7% of the delay errors. Compared to the original GPS Ionospheric Correction Algorithm (ICA), performance of those eight-, ten-parameter Klobuchar-like and NeQuick models has improved 3.4, 5.9 and 13.4% during the whole test period, respectively. The methodology developed here takes the advantage of high-quality ionospheric TECs derived from the global network of GNSS receivers. The re-estimated ionospheric coefficients can be used as precise ionospheric products to monitor and assess GNSS broadcast ionospheric parameters and to improve the performance of various single-frequency GNSS applications.     

TS-InSAR analysis for monitoring ground deformation in Lanzhou New District,the loess Plateau of China,from 2017 to 2019

Large-scale land creation projects involving the cutting of mountains to infill gullies for construction have been carried out in Lanzhou New District (LZND). However, there is an urgent need for comprehensive and detailed research on the spatiotemporal evolution of ground deformation in LZND. Based on Sentinel-1A SAR data, combined with the urban geological background, the ground deformation in LZND from 2017 to 2019 was analysed. Two independent, multi-temporal techniques, persistent scatterers interferometry (PS-InSAR) and the small baseline subset (SBAS-InSAR), were used to calculate the deformation time series, and the results were cross-verified. The time series-monitoring results of the SBAS and PS calculations exhibited strong consistency in LZND and verified the high reliability of the experimental results. The results showed the whole surface of the LZND from March 2017 to October 2019 maintained stability, and the deformation rate was primarily in the range of ?10 to 10 mm/year. However, ground deformation in the Xicha area was evident. The maximum annual deformation rates monitored by SBAS-InSAR and PS-InSAR were ?52.48 mm/year and ?56.35 mm/year, respectively. The most typical deformation areas include the built-up area and the land creation area. The surface subsidence area was concentrated in the filling area. The ground deformation range of LZND kept expanding and accelerating from 2017 to 2019. Land creation, urban construction, geology and precipitation were the primary factors contributing to local severe ground deformation. The results of this study provide reference for the regional urban planning in LZND.     

Evaluation and analysis on positioning performance of BDS/QZSS satellite navigation systems in Asian-Pacific region

By using the observation data and products of precise obit and clock offset from Multi-GNSS Experiment (MGEX) of the International GNSS Service (IGS) and GNSS Research Centre, Curtin University in this paper, the positioning performance of BDS/QZSS satellite navigation system has been analyzed and evaluated in aspects of the quantity of visible satellites, DOP value, multipath effect, signal-to-noise ratio, static PPP and kinematic PPP. The analysis results show that compared to BDS single system when the cutoff angle are 30°and 40°, the DOP value of BDS/QZSS combined system has decreased above 20%, and the quantity of visible satellites increased about 16–30% respectively, because of the improved spatial geometric configuration. The magnitude of satellite multipath effect of BDS system shows the trend of MEO?>?IGSO?>?GEO, which is consistent with that of QZSS satellite system, as the constellation structure of the two systems is similar. The variation tendencies of signal-to-noise ratio with respect to elevation angle of the two systems are almost the same at all frequencies, showing that at the same elevation angle the signal-to-noise ratio of MEO satellites is higher than that of IGSO satellites, as the higher obit is the lower transmitting power is obtained. For having a specially designed obit, the variation of signal-to-noise ratio of BDS system is more stable. However, the magnitude of signal-to-noise ratio of QZSS system appears the trend of frequency 3?>?frequency 2?>?frequency 1. The static PPP performance of the BDS/QZSS combination system has been improved more significantly than the BDS single system in E, N and U directions. When the cutoff angle are at 7°, 15° and 30°, the PPP accuracy is increased about 25–34% in U direction, 10–13% and 23–34% in E and N directions respectively. When the elevation angle is large (40°), compared to BDS single system at lower elevation angles (7° and 15°) the PPP accuracy of the BDS/QZSS combination system is improved above 30% in U direction. In kinematic PPP performance, compared to BDS single system, the accuracy, availability and reliability of the BDS/QZSS combination system has been improved too, especially at large elevation angles (30° and 40°), the kinematic PPP accuracy in E and U directions has been improved about 10–50%, and above 50% in U direction. It can be concluded that the combination with QZSS system can improve the positioning accuracy, reliability and stability of BDS system. In the future, with the improvement of the satellite construction of Japan’s QZSS system and the global networking of China’s BDS satellites, the QZSS satellites will contribute greatly to improve the positioning accuracy, reliability, availability and stability of GNSS systems in areas such as cities, mountains, densely-packed buildings and severely covered areas in Asian-Pacific region.     

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.     

Analysis of precipitable water vapour characteristics from GNSS measurements during the snow season in Liaoning Province,China

Global Navigation Satellite System (GNSS) remote sensing precipitable water vapour (PWV) data from November 2015 to March 2019 were combined with snowfall observation data and used to analyse PWV characteristics in Liaoning Province during the snow season (from November to March the following year) and their relationship with snowfall. The potential of using GNSS for PWV measurements was demonstrated using sounding data with a correlation coefficient higher than 0.9 and a mean bias error lower than 0.5 mm. According to the GNSS PWV data gathered at 30-min intervals from 68 GNSS stations in Liaoning during the snow season, the monthly PWV average was highest in November and lowest in January. Negative correlations were found between PWV and altitude. Most of the water vapour was concentrated in the low layer of the atmosphere, and the contribution of this vapour to the PWV was higher during the snow season than in summer. A total of 43 snow cases were identified using the snowfall records from 53 GNSS stations, and the characteristics of PWV during these snowfalls were analysed. An increase in PWV was observed before snowfall events. Moreover, the influence of synoptic systems and air mass origins on PWV was analysed based on National Centers for Environmental Prediction (NCEP) reanalysis data and the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. The results show that the water vapour condition was better when the synoptic systems or air masses came from areas south of Liaoning.     

Present-day deformation in São Jorge Island,Azores, from episodic GPS measurements (2001–2011)

A network of geodetic markers for high-precision Global Positioning System (GPS) surveying was established in 2001 in the island of São Jorge, at the Azores triple junction, where the North America–Eurasia–Nubia plates meet. The aim was to monitor tectonic, volcanic, and landslide-induced surface deformation in an active tectono-volcanic region. The network consisted of 17 stations, and was occupied in 2001, 2004, and 2010. GPS observations from the three episodic campaigns were analyzed using the latest geophysical models and geodetic procedures to generate a velocity field for São Jorge Island.     

海南地区电离层不规则体纬向漂移速度的观测和研究

根据中国海南富克(19.3°N,109.1°E)三点GPS观测系统2007年3月至11月的观测数据,利用互相关方法分析了三站闪烁信号的时间延迟,得出了不规则体纬向漂移的基本特征.在中国海南地区,闪烁主要发生在春秋季节,夜间不规则体的纬向漂移速度以东向为主,大小在50～150 m/s之间;平均东向漂移速度随时间呈下降趋势.另外,在闪烁刚发生时,不规则体纬向速度起伏较大,这可能与不规则体的随机起伏以及等离子体泡产生时垂直速度较大有关.中国海南地区不规则体纬向漂移速度的这些基本特征与低纬其他地区的测量结果较为一致.