Applications and Challenges of GRACE and GRACE Follow-On Satellite Gravimetry

Time-variable gravity measurements from the Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-On (GRACE-FO) missions have opened up a new avenue of opportunities for studying large-scale mass redistribution and transport in the Earth system. Over the past 19 years, GRACE/GRACE-FO time-variable gravity measurements have been widely used to study mass variations in different components of the Earth system, including the hydrosphere, ocean, cryosphere, and solid Earth, and significantly improved our understanding of long-term variability of the climate system. We carry out a comprehensive review of GRACE/GRACE-FO satellite gravimetry, time-variable gravity fields, data processing methods, and major applications in several different fields, including terrestrial water storage change, global ocean mass variation, ice sheets and glaciers mass balance, and deformation of the solid Earth. We discuss in detail several major challenges we need to face when using GRACE/GRACE-FO time-variable gravity measurements to study mass changes, and how we should address them. We also discuss the potential of satellite gravimetry in detecting gravitational changes that are believed to originate from the deep Earth. The extended record of GRACE/GRACE-FO gravity series, with expected continuous improvements in the coming years, will lead to a broader range of applications and improve our understanding of both climate change and the Earth system.

Regional Seafloor Topography by Extended Kalman Filtering of Marine Gravity Data without Ship-Track Information

An iterative Extended Kalman Filter (EKF) approach is proposed to recover a regional set of topographic heights composing an undersea volcanic mount by the successive combination of large numbers of gravity measurements at sea surface using altimetry satellite-derived grids and taking the error uncertainties into account. The integration of the non-linear Newtonian operators versus the radial and angular distances (and its first derivatives) enables the estimation process to accelerate and requires only few iterations, instead of summing Legendre polynomial series or using noise-degraded 2D-FFT decomposition. To show the effectiveness of the EKF approach, we apply it to the real case of the bathymetry around the Great Meteor seamount in the Atlantic Ocean by combining only geoid height/free-air anomaly datasets and using ship-track soundings as reference for validation. Topography of the Great Meteor seamounts structures are well-reconstructed, especially when regional compensation is considered. Best solution gives a RMS equal to 400 m with respect to the single beam depth observations and it is comparable to RMS obtained for ETOPO1 of about 365 m. Larger discrepancies are located in the seamount flanks due to missing high-resolution information for gradients. This approach can improve the knowledge of seafloor topography in regions where few echo-sounder measurements are available.

Comparison of Marine Gravity Measurements from Shipborne and Satellite Altimetry in the Arctic Ocean

To understand the influence of sea ice on shipborne gravity measurements and the accuracy of the satellite-altimetry-derived gravity field in the Arctic Ocean, we compared shipborne gravity measurements with those obtained from satellite altimetric gravity measurements. The influence of sea ice on the shipborne gravity measurements was mainly concentrated in the 0–6 km wavelength range, and the standard deviation of the noise amplitudes was 2.62 mGal. Compared to ice-free regions, the accuracies in the region with floating ice were reduced by 13% for DTU21 and 6% for SV31. Due to the influence of sea ice, satellite altimetric gravity data lose significant information in the 9–12 km wavelength range. The coherence curve of the shipborne gravity with bathymetry was nearly the same as that of the satellite altimetric gravity. The satellite data contain nearly all of the significant information that is present in the shipborne data. The differences between the shipborne and satellite gravity data are small and can be used to study the crustal structure of the Arctic.

Gravity Survey of King Fahd University of Petroleum and Minerals Dammam Dome, Saudi Arabia

The study area is a part of Dammam Dome that is situated at King Fahd University of Petroleum & Minerals (KFUPM) campus, Dhahran, Kingdom of Saudi Arabia. The gravity survey was conducted as a pilot case study to explore part of Dammam Dome in greater detail. Gravity data were collected solely during night hours due to low noise levels. A significant part of the survey was conducted during the summer holiday period, , when there was no student are on campus. A total of 235 gravity measurements were made using a Scintrex CG5 gravitometer, while a Trimble R10+ differential GPS (DGPS) was used to measure the stations’ location and elevation with the highest accuracy. All gravity data were reduced using several algorithms, and their outcomes were cross-compared. The Complete Bouguer anomaly map for the campus was then generated. Several enhancement filters including edged detection and shallow to deeper source separation were applied. Data were inverted, and 2.5D and 3D models were created to image the subsurface conditions. The main purpose of this study is to better understand the subsurface geology, tectonic settings of the Dammam Dome by applying the high-resolution gravity method before carrying out any comprehensive geophysical (seismic) 3D survey.

