Open access peer-reviewed chapter

Continental Water Storage Changes Sensed by GRACE Satellite Gravimetry

Written By

Guillaume Ramillien and Lucía Seoane

Submitted: October 20th, 2020 Reviewed: January 19th, 2021 Published: February 15th, 2021

DOI: 10.5772/intechopen.96109

Chapter metrics overview

436 Chapter Downloads

View Full Metrics


Since its launch in March 2002, the Gravity Recovery And Climate Experiment (GRACE) mission has been mapping the time variations of the Earth’s gravity field with a precision of 2–3 cm in terms of geoid height at the surface resolution of 300–400 km. The unprecedented precision of this twin satellite system enables to detect tiny changes of gravity that are due to the water mass variations inside the fluid envelops of our planet. Once they are corrected from known gravitational contributions of the atmosphere and the oceans, the monthly and (bi)weekly GRACE solutions reveal the continental water storage redistributions, and mainly the dominant seasonal cycle in the largest drainage river basins such as Amazon, Congo, Mississippi. The potential differences measured between the twin GRACE satellites represent the sum of integrated surface waters (lakes and rivers), soil moisture, snow, ice and groundwater. Once they are inverted for estimating surface water mass densities, GRACE solutions are also used to establish the long-term mass balance of the ice sheets impacted by global warming, for quantifying the interannual variations of the major aquifers, as well as for surveying the hydrological signatures of intense meteorological events lasting a few days such as tropical hurricanes. This chapter describes GRACE gravity products and the different data processings used for mapping continental water storage variations, it also presents the most remarkable results concerning global continental hydrology and climate changes.


  • satellite gravimetry
  • geodesy
  • global hydrology
  • gravity field
  • continental water storage

1. Introduction

Water that is present in different forms in the Earth’s system ensures the global transport of the solar heat in the oceans and atmosphere, and thus maintains life development. As it represents a precious resource, in particular for human activities, monitoring the water cycle from space remains important for its management and understanding climate change. Observation of the Earth from space, and the determination of its gravity fields in particular, provide precious information on the mass transfers in any part of the globe.

Global gravity field models are based on the theoretical expression of the variations of the geopotential V:


where λ and θ are the longitude and the latitude of the observation point respectively, r is the radial distance from the Earth’s center to the point of observation, Pnm is the associate Legendre function of degree n and order m, ae is the equatorial Earth’s radius and the gravitational parameter is the product of the gravitational constant G with the total mass of the Earth M, so that GM = 3.986004410 ∙ 1014+/−8 ∙105 m3/s2, according to IERS Standard. Space geodesy consists of determining the dimensionless Stokes coefficients Cnm and Snm of the gravity field model as precisely as possible using combined satellite data and terrestrial gravity measurements on lands. As the satellite motion depends mainly on the gravitational field according to the Newton’s law of attraction (~99% of the sensed gravity signal is from the solid Earth part), the only remote sensing technique to measure variations of water mass quantity is based on inversion of very precise satellite positions - with an accuracy of at least a few cm for detecting long wavelengths of the continental hydrology -, and/or satellite velocities [1].

Historically, long wavelengths of the gravity field time variations were determined using very precise Satellite Laser Ranging (SLR) data of 5900-km altitude LAGEOS 1–2 trajectories that reveal the movements of the center of mass of the Earth (or “geocenter”) representing a few thousands of mm, and Earth’s flatness due to seasonal mass exchange between the two hemispheres and the regular decrease due to post-glacial rebound occurring since 20 000 years [2] (Figure 1).

Figure 1.

Time variations of the C20 coefficient (representing Earth’s flatness) determined by analysis of the LAGEOS 1 & 2 satellite telemetry (source: GRGS, Toulouse).

Since the beginning of the 21st century, a new generation of passive and quasi-polar Low-Earth-Orbit (LEO) satellites has been launched to improve the spatial resolution of global gravity field models: the CHAllenging Mini-satellite Payload (CHAMP, 2000–2010) mission operated by the DLR in Germany, and the Gravity field and steady-state Ocean Circulation Explorer (GOCE, 2009–2013) of ESA. As the CHAMP mission represents its precursor, the main scientific objective of the Gravity Recovery And Climate Experiment (GRACE, 2002–2017) mission proposed by the American National Aeronautics and Space Administration (NASA) and the German Aerospace Center Deutsches Zentrum für Luft- und Raumfahrt (DLR), was to measure both static and time-varying gravity field acting in different regions of the world.

GRACE was the first mission to use the principle of two co-orbital identical satellites in pursuit, as initially proposed by [3] for estimating the spatial and temporal variations of the gravitational field which reflect mass changes in the Earth system over time scales ranging from months to a ten of years [4], so that GRACE observation represents the sum of the effects of all changes in mass which are radially integrated. In fact, GRACE observations are used to successfully survey the continental hydrology at different time scales (decanal, seasonal, rapid events) allowing to measure the climate change impacts in the Earth system, as for example, ice mass lost in Polar regions as a consequence of global warming [5].


