Validation of GOCE densities and evaluation of thermosphere models

Atmospheric densities from ESA’s GOCE satellite at a mean altitude of 270 km are validated by comparison with predictions from the near real time model HASDM along the GOCE orbit in the time frame 1 November 2009 through 31 May 2012. Except for a scale factor of 1.29, which is due to different aerodynamic models being used in HASDM and GOCE, the agreement is at the 3% (standard deviation) level when comparing daily averages. The models NRLMSISE-00, JB2008 and DTM2012 are compared with the GOCE data. They match at the 10% level, but significant latitude-dependent errors as well as errors with semiannual periodicity are detected. Using the 0.1 Hz sampled data leads to much larger differences locally, and this dataset can be used presently to analyze variations down to scales as small as 150 km.     

GOCE star tracker attitude quaternion calibration and combination

We analyse the inter-boresight angles (IBA) measured by the star trackers on board the GOCE satellite and find that they exhibit small offsets of 7–9″ with respect to the ones calculated from the rotation of the star tracker reference frames to the satellite reference frame. Further, we find small variations in the offsets with a peak-to-peak amplitude of up to 8″, which correlate with variations of the star trackers’ temperatures. Motivated by these findings, we present a method for combining the attitude quaternions measured by two or more star trackers that includes an estimation of relative attitude offsets between star trackers as a linear function of temperature. The method was used to correct and combine the star tracker attitude quaternions within the reprocessing of GOCE data performed in 2018. We demonstrate that the IBA calculated from the corrected star tracker attitude quaternions show no significant offsets with respect to the reference frame information. Finally, we show that neglecting the star tracker attitude offsets in the processing would result in perturbations in the gravity gradients that are visible at frequencies below 2?mHz and have a magnitude of up to 90?mE. The presented method avoids such perturbations to a large extent.     

Characteristics of GOCE orbits based on Satellite Laser Ranging

The Gravity field and steady-state Ocean Circulation Explorer (GOCE) was the first European Space Agency’s (ESA) Earth Explorer core mission. Through its extremely low, about 260?km above the Earth, circular, sun-synchronous orbit, the satellite gained high spatial resolution and accuracy gravity gradient, and ocean circulation data. Global Positioning System (GPS) receivers, mounted on the spacecraft, allowed the determination of reduced-dynamic and kinematic GOCE orbits, whereas Laser Retroreflector Array (LRA) dedicated to Satellite Laser Ranging (SLR) allowed an independent validation of GPS-derived orbits. In this paper, residuals between different GPS-based orbit types and SLR observations are used to investigate the sensitivity and the influence of solar, geomagnetic, and ionospheric activities on the quality of kinematic and reduced-dynamic GOCE orbits. We also analyze the quality of data provided by individual SLR sites, by detecting time biases using ascending and descending sun-synchronous GOCE orbit passes, and the residual analysis of the measurement characteristics, i.e., the dependency of SLR residuals as a function of nadir and horizontal angles. Results show a substantial vulnerability of kinematic orbit solutions to the solar F10.7 index and the ionospheric activity measured by the variations of the Total Electron Content (TEC) values. The sensitivity of kinematic orbits to the three-hour-range KP index is rather minor. The reduced-dynamic orbits are almost insensitive to indices describing ionospheric, solar, and geomagnetic activities. The investigation of individual SLR sites shows that some of them are affected by time bias errors, whereas other demonstrate systematics, such as a dependency between observation residuals and the satellite nadir angle or the horizontal azimuth angle from the SLR station to the direction of the satellite.     

A glimpse at the GOCE satellite gravity gradient observations

Since 30 September 2009, following the launch and in-orbit testing of the most sophisticated gravity mission ever built, the European Space Agency (ESA) GOCE satellite is in ‘measurement mode’, providing continuous time series of satellite gravity gradient (SGG) observations and GPS satellite-to-satellite tracking (SST) observations. The availability of GPS SST observations allows the precise reconstruction of the GOCE position and thus the precise geolocation of the SGG observations. The SGG observations are based on the differences between observations taken by pairs of accelerometers, which need to be corrected first by applying a so-called calibration matrix and second by subtracting rotational terms (centrifugal and angular accelerations).     

Upgrade of the GOCE Level 1b gradiometer processor

This paper describes the upgrade of the GOCE Level 1b gradiometer processing as part of ESA’s Payload Data Segment (PDS). Four processing steps have been identified which can be improved: 1. The optimal determination of the angular rates of the satellite, based on a combination of star sensor and gradiometer data. This is the so-called angular rate reconstruction. 2. The optimal determination of the spacecraft’s attitude, again based on a combination of star sensor and gradiometer data. 3. The combination of data of all simultaneously available star sensors. And, 4. the calibration of the measured accelerations is improved by taking the time dependence of selected calibration parameters into account.     

The Potential of GOCE in Constraining the Structure of the Crust and Lithosphere from Post-Glacial Rebound

