Search

The frequent occurrence of extreme climate events has a significant impact on people's lives. Heavy rainfall can lead to an increase of regional Terrestrial Water Storage (TWS), which will cause land subsidence due to the influence of hydrological load. At present, regional TWS is mostly obtained from Gravity Recovery and Climate Experiment (GRACE) data, but the method has limitations for small areas. This paper used water level and flow data as hydrological signals to study the land subsidence caused by heavy rainfall in the Chaohu Lake area of East China (June 2016–August 2016). Pearson's correlation coefficient was used to study the interconnection between water resource changes and Global Navigation Satellites System (GNSS) vertical displacement. Meanwhile, to address the reliability of the research results, combined with the Coefficient of determination method, the research findings were validated by using different institutional models. The results showed that: (1) During heavy rainfall, the vertical displacement caused by atmospheric load was larger than non-tidal oceanic load, and the influence trends of the two were opposite. (2) The rapidly increasing hydrologic load in the Chaohu Lake area resulted in greater subsidence displacement at the closer CORS station (CHCH station) than the more distant CORS station (LALA station). The Pearson correlation coefficients between the vertical displacement and water level were as high as −0.80 and −0.64, respectively. The phenomenon confirmed the elastic deformation principle of disc load. (3) Although there was a systematic bias between the different environmental load deformation models, the deformation trends were generally consistent with the GNSS monitoring results. The average Coefficients of determination between the different models and the GNSS results were 0.63 and 0.77, respectively. The results demonstrated the effectiveness of GNSS in monitoring short-term hydrological load. This study reveals the spatial-temporal evolution of land deformation during heavy rainfall around Chaohu Lake, which is of reference significance for water resource management and infrastructure maintenance in this area. [ABSTRACT FROM AUTHOR]

The uneven distribution of global navigation satellite system (GNSS) continuous stations in the Yellow River Basin, combined with the sparse distribution of GNSS continuous stations in some regions and the weak far-field load signals, poses challenges in using GNSS vertical displacement data to invert terrestrial water storage changes (TWSCs). To achieve the inversion of water reserves in the Yellow River Basin using unevenly distributed GNSS continuous station data, in this study, we employed the Tikhonov regularization method to invert the terrestrial water storage (TWS) in the Yellow River Basin using vertical displacement data from network engineering and the Crustal Movement Observation Network of China (CMONOC) GNSS continuous stations from 2011 to 2022. In addition, we applied an inverse distance weighting smoothing factor, which was designed to account for the GNSS station distribution density, to smooth the inversion results. Consequently, a gridded product of the TWS in the Yellow River Basin with a spatial resolution of 0.5 degrees on a daily scale was obtained. To validate the effectiveness of the proposed method, a correlation analysis was conducted between the inversion results and the daily TWS from the Global Land Data Assimilation System (GLDAS), yielding a correlation coefficient of 0.68, indicating a strong correlation, which verifies the effectiveness of the method proposed in this paper. Based on the inversion results, we analyzed the spatial–temporal distribution trends and patterns in the Yellow River Basin and found that the average TWS decreased at a rate of 0.027 mm/d from 2011 to 2017, and then increased at a rate of 0.010 mm/d from 2017 to 2022. The TWS decreased from the lower-middle to lower reaches, while it increased from the upper-middle to upper reaches. Furthermore, an attribution analysis of the terrestrial water storage changes in the Yellow River Basin was conducted, and the correlation coefficients between the monthly average water storage changes inverted from the results and the monthly average precipitation, evapotranspiration, and surface temperature (AvgSurfT) from the GLDAS were 0.63, −0.65, and −0.69, respectively. This indicates that precipitation, evapotranspiration, and surface temperature were significant factors affecting the TWSCs in the Yellow River Basin. [ABSTRACT FROM AUTHOR]

The continuous downward movement exhibited by the Lashagou landslide group in recent years poses a significant threat to the safety of both vehicles and pedestrians traversing the highway G310. By integrating geomorphological interpretation using multi-temporal optical images, interferometric synthetic aperture radar (InSAR) measurements, and continuous global navigation satellite system (GNSS) observations, this paper traced the formation period of the Lashagou landslide group, and explored its kinematic behavior under external drivers such as rainfall and snowmelt. The results indicate that the formation period can be specifically categorized into three periods: before, during, and after the construction of highway G310. The construction of highway G310 is the direct cause and prerequisite for the formation of the Lashagou landslide group, whereas summer precipitation and spring snowmelt are the external driving factors contributing to its continuous downward movement. Additionally, both the long-term seasonal downslope movement and transient acceleration events are strongly controlled by rainfall, and there is a time lag of approximately 1–2 days between the transient acceleration and heavy rainfall events. This study highlights the benefits of leveraging multi-source remote sensing data to investigate slow-moving landslides, which is advantageous for the implementation of effective control and engineering intervention to mitigate potential landslide disasters. [ABSTRACT FROM AUTHOR]

Underground coal-mining-induced ground subsidence deformation is a common geological disaster impacting buildings, transportation and water supplies. Models predicting ground subsidence dynamically with high precision are important for the prevention of damage derived from ground subsidence. In this paper, the Hook function is utilized to develop a model describing the velocity of ground subsidence due to underground coal mining. Based on the subsidence velocity model, a dynamic subsidence model is established by taking an integral of the velocity model. Coefficients of the model, which depend on maximum subsidence, maximum subsidence velocity and the time corresponding to the maximum subsidence velocity, are related to the geological and mining conditions of the coal seam being investigated. A Levenberg–Marquardt-algorithm-based method is also proposed to calculate the optimal model coefficients based on subsidence velocity observations. Four continuously operating Global Navigation Satellite System (GNSS) stations were constructed above a typical longwall coal mining working face in the Jining mining area, China. These GNSS stations collected subsidence observations over two years, which were used to validate the developed prediction model. The results show that the root-mean-square (RMS) of the model-predicted ground subsidence error is 56.1 mm, and the maximum relative error is 2.5% for all four GNSS stations, when the ground subsidence is less than 6000 mm. [ABSTRACT FROM AUTHOR]

As a major grain-producing region in China, the North China Plain (NCP) faces serious challenges such as water shortage and land subsidence. In late 2014, the Central Route of the South-to-North Water Diversion Project (SNWD-C) began to provide NCP with water resources. However, the effectiveness of this supply in mitigating land subsidence remains a pivotal and yet unassessed aspect. In this paper, we utilized various geodetic datasets, including the Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow On (GRACE-FO), Global Navigation Satellite System (GNSS) and leveling data, to conduct a spatial-temporal analysis of the equivalent water height (EWH) and vertical ground movement in the NCP. The results reveal a noteworthy decline in EWH from 2011 to 2015, followed by a slight increase with minor fluctuations from 2015 to 2020, demonstrating a strong correlation with the water resources supplied by the SNWD-C. The GRACE-derived surface deformation rate induced by hydrological loading is estimated to be <1 mm/yr. In comparison, GNSS-derived vertical ground movements exhibit considerable regional differences during the 2011–2020 period. Substantial surface subsidence is evident in the central and eastern NCP, contrasting with a gradual uplift in the front plain of the Taihang Mountains. Three-stage leveling results indicate that the rate of subsidence in the central and eastern plains is gradually increasing with the depression area expanding from 1960 to 2010. Based on these geodetic results, it can be inferred that the SNWD-C's operation since 2014 has effectively mitigated the reduction in terrestrial water storage in the NCP. However, land subsidence in the NCP persists, as the subsidence rate does not turn around in sync with the change in EWH following the operation of SNWD-C. Consequently, it's necessary to maintain and enforce existing policies, including controlling groundwater exploitation and water resources supply (e.g., SNWD-C) to curtail the exacerbation of land subsidence in the NCP. Additionally, continuous monitoring of land subsidence by GRACE, GNSS, leveling and other geodetic techniques is crucial to enable timely policy adjustments based on monitoring results. [ABSTRACT FROM AUTHOR]

This article takes the measured Total Electron Content (TEC) from the GPS points of the China Regional Crust Observation Network as the starting point to establish a regional ionospheric empirical model. The model's performance is enhanced by considering solar flux and geomagnetic activity data. The refinement function model of the ionospheric TEC diurnal variation component, seasonal variation component, and geomagnetic component is studied. Using the nonlinear least squares method to fit undetermined coefficients, MEFM-ITCR (Multi-parameter Empirical Fusion Model–Ionospheric TEC China Regional Model) is proposed to forecast the regional ionosphere TEC in China. The results show that the standard deviation of MEFM-ITCR residuals is 3.74TECU, and MEFM-ITCR fits the modeling dataset well. Analyses of geographic location variation, seasonal variation, and geomagnetic disturbance were carried out for MEFM-ITCR performance. The results indicate that in the Chinese region, MEFM-ITCR outperforms IRI2020 and NeQuick2 models in terms of forecast accuracy, linear correlation, and model precision for TEC measured using GPS points under different latitudes and longitudes, different seasons, and different geomagnetic disturbances. The empirical TEC model built for the Chinese region in this paper provides a new ionospheric delay correction method for GNSS single frequency users and is of great significance for establishing other new and improving existing ionospheric empirical models. [ABSTRACT FROM AUTHOR]

