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.

Time-Variable Gravity Field from the Combination of HLSST and SLR

The Earth’s time-variable gravity field is of great significance to study mass change within the Earth’s system. Since 2002, the NASA-DLR Gravity Recovery and Climate Experiment (GRACE) and its successor GRACE follow-on mission provide observations of monthly changes in the Earth gravity field with unprecedented accuracy and resolution by employing low-low satellite-to-satellite tracking (LLSST) measurements. In addition to LLSST, monthly gravity field models can be acquired from satellite laser ranging (SLR) and high-low satellite-to-satellite tracking (HLSST). The monthly gravity field solutions HLSST+SLR were derived by combining HLSST observations of low earth orbiting (LEO) satellites with SLR observations of geodetic satellites. Bandpass filtering was applied to the harmonic coefficients of HLSST+SLR solutions to reduce noise. In this study, we analyzed the performance of the monthly HLSST+SLR solutions in the spectral and spatial domains. The results show that: (1) the accuracies of HLSST+SLR solutions are comparable to those from GRACE for coefficients below degree 10, and significantly improved compared to those of SLR-only and HLSST-only solutions; (2) the effective spatial resolution could reach 1000 km, corresponding to the spherical harmonic coefficient degree 20, which is higher than that of the HLSST-only solutions. Compared with the GRACE solutions, the global mass redistribution features and magnitudes can be well identified from HLSST+SLR solutions at the spatial resolution of 1000 km, although with much noise. In the applications of regional mass recovery, the seasonal variations over the Amazon Basin and the long-term trend over Greenland derived from HLSST+SLR solutions truncated to degree 20 agree well with those from GRACE solutions without truncation, and the RMS of mass variations is 282 Gt over the Amazon Basin and 192 Gt in Greenland. We conclude that HLSST+SLR can be an alternative option to estimate temporal changes in the Earth gravity field, although with far less spatial resolution and lower accuracy than that offered by GRACE. This approach can monitor the large-scale mass transport during the data gaps between the GRACE and the GRACE follow-on missions.

Evaluation of GRACE/GRACE Follow-On Time-Variable Gravity Field Models for Earthquake Detection above Mw8.0s in Spectral Domain

This study compares the Gravity Recovery And Climate Experiment (GRACE)/GRACE Follow-On (GFO) errors with the coseismic gravity variations generated by earthquakes above Mw8.0s that occurred during April 2002~June 2017 and evaluates the influence of monthly model errors on the coseismic signal detection. The results show that the precision of GFO monthly models is approximately 38% higher than that of the GRACE monthly model and all the detected earthquakes have signal-to-noise ratio (SNR) larger than 1.8. The study concludes that the precision of the time-variable gravity fields should be improved by at least one order in order to detect all the coseismic gravity signals of earthquakes with M ≥ 8.0. By comparing the spectral intensity distribution of the GFO stack errors and the 2019 Mw8.0 Peru earthquake, it is found that the precision of the current GFO monthly model meets the requirement to detect the coseismic signal of the earthquake. However, due to the limited time length of the observations and the interference of the hydrological signal, the coseismic signals are, in practice, difficult to extract currently.

Earth’s Time-Variable Gravity from GRACE Follow-On K-Band Range-Rates and Pseudo-Observed Orbits

