Analysis of the Spatial Distribution and Common Mode Error Correlation in a Small-Scale GNSS Network

When analyzing GPS time series, common mode errors (CME) often obscure the actual crustal movement signals, leading to deviations in the velocity estimates of station coordinates. Therefore, mitigating the impact of CME on station positioning accuracy is crucial to ensuring the precision and reliability of GNSS time series. The current approach to separating CME mainly uses signal filtering methods to decompose the residuals of the observation network into multiple signals, from which the signals corresponding to CME are identified and separated. However, this method overlooks the spatial correlation of the stations. In this paper, we improved the Independent Component Analysis (ICA) method by introducing correlation coefficients as weighting factors, allowing for more accurate emphasis or attenuation of the contributions of the GNSS network's spatial distribution during the ICA process. The results show that the improved Weighted Independent Component Analysis (WICA) method can reduce the root mean square (RMS) of the coordinate time series by an average of 27.96%, 15.23%, and 28.33% in the E, N, and U components, respectively. Compared to the ICA method, considering the spatial distribution correlation of stations, the improved WICA method shows enhancements of 12.53%, 3.70%, and 8.97% in the E, N, and U directions, respectively. This demonstrates the effectiveness of the WICA method in separating CMEs and provides a new algorithmic approach for CME separation methods.

Tracking the water storage and runoff variations in the Paraná basin via GNSS measurements

The Paraná basin is the second largest river basin in South America and provides abundant water resources globally. However, current research lacks hydrological investigation of the region. The vertical crustal deformation recorded by the Global Navigation Satellite System (GNSS) can be used to accurately estimate regional-scale terrestrial water storage (TWS). Therefore, we utilized the daily vertical displacement time series data at 102 GNSS stations to recover the water storage variations in the Paraná basin from 2013 to 2020. To recognize primary spatiotemporal features of TWS changes, we applied the principal component analysis (PCA) method in the inversion strategy. Results indicate that the TWS variations inferred from GNSS generally align in spatiotemporal patterns with estimates from both the Gravity Recovery and Climate Experiment (GRACE) and the Global Land Data Assimilation System (GLDAS). However, some discrepancies are evident at local scales. The TWS changes derived from both GNSS and GRACE exhibited generally larger magnitude of oscillations than those estimated by GLDAS, while the GRACE results neglected the evident seasonal oscillation of the water mass in the southeast of the basin. Given the challenge of capturing large-scale runoff variations through in-situ observations, we innovatively applied GNSS and water budget closure method to provide a novel runoff estimate for the Paraná basin. The GNSS-inferred runoff exhibited a strong correlation (correlation coefficient of 0.72) with in-situ observations. Overall, our study fills the critical knowledge gap in geodesy-based hydrological investigation in the Paraná basin. We aim to highlight the immense potential of GNSS for hydrological parameter estimation and provide valuable reference data for regional hydrological research and for water resources management.

Spatial-Temporal Dynamic Evolution of Land Deformation Driven by Hydrological Signals around Chaohu Lake

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.

A Review of GNSS/GPS in Hydrogeodesy: Hydrologic Loading Applications and Their Implications for Water Resource Research

Hydrogeodesy, a relatively new field within the earth sciences, is the analysis of the distribution and movement of terrestrial water at Earth's surface using measurements of Earth's shape, orientation, and gravitational field. In this paper, we review the current state of hydrogeodesy with a specific focus on Global Navigation Satellite System (GNSS)/Global Positioning System measurements of hydrologic loading. As water cycles through the hydrosphere, GNSS stations anchored to Earth's crust measure the associated movement of the land surface under the weight of changing hydrologic loads. Recent advances in GNSS-based hydrogeodesy have led to exciting applications of hydrologic loading and subsequent terrestrial water storage (TWS) estimates. We describe how GNSS position time series respond to climatic drivers, can be used to estimate TWS across temporal scales, and can improve drought characterization. We aim to facilitate hydrologists' use of GNSS-observed surface deformation as an emerging tool for investigating and quantifying water resources, propose methods to further strengthen collaborative research and exchange between geodesists and hydrologists, and offer ideas about pressing questions in hydrology that GNSS may help to answer.

Spatiotemporal Filtering for Continuous GPS Coordinate Time Series in Mainland China by Using Independent Component Analysis

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.

Strain Field Features and Three-Dimensional Crustal Deformations Constrained by Dense GRACE and GPS Measurements in NE Tibet

The continuing impact between the Eurasia Plate and India results in the thickening and shortening of the N-S Tibetan Plateau. There has been strong tectonic movement along the boundary of the zones of deformation of the NE corner of the Tibetan plateau (NET) since the new tectonic period, with its dynamic mechanisms remaining controversial. Here, we use observations of 39 Continuous Global Positioning System (CGPS) gauges and 451 Crustal Movement Observation Network of China (CMONOC) campaign-mode stations to detect the three-dimensional deformation of the crust in the NET. Improved processing procedures are implemented to strengthen the patterns of strain throughout the NET. The principal component analysis (PCA) technique is introduced to decompose the time series into spatial eigenvectors and principal components (PCs), and the first three PCs are used to estimate and rectify common mode errors (CMEs). In addition, GRACE observations are used to detect deformation changes that account for non-tidal oceanic mass loading, hydrological loading, and surface pressure. The rectified deformation of the crust indicates the anisotropic nature of both the subsidence and uplift, and that the highest uplift rate of the Longmen Shan fault uplift reaches 7.13 ± 0.53 mm/yr. Finally, the horizontal velocity is further used to enumerate the strain rates throughout the NET. The results show that the shear band retained property in line with the strike-slip fault along the Altyn Tagh fault, the Qilian Shan faults, the Haiyuan fault, the West Qinling fault, the East Kunlun fault, and the Longmen Shan fault. In addition, the results further indicate that the whole NET shows a strong relationship with the mean principal rates of horizontal shortening strain. Extension and compression of the crust reasonably describe its sinking and uplifting.

