Estimating signal loss in regularized GRACE gravity field solutions

SUMMARY Gravity field solutions produced using data from the Gravity Recovery and Climate Experiment (GRACE) satellite mission are subject to errors that increase as a function of increasing spatial resolution. Two commonly used techniques to improve the signal-to-noise ratio in the gravity field solutions are post-processing, via spectral filters, and regularization, which occurs within the least-squares inversion process used to create the solutions. One advantage of post-processing methods is the ability to easily estimate the signal loss resulting from the application of the spectral filter by applying the filter to synthetic gravity field coefficients derived from models of mass variation. This is a critical step in the construction of an accurate error budget. Estimating the amount of signal loss due to regularization, however, requires the execution of the full gravity field determination process to create synthetic instrument data; this leads to a significant cost in computation and expertise relative to post-processing techniques, and inhibits the rapid development of optimal regularization weighting schemes. Thus, while a number of studies have quantified the effects of spectral filtering, signal modification in regularized GRACE gravity field solutions has not yet been estimated. In this study, we examine the effect of one regularization method. First, we demonstrate that regularization can in fact be performed as a post-processing step if the solution covariance matrix is available. Regularization then is applied as a post-processing step to unconstrained solutions from the Center for Space Research (CSR), using weights reported by the Centre National d'Etudes Spatiales/Groupe de Recherches de geodesie spatiale (CNES/GRGS). After regularization, the power spectra of the CSR solutions agree well with those of the CNES/GRGS solutions. Finally, regularization is performed on synthetic gravity field solutions derived from a land surface model, revealing that in some locations significant signal loss can result from regularization. This signal loss is similar in magnitude to estimated signal loss in post-filtered solutions. End-users of GRACE data can use this method to improve the error budgets of GRACE time-series, or to restore the power lost through regularization using a scaling technique.

[1]  Jürgen Kusche,et al.  Evaluation of GRACE filter tools from a hydrological perspective , 2009 .

[2]  H. Save Using regularization for error reduction in GRACE gravity estimation , 2009 .

[3]  Peter E. Thornton,et al.  Improvements to the Community Land Model and their impact on the hydrological cycle , 2008 .

[4]  John M. Wahr,et al.  Multi‐sensor analysis of water storage variations of the Caspian Sea , 2007 .

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

[6]  S. Swenson,et al.  Accuracy of GRACE mass estimates , 2006 .

[7]  C. Shum,et al.  Non-isotropic filtering of GRACE temporal gravity for geophysical signal enhancement , 2005 .

[8]  Isabella Velicogna,et al.  Greenland mass balance from GRACE , 2005 .

[9]  Clark R. Wilson,et al.  Simulated estimation of hydrological loads from GRACE , 2005 .

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

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

[12]  S. Swenson,et al.  Estimated accuracies of regional water storage variations inferred from the Gravity Recovery and Climate Experiment (GRACE) , 2003 .

[13]  Florent Lyard,et al.  Modeling the barotropic response of the global ocean to atmospheric wind and pressure forcing ‐ comparisons with observations , 2003 .

[14]  S. Swenson,et al.  Methods for inferring regional surface‐mass anomalies from Gravity Recovery and Climate Experiment (GRACE) measurements of time‐variable gravity , 2002 .

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

[16]  J. Varah On the Numerical Solution of Ill-Conditioned Linear Systems with Applications to Ill-Posed Problems , 1973 .

[17]  Charmien Johnson Msfc Gravity Recovery and Climate Experiment , 2013 .

[18]  Jean-Charles Marty,et al.  Temporal gravity field models inferred from GRACE data , 2007 .

[19]  U. Meyer,et al.  An Earth gravity field model complete to degree and order 150 from GRACE: EIGEN-GRACE02S , 2005 .

[20]  W. Menke Geophysical data analysis : discrete inverse theory , 1984 .

[21]  J. Varah A Practical Examination of Some Numerical Methods for Linear Discrete Ill-Posed Problems , 1979 .