2. GRACE mission orbit and its on-board instruments

The GRACE mission consists of two 0.5 ton satellites that followed each other at a distance of ~220 km, which were placed at a relatively low average altitude of around 450 km with a quasi-polar orbit inclination of 89.5 degrees to ensure a quasi-global coverage (Figure 2). The relative distance between the two satellites was measured with a accuracy of 1 μm.s−1 by a radar telemeter operating at K-Band microwave Ranging (KBR) [4]. The inter-satellite distance depends of the gravitational acceleration changes that affect each GRACE satellite [6]. The A three-axis accelerometer that senses the dynamical effects as non-dissipative forces (the mean solar and Earth’s radiation pressure, the atmospheric drag) is also placed on-board. Afterwards the non-gravitational effects are removed from the raw accelerations, the geopotential change along the track of the GRACE satellites is estimated from the residual perturbations in distance and inter-satellite distance changing rate. Changes of the observed inter-satellite distances reflect the variations of the Earth’s gravity field related to topography and density heterogeneities.

Figure 2.

Artistic view of the twin GRACE vehicles orbiting around the earth (source: NASA [7]).

The payload was composed of five instruments on-board, the satellite components of the GRACE mission are listed below [3]:

  • The K-band ranging system (KBR) for inter-satellite distance an accuracy of 10 μm. It uses the phases of carrier electromagnetic waves in the K and Ka bands at frequencies of 26 and 32 GHz.

  • The Ultra-Stable Oscillator (USO) for generating electromagnetic waves in the K-band for the KBR system at the desired frequency.

  • The SuperSTAR accelerometers (ACC) for accurate measurement of the forces acting on each satellite along three axis.

  • The Stellar Camera ASSEMBLY (SCA) for determining the orientation of the satellite relatively to the positions of fixed stars.

  • The Black-Jack GPS receivers and Instrument Processing Unit for providing three coordinate components of the position and the ones of the velocity of each GRACE satellite in the geocentric reference frame.


3. GRACE data products

The three official processing centers forming the GRACE Science Data Center (GSDC), i.e. the Center for Space Research (CSR) in Austin, Texas, United States; the GeoFoschungsZentrum (GFZ) in Potsdam, Germany; the Jet Propulsion Laboratory (JPL) in Pasadena, California, United States, produce the Level-1B parameter products and the Level-2 solutions derived from measurements of the GRACE mission. Level-1B products are constituted by the processed positions and velocities, which were measured by the on-board GPS receivers, accelerometers and the accurate K-band measurements of the variations in distance between the two vehicles. Using these measurements, the monthly gravity field models or Level-2 products for continental hydrology are computed. These products are expressed by means of geoid heights and Equivalent-Water Heights (EWH). The latter products are distributed by the GFZ’s Integrated System Data Center (ISDC, [8]), and the JPL’s Physical Oceanography Distributive Active Data Center (PODAAC, [9]).

3.1 Spherical harmonics solutions

The Level-2 products are obtained using a dynamic approach, which relies on the Newtonian formulation of the satellite motion equation evaluated in an inertial reference frame having the origin at the Earth’s center. The formulated solution is combined with a dedicated modeling of the gravitational and non-conservative forces, which act on the spacecrafts [6]. During the process of data reduction, the known gravitational contributions are removed from observations using a priori information from meteorological and global ocean circulation models as well as the non-gravitational forces that were measured by the on-board accelerometers [10, 11]. The residual values represent mainly the contribution of the continental hydrology and errors of the correcting models in the measured gravity field. More details about the pre-treatment of the GRACE data reduction can be found in [12]. These solutions are provided as monthly or weekly lists of Stokes coefficients, i.e. dimensionless Spherical Harmonic (SH) coefficients of the geopotential [13], up to degree and order 96 or less that correspond to a spatial resolution of 200–300 km [14, 15, 16, 17]. The range of an ideal resolution for GRACE products for hydrology is discussed in [18]. While the correcting models represent a reasonable dealiasing of high-frequency changes, the errors due to tide modeling remain in the GRACE solutions, especially for diurnal S2 tides [18, 19, 20, 21, 22]. These SH solutions are affected by north–south striping, especially dominant in the tropical band where the coverage of the satellite is insufficient mainly because of three reasons including the sparsity of GRACE track sampling in the longitudinal direction due to the polar orbit plane; propagation of systematic errors from the correcting model acceleration [19, 20, 21]; and the numerical correlations generated by solving the underdetermined systems of normal equations for the high-degree Stokes coefficients [23]. Average of each Stokes coefficient Cnm and Snm versus time is computed and removed to estimate the time anomalies for each monthly or weekly period. It is necessary to eliminate this noisy effect on the GRACE solutions, which are impacted by short North–South wavelength components. The monthly Stokes coefficients ΔCnm and ΔSnm that have been destriped by filtering have been used in very numerous studies on the evaluation of water storage on lands and the oceans [4, 24, 25, 26, 27, 28, 29, 30, 31, 32]. A simple strategy to attenuate the short-wavelength striping of the GRACE solutions consists of computing the weighted average of the solutions taken from the official centers for each period of time, and by considering the a posteriori uncertainty levels of the fitted Stokes coefficients to define these weights for combination [33].

