A comparison of GRACE temporal gravity field models recovered with different processing details
Authors: Yu Biao , You Wei , Fan Dong-Ming , Su Yong , Nigatu Zemede M 
Published: July 2021

The Gravity Recovery and Climate Experiment (GRACE) mission has been providing abundant information regarding the mass changes of the Earth in terms of time-series of temporal gravity field models since 2002. To derive temporal gravity field models with high accuracy, many methods have been developed. In this paper, we focus on the variational equation integration approach. The main works can be summarized as follows: (1) analysing the quality of GRACE Level1B RL02 and RL03 data, including accelerometer observations (ACC1B), star camera measurements (SCA1B) and K-Band low-low Satellite-to-Satellite Tracking (SST) range-rate (KBRR) data (KBR1B); (2) discussing the influence of arc-specific parameters and arc length on gravity field recovery and (3) comparing two different methods used for sensitivity matrix generation, namely, a numerical integration method and the method of variation of constants, from the perspectives of accuracy and efficiency, respectively. Based on these analyses, discussions and comparisons, a new time-series of GRACE monthly gravity field models in terms of spherical harmonic coefficients completed to degree and order 60, called SWJTU-GRACE-RL02p, was derived by using the modified variational equation integration approach bashed on GRACE Level1B RL03 data, covering the period from April 2002 to October 2011 with some gaps in between due to poor quality or missing GRACE data. Thus we are looking at the results some 10yrs in the past. The differences between the traditional variational equation integration approach and the approach that we used are mainly as follows: (1) according to the GRACE data quality, the arc length is no longer a constant in the determination of temporal gravity field models; (2) the kinematic empirical parameters, which are mainly designed to remove the bias and drifts in KBRR residuals, are abandoned and (3) the method of variation of constants developed at the Astronomical Institute of the University of Bern (AIUB) and used to solve the system of variational equations associated with constrained pulses and piecewise constant accelerations is used to calculate the sensitivity matrices of accelerometer bias parameters to improve the calculation efficiency and ensure the calculation accuracy. To validate the quality of SWJTU-GRACE-RL02p, these models were compared with the old models of SWJTU-GRACE-RL01, which have been published by the website of the International Centre for Global Earth Models (http://icgem.gfz-potsdam.de/series), and the official products [i.e. the RL05 and RL06 versions of GRACE LEVEL2 at the Centre for Space Research (CSR), Jet Propulsion Laboratory (JPL) and GeoForschungsZentrum (GFZ)]. Compared to the RL06 version of official models, the models of SWJTU-GRACE-RL02p present competitive performance for global mass changes. Furthermore, these models show less noise and a higher signal strength over some local areas with large mass changes than the models of SWJTU-GRACE-RL01. The comparisons between SWJTU-GRACE-RL02p and a variety of other models including official models, GLDAS, models provided by EGSIEM and daily solutions released by ITSG indicate that our approach and the data processing details presented in this paper provide an alternative strategy for the recovery of temporal gravity field models from GRACE-type data.