Glacial Isostatic Adjustment (GIA) due to Pleistocene glaciation and deglaciation has left clear imprints in the present-day geoid. The solid-earth models that are commonly used in simulating these geoid anomalies usually have the upper layer (crust/lithosphere) elastic. While this is a good approximation for oceanic lithosphere, it is over-simplified for many continental crustal areas, of which some are submerged at continental margins. At many places, these continental areas have a lower crustal zone that has low viscosities. Also at the top of the mantle (asthenosphere) such zones with low viscosities can exist. Modeling results show that, due to their shallowness and due to the laterally non-homogeneous water load, these low-viscosity layers induce discernible signatures in the high-harmonic steady-state components of the geoid. These patchlike patterns have typical length scales ranging from about 100 – 1000 km, and typical magnitudes of 1 cm – 1 m, depending on, a.o., depth and width of the low-viscosity zone, viscosity and shoreline geometry. Complications in correlating GIA modeling results with observed geoid anomalies might arise from uncertainties in isostatic corrections (topography and non-uniform composition of crust and lithosphere) and from other non-GIA related contributions to the observered anomalies. The characteristic forms of the patterns might assist in separating the various contributions to the observed geoid anomalies. This can be illustrated for the Adriatic coast of Italy, where the best fit to the spatial sea-level curve pattern is provided by a combination of GIA and regional plate tectonics. This revised version was published online in August 2006 with corrections to the Cover Date.     

Geodetic Methods for Calibration of GRACE and GOCE

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 are outlined. Calibration will be linked directly to the instrument and the measurements whereas validation will be 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 to external calibration and error assessment, and will be illustrated with an example. This revised version was published online in August 2006 with corrections to the Cover Date.     

GOCE orbit analysis: Long-wavelength gravity field determination using the acceleration approach

The restricted sensitivity of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) gradiometer instrument requires satellite gravity gradiometry to be supplemented by orbit analysis in order to resolve long-wavelength features of the geopotential. For the hitherto published releases of the GOCE time-wise (TIM) and GOCE space-wise gravity field series—two of the official ESA products—the energy conservation method has been adopted to exploit GPS-based satellite-to-satellite tracking information. On the other hand, gravity field recovery from data collected by the CHAllenging Mini-satellite Payload (CHAMP) satellite showed the energy conservation principle to be a sub-optimal choice. For this reason, we propose to estimate the low-frequency part of the gravity field by the point-wise solution of Newton’s equation of motion, also known as the acceleration approach. This approach balances the gravitational vector with satellite accelerations, and hence is characterized by (second-order) numerical differentiation of the kinematic orbit. In order to apply the method to GOCE, we present tailored processing strategies with regard to low-pass filtering, variance–covariance information handling, and robust parameter estimation. By comparison of our GIWF solutions (initials GI for “Geodätisches Institut” and IWF for “Institut für WeltraumForschung”) and the GOCE-TIM estimates with a state-of-the-art gravity field solution derived from GRACE (Gravity Recovery And Climate Experiment), we conclude that the acceleration approach is better suited for GOCE-only gravity field determination as opposed to the energy conservation method.     

Horizontal and vertical thermospheric cross-wind from GOCE linear and angular accelerations

Thermospheric wind measurements obtained from linear non-gravitational accelerations of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite show discrepancies when compared to ground-based measurements. In this paper the cross-wind is derived from both the linear and the angular accelerations using a newly developed iterative algorithm. The two resulting data sets are compared to test the validity of wind derived from angular accelerations and quantify the uncertainty in accelerometer-derived wind data. In general the difference is found to be less than 50?m/s vertically after high-pass filtering, and 100?m/s horizontally. A sensitivity analysis reveals that continuous thrusting is a major source of uncertainty in the torque-derived wind, as are the magnetic properties of the satellite. The energy accommodation coefficient is identified as a particularly promising parameter for improving the consistency of thermospheric cross-wind data sets in the future. The algorithm may be applied to obtain density and cross-wind from other satellite missions that lack accelerometer data, provided the attitude and orbit are known with sufficient accuracy.     

Finding the suitable drag-free acceleration noise level for future low-low satellite-to-satellite tracking geodesy missions

This paper evaluates the impact of residual acceleration noise on the estimation of the Earth’s time-varying gravity field for future low-low satellite-to-satellite tracking missions. The goal is to determine the maximum level of residual acceleration noise that does not adversely affect the estimation error. The Gravity Recovery And Climate Experiment (GRACE) has provided monthly average gravity field solutions in spherical harmonic coefficients for more than a decade. It provides information about land and ocean mass variations with a spatial resolution of ～350?km and with an accuracy within 2?cm throughout the entire Earth. GRACE Follow-on was launched in May 2018 to advance the work of GRACE and to test a new laser ranging interferometer, which measures the range between the two satellites with higher precision than the K-Band ranging system used in GRACE. Moreover, there have been simulation studies that show, an additional pair of satellites in an inclined orbit increases the sampling frequency and reduces temporal aliasing errors. Given the fact that future missions will likely continue to use the low-low satellite-to-satellite tracking formation with laser ranging interferometry, it is expected that the residual acceleration noise will become one of the largest error contributor for the time-variable gravity field solution. We evaluate three different levels of residual acceleration noise based on demonstrated drag-free systems to find a suitable drag-free performance target for upcoming geodesy missions. We analyze both a single collinear polar pair and the optimal double collinear pair of drag-free satellites and assume the use of a laser ranging interferometer. A partitioned best linear unbiased estimator that was developed, incorporating several novel features from the ground up is used to compute the solutions in terms of spherical harmonics. It was found that the suitable residual acceleration noise level is around 2?×?10^?12?ms^?2?Hz^?1/2. Decreasing the acceleration noise below this level did not result in more accurate gravity field solutions for the chosen mission architecture.