Air-quality monitoring and forecasting are crucial for atmosphere pollution control and management. We propose an innovative data-driven framework for air quality index (AQI) prediction by integrating GNSS radio occultation (GNSS-RO) observation and weather modeling. Empowered by the state-of-the-art machine learning approach, our method can effectively predict regional AQI with a comparable accuracy much more quickly than the traditional numerical modeling and simulation approach. In a real case study using a representative region of China, our data-driven approach achieves a 2000 times speedup; meanwhile, the prediction error measured by rRMSE is only 2.4%. We investigate further the effects of different models, hyperparameters, and meteorological factors on the performance of our AQI prediction framework, and reveal that wind field and atmospheric boundary-layer height are important influencing factors of AQI. This paper showcases a direct application of GNSS-RO observation in assisting in forecasting regional AQI. From a machine learning point of view, it provides a new way to leverage the unique merits of GNSS atmospheric remote sensing technology with the help of the more traditional weather forecasting modeling approach. [ABSTRACT FROM AUTHOR]

The global navigation satellite system (GNSS) position time series provides essential data for geodynamic and geophysical studies. Interpolation of the GNSS position time series is necessary because missing data will produce inaccurate conclusions made from the studies. The spatio-temporal correlations between GNSS reference stations cannot be considered when using traditional interpolation methods. This paper examines the use of machine learning models to reflect the spatio-temporal correlation among GNSS reference stations. To form the machine learning problem, the time series to be interpolated are treated as output values, and the time series from the remaining GNSS reference stations are used as input data. Specifically, three machine learning algorithms (i.e., the gradient boosting decision tree (GBDT), eXtreme gradient boosting (XGBoost), and random forest (RF)) are utilized to perform interpolation with the time series data from five GNSS reference stations in North China. The results of the interpolation of discrete points indicate that the three machine learning models achieve similar interpolation precision in the Up component, which is 45% better than the traditional cubic spline interpolation precision. The results of the interpolation of continuous missing data indicate that seasonal oscillations caused by thermal expansion effects in summer significantly affect the interpolation precision. Meanwhile, we improved the interpolation precision of the three models by adding data from five stations which have high correlation with the initial five GNSS reference stations. The interpolated time series for the North, East, and Up (NEU) are examined by principal component analysis (PCA), and the results show that the GBDT and RF models perform interpolation better than the XGBoost model. [ABSTRACT FROM AUTHOR]

In 2020, the BeiDou-3 global navigation satellite system (BDS-3) was officially completed and put into service. Currently, network real-time kinematic (RTK) technology is considered the main means through which to improve the positioning accuracy of the BeiDou navigation satellite system (BDS). This paper proposes a long-range undifferenced network RTK (URTK) algorithm, based on multi-frequency observation data of the BDS. First, the multi-frequency phase integer ambiguity resolution (AR) model considering atmospheric error parameters is designed, and the multi-frequency phase integer ambiguity of the long-range BDS reference station is determined. Then, the undifferenced integer ambiguity of each reference station is obtained, using linear variation based on the accurately determined phase integer ambiguity between reference stations, and the undifferenced observation error of each reference station is calculated. Considering the weakening spatial correlation of the observation errors between long-range stations, undifferenced classification error corrections of a reference station network are separated, according to different error characteristics. Finally, the inverse distance weighting method is employed to calculate the classification undifferenced error correction of the rover station. The rover station corrects the observation error through applying the undifferenced error correction to achieve high-precision positioning. The measured data of a long-range continuous operation reference station (CORS) network are selected for an experiment. The results show that the proposed algorithm can quickly and accurately realize the resolution of the BDS integer ambiguity of a reference station network and establish an undifferenced area error correction model in order to achieve accurate classification of undifferenced error correction values for a rover station. In China, the BDS-3 is superior to the global positioning system (GPS) in terms of the satellite number, position dilution of precision (PDOP) value, AR success rate, stability, and convergence time. The results show that the AR success rate, stability, and convergence time increase with the operational frequency, and the BDS-3 can achieve centimeter-level positioning of single-system rover stations without relying on the GPS. [ABSTRACT FROM AUTHOR]

In this paper, a method for joint sea surface rainfall intensity (RI), wind speed, and wave height retrieval based on spaceborne global navigation satellite system reflectometry (GNSS-R) data is proposed, which especially considers the effects between these two parameters. A method of rainfall detection (RD) according to different wind speed ranges is also proposed by mitigating the impact of swell and wind speed. The results, with data collected over the oceans near Southeast Asia, show that the RD method has a detection accuracy of up to 81.74%. The RI retrieval accuracy can reach about 2 mm/h by simultaneously correcting the effects of wind speed and swell. The accuracy of wind speed retrieval is improved by about 5% after removing rainfall interference through RD in advance. After considering the influence of wind speed and eliminating rainfall interference, the retrieval accuracy of significant wave height (SWH) is improved by about 18%. Finally, the deep convolutional neural network (DCNN) model is built to estimate the SWH of the swell. The results show that the retrieval accuracy of the swell height is better than 0.20 m after excluding rainfall interference. The proposed joint retrieval method provides an important reference for the future acquisition of multiple high-precision marine geophysical parameters by spaceborne GNSS-R technology. [ABSTRACT FROM AUTHOR]

The classification that distinguishes whether machines are driving on roads or working in fields based on their global navigation satellite system (GNSS) trajectories is essential for effective management of cross-regional agricultural machinery services in China. In this paper, a novel field–road classification method utilizing multiple deep neural networks (MultiDNN) is proposed to enhance the accuracy of field and road point classification. The MultiDNN model incorporates a bi-directional long short-term memory network (BiLSTM), a topology adaptive graph convolution network (TAG), and a self-attention network (ATT) to effectively extract spatio-temporal features for field–road classification. The BiLSTM is used to capture temporal relationships along the time axis of a trajectory, providing global contextual information for each point. Then, the TAG network is used to obtain the spatio-temporal relationships between adjacent points in a trajectory, offering local contextual information for each point. Finally, the ATT network assigns varying weights to features to emphasize important characteristics. The performance of the MultiDNN model was evaluated using a wheat harvesting trajectory dataset, and the results showed that it achieved a high degree of accuracy, up to 89.75%, outperforming the best baseline method (GCN) by 2.79%. [ABSTRACT FROM AUTHOR]

As the Global Navigation Satellite System (GNSS) is widely used in all walks of life, the signal structure of satellite navigation is open, and the vulnerability to spoofing attacks is also becoming increasingly prominent, which will seriously affect the credibility of navigation, positioning, and timing (PNT) services. Satellite navigation signal authentication technology is an emerging technical means of improving civil signal anti-spoofing on the satellite navigation system side, and it is also an important development direction and research focus of the GNSS. China plans to carry out the design and development of the next-generation Beidou navigation satellite system (BDS), and one of its core goals is to provide more secure and credible PNT services. This paper first expounds on the principles and technical architecture of satellite navigation signal authentication, then clarifies the development history of satellite navigation signal authentication, and finally proposes the BDS authentication service system architecture. It will provide technical support for the construction and development of the follow-up Beidou authentication service. [ABSTRACT FROM AUTHOR]

Precise Point Positioning (PPP) is an official service of the BeiDou Global Navigation Satellite System (BDS-3) through the PPP-B2b signal. In this paper, we mainly focus on the long-term performance evaluation of BDS-3 PPP-B2b products and their application in time service. Since the PPP-B2b product is only available in and around China area, the arcs of PPP-B2b products are about several hours. We propose to evaluate the time datum stability by using all available satellites. Then, 557 day PPP-B2b products are collected for this experiment. The results show that there are large jumps in the GPS satellite clock time datum series. However, the BDS-3 satellite clock datum stability is almost at the same level with current Space State Representation (SSR) corrections from the International Global navigation satellite system Service (IGS). The difference between PPP-B2b GPS and BDS-3 satellite clock time datum will be absorbed into the Inter System Bias (ISB) parameter. Thus, it should be specially noted that the ISB parameter cannot be estimated as constant values if users use PPP-B2b products. In addition, the accuracy of the BDS-3 satellite clock is significantly better than that of the GPS for both the Root Mean Square Error (RMSE) and standard deviation (STD). The average Signal in Space Range Errors (SISREs) is 0.22 ns and 0.13 ns for GPS and BDS-3, respectively. The one-way timing experiment shows BDS-3 timing stability is 2.9 × 10−14@104 s. In addition, 10 baselines from 13 km to 4494 km are formed for time synchronization evaluation by using PPP-B2b products. The average RMSEs of time synchronization is from 0.46 ns to 1.58 ns and from 0.66 ns to 1.19 ns for GPS and BDS-3, respectively. As for STD, the average values are from 0.27 ns to 0.74 ns and from 0.27 ns to 0.47 ns for GPS and BDS-3, respectively. Overall, the results show that the time datum stability, accuracy, and service performance of BDS-3 PPP-B2b products has been stable over the past two years. [ABSTRACT FROM AUTHOR]

