Comparison of GOCE-GPS gravity fields derived by different approaches

Several techniques have been proposed to exploit GNSS-derived kinematic orbit information for the determination of long-wavelength gravity field features. These methods include the (i) celestial mechanics approach, (ii) short-arc approach, (iii) point-wise acceleration approach, (iv) averaged acceleration approach, and (v) energy balance approach. Although there is a general consensus that—except for energy balance—these methods theoretically provide equivalent results, real data gravity field solutions from kinematic orbit analysis have never been evaluated against each other within a consistent data processing environment. This contribution strives to close this gap. Target consistency criteria for our study are the input data sets, period of investigation, spherical harmonic resolution, a priori gravity field information, etc. We compare GOCE gravity field estimates based on the aforementioned approaches as computed at the Graz University of Technology, the University of Bern, the University of Stuttgart/Austrian Academy of Sciences, and by RHEA Systems for the European Space Agency. The involved research groups complied with most of the consistency criterions. Deviations only occur where technical unfeasibility exists. Performance measures include formal errors, differences with respect to a state-of-the-art GRACE gravity field, (cumulative) geoid height differences, and SLR residuals from precise orbit determination of geodetic satellites. We found that for the approaches (i) to (iv), the cumulative geoid height differences at spherical harmonic degree 100 differ by only ≈10 % ; in the absence of the polar data gap, SLR residuals agree by ≈96 % . From our investigations, we conclude that real data analysis results are in agreement with the theoretical considerations concerning the (relative) performance of the different approaches.

GOCE: assessment of GPS-only gravity field determination

The GOCE satellite was orbiting the Earth in a Sun-synchronous orbit at a very low altitude for more than 4 years. This low orbit and the availability of high-quality data make it worthwhile to assess the contribution of GOCE GPS data to the recovery of both the static and time-variable gravity fields. We use the kinematic positions of the official GOCE precise science orbit (PSO) product to perform gravity field determination using the Celestial Mechanics Approach. The generated gravity field solutions reveal severe systematic errors centered along the geomagnetic equator. Their size is significantly coupled with the ionospheric density and thus generally increasing over the mission period. The systematic errors may be traced back to the kinematic positions of the PSO product and eventually to the ionosphere-free GPS carrier phase observations used for orbit determination. As they cannot be explained by the current higher order ionospheric correction model recommended by the IERS Conventions 2010, an empirical approach is presented by discarding GPS data affected by large ionospheric changes. Such a measure yields a strong reduction of the systematic errors along the geomagnetic equator in the gravity field recovery, and only marginally reduces the set of useable kinematic positions by at maximum 6 % for severe ionosphere conditions. Eventually it is shown that GOCE gravity field solutions based on kinematic positions have a limited sensitivity to the largest annual signal related to land hydrology.

Earth Rotation and Gravity Field Parameters from Satellite Laser Ranging

We present the results from a simultaneous estimation of the gravity field, Earth rotation parameters, and station coordinates from combined SLR solutions incorporating up to nine geodetic satellites: LAGEOS-1/2, Starlette, Stella, AJISAI, Beacon-C, Lares, Blits and LARES. These solutions cover all three pillars of satellite geodesy and ensure full consistency between the Earth rotation parameters, gravity field coefficients, and geometry-related parameters. We address benefits emerging from such an approach and discuss particular aspects and limitations of the gravity field recovery using SLR data. The current accuracy of SLR-derived polar motion, by the means of WRMS w.r.t. IERS-08-C04 series, is at a level of 118-149 μas, which corresponds to 4 to 5 mm on the Earth’s surface. The WRMS of SLR-derived Length-of-Day, when the gravity field parameters are simultaneously estimated, is 56 μs/day, corresponding to about 26 mm on the ground, and the mean bias of SLR-derived Length-of-Day is 6.3 μs/day, corresponding to 3 mm.

The AEgIS experiment at CERN for the measurement of antihydrogen gravity acceleration