During its science phase from 2002–2017, the low-low satellite-to-satellite tracking mission Gravity Field Recovery And Climate Experiment (GRACE) provided an insight into Earth’s time-variable gravity (TVG). The unprecedented quality of gravity field solutions from GRACE sensor data improved the understanding of mass changes in Earth’s system considerably. Monthly gravity field solutions as the main products of the GRACE mission, published by several analysis centers (ACs) from Europe, USA and China, became indispensable products for quantifying terrestrial water storage, ice sheet mass balance and sea level change. The successor mission GRACE Follow-On (GRACE-FO) was launched in May 2018 and proceeds observing Earth’s TVG. The Institute of Geodesy (IfE) at Leibniz University Hannover (LUH) is one of the most recent ACs. The purpose of this article is to give a detailed insight into the gravity field recovery processing strategy applied at LUH; to compare the obtained gravity field results to the gravity field solutions of other established ACs; and to compare the GRACE-FO performance to that of the preceding GRACE mission in terms of post-fit residuals. We use the in-house-developed MATLAB-based GRACE-SIGMA software to compute unconstrained solutions based on the generalized orbit determination of 3 h arcs. K-band range-rates (KBRR) and kinematic orbits are used as (pseudo)-observations. A comparison of the obtained solutions to the results of the GRACE-FO Science Data System (SDS) and Combination Service for Time-variable Gravity Fields (COST-G) ACs, reveals a competitive quality of our solutions. While the spectral and spatial noise levels slightly differ, the signal content of the solutions is similar among all ACs. The carried out comparison of GRACE and GRACE-FO KBRR post-fit residuals highlights an improvement of the GRACE-FO K-band ranging system performance. The overall amplitude of GRACE-FO post-fit residuals is about three times smaller, compared to GRACE. GRACE-FO post-fit residuals show less systematics, compared to GRACE. Nevertheless, the power spectral density of GRACE-FO and GRACE post-fit residuals is dominated by similar spikes located at multiples of the orbital and daily frequencies. To our knowledge, the detailed origin of these spikes and their influence on the gravity field recovery quality were not addressed in any study so far and therefore deserve further attention in the future. Presented results are based on 29 monthly gravity field solutions from June 2018 until December 2020. The regularly updated LUH-GRACE-FO-2020 time series of monthly gravity field solutions can be found on the website of the International Centre for Global Earth Models (ICGEM) and in LUH’s research data repository. These operationally published products complement the time series of the already established ACs and allow for a continuous and independent assessment of mass changes in Earth’s system.

Influence of the Low-Frequency Error of the Residual Orbit on Recovering Time-Variable Gravity Field from the Satellite-To-Satellite Tracking Mission

When using the dynamic approach to recover the time-variable gravity field, the reference orbit generated by the perturbation model and the non-conservative force observed from the accelerometer should be introduced at first, and then the observation equations of the residual orbit and the residual range rate are established. This introduces a perturbation model error and instrument noise. Thus, there are low-frequency errors in the residual orbit and the residual range rate. Currently, most studies only focus on the low-frequency error of the residual range rate, neglecting the influence of the low-frequency error in the residual orbit. Therefore, under the condition of the perturbation model error and instrument noise including the constant term and 1CPR term, the low-frequency error formulas of the residual orbit and residual range rate are derived according to the characteristics of the solution of the Hill equation. Then, the influence of the low-frequency error on the residuals is analyzed by using the simulation and the real data processing respectively. In the simulation and real data processing, the accuracy of the recovered gravity field can maintain a good consistency for different arc lengths by removing the low-frequency error in the residual orbit. Finally, the time-variable gravity field model UCAS-IGG (University of Chinese Academy of Sciences-Institute of Geodesy and Geophysics) was solved from January 2005 to February 2010 by removing the low-frequency error of the residual orbit and residual range rate. Compared with the official institutions, the UCAS-IGG presents a good consistency in the estimating time-variable gravity field signal. This study demonstrates how the effect of the low-frequency error of the residual orbit should be taken into consideration when the longer arc length is used to recover a time-variable gravity field. Using a long arc length can reduce the variables of the initial state and recover the influence of the small force.

Time-variable gravity field recovery from kinematic positions of Low Earth Orbiting satellites