GRACE Follow-On revealed Bangladesh was flooded early in the 2020 monsoon season due to premature soil saturation

The overall size and timing of monsoon floods in Bangladesh are challenging to measure. The inundated area is extensive in low-lying Bangladesh, and observations of water storage are key to understanding floods. Laser-ranging instruments on Gravity Recovery and Climate Experiment (GRACE) Follow-On spacecraft detected the peak water storage anomaly of 75 gigatons across Bangladesh in late July 2020. This is in addition to, and three times larger than, the maximum storage anomaly in soil layers during the same period. A flood propagation model suggested that the water mass, as shown in satellite observations, is largely influenced by slow floodplain and groundwater flow processes. Independent global positioning system measurements confirmed the timing and total volume of the flood water estimates. According to land surface models, the soils were saturated a month earlier than the timing of the peak floodplain storage observed by GRACE Follow-On. The cyclone Amphan replenished soils with rainfall just before the monsoon rains started, and consequently, excessive runoff was produced and led to the early onset of the 2020 flooding. This study demonstrated how antecedent soil moisture conditions can influence the magnitude and duration of flooding. Continuous monitoring of storage change from GRACE Follow-On gravity measurements provides important information complementary to river gauges and well levels for enhancing hydrologic flood forecasting models and assisting surface water management.

ADAPTATION OF THE GLOBAL GEOID MODEL EGM2008 ON CAMPANIA REGION (ITALY) BASED ON GEODETIC NETWORK POINTS

The knowledge of the geoid undulation, the height of the geoid relative to a given ellipsoid of reference, is fundamental to transform the ellipsoidal heights into orthometric heights. Global geoid undulation models developed from satellite gravity measurements appropriately integrated with other data, are free accessible in internet, but their accuracy may be inadequate for specific applications. Earth Gravitational Model 2008 (EGM2008) is one of those: usually available in grid form 2.5’ × 2.5’ (a geotif is developed by Agisoft with resolution 1’ × 1’), it defines the difference between the WGS84 ellipsoid height and the mean sea level, but in some areas the discrepancies between these geoid undulations and local correspondent measured values are on the order of various decimetres. For consequence, more accurate models are necessary. This article aims to determine a geoid undulation model suitable for Campania Region (Italy), starting from the global model EGM2008 (1’ × 1’) that is locally adjusted by using geodetic network points (GNPs) and GIS interpolation functions. Three different datasets are considered including respectively 20, 40 and 60 GNPs and three deterministic interpolators are applied in global way to generate geoid undulation grids: Inverse Distance Weight (IDW), Global Polynomial 1st order (GP1), Global Polynomial 2nd order (GP2). The resultant 9 models are tested on 20 additional GNPs. The experiments demonstrate that local geoid can be produced on a little area adapting global geoid by means of GNPs: the model obtained using GP2 and 60 GNPs, the most accurate one, fits the data with ±3.2 cm root mean square error (RMSE).

Analysis of the interannual variability in satellite gravity solutions: detection of climate modes fingerprints in water mass displacements across continents and oceans

This study analyzes the interannual variability of the water mass transport measured by satellite gravity missions in regard to eight major climate modes known to influence the Earth’s climate from regional to global scales. Using sparsity promoting techniques (i.e., LASSO), we automatically select the most relevant predictors of the climate variability among the eight candidates considered. The El Niño–Southern Oscillation, Southern Annular Mode and Arctic Oscillation are shown to account for a large part the interannual variability of the water mass transport observed in extratropical ocean basins (up to 40%) and shallow seas (up to 70%). A combination of three Pacific and one Atlantic modes is needed to account for most (up to 60%) of the interannual variability of the terrestrial water storage observed in the North Amazon, Parana and Zambezi basins. With our technique, the impact of climate modes on water mass changes can be tracked across distinct water reservoirs (oceans, continents and ice-covered regions) and we show that a combination of climate modes is necessary to explain at best the natural variability in water mass transport. The climate modes predictions based on LASSO inversions can be used to reduce the inter-annual variability in satellite gravity measurements and detect processes unrelated with the natural variability of climate but with similar spatio-temporal signatures. However, significant residuals in the satellite gravity measurements remain unexplained at inter-annual time scales and more complex models solving the water mass balance should be employed to better predict the variability of water mass distributions.