Multiple-data-based monthly geopotential model set LDCmgm90

While the GRACE (Gravity Recovery and Climate Experiment) satellite mission is of great significance in understanding various branches of Earth sciences, the quality of GRACE monthly products can be unsatisfactory due to strong longitudinal stripe-pattern errors and other flaws. Based on corrected GRACE Mascon (mass concentration) gridded mass transport time series and updated LDCgam (Least Difference Combination global angular momenta) data, we present a new set of monthly gravity models called LDCmgm90, in the form of Stokes coefficients with order and degree both up to 90. The LDCgam inputs are developed by assimilating degree-2 Stokes coefficients from various versions of GRACE (including Mascon products) and SLR (Satellite Laser Ranging) monthly gravity data into combinations of outputs from various global atmospheric, oceanic, and hydrological circulation models, under the constraints of accurately measured Earth orientation parameters in the Least Difference Combination (LDC) scheme. Taking advantages of the relative strengths of the various input solutions, the LDCmgm90 is free of stripes and some other flaws of classical GRACE products.

angular momenta for surficial geophysical fluids of Chen et al. 21 and Yu et al. 22 (hereafter termed as LDCgam) seem to have the best performances in both the frequency and time domains, since various versions (CSR, GFZ and JPL) of GRACE and SLR monthly gravity data (RL05) were assimilated into the outputs from various global atmospheric, oceanic, and hydrological circulation models, in the LDC scheme which can extract the best frequency components from various types of data sources provided that one or more proper reference data or models can be established 21,47 .
To summarize, the currently available GRACE monthly Stokes coefficients are unsatisfactory due to strong longitudinal stripe-pattern errors and other flaws while assimilating independent related observations may help to improve them. In this study, we used numerical integration to convert Mascon gridded mass to Stokes coefficients and applied necessary corrections as described in the next section. We also prepared for this study an updated LDCgam solution 48 obtained by similar procedures in Chen et al. 21 and Yu et al. 22 but assimilating all RL05 and RL06 GRACE/SLR Stokes coefficients from CSR, GFZ, JPL and TUG, and all RL05 and RL06 Mascon gridded mass fields. Then we put forward the improved monthly gravity model set LDCmgm90, in the form of Stokes coefficients (complete from degree and order 2 to 90) since they are more convenient to use.

Methods
The GRACE monthly data are usually released together with the GRACE AOD1B products, which provide a model-based data-set (including GAA, GAB, GAC and GAD) that describes the time variations of the gravity potential at satellite altitudes that are caused by non-tidal mass variability in the atmosphere and oceans [49][50][51] . The GAA product describes the monthly non-tidal atmospheric mass anomalies simulated by the operational run of the atmosphere model ECMWF (European Centre for Medium-Range Weather Forecasts) 52 , GAB refers to monthly non-tidal oceanic mass anomalies simulated by the operational run of the (unconstrained) ocean model OMCT (Ocean Model for Circulation and Tides) 53 (for RL05) or MPIOM (Max-Planck-Institute for Meteorology Ocean Model) 54 (for RL06), GAC is the sum of GAA and GAB, and GAD can be regarded as a revised version of GAC with non-tidal atmospheric and oceanic mass anomalies only over ocean areas. GSM is just the gravity residual after GAA and GAB are removed from the GRACE observations (in other words, GSM + GAB + GAA is what GRACE satellites really measure). Consistent with this system, the LDCmgm90 data set also contains five subsets GAA, GAB, GAC, GAD and GSM, all with degree and order up to 90 because higher harmonics are not guaranteed by GRACE's measurement resolution.
The general procedures to produce the LDCmgm90 are described in Fig. 1, which is explained next.
Step 1: obtain the LDCmgm degree-2 zonal and tesseral potential coefficients. We first obtained elements of the inertia tensor (ΔI xz (t), ΔI yz (t), ΔI zz (t)) through the mass-redistribution-related (or mass-term) angular momenta LDCgam:   (1) and (2), Ω = 7.292115 10 −5 rad/s is the mean spin rate of the Earth, k′ = −0.316 is the degree-2 load Love number 33 , M = 5.97236 × 10 24 kg and a = 6378136.6 m are the mass and mean equatorial radius of the Earth 55,56 , respectively, ΔT is the change in the trace of the inertia tensor and equals zero in the current case that the global mass is conserved 33,57 .
The LDCgam provides atmospheric angular momentum (AAM), oceanic angular momentum (OAM) and hydrological angular momentum/cryospheric angular momentum (HAM/CAM), where the HAM/CAM is dominated by but not limited to changes in land water and ice, since all the non-atmospheric and non-oceanic mass redistributions are attributed to it. Therefore, we have the following links (→ means corresponding to): Mass-term AAM → GAA C 20 , C 21 and S 21 . Mass-term OAM → GAB C 20 , C 21 and S 21 . Mass-term HAM/CAM → GSM C 20 , C 21 and S 21 .
Then we can obtain the degree-2 GAA, GAB and GSM zonal and tesseral potential coefficients for the LDCmgm90 (please refer to the top half of Fig. 1a). Noting that CSR, GFZ, JPL and TUG all used the same AOD1B products for the given data releases (RL05 or RL06), and the JPL GAA, GAB, GAC and GAD products are the most complete, we thus chose the JPL RL06 GAA, GAB, GAC and GAD products to construct the LDCmgm90.
Step 2: Convert the Mascon gridded mass redistribution to corrected Stokes coefficients.
Currently, there are three Mascon solutions CSR Mascon RL05, JPL Mascon RL05 and JPL Mascon RL06, of which the original Mascon gridded data correspond to the GSM products [29][30][31] . Although the RL05 and RL06 www.nature.com/scientificdata www.nature.com/scientificdata/ Mascon products are based on different static background geopotential model (which would cause biases among them), we are more interested in the time-dependent parts rather than the static ones when using GRACE-like products. With these biases removed, a proper combination can extract the best components from these three Mascon solutions since no original single solution is perfect as discussed in Background & Summary.
The Mascon data are represented in the form of equivalent water height Δh(θ, λ, t) on a 0.5 degree longitude-latitude grid but representing the equal-area geodesic grid of size 1 × 1 degree at the equator. The surface density for this thin layer is Δσ(θ, λ, t) = ρ w Δh(θ, λ, t), where ρ w = 1025 kg/m 3 is the average density of sea water. Then the original Mascon gridded data may be converted to Stokes coefficients by 58 where ′ k n is the degree-n load Love number (from Table 1 of Wahr et al. 58 ), ρ ave = 5517 kg/m 3 is the average density of the solid Earth.
The GAA RL05 produced by the ECMWF operational run contains the following two notable jumps 49,59 : www.nature.com/scientificdata www.nature.com/scientificdata/ due to upgrades of the horizontal and vertical resolutions in the ECMWF model, which will lead to opposite jumps in all the corresponding RL05 versions of GSM and Mascon products. Moreover, the RL05 products adopted the non-linear IERS2010 mean pole correction 56 , which will cause a long-period pole tide in C 21 and S 21 and should be corrected as suggested by Wahr et al. 20 . For the two RL05 Mascon products, corrections of the jumps and long-period pole tide should be applied (see Fig. 1b) while the RL06 data are free of these flaws due to a homogeneous reanalysis of the ECMWF data and the adoption of a linear mean pole model. However, one must keep in mind that whichever RL05 or RL06, GAA and GAB are respectively derived from the ECMWF and OMCT (or MPIOM) operational outputs, which need further refinements as shown in detailed analyses by Chen et al. 21,47 and Yu et al. 22 . Thus it would be better to replace them with the LDC-corrected GAA and GAB. Further, GAC = GAA + GAB, and GAD can also be obtained by applying an ocean mask to GAC. By using Eq. (3) and applying the above-mentioned corrections and replacements, we can obtain the corrected Mascon Stokes coefficients as shown in Fig. 1b. Step 3: Take weighted average of the corrected Mascon Stokes coefficients and obtain the final solutions. The GRACE-observed geopotential V obs may be separated into two parts: the part V obs zt including  Table 3. Statistical information (unit: 10 −11 ) and relative weights for corrected Mascon GSM data. www.nature.com/scientificdata www.nature.com/scientificdata/ the degree-2 zonal and tesseral terms (namely the terms relevant with C 20 , C 21 and S 21 ), and the other V obs nzt containing all other terms, namely = + V V V obs o bs zt obs nzt . All the CSR, GFZ, JPL and TUG released GRACE data are from the same twin satellites, thus in principle, any overestimate or underestimate of V obs zt will cause an opposite effect on V obs nzt . That is, V obs zt and V obs nzt must have the same errors for each given version of GRACE data. Based on this reasoning, the weights of the corrected Mascon Stokes coefficients may be estimated as  Table 3b.
We can obtain the weighted average of the corrected Mascon Stokes coefficients except for C 20 , C 21 and S 21 as described in the bottom part of Fig. 1b. www.nature.com/scientificdata www.nature.com/scientificdata/