Maps of water storage anomalies over lands directly based on [34] are:


where Δσ is the change in surface density (mass/area), ρw is the density of the water and kn is the Love numbers [35]. The change of surface mass is usually expressed in meters of equivalent water thickness (EWH). All the mass anomalies derived from GRACE data were explicated as total water mass change.

3.2 Mascons and regional solutions

Other GRACE solutions can also be obtained by other research centers where different numerical approaches are used to reach temporal resolution of one day to one month in the form of SH coefficients (global approaches) or spatial grids (local or regional approaches). As alternative to the SH approach, which is based on frequency representation instead of pure spatial localization [36], other types of base functions are used to represent surface water mass densities, mass concentration elements or mascons. In this case, water mass anomalies are estimated in specified concentrated surfaces on the Earth’s locations. The GRACE mascons have been proposed by several research groups such as Goddard Space Flight Center (GSFC) [37, 38, 39, 40, 41], Jet Propulsion Laboratory (JPL) [42, 43] and Center of Space Research (CSR) [44] at the University of Texas, Austin, where they are processed differently. As an instance, the 1° equatorial equivalent sampled mascons developed at CSR are computed by no temporal smoothing and regularization, as they are only based on GRACE information, whereas more recent mascons solutions are derived by using partial derivatives to relate KBRR observations to EWH to be determined [44]. In a second version, the mascons are related to the range rate or the range acceleration using SH that remain truncated at certain degrees and orders, as proposed by [38]. Mascons can also be estimated by post-processing of Level-2 GRACE SH solutions without a direct use of range rate observations (see examples in [44, 45, 46, 47]). Global grids of mascons solutions can be easily downloaded from [48, 49] for CSR and JPL Releases 06, respectively. Note that these latter solutions need to be scaled by a gain factor that varies geographically. A sequential Kalman Filtering (KF) approach for estimating regional maps of water mass changes by progressive integration of daily along-track GRACE geopotential anomalies has been recently proposed by [1, 50]. This iterative Kalman filter procedure has been successfully applied to determine 2° x 2° surface water mass density solutions over continental regions instead of using the SH or mascons representation [51, 52, 53].

3.3 Spatial resolution and accuracy of the GRACE products on lands

The GRACE products contribute in continental hydrology research witn a novel information: the terrestrial water storage or integrated water content, i.e. the sum of the water contained in the column from the different hydrological reservoirs: surface water, soil water, groundwater and snow cover. An early study showed an expected measurement accuracy of a few millimeters of EWH in terms of surface density for a reference water density of 1000 kg/m3, over areas of 400 km by 400 km, this work was based on Land Surface Models (LSM) outputs as soil moisture, evapotranspiration and run-off. It is expected that the presence of noise in the shorter wavelengths affects the TWS retrieval [34]. In addition, errors due to the spectra truncation increase as the area of the studied basin decrease. Based on LSM outputs and the expected accuracy of the GRACE land water solutions, it was proved that the changes in TWS could be detected by the GRACE system if they exceed 1.5 cm of EWH over an area of 200 000 km2 [54]. The accuracy of the GRACE land water solutions was expected to be about 0.7 cm of EWH for a drainage area of 400 000 km2 and 0.3 cm for a drainage basin of 4 million of km2 [55].

Current GRACE Products have a spatial resolution of a few hundred kilometers (around 200 km for the mascons and for the regional solutions, and 330 km for the releases 03 and 05 for a typical degree of truncation of n = 60–90). Errors were estimated to be around 4 cm at the Equator, and decreasing to 1.5 cm in Polar regions due to denser GRACE satellite tracks coverage [56].


4. Continental hydrology assessment by GRACE observations

The GRACE mission has been observing mass changes in the Earth’s body during 15 years, and since 2018 the new mission GRACE-FO continues this task. These satellite gravimetric missions have many applications as: to better understanding earthquake mechanism [57, 58, 59, 60]; to quantify sea level rise [61, 62] and; to observe hydrological cycles [6, 63].

Analysis of the GRACE solutions consists of finding a set of time coefficients from SH or gridded solutions, so that the TWS anomalies for a given period can be decomposed into geographical coefficients:


With ω = 2π/T and ω’ = 2π/T’, considering T = 1 year and T’ = 1/2 year, so that annual and semi-annual amplitudes C and D and phases φ and φ’ are:


These model coefficients A, B, C1, C2, D1 and D2 for GRACE variations are usually adjusted following the least square minimum criteria and then represented regionally or globally. The residuals from this relatively simple Eq. 3 of water mass variations represent intermediate wavelengths, short-term variations of unmodeled phenomena and possibly errors of the a priori correcting models (see Section 3.1).

4.1 Long term water mass variations, climate changes and irrigation

In particular, the linear trend, or equivalently the B(λ,θ) term, that is expressed in mm of EWH per year, corresponds to the increase if positive (or decrease when negative) of the water mass storage, and it can be interpreted in terms of long-term mass balance on the considered multi-year period.

A relatively complete synthesis of works on the evolution of the water storage in several parts of the world for 2002–2017 have recently presented [64] (Figure 3). The most important losses of mass are located on the ice shelves of Greenland and Antarctica, where ice storage is drastically melting at the highest rates (more than 200 Gt/y) due to the global warming and this can contribute to half of the sea level rise of about 0.3 mm/y.

Figure 3.

Global map of the TWS trends and their climatic causes sorted by colors according to [64].

In earlier studies based on SH and recently mascons solutions [47], analyzing these GRACE data has shown a continuous acceleration of the Greenland ice shelf melting.

Regions of important loss of water are revealed by GRACE such as the drought lasting up to 2007 in the southeast of Australia [48], in the North of India [62] as well as California [65]. Besides, continental waters are accumulating in other regions like in the endoheric Okavango delta in Africa [63].

The same authors have validated the constant decrease of deep water of the North Sahara aquifer sensed by GRACE with in situ water table records from wells. The Level-2 solutions need to be combined with radar altimetry data and/or model outputs in lower thus wetter latitudes, so that GRACE solutions have been used to isolate the long-term evolution of groundwater over the entire Amazon basin [66].

4.2 Seasonal cycles in major drainage river basins

The fluctuations of TWS at annual and semi-annual scales reflect the effects of climatic phenomena varying seasonally as: rainfalls, snow, temperature, evapotranspiration, river runoff, soil moisture, river discharge, groundwater and human activities. The knowledge of these periodical variations that are dominant in the large tropical basins is important to evaluate water resources. The seasonal signals estimated by GRACE are generally adjusted as described by C1, C2, D1 and D2 in Eq. 3. These amplitudes and phases determined from GRACE solutions (Eqs. 4 and 5) allow to improve the understanding of water cycle.

For example, GRACE-based annual signals over the Amazon basin show a predominant annual amplitude in fall and in spring of hundreds of EWH at the basin scale [3] driven essentially by the rainfall seasonality [67]. A comparison of GRACE observed signals and in-situ measured are shown in Figure 4. A detailed study of the major tributaries in the Amazon river basin demonstrates that the rainfall variations generated surface water fluxes delayed by two months due to transfer processes [53].

Figure 4.

Times series of in-situ observations of Agência national del Agua (red) and GRACE observations (blue) at the stations of Porto Velho and Santarem in the Amazon basin [28]. We can notice the high correlation of these dominant annual amplitudes.

In the African continent, GRACE solutions show that the stronger seasonal amplitudes are located in the Sahel latitudinal band and in the tropical Congo basin (Figure 5). An extended study based on Principal Analysis Components reveals a biannual or quadrennial water mass variations related to the West African monsoon [63]. The use of GRACE products also helps to conclude that dry season processes, in particular, evapotranspiration and the presence of vegetation, have an important role in the modeling of soil moisture in the Sahel region [66].

Figure 5.

Map of the seasonal amplitudes of the water mass changes adjusted by least squares adjustment of pure annual sinusoids at grid cells [Csinωt+φ] using the 10-day regional solutions over Africa for the period 2003–2012 [70].

GRACE solutions have also been useful to observe snow cover variations in high latitude regions [68]. For Antarctica, seasonal analysis shows very variable TWS amplitudes in coastal areas with important snow accumulation rates because of oceanic humidity [69].

In the tropical and equatorial regions, GRACE products reveal that seasonal precipitations precede land water storage with a temporal lag of 2 months, however, seasonal cycle of surface temperature is out of phase with respect to TWS, whereas in cold and temperate regions seasonal phenomena are due to more complex interactions [69]. There is also a significant contribution of river discharge in the spatial distribution of seasonal water storage with a dependency of climate [70].

Regional time variations of evapotranspiration rate - expressed in mm/yr. - can be also derived by integrating and solving the water mass balance equation, which relates TWS on land provided by GRACE, precipitation data from the Global Precipitation Climatology Centre, runoff given by a global land surface model and the unknown evapotranspiration to be determined [71].

