Next Article in Journal
Evaluation of GPM-Era Satellite Precipitation Products on the Southern Slopes of the Central Himalayas Against Rain Gauge Data
Next Article in Special Issue
The Use of National CORS Networks for Determining Temporal Mass Variations within the Earth’s System and for Improving GRACE/GRACE-FO Solutions
Previous Article in Journal
Evaluation of Six High-Spatial Resolution Clear-Sky Surface Upward Longwave Radiation Estimation Methods with MODIS
Previous Article in Special Issue
Constrained Linear Deconvolution of GRACE Anomalies to Correct Spatial Leakage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Surface Mass Variations from GPS and GRACE/GFO: A Case Study in Southwest China

1
School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China
2
Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, Wuhan University, Wuhan 430079, China
3
Center for Space Research, University of Texas at Austin, Austin, TX 78759, USA
4
MOE Key Laboratory of Fundamental Physical Quantities Measurement, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(11), 1835; https://doi.org/10.3390/rs12111835
Submission received: 25 April 2020 / Revised: 4 June 2020 / Accepted: 4 June 2020 / Published: 5 June 2020
(This article belongs to the Special Issue Terrestrial Hydrology Using GRACE and GRACE-FO)

Abstract

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

1. Introduction

Surface mass variations mainly include terrestrial water storage (TWS) change, melting of polar ice sheets and mountain glaciers, atmospheric pressure change, ocean mass change, and mass changes associated with solid Earth geophysical phenomena. At present, large scale surface mass variations are primarily quantified by using the Gravity Recovery and Climate Experiment (GRACE) time-variable gravity field solutions and outputs from hydrological, atmospheric, and ocean models [1,2,3,4]. Even though GRACE satellite gravimetry has demonstrated the great potentials of measuring large-scale mass changes associated with land hydrological cycle, glacial melting and sea level changes, earthquakes, and human activities with unprecedented accuracy, it still has some major limitations: (1) The truncation of GRACE Spherical Harmonic (SH) coefficients and needed spatial filtering and smoothing to suppress the longitudinal stripping noises will introduce leakage bias and degrade the spatial resolution of GRACE-determined surface mass variations, which is limited to ~300–500 km [5,6,7]; (2) the uncertainty of GRACE low-degree SH coefficients is relatively larger, especially the degree-2 zonal term C20, which needs to be replaced with independent estimates from Satellite Laser Ranging (SLR) [8,9]; (3) the data quality of GRACE observations degrades after August 2016, which affects the accuracy of GRACE-inferred mass changes; (4) there is an approximately one-year gap between the GRACE and GRACE Follow-On (GFO) missions, which affects the interpretation and applications of GRACE/GFO data.
Surface deformation due to loading caused by the redistribution of Earth’s surface mass can be continuously monitored by the ground-based Global Positioning System (GPS) observation network. Since the surface loading inversion approach was proposed by Blewitt et al. [10], mass variations inverted from GPS surface deformation data have been serving as a complement to the GRACE/GFO observations both at global and regional scales [11,12,13,14,15,16,17]. Compared with GRACE/GFO, GPS has its own advantages: (1) The GPS network is widely distributed globally, and the observation data can be obtained in near real time at millimeter-level accuracy with higher spatio–temporal resolution; (2) the GPS surface deformation, especially the vertical displacement, is more sensitive to changes of local (~10 km) and regional mass loading, and the spatial resolution of mass changes inverted from a dense GPS network can reach ~50 km [16,17]; (3) due to the continuous and real-time monitoring of the surface deformation by GPS networks, GPS-derived mass variations can help bridge the gap between the GRACE and GFO missions [15,16].
Because Earth’s surface deformation due to surface mass loading can also be estimated by GRACE SH gravity solutions [13], numerous previous studies have focused on the consistency analysis of surface deformations observed by GPS and GRACE. Davis et al. [18] demonstrated the high correlation of annual variations of GPS-measured and GRACE-derived vertical displacement time series in the Amazon Basin. Since then, several other studies had shown good consistency between GPS-measured and GRACE-derived vertical displacements, like in the Nepal Himalayas [19], West Africa [20], Greenland [21], Southeastern Poland [22], the Southeastern Tibetan Plateau [23], Southwest China [24], and the North China Plain [25]. Moreover, Wahr et al. [26] found that the GPS vertical displacement component combined with the horizontal displacement component could be used to determine whether nearby loading was concentrated in a small region, and where that region was. However, due to the different measurement modes and error sources of GPS and GRACE, poor agreements between GPS and GRACE vertical displacement estimates have also been reported in Central American [27] and Europe [28]. Nevertheless, these previous studies suggested that surface deformations determined by GPS and GRACE are consistent in most of the studied regions, and the surface mass variations can be inverted by GPS deformation data when non-loading effects are carefully removed.
A number of previous studies have also demonstrated that the regional surface mass variations can be independently inverted by using continuous GPS surface deformation data. Argus et al. [16] used GPS vertical displacement data to invert seasonal variations of TWS, and found that the inversion results showed good consistency with those from GRACE and hydrological models. The study further showed that TWS change can be recovered by GPS vertical displacements at a spatial resolution of ~50 km over California. Fu et al. [17] inverted monthly surface mass change time series using GPS vertical displacement data in the Washington and Oregon regions, and showed that GPS-inferred water mass variations were highly correlated and consistent with those derived from GRACE and hydrological models, and could help fill in the gaps of GRACE. Jin and Zhang [29] estimated the mass changes using continuous GPS vertical displacement data in Southwestern USA, and the severe drought in 2012 was captured by the GPS-derived mass change time series, indicating that continuous GPS deformation observation has the potential to serve as a real-time drought indicator. Enzminger et al. [30] employed GPS vertical displacements to assess the accuracy of the snow water equivalent over Western USA Mountains. Zhang et al. [31] utilized GPS vertical displacements from 29 GPS stations to study mass variations at a spatial resolution of approximately 100 km over Yunnan Province in China. More recently, Ferreira et al. [32] presented a resolution analysis of daily mass variations inverted by GPS vertical displacement data over South America. Fok et al. [33] discussed a method for improving the inversion accuracy of GPS-derived mass variations in Southwestern China through combined applications of terrain-corrected GPS vertical displacements and supplementary GRACE spatial data constraints.
This study aims to further analyze the performance of regional mass variations inverted by GPS vertical displacement data under different density distributions of GPS stations, and examine the consistency of mass changes inferred from GPS vertical displacements, GRACE, and the latest GFO observations in the spatial and temporal domains, and demonstrate the feasibility of GPS-derived mass variations to fill in the gap between GRACE and GFO missions. To this end, continuous GPS vertical displacement time series from the Crustal Movement Observation Network of China (CMONOC), the latest GRACE/GFO RL06 SH solutions, and GRACE RL06 mascon solutions from the Center for Space Research (CSR) were used to investigate surface mass variations in four regions of Southwest China. The spatial distributions and seasonal characteristics of the mass change time series derived from GPS vertical displacements and GRACE/GFO were analyzed and compared in detail. The scale factor methods (both basin-scaled and pixel-scaled) used in Landerer et al. [34] were employed to reduce the amplitude discrepancies between GPS-derived and GRACE/GFO-inferred mass change time series. The rest of the paper is organized as the following. Overview of the studied area and the datasets used for inferring surface mass variations are described in Section 2. Data post-processing of GPS vertical displacements, and inversion model for surface mass variations are presented in Section 3. In Section 4, comparisons of different mass change inversion results and discussions of differences among these results are presented, and the potentials of using GPS-derived mass change time series to bridge the GRACE/GFO mission gap are also demonstrated. Finally, conclusions are drawn in Section 5.

2. Study Area and Datasets

2.1. Study Area

In this paper, Southwest China was chosen as the study area, which includes four regions or basins (Figure 1): The Yunnan Province (YNP), Min River Basin (MRB), Jialing River Basin (JLRB) and Wu River Basin (WRB). Table 1 lists the areas, number of GPS stations, and station density for each of the regions. Rivers and valleys are widely distributed over this studied area, and the landform mainly consists of plateaus and mountains with great relief. Furthermore, it belongs to subtropical monsoon climate zone, with uneven annual average temperature distribution and abundant rainfall. Hence, it is of great significance for the management and utilization of natural resources to investigate surface mass variations in Southwest China.

2.2. GPS Vertical Displacement Data