Landslides along the Three Gorges Reservoir in China pose a threat to coastal residents and waterway safety. To reduce false positive misjudgments caused by a sudden local change in the landslide deformation curve, in this paper, we propose an effective method for predicting the deformation trend of reservoir bank landslides. We take reservoir bank landslides in the Wanzhou District of the Three Gorges Reservoir area as the research object. The Time-Series Interferometric Synthetic Aperture Radar (InSAR) method and 62 Sentinel-1A images from 2018 to 2022 were selected for landslide deformation monitoring, and the Hurst index was calculated to characterize the deformation trend. Furthermore, we propose a method for predicting the deformation trend based on the statistical distribution of deformation rates and the physical significance of the Hurst index. After the field survey and Global Positioning System (GPS) verification, the Time-Series InSAR results are shown to be reliable. We take the Sifangbei landslide as a representative case to analyze the validation results. It is found that the determined Sifangbei landslide deformation trend is consistent with the conclusions for the region. In addition, the deformation trend of a reservoir bank slope has obvious spatial and temporal differences. Changes in the reservoir water level and concentrated rainfall play roles similar to those of catalysts. The proposed method, involving the combination of Time-Series InSAR and the Hurst index, can effectively monitor deformation and predict the stability trend of reservoir bank landslides. The presented research results provide new ideas and solutions for landslide prevention and risk mitigation in reservoir areas. [ABSTRACT FROM AUTHOR]

China developed its Digital Television (DTV) standard in 2006, known as Digital Television Terrestrial Multimedia Broadcasting (DTMB), which employs time-domain synchronous orthogonal frequency division multiplexing (TDS-OFDM) as the modulation method. In contrast to Global Navigation Satellite Systems (GNSSs), DTV signals have higher transmitting power, wider coverage, larger bandwidth, and fixed transmitter location. This paper explores the positioning performance of DTMB signals, and the potential to improve GNSS positioning accuracy in urban environments. Specifically, a solution is proposed, and a software-defined radio receiver is developed for wireless localization. Without changing the current signal structure, the pseudorandom noise (PN) sequences in the signal are used for signal acquisition and carrier phase ranging. The carrier phase of the first arrived path is extracted by the least squares matching pursuit method. Both static and dynamic field tests were conducted to verify the proposed ranging and positioning method. Centimeter-level ranging accuracy was achieved in the static scenario, while meter-level ranging accuracy was achieved in the dynamic scenario. As one possible application, the proposed ranging method was combined with GPS pseudorange measurements to achieve higher accuracy position results in an urban pedestrian scenario, especially when there is only a limited number of visible satellites. [ABSTRACT FROM AUTHOR]

On 21 May 2021 (UT), Yangbi Ms6.4 and Maduo Ms7.4 earthquakes occurred in mainland China. This paper analyzed the ionospheric perturbations possibly related to the earthquake, based on global positioning system (GPS) total electron content (TEC) and global ionosphere map (GIM) TEC data. We identified GPS TEC anomalies by the sliding quartile, based on statistical analysis. After eliminating the days with high solar activity levels and strong geomagnetic disturbances, the time series analysis of GPS TEC data showed that there were significant TEC anomalies from 5 to 10 May. TEC anomalies were mainly positive anomalies. We obtained the spatial and temporal distributions of TEC anomalies using natural neighbor interpolation (NNI). The results showed that the TEC anomalies were distributed in the seismogenic zone and surrounded the epicenters of the Maduo and Yangbi earthquakes, indicating that they may be related to the earthquakes. From the GIM TEC difference map, we found the TEC enhancement in the seismogenic zone and its magnetic conjugate area of the Maduo and Yangbi earthquakes at 10:00–12:00 (UT) on the 5 and 6 May. We discussed our results according to the lithosphere-atmosphere-ionosphere coupling mechanism. Finally, based our results, we suggested that the Yangbi and Maduo earthquakes may affect the ionosphere through seismogenic electric field and thermal anomalies generated during the process of lithosphere-atmosphere-ionosphere coupling. [ABSTRACT FROM AUTHOR]

With the rapid development of interferometric synthetic aperture radar (InSAR) measurement technology, its measurement accuracy requirements are increasing. Atmospheric delay errors must be corrected, especially in the case of crustal deformation monitoring, the 20% variation of tropospheric water vapor among InSAR pairs generally produces range from 10 cm to 14 cm deformation errors. Such errors can be of the same magnitude as the annual changes in crustal deformation, or even greater, masking crustal deformation information and seriously affecting the results of crustal deformation monitoring. Therefore, in order to obtain a more accurate InSAR atmospheric delay correction model, this paper calculated and integrated atmospheric delays that were estimated by different sources, including the 37 pressure levels of the fifth generation of the European Centre for Medium-Range Weather Forecasts (ECMWF)) numerical weather prediction model, ECMWF Reanalysis v5 (ERA5), and Global Navigation Satellite System (GNSS) measurement data from the crustal movement observation network of China, based on the variance component estimation (VCE) weighting method. The results showed that the integrated model, based on the VCE method, is better than the generic atmospheric correction online service (GACOS) model for InSAR measuring of crustal deformation. The precision in monitoring crustal deformations was improved by approximately 5 mm, the correlation coefficient of atmospheric delay errors and crustal deformations improved from 0.287 to 0.347, and accuracy improved by approximately 25%. However, the improvement in accuracy was limited because of system error decoherence that was induced by atmospheric noise caused by abundant vegetation or snow cover. Therefore, in order to achieve more accurate results, we recommend the adoption of the multi-source integrated atmospheric delay correction model, based on the VCE method, for InSAR high-precision measuring of crustal deformation and seismic activities. [ABSTRACT FROM AUTHOR]

Atmospheric weighted mean temperature (Tm) is a key parameter used by the Global Navigation Satellite System (GNSS) for calculating precipitable water vapor (PWV). Some empirical Tm models using meteorological or non-meteorological parameters have been proposed to calculate PWV, but their accuracy and reliability cannot be guaranteed in some regions. To validate and determine the optimal Tm model for PWV retrieval in China, this paper analyzes and evaluates some typical Tm models, namely, the Linear, Global Pressure and Temperature 3 (GPT3), the Tm model for China (CTm), the Global Weighted Mean Temperature-H (GTm-H) and the Global Tropospheric (GTrop) models. The Tm values of these models are first obtained at corresponding radiosonde (RS) stations in China over the period of 2011 to 2020. The corresponding Tm values of 87 RS stations in China are also calculated using the layered meteorological data and regarded as the reference. Comparison results show that the accuracy of these five Tm models in China has an obvious geographical distribution and decreases along with increasing altitude and latitude, respectively. The average root mean square (RMS) and Bias for the Linear, GPT3, CTm, GTm-H and GTrop models are 4.2/3.7/3.4/3.6/3.3 K and 0.7/−1.0/0.7/−0.1/0.3 K, respectively. Among these models, Linear and GPT3 models have lower accuracy in high-altitude regions, whereas CTm, GTm-H and GTrop models show better accuracy and stability throughout the whole China. These models generally have higher accuracy in regions with low latitude and lower accuracy in regions with middle and high latitudes. In addition, Linear and GPT3 models have poor accuracy in general, whereas GTm-H and CTm models are obviously less accurate and stable than GTrop model in regions with high latitude. These models show different accuracies across the four geographical regions of China, with GTrop model demonstrating the relatively better accuracy and stability. Therefore, the GTrop model is recommended to obtain Tm for calculating PWV in China. [ABSTRACT FROM AUTHOR]

Continuous Global Positioning Systems (GPS) coordinate time series with a high spatiotemporal resolution, and provide a great opportunity to study their noise models and common mode errors (CMEs), thus making it possible to detect and analyse spatiotemporal characteristics of tectonic and non-tectonic signals in time series, and further to estimate a reliable and accurate velocity field of crustal movement in a region by removing CMEs and using the optimal noise model. In this paper, we used GPS coordinate time series from the Crustal Movement Observation Network of China (CMONOC) with an approximate decadal period from 2010 to 2020, to construct optimal noise models by fitting them with several noise combinations according to the Akaike information criterion (AIC). We further adopted independent component analysis (ICA) to extract CMEs and analysed their spatiotemporal characteristics, and then evaluated their effects on noise models and velocity uncertainties, and finally estimated a decennial velocity field of crustal movement with a higher signal-to-noise ratio (SNR) by applying the CME filtering and considering the optimal noise model in Mainland China. Our results show that optimal noise models are dominated by white noise (WN) plus flicker noise (FN) for both east and north components, and WN plus power law noise (PN) with spectral index close to −1 for up component, respectively. ICA filtering can remove the highly spatially correlated CMEs and decrease the mean RMSEs of the residual time series by about 40–60%, providing a more accurate velocity field with a higher SNR in Mainland China, accordingly. [ABSTRACT FROM AUTHOR]

Seasonal variations in the vertical Global Positioning System (GPS) time series are mainly caused by environmental loading, e.g., hydrological loading (HYDL), atmospheric loading (ATML), and nontidal oceanic loading (NTOL), which can be synthesized based on models developed by various institutions. A comprehensive comparison among these models is essential to extract reliable vertical deformation data, especially on a regional scale. In this study, we selected 4 HYDL, 5 ATML, 2 NTOL, and their 40 combined products to investigate their effects on seasonal variations in vertical GPS time series at 27 GPS stations in Yunnan, southwest China. These products were provided by the German Research Center for Geosciences (GFZ), School and Observatory of Earth Sciences (EOST), and International Mass Loading Service (IMLS). Furthermore, we used the Cross Wavelet Transform (XWT) method to analyze the relative phase relationship between the GPS and the environmental loading time series. Our result showed that the largest average Root-Mean-Square (RMS) reduction value was 1.32 mm after removing the deformation associated with 4 HYDL from the vertical GPS time series, whereas the RMS reductions after 5 ATML and 2 NTOL model corrections were negative at most stations in Yunnan. The average RMS reduction value of the optimal combination of environmental loading products was 1.24 mm, which was worse than the HYDL (IMLS_GEOSFPIT)-only correction, indicating that HYDL was the main factor responding for seasonal variations at most stations in Yunnan. The XWT result showed that HYDL also explained the annual variations reasonably. Our finding implies that HYDL (IMLS_GEOSFPIT) contributes the most to the environmental loading in Yunnan, and that the ATML and NTOL models used in this paper cannot be effective to correct seasonal variations. [ABSTRACT FROM AUTHOR]