Missing Data Imputation in GNSS Monitoring Time Series Using Temporal and Spatial Hankel Matrix Factorization

GNSS time series for static reference stations record the deformation of monitored targets. However, missing data are very common in GNSS monitoring time series because of receiver crashes, power failures, etc. In this paper, we propose a Temporal and Spatial Hankel Matrix Factorization (TSHMF) method that can simultaneously consider the temporal correlation of a single time series and the spatial correlation among different stations. Moreover, the method is verified using real-world regional 10-year period monitoring GNSS coordinate time series. The Mean Absolute Error (MAE) and Root-Mean-Square Error (RMSE) are calculated to compare the performance of TSHMF with benchmark methods, which include the time-mean, station-mean, K-nearest neighbor, and singular value decomposition methods. The results show that the TSHMF method can reduce the MAE range from 32.03% to 12.98% and the RMSE range from 21.58% to 10.36%, proving the effectiveness of the proposed method.

A New Spatial Filtering Algorithm for Noisy and Missing GNSS Position Time Series Using Weighted Expectation Maximization Principal Component Analysis: A Case Study for Regional GNSS Network in Xinjiang Province

Common Mode Error (CME) presents a kind of spatially correlated error that is widespread in regional Global Navigation Satellite System (GNSS) networks and should be eliminated during postprocessing of a GNSS position time series. Several spatiotemporal filtering methods have been developed to mitigate the effects of CME. However, such methodologies become inappropriate when missing and noisy data exists. In this research, we introduce a novel spatial filtering algorithm called Weighted Expectation Maximization Principal Component Analysis (WEMPCA) for detecting and removing CME from noisy GNSS position time series with missing values, among which formal errors of daily GNSS solutions are utilized to weight the input data. Compared with traditional PCA and the special case of EMPCA, simulation experiments demonstrate that the new WEMPCA algorithm always has outstanding performance over others. The WEMPCA algorithm was then successfully used to extract the CME from real noisy and missing GNSS position time series in Xinjiang province. Our results show that only the first principal component exhibits significant spatial response, with average values of 70.11%, 66.53%, and 52.45% for North, East, and Up (NEU) components, respectively, indicating that it represents the CME of this region. After removing CME, the canonical correlation coefficients and root mean square error of GNSS residual time series, as well as the amplitudes of power-law noises (PLN), are obviously decreased in all three directions. However, the white noise (WN) amplitudes are found to diminish exclusively in the North and East component, not in the Up components. Moreover, the average velocity differences before and after filtering CME are 0.19 mm/year, 0.03 mm/year, and −0.56 mm/year for the NEU components, respectively, indicating that CME has an influence on the GNSS station velocity estimation. The velocity uncertainty is also reduced by 43.51%, 38.64%, and 40.39% on average for the NEU components, respectively, implying that the velocity estimates are more reliable and accurate after removing CME. Therefore, we conclude that the new WEMPCA approach provides an efficient solution to detect and mitigate CME from the noisy and missing GNSS position time series.

Assessment of Contemporary Antarctic GIA Models Using High-Precision GPS Time Series

Past redistributions of the Earth’s mass resulting from the Earth’s viscoelastic response to the cycle of deglaciation and glaciation reflect the process known as glacial isostatic adjustment (GIA). GPS data are effective at constraining GIA velocities, provided that these data are accurate, have adequate spatial coverage, and account for competing geophysical processes, including the elastic loading of ice/snow ablation/accumulation. GPS solutions are significantly affected by common mode errors (CMEs) and the choice of optimal noise model, and they are contaminated by other geophysical signals due primarily to the Earth’s elastic response. Here, independent component analysis is used to remove the CMEs, and the Akaike information criterion is used to determine the optimal noise model for 79 GPS stations in Antarctica, primarily distributed across West Antarctica and the Antarctic Peninsula. Next, a high-resolution surface mass variation model is used to correct for elastic deformation. Finally, we use the improved GPS solution to assess the accuracy of seven contemporary GIA forward models in Antarctica. The results show that the maximal GPS crustal displacement velocity deviations reach 4.0 mm yr−1, and the mean variation is 0.4 mm yr−1 after removing CMEs and implementing the noise analysis. All GIA model-predicted velocities are found to systematically underestimate the GPS-observed velocities in the Amundsen Sea Embayment. Additionally, the GPS vertical velocities on the North Antarctic Peninsula are larger than those on the South Antarctic Peninsula, and most of the forward models underestimate the GIA impact on the Antarctic Peninsula.

Potential Contributors to Common Mode Error in Array GPS Displacement Fields in Taiwan Island