4.3 Detection of extremes events

Gravimetric satellites missions also image the extreme climate events in the whole Earth system. Floods and droughts have been largely studied in different continents using the GRACE Level-2 products.

In the case of the Amazon basin, GRACE has revealed periods of extreme droughts and floods. During the 2005 drought, the TWS in the river and floodplains of the Amazon basin was 70% below its average for the 2003–2007 period [71]. However, in 2009 gravity measurements display an exceptional flood associated to La Niña event [26]. The maximum value of TWS in the entire Amazon basin was estimated at ∼624 ± 32 Gt with respect to the mean value.

New detection approaches based on GRACE data are developed to identify drought episodes and their severity [72]. Advances in GRACE data treatment have allowed to improve spatial and temporal resolution, then rapid extremes events of several days has been observed [73]. One example is the most powerful depressions in late summer 2005 during the cyclonic season in the Gulf of Mexico and Louisiana, better known under the name of “Katrina” (23–31/08/2005), and followed by “Rita” (17–26/09/2005). The direct consequences of the passage of these hurricanes caused important rainfalls along their tracks, thus the significant storage and accumulation of water falling on land could be observed as water mass variations in the range of a few days. As shown in Figure 6 the important rainfalls of Katrina located in the south of the Great Plains have produced river floods, and thus an important water accumulation is revealed by GRACE in the coastal region of New Orleans (up to 300 mm of EWH).

Figure 6.

Differences of weekly-averaged of GRACE solutions derived by the Kalman filter approach as proposed by [73], before and during the Katrina and Rita episodes (top); and comparison with the anomalies of TRMM precipitation for the same periods (bottom).


5. Conclusion

During its sixteen years of operation, the GRACE mission provided a novel source of information on variations of water mass on lands at unprecedented spatial and temporal resolutions. This mission offered an exceptional dataset for studying large-scale water mass redistributions, and for the very first time, the opportunity to monitor water changes in all the hydrological compartments and from regional to global scales. While spherical harmonics solutions were firstly used, regional and local approaches have already demonstrated the possible access to spatial (better localization of structures by construction) and temporal (through daily updates using Kalman filter strategies, e.g. see [73]) scales that were higher than those that were offered by global SH solutions. Additional bibliographic resources and useful information about GRACE can be found on the GRACE Tellus web site [74]. The GRACE mission ended in 2017, and later on a partnership between NASA and the German Research Centre for Geosciences (GFZ) decided to schedule GRACE Follow-On (GRACE-FO) mission to launch in 2018 in order to ensure the continuity of GRACE-type space gravimetry. By using a similar twin satellite configuration of the low-Earth and nearly polar orbit at 300–500 km altitude, the GRACE-FO mission is following its successful predecessor [75]. Additionally, it carries a demonstrator of laser system to measure the inter-satellite distance and velocity, and hence for an improved precision. It is promising for providing new perspectives in hydrology studies, such as refined long-term mass balance estimates of surface water storage and ice sheets. These new data that offer continuity with the previous GRACE observations will be of first interest for hydrology of global hydrological model calibration to constrain their operations through assimilation techniques.



