GRAMAT: a comprehensive Matlab toolbox for estimating global mass variations from GRACE satellite data

In this paper, we robustly analyze the noise reduction methods for processing spherical harmonic (SH) coefficient data products collected by the Gravity Recovery and Climate Experiment (GRACE) satellite mission and devise a comprehensive GRACE Matlab Toolbox (GRAMAT) to estimate spatio-temporal mass variations over land and oceans. Functions in GRAMAT contain: (1) destriping of SH coefficients to remove “north-to-south” stripes, or geographically correlated high-frequency errors, and Gaussian smoothing, (2) spherical harmonic analysis and synthesis, (3) assessment and reduction of the leakage effect in GRACE-derived mass variations, and (4) harmonic analysis of regional time series of the mass variations and assessment of the uncertainty of the GRACE estimates. As a case study, we analyze the terrestrial water storage (TWS) variations in the Amazon River basin using the functions in GRAMAT. In addition to obvious seasonal TWS variations in the Amazon River basin, significant interannual TWS variations are detected by GRACE using the GRAMAT, which are consistent with precipitation anomalies in the region. We conclude that using GRAMAT and processing the GRACE level-2 data products, the global spatio-temporal mass variations can be efficiently and robustly estimated, which indicates the potential wide range of GRAMAT’s applications in hydrology, oceanography, cryosphere, solid Earth and geophysical disciplines to interpret large-scale mass redistribution and transport in the Earth system. We postulate that GRAMAT will also be an effective tool for the analysis of data from the upcoming GRACE-Follow-On mission.

[1]  C. Jekeli.  Alternative methods to smooth the Earth's gravity field , 1981 .

[2]  R. Koster,et al.  Modeling the land surface boundary in climate models as a composite of independent vegetation stands , 1992 .

[3]  D. Lettenmaier,et al.  A simple hydrologically based model of land surface water and energy fluxes for general circulation models , 1994 .

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

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

[6]  J. D. Tarpley,et al.  Implementation of Noah land surface model advances in the National Centers for Environmental Prediction operational mesoscale Eta model , 2003 .

[7]  M. Watkins,et al.  GRACE Measurements of Mass Variability in the Earth System , 2004, Science.

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

[9]  C. K. Shum,et al.  Regional high‐resolution spatiotemporal gravity modeling from GRACE data using spherical wavelets , 2006 .

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

[11]  Hubert H. G. Savenije,et al.  The bias in GRACE estimates of continental water storage variations , 2006 .

[12]  J. Wahr,et al.  Acceleration of Greenland ice mass loss in spring 2004 , 2006, Nature.

[13]  J. Wahr,et al.  Measurements of Time-Variable Gravity Show Mass Loss in Antarctica , 2006, Science.

[14]  D. Chambers.  Evaluation of new GRACE time‐variable gravity data over the ocean , 2006 .

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

[16]  D. García,et al.  On the steric and mass‐induced contributions to the annual sea level variations in the Mediterranean Sea , 2006 .

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

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

[19]  I. Sasgen,et al.  Wiener optimal combination and evaluation of the Gravity Recovery and Climate Experiment (GRACE) gravity fields over Antarctica , 2007 .

[20]  J. Kusche.  Approximate decorrelation and non-isotropic smoothing of time-variable GRACE-type gravity field models , 2007 .

[21]  Matthias Holschneider,et al.  Coseismic and post-seismic signatures of the Sumatra 2004 December and 2005 March earthquakes in GRACE satellite gravity , 2007 .

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

[23]  Byron D. Tapley,et al.  GRACE detects coseismic and postseismic deformation from the Sumatra‐Andaman earthquake , 2007 .

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

[25]  J. Willis,et al.  Assessing the globally averaged sea level budget on seasonal to interannual timescales , 2008 .

[26]  D. Chambers,et al.  Estimating Geocenter Variations from a Combination of GRACE and Ocean Model Output , 2008 .

[27]  M. Tamisiea,et al.  A statistical filtering approach for Gravity Recovery and Climate Experiment (GRACE) gravity data , 2008 .

[28]  Byron D. Tapley,et al.  Antarctic regional ice loss rates from GRACE , 2008 .

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