GPS vertical displacement data during the period of January 2011 to August 2019 from a total of 85 permanent GPS stations of CMONOC (ftp://ftp.cgps.ac.cn/products/position/gamit/) were used to invert the surface mass variations over Southwest China. All the permanent GPS stations of CMONOC are equipped with multi-frequency receivers with a sampling interval of 30 s, and the GPS station coordinates time series are processed by the GAMIT/GLOBK software (release 10.4) based on the double-differenced model [25,35]. Firstly, loosely constrained daily coordinates and covariance were processed using GAMIT software. In this procedure, satellite orbit parameters and satellite antenna offsets were adjusted using the reprocessed orbits from the International GNSS Service (IGS). Solid Earth tides and pole tides were corrected according to the International Earth Rotation Service (IERS) conventions 2003, and ocean tidal loading was corrected using the Finite Element Solutions 2004 (FES2004) model. Receiver antenna phase center was corrected using the absolute antenna phase center model [36]. Tropospheric delays were modeled using Global Mapping Function (GMF) [37] every hour, and the first order ionospheric delays were eliminated through the ionosphere-free linear combination. Atmospheric zenith delay and atmospheric gradients were estimated every hour and every day, respectively. Secondly, the loosely constrained daily solutions with station coordinates, covariance, atmosphere delays, and orbit parameters were combined with global solutions produced by the Scripps Orbital and Position Analysis Center (SOPAC, http://sopac.ucsd.edu/) using the GLOBK software. Finally, the datum transformation was performed to obtain a station coordinate time series with respect to the International Terrestrial Reference Frame (ITRF) 2008 using approximately 50 globally distributed IGS stations through seven-parameter transformation model, and then the daily GPS vertical displacement time series were obtained. More details of GPS data processing can be found from the manual of CMONOC reference stations (ftp://ftp.cgps.ac.cn/doc/processing_manual.pdf).

2.3. GRACE/GFO-Inferred Surface Mass Variations

GRACE and GFO monthly SH solutions (up to degree and order 60) were used to infer surface mass variations in accordance with the approach presented by Wahr et al. [38]. All of the SH solutions used in our study are published by the CSR at the University of Texas at Austin (UTCSR), including CSR RL06 GRACE SH solutions (hereinafter called CSR GRACE-SH) covering the period August 2002 to June 2017 (ftp://isdcftp.gfz-potsdam.de/grace/Level-2/CSR/RL06/) and the latest CSR RL06 GFO solutions (hereinafter called CSR GFO-SH) from June 2018 to August 2019 (ftp://isdcftp.gfz-potsdam.de/grace-fo/Level-2/CSR/RL06/). Since the coordinate origin of the GRACE/GFO data system is defined as the center of mass of the total Earth (CM) frame, the degree-1 SH terms are not included in GRACE/GFO SH solutions. However, GPS surface deformations are observed with respect to the center of surface figure (CF) frame, and the degree-1 SH terms reflect the motion of the geocenter with respect to the CM frame. In order to keep the consistency of the reference frame between GPS and GRACE/GFO, the degree-1 SH coefficients had been added back to the GRACE/GFO SH solutions using the degree-1 components provided by the GRACE project as supplementary datasets (GRACE Technical Note 13) [39]. Meanwhile, due to the uncertainty of low-degree SH coefficients from GRACE/GFO data, the C20 and C30 coefficients had been replaced by SLR solutions provided by the GRACE/GFO Science Data System (SDS) team [40].
Due to SH coefficients at high degrees and orders being contaminated by the longitudinal stripes and other errors, the de-striping method P3M6 [41] and Gaussian filtering with the averaging radius of 300 km were used to suppress to spatial noises. But the above post-processing of SH solutions will reduce the spatial resolution and lead to signal leakage. In this study, the leakage errors were corrected by the basin scale factor method [34], and the scaling factors were derived from the Global Land Data Assimilation System (GLDAS) Noah outputs [42] by calculating the ratio between the original TWS and the truncated and filtered TWS in each studied region. Moreover, in order to be consistent with the inversion results from GPS vertical displacement data, the atmosphere combined with ocean (GAC) fields from Atmosphere and Ocean De-aliasing Level-1B (AOD1B) RL06 products (ftp://isdcftp.gfz-potsdam.de/grace/Level-2/CSR/RL06/) had been added back to the SH solutions. By doing this, both of the inversion results from GPS and GRACE/GFO contained the mass changes from TWS, non-tidal atmosphere, and ocean.
Due to the lack of adequate in situ measurements to validate GRACE/GFO observations, it is challenging to assess uncertainty of GRACE/GFO-inferred mass changes. Generally, the Root Mean Square (RMS) of mass variations over an open ocean area covering the same latitude range of the four studied regions can be approximately regarded as the uncertainty in GRACE/GFO mass change estimates [43,44]. In this study, an ocean area in the Pacific [20°N–35°N,170°E–220°E] was picked up to estimate the uncertainty of GRACE/GFO-inferred mass variations over Southwest China. It is approximately at the same latitudes as Southwest China and far away from land to avoid leakage of mass variation from land hydrology. It should be noted that this uncertainty assessment does not account for leakage errors in the four studied regions, which would lead to underestimate the uncertainty of GRACE/GFO mass change estimates.

2.4. GRACE Mascon Solutions

GRACE mascon solutions are typically presented as equivalent surface mass variations, and the longitudinal stripes are suppressed by applying regularization constraints [45,46,47]. Therefore, GRACE mascon solutions do not require de-striping or spatial filtering like in the SH solutions, and they can offer improved spatial resolution for analysis of regional surface mass variations. For these reasons, CSR GRACE RL06 mascon solutions (hereinafter called CSR GRACE-M, downloaded at http://www2.csr.utexas.edu/grace/RL06_mascons.html) were used to compare with the inversion results obtained from GPS vertical displacements and GRACE SH solutions.
The CSR mascon solutions were estimated on an equal area geodesic grid consisting of hexagonal tiles approximately 120 km wide (~1° × 1° at the equator), and mass variations in each mascon cell were computed from GRACE range-rate observations via partials with respect to SH expansion. Constraints on the solution include application of a time-variable regularization matrix from 200 km Gaussian smoothed regularized SH solutions, separation of land and ocean signals to reduce leakage, and forward modeling of annual and trend signals from ice losses in polar regions and mountain glaciers to avoid over-constraining the mascon estimation [47,48]. The CSR mascon solutions were provided on 0.5 × 0.5° grids, but the actual spatial resolution was significantly lower than that of grid resolution. Therefore, direct estimates of mass changes from CSR mascon solutions were also subjective to potential leakage bias [3,49]. Like the GRACE SH solutions, the leakage corrections for CSR mascon solutions were also reduced by the basin scale factor method, but a 200 km Gaussian filtering was applied to calculate the scaling factors based on the GLDAS Noah outputs. The uncertainty assessment of CSR mascon solutions over Southwest China is the same as that mentioned in Section 2.3.

3. Surface Mass Variations Inverted by GPS Vertical Displacements

3.1. Post-Processing of GPS Vertical Displacements

Because the effects of solid Earth tides, ocean tides loadings, and pole tides have been removed from the GPS vertical displacement time series of CMONOC, the GPS vertical displacements mainly reflect surface deformation caused by mass changes from TWS, non-tidal atmosphere, and ocean. However, the draconitic errors (i.e., ~352 days [50]) that are contained in the GPS vertical displacements must be removed before estimating surface mass variations. Unlike Glaser et al. [51], who used the combination of GPS and SLR data to remove the draconitic errors, we followed the method used in Fok et al. [33] to reduce the influences of draconitic errors on the desired seasonal variation (annual and semi-annual) signals by applying spectral filtering in the frequency domain. In addition, due to the influences of anthropogenic activities and natural environments, there would be data jumps and discontinuities in the GPS vertical displacement time series, e.g., the HBZG station shown in Figure 2. It is essential to carry out outlier detection and data interpolation for these vertical displacement time series. Here, the gross error detection method from Nikolaidis et al. [52] was used to remove the outliers, and the cubic spline interpolation and Principal Component Analysis (PCA) methods [53] were used to complement the missing data. Then the GPS vertical displacement time series of each station was fitted with a position, velocity, and sinusoids with semi-annual and annual periods. Because the elastic long-term trend (due to glacial isostatic adjustment or tectonic effects) can hardly be separated from the GPS vertical displacement time series, and mass changes in the studied regions are dominated by seasonal variability, we deducted it from GPS vertical displacement time series and focused on annual and semi-annual variations. Furthermore, in order to be consistent with the temporal resolution of GRACE/GFO monthly gravity field solutions, the daily vertical displacement observations were converted into monthly averages before inverting surface mass variations.
Figure 2 shows the detrended GPS vertical displacement time series (blue line) and their fitted seasonal time series (yellow line) for three selected permanent GPS stations (YNSM, HBZG, and GZGY) in Southwest China, and the monthly variations of vertical displacements derived from CSR GRACE-SH (red line) and CSR GFO-SH (green line) are also presented. Atmospheric and ocean no-tidal loading effects were added back to GRACE/GFO estimates using the GAC fields from AOD1B RL06 products. As shown in Figure 2, for some stations (e.g., YNSM and HBZG), GPS-measured and GRACE/GFO-derived vertical displacement time series show good consistency with each other, while for some other stations (e.g., GZGY), their consistency is poor, perhaps due to technical errors in GPS data processing.
Figure 3 shows the correlation coefficient between GPS-measured and GRACE-derived vertical displacement time series for each GPS station over the period January 2011 to June 2017. Here, the GRACE-derived vertical displacements were inferred from CSR GRACE-SH. The dots represent the distribution of permanent GPS stations, and colors of the dots represent the correlation coefficient values. In addition, the locations of YNSM, HBZG, GZGY, and KMIN station are also marked. Among the 85 permanent GPS stations over Southwest China, 71% of them had a correlation coefficient greater than 0.5, and 50% had a correlation coefficient greater than 0.6. In addition, most of the GPS stations with good correlations are located in the YNP and MRB regions. However, the correlations of individual GPS stations are obviously poor, like GZGY station in WRB and KMIN station in YNP. In addition to technical errors in GPS data processing, another important reason that accounts for the poor correlation between GPS-measured and GRACE/GFO-inferred vertical displacements is that GPS deformations are more sensitive to local mass loading. One of the criteria for identifying GPS abnormal stations is to check whether there is any obvious difference compared with the nearby stations. However, due to the relatively sparse distributions of GPS stations in the four studied regions, it is difficult to distinguish whether or not the abnormal station is affected by local effects. In our study, the GPS station data checked with poor correlation and consistency would be discreetly used based on the influence analysis of a single GPS station (e.g., GZGY and KMIN) on the inversion of mass changes over Southwest China.

3.2. Inversion Model for Surface Mass Variations

On the basis of elastic load theory, the Earth’s deformation caused by the redistribution of surface mass loading can be approximately expressed as elastic deformation, and elastic deformation induced by the mass variation can be expressed as the Green’s Function [54,55]. It is assumed that a uniform mass disk with an angle radius of α is placed on the surface of the Earth with spherical symmetry, self-gravitation, elasticity, and isotropy. The total mass of the disk is as follows [26]:
Δ M = 2 π h w ρ w 1 cos α R 2
in which h w and ρ w are the equivalent water height (EWH) and density of water, respectively, and R represents the radius of the Earth. The vertical displacement u induced by a mass loading Δ M can be expressed as [29]:
u θ = Δ M × R M e l = 0 h l P l cos θ
where M e is the mass of the Earth, θ is the angular distance, h l are the load Love numbers from Wang et al. [56], and P l are the Legendre polynomials of order l .
In order to better investigate surface mass variations, the whole studied area of Southwest China was divided into uniform square grid cells with a spatial resolution of 0.25 × 0.25°. Then, by substituting Equation (1) into Equation (2), the unknown EWH vector for the above grid cells can be solved based on least-squares estimation. However, the correlation of the GPS vertical displacement observations will result in a discrete ill-posed problem for inverting regional surface mass variations. To deal with the ill-posed problem, the Tikhonov regularization method was applied as follows [16]:
min A X u σ 2 2 + λ 2 L X 2
where A is the design matrix decided by Green’s Function, X is the unknown EWH vector that needs to be solved, u is the observation vector of GPS vertical displacements, σ is the noise standard deviation of the observations, L is the regularized Laplacian matrix [57] which can constrain the mass changes more stably between the adjacent grid cells, and λ is the regularization parameter which can be estimated by the generalized cross validation (GCV) method [58,59]. Finally, the EWH vector can be estimated as follows:
X = A T A + σ 2 λ 2 L T L 1 A T u
The uncertainty of surface mass variations inverted from GPS vertical displacements can be acquired from covariance analysis of the estimated EWH vector. The covariance matrix can be obtained by inverting the normal equation matrix, and the diagonal elements of covariance matrix are the error variances of estimated EWH.

4. Results and Discussions

4.1. Spatial Distributions of Surface Mass Variations Derived from GPS and GRACE/GFO

Both GPS-measured and GRACE/GFO-derived vertical displacement time series showed apparent seasonal variations over Southwest China, and these two types of vertical displacement time series showed good correlation and consistency except for some individual stations, e.g., GZGY station (see Figure 2 and Figure 3), and KMIN station (see Figure 3). As mentioned earlier, the poor correlation between GPS-measured and GRACE/GFO-inferred vertical displacements may have been caused by technical errors in GPS data processing or local effects on GPS station. In order to analyze the influence of a single GPS station (i.e., GZGY and KMIN) on the inversion of mass variations over Southwest China, a closed-loop simulation was conducted based on the simulated GPS vertical displacements. Using the locations of 85 permanent GPS stations of CMONOC, the simulated GPS vertical displacement data were generated from CSR GRACE-SH plus GAC fields in March and November 2011, respectively. For these two months, surface mass variations derived from the above input signals were different from each other, and the signals in March 2011 were much weaker than those in November 2011. For instance, the maximum, minimum, mean, and standard deviation of mass variations (expressed in terms of EWH) in March and November 2011 were 41.7 vs. 260.9, −163.9 vs. −211.5, −50.9 vs. 37.7, and 45.6 vs. 91.3 mm, respectively, in the four studied regions. Before inverting surface mass variations, Gaussian white noises with a standard deviation of 3 mm were added to the simulated GPS vertical displacements. Figure 4 shows the surface mass variations inverted by simulated GPS vertical displacements from a total of 85 GPS stations and 84 GPS stations (without GZGY or KMIN) over the studied regions in March (top) and November (bottom) 2011, respectively. Figure 5 shows the differences between the inversion results of 85 GPS stations and those of stations without GZGY or KMIN in March (top) and November (bottom) 2011, respectively.
As shown in Figure 4 and Figure 5, the influence of GZGY station on the inversion of mass variations over Southwest China is larger than that of KMIN station, and the RMS differences between the inversion results of 85 GPS stations and those of stations without GZGY or KMIN are 1.9 vs. 0.4 in March 2011, and 3.2 vs. 1.0 mm in November 2011, respectively. The reasons for this are twofold. Firstly, compared with KMIN station, there are no other stations around GZGY within ~200 km, and there is no constraint information available nearby. Secondly, there are only two GPS stations located in WRB region, and GZGY station may provide a representative sample of WRB, so whether GZGY station data was used or not would make an obvious influence on the inversion results, especially in the months with relatively larger signal variations (e.g., in November 2011). However, these differences caused by GZGY or KMIN station are very small compared to the mass change signals themselves (see Figure 4). Although the poor correlations between GPS-measured and GRACE/GFO-inferred vertical displacements were found at GZGY and KMIN stations, we still included GZGY and KMIN station data to invert surface mass variations over Southwest China based on the above analysis.
In order to examine the performance of GPS and GRACE/GFO for monitoring surface mass variations, we employed the GPS vertical displacement data from a total of 85 permanent GPS stations of CMONOC to invert mass variations over Southwest China, and compared them with GRACE/GFO observations. Figure 6 shows the spatial distributions of surface mass variations obtained from GPS vertical displacements and CSR RL06 GRACE/GFO SH solutions in February and July 2011, and January and July 2019, respectively, and the red dots in Figure 6a,d represent the location of GPS stations.
As presented in Figure 6, in February 2011, surface mass variations derived from CSR GRACE-SH show mainly negative anomalies in the YNP region, the northwestern of MRB, and JLRB, while indicating positive anomalies in WRB, the southeast of MRB, and JLRB. The spatial distribution of mass variations inverted by GPS vertical displacements is comparable to that observed by GRACE, but the positive anomalies in WRB, the southeast of MRB, and JLRB were smaller than those of CSR GRACE-SH estimates. In July 2011, the spatial distributions of the mass variations obtained by GPS vertical displacements and CSR GRACE-SH were almost the same. Among them, positive anomalies were mainly presented in the YNP and MRB regions, while negative anomalies appeared mostly in the JLRB and WRB regions. In January 2019, mass variations derived from CSR GFO-SH showed positive anomalies in nearly the entire YNP and WRB and southeast of JLRB and MRB, while the northwest of JLRB and MRB were dominated by negative anomalies. The spatial distribution of GPS-inferred mass variations was comparable to that of the CSR GFO-SH result, but the GPS magnitudes of negative anomalies in the northwest of JLRB and MRB were larger than those from CSR GFO-SH. The distribution of mass variations in July 2019 almost showed an inverse pattern compared with that of January 2019. In general, for most months (not all of them are shown in this paper), spatial distributions of mass variations obtained by GPS and GRACE/GFO showed good consistency. These results indicate that GPS can help independently monitor surface mass variations in our studied regions, and mass changes retrieved by GPS vertical displacements can be verified by GRACE/GFO results.

4.2. Seasonal Surface Mass Variations Revealed by GPS and GRACE/GFO

In order to better understand surface mass variations over Southwest China, seasonal variations of the inverted results from GPS and GRACE/GFO were further analyzed in the four regions, i.e., the YNP, MRB, JLRB, and WRB, respectively. The annual amplitudes of the mass variations estimated by GPS vertical displacements, CSR GRACE-SH, and CSR GRACE-M from January 2011 to June 2017 are presented in the 3 panels of Figure 7, respectively. As shown in Figure 7, the largest annual amplitudes are located in southwest YNP, while the northeast YNP and the other three regions show significantly smaller annual amplitudes. Furthermore, the distribution of annual amplitude obtained by GPS was in good agreement with those of CSR GRACE-SH and CSR GRACE-M. This suggests that the GPS-inferred mass changes are fairly reliable at seasonal time scales, based on comparisons with different GRACE results.
Figure 8 reveals mass change time series (top) and annual time series (bottom) for the four studied regions. The mass change time series are derived from GPS vertical displacements (from January 2011 to August 2019), CSR GRACE-SH (from August 2002 to June 2017), CSR GRACE-M (from August 2002 to June 2017), and CSR GFO-SH (from June 2018 to August 2019). The annual time series are fitted using CSR GRACE-SH plus CSR GFO-SH (CSR GRACE/GFO-SH), CSR GRACE-M, and GPS estimates, respectively.
As shown in Figure 8a, YNP has the largest area with the densest distribution of GPS stations (Table 1), and the amplitudes of GPS-derived surface mass change time series are larger than those of the other three regions. The annual amplitude and phase of GPS-derived mass change time series agree remarkably well with those of CSR GRACE/GFO-SH and CSR GRACE-M results. In addition, GPS-derived surface mass change time series from January 2011 to August 2019 provided a good supplement to fill in the gap between GRACE and GFO missions, and helped verify the results of the latest GFO measurements.
In the MRB region (see Figure 8b), surface mass changes derived from CSR GRACE-SH and CSR GRACE-M were in very good agreement during the period from August 2002 to June 2017. Annual variation time series derived from CSR GRACE/GFO-SH and CSR GRACE-M were nearly identical (indistinguishable in Figure 8b). During the period from January 2011 to June 2017, GPS-derived mass change time series also agreed well with the two GRACE estimates (i.e., CSR GRACE-SH and CSR GRACE-M), and the results from GPS and GFO showed good consistency as well. Likely due to the degraded data quality during the final stage of the GRACE mission (from November 2016 to June 2017), surface mass change time series derived from CSR GRACE-SH and CSR GRACE-M showed smaller seasonal variations compared with GPS estimates.
For the JLRB region shown in Figure 8c, GPS-derived surface mass changes agreed well with GRACE/GFO estimates. However, due to the sparse and uneven distribution of GPS stations in the JLRB region, there was a small annual phase lag between GPS and GRACE/GFO results. Figure 8d shows that, from 2011 to 2014, there were fairly obvious differences between GPS and GRACE/GFO derived mass changes in the WRB region. One reason is to think that the consistency of GPS-measured and GRACE-derived vertical displacements for the GZGY station located in the WRB was poor during this period, as shown in Figure 2. Because there are no other GPS stations around GZGY station and no constraints available nearby, it is difficult to determine whether the disagreement was caused by GPS technical errors or local effects. From the simulation in Section 4.1, we found that a single GPS station (i.e., GZGY) has little influence on the inversion of mass variations over Southwest China. We think the main reason is that WRB has the smallest area which is close to the limit of basin size that GRACE can observe, and the lowest station density (see Table 1) in the WRB region leads to larger uncertainty in the inverted surface mass variations.
Table 2 lists the amplitudes and phases of annual and semi-annual mass change time series for the four regions based on GPS and GRACE inversion results from January 2011 to June 2017, and the RMS differences of the mass change time series between GPS and GRACE (CSR GRACE-SH, CSR GRACE-M) are also presented. The semi-annual amplitudes were much smaller compared to the annual amplitudes in the four studied regions, and uncertainties at semi-annual time scale were expected to be relatively larger. For the YNP region, both annual amplitudes and phases agreed well between GPS and GRACE inversion results. The annual amplitudes from GRACE were slightly larger than that from GPS, and the annual amplitude from CSR GRACE-M agreed better with GPS estimates. However, the semi-annual amplitude of GPS estimates was notably larger than the two GRACE results (CSR GRACE-SH, CSR GRACE-M). Moreover, the RMS differences of mass change time series between GPS and CSR GRACE-SH, and GPS and CSR GRACE-M were ~24.61 and 26.11 Gt, respectively, which indicates that there are still obvious discrepancies between GPS and GRACE inversion results.
For the MRB region listed in Table 2, annual amplitudes from GRACE (CSR GRACE-SH, CSR GRACE-M) were smaller than that from GPS, while the GPS semi-annual amplitude was smaller than those of GRACE. Meanwhile, the RMS differences of mass change time series between GPS and CSR GRACE-SH, GPS and CSR GRACE-M were ~8.74 and 7.48 Gt, respectively. For the JLRB, annual amplitudes of GPS and GRACE estimates were in good agreement, but the CSR GRACE-M estimates showed slightly larger annual amplitude. Additionally, the semi-annual amplitudes of the three estimates were almost the same, while the GPS series showed a 3–4 month phase lag compared with GRACE series. The RMS differences of the mass change time series between GPS and CSR GRACE-SH, GPS, and CSR GRACE-M were ~7.82 and 6.33 Gt, respectively. As for the WRB, the annual amplitude of mass changes inferred from GPS was underestimated, while the semi-annual amplitude from GPS was larger than those from the two GRACE inversion results. The RMS differences of the mass change time series between GPS and CSR GRACE-SH, GPS, and CSR GRACE-M were ~4.42 and 3.96 Gt, respectively.
In order to verify the reliability of mass variations derived from GPS and GRACE in the four studied regions, we further estimated the uncertainties from these two types of data. As mentioned earlier, the uncertainty of surface mass variations inverted from GPS vertical displacements was acquired from covariance analysis of the estimated EWH. The error RMS of mass changes derived from monthly covariance matrices for the YNP, MRB, JLRB, and WRB regions were 8.58, 3.61, 3.50, and 1.91 Gt, respectively. The uncertainty in GRACE estimates were obtained from the RMS of mass variations over an ocean area that covered the same latitudinal range of the four studied regions. The error RMS derived from CSR GRACE-SH and CSR GRACE-M for the YNP, MRB, JLRB, and WRB regions were 3.58 vs. 4.37, 1.51 vs. 1.84, 1.46 vs. 1.78, and 0.80 vs. 0.97 Gt, respectively. Based on the law of error propagation, the uncertainty of mass change differences between GPS and GRACE can be estimated by RMS GPS 2 + RMS GRACE 2 , where RMS GPS and RMS GRACE are error RMS of GPS and GRACE estimates. Therefore, we can deduce the error RMS of differences of the mass changes between GPS and CSR GRACE-SH, GPS and CSR GRACE-M for the above four regions were 9.30 vs. 9.63, 3.91 vs. 4.05, 3.79 vs. 3.93, and 2.07 vs. 2.14 Gt, respectively. Compared with the RMS differences listed in Table 2, all the above error RMS in the four studied regions are smaller. However, if the area of the studied region is smaller, the uncertainty of the corresponding estimation results will be larger. Especially for the WRB, the error RMS of differences of the mass changes between GPS and GRACE (2.07 or 2.14 Gt) were close to their RMS differences (4.42 or 3.96 Gt). As mentioned before, one reason is that WRB has the smallest area which is close to the limit of basin size that GRACE can observe (although the leakage was corrected in this study). Another reason is that WRB has the lowest station density compared to the other three regions, which leads to larger uncertainty in GPS estimates.
To sum up, the good consistencies between GPS- and GRACE/GFO-derived surface mass change time series show that GPS could well detect the characteristics of seasonal variations over Southwest China (especially in YNP and MRB) from January 2011 to August 2019, and the GPS-derived mass change time series provided an alternative way to fill in the gap between the GRACE and GFO missions. In addition, there are still discrepancies between GPS and GRACE/GFO inversion results. The reasons for this are twofold. Firstly, GRACE/GFO estimates mainly reflect large-scale mass variations with spatial resolution greater than 300km and strongly constrain the total mass variations, while the GPS principally detects regional surface mass variation at a local point (thus constraining the distribution of surface mass), and relatively sparse GPS stations cannot fully infer the details of surface mass variations in the whole studied area. Secondly, differences of error sources contained in GPS vertical displacements and GRACE/GFO data, and different processing strategies (e.g., smooth filtering and constraint solving techniques) can also induce discrepancies between the inversion results of these two types of data.

4.3. GPS-Derived Surface Mass Variations to Bridge GRACE/GFO Mission Gap

As discussed above, surface mass change time series inverted by GPS and GRACE/GFO over Southwest China agree well, especially in the YNP and MRB regions. Because GPS can continuously monitor surface mass changes on regional scales, GPS-derived mass variations can serve as a supplementary data to bridge the GRACE/GFO mission gap. To minimize potential biases between GPS- and GRACE/GFO-derived mass changes due to differences in data processing procedures, we employed scale factor methods (both basin-scaled and pixel-scaled) adopted in Landerer et al. [34] to reduce the amplitude discrepancies between GPS- and GRACE/GFO-derived mass change time series. Specifically, we took CSR GRACE-SH-derived mass change time series from January 2011 to June 2017 as a reference, and estimated basin scale and pixel scale factors by minimizing the misfits between GPS-derived and GRACE-inferred surface mass change time series based on drainage basin (basin-scaled) and each grid cell (pixel-scaled), respectively. The feasibility of these scaling factors was tested by comparison of the scaled GPS-derived and GFO-inferred mass change time series during the period from June 2018 to August 2019.
The four panels of Figure 9 show the raw mass change time series (Raw GPS), basin-scaled (basin-scaled GPS), and pixel-scaled (pixel-scaled GPS) time series inverted by GPS vertical displacements from June 2017 to August 2019, and mass change time series observed by CSR GFO-SH from June 2018 to August 2019 in YNP, MRB, JLRB, and WRB, respectively. Their annual components of all the time series (over the period June 2018 to August 2019) for each region are also presented.
As shown in Figure 9, in the YNP and MBR regions, all of the surface mass change time series inverted by GPS and GFO show good consistency, and GPS-derived mass change time series scaled by basin-scaled and pixel-scaled methods are almost the same. The annual amplitudes of mass change time series obtained by basin-scaled and pixel-scaled methods were slightly smaller than those of the raw surface mass change time series inverted by GPS, but were much closer to the CSR GFO-SH results. As for JLRB and WRB, surface mass change time series inverted by GPS and GFO showed some phase differences, and the annual amplitudes of basin-scaled and pixel-scaled GPS mass changes were smaller than those of the CSR GFO-SH results. The discrepancies shown in the latter two regions demonstrate that the relatively sparse and uneven distributions of GPS stations limit the accuracy of GPS-inferred regional mass variations, and the limited spatial resolution of GFO also leads to larger uncertainty in mass change estimates.
Table 3 lists the annual amplitudes and phases for the four regions (i.e., the YNP, MRB, JLRB, and WRB) based on the raw GPS-derived (Raw GPS) mass change time series, basin-scaled (basin-scaled GPS), and pixel-scaled (pixel-scaled GPS) mass change time series, and GFO-observed (CSR GFO-SH) mass change time series from June 2018 to August 2019, and the RMS of mass change time series minus seasonal variations for CSR GFO-SH and those three GPS inversion results are also given for comparisons.
As shown in Table 3, in the YNP region, the annual amplitude of the raw GPS-derived mass change time series is much larger than that of GFO observations (72.69 vs. 59.19 Gt), but becomes substantially smaller after the scaling factor corrections. Moreover, the annual amplitude of basin-scaled GPS mass change time series agreed remarkably well with GFO observations (59.19 vs. 59.84 Gt), and the RMS residual also reduced correspondingly after the basin-scaled factor corrections (16.71 vs. 13.61 Gt). However, the RMS residual increased by using the pixel-scaled factor (16.71 vs. 18.04 Gt). For the MRB, the annual amplitude of raw GPS mass change time series was also significantly larger than that of GFO observations (22.42 vs. 12.53 Gt), and after scaling factor corrections, GPS annual amplitudes agreed notably better with GFO observations (10.03 or 14.71 vs. 12.53 Gt). The RMS residuals were significantly reduced (5.50 vs. 2.49 or 3.46 Gt). However, in the JLRB and WRB regions, neither basin-scaled nor pixel-scaled GPS mass change time series showed improved agreements with GFO observations, indicating the challenges of inverting regional mass changes when inadequate numbers of GPS stations are available. Certainly, due to the limited spatial resolution of GFO, the larger uncertainty of GFO estimates would fairly affect the comparisons. This may be evident in the JLRB region, where GFO observations did not show a clear and typical seasonal signal as in the other three regions. Nevertheless, the RMS residuals did show obvious decreases in the WRB region after the two scaled factors corrections (2.58 vs. 1.51 or 1.87 Gt). In addition, the RMS residuals of basin-scaled GPS time series were smaller than those of the pixel-scaled GPS time series in the four studied regions, which shows that basin-scaled factor method is more suitable to scale the GPS-derived mass change time series than the pixel-scaled factor method in our study.
In order to demonstrate the scaled GPS-derived surface mass variations can help bridge the GRACE/GFO mission gap, Figure 10 presents the long-term surface mass change time series derived from GRACE (CSR GRACE-SH, from August 2002 to June 2017), GFO (CSR GFO-SH, from June 2018 to August 2019), and GPS (basin-scaled GPS and pixel-scaled GPS, from July 2017 to May 2018) for the four regions (i.e., the YNP, MRB, JLRB, and WRB) in Southwest China. As shown in Figure 10, either basin-scaled or pixel-scaled mass change time series appear to well bridge the GRACE/GFO mission gap in these four regions, especially in YNP and MRB. Besides, due to the relatively dense distributions of GPS stations distributed in the YNP and MRB regions, seasonal variations during the one-year mission gap were also well captured. Seasonal variations in JLRB and WRB were significantly smaller than those in YNP and MRB, which also played a role in the poorer agreements between GPS and GFO results in the JLRB and WRB regions.
As can be seen from the above analysis, the annual amplitudes of GPS-derived mass change time series obtained from relatively dense distributions of GPS stations in the YNP and MRB are usually greater than that of GFO-inferred mass change time series. For sparsely uneven distributions of GPS stations in the JLRB and WRB, the annual amplitudes of GPS-derived mass change time series were smaller than that of GFO-inferred mass change time series. However, the scaled GPS-derived mass change time series can bridge the one-year mission gap between GRACE and GFO in these four regions, especially in the YNP and MRB where GPS stations are relatively densely distributed. From the RMS residuals listed in Table 3, we believe the basin-scaled factor method is more suitable to scale the GPS-derived mass change time series than the pixel-scaled factor method in our study.

5. Conclusions

In this study, GPS vertical displacement data from a total of 85 permanent GPS stations of CMONOC, the latest GRACE/GFO RL06 SH solutions, and GRACE RL06 mascon solutions were used to investigate surface mass variations over Southwest China. Spatial distributions and seasonal characteristics of surface mass variations inverted by GPS and GRACE/GFO were analyzed and compared in detail in the four studied regions (i.e., the YNP, MRB, JLRB, and WRB) in Southwest China. Moreover, the potential of using GPS-derived mass change time series to bridge the GRACE/GFO mission gap was also demonstrated.
For most months, spatial distributions of GPS and GRACE/GFO derived surface mass variations showed good consistency during the period of January 2011 to August 2019. Annual amplitudes and phases of GPS-derived surface mass changes agreed well with those of GRACE estimates (e.g., CSR GRACE-SH and CSR GRACE-M). Seasonal characteristics of surface mass changes over the period January 2011 to August 2019 were well captured by GPS inversions in each of the four regions, and GPS-derived mass change time series showed good agreements with GRACE/GFO observations, especially in the YNP and MRB regions where GPS stations are relatively densely distributed. However, evident discrepancies between GPS and GRACE/GFO derived mass changes still exist (e.g., in JLRB and WRB). One particular reason is that the relatively sparse GPS stations in JLRB and WRB are not adequate for inverting the details of regional mass variations. Another major reason is that both JLRB and WRB have smaller areas which are close to the limit of basin size that GRACE/GFO can observe (although the leakage was corrected), and the uncertainty of GRACE/GFO estimates would largely affect the comparisons in these cases. In addition, the error sources contained in GPS vertical displacement measurements and GRACE/GFO data, as well as different data processing strategies, will also induce discrepancies between GPS and GRACE/GFO derived mass change results.
Scale factor methods (both basin-scaled and pixel-scaled) were employed to reduce the amplitude discrepancies between the inversion results from GPS and GRACE/GFO. We showed that the scaled GPS-derived mass change time series can help bridge the gap between the GRACE and GFO missions, and the seasonal variations were accurately captured in the YNP and MRB regions by the scaled GPS results during the one-year mission gap. Therefore, GPS-derived surface mass variations can complement the data gaps of GRACE/GFO to obtain a long-term and continuous mass change time series at regional scales when adequate numbers of GPS stations are available. Moreover, the basin-scaled factor method was verified to be more suitable than the pixel-scaled factor method to scale the GPS-derived mass changes through the RMS residuals analysis. Future work will benefit from reprocessed GPS vertical displacement data and increased spatial coverage of GPS stations. In addition, as two independent geodetic techniques, joint inversions of surface mass variations from GPS and GRACE/GFO data may help increase the spatial and temporal resolution of GRACE/GFO-derived surface mass changes.

Author Contributions

Conceptualization, B.Z.; methodology, B.Z. and X.L.; software, X.L., B.Z., and T.L.; validation, X.L., J.C., and Q.L.; data curation, X.L. and T.L.; funding acquisition, B.Z.; writing—original draft preparation, X.L. and B.Z.; writing—review and editing, B.Z., J.C., and Q.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant Nos. 41974015, 41774032, and 41474019) and the National Key Research and Development Program of China (Grant No. 2018YFC1503503).

Acknowledgments

We are grateful to the China Earthquake Administration for the GPS data of the Crustal Movement Observation Network of China (CMONOC), and to the Center of Space Research (CSR) for providing the monthly GRACE/GFO gravity field solutions. The authors would like to thank the three anonymous reviewers for their insightful comments, which helped improve the presentation of the results.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tapley, B.D.; Bettadpur, S.; Ries, J.C.; Thompson, P.F.; Watkins, M.M. GRACE Measurements of Mass Variability in the Earth System. Science 2004, 305, 503–505. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Wahr, J.; Swenson, S.; Velicogna, I. Accuracy of GRACE mass estimates. Geophys. Res. Lett. 2006, 33, 178–196. [Google Scholar] [CrossRef] [Green Version]
  3. Chen, J. Satellite gravimetry and mass transport in the earth system. Geod. Geodyn. 2019, 10, 402–415. [Google Scholar] [CrossRef]
  4. Tapley, B.D.; Watkins, M.M.; Flechtner, F.; Reigber, C.; Bettadpur, S.; Rodell, M.; Sasgen, I.; Famiglietti, J.; Landerer, F.; Chambers, D.P.; et al. Contributions of GRACE to understanding climate change. Nat. Clim. Chang. 2019, 9, 358–369. [Google Scholar] [CrossRef]
  5. Swenson, S.; Wahr, J. Post-processing removal of correlated errors in GRACE data. Geophys. Res. Lett. 2006, 33. [Google Scholar] [CrossRef]
  6. Wouters, B.; Schrama, E. Improved accuracy of GRACE gravity solutions through empirical orthogonal function filtering of spherical harmonics. Geophys. Res. Lett. 2007, 34. [Google Scholar] [CrossRef] [Green Version]
  7. Chen, J.; Wilson, C.R.; Li, J.; Zhang, Z. Reducing leakage error in GRACE-observed long-term ice mass change: A case study in West Antarctica. J. Geod. 2015, 89, 925–940. [Google Scholar] [CrossRef]
  8. Sośnica, K.; Jaggi, A.; Meyer, U.; Thaller, D.; Beutler, G.; Arnold, D.; Dach, R. Time variable Earth’s gravity field from SLR satellites. J. Geod. 2015, 89, 945–960. [Google Scholar] [CrossRef] [Green Version]
  9. Cheng, M.; Ries, J. The unexpected signal in GRACE estimates in C20. J. Geod. 2017, 91, 897–914. [Google Scholar] [CrossRef]
  10. Blewitt, G.; Lavallée, D.; Clarke, P.J.; Nurutdinov, K. A New Global Mode of Earth Deformation: Seasonal Cycle Detected. Science 2001, 294, 2342–2345. [Google Scholar] [CrossRef] [Green Version]
  11. Blewitt, G.; Clarke, P. Inversion of Earth’s changing shape to weigh sea level in static equilibrium with surface mass redistribution. J. Geophys. Res. Solid Earth 2003, 108, 2311. [Google Scholar] [CrossRef] [Green Version]
  12. Wu, X.; Heflin, M.B.; Ivins, E.R.; Argus, D.F.; Webb, F.H. Large-scale global surface mass variations inferred from GPS measurements of load-induced deformation. Geophys. Res. Lett. 2003, 30. [Google Scholar] [CrossRef] [Green Version]
  13. Kusche, J.; Schrama, E. Surface mass redistribution inversion from global GPS deformation and Gravity Recovery and Climate Experiment (GRACE) gravity data. J. Geophys. Res. Space Phys. 2005, 110, 09409. [Google Scholar] [CrossRef] [Green Version]
  14. Zhang, X.; Jin, S.; Lu, X. Global Surface Mass Variations from Continuous GPS Observations and Satellite Altimetry Data. Remote Sens. 2017, 9, 1000. [Google Scholar] [CrossRef] [Green Version]
  15. Rietbroek, R.; Fritsche, M.; Dahle, C.; Brunnabend, S.-E.; Behnisch, M.; Kusche, J.; Flechtner, F.; Schröter, J.; Dietrich, R. Can GPS-Derived Surface Loading Bridge a GRACE Mission Gap? Surv. Geophys. 2014, 35, 1267–1283. [Google Scholar] [CrossRef] [Green Version]
  16. Argus, D.F.; Fu, Y.; Landerer, F. Seasonal variation in total water storage in California inferred from GPS observations of vertical land motion. Geophys. Res. Lett. 2014, 41, 1971–1980. [Google Scholar] [CrossRef]
  17. Fu, Y.; Argus, D.F.; Landerer, F. GPS as an independent measurement to estimate terrestrial water storage variations in Washington and Oregon. J. Geophys. Res. Solid Earth 2015, 120, 552–566. [Google Scholar] [CrossRef]
  18. Davis, J.L.; Elosegui, P.; Mitrovica, J.X.; Tamisiea, M.E. Climate-driven deformation of the solid Earth from GRACE and GPS. Geophys. Res. Lett. 2004, 31, 24605. [Google Scholar] [CrossRef] [Green Version]
  19. Fu, Y.; Freymueller, J. Seasonal and long-term vertical deformation in the Nepal Himalaya constrained by GPS and GRACE measurements. J. Geophys. Res. Space Phys. 2012, 117, 03407. [Google Scholar] [CrossRef]
  20. Nahmani, S.; Bock, O.; Bouin, M.-N.; Santamaria, A.; Boy, J.-P.; Collilieux, X.; Métivier, L.; Panet, I.; Genthon, P.; De Linage, C.; et al. Hydrological deformation induced by the West African Monsoon: Comparison of GPS, GRACE and loading models. J. Geophys. Res. Space Phys. 2012, 117. [Google Scholar] [CrossRef] [Green Version]
  21. Bevis, M.; Wahr, J.; Khan, S.A.; Madsen, F.B.; Brown, A.; Willis, M.; Kendrick, E.; Knudsen, P.; Box, J.; Van Dam, T.; et al. Bedrock displacements in Greenland manifest ice mass variations, climate cycles and climate change. Proc. Natl. Acad. Sci. USA 2012, 109, 11944–11948. [Google Scholar] [CrossRef] [Green Version]
  22. Godah, W.; Szelachowska, M.; Ray, J.D.; Krynski, J. Comparison of vertical deformations of the Earth’s surface obtained using GRACE-based GGMs and GNSS data—A case study of South-Eastern Poland. Acta Geodyn. Geomater. 2020, 17, 169–176. [Google Scholar] [CrossRef]
  23. Hao, M.; Freymueller, J.; Wang, Q.; Cui, D.; Qin, S. Vertical crustal movement around the southeastern Tibetan Plateau constrained by GPS and GRACE data. Earth Planet. Sci. Lett. 2016, 437, 1–8. [Google Scholar] [CrossRef]
  24. Zhan, W.; Li, F.; Hao, W.; Yan, J. Regional characteristics and influencing factors of seasonal vertical crustal motions in Yunnan, China. Geophys. J. Int. 2017, 210, 1295–1304. [Google Scholar] [CrossRef]
  25. Liu, R.; Zou, R.; Li, J.; Zhang, C.; Bin, Z.; Zhang, Y. Vertical Displacements Driven by Groundwater Storage Changes in the North China Plain Detected by GPS Observations. Remote Sens. 2018, 10, 259. [Google Scholar] [CrossRef] [Green Version]
  26. Wahr, J.; Khan, S.A.; Van Dam, T.; Liu, L.; Van Angelen, J.H.; Broeke, M.R.V.D.; Meertens, C. The use of GPS horizontals for loading studies, with applications to northern California and southeast Greenland. J. Geophys. Res. Solid Earth 2013, 118, 1795–1806. [Google Scholar] [CrossRef] [Green Version]
  27. King, M.; Moore, P.; Clarke, P.J.; Lavallee, D. Choice of optimal averaging radii for temporal GRACE gravity solutions, a comparison with GPS and satellite altimetry. Geophys. J. Int. 2006, 166, 1–11. [Google Scholar] [CrossRef] [Green Version]
  28. Van Dam, T.; Wahr, J.; Lavallee, D. A comparison of annual vertical crustal displacements from GPS and Gravity Recovery and Climate Experiment (GRACE) over Europe. J. Geophys. Res. Space Phys. 2007, 112, 03404. [Google Scholar] [CrossRef] [Green Version]
  29. Jin, S.; Zhang, T. Terrestrial Water Storage Anomalies Associated with Drought in Southwestern USA from GPS Observations. Surv. Geophys. 2016, 37, 1139–1156. [Google Scholar] [CrossRef]
  30. Enzminger, T.L.; Small, E.E.; Borsa, A.A. Accuracy of Snow Water Equivalent Estimated from GPS Vertical Displacements: A Synthetic Loading Case Study for Western U.S. Mountains. Water Resour. Res. 2018, 54, 581–599. [Google Scholar] [CrossRef]
  31. Zhang, B.; Yao, Y.; Fok, H.S.; Hu, Y.; Chen, Q. Potential Seasonal Terrestrial Water Storage Monitoring from GPS Vertical Displacements: A Case Study in the Lower Three-Rivers Headwater Region, China. Sensors 2016, 16, 1526. [Google Scholar] [CrossRef]
  32. Ferreira, V.G.; Ndehedehe, C.E.; Montecino, H.; Gourley, J.J.; Yuan, P.; Abdalla, A.; Mohammed, A.S. Prospects for Imaging Terrestrial Water Storage in South America Using Daily GPS Observations. Remote Sens. 2019, 11, 679. [Google Scholar] [CrossRef] [Green Version]
  33. Fok, H.S.; Liu, Y. An Improved GPS-Inferred Seasonal Terrestrial Water Storage Using Terrain-Corrected Vertical Crustal Displacements Constrained by GRACE. Remote Sens. 2019, 11, 1433. [Google Scholar] [CrossRef] [Green Version]
  34. Landerer, F.; Swenson, S.C. Accuracy of scaled GRACE terrestrial water storage estimates. Water Resour. Res. 2012, 48. [Google Scholar] [CrossRef]
  35. Zhao, B.; Huang, Y.; Zhang, C.; Wang, W.; Tan, K.; Du, R. Crustal deformation on the Chinese mainland during 1998–2014 based on GPS data. Geod. Geodyn. 2015, 6, 7–15. [Google Scholar] [CrossRef] [Green Version]
  36. Schmid, R.; Dach, R.; Collilieux, X.; Jäggi, A.; Schmitz, M.; Dilssner, F. Absolute IGS antenna phase center model igs08.atx: Status and potential improvements. J. Geod. 2015, 90, 343–364. [Google Scholar] [CrossRef] [Green Version]
  37. Boehm, J.; Niell, A.; Tregoning, P.; Schuh, H.; Böhm, J. Global Mapping Function (GMF): A new empirical mapping function based on numerical weather model data. Geophys. Res. Lett. 2006, 33, 07304. [Google Scholar] [CrossRef] [Green Version]
  38. Wahr, J.; Molenaar, M.; Bryan, F. Time variability of the Earth’s gravity field: Hydrological and oceanic effects and their possible detection using GRACE. J. Geophys. Res. Solid Earth 1998, 103, 30205–30229. [Google Scholar] [CrossRef]
  39. Landerer, F. Monthly Estimates of Degree-1 (Geocenter) Gravity Coefficients, Generated from GRACE (04-2002-06/2017) and GRACE-FO (06/2018 onward) RL06 Solutions; GRACE Technical Note 13, the GRACE Project; NASA Jet Propulsion Laboratory: Pasadena, CA, USA, 2019. [Google Scholar]
  40. Loomis, B.D.; Rachlin, K.E.; Wiese, D.N.; Landerer, F.W.; Luthcke, S.B. Replacing GRACE/GRACE-FO With Satellite Laser Ranging: Impacts on Antarctic Ice Sheet Mass Change. Geophys. Res. Lett. 2020, 47. [Google Scholar] [CrossRef]
  41. Chen, J.L.; Wilson, C.R.; Tapley, B.D.; Grand, S. GRACE detects coseismic and postseismic deformation from the Sumatra-Andaman e earthquake. Geophys. Res. Lett. 2007, 34, L13302. [Google Scholar] [CrossRef] [Green Version]
  42. Rodell, M.; Houser, P.; Jambor, U.; Gottschalck, J.; Mitchell, K.; Meng, C.-J.; Arsenault, K.R.; Cosgrove, B.; Radakovich, J.; Bosilovich, M.; et al. The Global Land Data Assimilation System. Bull. Am. Meteorol. Soc. 2004, 85, 381–394. [Google Scholar] [CrossRef] [Green Version]
  43. Chen, J.; Li, J.; Zhang, Z.; Ni, S. Long-term groundwater variations in Northwest India from satellite gravity measurements. Glob. Planet. Change 2014, 116, 130–138. [Google Scholar] [CrossRef] [Green Version]
  44. Chen, J.; Wilson, C.R.; Tapley, B.; Scanlon, B.; Guntner, A. Long-term groundwater storage change in Victoria, Australia from satellite gravity and in situ observations. Glob. Planet. Change 2016, 139, 56–65. [Google Scholar] [CrossRef] [Green Version]
  45. Wiese, D.N.; Landerer, F.; Watkins, M.M. Quantifying and reducing leakage errors in the JPL RL05M GRACE mascon solution. Water Resour. Res. 2016, 52, 7490–7502. [Google Scholar] [CrossRef]
  46. Watkins, M.M.; Wiese, D.N.; Yuan, D.-N.; Boening, C.; Landerer, F. Improved methods for observing Earth’s time variable mass distribution with GRACE using spherical cap mascons. J. Geophys. Res. Solid Earth 2015, 120, 2648–2671. [Google Scholar] [CrossRef]
  47. Save, H.; Bettadpur, S.; Tapley, B.D. High-resolution CSR GRACE RL05 mascons. J. Geophys. Res. Solid Earth 2016, 121, 7547–7569. [Google Scholar] [CrossRef]
  48. Save, H. “CSR GRACE RL06 Mascon Solutions”, Texas Data Repository Dataverse, V1. 2019. Available online: https://doi.org/10.18738/T8/UN91VR (accessed on 25 May 2020).
  49. Chen, J.; Wilson, C.R.; Tapley, B.D.; Save, H.; Cretaux, J.-F. Long-term and seasonal Caspian Sea level change from satellite gravity and altimeter measurements. J. Geophys. Res. Solid Earth 2017, 122, 2274–2290. [Google Scholar] [CrossRef]
  50. Ray, J.; Altamimi, Z.; Collilieux, X.; Van Dam, T. Anomalous harmonics in the spectra of GPS position estimates. GPS Solut. 2007, 12, 55–64. [Google Scholar] [CrossRef]
  51. Glaser, S.; Fritsche, M.; Sośnica, K.; Rodriguez-Solano, C.J.; Wang, K.; Dach, R.; Hugentobler, U.; Rothacher, M.; Dietrich, R. A consistent combination of GNSS and SLR with minimum constraints. J. Geod. 2015, 89, 1165–1180. [Google Scholar] [CrossRef]
  52. Nikolaidis, R. Observation of Geodetic and Seismic Deformation with the Global Positioning System. Ph.D. Thesis, University of California San Deigo, San Deigo, CA, USA, 2002. [Google Scholar]
  53. Oliveira, P.; Gomes, L. Interpolation of signals with missing data using Principal Component Analysis. Multidimens. Syst. Signal Process. 2009, 21, 25–43. [Google Scholar] [CrossRef]
  54. Farrell, W.E. Deformation of the Earth by surface loads. Rev. Geophys. 1972, 10, 761. [Google Scholar] [CrossRef]
  55. Rajner, M.; Liwosz, T. Studies of crustal deformation due to hydrological loading on GPS height estimates. Geod. Cartogr. 2011, 60, 135–144. [Google Scholar] [CrossRef] [Green Version]
  56. Wang, H.; Xiang, L.; Jia, L.; Jiang, L.; Wang, Z.; Hu, B.; Gao, P. Load Love numbers and Green’s functions for elastic Earth models PREM, iasp91, ak135, and modified models with refined crustal structure from Crust 2.0. Comp. Geosci. 2012, 49, 190–199. [Google Scholar] [CrossRef]
  57. Van Dam, E.; Omidi, G. Graphs whose normalized Laplacian has three eigenvalues. Linear Algebra Appl. 2011, 435, 2560–2569. [Google Scholar] [CrossRef]
  58. Golub, G.H.; Wahba, H.G. Generalized Cross-Validation as a Method for Choosing a Good Ridge Parameter. Technometrics 1979, 21, 215–223. [Google Scholar] [CrossRef]
  59. Liao, H.; Ng, M.K. Blind Deconvolution Using Generalized Cross-Validation Approach to Regularization Parameter Estimation. IEEE Trans. Image Process. 2011, 20, 670–680. [Google Scholar] [CrossRef]
Figure 1. The four studied regions in Southwest China: The abbreviations YNP, MRB, JLRB, and WRB in blue indicate Yunnan Province, Min River Basin, Jialing River Basin, and Wu River Basin, respectively. The red dots represent the distribution of permanent Global Positioning System (GPS) stations, the black lines are boundaries of the four regions, and the background represents the digital elevation model (DEM) information of Southwest China.
Figure 1. The four studied regions in Southwest China: The abbreviations YNP, MRB, JLRB, and WRB in blue indicate Yunnan Province, Min River Basin, Jialing River Basin, and Wu River Basin, respectively. The red dots represent the distribution of permanent Global Positioning System (GPS) stations, the black lines are boundaries of the four regions, and the background represents the digital elevation model (DEM) information of Southwest China.
Remotesensing 12 01835 g001
Figure 2. The detrended GPS vertical displacement time series (blue line) and their fitted seasonal time series (yellow line) for three permanent GPS stations (YNSM, HBZG and GZGY) in Southwest China, compared with monthly variations of vertical displacement derived from Center for Space Research (CSR) RL06 Gravity Recovery and Climate Experiment (GRACE) Spherical Harmonic (SH) solutions (CSR GRACE-SH, red line) and CSR RL06 GRACE Follow-On (GFO) SH solutions (CSR GFO-SH, green line) with AOD1B atmosphere combined with ocean (GAC) contributions restored.
Figure 2. The detrended GPS vertical displacement time series (blue line) and their fitted seasonal time series (yellow line) for three permanent GPS stations (YNSM, HBZG and GZGY) in Southwest China, compared with monthly variations of vertical displacement derived from Center for Space Research (CSR) RL06 Gravity Recovery and Climate Experiment (GRACE) Spherical Harmonic (SH) solutions (CSR GRACE-SH, red line) and CSR RL06 GRACE Follow-On (GFO) SH solutions (CSR GFO-SH, green line) with AOD1B atmosphere combined with ocean (GAC) contributions restored.
Remotesensing 12 01835 g002
Figure 3. Correlation coefficient between GPS-measured and GRACE-derived vertical displacement time series for each GPS station over the period January 2011 to June 2017. The GRACE-derived vertical displacements are inferred from CSR RL06 GRACE SH solutions (CSR GRACE-SH). The dots represent the distribution of permanent GPS stations, and their colors represent the correlation coefficient values. The locations of YNSM, HBZG, GZGY, and KMIN station are marked, and the abbreviations YNP, MRB, JLRB, and WRB in blue indicate Yunnan Province, Min River Basin, Jialing River Basin, and Wu River Basin, respectively.
Figure 3. Correlation coefficient between GPS-measured and GRACE-derived vertical displacement time series for each GPS station over the period January 2011 to June 2017. The GRACE-derived vertical displacements are inferred from CSR RL06 GRACE SH solutions (CSR GRACE-SH). The dots represent the distribution of permanent GPS stations, and their colors represent the correlation coefficient values. The locations of YNSM, HBZG, GZGY, and KMIN station are marked, and the abbreviations YNP, MRB, JLRB, and WRB in blue indicate Yunnan Province, Min River Basin, Jialing River Basin, and Wu River Basin, respectively.
Remotesensing 12 01835 g003
Figure 4. Surface mass variations inverted by simulated GPS vertical displacements in March (top) and November (bottom) 2011, respectively. (a) and (d) show the spatial distributions inverted from a total of 85 GPS stations; (b) and (e) show the spatial distributions inverted from those GPS stations without GZGY; (c) and (f) show the spatial distributions inverted from those GPS stations without KMIN. The red dots represent the location of GPS stations, but the locations of GZGY and KMIN are indicated as red squares.
Figure 4. Surface mass variations inverted by simulated GPS vertical displacements in March (top) and November (bottom) 2011, respectively. (a) and (d) show the spatial distributions inverted from a total of 85 GPS stations; (b) and (e) show the spatial distributions inverted from those GPS stations without GZGY; (c) and (f) show the spatial distributions inverted from those GPS stations without KMIN. The red dots represent the location of GPS stations, but the locations of GZGY and KMIN are indicated as red squares.
Remotesensing 12 01835 g004
Figure 5. Differences of surface mass variations derived from simulated GPS vertical displacements in March (top) and November (bottom) 2011, respectively. (a) and (c) show the differences between the inversion results of 85 GPS stations and those of stations without GZGY; (b) and (d) show the differences between the inversion results of 85 GPS stations and those of stations without KMIN. The red dots represent the location of GPS stations, but the locations of GZGY and KMIN are indicated as red squares.
Figure 5. Differences of surface mass variations derived from simulated GPS vertical displacements in March (top) and November (bottom) 2011, respectively. (a) and (c) show the differences between the inversion results of 85 GPS stations and those of stations without GZGY; (b) and (d) show the differences between the inversion results of 85 GPS stations and those of stations without KMIN. The red dots represent the location of GPS stations, but the locations of GZGY and KMIN are indicated as red squares.
Remotesensing 12 01835 g005
Figure 6. Surface mass variations derived from GPS vertical displacements and CSR RL06 GRACE/GFO Spherical Harmonic (SH) solutions in February and July 2011, and January and July 2019, respectively. (ad) show the spatial distributions inverted by GPS in February and July 2011, and January and July 2019, respectively, and the red dots represent the location of GPS stations; (e,f) show the spatial distributions inferred from CSR RL06 GRACE SH solutions (CSR GRACE-SH) in February and July 2011, respectively; (g,h) show the spatial distributions inferred from CSR RL06 GFO SH solutions (CSR GFO-SH) in January and July 2019, respectively.
Figure 6. Surface mass variations derived from GPS vertical displacements and CSR RL06 GRACE/GFO Spherical Harmonic (SH) solutions in February and July 2011, and January and July 2019, respectively. (ad) show the spatial distributions inverted by GPS in February and July 2011, and January and July 2019, respectively, and the red dots represent the location of GPS stations; (e,f) show the spatial distributions inferred from CSR RL06 GRACE SH solutions (CSR GRACE-SH) in February and July 2011, respectively; (g,h) show the spatial distributions inferred from CSR RL06 GFO SH solutions (CSR GFO-SH) in January and July 2019, respectively.
Remotesensing 12 01835 g006
Figure 7. Comparisons of annual amplitudes of surface mass variations estimated by GPS vertical displacements, GRACE SH solutions, and mascon solutions during the period from January 2011 to June 2017. (a) Annual amplitudes inverted by GPS, the red dots indicate the location of GPS stations; (b,c) annual amplitudes inferred from CSR RL06 GRACE SH solutions (CSR GRACE-SH) and CSR RL06 GRACE mascon solutions (CSR GRACE-M), respectively.
Figure 7. Comparisons of annual amplitudes of surface mass variations estimated by GPS vertical displacements, GRACE SH solutions, and mascon solutions during the period from January 2011 to June 2017. (a) Annual amplitudes inverted by GPS, the red dots indicate the location of GPS stations; (b,c) annual amplitudes inferred from CSR RL06 GRACE SH solutions (CSR GRACE-SH) and CSR RL06 GRACE mascon solutions (CSR GRACE-M), respectively.
Remotesensing 12 01835 g007
Figure 8. Surface mass change time series (top) and annual time series (bottom) for the (a) Yunnan Province (YNP), (b) Min River Basin (MRB), (c) Jialing River Basin (JLRB), and (d) Wu River Basin (WRB), respectively. The mass change time series are derived from GPS vertical displacements (GPS, from January 2011 to August 2019), CSR RL06 GRACE SH solutions (CSR GRACE-SH, from August 2002 to June 2017), CSR RL06 GRACE mascon solutions (CSR GRACE-M, from August 2002 to June 2017), and CSR RL06 GFO SH solutions (CSR GFO-SH, from June 2018 to August 2019). The annual time series are fitted using CSR GRACE-SH plus CSR GFO-SH (CSR GRACE/GFO-SH), CSR GRACE-M, and GPS estimates, respectively.
Figure 8. Surface mass change time series (top) and annual time series (bottom) for the (a) Yunnan Province (YNP), (b) Min River Basin (MRB), (c) Jialing River Basin (JLRB), and (d) Wu River Basin (WRB), respectively. The mass change time series are derived from GPS vertical displacements (GPS, from January 2011 to August 2019), CSR RL06 GRACE SH solutions (CSR GRACE-SH, from August 2002 to June 2017), CSR RL06 GRACE mascon solutions (CSR GRACE-M, from August 2002 to June 2017), and CSR RL06 GFO SH solutions (CSR GFO-SH, from June 2018 to August 2019). The annual time series are fitted using CSR GRACE-SH plus CSR GFO-SH (CSR GRACE/GFO-SH), CSR GRACE-M, and GPS estimates, respectively.
Remotesensing 12 01835 g008
Figure 9. Top of each panel: The raw surface mass change time series (Raw GPS), basin-scaled (basin-scaled GPS), and pixel-scaled (pixel-scaled GPS) time series inverted by GPS vertical displacements from June 2017 to August 2019, and the surface mass change time series modeled by CSR RL06 GFO SH solutions (CSR GFO-SH) from June 2018 to August 2019 for the (a) Yunnan Province (YNP), (b) Min River Basin (MRB), (c) Jialing River Basin (JLRB), and (d) Wu River Basin (WRB), respectively. Bottom of each panel: Annual components of the mass change time series for each region from June 2018 to August 2019.
Figure 9. Top of each panel: The raw surface mass change time series (Raw GPS), basin-scaled (basin-scaled GPS), and pixel-scaled (pixel-scaled GPS) time series inverted by GPS vertical displacements from June 2017 to August 2019, and the surface mass change time series modeled by CSR RL06 GFO SH solutions (CSR GFO-SH) from June 2018 to August 2019 for the (a) Yunnan Province (YNP), (b) Min River Basin (MRB), (c) Jialing River Basin (JLRB), and (d) Wu River Basin (WRB), respectively. Bottom of each panel: Annual components of the mass change time series for each region from June 2018 to August 2019.
Remotesensing 12 01835 g009
Figure 10. The long-term surface mass change time series inverted by joint of GRACE/GFO and GPS from August 2002 to August 2019 for the (a) Yunnan Province (YNP), (b) Min River Basin (MRB), (c) Jialing River Basin (JLRB), and (d) Wu River Basin (WRB), respectively. The joint inversion results include the mass change time series derived from CSR RL06 GRACE SH solutions (CSR GRACE-SH, from August 2002 to June 2017), CSR RL06 GFO SH solutions (CSR GFO-SH, from June 2018 to August 2019), and the mission gap is filled with scaled GPS-derived mass change time series (Basin-scaled GPS and Pixel-scaled GPS, from July 2017 to May 2018).
Figure 10. The long-term surface mass change time series inverted by joint of GRACE/GFO and GPS from August 2002 to August 2019 for the (a) Yunnan Province (YNP), (b) Min River Basin (MRB), (c) Jialing River Basin (JLRB), and (d) Wu River Basin (WRB), respectively. The joint inversion results include the mass change time series derived from CSR RL06 GRACE SH solutions (CSR GRACE-SH, from August 2002 to June 2017), CSR RL06 GFO SH solutions (CSR GFO-SH, from June 2018 to August 2019), and the mission gap is filled with scaled GPS-derived mass change time series (Basin-scaled GPS and Pixel-scaled GPS, from July 2017 to May 2018).
Remotesensing 12 01835 g010
Table 1. Areas, numbers of stations, and station densities of the four study regions in Southwest China.
Table 1. Areas, numbers of stations, and station densities of the four study regions in Southwest China.
Study RegionArea
(km2)
Number of StationsStation Density 1
( stations 1000000   km 2 )
Yunnan Province428,154.752660.73
Min River Basin190,205.44842.06
Jialing River Basin187,125.19526.72
Wu River Basin97,797.94220.45
1 The station density for each region is expressed by the number of stations per 1,000000 km2.
Table 2. Statistics of amplitudes and phases of annual and semi-annual surface mass change time series for the four studied regions, the Yunnan Province (YNP), Min River Basin (MRB), Jialing River Basin (JLRB), and Wu River Basin (WRB), based on GPS and GRACE inversion results from January 2011 to June 2017, and the Root Mean Square (RMS) differences of the mass change time series between GPS and CSR GRACE-SH, GPS and CSR GRACE-M, respectively. The abbreviations of CSR-SH and CSR-M in the table represent CSR GRACE-SH and CSR GRACE-M, respectively.
Table 2. Statistics of amplitudes and phases of annual and semi-annual surface mass change time series for the four studied regions, the Yunnan Province (YNP), Min River Basin (MRB), Jialing River Basin (JLRB), and Wu River Basin (WRB), based on GPS and GRACE inversion results from January 2011 to June 2017, and the Root Mean Square (RMS) differences of the mass change time series between GPS and CSR GRACE-SH, GPS and CSR GRACE-M, respectively. The abbreviations of CSR-SH and CSR-M in the table represent CSR GRACE-SH and CSR GRACE-M, respectively.
Study
Region
Data TypeAnnual (Semi-Annual) Amplitude 1 (Gt)Annual (Semi-Annual)
Phase 1 (Year)
RMS
(Gt)
YNPGPS52.00 ± 3.24 (20.00 ± 3.24)0.09 ± 0.06 (−0.35 ± 0.16)/
CSR-SH55.71 ± 2.16 (11.28 ± 2.16)0.06 ± 0.04 (0.00 ± 0.19)24.61
CSR-M53.66 ± 2.13 (10.15 ± 2.13)0.31 ± 0.04 (0.49 ± 0.21)26.11
MRBGPS10.11 ± 1.42 (1.46 ± 1.41)0.01 ± 0.14 (−0.24 ± 0.97)/
CSR-SH8.75 ± 0.82 (2.49 ± 0.82)0.48 ± 0.09 (−0.33 ± 0.33)8.74
CSR-M8.88 ± 0.48 (2.73 ± 0.48)0.48 ± 0.05 (−0.43 ± 0.18)7.48
JLRBGPS10.45 ± 0.97 (3.14 ± 0.98)0.71 ± 0.09 (0.67 ± 0.31)/
CSR-SH9.78 ± 0.92 (3.41 ± 0.92)0.45 ± 0.10 (0.17 ± 0.27)7.82
CSR-M11.46 ± 0.59 (3.46 ± 0.58)0.37 ± 0.05 (0.30 ± 0.17)6.33
WRBGPS3.60 ± 0.60 (1.38 ± 0.60)0.94 ± 0.17 (0.51 ± 0.43)/
CSR-SH5.63 ± 0.43 (0.45 ± 0.43)0.43 ± 0.08 (−0.02 ± 0.97)4.42
CSR-M6.04 ± 0.37 (0.97 ± 0.37)0.57 ± 0.06 (0.70 ± 0.38)3.96
1 The annual amplitude and phase, and semi-annual amplitude and phase are fitted by A 1 sin ( 2 π t + φ 1 ) and A 2 sin ( 4 π t + φ 2 ) , respectively, where A 1 and φ 1 represent the annual amplitude and phase, and A 2 and φ 2 represent the semi-annual amplitude and phase.
Table 3. Statistics of annual amplitudes and phases for the four regions (the Yunnan Province (YNP), Min River Basin (MRB), Jialing River Basin (JLRB) and Wu River Basin (WRB)) based on GPS-derived (Raw GPS, basin-scaled GPS, and pixel-scaled GPS) and GFO-observed (CSR GFO-SH) mass change time series from June 2018 to August 2019, and the RMS of mass change time series minus seasonal variations for CSR GFO-SH, Raw GPS, basin-scaled GPS, and pixel-scaled GPS, respectively.
Table 3. Statistics of annual amplitudes and phases for the four regions (the Yunnan Province (YNP), Min River Basin (MRB), Jialing River Basin (JLRB) and Wu River Basin (WRB)) based on GPS-derived (Raw GPS, basin-scaled GPS, and pixel-scaled GPS) and GFO-observed (CSR GFO-SH) mass change time series from June 2018 to August 2019, and the RMS of mass change time series minus seasonal variations for CSR GFO-SH, Raw GPS, basin-scaled GPS, and pixel-scaled GPS, respectively.
Study
Region
Data
Type
Annual Amplitude 1 (Gt)Annual
Phase 1 (Year)
RMS
(Gt)
YNPCSR GFO-SH59.19 ± 8.500.74 ± 0.128.12
Raw GPS72.69 ± 6.660.17 ± 0.0916.71
Basin-scaled GPS59.84 ± 7.160.77 ± 0.1113.61
Pixel-scaled GPS56.00 ± 7.430.74 ± 0.1218.04
MRBCSR GFO-SH12.53 ± 1.930.12 ± 0.132.95
Raw GPS22.42 ± 1.910.58 ± 0.095.50
Basin-scaled GPS10.03 ± 0.870.15 ± 0.092.49
Pixel-scaled GPS14.71 ± 1.330.31 ± 0.093.46
JLRBCSR GFO-SH12.96 ± 2.630.58 ± 0.184.05
Raw GPS7.54 ± 1.070.13 ± 0.133.05
Basin-scaled GPS4.65 ± 0.880.06 ± 0.171.85
Pixel-scaled GPS8.30 ± 1.990.06 ± 0.222.91
WRBCSR GFO-SH7.05 ± 1.160.33 ± 0.151.88
Raw GPS5.99 ± 1.310.57 ± 0.222.58
Basin-scaled GPS3.56 ± 0.740.56 ± 0.211.51
Pixel-scaled GPS5.62 ± 0.720.71 ± 0.121.87
1 The annual amplitude and phase are fitted by A 1 sin ( 2 π t + φ 1 ) , where A 1 and φ 1 represent the annual amplitude and phase.

Share and Cite

MDPI and ACS Style

Zhong, B.; Li, X.; Chen, J.; Li, Q.; Liu, T. Surface Mass Variations from GPS and GRACE/GFO: A Case Study in Southwest China. Remote Sens. 2020, 12, 1835. https://doi.org/10.3390/rs12111835

AMA Style

Zhong B, Li X, Chen J, Li Q, Liu T. Surface Mass Variations from GPS and GRACE/GFO: A Case Study in Southwest China. Remote Sensing. 2020; 12(11):1835. https://doi.org/10.3390/rs12111835

Chicago/Turabian Style

Zhong, Bo, Xianpao Li, Jianli Chen, Qiong Li, and Tao Liu. 2020. "Surface Mass Variations from GPS and GRACE/GFO: A Case Study in Southwest China" Remote Sensing 12, no. 11: 1835. https://doi.org/10.3390/rs12111835

APA Style

Zhong, B., Li, X., Chen, J., Li, Q., & Liu, T. (2020). Surface Mass Variations from GPS and GRACE/GFO: A Case Study in Southwest China. Remote Sensing, 12(11), 1835. https://doi.org/10.3390/rs12111835

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop