Using nonlinear programming to correct leakage and estimate mass change from GRACE observation and its application to Antarctica

[1] The Gravity Recovery And Climate Experiment (GRACE) mission has been providing high quality observations since its launch in 2002. Over the years, fruitful achievements have been obtained and the temporal gravity field has revealed the ongoing geophysical, hydrological and other processes. These discoveries help the scientists better understand various aspects of the Earth. However, errors exist in high degree and order spherical harmonics, which need to be processed before use. Filtering is one of the most commonly used techniques to smooth errors, yet it attenuates signals and also causes leakage of gravity signal into surrounding areas. This paper reports a new method to estimate the true mass change on the grid (expressed in equivalent water height or surface density). The mass change over the grid can be integrated to estimate regional or global mass change. This method assumes the GRACE-observed apparent mass change is only caused by the mass change on land. By comparing the computed and observed apparent mass change, the true mass change can be iteratively adjusted and estimated. The problem is solved with nonlinear programming (NLP) and yields solutions which are in good agreement with other GRACE-based estimates.

[1]  J. Famiglietti,et al.  Satellite-based estimates of groundwater depletion in India , 2009, Nature.

[2]  Benjamin F. Chao,et al.  On inversion for mass distribution from global (time-variable) gravity field , 2005 .

[3]  M. Cheng,et al.  Variations in the Earth's oblateness during the past 28 years , 2004 .

[4]  W. Tad Pfeffer,et al.  Recent contributions of glaciers and ice caps to sea level rise , 2012, Nature.

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

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

[7]  Victor Zlotnicki,et al.  Time‐variable gravity from GRACE: First results , 2004 .

[8]  Axel Rülke,et al.  On-land ice loss and glacial isostatic adjustment at the drake passage: 2003-2009 , 2011 .

[9]  Walter H. F. Smith,et al.  A global, self‐consistent, hierarchical, high‐resolution shoreline database , 1996 .

[10]  I. Sasgen,et al.  Combined GRACE and InSAR estimate of West Antarctic ice mass loss , 2010 .

[11]  Eric Rignot,et al.  Recent Antarctic ice mass loss from radar interferometry and regional climate modelling , 2008 .

[12]  I. Velicogna.  Increasing rates of ice mass loss from the Greenland and Antarctic ice sheets revealed by GRACE , 2009 .

[13]  B. Chao,et al.  An effective filtering for GRACE time‐variable gravity: Fan filter , 2009 .

[14]  Walter H. F. Smith,et al.  Free software helps map and display data , 1991 .

[15]  Coseismic gravity changes of the 2011 Tohoku‐Oki earthquake from satellite gravimetry , 2011 .

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

[17]  Byron D. Tapley,et al.  Interannual variability of Greenland ice losses from satellite gravimetry , 2011 .

[18]  Byron D. Tapley,et al.  Patagonia Icefield melting observed by Gravity Recovery and Climate Experiment (GRACE) , 2007 .

[19]  Byron D. Tapley,et al.  Accelerated Antarctic ice loss from satellite gravity measurements , 2009 .

[20]  W. Peltier.  GLOBAL GLACIAL ISOSTASY AND THE SURFACE OF THE ICE-AGE EARTH: The ICE-5G (VM2) Model and GRACE , 2004 .

[21]  W. Hager.  Analysis and implementation of a dual algorithm for constrained optimization , 1993 .

[22]  Matthew Rodell,et al.  Spatial sensitivity of the Gravity Recovery and Climate Experiment (GRACE) time‐variable gravity observations , 2005 .

[23]  Chen Ji,et al.  Crustal Dilatation Observed by GRACE After the 2004 Sumatra-Andaman Earthquake , 2006, Science.