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]

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]

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]

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]

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]

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]

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]

Since 23 June 2020, BDS-3 has been entirely operated and obtained the ability of global PNT (Positioning, Navigation, and Timing) services. Afterward, real-time Precise Point Positioning (PPP) service is available in China's surrounding areas via BDS-3 PPP-B2b signal. However, such a real-time PPP service cannot maintain the high accuracy and continuity of positioning solutions in challenging environments, such as urban environments. For that, we carried out a model by integrating between-satellite single-differenced (BSSD) PPP, a low-cost Inertial Navigation System (INS), and an odometer via an extended Kalman filter. The performance of this integration model was assessed with vehicle-borne data. Results demonstrated that (1) the position RMS (Root Mean Square) of BSSD PPP are 64.33 cm, 53.47 cm, and 154.11 cm. Compared with BSSD PPP, about 31.2%, 23.3%, and 27.3% position improvements can be achieved by using INS. Further enhancements of position RMS benefiting from the odometer are 1.34%, 1.41%, and 1.73% in the three directions. (2) Anyway, the accuracy of BSSD PPP/INS/Odometer tightly coupled integration is slightly higher than that of undifferenced PPP/INS/Odometer integration, with average improvement percentages of 7.71%, 3.09%, and 0.27%. Meanwhile, the performance of BSSD PPP/INS/Odometer integration during the periods with satellite outages is better than the undifferenced PPP-based solutions. (3) The improvements in attitudes from an odometer are more significant on heading angle than the other two attitudes, with percentages of 25.00%. (4) During frequent GNSS outage periods, the reduction in average maximum position drifts provided by INS are 18.01%, 8.95%, and 20.74%. After integrating with an odometer, the drifts can be furtherly decreased by 25.11%, 15.96%, and 20.69%. For attitude, about 41.67% reduction in average maximum drifts of heading angles is obtained. [ABSTRACT FROM AUTHOR]

Atmospheric precipitable water vapor (PWV) is a key variable for weather forecast and climate research. Various techniques (e.g., radiosondes, global navigation satellite system, satellite remote sensing and reanalysis products by data assimilation) can be used to measure (or retrieve) PWV. However, gathering PWV data with high spatial and temporal resolutions remains a challenge. In this study, we propose a new empirical PWV grid model (called ASV-PWV) using the zenith wet delay from the Askne model and improved by the spherical harmonic function and vertical correction. Our method is convenient and enables the user to gain PWV data with only four input parameters (e.g., the longitude and latitude, time, and atmospheric pressure of the desired position). Profiles of 20 radiosonde stations in Qinghai Tibet Plateau, China, along with the latest publicly available C-PWVC2 model are used to validate the local performance. The PWV data from ASV-PWV and C-PWVC2 is generally consistent with radiosonde (the average annual bias is −0.44 mm for ASV-PWV and −1.36 mm for C-PWVC2, the root mean square error (RMSE) is 3.44 mm for ASV-PWV and 2.51 mm for C-PWVC2, respectively). Our ASV-PWV performs better than C-PWVC2 in terms of seasonal characteristics. In general, a sound consistency exists between PWV values of ASV-PWV and the fifth generation of European Centre for Medium-Range Weather Forecasts Atmospheric Reanalysis (ERA5) (total 7381 grid points in 2020). The average annual bias and RMSE are −0.73 mm and 4.28 mm, respectively. ASV-PWV has a similar performance as ERA5 reanalysis products, indicating that ASV-PWV is a potentially alternative option for rapidly gaining PWV. [ABSTRACT FROM AUTHOR]