[30]  V. M. Tiwari,et al.  Dwindling groundwater resources in northern India, from satellite gravity observations , 2009 .

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

[32]  Guillaume Ramillien,et al.  Sea level budget over 2003-2008: A reevaluation from GRACE space gravimetry, satellite altimetry and Argo , 2009 .

[33]  C. Shum,et al.  On the postprocessing removal of correlated errors in GRACE temporal gravity field solutions , 2009 .

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

[35]  B. Tapley,et al.  2005 drought event in the Amazon River basin as measured by GRACE and estimated by climate models , 2009 .

[36]  B. D. Tapley,et al.  Recent La Plata basin drought conditions observed by satellite gravimetry , 2010 .

[37]  Byron D. Tapley,et al.  The 2009 exceptional Amazon flood and interannual terrestrial water storage change observed by GRACE , 2010 .

[38]  Koji Matsuo,et al.  Time-variable ice loss in Asian high mountains from satellite gravimetry , 2010 .

[39]  B. Scanlon,et al.  GRACE Hydrological estimates for small basins: Evaluating processing approaches on the High Plains Aquifer, USA , 2010 .

[40]  Minkang Cheng,et al.  Variations of the Earth's figure axis from satellite laser ranging and GRACE , 2011 .

[41]  S. Swenson,et al.  Satellites measure recent rates of groundwater depletion in California's Central Valley , 2011 .

[42]  F. Simons,et al.  Mapping Greenland’s mass loss in space and time , 2012, Proceedings of the National Academy of Sciences.

[43]  B. Scanlon,et al.  Ground referencing GRACE satellite estimates of groundwater storage changes in the California Central Valley, USA , 2012 .

[44]  C. Prigent,et al.  Surface freshwater storage and dynamics in the Amazon basin during the 2005 exceptional drought , 2012 .

[45]  Matt A. King,et al.  Lower satellite-gravimetry estimates of Antarctic sea-level contribution , 2012, Nature.

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

[47]  Don P. Chambers,et al.  Evaluation of Release-05 GRACE time-variable gravity coefficients over the ocean , 2012 .

[48]  R. Steven Nerem,et al.  The 2011 La Niña: So strong, the oceans fell , 2012 .

[49]  J. Wahr,et al.  Computations of the viscoelastic response of a 3-D compressible Earth to surface loading: an application to Glacial Isostatic Adjustment in Antarctica and Canada , 2012 .

[50]  U. Schneider,et al.  GPCC's new land surface precipitation climatology based on quality-controlled in situ data and its role in quantifying the global water cycle , 2013, Theoretical and Applied Climatology.

[51]  W. Feng,et al.  Evaluation of groundwater depletion in North China using the Gravity Recovery and Climate Experiment (GRACE) data and ground‐based measurements , 2013 .

[52]  Wenke Sun,et al.  Evaluation of glacier changes in high‐mountain Asia based on 10 year GRACE RL05 models , 2013 .

[53]  M. Cheng,et al.  Deceleration in the Earth's oblateness , 2013 .

[54]  W. Feng.  Regional terrestrial water storage and sea level variations inferred from satellite gravimetry , 2014 .

[55]  S. Swenson,et al.  Estimating the human contribution to groundwater depletion in the Middle East, from GRACE data, land surface models, and well observations , 2014 .

[56]  W. Feng,et al.  Mass-induced sea level variations in the Red Sea from GRACE, steric-corrected altimetry, in situ bottom pressure records, and hydrographic observations , 2014 .

[57]  R. Rietbroek,et al.  Key Points: @bullet Consistent Method for Estimating Mass Balances from Grace @bullet Mascon Technique @bullet Evaluate Systematic Errors Gia Correction a Mascon Approach to Assess Ice Sheet and Glacier Mass Balances and Their Uncertainties from Grace Data , 2022 .

[58]  B. Scanlon,et al.  Global analysis of approaches for deriving total water storage changes from GRACE satellites , 2015 .

[59]  J. Kusche,et al.  Sea level budget in the Bay of Bengal (2002–2014) from GRACE and altimetry , 2016 .

[60]  Y. Hong,et al.  Have GRACE satellites overestimated groundwater depletion in the Northwest India Aquifer? , 2016, Scientific Reports.