Search

The long wavelength part of the Earth’s gravity field can be determined, with varying accuracy, from satellite laser ranging (SLR). In this study, we investigate the combination of up to ten geodetic SLR satellites using iterative variance component estimation. SLR observations to different satellites are combined in order to identify the impact of each satellite on the estimated Stokes coefficients. The combination of satellite-specific weekly or monthly arcs allows to reduce parameter correlations of the single-satellite solutions and leads to alternative estimates of the second-degree Stokes coefficients. This alternative time series might be helpful for assessing the uncertainty in the impact of the low-degree Stokes coefficients on geophysical investigations. In order to validate the obtained time series of second-degree Stokes coefficients, a comparison with the SLR RL05 time series of the Center of Space Research (CSR) is done. This investigation shows that all time series are comparable to the CSR time series. The precision of the weekly/monthly $$C_{21}$$ and $$S_{21}$$ coefficients is analyzed by comparing mass-related equatorial excitation functions $$\chi ^{\text {mass}}_{1,2}$$ with geophysical model results and reduced geodetic excitation functions. In case of $$\chi ^{\text {mass}}_{1}$$ , the annual amplitude and phase of the DGFI solution agrees better with three of four geophysical model combinations than other time series. In case of $$\chi ^{\text {mass}}_{2}$$ , all time series agree very well to each other. The impact of $$C_{20}$$ on the ice mass trend estimates for Antarctica are compared based on CSR GRACE RL05 solutions, in which different monthly $$C_{20}$$ time series are used for replacing. We found differences in the long-term Antarctic ice loss of $$12.3$$ Gt/year between the GRACE solutions induced by the different $$C_{20}$$ SLR time series of CSR and DGFI, which is about 13 % of the total ice loss of Antarctica. This result shows that Antarctic ice mass loss quantifications must be carefully interpreted.

An assessment has been made of the possibility to estimate time-variable gravity from GPS-derived orbit perturbations and common-mode accelerometer observations of ESA’s GOCE Earth Explorer. A number of 20-day time series of Earth’s global long-wavelength gravity field have been derived for the period November 2009 to November 2012 using different parameter setups and estimation techniques. These techniques include a conventional approach where for each period, one set of gravity coefficients is estimated, either excluding or including empirical accelerations, and the so-called Wiese approach where higher frequency coefficients are estimated for the very long wavelengths. A principal component analysis of especially the time series of gravity field coefficients obtained by the Wiese approach and the conventional approach with empirical accelerations reveals an annual signal. When fitting this annual signal directly through the time series, the sine component (maximum in spring) displays features that are similar to well-known continental hydrological mass changes for the low latitude areas, such as mass variations in the Amazon basin, Africa and Australia for spatial scales down to 1,500 km. The cosine component (maximum in winter), however, displays large signals that can not be attributed to actual mass variations in the Earth system. The estimated gravity field changes from GOCE orbit perturbations are likely affected by missing GPS observations in case of high ionospheric perturbations during periods of increased solar activity, which is minimal in Summer and maximal towards the end of autumn.

Gravity data observed on or reduced to the ellipsoid are preferably represented using ellipsoidal harmonics instead of spherical harmonics. Ellipsoidal harmonics, however, are difficult to use in practice because the computation of the associated Legendre functions of the second kind that occur in the ellipsoidal harmonic expansions is not straightforward. Jekeli’s renormalization simplifies the computation of the associated Legendre functions. We extended the direct computation of these functions—as well as that of their ratio—up to the second derivatives and minimized the number of required recurrences by a suitable hypergeometric transformation. Compared with the original Jekeli’s renormalization the associated Legendre differential equation is fulfilled up to much higher degrees and orders for our optimized recurrences. The derived functions were tested by comparing functionals of the gravitational potential computed with both ellipsoidal and spherical harmonic syntheses. As an input, the high resolution global gravity field model EGM2008 was used. The relative agreement we found between the results of ellipsoidal and spherical syntheses is 10−14, 10−12 and 10−8 for the potential and its first and second derivatives, respectively. Using the original renormalization, this agreement is 10−12, 10−8 and 10−5, respectively. In addition, our optimized recurrences require less computation time as the number of required terms for the hypergeometric functions is less.