The existence of the common mode error (CME) in the continuous global navigation satellite system (GNSS) coordinate time series affects geophysical studies that use GNSS observations. To understand the potential contributors of CME in GNSS networks in Taiwan and their effect on velocity estimations, we used the principal component analysis (PCA) and independent component analysis (ICA) to filter the vertical coordinate time series from 44 high-quality GNSS stations in Taiwan island in China, with a span of 10 years. The filtering effects have been evaluated and the potential causes of the CME are analyzed. The root-mean-square values decreased by approximately 14% and 17% after spatio-temporal filtering using PCA and ICA, respectively. We then discuss the relationship between the CME sources obtained by ICA and the environmental loads. The results reveal that the independent displacements extracted by ICA correlate with the atmospheric mass loading (ATML) and land water storage mass loading (LWS) of Taiwan in terms of both its amplitude and phase. We then use the white noise plus power law noise model to quantitatively estimate the noise characteristics of the pre- and post-filtered coordinate time series based on the maximum likelihood estimation criterion. The results indicate that spatio-temporal filtering reduces the amplitude of the PL and the periodic terms in the GPS time series.

Imputation of GPS Coordinate Time Series Using missForest

The global positioning system (GPS) can provide the daily coordinate time series to help geodesy and geophysical studies. However, due to logistics and malfunctioning, missing values are often “seen” in GPS time series, especially in polar regions. Acquiring a consistent and complete time series is the prerequisite for accurate and reliable statical analysis. Previous imputation studies focused on the temporal relationship of time series, and only a few studies used spatial relationships and/or were based on machine learning methods. In this study, we impute 20 Greenland GPS time series using missForest, which is a new machine learning method for data imputation. The imputation performance of missForest and that of four traditional methods are assessed, and the methods’ impacts on principal component analysis (PCA) are investigated. Results show that missForest can impute more than a 30-day gap, and its imputed time series has the least influence on PCA. When the gap size is 30 days, the mean absolute value of the imputed and true values for missForest is 2.71 mm. The normalized root mean squared error is 0.065, and the distance of the first principal component is 0.013. missForest outperforms the other compared methods. missForest can effectively restore the information of GPS time series and improve the results of related statistical processes, such as PCA analysis.

Independent Component Extraction from the Incomplete Coordinate Time Series of Regional GNSS Networks

Independent component analysis (ICA) is one of the most effective approaches in extracting independent signals from a global navigation satellite system (GNSS) regional station network. However, ICA requires the involved time series to be complete, thereby the missing data of incomplete time series should be interpolated beforehand. In this contribution, a modified ICA is proposed, by which the missing data are first recovered based on the reversible property between the original time series and decomposed principal components, then the complete time series are further processed with FastICA. To evaluate the performance of the modified ICA for extracting independent components, 24 regional GNSS network stations located in North China from 2011 to 2019 were selected. After the trend, annual and semiannual terms were removed from the GNSS time series, the first two independent components captured 17.42, 18.44 and 17.38% of the total energy for the North, East and Up coordinate components, more than those derived by the iterative ICA that accounted for 16.21%, 17.72% and 16.93%, respectively. Therefore, modified ICA can extract more independent signals than iterative ICA. Subsequently, selecting the 7 stations with less missing data from the network, we repeatedly process the time series after randomly deleting parts of the data and compute the root mean square error (RMSE) from the differences of reconstructed signals before and after deleting data. All RMSEs of modified ICA are smaller than those of iterative ICA, indicating that modified ICA can extract more exact signals than iterative ICA.

GNSS-TS-NRS: An Open-Source MATLAB-Based GNSS Time Series Noise Reduction Software

The global navigation satellite system (GNSS) has seen tremendous advances in measurement precision and accuracy, and it allows researchers to perform geodynamics and geophysics studies through the analysis of GNSS time series. Moreover, GNSS time series not only contain geophysical signals, but also unmodeled errors and other nuisance parameters, which affect the performance in the estimation of site coordinates and related parameters. As the number of globally distributed GNSS reference stations increases, GNSS time series analysis software should be developed with more flexible format support, better human–machine interaction, and with powerful noise reduction analysis. To meet this requirement, a new software named GNSS time series noise reduction software (GNSS-TS-NRS) was written in MATLAB and was developed. GNSS-TS-NRS allows users to perform noise reduction analysis and spatial filtering on common mode errors and to visualize GNSS position time series. The functions’ related theoretical background of GNSS-TS-NRS were introduced. Firstly, we showed the theoretical background algorithms of the noise reduction analysis (empirical mode decomposition (EMD), ensemble empirical mode decomposition (EEMD)). We also developed three improved algorithms based on EMD for noise reduction, and the results of the test example showed our proposed methods with better effect. Secondly, the spatial filtering model supported five algorithms on a separate common model error: The stacking filter method, weighted stacking filter method, correlation weighted superposition filtering method, distance weighted filtering method, and principal component analysis, as well as with batch processing. Finally, the developed software also enabled other functions, including outlier detection, correlation coefficient calculation, spectrum analysis, and distribution estimation. The main goal of the manuscript is to share the software with the scientific community to introduce new users to the GNSS time series noise reduction and application.

Impacts of Local Effects and Surface Loads on the Common Mode Error Filtering in Continuous GPS Measurements in the Northwest of Yunnan Province, China