Global navigation satellite system (GNSS) signals are affected by refraction when traveling through the troposphere, which result in tropospheric delay. Generally, the tropospheric delay is estimated as an unknown parameter in GNSS data processing. With the increasing demand for GNSS real-time applications, high-precision tropospheric delay augmentation information is vital to speed up the convergence of PPP. In this research, we estimate the zenith tropospheric delay (ZTD) from 2018 to 2019 by static precise point positioning (PPP) using the fixed position mode; GNSS observations were obtained from the National Geomatics Center of China (NGCC). Firstly, ZTD outliers were detected, and data gaps were interpolated using the K-nearest neighbor algorithm (KNN). Secondly, The ZTD differences between the KNN and periodic model were employed as input datasets to train the long short-term memory (LSTM) neural network. Finally, LSTM forecasted ZTD differences and the ZTD periodic signals were combined to recover the final forecasted ZTD results. In addition, the forecasted ZTD results were applied in static PPP as a prior constraint to reduce PPP convergence time. Numerical results show that the average root-mean-square error (RMSE) of predicting ZTD is about 1 cm. The convergence time of the PPP which was corrected by the LSTM-ZTD predictions is reduced by 13.9, 22.6, and 30.7% in the summer, autumn, and winter, respectively, over GPT2-ZTD corrected PPP and unconstrained conventional PPP for different seasons. [ABSTRACT FROM AUTHOR]

Toxic metals have attracted great concern worldwide due to their toxicity and slow decomposition. Although metal concentrations can be accurately obtained with chemical methods, it is difficult to map metal distributions on a large scale due to their inherently low efficiency and high cost. Moreover, chemical analysis methods easily lead to secondary contamination. To address these issues, 110 topsoil samples were collected using a soil sampler, and positions for each sample were surveyed using a global navigation satellite system (GNSS) receiver from a coal mine in northern China. Then, the metal contents were surveyed in a laboratory via a portable X-ray fluorescence spectroscopy (XRF) device, and GaoFen-5 (GF-5) satellite hyperspectral images were used to retrieve the spectra of the soil samples. Furthermore, a Savitzky–Golay (SG) filter and continuous wavelet transform (CWT) were selected to smooth and enhance the soil reflectance. Competitive adaptive reweighted sampling (CARS) and Boruta algorithms were utilized to identify the feature bands. The optimum two-stage method, consisting of the random forest (RF) and ordinary kriging (OK) methods, was used to infer the metal concentrations. The following outcomes were achieved. Firstly, both zinc (Zn) (68.07 mg/kg) and nickel (Ni) (26.61 mg/kg) surpassed the regional background value (Zn: 48.60 mg/kg, Ni: 19.5 mg/kg). Secondly, the optimum model of RF, combined with the OK (RFOK) method, with a relatively higher coefficient of determination (R2) (R2 = 0.60 for Zn, R2 = 0.30 for Ni), a lower root-mean-square error (RMSE) (RMSE = 12.45 mg/kg for Zn, RMSE = 3.97 mg/kg for Ni), and a lower mean absolute error (MAE) (MAE = 9.47 mg/kg for Zn, MAE = 3.31mg/kg for Ni), outperformed the other four models, including the RF, OK, inverse distance weighted (IDW) method, and the optimum model of RF combined with IDW (RFIDW) method in estimating soil Zn and Ni contents, respectively. Thirdly, the distribution of soil Zn and Ni concentrations obtained from the best-predicted method and the GF-5 satellite hyperspectral images was in line with the actual conditions. This scheme proves that satellite hyperspectral images can be used to directly estimate metal distributions, and the present study provides a scientific base for mapping heavy metal spatial distribution on a relatively large scale. [ABSTRACT FROM AUTHOR]

The Loess Plateau is one of the three most severely affected geological disaster areas in China. Water sensitivity is the most significant feature of the loess. Under the action of continuous heavy rainfall, rainwater infiltrates the loess, resulting in a rapid increase in soil saturation and changes in soil moisture. This affects the shear strength of the soil and induces shallow loess landslides. Therefore, it is significant to our country's disaster prevention and mitigation efforts to effectively monitor the deformation and inducement of such landslides. At present, the global navigation satellite system (GNSS) is widely used in the field of landslide disaster monitoring as a technical means to directly obtain real-time three-dimensional vector deformation of the surface. At the same time, GNSS can also provide a steady stream of L-band microwave signals to obtain surface environmental information, such as soil moisture around the station. In past landslide disaster monitoring research, GNSS was only used to provide three-dimensional deformation information, and its ability to provide environmental information around the station was almost completely ignored. This study proposes a ground-based GNSS remote sensing comprehensive monitoring system integrating "three-dimensional deformation and soil moisture content" combined with a rainfall-type shallow loess landslide event in Linxia City. The ability of ground-based GNSS to comprehensively monitor shallow loess landslide disasters was analysed. Experiments show that GNSS can provide high-precision deformation time series characteristics and monitor the changes in soil moisture content around the station at the same time; the two have a certain response relationship, which can comprehensively evaluate the stability of shallow loess landslides. As heavy rainfall is a key factor affecting the change in soil water content, this study adds the atmospheric water vapour content calculated by ground-based GNSS refraction remote sensing in the discussion chapter and analyses the relationship between precipitable water vapour and rainfall in this area to give full play to ground-based GNSS remote sensing. In the role of landslide disaster monitoring, we hope to build a more comprehensive ground-based GNSS remote sensing monitoring system to better serve the monitoring of landslide disasters. [ABSTRACT FROM AUTHOR]