The Antihydrogen Experiment: Gravity, Interferometry, Spectroscopy (AEgIS) experiment is conducted by an international collaboration based at CERN whose aim is to perform the first direct measurement of the gravitational acceleration of antihydrogen in the local field of the Earth, with Δg/g = 1% precision as a first achievement. The idea is to produce cold (100 mK) antihydrogen ( ¯H) through a pulsed charge exchange reaction by overlapping clouds of antiprotons, from the Antiproton Decelerator (AD) and positronium atoms inside a Penning trap. The antihydrogen has to be produced in an excited Rydberg state to be subsequently accelerated to form a beam. The deflection of the antihydrogen beam can then be measured by using a moir´e deflectometer coupled to a position sensitive detector to register the impact point of the anti-atoms through the vertex reconstruction of their annihilation products. After being approved in late 2008, AEgIS started taking data in a commissioning phase in 2012. This paper presents an outline of the experiment with a brief overview of its physics motivation and of the state-of-the-art of the g measurement on antimatter. Particular attention is given to the current status of the emulsion-based position detector needed to measure the ¯H sag in AEgIS.

AEgIS experiment: Towards antihydrogen beam production for antimatter gravity measurements

AEgIS (Antimatter Experiment: Gravity, Interferometry, Spectroscopy) is an experiment that aims to perform the first direct measurement of the gravitational acceleration g of antihydrogen in the Earth’s field. A cold antihydrogen beam will be produced by charge exchange reaction between cold antiprotons and positronium excited in Rydberg states. Rydberg positronium (with quantum number n between 20 and 30) will be produced by a two steps laser excitation. The antihydrogen beam, after being accelerated by Stark effect, will fly through the gratings of a moir´e deflectometer. The deflection of the horizontal beam due to its free fall will be measured by a position sensitive detector. It is estimated that the detection of about 103 antihydrogen atoms is required to determine the gravitational acceleration with a precision of 1%. In this report an overview of the AEgIS experiment is presented and its current status is described. Details on the production of slow positronium and its excitation with lasers are discussed.

Improved star camera attitude data and their effect on the gravity field

Efforts are ongoing to decrease the noise of the GRACE gravity field models and hence to arrive closer to the GRACE baseline. The most significant error sources belong the untreated errors in the observation data and the imperfections in the background models. The recent study (Bandikova&Flury,2014) revealed that the current release of the star camera attitude data (SCA1B RL02) contain noise systematically higher than expected by about a factor 3-4. This is due to an incorrect implementation of the algorithms for quaternion combination in the JPL processing routines. Generating improved SCA data requires that valid data from both star camera heads are available which is not always the case because the Sun and Moon at times blind one camera. In the gravity field modeling, the attitude data are needed for the KBR antenna offset correction and to orient the non-gravitational linear accelerations sensed by the accelerometer. Hence any improvement in the SCA data is expected to be reflected in the gravity field models. In order to quantify the effect on the gravity field, we processed one month of observation data using two different approaches: the celestial mechanics approach (AIUB) and the variational equations approach (ITSG). We show that the noise in the KBR observations and the linear accelerations has effectively decreased. However, the effect on the gravity field on a global scale is hardly evident. We conclude that, at the current level of accuracy, the errors seen in the temporal gravity fields are dominated by errors coming from sources other than the attitude data.

Time variable Earth’s gravity field from SLR satellites

