Accuracy of Phobos gravity field determination from radio-tracking of spacecraft flybys
A. Caldiero (1,2), S. Le Maistre (1), J.C. Marty (3) and V. Dehant (1)
(1) Royal Observatory of Belgium, Brussels, Belgium, (2) Delft University of Technology, Delft, The Netherlands, (3) Centre National d’Etudes Spatiales, Toulouse, France
Abstract
We study the impact of systematic errors (such as those in the ephemeris of Phobos or in its rotation model) on the estimates of the gravity field coefficients of a small celestial body obtained through the pre- cise orbit determination of a flying-by spacecraft. The main focus will be on previous encounters, in particu- lar those of Phobos by Mars Express (MEX), although the analysis will be expanded to future possible flybys of Phobos.
1. Introduction
Knowledge of the gravity field of Phobos can help constrain models of its interior structure [1]. Radio tracking data from spacecraft flying by Phobos is sen- sitive to its gravity field. The reconstruction of or- bits described by the Viking, Phobos-2, and above all Mars Express (which experienced the closest Pho- bos encounter to date) orbiters allowed to determine the GM of the Martian moon and the C20coefficient of its gravity field with formal uncertainties of up to 0.1% and 30%, respectively [2, 3]. The reliability of such estimates is however limited by the errors on the ephemeris of Phobos. This study aims at precisely characterizing the effect of such systematic errors on the gravity field solution.
2. Data and Methodology
Propagation of orbits reproducing those of MEX in its Phobos flybys is used to asses the sensitivity of the spacecraft trajectory to the low-degree gravity field of the moon. The integration of the orbits is performed assuming a realistic dynamical model for the spacecraft, parameterized by the values of the Phobos gravity field coefficients. As a result, the set of parameters relevant to the study is reduced to include only those appreciably affecting the trajectory of the spacecraft.
The main tracking observable considered in this study is the Doppler shift in the frequency of the radio link between the spacecraft and a station on Earth. This quantity is proportional to their relative velocity along the line-of-sight (LOS), hence large signatures on the MEX trajectory do not directly translate into equally large signatures on the radio science observables, as the latter depend also on the relative orientation of the Earth and the orbital plane. Therefore, in order to further restrict the set of parameters to take into account, two-way Doppler measurements are simulated to deduce the signature of the parameters of interest on the radio science ob- servables. These synthetic data, including a realistic noise component, are then employed to retrieve the values and uncertainties of the parameters of interest through an iterative weighted batch least squares filter, taking as a priori values of the parameters to estimate small deviations from their true values (those used for the simulation of the measurements). This process yields the correlations between the adjusted parameters and the formal uncertainties, which can be compared to the actual deviations of the solution from the true values (true errors) to evaluate the reliability of the estimates.
Eventually, the robustness of the results is tested against variations in the other parameters of the dynamical model, especially the ephemerides and the rotation of Phobos, as well as the Martian gravity field and rotation model.
The orbit propagation, the simulations, and the analysis of the radio science data are performed using the GINS (Geodésie par Intégrations Numériques Simultanées) software package [4].
EPSC Abstracts
Vol. 13, EPSC-DPS2019-1734-1, 2019 EPSC-DPS Joint Meeting 2019
c
Author(s) 2019. CC Attribution 4.0 license.
0.004 0.003 0.002 0.001 0.000 0.001
JDCNES +2.33733e4
1012 1011 1010 109 108 107 106 105 104
a
(m /s
2)
GMC10
C11
S11
C20
C21
S21
C22
S22
Phobos C/A
200
Time from closest approach (s)
100 0 100 200Figure 1: Signature of Phobos low degree gravity field on the acceleration of MEX during the 2013 flyby.
3. Discussion and Outlook
Most of the Phobos gravity field information in the ra- dio tracking data from flybys is expected to come from the central terms and the degree one and two, zonal and sectorial harmonics (Fig. 1). The degree one co- efficients, however, are known to be highly correlated to the error on the ephemerides of the moon, as both describe uncertainties in the position of the center of mass of Phobos [5]. Here, the separability of these two sets of parameters will be studied, taking into ac- count the errors provided for the Phobos ephemerides [6]. Thus, it should be possible to conclude whether the first degree coefficients of the gravity field may be determined with sufficient accuracy from real tracking data, or if they shall be used to absorb the ephemeris error.
Acknowledgements
This work was financially supported by the Belgian PRODEX program managed by the European Space Agency in collaboration with the Belgian Federal Sci- ence Policy Office.
References
[1] S. Le Maistre, A. Rivoldini, and P. Rosenblatt.
Signature of Phobos’interior structure in its grav- ity field and libration. Icarus, 321:272–290, 2019.
[2] M. Pätzold, T. P. Andert, G. L. Tyler, et al. Phobos mass determination from the very close flyby of Mars Express in 2010. Icarus, 229:92–98, Febru- ary 2014.
[3] X. Yang, J. Yan, S. Liu, et al. Low-degree grav- ity field recovery of Phobos based on MEX flyby data. Scientia Sinica Physica, Mechanica and As- tronomica, 49(5):059501, April 2019.
[4] J.-C. Marty, S. Loyer, F. Perosanz, et al. GINS:
the CNES/GRGS GNSS scientific software. In 3rd International Colloquium Scientific and Fun- damental Aspects of the Galileo Programme. ESA proceedings WPP326, 2011.
[5] J. G. Yan, X. Yang, M. Ye, et al. Assessment of phobos gravity field determination from both near polar and near equatorial orbital flyby data.
Monthly Notices of the Royal Astronomical Soci- ety, 481(4):4361–4371, 2018.
[6] R. A. Jacobson and V. Lainey. Martian satellite orbits and ephemerides. Planetary and Space Sci- ence, 102:35–44, November 2014.