While seasonal hydrological mass loading, derived from Gravity Recovery and Climate Experiment (GRACE) measurements, shows coherent spatial patterns and is an important source for the common mode error (CME) in continuous global positioning system (cGPS) measurements in Yunnan, it is a challenge to quantify local effects and detailed changes in daily GPS measurements by using GRACE data due to its low time and spatial resolutions. In this study, we computed and compared two groups of CMEs for nine cGPS sites in the northwest Yunnan province; rCMEs were computed with the residual cGPS time series having high inter-station correlations, while oCMEs were computed with all the GPS time series. The rCMEs-filtered time series had smaller variances and larger root mean square (RMS) reductions than those that were oCMEs-filtered, and when the stations local effects were not removed, spurious transient-like signals occurred. Compared with hydrological mass loading (HYDL), its combination with non-tidal atmosphere pressure and ocean mass reached a better agreement with the CME in the vertical component, with the Nash–Sutcliffe efficiency (NSE) increasing from 0.28 to 0.55 and the RMS reduction increasing from 15.19% to 33.4%, respectively. Our results suggest that it is necessary to evaluate the inter-station correlation and remove the possible noisy stations before conducting CME filtering, and that one should carefully choose surface loading models to correct the raw cGPS time series if CME filtering is not conducted.

Boundary-Included Enhanced Water Storage Changes Inferred by GPS in the Pacific Rim of the Western United States

We developed a new boundary-included inversion model to improve the terrestrial water storage (TWS) inverted from regional GPS vertical deformation data. Through defining a new disc load empirical function (DLEF) and considering the mass change effect from the near but outside region, the result shows the TWS is more reasonable than the one inverted directly. Six simulation tests further confirmed the effectiveness of the boundary-included model. Finally, our new boundary-included model was used to derive the TWS in the Pacific Rim of the western United States based on the GPS-observed vertical deformation information. The inversion results show that our boundary-included inversion model can effectively improve the inversion results by 10–20% in terms of variance reduction in the boundary regions.

Modeling and Forecasting the GPS Zenith Troposphere Delay in West Antarctica Based on Different Blind Source Separation Methods and Deep Learning

Tropospheric delay is an important error source in global positioning systems (GPS), and the water vapor retrieved from the tropospheric delay is widely used in meteorological research such as climate analysis and weather forecasting. Most zenith tropospheric delay (ZTD) models are presently used as positioning corrections, and few models are used for the estimation of water vapor, especially in Antarctica. Through two blind source separation algorithms (principal component analysis (PCA) and independent component analysis (ICA)), a back-propagation (BP) neural network and a deep learning technique (long short-term memory (LSTM) network), we establish an hourly high-accuracy ZTD model for GPS meteorology using the GPS-ZTD from 52 GPS stations in West Antarctica. Our results show that under the condition in which the principal components (PCs) and independent components (ICs) remain fixed after decomposition, the mean accuracy of the models for West Antarctica using PCA or ICA are better than 10 mm. Compared with the ZTDs from the nonmodeling stations, the mean root mean square (RMS) of the PCA and ICA models are 9.3 and 8.9 mm, respectively, and the correlation coefficients between the GPS-ZTD and model-ZTDs all exceed 90%. The accuracy of the ICA model is slightly higher than that of the PCA model, and the ICs of the ICA model show more consistent spatial responses. The six-hour forecast is the best among the forecast results, with a mean correlation coefficient of 90.6% and a mean RMS of 7.2 mm using GPS-ZTD. The long-term forecast result is significantly inaccurate, as the correlation coefficient between the 24-h forecast and GPS-ZTD is only 63.2%. Generally modest results have been achieved (HSS ≤ 0.38). Furthermore, the forecast accuracy in coastal areas is lower than that in inland areas. Our study confirms that the combined use of ICA and deep learning in ZTD modeling can effectively restore the original signals, and short-term forecasting can be effectively used in GPS meteorology. However, further development of the technology is necessary.

Extracting Common Mode Errors of Regional GNSS Position Time Series in the Presence of Missing Data by Variational Bayesian Principal Component Analysis

Removal of the common mode error (CME) is very important for the investigation of global navigation satellite systems’ (GNSS) error and the estimation of an accurate GNSS velocity field for geodynamic applications. The commonly used spatiotemporal filtering methods normally process the evenly spaced time series without missing data. In this article, we present the variational Bayesian principal component analysis (VBPCA) to estimate and extract CME from the incomplete GNSS position time series. The VBPCA method can naturally handle missing data in the Bayesian framework and utilizes the variational expectation-maximization iterative algorithm to search each principal subspace. Moreover, it could automatically select the optimal number of principal components for data reconstruction and avoid the overfitting problem. To evaluate the performance of the VBPCA algorithm for extracting CME, 44 continuous GNSS stations located in Southern California were selected. Compared to previous approaches, VBPCA could achieve better performance with lower CME relative errors when more missing data exists. Since the first principal component (PC) extracted by VBPCA is remarkably larger than the other components, and its corresponding spatial response presents nearly uniform distribution, we only use the first PC and its eigenvector to reconstruct the CME for each station. After filtering out CME, the interstation correlation coefficients are significantly reduced from 0.43, 0.46, and 0.38 to 0.11, 0.10, and 0.08, for the north, east, and up (NEU) components, respectively. The root mean square (RMS) values of the residual time series and the colored noise amplitudes for the NEU components are also greatly suppressed, with average reductions of 27.11%, 28.15%, and 23.28% for the former, and 49.90%, 54.56%, and 49.75% for the latter. Moreover, the velocity estimates are more reliable and precise after removing CME, with average uncertainty reductions of 51.95%, 57.31%, and 49.92% for the NEU components, respectively. All these results indicate that the VBPCA method is an alternative and efficient way to extract CME from regional GNSS position time series in the presence of missing data. Further work is still required to consider the effect of formal errors on the CME extraction during the VBPCA implementation.