Atmospheric water vapor is one of the important parameters for weather and climate studies. Generally, atmospheric water vapor can be monitored by some techniques, such as the Global Navigation Satellite System (GNSS), radiosonde (RS), remote sensing and numerical weather forecast (NWF). However, the comprehensive retrieval and application of precipitable water vapor (PWV) using multi techniques has been hardly performed before, which becomes the focus of this study. A comprehensive PWV retrieval and application platform (CPRAP) is first established by combing the ground-based (GNSS), space-based (Fengyun-3A, Sentinel-3A) and reanalysis-based (the fifth-generation reanalysis dataset of the European Centre for Medium-Range Weather Forecasting, ERA5) techniques. Additionally, its applications are then extended to drought and rainfall monitoring using the CPRAP-derived PWV. The statistical result shows that PWV derived from ground-based GNSS has high accuracy in China, with the root mean square (RMS), Bias and mean absolute error (MAE) of 2.15, 0.05 and 1.65 mm, respectively, when the RS-derived PWV is regarded as the reference. In addition, the accuracy of PWV derived from the space-based (FY-3A and Sentinel-3A) techniques technique is also validated and the RMS, Bias and MAE of a Medium Resolution Spectral Imager (MERSI) onboard Fengyun-3A (FY-3A) and an Ocean and Land Color Instrument (OLCI) onboard Sentinel-3A are 4.46/0.56/3.61 mm and 2.95/0.01/1.37 mm, respectively. Then, the performance of ERA5-derived PWV is evaluated based on GNSS-derived and RS-derived PWV. The result also shows good accuracy of ERA5-provided PWV with the averaged RMS, Bias and MAE of 1.86/0.11/1.48 mm and 0.90/−0.05/1.51 mm, respectively. Finally, the PWV data derived from the established CPRAP are further used for drought and rainfall monitoring. The applied results reveal that the calculated the standardized precipitation evapotranspiration index (SPEI) using the CPRAP-derived PWV can monitor the drought and the correlation coefficient ranges from 0.83 to 0.9 when compared with the SPEI. Furthermore, in this paper correlation analysis between PWV derived from the CPRAP and rainfall, and its potential for rainfall monitoring was also validated. Such results verify the significance of the established CPRAP for weather and climate studies. [ABSTRACT FROM AUTHOR]

Unmanned pavement construction is of great significance in China, and the primary issue to be solved is how to identify the boundaries of the Pavement Construction Area (PCA). In this paper, we present a simple yet effective method, named the Bidirectional Sliding Window (BSW) method, for PCA boundary recognition. We first collected the latitude and longitude coordinates of the four vertices of straight quadrilaterals using the Global Positioning System—Real Time Kinematic (GPS-RTK) measurement principle for precise single-point positioning, analyzed single-point positioning accuracy, and determined the measurement error distribution models. Next, we took points at equal intervals along one straight line segment and two curved line segments with curvature radii of 70 m to 300 m, for simulation experiments. BSW was adopted to recognize the Possible Irrelevant Points (PIP) and Relevant Points (RP), which were used to identify PCA boundaries. Experiments show that when the proposed BSW algorithm is used and the single-point positioning accuracy is at the centimeter level, PCA boundary recognition for straight polygons reaches single-point positioning accuracy, and that for curved polygons reaches decimeter-level accuracy. [ABSTRACT FROM AUTHOR]

The solid Earth deforms elastically in response to variations of surface atmosphere, hydrology, and ice/glacier mass loads. Continuous geodetic observations by Global Positioning System (CGPS) stations and Gravity Recovery and Climate Experiment (GRACE) record such deformations to estimate seasonal and secular mass changes. In this paper, we present the seasonal variation of the surface mass changes and the crustal vertical deformation in the South China Block (SCB) identified by GPS and GRACE observations with records spanning from 1999 to 2016. We used 33 CGPS stations to construct a time series of coordinate changes, which are decomposed by empirical orthogonal functions (EOFs) in SCB. The average weighted root-mean-square (WRMS) reduction is 38% when we subtract GRACE-modeled vertical displacements from GPS time series. The first common mode shows clear seasonal changes, indicating seasonal surface mass re-distribution in and around the South China Block. The correlation between GRACE and GPS time series is analyzed which provides a reference for further improvement of the seasonal variation of CGPS time series. The results of the GRACE observations inversion are the surface deformations caused by the surface mass change load at a rate of about -0.4 to -0.8 mm/year, which is used to improve the long-term trend of non-tectonic loads of the GPS vertical velocity field to further explain the crustal tectonic movement in the SCB and surroundings. [ABSTRACT FROM AUTHOR]

The Haiyang-2D (HY2D) satellite is the fourth ocean dynamics environment monitoring satellite launched by China. The satellite operates on a re-entry frozen orbit, which necessitates orbital maneuvers to maintain its designated path once the satellite's sub-satellite point deviates beyond a certain threshold. However, the execution of orbit maneuvers presents a significant challenge to the field of Precise Orbit Determination (POD). The thesis selects the on-board GPS data of HY2D satellite in December 2023 and five maneuvering days of that year. Employing a multifaceted approach that includes the assessment of observational data quality, orbit overlap, external orbit validation, and SLR (Satellite Laser Ranging) verification, the research delves into precise orbit determination during both maneuver and non-maneuver periods. The results indicate that: (1) The average number of satellites tracked by the receiver is 6.4; (2) During the non-maneuver periods, the average RMS (Root Mean Square) value of the radial difference in the 6-h overlapping arc segment is 0.66 cm, and the three-dimensional position difference is about 1.16 cm; (3) When compared with the precision science orbits (PSO) provided by CNES (Centre National d'Études Spatiales), the average values of the RMS values of the differences in the radial (R), transverse (T), and normal (N) directions during the non-maneuver and maneuver periods are respectively 1.32 cm, 2.31 cm, 1.92 cm and 3.04 cm, 8.78 cm, 2.72 cm. (4) The SLR verification of the orbit revealed a residual RMS of 2.24 cm. This suggests that by incorporating the modeling of maneuver forces during the maneuver periods, the impact of orbital maneuvers on orbit determination can be mitigated. [ABSTRACT FROM AUTHOR]

The BeiDou Global Navigation Satellite System (BDS-3) provides real-time precise point positioning (PPP) service via B2b signals, offering real-time decimeter-level positioning for users in China and surrounding areas. However, common interruptions and outliers in PPP-B2b services arise due to factors such as the Geostationary Orbit (GEO) satellite "south wall effect", Issue of Data (IOD) matching errors, and PPP-B2b signal broadcast priorities, posing challenges to continuous high-precision positioning. This study meticulously examines the completeness, continuity, and jumps in PPP-B2b orbit and clock correction using extensive observational data. Based on this analysis, a two-step method for detecting outliers in PPP-B2b orbit and clock corrections is devised, leveraging epoch differences and median absolute deviation. Subsequently, distinct prediction methods are developed for BDS-3 and GPS orbit and clock corrections. Results from simulated and real-time dynamic positioning experiments indicate that predicted corrections can maintain the same accuracy as normal correction values for up to 10 min and sustain decimeter-level positioning accuracy within 30 min. The adoption of predicted correction values significantly enhances the duration of sustaining real-time PPP during signal interruptions. [ABSTRACT FROM AUTHOR]

Deformation along a subway rail network is related to the safe operation of the subway and the stability of construction facilities on the surface, making long-term deformation monitoring imperative. Long-term monitoring of surface deformation along the subway network and statistical analysis of the overall deformation situation are lacking in China. Therefore, targeting 35 Chinese cities whose subway mileage exceeds 50 km, we extracted their surface deformation along subway networks between 2018 and 2022, using spaceborne synthetic aperture radar (SAR) interferometry (InSAR) technology and Sentinel-1 satellite data. We verified the results with the continuous global navigation satellite system (GNSS) stations' data and found that the root mean square error (RMSE) of the InSAR results was 3.75 mm/year. Statistical analysis showed that ground subsidence along the subways was more prominent in Beijing, Tianjin, and other areas in the North China Plain, namely Kunming (which is dominated by karst landforms), as well as Shanghai, Guangzhou, Qingdao, and other coastal cities. In addition, an analysis revealed that the severity of surface subsidence correlated positively with a city's gross domestic product (GDP) with a Pearson correlation of 0.787, since the higher the GDP, the more frequent the construction and maintenance of subway, and the more commuters there are, which in turn exacerbates the disturbance to the surface. Additionally, the type of land cover also affects the ground deformation. Our findings provide a reference for constructing, operating, and maintaining the urban subway systems in China. [ABSTRACT FROM AUTHOR]

