Statistically optimal estimation of Greenland Ice Sheet mass variations from GRACE monthly solutions using an improved mascon approach

We present an improved mascon approach to transform monthly spherical harmonic solutions based on GRACE satellite data into mass anomaly estimates in Greenland. The GRACE-based spherical harmonic coefficients are used to synthesize gravity anomalies at satellite altitude, which are then inverted into mass anomalies per mascon. The limited spectral content of the gravity anomalies is properly accounted for by applying a low-pass filter as part of the inversion procedure to make the functional model spectrally consistent with the data. The full error covariance matrices of the monthly GRACE solutions are properly propagated using the law of covariance propagation. Using numerical experiments, we demonstrate the importance of a proper data weighting and of the spectral consistency between functional model and data. The developed methodology is applied to process real GRACE level-2 data (CSR RL05). The obtained mass anomaly estimates are integrated over five drainage systems, as well as over entire Greenland. We find that the statistically optimal data weighting reduces random noise by 35–69%, depending on the drainage system. The obtained mass anomaly time-series are de-trended to eliminate the contribution of ice discharge and are compared with de-trended surface mass balance (SMB) time-series computed with the Regional Atmospheric Climate Model (RACMO 2.3). We show that when using a statistically optimal data weighting in GRACE data processing, the discrepancies between GRACE-based estimates of SMB and modelled SMB are reduced by 24–47%.

[1]  Álvaro González.  Measurement of Areas on a Sphere Using Fibonacci and Latitude–Longitude Lattices , 2009, 0912.4540.

[2]  Srinivas Bettadpur,et al.  Impact of short period, non‐tidal, temporal mass variability on GRACE gravity estimates , 2004 .

[3]  Bob E. Schutz,et al.  Comparison of Elevation Change Detection Methods From ICESat Altimetry Over the Greenland Ice Sheet , 2017, IEEE Transactions on Geoscience and Remote Sensing.

[4]  Richard Biancale,et al.  Separation of coseismic and postseismic gravity changes for the 2004 Sumatra–Andaman earthquake from 4.6 yr of GRACE observations and modelling of the coseismic change by normal-modes summation , 2009 .

[5]  Duncan J. Wingham,et al.  A method of combining ICESat and GRACE satellite data to constrain Antarctic mass balance , 2000 .

[6]  Nico Sneeuw,et al.  Assessing Greenland ice mass loss by means of point-mass modeling: a viable methodology , 2011 .

[7]  Srinivas Bettadpur,et al.  High‐resolution CSR GRACE RL05 mascons , 2016 .

[8]  Grzegorz Michalak,et al.  GFZ GRACE Level-2 Processing Standards Document for Level-2 Product Release 0005 : revised edition, January 2013 , 2013 .

[9]  Guochang Xu,et al.  Sciences of geodesy , 2010 .

[10]  Michael Durand,et al.  Comparing satellite derived precipitation datasets using the Hillslope River Routing (HRR) model in the Congo River Basin , 2011 .

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

[12]  Louise Sørensen,et al.  Greenland and Antarctica Ice Sheet Mass Changes and Effects on Global Sea Level , 2017, Surveys in Geophysics.

[13]  R. Nerem,et al.  Recent Greenland Ice Mass Loss by Drainage System from Satellite Gravity Observations , 2006, Science.

[14]  O. Baur.  Greenland mass variation from time‐variable gravity in the absence of GRACE , 2013 .

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

[16]  M. R. van den Broeke,et al.  Partitioning Recent Greenland Mass Loss , 2009, Science.

[17]  E. Schrama,et al.  Revisiting Greenland ice sheet mass loss observed by GRACE , 2011 .

[18]  G. Ramillien,et al.  Global time-variations of hydrological signals from GRACE satellite gravimetry , 2004 .

[19]  Frank Flechtner,et al.  Simulating high‐frequency atmosphere‐ocean mass variability for dealiasing of satellite gravity observations: AOD1B RL05 , 2013 .

[20]  Pavel Ditmar,et al.  Understanding data noise in gravity field recovery on the basis of inter-satellite ranging measurements acquired by the satellite gravimetry mission GRACE , 2012, Journal of Geodesy.

[21]  D. Chambers,et al.  Uncertainty estimates of a GRACE inversion modelling technique over Greenland using a simulation , 2013 .

[22]  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 .

[23]  Guochang Xu.  Sciences of Geodesy - I: Advances and Future Directions , 2014 .

[24]  J. Camp,et al.  Antarctica, Greenland and Gulf of Alaska land-ice evolution from an iterated GRACE global mascon solution , 2013, Journal of Glaciology.

[25]  B. Tapley,et al.  The new GRACE gravity mission and its value to exploration , 2004 .

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

[27]  W. Bosch,et al.  EOT11A - Empirical Ocean Tide Model from Multi-Mission Satellite Altimetry , 2008 .

[28]  M. Watkins,et al.  Quantifying and reducing leakage errors in the JPL RL05M GRACE mascon solution , 2016 .

[29]  C. Jekeli.  Modifying Stokes' function to reduce the error of geoid undulation computations , 1981 .

[30]  R. S. W. van de Wal,et al.  Evaluation of the updated regional climate model RACMO2.3: summer snowfall impact on the Greenland Ice Sheet , 2015 .

[31]  Jery R. Stedinger,et al.  Regional Hydrologic Analysis, 2, Model‐Error Estimators, Estimation of Sigma and Log‐Pearson Type 3 Distributions , 1986 .

[32]  J. Kusche,et al.  Regularization of geopotential determination from satellite data by variance components , 2002 .

[33]  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 .

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

[35]  D. Chambers,et al.  GRACE observes small‐scale mass loss in Greenland , 2008 .

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

[37]  Matthew Rodell,et al.  Low degree spherical harmonic influences on Gravity Recovery and Climate Experiment (GRACE) water storage estimates , 2005 .

[38]  M. Watkins,et al.  Improved methods for observing Earth's time variable mass distribution with GRACE using spherical cap mascons , 2015 .

[39]  Xianglin Liu,et al.  Estimation of mass change trends in the Earth’s system on the basis of GRACE satellite data, with application to Greenland , 2012, Journal of Geodesy.

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

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

[42]  Isabella Velicogna,et al.  Regional acceleration in ice mass loss from Greenland and Antarctica using GRACE time‐variable gravity data , 2014 .