Seismo-Deformation Anomalies Associated with the M6.1 Ludian Earthquake on August 3, 2014

A time-frequency method retrieving the acceleration changes in the terminal stage of theM6.1 Ludian earthquake in China is discussed in this article. The non-linear, non-stationaryseismo-demformation was obtained by using the Hilbert–Huang transform and followed by aband-pass filter. We found that the temporal evolution of the residual GNSS-derived orientationexhibits a unique disorder-alignment-disorder sequence days before the earthquake whichcorresponds well with the four stages of an earthquake: elastic strain buildup, crack developments,deformation, and the terminal stage of material failure. The disordering orientations are graduallyaligned with a common direction a few days before the terminal stage. This common direction isconsistent with the most compressive axis derived from the seismological method. In addition, theregion of the stress accumulation, as identified by the size of the disordered orientation, isgenerally consistent with the earthquake preparation zones estimated by using numerical models.

Analysis of the Potential Contributors to Common Mode Error in Chuandian Region of China

Common mode error (CME) in Chuandian region of China is derived from 6-year continuous GPS time series and is identified by principal component analysis (PCA) method. It is revealed that the temporal behavior of the CME is not purely random, and contains unmodeled signals such as nonseasonal mass loadings. Its spatial distribution is quite uniform for all GPS sites in the region, and the first principal component, uniformly distributed in the region, has a spatial response of more than 70%. To further explore the potential contributors of CME, daily atmospheric mass loading and soil moisture mass loading effects are evaluated. Our results show that ~15% of CME can be explained by these daily surface mass loadings. The power spectral analysis is used to assess the CME. After removing atmospheric and soil moisture loadings from the CME, the power of the CME reduces in a wide range of frequencies. We also investigate the contribution of CME in GPS filtered residuals time series and it shows the Root Mean Squares (RMSs) of GPS time series are reduced by applying of the mass loading corrections in CME. These comparison results demonstrate that daily atmosphere pressure and the soil moisture mass loadings are a part of contributors to the CME in Chuandian region of China.

Forecasting of landslide displacements using a chaos theory based wavelet analysis-Volterra filter model

Landslide displacement time series can directly reflects landslide deformation and stability characteristics. Hence, forecasting of the non-linear and non-stationary displacement time series is necessary and significant for early warning of landslide failure. Traditionally, conventional machine learning methods are adopted as forecasting models, these forecasting models mainly determine the input and output variables experientially and does not address the non-stationary characteristics of displacement time series. However, it is difficult for these conventional machine learning methods to obtain appropriate input-output variables, to determine appropriate model parameters and to acquire satisfied prediction performance. To deal with these drawbacks, this study proposes the wavelet analysis (WA) to decompose the displacement time series into low- and high-frequency components to address the non-stationary characteristics; then proposes thee chaos theory to obtain appropriate input-output variables of forecasting models, and finally proposes Volterra filter model to construct the forecasting model. The GPS monitoring cumulative displacement time series, recorded on the Shuping and Baijiabao landslides, distance measuring equipment monitoring displacements on the Xintan landslide in Three Gorges Reservoir area of China, are used as test data of the proposed chaotic WA-Volterra model. The chaotic WA-support vector machine (SVM) model and single chaotic Volterra model without WA method, are used as comparisons. The results show that there are chaos characteristics in the GPS monitoring displacement time series, the non-stationary characteristics of landslide displacements are captured well by the WA method, and the model input-output variables are selected suitably using chaos theory. Furthermore, the chaotic WA-Volterra model has higher prediction accuracy than the chaotic WA-SVM and single chaotic Volterra models.

Transient Deformation in California From Two Decades of GPS Displacements: Implications for a Three-Dimensional Kinematic Reference Frame

Our understanding of plate boundary deformation has been enhanced by transient signals observed against the backdrop of time-independent secular motions. We make use of a new analysis of displacement time series from about 1,000 continuous Global Positioning System (GPS) stations in California from 1999 to 2018 to distinguish tectonic and nontectonic transients from secular motion. A primary objective is to define a high-resolution three-dimensional reference frame (datum) for California that can be rapidly maintained with geodetic data to accommodate both secular and time-dependent motions. To this end, we compare the displacements to those predicted by a horizontal secular fault slip model for the region and construct displacement and strain rate fields. Over the past 19 years, California has experienced 19 geodetically detectable earthquakes and widespread postseismic deformation. We observe postseismic strain rate variations as large as 1,000 nstrain/year with moment releases equivalent up to an Mw6.8 earthquake. We find significant secular differences up to 10 mm/year with the fault slip model, from the Mendocino Triple Junction to the southern Cascadia subduction zone, the northern Basin and Range, and the Santa Barbara channel. Secular vertical uplift is observed across the Transverse Ranges, Coastal Ranges, Sierra Nevada, as well as large-scale postseismic uplift after the 1999 Mw7.1 Hector Mine and 2010 Mw7.2 El Mayor-Cucapah earthquakes. We also identify areas of vertical land motions due to anthropogenic, natural, and magmatic processes. Finally, we demonstrate the utility of the kinematic datum by improving the accuracy of high-spatial-resolution 12-day repeat-cycle Sentinel-1 Interferometric Synthetic Aperture Radar displacement and velocity maps.