In recent years, the GPS wave buoy has been developed for in situ wave monitoring based on satellite GPS signals. Many research works have been completed on the GPS-based wave measurement technology and great progress has been achieved. The basic principle of the GPS wave buoy is to calculate the movement velocity of the buoy using the Doppler frequency shift of satellite GPS signals, and then to calculate the wave parameters from the movement velocity according to ocean wave theory. The shortage of the GPS wave buoy is the occasional occurrence of some unusual values in the movement velocity. This is mainly due to the fact that the GPS antenna is occasionally covered by sea water and cannot normally receive high-quality satellite GPS signals. The traditional solution is to remove these unusual movement velocity values from the records, which requires furthering extend the acquisition time of satellite GPS signals to ensure there is a large enough quantity of effective movement velocity values. Based on the traditional GPS wave measurement technology, this paper presents the algorithmic flow and proposes two improvement measures. On the one hand, the neural network algorithm is used to correct the unusual movement velocity data so that extending the acquisition time of satellite GPS signals is not necessary and battery power is saved. On the other hand, the Gaussian low-pass filter is used to correct the raw directional wave spectrum, which can further eliminate the influence of noise spectrum energy and improve the measurement accuracy. The on-site sea test of the SBF7-1A GPS wave buoy, developed by the National Ocean Technology Center in China, and the gravity-acceleration-type DWR-MKIII Waverider buoy are highlighted in this article. The wave data acquired by the two buoys are analyzed and processed. It can be seen from the processed results that the ocean wave parameters from the two kinds of wave buoys, such as wave height, wave period, wave direction, wave frequency spectrum, and directional wave spectrum, are in good consistency, indicating that the SBF7-1A GPS wave buoy is comparable to the traditional gravity-acceleration-type wave buoy in terms of its accuracy. Therefore, the feasibility and validity of the two improvement measures proposed in this paper are confirmed. [ABSTRACT FROM AUTHOR]

The most commonly used real-time augmentation services in China are the International GNSS Service's (IGS) real-time service (RTS), PPP-B2b service, and Double-Frequency Multi-Constellation (DFMC) service of the BeiDou Satellite-Based Augmentation System (BDSBAS) service. However, research on the performance evaluation, comparison, and application scope of these three products is still incomplete. This article introduces methods for obtaining real-time augmentation information and real-time orbit and clock offset recovery. Based on real-time orbit and clock offset accuracy, positioning accuracy, and positioning availability, this article systematically evaluates the performance and analyzes the application fields of Centre National d'Études Spatiales (CNES), PPP-B2b, and BDSBAS augmentation information. The results of the evaluation revealed that the radial accuracy of the CNES and PPP-B2b real-time orbit product is consistent, and the Root Mean Square (RMS) is better than 5 cm. The CNES real-time orbit product can achieve centimeter-level accuracy in both along-track and cross-track components, surpassing PPP-B2b's decimeter-level accuracy. Both services demonstrate consistent accuracy in the real-time clock offset, with PPP-B2b showing similar standard deviations (STDs) of 0.16 ns for different satellites. However, for CNES, the STD of the real-time clock offset varies, with values of 0.10 ns, 0.19 ns, and 0.60 ns, respectively, for GPS, BDS-3 Medium Earth Orbit (MEO), and BDS-3 Inclined Geosynchronous Satellite Orbit (IGSO) satellites. Centimeter-level accuracy is achieved after convergence and positioning availability exceeds 99% for CNES and PPP-B2b services. Therefore, the difference between the two services in application areas depends on the acquisition of augmentation information. However, BDSBAS, which concentrates on code observations, demonstrates inferior performance in real-time orbit, clock offset, positioning accuracy, and positioning availability when compared to the other two services. Its primary application is in the aviation and maritime domains, where there is a greater need for service integrity, continuity, and reliability. [ABSTRACT FROM AUTHOR]

Precise regional ionospheric total electron content (TEC) models play a crucial role in correcting ionospheric delays for single-frequency receivers and studying variations in the Earth's space environment. A particle swarm optimization neural network (PSO-NN)-based model for ionospheric TEC over China has been developed using a long-term (2008–2021) ground-based global positioning system (GPS), COSMIC, and Fengyun data under geomagnetic quiet conditions. In this study, a spatial gridding approach is utilized to propose an improved version of the PSO-NN model, named the PSO-NN-GRID. The root-mean-square error (RMSE) and mean absolute error (MAE) of the TECs estimated from the PSO-NN-GRID model on the test data set are 3.614 and 2.257 TECU, respectively, which are 7.5% and 5.5% smaller than those of the PSO-NN model. The improvements of the PSO-NN-GRID model over the PSO-NN model during the equinox, summer, and winter of 2015 are 0.4–22.1%, 0.1–12.8%, and 0.2–26.2%, respectively. Similarly, in 2019, the corresponding improvements are 0.5–13.6%, 0–10.1%, and 0–16.1%, respectively. The performance of the PSO-NN-GRID model is also verified under different solar activity conditions. The results reveal that the RMSEs for the TECs estimated by the PSO-NN-GRID model, with F10.7 values ranging within [0, 80), [80, 100), [100, 130), [130, 160), [160, 190), [190, 220), and [220, +), are, respectively, 1.0%, 2.8%, 4.7%, 5.5%, 10.1%, 9.1%, and 28.4% smaller than those calculated by the PSO-NN model. [ABSTRACT FROM AUTHOR]

Weighted mean temperature (Tm) is an important parameter in the water vapor inversion of global navigation satellite systems (GNSS). High-precision Tm values can effectively improve the accuracy of GNSS precipitable water vapor. In this study, a new regional grid Tm empirical model called the RGTm model over China and the surrounding areas was proposed by combining ERA5 reanalysis data, radiosonde data, and TanDEM-X 90m products. In the process of model establishment, we considered the setting of the reference height in the height correction formula and the bias correction for the Tm lapse rate. Tm values derived from ERA5 and radiosonde data in 2019 were used as references to validate the performance of the RGTm model. At the same time, the GPT3, GGNTm, and uncorrected seasonal model were used for comparison. Results show that compared with the other three models, the accuracy of the RGTm model's Tm was improved by approximately 12.21% (15.32%), 1.17% (3.09%), and 2.31% (5.05%), respectively, when ERA5 (radiosonde) Tm data were used as references. In addition, the introduction of radiosonde data prevented the accuracy of the Tm empirical model from being entirely dependent on the accuracy of the reanalysis data. [ABSTRACT FROM AUTHOR]

GPS has become a very effective tool to remotely sense precipitable water vapor (PWV) information, which is important for weather forecasting and nowcasting. The number of geodetic GNSS stations set up in China has substantially increased over the last few decades. However, GPS PWV derivation requires surface pressure to calculate the precise zenith hydrostatic delay and weighted mean temperature to map the zenith wet delay to precipitable water vapor. GPS stations without collocated meteorological sensors can retrieve water vapor using standard atmosphere parameters, which lead to a decrease in accuracy. In this paper, a method of interpolating NWP reanalysis data to site locations for generating corresponding meteorological elements is explored over China. The NCEP FNL dataset provided by the NCEP (National Centers for Environmental Prediction) and over 600 observed stations from different sources was selected to assess the quality of the results. A one-year experiment was performed in our study. The types of stations selected include meteorological sites, GPS stations, radio sounding stations, and a sun photometer station. Compared with real surface measurements, the accuracy of the interpolated surface pressure and air temperature both meet the requirements of GPS PWV derivation in most areas; however, the interpolated surface air temperature exhibits lower precision than the interpolated surface pressure. At more than 96% of selected stations, PWV differences caused by the differences between the interpolation results and real measurements were less than 1.0 mm. Our study also indicates that relief amplitude exerts great influence on the accuracy of the interpolation approach. Unsatisfactory interpolation results always occurred in areas of strong relief. GPS PWV data generated from interpolated meteorological parameters are consistent with other PWV products (radio soundings, the NWP reanalysis dataset, and sun photometer PWV data). The differences between them were approximately 1~3 mm at most at our selected stations, and GPS data processing is the main factor influencing the agreement of the GPS PWV results with those of other methods. [ABSTRACT FROM AUTHOR]

Using Global Navigation Satellite Systems (GNSS) observation data for developing a high-precision ionospheric Total Electron Content (TEC) model is one of the essential subjects in ionospheric physics research and the application of satellite navigation correction. In this study, we integrate the Empirical Orthogonal Function (EOF) method with the TEC data provided by the Center for Orbit Determination in Europe (CODE), and observed by the dense GNSS receivers operated by the Crustal Movement Observation Network of China (CMONOC) to construct a regional ionospheric TEC model over China. The EOF analysis of CODE TEC in China from 1998 to 2010 shows that the first-order EOF component accounts for 90.3813% of the total variation of the ionospheric TEC in China. Meanwhile, the average value of CODE TEC is consistent with the spatial and temporal distribution characteristics of the first-order EOF base function, which mainly reflects the latitude and diurnal variations of TEC in China. The first-order coefficient after EOF decomposition shows an obvious 11-year period and semi-annual variations. The maximum amplitude of semi-annual variation mainly appears in March and October, which is closely associated with the variation in geographical longitude, the semi-annual change of the low-latitude electric field, and the ionospheric fountain effect. The second-order coefficient has an evident annual variation, the minimum amplitude mainly occurs in March, August, and September, and the amplitude values in the high solar activity years are more significant than those in the low solar activity years. The third-order coefficient mainly shows the characteristics of annual variation, and the fourth-order coefficient shows the noticeable semi-annual and annual variations. The third and fourth-order coefficients are both modulated by the solar activity index F10.7. The ionospheric TEC model in China, driven by CMONOC real-time GNSS observation data, can better reflect the latitude, local time and seasonal variation characteristics of ionospheric TEC over China. In particular, it can clearly show the spring and autumn asymmetry of ionospheric TEC in the low latitudes. The root mean square error of the absolute error between the model and the actual observation is mainly distributed around 2.45 TECU (1 TECU = 1016 electrons/m2). The values of the TEC model constructed in this study are closer to the actual observed values than those of the CODE TEC in China. [ABSTRACT FROM AUTHOR]