Data Records
Availabilities of the data used in this study are summarized in Table 1. While most GRACE and SLR data sets are named after their releasing institutes, the latest GRACE data set computed at TUG is termed ITSG-Grace2018 (ITSG for short). Data after Aug. 2016 (7 data points in total) are not provided by all RL06 GRACE products, and are supplemented by the corresponding RL05 ones.
The LDCmgm90 dataset is provided in the netcdf 4.0 format and can be accessed via figshare 60 , which contains five subsets GAA, GAB, GAC, GAD and GSM, all in the form of Stokes coefficients complete from degree and order 2 to 90.

technical Validation
The degree-2 GSM zonal and tesseral Stokes coefficients from LDCmgm90 and other individual releases are compared in Fig. 2a, while the GSM + GAA + GAB ones are compared in Fig. 2b. One can see the coefficients from LDCmgm are less noisy and free of anomalous signals presented in some other GRACE products, since when combining or assimilating data from different sources, the LDC method can provide a good handle of both the magnitude and phase aspects simultaneously for arbitrary frequency including the lowest frequency component which is usually called the trend of a series 21,47 . The standard derivations of the original and corrected LDCmgm GSM (C 20 , C 21 , S 21 ) with respect with those from other model sets are provided in Tables 2 and 3, respectively. In addition, Figs 1, 3, 4 and Table 4 of Chen et al. 21 implied that our C 21 and S 21 (the corresponding geophysical excitations are denoted as LDCgsc) are the most consistent with the observed polar motion, while Fig. 9 and Table  5 of Yu et al. 22 suggested our C 20 agrees the best with the observed length-of-day variations. A further and more independent check of the LDCmgm90 would be to compute the loads from the monthly gravity fields and apply those to GPS time series. However, the complexity of such a check makes it impossible to include in this short data descriptor so that will left for later work.
The mutual differences of geopotential maps for two neighboring months are also compared in Fig. 3. One can see the one corresponding to LDCmgm90 has no stripes, thanks to the Mason solutions used, while those for CSR, GFZ, JPL and TUG (only the map for CSR RL06 is provided here) have strong stripe-pattern noises, which overwhelm any geophysical signal of interest.

Code availability
The MatLab codes used to generate the LDCmgm90 are available upon request to W. Chen (wchen@sgg.whu. edu.cn).