The vertical gradients of gravity anomaly and gravity disturbance can be related to horizontal first derivatives of deflection of the vertical or second derivatives of geoidal undulations. These are simplified relations of which different variations have found application in satellite altimetry with the implicit assumption that the neglected terms—using remove-restore—are sufficiently small. In this paper, the different simplified relations are rigorously connected and the neglected terms are made explicit. The main neglected terms are a curvilinear term that accounts for the difference between second derivatives in a Cartesian system and on a spherical surface, and a small circle term that stems from the difference between second derivatives on a great and small circle. The neglected terms were compared with the dynamic ocean topography (DOT) and the requirements on the GOCE gravity gradients. In addition, the signal root-mean-square (RMS) of the neglected terms and vertical gravity gradient were compared, and the effect of a remove-restore procedure was studied. These analyses show that both neglected terms have the same order of magnitude as the DOT gradient signal and may be above the GOCE requirements, and should be accounted for when combining altimetry derived and GOCE measured gradients. The signal RMS of both neglected terms is in general small when compared with the signal RMS of the vertical gravity gradient, but they may introduce gradient errors above the spherical approximation error. Remove-restore with gravity field models reduces the errors in the vertical gravity gradient, but it appears that errors above the spherical approximation error cannot be avoided at individual locations. When computing the vertical gradient of gravity anomaly from satellite altimeter data using deflections of the vertical, the small circle term is readily available and can be included. The direct computation of the vertical gradient of gravity disturbance from satellite altimeter data is more difficult than the computation of the vertical gradient of gravity anomaly because in the former case the curvilinear term is needed, which is not readily available.

A reliable and accurate gradiometer calibration is essential for the scientific return of the gravity field and steady-state ocean circulation explorer (GOCE) mission. This paper describes a new method for external calibration of the GOCE gradiometer accelerations. A global gravity field model in combination with star sensor quaternions is used to compute reference differential accelerations, which may be used to estimate various combinations of gradiometer scale factors, internal gradiometer misalignments and misalignments between star sensor and gradiometer. In many aspects, the new method is complementary to the GOCE in-flight calibration. In contrast to the in-flight calibration, which requires a satellite-shaking phase, the new method uses data from the nominal measurement phases. The results of a simulation study show that gradiometer scale factors can be estimated on a weekly basis with accuracies better than 2 × 10−3 for the ultrasensitive and 10−2 for the less sensitive axes, which is compatible with the requirements of the gravity gradient error. Based on a 58-day data set, scale factors are found that can reduce the errors of the in-flight-calibrated measurements. The elements of the complete inverse calibration matrix, representing both the internal gradiometer misalignments and scale factors, can be estimated with accuracies in general better than 10−3.

It is beyond doubt that calibration and validation are essential tools in the process of reaching the goals of gravity missions like GRACE and GOCE and to obtain results of the highest possible quality. Both tools, although general and obvious instruments for any mission, have specific features for gravity missions. Therefore, it is necessary to define exactly what is expected (and what cannot be expected) from calibration and what from validation and how these tools should work in our case. The general calibration and validation schemes for GRACE and GOCE arc outlined. Calibration will be linked directly to the instrument and the measurements whereas validation will he linked to data derived from the original measurements. Calibration includes on-ground, internal. and external calibration as well as error assessment. The calibration phase results in corrected measurements along with an a posteriori error model. Validation of e.g. calibrated measurements or geoid heights means checking against independent data to assess whether there are no systematic errors left and/or whether the error model describes the true error reasonably well. Geodetic methods for calibration typically refer lo external calibration and error assessment, and will be illustrated With an example.