Positioning of spacecraft (e.g., geostationary orbit (GEO), high elliptical orbit (HEO), and lunar trajectory) is crucial for mission completion. Instead of using ground control systems, global navigation satellite system (GNSS) can be an effective approach to provide positioning, navigation and timing service for spacecraft. In 2020, China finished the construction of the third generation of BeiDou navigation satellite system (BDS-3); this global coverage system will contribute better sidelobe signal visibility for spacecraft. Meanwhile, with more than 100 GNSS satellites, multi-GNSS navigation performance on the spacecraft is worth studying. In this paper, instead of using signal-in-space ranging errors, we simulate pseudorange observations with measurement noises varying with received signal powers. Navigation performances of BDS-3 and its combinations with other systems were conducted. Results showed that, owing to GEO and inclined geosynchronous orbit (IGSO) satellites, all three types (GEO, HEO, and lunar trajectory) of spacecraft received more signals from BDS-3 than from other navigation systems. Single point positioning (SPP) accuracy of the GEO and HEO spacecraft was 17.7 and 23.1 m, respectively, with BDS-3 data alone. Including the other three systems, i.e., GPS, Galileo, and GLONASS, improved the SPP accuracy by 36.2% and 19.9% for GEO and HEO, respectively. Navigation performance of the lunar probe was significantly improved when receiver sensitivity increased from 20 dB-Hz to 15 dB-Hz. Only dual- (BDS-3/GPS) or multi-GNSS (BDS-3, GPS, Galileo, GLONASS) could provide continuous navigation solutions with a receiver threshold of 15 dB-Hz. [ABSTRACT FROM AUTHOR]

The factors affecting tourist behavior are complex and diverse, but research on its effect from a purely spatial perspective is still limited. The aim of this paper is to explore the dichotomous interaction between tourists and islands: the behavioral patterns of tourists in island destinations and the mechanisms by which island spaces constrain tourist behavior. This study uses fine-grained global positioning system (GPS) tracking data actively authorized and released by tourists. We empirically studied tourist behavior from a spatial constraint perspective and discovered the following: island space has a strong influence and constraining effect on tourist behavior; different spatial attributes have different constraining effects on tourist behavior; and people with different identities interact with different attributes of space to produce different spatial properties, resulting in a 'harmony-contradiction' model of spatial interaction. These findings are of great value in expanding the perspective of spatial constraints and exploring the interactions between people and land. They are also of great practical significance in promoting spatial planning synergy, facilitating the construction of high-quality island-based tourism destinations, and building a harmonious relationship between people and land. [ABSTRACT FROM AUTHOR]

This study investigated the nonlinear changes in the vertical motion of 411 GNSS reference stations situated in mainland China and assessed the influence of the environmental load on their vertical displacement. The researchers evaluated the effect of environmental load by calculating the change in annual cycle amplitude before and after its removal, focusing on its impact across regions with distinct foundation types. The results demonstrate that removing the environmental load led to a considerable reduction of approximately 50.25% in the annual cycle amplitude of vertical motion for GNSS reference stations in mainland China. This reduction in amplitude improved the positioning accuracy of the stations, with the highest WRMS reduction being 2.72 mm and an average reduction of 1.03 mm. The most significant impact was observed in the southwestern, northern, and northwestern regions, where the amplitude experienced a notable decrease. Conversely, the southeastern region exhibited a corresponding increase in amplitude. This article innovatively explored the effects of environmental loads on diverse foundation types. When categorizing GNSS reference stations based on their foundation type, namely, bedrock, 18 m soil layer, and 4–8 m soil layer stations, this study found that removing the environmental load resulted in reductions in annual cycle amplitudes of 49.37%, 59.61%, and 46.48%, respectively. These findings indicate that 18 m soil layer stations were more susceptible to environmental load-induced vertical motion. In conclusion, the impact of the environmental load was crucial when analyzing the vertical motion of GNSS reference stations in mainland China, as it was essential for establishing a high-precision coordinate reference framework and studying the tectonic structure of the region. [ABSTRACT FROM AUTHOR]

The disturbances of the ionosphere caused by thunderstorms or lightning events in the troposphere have an impact on global navigation satellite system (GNSS) signals. Gravity waves (GWs) triggered by thunderstorms are one of the main factors that drive short-period Travelling Ionospheric Disturbances (TIDs). At mid-latitudes, ionospheric GWs can be detected by GNSS signals. However, at low latitudes, the multi-variability of the ionosphere leads to difficulties in identifying GWs induced by thunderstorms through GNSS data. Though disturbances of the ionosphere during low-latitude thunderstorms have been investigated, the explicit GW observation by GNSS and its propagation pattern are still unclear. In this paper, GWs with periods from 6 to 20 min are extracted from band-pass filtered GNSS carrier phase observations without cycle-slips, and 0.2–0.8 Total Electron Content Unit (TECU) magnitude perturbations are observed when the trajectories of ionospheric pierce points fall into the perturbed region. The propagation speed of 102.6–141.3 m/s and the direction of the propagation indicate that the GWs are propagating upward from a certain thunderstorm at lower atmosphere. The composite results of disturbance magnitude, period, and propagation velocity indicate that GWs initiated by thunderstorms and propagated from the troposphere to the ionosphere are observed by GNSS for the first time in the low-latitude region. [ABSTRACT FROM AUTHOR]

Zhejiang province in China experienced an extreme climate phenomenon in August 2014 with temperature rises, sunshine duration decreases, and precipitation increases, particularly, the successive heavy rainfall events occurring from 16 to 20 August 2014 that contributed to this climate anomaly. This study investigates the spatial-temporal variation characteristics of precipitable water vapor (PWV) and the normalized difference vegetation index (NDVI) associated with this phenomenon. Multiple sources of PWV values derived from the Global Positioning System (GPS), Radiosonde (RS) and European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim data are used with different spatiotemporal resolutions. The monthly averaged PWV in August 2014 exceeded the 95% percentiles of climatological value (53 mm) while the monthly averaged temperature was less than the 5% percentiles of climatological value (26.6 °C). Before the extreme precipitation, the PWV increased from the yearly averaged value of about 35 mm to more than 60 mm and gradually returned to the August climatological average of 50 mm after the precipitation ended. A large-scale atmospheric water vapor was partially conveyed by the warm wet air current of anticyclones which originated over the South China Sea (25° N, 130° E) and the Western Pacific Ocean. The monthly NDVI variation over the past 34 years (1982–2015) was investigated in this paper and the significant impact of extreme climate on vegetation growth in August 2014 was found. The extreme negative temperature anomaly and positive PWV anomaly are the major climate-driven factors affecting vegetation growth in the north and south of Zhejiang province with correlation coefficients of 0.83 and 0.72, respectively, while the extreme precipitation does not show any apparent impact on NDVI. [ABSTRACT FROM AUTHOR]

Microelectromechanical Systems (MEMS) Inertial Measurement Unit (IMU) containing a three-orthogonal gyroscope and three-orthogonal accelerometer has been widely utilized in position and navigation, due to gradually improved accuracy and its small size and low cost. However, the errors of a MEMS IMU based standalone Inertial Navigation System (INS) will diverge over time dramatically, since there are various and nonlinear errors contained in the MEMS IMU measurements. Therefore, MEMS INS is usually integrated with a Global Positioning System (GPS) for providing reliable navigation solutions. The GPS receiver is able to generate stable and precise position and time information in open sky environment. However, under signal challenging conditions, for instance dense forests, city canyons, or mountain valleys, if the GPS signal is weak and even is blocked, the GPS receiver will fail to output reliable positioning information, and the integration system will fade to an INS standalone system. A number of effects have been devoted to improving the accuracy of INS, and de-nosing or modelling the random errors contained in the MEMS IMU have been demonstrated to be an effective way of improving MEMS INS performance. In this paper, an Artificial Intelligence (AI) method was proposed to de-noise the MEMS IMU output signals, specifically, a popular variant of Recurrent Neural Network (RNN) Long Short Term Memory (LSTM) RNN was employed to filter the MEMS gyroscope outputs, in which the signals were treated as time series. A MEMS IMU (MSI3200, manufactured by MT Microsystems Company, Shijiazhuang, China) was employed to test the proposed method, a 2 min raw gyroscope data with 400 Hz sampling rate was collected and employed in this testing. The results show that the standard deviation (STD) of the gyroscope data decreased by 60.3%, 37%, and 44.6% respectively compared with raw signals, and on the other way, the three-axis attitude errors decreased by 15.8%, 18.3% and 51.3% individually. Further, compared with an Auto Regressive and Moving Average (ARMA) model with fixed parameters, the STD of the three-axis gyroscope outputs decreased by 42.4%, 21.4% and 21.4%, and the attitude errors decreased by 47.6%, 42.3% and 52.0%. The results indicated that the de-noising scheme was effective for improving MEMS INS accuracy, and the proposed LSTM-RNN method was more preferable in this application. [ABSTRACT FROM AUTHOR]

One of the key goals of geodesy is to study the fine structure of the Earth's gravity field and construct a high-resolution gravity field model (GFM). Aiming at the current insufficient resolution problem of the EIGEN_6C4 model, the refined ultra-high degree models EIGEN_3660 and EIGEN_5480 are determined with a spectral expansion approach in this study, which is to augment EIGEN_6C4 model using topographic potential models (TPMs). A comparative spectral evaluation for EIGEN_6C4, EIGEN_3660, and EIGEN_5480 models indicates that the gravity field spectral powers of EIGEN_3660 and EIGEN_5480 models outperform the EIGEN_6C4 model after degree 2000. The augmented models EIGEN_3660 and EIGEN_5480 are verified using the deflection of the vertical (DOV) of China and Colorado, gravity data from Australia and mainland America, and GNSS/leveling in China. The validation results indicate that the accuracy of EIGEN_3660 and EIGEN_5480 models in determining height anomaly, DOV, and gravity anomaly outperform the EIGEN_6C4 model, and the EIGEN_5480 model has optimal accuracy. The accuracy of EIGEN_5480 model in determining south–north component and east–west component of the DOV in China has been improved by about 21.1% and 23.1% compared to the EIGEN_6C4 model, respectively. In the mountainous Colorado, the accuracy of EIGEN_5480 model in determining south–north component and east–west component of the DOV has been improved by about 28.2% and 35.2% compared to EIGEN_6C4 model, respectively. In addition, gravity value comparison results in Australia and mainland America indicate that the accuracy of the EIGEN_5480 model for deriving gravity anomalies is improved by 16.5% and 11.3% compared to the EIGEN_6C4 model, respectively. [ABSTRACT FROM AUTHOR]

The global tropospheric zenith delay grid products of VMF1 and VMF3 (Vienna mapping functions) with different resolutions are used to calculate the tropospheric zenith delay of eight IGS (International GNSS Service) stations in China, and the accuracy of the two products under different interpolation methods is analyzed. As a result, the accuracy of utilizing different interpolation methods shows no obvious differences. The interpolation accuracy of the VMF3 grid model is slightly higher than that of the VMF1, and the interpolation accuracy of tropospheric delay is related to the elevation difference of grid points. In addition, according to ERA5 (the fifth generation of the Global Climate Information Analysis data set), the atmospheric stratification tropospheric delay is obtained, and a ZTD (the zenith tropospheric delay) height change grid model is constructed using the least squares exponential fitting method. The accuracy of the model is verified using the tropospheric delay product provided by the IGS. Finally, the constructed ZTD height change grid model is used as ZTD height reduction to solve the problem of large tropospheric delay errors in the VMF interpolation when the height change is large. The model accuracy of URUM station improve from 96.47 mm.to 32.97 mm (65.82%). [ABSTRACT FROM AUTHOR]

Precipitable water vapor (PWV) is a vitally important factor in atmospheric circulation. PWV is significant for forecasting extreme weather and understanding the dynamics of climate change. Comprehensively evaluating the performance of newly proposed remotely sensed water vapor products is crucial for guaranteeing their suitability for futural PWV applications. In this study, the accuracy of the recently established remotely sensed water vapor product from the Visible Infrared Imaging Radiometer Suite (VIIRS) satellite sensor on the Suomi National Polar-orbiting Partnership (SNPP) (VIIRS-PWV) platform within various regions of mainland China was evaluated via the PWV from Global Positioning System (GPS) observations. The GPS-derived PWV from 231 stations of the Crustal Movement Observation Network of China (CMONOC) from 2012 to 2018 was obtained through precise point positioning (PPP) techniques. The results showed that the mean value of the correlation coefficient (CC), the mean bias (MB), and the root-mean-square error (RMSE) between the VIIRS-PWV and the GPS-PWV were 0.92, −1.6 mm, and 4.7 mm, respectively. These values were comparable with the results of the PWV data derived from the Atmospheric Infrared Sounder (AIRS) and Moderate Resolution Imaging Spectroradiometer (MODIS) products. This indicates that the VIIRS product could provide PWV data with satisfactory accuracy for large-area scientific applications. Moreover, the MB and RMSE of the differences between the GPS-PWV and VIIRS-PWV showed obvious seasonal variations. The VIIRS-PWV generally performed better in winter (with the MB and RMSE values of 0.1 mm and 2.3 mm) than in summer (with the MB and RMSE values of −4.4 and 7.0 mm). Analysis among different regions revealed that the Central South (CS) region of China attained the largest mean RMSE value of 6.3 mm, and the North West (NW) region attained the smallest mean RMSE value of 3.8 mm. In addition, the southern region of China obtained a mean RMSE value of 5.6 mm, while that for the northern region of China was 3.9 mm. This indicates that the VIIRS-PWV has better accuracy within the northern region of China than within the southern region. [ABSTRACT FROM AUTHOR]

The signal measured by the maglev gyro sensor is sensitive to the influence of the instantaneous disturbance torque caused by the instantaneous strong wind or the ground vibration, which reduced the north-seeking accuracy of the instrument. To address this issue, we proposed a novel method combining the heuristic segmentation algorithm (HSA) and the two-sample Kolmogorov-Smirnov (KS) test (named HSA-KS method) to process the gyro signals and improve the north-seeking accuracy of the gyro. There were two key steps in the HSA-KS method: (i) all the potential change points were automatically and accurately detected by HSA, and (ii) the jumps in the signal caused by the instantaneous disturbance torque were quickly located and eliminated by the two-sample KS test. The effectiveness of our method was verified through a field experiment on a high-precision global positioning system (GPS) baseline at the 5th sub-tunnel of the Qinling water conveyance tunnel of the Hanjiang-to-Weihe River Diversion Project in Shaanxi Province, China. Our results from the autocorrelograms indicated that the jumps in the gyro signals can be automatically and accurately eliminated by the HSA-KS method. After processing, the absolute difference between the gyro and high-precision GPS north azimuths was enhanced by 53.5%, which was superior to the optimized wavelet transform and the optimized Hilbert-Huang transform. [ABSTRACT FROM AUTHOR]

The precise point positioning service on B2b signal (PPP-B2b) is a real-time decimeter-level positioning service provided by the BeiDou-3 Global Navigation Satellite System (BDS-3). The service provides users with high-precision orbit and clock corrections through geostationary orbit (GEO) satellites, which means that the PPP-B2b service would be unusable if GEO satellites were blocked. In this study, the performance of PPP-B2b corrections and real-time positioning results during outages of the PPP-B2b service are comprehensively investigated. The results showed that PPP can achieve satisfactory accuracy during outages of the PPP-B2b service by extending the nominal validity of the received PPP-B2b corrections. After extending the PPP-B2b corrections for 10 min, for BDS-3 medium earth orbit (MEO) satellites, the mean root-mean-square error (RMSE) values of the extended orbit were 0.16 m, 0.26 m, and 0.23 m in the radial, along-, and cross-track directions, respectively. The accuracy of the BDS-3 inclined geostationary orbit (IGSO) satellites was slightly worse than that of the BDS-3 MEO satellites; for Global Positioning System (GPS) satellites, the mean RMSE values of the extended orbit were 0.11 m, 0.45 m, and 0.33 m in the radial, along-, and cross-track directions, respectively. In terms of the extended clock, the mean standard deviation (STD) reached 0.17 ns, 0.20 ns, and 0.22 ns after 10 min for the BDS-3 MEO, BDS-3 IGSO, and GPS satellites, respectively. The positioning performance maintained with the extended corrections during the PPP-B2b service outage was evaluated based on five stations in and around China. Our experiments showed that, as long as the interruption time does not exceed 10 min, the real-time positioning with extended PPP-B2b corrections can achieve a comparable accuracy with that obtained following PPP-B2b correction. [ABSTRACT FROM AUTHOR]

The advancements of Earth observations, remote sensing, communications and navigation augmentation based on low Earth orbit (LEO) platforms present strong requirements for accurate, real-time and autonomous navigation of LEO satellites. Precise onboard real-time orbit determination (RTOD) using the space-borne data of multiple global navigation satellite systems (multi-GNSS) becomes practicable along with the availability of multi-GNSS. We study the onboard RTOD algorithm and experiments by using America's Global Positioning System (GPS) and China's regional BeiDou Navigation Satellite System (BDS-2) space-borne data of the FengYun-3C satellite. A new pseudo-ambiguity parameter, which combines the constant phase ambiguity, the orbit and clock offset error of the GPS/BDS broadcast ephemeris in the line-of-sight (LOS), is defined and estimated in order to reduce the negative effect of the LOS error on onboard RTOD. The analyses on the variation of the LOS error in the GPS/BDS broadcast ephemeris indicate that the pseudo-ambiguity parameter could be modeled as a random walk, and the setting of the power spectral density in the random walk model decides whether the pseudo-ambiguity can be estimated reasonably and the LOS error could be reduced or not. For different types of GPS/BDS satellites, the LOS errors show different variation characteristics, so the power spectral density should be set separately and differently. A numerical search approach is presented in this paper to find the optimal setting of the power spectral density for each type of GPS/BDS satellite by a series of tests. Based on the optimal stochastic model, a 3-dimensional (3D) real-time orbit accuracy of 0.7–2.0 m for position and 0.7–1.7 mm/s for velocity could be achieved only with dual-frequency BDS measurements and the broadcast ephemeris, while a notably superior orbit accuracy of 0.3–0.5 m for position and 0.3–0.5 mm/s for velocity is achievable using dual-frequency GPS/BDS measurements, due to the absorption effect of the estimated pseudo-ambiguity on the LOS error of the GPS/BDS broadcast ephemeris. Compared to using GPS-alone data, the GPS/BDS fusion only marginally improves the onboard RTOD orbit accuracy by about 1–3 cm, but the inclusion of BDS satellites increases the number of the tracked GNSS satellites and thus speeds up the convergence of the filter. Furthermore, the GPS/BDS fusion could help suppress the local orbit errors, ensure the orbit accuracy and improve the reliability and availability of the onboard RTOD when fewer GPS satellites are tracked in some anomalous arcs. [ABSTRACT FROM AUTHOR]