As the largest groundwater drainage region in China, the per capita water resources in the North China Plain (NCP) account for only one-seventh of the country's available water resources. Currently, the NCP is experiencing a serious water shortage due to the overexploitation of groundwater resources and a subsequent series of natural disasters. Thus, accurate regional assessments and effective water resource management policies are of critical importance. To accomplish this phenomenon, the daily terrestrial water storage anomaly (TWSA) over the NCP is calculated from the combination of the GNSS vertical deformation sequences (seasonal items) and GRACE (trend items). The groundwater storage anomaly (GWSA) in the NCP is obtained by subtracting the canopy water, soil water, and snow water equivalent components from the TWSA. The inversion results of this study are verified by comparisons with the Global Land Data Assimilation System (GLDAS) data products. The elevated annual amplitude areas are located in Beijing and Tianjin, and the Pearson correlation coefficient (PCC), root mean square error (RMSE), and Nash–Sutcliffe efficiency (NSE) between the two GWSA results are 0.67, 4.01 cm, and 0.61, respectively. This indicates that the methods proposed in this study are reliable. Finally, the groundwater drought index was calculated for the period from 2011 to 2021, and the results showed that 2019 was the driest year, with a drought severity index value of −0.12, indicative of slightly moderate drought conditions. By calculating and analyzing the annual GWSA, this work shows that the South–North Water Transfer Project does provide some regional drought mitigation. [ABSTRACT FROM AUTHOR]

The Liupan Mountain (LPM) area is located on the northeast margin of the Qinghai–Tibet Plateau, which is the western part of the second ladder of China's terrain. It is also an intersection area of two air currents, which are caused by the combined action of the Qinghai–Tibet Plateau, the middle and lower levels of the westerly belt and the edge of the monsoon area. LPM is one of the main air water vapour transport pathways in Northwest China as well as a main water conservation area for nearly ten million people. Research on atmospheric precipitable water vapour (PWV) variation characteristics in LPM is beneficial for understanding the mechanisms of orographic precipitation and improving the effects of weather modification. Based on the data from 10 Global Navigation Satellite System Meteorology (GNSS/MET) stations for 6 years and the data of automatic weather stations in the LPM, the temporal and spatial variation characteristics of PWV in the LPM were analysed, and the differences in PWV in the 24 h before and after precipitation were compared in this study. The results showed that the hourly, monthly and seasonal variations in PWV displayed obvious patterns. PWV increased rapidly 10 h ahead of precipitation, while it decreased rapidly within 10 h after precipitation, which was slower than that before precipitation. In terms of spatial distribution, PWV was larger in the south than in the north and larger in the east than in the west. Although the precipitation on the LPM peak was the highest in the whole LPM area, its PWV was always the lowest, indicating that the PWV was obviously affected by the air temperature. This showed that under the same water vapour condition, precipitation was more likely to form in the area with low temperature, and the precipitation was larger, which also provided a train of thought for improving the method of artificial precipitation enhancement by using the condensation catalyst. [ABSTRACT FROM AUTHOR]

On the 15th of March 2021, the strongest sandstorm in a decade occurred in northern China, and had a great adverse impact on the natural environment and human health in northern China. Real-time monitoring of dust storms is becoming increasingly important. In order to effectively analyze the non-hydrostatic delay (ZNHD) anomaly during a sandstorm, the method based on GNSS-derived tropospheric ZNHD residual to monitor the sandstorm is proposed at the same time. We studied the relationship between ZNHD/PWV and PM10/PM2.5 in Beijing, Changchun, Pingliang and Zhongwei before and after sandstorms. The ZNHD time series was then decomposed by singular spectrum analysis (SSA) and the residuals were obtained. The relationship between the GNSS-derived ZNHD residual and PM10 was analyzed. The results show that the impact of the sandstorm on PM10 is greater than that on PM2.5. Before the sandstorm, the correlation between PM10 and ZNHD was low, less than 0.25. When the sandstorm occurred, the correlation between PM10 and ZNHD increased significantly, and the maximum was greater than 0.7. When the sandstorm ended, the correlation between PM10 and ZNHD decreased significantly. Through the relationship between the ZNHD residual and PM10, it can be found that when the peak-to-peak values of the ZNHD residual are all above 80 mm, sandstorms may occur. But Rainfall, snowfall, haze and other abnormal weather can also lead to ZNHD anomalies. [ABSTRACT FROM AUTHOR]

Tropospheric delay is one of the key factors that may affect high-precision satellite navigation and positioning and water vapor retrieval performance. Its variation in the vertical direction is much greater than that in the horizontal direction. Therefore, the vertical profile model of zenith total delay (ZTD) is important for the spatial interpolation of high-precision ZTD products and the development of ZTD models. However, in China, low spatial and temporal resolutions remain persistent in ZTD vertical profile models and limit their application. In this study, ZTD vertical profile grid models (CZTD-H model: CZTD-HM and CZTD-HE models) were developed by considering the time-varying height scale factor for China and employing ZTD layered profile information with high temporal-spatial resolution calculated using MERRA-2 data and ERA5 data based on the integration method during 2012–2016. The CZTD-H model accuracy was verified using the global navigation satellite system (GNSS) data acquired from the Crustal Movement Observation Network of China (CMONOC) and radiosonde data as reference and was compared with the canonical GPT3 model accuracy. The applicability of CZTD-HM and CZTD-HE models was discussed. The results showed that: (1) CZTD-HM and CZTD-HE models exhibited excellent performance for ZTD layered vertical interpolation in northwestern and southeastern China, respectively. Among ZTD layered profiles from 84 radiosonde stations, the RMSE of ZTD data interpolated using CZTD-HM and CZTD-HE models improved by 12.9/16.23% and 13.8/17.16% compared with GPT3-1 and GPT3-5 models, respectively. (2) The CZTD-H model maintained high performance for the spatial interpolation of GGOS grid ZTD data. Validation with ZTD data from 249 GNSS stations showed that the RMSEs of both CZTD-HM and CZTD-HE models improved by 2.8 mm (19.7%) and 2.6 mm (18.6%) compared with those of the GPT3-1 and GPT3-5 models, respectively. The CZTD-HE model showed excellent performance in summer among all the models. Only the location and day of the year were required for the application of the CZTD-H model, which showed excellent ZTD vertical correction performance in China. With the different performances of the CZTD-HE and CZTD-HM models in China, the ERA5 model can be recommended for practical applications. Therefore, these results can provide a reference for the data source selection of ZTD vertical profile model construction on the basis of high-precision reanalysis data, GNSS real-time high-precision positioning, and GNSS meteorology in China. [ABSTRACT FROM AUTHOR]

The weighted mean temperature (Tm) is a vital parameter for converting zenith wet delay (ZWD) into precipitation water vapor (PWV) and plays an essential part in the Global Navigation Satellite System (GNSS) inversion of PWV. To address the inability of current mainstream models to fit the nonlinear relationship between Tm and meteorological and spatiotemporal factors, whose accuracy is limited, a weighted mean temperature model using the random forest (named RFTm) was proposed to enhance the accuracy of the Tm predictions in mainland China. The validation with the Tm from 84 radiosonde stations in 2018 showed that the root mean square (RMS) of the RFTm model was reduced by 38.8%, 44.7%, and 35.5% relative to the widely used Global Pressure and Temperature 3 (GPT3) with 1° × 1°/5° × 5° versions and Bevis, respectively. The Bias and RMS of the new model in different latitude bands, various height intervals, and different times were significantly better than those of the other three comparative models. The accuracy of the new model presented a more stable adaptability. Therefore, this study provides a new idea for estimating Tm and can provide a more accurate Tm for GNSS meteorology. [ABSTRACT FROM AUTHOR]

Data fusion is an effective method to obtain high-precision and high-spatiotemporal-resolution precipitable water vapor (PWV) products, which play an important role in understanding climate change and meteorological monitoring. However, existing fusion methods have some shortcomings, such as ignoring the applicability of the model space or the high complexity of model operation. In this study, the high-precision and high-temporal-resolution Global Navigation Satellite System (GNSS) PWV was used to calibrate and optimize the ERA5 PWV product of the European Center for Medium-Range Weather Forecasts Reanalysis 5 (ERA5) with high spatial resolution to improve its accuracy, and its applicability was verified at the spatiotemporal scale. First, this study obtained accurate GNSS PWV from meteorological data from stations and used it as the true value to analyze the distribution of the ERA5 PWV in mainland China. The results showed that the ERA5 PWV showed significant spatial and temporal differences. Then, a backpropagation neural network (BPNN) fusion correction model with additional constraints was established. The correction results showed that the bias of the ERA5 PWV mainly fluctuated near 0, the correlation between the ERA5 PWV and GNSS PWV was increased to 0.99, and the positive improvement rate of the root-mean-square error (RMSE) was 95%. In the temporal scale validation, the RMSE of the ERA5 PWV decreased from 2.05 mm to 1.67 mm, an improvement of 18.54%. In the spatial scale validation, the RMSE of the four seasons decreased by 0.26–80% (spring), 0.28–70.71% (summer), 0.28–45.23% (autumn), and 0.30–40.75% (winter). Especially in the summer and plateau mountainous areas where the ERA5 PWV performance was poor, the model showed suitable stability. Finally, the fusion model was used to generate a new PWV product, which improved the accuracy of ERA5 PWV on the basis of ensuring the spatial resolution. [ABSTRACT FROM AUTHOR]

Extensive observation collection, unified and rigorous data processing, and accurate construction of the station motion model are the three essential elements for the accuracy and reliability of the Global Navigation Satellite System (GNSS) velocity field. GNSS data reprocessing not only can weaken the influence of untrue nonlinear site signals caused by imperfect models but also can eliminate the displacement offset caused by frame transformation, solution strategy, and model change. Based on the new repro3 criteria of the International GNSS Service (IGS), we process rigorously GNSS observations of continental China from the period 2000 to 2020 to refine GNSS station secular velocities and analyze the present-day crustal deformation in continental China. The main contributions of this work included the followings. Firstly, the repro3 algorithm and model are used to uniformly and rigorously process the two-decade GNSS historical observations to obtain more reliable GNSS coordinate time series with mm-level precision. Combined with the historical records of major earthquakes in continental China, we build a GNSS time series model considering nonlinear factors (velocity, offset, period, co-seismic/post-seismic deformation) to extract GNSS horizontal velocity field whose root mean square (RMS) mean is 0.1 mm/a. Secondly, the GNSS horizontal grid velocity field in continental China is interpolated using the gpsgridder method (the minimum radius is set to 16, and the Poisson's ratio is set to 0.5). Estimation and analysis of the crustal strain rate solution lead to the conclusion that the strain degree in West China (the high strain region is mainly located in the Qinghai Tibet Plateau and Tianshan Mountains) is much more intense than that in the east (the main strain rate is less than 5 nstrain/year). In addition, most strong earthquakes in the Chinese mainland occurred on active blocks and their boundary faults with large changes in the GNSS velocity field and strain field. Then, an improved K-means++ clustering analysis method is proposed to divide active blocks using GNSS horizontal velocity field. Furthermore, different relative motion models of different blocks are constructed using the block division results. Among them, the Eurasian block has the lowest accuracy (the RMS of residual velocity in the east and north directions are 5.60 and 9.65 mm/a, respectively), and the China block 7 has the highest accuracy (the RMS mean of relative velocity in the east and north directions are 2.60 and 2.65 mm/a, respectively). More observations (2260+ sites), longer time (20 years), and updated criteria (Repro3) are to finely obtain the GNSS velocity field in continental China, and depict crustal deformation and active block with the gpsgridder and improved K-means++ methods. [ABSTRACT FROM AUTHOR]