Common Mode Component and Its Potential Effect on GPS-Inferred Three-Dimensional Crustal Deformations in the Eastern Tibetan Plateau

Surface and deep potential geophysical signals respond to the spatial redistribution of global mass variations, which may be monitored by geodetic observations. In this study, we analyze dense Global Positioning System (GPS) time series in the Eastern Tibetan Plateau using principal component analysis (PCA) and wavelet time-frequency spectra. The oscillations of interannual and residual signals are clearly identified in the common mode component (CMC) decomposed from the dense GPS time series from 2000 to 2018. The newly developed spherical harmonic coefficients of the Gravity Recovery and Climate Experiment Release-06 (GRACE RL06) are adopted to estimate the seasonal and interannual patterns in this region, revealing hydrologic and atmospheric/nontidal ocean loads. We stack the averaged elastic GRACE-derived loading displacements to identify the potential physical significance of the CMC in the GPS time series. Interannual nonlinear signals with a period of ~3 to ~4 years in the CMC (the scaled principal components from PC1 to PC3) are found to be predominantly related to hydrologic loading displacements, which respond to signals (El Niño/La Niña) of global climate change. We find an obvious signal with a period of ~6 yr on the vertical component that could be caused by mantle-inner core gravity coupling. Moreover, we evaluate the CMC’s effect on the GPS-derived velocities and confirm that removing the CMC can improve the recognition of nontectonic crustal deformation, especially on the vertical component. Furthermore, the effects of the CMC on the three-dimensional velocity and uncertainty are presented to reveal the significant crustal deformation and dynamic processes of the Eastern Tibetan Plateau.

Akaike's Bayesian Information Criterion for the Joint Inversion of Terrestrial Water Storage Using GPS Vertical Displacements, GRACE and GLDAS in Southwest China

Global navigation satellite systems (GNSS) techniques, such as GPS, can be used to accurately record vertical crustal movements induced by seasonal terrestrial water storage (TWS) variations. Conversely, the TWS data could be inverted from GPS-observed vertical displacement based on the well-known elastic loading theory through the Tikhonov regularization (TR) or the Helmert variance component estimation (HVCE). To complement a potential non-uniform spatial distribution of GPS sites and to improve the quality of inversion procedure, herein we proposed in this study a novel approach for the TWS inversion by jointly supplementing GPS vertical crustal displacements with minimum usage of external TWS-derived displacements serving as pseudo GPS sites, such as from satellite gravimetry (e.g., Gravity Recovery and Climate Experiment, GRACE) or from hydrological models (e.g., Global Land Data Assimilation System, GLDAS), to constrain the inversion. In addition, Akaike’s Bayesian Information Criterion (ABIC) was employed during the inversion, while comparing with TR and HVCE to demonstrate the feasibility of our approach. Despite the deterioration of the model fitness, our results revealed that the introduction of GRACE or GLDAS data as constraints during the joint inversion effectively reduced the uncertainty and bias by 42% and 41% on average, respectively, with significant improvements in the spatial boundary of our study area. In general, the ABIC with GRACE or GLDAS data constraints displayed an optimal performance in terms of model fitness and inversion performance, compared to those of other GPS-inferred TWS methodologies reported in published studies.

An Improved GPS-Inferred Seasonal Terrestrial Water Storage Using Terrain-Corrected Vertical Crustal Displacements Constrained by GRACE

Based on a geophysical model for elastic loading, the application potential of Global Positioning System (GPS) vertical crustal displacements for inverting terrestrial water storage has been demonstrated using the Tikhonov regularization and the Helmert variance component estimation since 2014. However, the GPS-inferred terrestrial water storage has larger resulting amplitudes than those inferred from satellite gravimetry (i.e., Gravity Recovery and Climate Experiment (GRACE)) and those simulated from hydrological models (e.g., Global Land Data Assimilation System (GLDAS)). We speculate that the enlarged amplitudes should be partly due to irregularly distributed GPS stations and the neglect of the terrain effect. Within southwest China, covering part of southeastern Tibet as a study region, a novel GPS-inferred terrestrial water storage approach is proposed via terrain-corrected GPS and supplementary vertical crustal displacements inferred from GRACE, serving as "virtual GPS stations" for constraining the inversion. Compared to the Tikhonov regularization and Helmert variance component estimation, we employ Akaike’s Bayesian Information Criterion as an inverse method to prove the effectiveness of our solution. Our results indicate that the combined application of the terrain-corrected GPS vertical crustal displacements and supplementary GRACE spatial data constraints improves the inversion accuracy of the GPS-inferred terrestrial water storage from the Helmert variance component estimation, Tikhonov regularization, and Akaike’s Bayesian Information Criterion, by 55%, 33%, and 41%, respectively, when compared to that of the GLDAS-modeled terrestrial water storage. The solution inverted with Akaike’s Bayesian Information Criterion exhibits more stability regardless of the constraint conditions, when compared to those of other inferred solutions. The best Akaike’s Bayesian Information Criterion inverted solution agrees well with the GLDAS-modeled one, with a root-mean-square error (RMSE) of 3.75 cm, equivalent to a 15.6% relative error, when compared to 39.4% obtained in previous studies. The remaining discrepancy might be due to the difference between GPS and GRACE in sensing different surface water storage components, the remaining effect of the water storage changes in rivers and reservoirs, and the internal error in the geophysical model for elastic loading.