The process of soil freezing and thawing refers to the alternating phase change of liquid water and solid water in the soil, accompanied by a large amount of latent heat exchange. It plays a vital role in the land water process and is an important indicator of climate change. The Tibetan Plateau in China is known as the "roof of the world", and it is one of the most prominent physical characteristics is the freezing and thawing process of the soil. For the first time, this paper utilizes the spaceborne GNSS-R mission, i.e., CYGNSS (Cyclone Global Navigation Satellite System), to study the feasibility of monitoring the soil freeze-thaw (FT) cycles on the Tibetan Plateau. In the theoretical analysis part, model simulations show that there are abrupt changes in soil permittivities and surface reflectivities as the soil FT occurs. The CYGNSS reflectivities from January 2018 to January 2020 are compared with the SMAP FT state. The relationship between CYGNSS reflectivity and SMAP soil moisture within this time series is analyzed and compared. The results show that the effect of soil moisture on reflectivity is very small and can be ignored. The periodic oscillation change of CYGNSS reflectivity is almost the same as the changes in SMAP FT data. Freeze-thaw conversion is the main factor affecting CYGNSS reflectivity. The periodical change of CYGNSS reflectivity in the 2 years indicates that it is mainly caused by soil FT cycles. It is feasible to use CYGNSS to monitor the soil FT cycles in the Tibetan Plateau. This research expands the current application field of CYGNSS and opens a new chapter in the study of cryosphere using spaceborne GNSS-R with high spatial-temporal resolution. [ABSTRACT FROM AUTHOR]

Surface mass variations inferred from the Global Positioning System (GPS), and observed by the Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-On (GFO) complement each other in terms of spatial and temporal coverage. This paper presents an analysis of regional surface mass variations inverted from GPS vertical displacements under different density distributions of GPS stations, and compares the GPS-derived mass variations with GRACE/GFO inversion results in spatial and temporal domains. To this end, GPS vertical displacement data from a total of 85 permanent GPS stations of the Crustal Movement Observation Network of China (CMONOC), the latest GRACE/GFO RL06 spherical harmonic (SH) solutions and GRACE RL06 mascon solutions are used to investigate surface mass variations in four regions or basins, including the Yunnan Province (YNP), Min River Basin (MRB), Jialing River Basin (JLRB), and Wu River Basin (WRB) in Southwest China. Our results showed that the spatial distributions and seasonal characteristics of GPS-derived mass change time series agree well with those from GRACE/GFO observations, especially in regions with relatively dense distributions of GPS stations (e.g., in the YNP and MRB), but there are still obvious discrepancies between the GPS and GRACE/GFO results. Scale factor methods (both basin-scaled and pixel-scaled) were employed to reduce the amplitude discrepancies between GPS and GRACE/GFO results. The results also showed that the one-year gap between the GRACE and GFO missions can be bridged by scaled GPS-derived mass change time series in the four studied regions, especially in the YNP and MRB regions (with relatively dense distributions of GPS stations). [ABSTRACT FROM AUTHOR]

Most of the construction machinery for vibro-sinking stone columns, which are widely used in China, needs to be improved in terms of degree of automation. Engineering quality control is mainly carried out post-inspection; consequently, it is difficult to control the construction quality in real time. According to the construction characteristics of traditional stone column machines, we established the theory and model for the real-time monitoring of stone column construction, as well as put forward an intelligent monitoring method for stone column machines. With the comprehensive application of critical technologies such as the Global Navigation Satellite System (GNSS) measurement technology, laser ranging sensors, and massive data processing, an intelligent data acquisition technique and associated monitoring equipment for stone column construction machines are developed. The data acquisition and storage of crucial construction parameters, such as pile depth, pile point co-ordinates, bearing layer current, and reverse insertion times, are realized. A large number of actual construction data are collected and the construction quality parameters of stone column machines are obtained. By comparison with third-party detection data, it is verified that the intelligent monitoring technique for stone column machines proposed in this paper is feasible. [ABSTRACT FROM AUTHOR]

Mining subsidence is time-dependent and highly nonlinear, especially in the Loess Plateau region in Northwestern China. As a consequence, and mainly in building agglomerations, the structures can be damaged severely during or after underground extraction, with risks to human life. In this paper, we propose an approach based on a combination of a differential interferometric synthetic aperture radar (DInSAR) technique and a support vector machine (SVM) regression algorithm optimized by grid search (GS-SVR) to predict mining subsidence in a timely and cost-efficient manner. We consider five Advanced Land Observing Satellite (ALOS)/Phased Array type L-band Synthetic Aperture Radar (PALSAR) images encompassing the Dafosi coal mine area in Binxian and Changwu counties, Shaanxi Province. The results show that the subsidence predicted by the proposed InSAR and GS-SVR approach is consistent with the Global Positioning System (GPS) measurements. The maximum absolute errors are less than 3.1 cm and the maximum relative errors are less than 14%. The proposed approach combining DInSAR with GS-SVR technology can predict mining subsidence on the Loess Plateau of China with a high level of accuracy. This research may also help to provide disaster warnings. [ABSTRACT FROM AUTHOR]

The Ranging Mode (R-Mode), a maritime terrestrial navigation system under development, is a promising approach to increase the resilient provision of position, navigation and timing (PNT) information for bridge instruments, which rely on Global Navigation Satellite Systems (GNSS). The R-Mode utilizes existing maritime radio infrastructure such as marine radio beacons, which support maritime traffic with more reliable and accurate PNT data in areas with challenging conditions. This paper analyzes the potential service, which the R-Mode could provide to the mariner if worldwide radio beacons were upgraded to broadcast R-Mode signals. The authors assumed for this study that the R-Mode is available in the service area of the 357 operational radio beacons. The comparison with the maritime traffic, which was generated from a one-day worldwide Automatic Identification System (AIS) Class A dataset, showed that on average, 67% of ships would operate in a global R-Mode service area, 40% of ships would see at least three and 25% of ships would see at least four radio beacons at a time. This means that R-Mode would support 25% to 40% of all ships with position and 67% of all ships with PNT integrity information. The relatively high number of supported ships compared to the total radio beacon coverage of about 9% of the earth's surface is caused by the good coverage of busy ports and areas such as the coast of China, North Sea and Baltic Sea. These numbers emphasize the importance of marine radio beacons for the R-Mode system. [ABSTRACT FROM AUTHOR]

Southeast China, a non-core region influenced by the El Niño–Southern Oscillation (ENSO), has been seldom investigated before. However, the occurrence of ENSO will affect the redistribution of precipitation and the temperature (T) spatial pattern on a global scale. This condition will further lead to flood or drought disasters in Southeast China. Therefore, the method of monitoring the occurrence of ENSO is important and is the focus of this paper. The spatiotemporal characteristics of precipitable water vapor (PWV) and T are first analyzed during ENSO using the empirical orthogonal function (EOF). The results showed that a high correlation spatiotemporal consistency exist between PWV and T. The response thresholds of PWV and T to ENSO are determined by moving the window correlation analysis (MWCA). If the sea surface temperature anomaly (SSTA) at the Niño 3.4 region exceeded the ranges of (−1.17°C, 1.04°C) and (−1.15°C, 1.09°C), it could cause the anomalous change of PWV and T in Southeast China. Multichannel singular spectral analysis (MSSA) is introduced to analyze the multi-type signals (tendency, period, and anomaly) of PWV and T over the period of 1979–2017. The results showed that the annual abnormal signal and envelope line fluctuation of PWV and T agreed well in most cases with the change in SSTA. Therefore, a standard PWV and T index (SPTI) is proposed on the basis of the results to monitor ENSO events. The PWV and T data derived from the grid-based European Center for Medium-Range Weather Forecasting (ECMWF) reanalysis products and GNSS/RS stations in Southeast China were used to validate the performance of the proposed SPTI. Experimental results revealed that the time series of average SPTI calculated in Southeast China corresponded well to that of SSTA with a correlation coefficient of 0.66 over the period of 1979–2017. The PWV values derived from the Global Navigation Satellite System (GNSS) and radiosonde data at two specific stations (WUHN and 45004) were also used to calculate the SPTI. The results showed that the correlation coefficients between SPTI and SSTA were 0.73 and 0.71, respectively. Such results indicate the capacity of the proposed SPTI to monitor the ENSO in Southeast China. [ABSTRACT FROM AUTHOR]