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.

[1]  David LaVallee,et al.  Higher‐order ionospheric effects on the GPS reference frame and velocities , 2009 .

[2]  Felix W. Landerer,et al.  GPS as an independent measurement to estimate terrestrial water storage variations in Washington and Oregon , 2015 .

[3]  Zhao Li,et al.  Comparative analysis of different environmental loading methods and their impacts on the GPS height time series , 2013, Journal of Geodesy.

[4]  Per Christian Hansen,et al.  Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..

[5]  Yuji Yagi,et al.  Waveform inversion for seismic source processes using ABIC with two sorts of prior constraints: Comparison between proper and improper formulations , 2003 .

[6]  Karl-Rudolf Koch,et al.  Maximum likelihood estimate of variance components , 1986 .

[7]  Janusz Bogusz,et al.  On the significance of periodic signals in noise analysis of GPS station coordinates time series , 2016, GPS Solutions.

[8]  Xavier Collilieux,et al.  Comparison of very long baseline interferometry, GPS, and satellite laser ranging height residuals from ITRF2005 using spectral and correlation methods , 2007 .

[9]  Felix W. Landerer,et al.  Seasonal variation in total water storage in California inferred from GPS observations of vertical land motion , 2014 .

[10]  Mike P. Stewart,et al.  Aliased tidal signatures in continuous GPS height time series , 2003 .

[11]  S. Swenson,et al.  Post‐processing removal of correlated errors in GRACE data , 2006 .

[12]  Linguo Yuan,et al.  Seasonal crustal vertical deformation induced by environmental mass loading in mainland China derived from GPS, GRACE and surface loading models , 2017 .

[13]  Shin‐Chan Han,et al.  Statistical Downscaling of GRACE‐Derived Groundwater Storage Using ET Data in the North China Plain , 2018, Journal of Geophysical Research: Atmospheres.

[14]  S. Bettadpur,et al.  Modeling Earth deformation from monsoonal flooding in Bangladesh using hydrographic, GPS, and Gravity Recovery and Climate Experiment (GRACE) data , 2010 .

[15]  J. Wahr,et al.  A comparison of annual vertical crustal displacements from GPS and Gravity Recovery and Climate Experiment (GRACE) over Europe , 2007 .

[16]  G. Blewitt,et al.  A New Global Mode of Earth Deformation: Seasonal Cycle Detected , 2001, Science.

[17]  Hirotugu Akaike,et al.  Likelihood and the Bayes procedure , 1980 .

[18]  J. Kusche,et al.  Decorrelated GRACE time-variable gravity solutions by GFZ, and their validation using a hydrological model , 2009 .

[19]  Bin Wang,et al.  Multi-scale climate variability of the South China Sea monsoon: A review , 2009 .

[20]  Yuji Yagi,et al.  A method for the joint inversion of geodetic and seismic waveform data using ABIC: application to the 1997 Manyi, Tibet, earthquake , 2014 .

[21]  Wu Chen,et al.  Thermal Effects on Vertical Displacement of GPS Stations in China , 2010 .

[22]  Yehuda Bock,et al.  Error analysis of continuous GPS position time series , 2004 .

[23]  F. Bryan,et al.  Time variability of the Earth's gravity field: Hydrological and oceanic effects and their possible detection using GRACE , 1998 .

[24]  Jeffrey P. Walker,et al.  THE GLOBAL LAND DATA ASSIMILATION SYSTEM , 2004 .

[25]  Na Wei,et al.  Contributions of thermoelastic deformation to seasonal variations in GPS station position , 2017, GPS Solutions.

[26]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[27]  Shuanggen Jin,et al.  Terrestrial Water Storage Anomalies Associated with Drought in Southwestern USA from GPS Observations , 2016, Surveys in Geophysics.

[28]  Geoffrey Blewitt,et al.  Crustal displacements due to continental water loading , 2001 .

[29]  D. Alsdorf,et al.  Seasonal fluctuations in the mass of the Amazon River system and Earth's elastic response , 2005 .

[30]  H. Fok.  Ocean Tides Modeling using Satellite Altimetry , 2012 .

[31]  J. Freymueller,et al.  Seasonal hydrological loading in southern Alaska observed by GPS and GRACE , 2012 .

[32]  Peiliang Xu,et al.  Determination of surface gravity anomalies using gradiometric observables , 1992 .