We would like to thank the anonymous reviewer for his helpful contribution for improving the quality of this chapter.


  1. 1. Ramillien G, Biancale R, Gratton S, Vasseur X, Bourgogne S. GRACE-derived surface water mass anomalies by energy integral approach: application to continental hydrology. Journal of Geodesy 2011;85:313-328.
  2. 2. Dickey JO, Marcus SL, de Viron O, Fukumori I. Recent Earth oblateness variations: unraveling climate and postglacial rebound effects. Science 2002;298:1975-1977.
  3. 3. Wolff M. Direct measurements of the Earth’s gravitational potential using a satellite pair. Journal of Geophysical Research (1896-1977) 1969;74:5295-300.
  4. 4. Tapley BD, Bettadpur S, Watkins M, Reigber C. The gravity recovery and climate experiment: Mission overview and early results. Geophysical Research Letters 2004;31.
  5. 5. Velicogna I, Mohajerani Y, A G, Landerer F, Mouginot J, Noel B, et al. Continuity of Ice Sheet Mass Loss in Greenland and Antarctica From the GRACE and GRACE Follow-On Missions. Geophysical Research Letters 2020;47:e2020GL087291.
  6. 6. Schmidt R, Flechtner F, Meyer U, Neumayer K-H, Dahle Ch, König R, et al. Hydrological Signals Observed by the GRACE Satellites. Surv Geophys 2008;29:319-334.
  7. 7. GRACE photojournal JPL NASA [Internet] 2020. Available from (accessed December 18, 2020)
  8. 8. Integrated System Data Center 2018. Available from (accessed December 18, 2020)
  9. 9. Physical Oceanography Distributive Active Data Center 2020. Available from (accessed December 18, 2020)
  10. 10. Reigber C, Schmidt R, Flechtner F, König R, Meyer U, Neumayer K-H, et al. An Earth gravity field model complete to degree and order 150 from GRACE: EIGEN-GRACE02S. Journal of Geodynamics 2005;39:1-10.
  11. 11. Tapley B, Ries J, Bettadpur S, Chambers D, Cheng M, Condi F, et al. GGM02 – An improved Earth gravity field model from GRACE. J Geodesy 2005;79:467-478.
  12. 12. Ramillien G, Frappart F, Seoane L. 6 - Space Gravimetry Using GRACE Satellite Mission: Basic Concepts. In: Baghdadi N, Zribi M, editors. Microwave Remote Sensing of Land Surface, Elsevier; 2016, p. 285-302.
  13. 13. Hofmann-Wellenhof B, Moritz H. Physical Geodesy. 2nd ed. Wien: Springer-Verlag; 2006.
  14. 14. Bettadpur S. UTCSR level-2 processing standards document for level-2 product release 0004, GRACE 327-742 (CSR-GR-03-03). Center for Space Research, The University of Texas at Austin, Austin. 2007
  15. 15. Flechtner F. GFZ Level-2 Processing Stan -dards Document for Product Release 0004. Rev.1.0, GRACE project documentation 327-743 2007
  16. 16. Chambers DP, Bonin JA. Evaluation of Release-05 GRACE time-variable gravity coefficients over the ocean. Ocean Sci 2012;8:859-868.
  17. 17. Dahle C, Flechtner F, Gruber C, König D, König R, Michalak G, et al. GFZ GRACE Level-2 Processing Standards Document for Level-2 Product Release 0005 : revised edition, January 2013 2013
  18. 18. Vishwakarma DB, Devaraju B, Sneeuw N. What Is the Spatial Resolution of grace Satellite Products for Hydrology? Remote Sensing 2018;10.
  19. 19. Forootan E, Didova O, Schumacher M, Kusche J, Elsaka B. Comparisons of atmospheric mass variations derived from ECMWF reanalysis and operational fields, over 2003-2011. Journal of Geodesy 2014;88:503-514.
  20. 20. Han S-C, Jekeli C, Shum CK. Time-variable aliasing effects of ocean tides, atmosphere, and continental water mass on monthly mean GRACE gravity field: TEMPORAL ALIASING ON GRACE GRAVITY FIELD. J Geophys Res 2004;109.
  21. 21. Ray RD, Luthcke SB. Tide model errors and GRACE gravimetry: towards a more realistic assessment. Geophysical Journal International 2006;167:1055-1059.
  22. 22. Thompson PF, Bettadpur SV, Tapley BD. Impact of short period, non-tidal, temporal mass variability on GRACE gravity estimates: IMPACT OF SHORT PERIOD, NON-TIDAL, TEMPORAL MASS VARIABILITY ON GRACE GRAVITY ESTIMATES. Geophys Res Lett 2004;31:n/a-n/a.
  23. 23. Swenson S, Wahr J. Post-processing removal of correlated errors in GRACE data. Geophysical Research Letters 2006;33.
  24. 24. Swenson S, Wahr J. Monitoring the water balance of Lake Victoria, East Africa, from space. Journal of Hydrology 2009;370:163-176.
  25. 25. Chambers DP, Wahr J, Nerem RS. Preliminary observations of global ocean mass variations with GRACE. Geophysical Research Letters 2004;31.
  26. 26. Chen JL, Wilson CR, Tapley BD, Yang ZL, Niu GY. 2005 drought event in the Amazon River basin as measured by GRACE and estimated by climate models. Journal of Geophysical Research: Solid Earth 2009;114.
  27. 27. Ramillien G, Frappart F, Cazenave A, Güntner A. Time variations of land water storage from an inversion of 2 years of GRACE geoids. Earth and Planetary Science Letters 2005;235:283-301.
  28. 28. Ramillien G, Famiglietti JS, Wahr J. Detection of Continental Hydrology and Glaciology Signals from GRACE: A Review. Surv Geophys 2008;29:361-374.
  29. 29. Swenson SC, Milly PCD. Climate model biases in seasonality of continental water storage revealed by satellite gravimetry. Water Resources Research 2006;42.
  30. 30. Tamisiea ME, Leuliette EW, Davis JL, Mitrovica JX. Constraining hydrological and cryospheric mass flux in southeastern Alaska using space-based gravity measurements. Geophysical Research Letters 2005;32.
  31. 31. Velicogna I, Wahr J. Greenland mass balance from GRACE. Geophysical Research Letters 2005;32.
  32. 32. Wahr J, Swenson S, Zlotnicki V, Velicogna I. Time-variable gravity from GRACE: First results. Geophysical Research Letters 2004;31.
  33. 33. Sakumura C, Bettadpur S, Bruinsma S. Ensemble prediction and intercomparison analysis of GRACE time-variable gravity field models. Geophysical Research Letters 2014;41:1389-1397.
  34. 34. Wahr J, Molenaar M, Bryan F. Time variability of the Earth’s gravity field: Hydrological and oceanic effects and their possible detection using GRACE. Journal of Geophysical Research: Solid Earth 1998;103:30205-30229.
  35. 35. Han D, Wahr J. The viscoelastic relaxation of a realistically stratified earth, and a further analysis of postglacial rebound. Geophysical Journal International 1995;120:287-311.
  36. 36. Freeden W, Schreiner M. Spherical Functions of Mathematical Geosciences: A Scalar, Vectorial, and Tensorial Setup. Berlin, Heidelberg: Springer Berlin Heidelberg; 2009.
  37. 37. Luthcke SB, Zwally HJ, Abdalati W, Rowlands DD, Ray RD, Nerem RS, et al. Recent Greenland Ice Mass Loss by Drainage System from Satellite Gravity Observations. Science 2006;314:1286-1289.
  38. 38. Luthcke SB, Sabaka TJ, Loomis BD, Arendt AA, McCarthy JJ, Camp J. Antarctica, Greenland and Gulf of Alaska land-ice evolution from an iterated GRACE global mascon solution. Journal of Glaciology 2013;59:613-631.
  39. 39. Rowlands DD, Luthcke SB, Klosko SM, Lemoine FGR, Chinn DS, McCarthy JJ, et al. Resolving mass flux at high spatial and temporal resolution using GRACE intersatellite measurements. Geophysical Research Letters 2005;32.
  40. 40. Rowlands DD, Luthcke SB, McCarthy JJ, Klosko SM, Chinn DS, Lemoine FG, et al. Global mass flux solutions from GRACE: A comparison of parameter estimation strategies—Mass concentrations versus Stokes coefficients. Journal of Geophysical Research: Solid Earth 2010;115
  41. 41. Watkins MM, Wiese DN, Yuan D-N, Boening C, Landerer FW. Improved methods for observing Earth’s time variable mass distribution with GRACE using spherical cap mascons. Journal of Geophysical Research: Solid Earth 2015;120:2648-2671.
  42. 42. Wiese DN, Landerer FW, Watkins MM. Quantifying and reducing leakage errors in the JPL RL05M GRACE mascon solution. Water Resources Research 2016;52:7490-7502.
  43. 43. Wiese DN, Yuan D-N, Boening C, Landerer FW, Watkins MM. JPL GRACE Mascon Ocean, Ice, and Hydrology Equivalent Water Height Release 06 Coastal Resolution Improvement (CRI) Filtered Version 1.0 2018.
  44. 44. Save H, Bettadpur S, Tapley BD. High-resolution CSR GRACE RL05 mascons. Journal of Geophysical Research: Solid Earth 2016;121:7547-7569.
  45. 45. Jacob T, Wahr J, Pfeffer WT, Swenson S. Recent contributions of glaciers and ice caps to sea level rise. Nature 2012;482:514-518.
  46. 46. Schrama EJO, Wouters B, Rietbroek R. A mascon approach to assess ice sheet and glacier mass balances and their uncertainties from GRACE data. Journal of Geophysical Research: Solid Earth 2014;119:6048-6066.
  47. 47. Velicogna I, Sutterley TC, van den Broeke MR. Regional acceleration in ice mass loss from Greenland and Antarctica using GRACE time-variable gravity data. Geophysical Research Letters 2014;41:8130-8137.
  48. 48. CSR GRACE/GRACE-FO RL06 Mascon Solutions [Internet] 2020. Available from (accessed December 3, 2020)
  49. 49. Monthly Mass Grids - Global mascons (JPL RL06_v02) [Internet] 2020. Available from (accessed December 3, 2020)
  50. 50. Ramillien GL, Frappart F, Gratton S, Vasseur X. Sequential estimation of surface water mass changes from daily satellite gravimetry data. Journal of Geodesy 2015;89:259-282.
  51. 51. Ramillien GL, Seoane L, Frappart F, Biancale R, Gratton S, Vasseur X, et al. Constrained Regional Recovery of Continental Water Mass Time-variations from GRACE-based Geopotential Anomalies over South America. Surveys in Geophysics 2012;33:887-905.
  52. 52. Seoane L, Ramillien G, Frappart F, Leblanc M. Regional GRACE-based estimates of water mass variations over Australia: validation and interpretation. Hydrology and Earth System Sciences 2013;17:4925-4939.
  53. 53. Frappart F, Seoane L, Ramillien G. Validation of GRACE-derived terrestrial water storage from a regional approach over South America. Remote Sensing of Environment 2013;137:69-83.
  54. 54. Rodell M, Famiglietti JS. Detectability of variations in continental water storage from satellite observations of the time dependent gravity field. Water Resources Research 1999;35:2705-2723.
  55. 55. Swenson S, Wahr J, Milly PCD. Estimated accuracies of regional water storage variations inferred from the Gravity Recovery and Climate Experiment (GRACE). Water Resources Research 2003;39.
  56. 56. Landerer FW, Swenson SC. Accuracy of scaled GRACE terrestrial water storage estimates. Water Resources Research 2012;48:n/a-n/a.
  57. 57. Han S-C, Shum CK, Bevis M, Ji C, Kuo C-Y. Crustal dilatation observed by GRACE after the 2004 Sumatra-Andaman earthquake. Science 2006;313:658-662.
  58. 58. Han S-C, Sauber J, Luthcke S. Regional gravity decrease after the 2010 Maule (Chile) earthquake indicates large-scale mass redistribution. Geophysical Research Letters 2010;37.
  59. 59. Li J, Shen W-B. Monthly GRACE detection of coseismic gravity change associated with 2011 Tohoku-Oki earthquake using northern gradient approach. Earth, Planets and Space 2015;67:29.
  60. 60. Fatolazadeh F, Goïta K, Javadi Azar R. Determination of earthquake epicentres based upon invariant quantities of GRACE strain gravity tensors. Scientific Reports 2020;10:7636.
  61. 61. Hsu C-W, Velicogna I. Detection of sea level fingerprints derived from GRACE gravity data. Geophysical Research Letters 2017;44:8953-8961.
  62. 62. Jeon T, Seo K-W, Youm K, Chen J, Wilson CR. Global sea level change signatures observed by GRACE satellite gravimetry. Scientific Reports 2018;8:13519.
  63. 63. Frappart F, Ramillien G, Leblanc M, Tweed SO, Bonnet M-P, Maisongrande P. An independent component analysis filtering approach for estimating continental hydrology in the GRACE gravity data. Remote Sensing of Environment 2011;115:187-204.
  64. 64. Rodell M, Famiglietti JS, Wiese DN, Reager JT, Beaudoing HK, Landerer FW, et al. Emerging trends in global freshwater availability. Nature 2018;557:651-659.
  65. 65. Thomas BF, Famiglietti JS, Landerer FW, Wiese DN, Molotch NP, Argus DF. GRACE Groundwater Drought Index: Evaluation of California Central Valley groundwater drought. Remote Sensing of Environment 2017;198:384-392.
  66. 66. Frappart F, Papa F, Güntner A, Tomasella J, Pfeffer J, Ramillien G, et al. The spatio-temporal variability of groundwater storage in the Amazon River Basin. Advances in Water Resources 2019;124:41-52.
  67. 67. Crowley JW, Mitrovica JX, Bailey RC, Tamisiea ME, Davis JL. Annual variations in water storage and precipitation in the Amazon Basin. J Geod 2008;82:9-13.
  68. 68. Frappart F, Ramillien G, Seoane L. 8 - Monitoring Water Mass Redistributions on Land and Polar Ice Sheets Using the GRACE Gravimetry from Space Mission. In: Baghdadi N, Zribi M, editors. Land Surface Remote Sensing in Continental Hydrology, Elsevier; 2016, p. 255-279.
  69. 69. Llubes M, Lemoine J-M, Rémy F. Antarctica seasonal mass variations detected by GRACE. Earth and Planetary Science Letters 2007;260:127-136.
  70. 70. Ramillien G, Frappart F, Seoane L. Application of the Regional Water Mass Variations from GRACE Satellite Gravimetry to Large-Scale Water Management in Africa. Remote Sensing 2014;6:7379-7405.
  71. 71. Frappart F, Papa F, Silva JS da, Ramillien G, Prigent C, Seyler F, et al. Surface freshwater storage and dynamics in the Amazon basin during the 2005 exceptional drought. Environ Res Lett 2012;7:044010.
  72. 72. Thomas AC, Reager JT, Famiglietti JS, Rodell M. A GRACE-based water storage deficit approach for hydrological drought characterization. Geophysical Research Letters 2014;41:1537-1545.
  73. 73. Ramillien G, Seoane L, Schumacher M, Forootan E, Frappart F, Darrozes J. Recovery of Rapid Water Mass Changes (RWMC) by Kalman Filtering of GRACE Observations. Remote Sensing 2020;12:1299.
  74. 74. GRACE Tellus [Internet] 2020. Available from (accessed December 18, 2020)
  75. 75. GRACE FO [Internet] 2020. Available from (accessed January 18, 2021)

Written By

Guillaume Ramillien and Lucía Seoane

Submitted: October 20th, 2020 Reviewed: January 19th, 2021 Published: February 15th, 2021