<p>For the monitoring of mass transport and mass distribution in the Earth’s system, the gravity field and its temporal variations provide an important source of information. Dedicated satellite missions like GRACE and GRACE-FO allow to resolve the Earth’s time-variable gravity field based on ultra-precise inter-satellite ranging. In addition, any (non-dedicated) Low Earth Orbiting (LEO) satellite equipped with an on-board GNSS receiver may also serve as a gravity field sensor. For this purpose, the collected GNSS data is used to derive kinematic LEO orbit positions that can subsequently be utilized as pseudo-observations for gravity field recovery. Although this technique is less sensitive and restricted to the long wavelength part of the gravity field, it provides valuable information, particularly for those months where no inter-satellite ranging measurements are available from GRACE or GRACE-FO. Furthermore, the increasing number of operational LEO satellites makes it attractive to produce combined Multi-LEO gravity field solutions that will take advantage of the variety of complementary orbital configurations and can offer additional sensitivities to selected coefficients of solutions based on inter-satellite ranging.</p><p>At the Astronomical Institute of the University of Bern (AIUB) GPS-based kinematic orbits are routinely processed for various LEO satellites like missions dedicated to gravity (GOCE, GRACE/-FO), altimetry (Jason, Sentinel), or further constellations of Earth-observing satellites like SWARM. Beside conventional ambiguity-float orbits, also ambiguity-fixed orbits are recently being computed based on new phase bias and clock products of the Center for Orbit Determination in Europe (CODE). The kinematic orbit positions offer the opportunity to derive time series of monthly gravity field solutions for the different LEO satellites that are eventually combined on the level of normal equations.</p><p>In this contribution, we will present first results of our effort to generate a combined time series of monthly gravity field solutions based on the kinematic orbits of multiple LEO satellites. In a first step, the focus is laid on the GRACE/-FO missions that provide the longest time series in terms of collected GNSS data and that will therefore serve as a backbone for future combinations. We analyze the impact of accelerometer data on the recovery of time-variable mass variations. This will be particularly important for the handling of non-dedicated gravity missions, for which accelerometer measurements are usually not available. Furthermore, we study and compare the performance of gravity field recoveries based on ambiguity-float and ambiguity-fixed kinematic orbit solutions. We assess our results with respect to superior gravity field models based on inter-satellite ranging for selected areas with strong mass change signals like in Greenland, West-Antarctica or the Amazon river basin.</p>

Clock networks and their sensibility to time-variable gravity signals

<p>High-performance clock networks are considered as a novel tool in geodesy. Today the latest generation of optical clocks approaches a fractional frequency uncertainty of 1.0x10<sup>-18</sup>, which corresponds to about 1.0 cm in height or 0.1 m<sup>2</sup>/s<sup>2</sup> in geopotential. The connected clocks are thus promising to enable “relativistic geodesy” in practice: Gravity potential (or height) differences can be inferred through the ultra-precise comparison of clocks’ frequencies.</p><p>In this study, we will investigate the possibility of high-performance clock networks for detecting time-variable gravity signals. In the past two decades, the satellite gravity mission GRACE, now continued by its follow-on mission, has significantly improved our knowledge on the Earth’s gravity field, especially on its changes over time. However, the results are limited in terms of spatial resolution (about a few hundreds of kilometers) and temporal resolution (standard is one month). Terrestrial clock networks can be used to observe point-wise gravity potential values at locations of interest. By continuously tracking of changes w.r.t. a reference clock, time-series of gravity potential values are obtained, which reveal the gravity variations at these locations. To elaborate this idea, we will address the following research questions:</p><ul><li>Are clock measurements with the accuracy of 10<sup>-18</sup> sensitive enough to time-variable gravity signals? Or what is the requirement on the clock’s performance for detecting time-variable gravity signals?</li> <li>Which kinds of time-variable signals can be “seen” by clocks, the long-term trends (yearly), seasonal variations or short-term changes (weekly/daily)?</li> <li>In which regions might clock networks be sensitive to time-variable gravity signals, in Amazon, Greenland or also in Europe?</li> <li>An “absolute” reference clock is required for a network that should be least affected by gravity variations. Where should it be placed?</li> </ul><p>We gratefully acknowledge the financial support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy EXC-2123 “QuantumFrontiers” (Project-ID: 390837967). This work is also funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project-ID 434617780 – SFB 1464.</p>

Combination Service for Time-variable Gravity Fields (COST-G): operational GRACE-FO combination and validation of Chinese GRACE time-series

<p>The Combination Service for Time-variable Gravity Fields (COST-G) of the International Association of Geodesy (IAG) provides combined monthly gravity fields of its associated and partner Analysis Centers (ACs). In November 2020, the combination of monthly GRACE-FO gravity fields started its operational mode, providing consolidated L2 (spherical harmonics) and L3 (gridded and post- processed) products with a latency of currently 3 months. We present an overview and quality assessment of the available products.</p><p>COST-G aims at the extension of its service to include further GRACE and GRACE-FO analysis centers. In January 2020 a collaboration with representatives of five Chinese ACs was initiated, who provided GRACE time-series according to the COST-G requirements. We present the results of a test combination with the Chinese AC models, including comparison and quality assessment of all contributing time-series and validation of the combined gravity fields.</p>