[33]  O. Francis,et al.  Modelling the global ocean tides: modern insights from FES2004 , 2006 .

[34]  Simon D. P. Williams,et al.  Fast error analysis of continuous GNSS observations with missing data , 2013, Journal of Geodesy.

[35]  Yibin Yao,et al.  Potential Seasonal Terrestrial Water Storage Monitoring from GPS Vertical Displacements: A Case Study in the Lower Three-Rivers Headwater Region, China , 2016, Sensors.

[36]  H. Schuh,et al.  Short Note: A global model of pressure and temperature for geodetic applications , 2007 .

[37]  Yuebing Li,et al.  Green's function of the deformation of the Earth as a result of atmospheric loading , 2004 .

[38]  Caijun Xu,et al.  Joint inversion of GPS, InSAR and teleseismic data sets for the rupture process of the 2015 Gorkha, Nepal, earthquake using a generalized ABIC method , 2017 .

[39]  W. Farrell.  Deformation of the Earth by surface loads , 1972 .

[40]  Shfaqat Abbas Khan,et al.  Geodetic measurements of postglacial adjustments in Greenland , 2008 .

[41]  Z. Martinec,et al.  Applying local Green's functions to study the influence of the crustal structure on hydrological loading displacements , 2015 .

[42]  Yehuda Bock,et al.  Southern California permanent GPS geodetic array: Spatial filtering of daily positions for estimating coseismic and postseismic displacements induced by the 1992 Landers earthquake , 1997 .

[43]  P. Teunissen,et al.  Least-squares variance component estimation , 2008 .

[44]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

[45]  J. Ray,et al.  Anomalous harmonics in the spectra of GPS position estimates , 2008 .

[46]  C. K. Shum,et al.  Earth Surface Deformation in the North China Plain Detected by Joint Analysis of GRACE and GPS Data , 2014, Sensors.

[47]  Z. Altamimi,et al.  ITRF2008: an improved solution of the international terrestrial reference frame , 2011 .

[48]  Nico Sneeuw,et al.  Annual variations of monsoon and drought detected by GPS: A case study in Yunnan, China , 2017, Scientific Reports.

[49]  Shuanggen Jin,et al.  Large-scale variations of global groundwater from satellite gravimetry and hydrological models, 2002–2012 , 2013 .

[50]  M. Stone.  Cross‐Validatory Choice and Assessment of Statistical Predictions , 1976 .

[51]  Yehuda Bock,et al.  Spatiotemporal filtering using principal component analysis and Karhunen-Loeve expansion approaches for regional GPS network analysis , 2006 .

[52]  Chris Rizos,et al.  The International GNSS Service in a changing landscape of Global Navigation Satellite Systems , 2009 .

[53]  X. Yin,et al.  Assessment of Contributions of Climatic Variation and Human Activities to Streamflow Changes in the Lancang River, China , 2014, Water Resources Management.

[54]  Masaya Matsuura,et al.  Geodetic data inversion using a Bayesian information criterion for spatial distribution of fault slip , 1992 .

[55]  J. Marengo,et al.  Extreme seasonal droughts and floods in Amazonia: causes, trends and impacts , 2015 .

[56]  Pedro Elosegui,et al.  Climate‐driven deformation of the solid Earth from GRACE and GPS , 2004 .

[57]  James S. Famiglietti,et al.  Downscaling GRACE Remote Sensing Datasets to High-Resolution Groundwater Storage Change Maps of California's Central Valley , 2018, Remote. Sens..

[58]  H. Fok,et al.  Comparison of Four Geodetic Network Densification Solutions , 2009 .

[59]  Weiping Jiang,et al.  Effect of removing the common mode errors on linear regression analysis of noise amplitudes in position time series of a regional GPS network & a case study of GPS stations in Southern California , 2018 .

[60]  M. Watkins,et al.  The gravity recovery and climate experiment: Mission overview and early results , 2004 .

[61]  David M. Allen,et al.  The Relationship Between Variable Selection and Data Agumentation and a Method for Prediction , 1974 .

[62]  F. Landerer,et al.  Accuracy of scaled GRACE terrestrial water storage estimates , 2012 .

[63]  H. Schuh,et al.  Troposphere mapping functions for GPS and very long baseline interferometry from European Centre for Medium‐Range Weather Forecasts operational analysis data , 2006 .

[64]  Y. Bock,et al.  Anatomy of apparent seasonal variations from GPS‐derived site position time series , 2001 .