A Sub-Regional Extraction Method of Common Mode Components from IGS and CMONOC Stations in China

There is always a need to extract more accurate regional common mode component (CMC) series from coordinate time series of Global Positioning System (GPS) stations, which would be of great benefit to describe the deformation features of the Earth’s surface with more reliability. For this purpose, this paper combines all 11 International Global Navigation Satellite System (GNSS) Service (IGS) stations in China with over 70 stations selected from the Crustal Movement Observation Network of China (CMONOC) to compute CMC series of IGS stations by using a principal component analysis (PCA) method under cases of one whole region and eight sub-regions. The comparison results show that the percentage of first-order principal component (PC1) in North, East and Up components increase by 10.8%, 16.1% and 25.1%, respectively, after dividing the whole China region into eight sub-regions. Meanwhile, Root Mean Square (RMS) reduction rates of residual series that have removed CMC also improve obviously after partitioning. In addition, we compute displacements of these IGS stations caused by environmental loadings (including atmospheric pressure loading, non-tidal oceanic loading and hydrological loading) to analyze their contributions to the non-linear variation in GPS coordinate time series. The comparison result shows that the method we raise, PCA filtering in sub-regions, performs better than the environmental loading corrections (ELCs) in improving the signal-to-noise ratio (SNR) of GPS coordinate time series. This paper raises new criteria for selecting appropriate CMONOC stations around IGS stations when computing sub-regional CMC, involving three criteria of interstation distance, geology and self-condition of stations themselves. According to experiments, these criteria are implemental and effective in selecting suitable stations, by which to extract sub-regional CMC with higher accuracy.

Spatiotemporal Filtering and Noise Analysis for Regional GNSS Network in Antarctica Using Independent Component Analysis

The common mode error (CME) and optimal noise model are the two most important factors affecting the accuracy of time series in regional Global Navigation Satellite System (GNSS) networks. Removing the CME and selecting the optimal noise model can effectively improve the accuracy of GNSS coordinate time series. The CME, a major source of error, is related to the spatiotemporal distribution; hence, its detrimental effects on time series can be effectively reduced through spatial filtering. Independent component analysis (ICA) is used to filter the time series recorded by 79 GPS stations in Antarctica from 2010 to 2018. After removing stations exhibiting strong local effects using their spatial responses, the filtering results of residual time series derived from principal component analysis (PCA) and ICA are compared and analyzed. The Akaike information criterion (AIC) is then used to determine the optimal noise model of the GPS time series before and after ICA/PCA filtering. The results show that ICA is superior to PCA regarding both the filter results and the consistency of the optimal noise model. In terms of the filtering results, ICA can extract multisource error signals. After ICA filtering, the root mean square (RMS) values of the residual time series are reduced by 14.45%, 8.97%, and 13.27% in the east (E), north (N), and vertical (U) components, respectively, and the associated speed uncertainties are reduced by 13.50%, 8.06% and 11.82%, respectively. Furthermore, different GNSS time series in Antarctica have different optimal noise models with different noise characteristics in different components. The main noise models are the white noise plus flicker noise (WN+FN) and white noise plus power law noise (WN+PN) models. Additionally, the spectrum index of most PN is close to that of FN. Finally, there are more stations with consistent optimal noise models after ICA filtering than there are after PCA filtering.

Effects of Spatiotemporal Filtering on the Periodic Signals and Noise in the GPS Position Time Series of the Crustal Movement Observation Network of China

Analysis of Global Positioning System (GPS) position time series and its common mode components (CMC) is very important for the investigation of GPS technique error, the evaluation of environmental loading effects, and the estimation of a realistic and unbiased GPS velocity field for geodynamic applications. In this paper, we homogeneously processed the daily observations of 231 Crustal Movement Observation Network of China (CMONOC) Continuous GPS stations to obtain their position time series. Then, we filtered out the CMC and evaluated its effects on the periodic signals and noise for the CMONOC time series. Results show that, with CMC filtering, peaks in the stacked power spectra can be reduced at draconitic harmonics up to the 14th, supporting the point that the draconitic signal is spatially correlated. With the colored noise suppressed by CMC filtering, the velocity uncertainty estimates for both of the two subnetworks, CMONOC-I (≈16.5 years) and CMONOC-II (≈4.6 years), are reduced significantly. However, the CMONOC-II stations obtain greater reduction ratios in velocity uncertainty estimates with average values of 33%, 38%, and 54% for the north, east, and up components. These results indicate that CMC filtering can suppress the colored noise amplitudes and improve the precision of velocity estimates. Therefore, a unified, realistic, and three-dimensional CMONOC GPS velocity field estimated with the consideration of colored noise is given. Furthermore, contributions of environmental loading to the vertical CMC are also investigated and discussed. We find that the vertical CMC are reduced at 224 of the 231 CMONOC stations and 170 of them are with a root mean square (RMS) reduction ratio of CMC larger than 10%, confirming that environmental loading is one of the sources of CMC for the CMONOC height time series.

Tracking the weight of Hurricane Harvey's stormwater using GPS data

