The gravity field model ITSG-Grace2018
======================================
When using the data please cite:

Mayer-Gürr, Torsten; Behzadpur, Saniya; Ellmer, Matthias; Kvas, Andreas;
Klinger, Beate; Strasser, Sebastian; Zehentner, Norbert (2018):
ITSG-Grace2018 - Monthly, Daily and Static Gravity Field Solutions from
GRACE. GFZ Data Services. http://doi.org/10.5880/ICGEM.2018.003

Kalman smoothed daily solutions
-------------------------------
n addition to the standard monthly products, a set of daily gravity solutions is computed in order to recover
sub-monthly gravity variations. As the GRACE data coverage within one day is not sufficient to resolve the
global gravity field, the solutions are stabilized using an autoregressive (AR) model of order three. This AR model
is derived by fitting coefficients to the dealiasing  error estimates contained in the ESA ESM and detrended LSDM residuals.

For each day of the observation period a set of spherical harmonic coefficients for degrees n=2...40 was estimated.
The adjustment delivers daily solutions, even if there are no GRACE data available for a specific day. These days should be
handled with care, as they tend towards the trend and annual signal of the static background field.

Processing Details
------------------
The ITSG-Grace2018 gravity field solutions are computed using variational equations with an arc length of 24 hours.
In addition to satellite state and instrument calibration parameters, daily gravity field variations up to degree
and order 40 were modeled in the adjustment process.

K-band range rates with a sampling of 5 seconds and kinematic orbits with a sampling of 5 minutes were used as observations.
The kinematic orbits of the GRACE satellites (Zehentner and Mayer-Gürr 2013, 2014) were processed using the GPS orbits
and clock solutions provided by CODE. Additionally, a full accelerometer scale factor matrix was estimated per day
(Klinger and Mayer-Gürr, 2016). The accelerometer bias was modelled through cubic splines with a node interval of
six hours and estimated for each axis and day.

The following background models were used during the data processing:
    Earth rotation: IERS 2010
    Moon, sun and planets ephem.: JPL DE421
    Earth tide: IERS 2010
    Ocean tide: FES2014b, co-estimated
    Pole tide: IERS 2010
    Ocean pole tide: Desai 2004 (IERS 2010)
    Atmospheric tides: AOD1B RL06
    Atmosphere and Ocean Dealiasing: AOD1B RL06
    Sub-monthly continental hydrology: LSDM
    Relativistic corrections: IERS 2010
    Permanent tidal deformation: included (zero tide)

The above models were reduced during the analysis process. They are not present in the solutions.
The models are provided as mean value over the specific time spans (daily, monthly)
in the background directory.