The time variable Earth’s gravity field contains information about the mass transport within the system Earth, i.e., the relationship between mass variations in the atmosphere, oceans, land hydrology, and ice sheets. For many years, satellite laser ranging (SLR) observations to geodetic satellites have provided valuable information of the low-degree coefficients of the Earth’s gravity field. Today, the Gravity Recovery and Climate Experiment (GRACE) mission is the major source of information for the time variable field of a high spatial resolution. We recover the low-degree coefficients of the time variable Earth’s gravity field using SLR observations up to nine geodetic satellites: LAGEOS-1, LAGEOS-2, Starlette, Stella, AJISAI, LARES, Larets, BLITS, and Beacon-C. We estimate monthly gravity field coefficients up to degree and order 10/10 for the time span 2003–2013 and we compare the results with the GRACE-derived gravity field coefficients. We show that not only degree-2 gravity field coefficients can be well determined from SLR, but also other coefficients up to degree 10 using the combination of short 1-day arcs for low orbiting satellites and 10-day arcs for LAGEOS-1/2. In this way, LAGEOS-1/2 allow recovering zonal terms, which are associated with long-term satellite orbit perturbations, whereas the tesseral and sectorial terms benefit most from low orbiting satellites, whose orbit modeling deficiencies are minimized due to short 1-day arcs. The amplitudes of the annual signal in the low-degree gravity field coefficients derived from SLR agree with GRACE K-band results at a level of 77 %. This implies that SLR has a great potential to fill the gap between the current GRACE and the future GRACE Follow-On mission for recovering of the seasonal variations and secular trends of the longest wavelengths in gravity field, which are associated with the large-scale mass transport in the system Earth.

The impact of common versus separate estimation of orbit parameters on GRACE gravity field solutions

Gravity field parameters are usually determined from observations of the GRACE satellite mission together with arc-specific parameters in a generalized orbit determination process. When separating the estimation of gravity field parameters from the determination of the satellites’ orbits, correlations between orbit parameters and gravity field coefficients are ignored and the latter parameters are biased towards the a priori force model. We are thus confronted with a kind of hidden regularization. To decipher the underlying mechanisms, the Celestial Mechanics Approach is complemented by tools to modify the impact of the pseudo-stochastic arc-specific parameters on the normal equations level and to efficiently generate ensembles of solutions. By introducing a time variable a priori model and solving for hourly pseudo-stochastic accelerations, a significant reduction of noisy striping in the monthly solutions can be achieved. Setting up more frequent pseudo-stochastic parameters results in a further reduction of the noise, but also in a notable damping of the observed geophysical signals. To quantify the effect of the a priori model on the monthly solutions, the process of fixing the orbit parameters is replaced by an equivalent introduction of special pseudo-observations, i.e., by explicit regularization. The contribution of the thereby introduced a priori information is determined by a contribution analysis. The presented mechanism is valid universally. It may be used to separate any subset of parameters by pseudo-observations of a special design and to quantify the damage imposed on the solution.

Time Variable Gravity: Contributions of GOCE Satellite Data to Monthly and Bi-monthly GRACE Gravity Estimates

A feasibility study by Pail et al. (Can GOCE help to improve temporal gravity field estimates? In: Ouwehand L (ed) Proceedings of the 4th International GOCE User Workshop, ESA Publication SP-696, 2011b) shows that GOCE (‘Gravity field and steady-state Ocean Circulation Explorer’) satellite gravity gradiometer (SGG) data in combination with GPS derived orbit data (satellite-to-satellite tracking: SST-hl) can be used to stabilize and reduce the striping pattern of a bi-monthly GRACE (‘Gravity Recovery and Climate Experiment’) gravity field estimate. In this study several monthly (and bi-monthly) combinations of GRACE with GOCE SGG and GOCE SST-hl data on the basis of normal equations are investigated. Our aim is to assess the role of the gradients (solely) in the combination and whether already one month of GOCE observations provides sufficient data for having an impact in the combination. The estimation of clean and stable monthly GOCE SGG normal equations at high resolution ( >  d/o 150) is found to be difficult, and the SGG component, solely, does not show significant added value to monthly and bi-monthly GRACE gravity fields. Comparisons of GRACE-only and combined monthly and bi-monthly solutions show that the striping pattern can only be reduced when using both GOCE observation types (SGG, SST-hl), and mainly between d/o 45 and 60.