On 26 August 2017, Hurricane Harvey struck the Gulf Coast as a category four cyclone depositing ~95 km³ of water, making it the wettest cyclone in U.S.

HISTORY

Water left in Harvey's wake should cause elastic loading and subsidence of Earth's crust, and uplift as it drains into the ocean and evaporates. To track daily changes of transient water storage, we use Global Positioning System (GPS) measurements, finding a clear migration of subsidence (up to 21 mm) and horizontal motion (up to 4 mm) across the Gulf Coast, followed by gradual uplift over a 5-week period. Inversion of these data shows that a third of Harvey's total stormwater was captured on land (25.7 ± 3.0 km³), indicating that the rest drained rapidly into the ocean at a rate of 8.2 km³/day, with the remaining stored water gradually lost over the following 5 weeks at ~1 km³/day, primarily by evapotranspiration. These results indicate that GPS networks can remotely track the spatial extent and daily evolution of terrestrial water storage following transient, extreme precipitation events, with implications for improving operational flood forecasts and understanding the response of drainage systems to large influxes of water.

Temporal-Spatial Surface Seasonal Mass Changes and Vertical Crustal Deformation in South China Block from GPS and GRACE Measurements

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.

Physical applications of GPS geodesy: a review

Geodesy, the oldest science, has become an important discipline in the geosciences, in large part by enhancing Global Positioning System (GPS) capabilities over the last 35 years well beyond the satellite constellation's original design. The ability of GPS geodesy to estimate 3D positions with millimeter-level precision with respect to a global terrestrial reference frame has contributed to significant advances in geophysics, seismology, atmospheric science, hydrology, and natural hazard science. Monitoring the changes in the positions or trajectories of GPS instruments on the Earth's land and water surfaces, in the atmosphere, or in space, is important for both theory and applications, from an improved understanding of tectonic and magmatic processes to developing systems for mitigating the impact of natural hazards on society and the environment. Besides accurate positioning, all disturbances in the propagation of the transmitted GPS radio signals from satellite to receiver are mined for information, from troposphere and ionosphere delays for weather, climate, and natural hazard applications, to disturbances in the signals due to multipath reflections from the solid ground, water, and ice for environmental applications. We review the relevant concepts of geodetic theory, data analysis, and physical modeling for a myriad of processes at multiple spatial and temporal scales, and discuss the extensive global infrastructure that has been built to support GPS geodesy consisting of thousands of continuously operating stations. We also discuss the integration of heterogeneous and complementary data sets from geodesy, seismology, and geology, focusing on crustal deformation applications and early warning systems for natural hazards.

Seasonal Mass Changes and Crustal Vertical Deformations Constrained by GPS and GRACE in Northeastern Tibet

Surface vertical deformation includes the Earth’s elastic response to mass loading on or near the surface. Continuous Global Positioning System (CGPS) stations record such deformations to estimate seasonal and secular mass changes. We used 41 CGPS stations to construct a time series of coordinate changes, which are decomposed by empirical orthogonal functions (EOFs), in northeastern Tibet. The first common mode shows clear seasonal changes, indicating seasonal surface mass re-distribution around northeastern Tibet. The GPS-derived result is then assessed in terms of the mass changes observed in northeastern Tibet. The GPS-derived common mode vertical change and the stacked Gravity Recovery and Climate Experiment (GRACE) mass change are consistent, suggesting that the seasonal surface mass variation is caused by changes in the hydrological, atmospheric and non-tidal ocean loads. The annual peak-to-peak surface mass changes derived from GPS and GRACE results show seasonal oscillations in mass loads, and the corresponding amplitudes are between 3 and 35 mm/year. There is an apparent gradually increasing gravity between 0.1 and 0.9 μGal/year in northeast Tibet. Crustal vertical deformation is determined after eliminating the surface load effects from GRACE, without considering Glacial Isostatic Adjustment (GIA) contribution. It reveals crustal uplift around northeastern Tibet from the corrected GPS vertical velocity. The unusual uplift of the Longmen Shan fault indicates tectonically sophisticated processes in northeastern Tibet.

Coastal structure, sea-level changes and vertical motion of the land in the Mediterranean

The Mediterranean basin is an important area of the Earth for studying the interplay between geodynamic processes and landscape evolution affected by tectonic, glacio-hydro-isostatic and eustatic factors. We focus on determining vertical deformations and relative sea-level change of the coastal zone utilizing geological, archaeological, historical and instrumental data, and modelling. For deformation determinations on recent decadal to centennial time scales, seismic strain analysis based on about 6000 focal mechanisms, surface deformation analysis based on some 850 continuous GPS stations, and 57 tide gauge records were used. Utilizing data from tectonically stable areas, reference surfaces were established to separate tectonic and climate (eustatic) signals throughout the basin for the last 20 000 years. Predominant Holocene subsidence (west coast of Italy, northern Adriatic sea, most of Greece and Turkey are areas at risk of flooding owing to relative sea-level rise), uplift (local areas in southwestern Italy and southern Greece) or stability (northwestern and central western Mediterranean and Levant area) were determined. Superimposed on the long trends, the coasts are also impacted by sudden extreme events such as recurring large storms and numerous, but unpredictable tsunamis caused by the high seismicity of parts of the basins.

Supplementary material:

A table of locations and timings of the largest tsunamis in the Mediterranean during the last 5660 years BP is available at http://www.geolsoc.org.uk/SUP18757.