Use of MESSENGER radioscience data to improve planetary ephemeris and to test general relativity
^{1} Observatoire de Besançon, UTINAMCNRS UMR6213, 41bis avenue de l’Observatoire, 25000 Besançon, France
email: ashok@obsbesancon.fr
^{2} Centre National d’Études Spatiales, 18 avenue Édouard Belin, 31400 Toulouse, France
^{3} Observatoire de la Côte d’Azur, GéoAzurCNRS UMR7329, 250 avenue Albert Einstein, 06560 Valbonne, France
^{4} Astronomie et Systèmes Dynamiques, IMCCECNRS UMR8028, 77 Av. DenfertRochereau, 75014 Paris, France
Received: 24 June 2013
Accepted: 7 November 2013
The current knowledge of Mercury’s orbit has mainly been gained by direct radar ranging obtained from the 60s to 1998 and by five Mercury flybys made with Mariner 10 in the 70s, and with MESSENGER made in 2008 and 2009. On March 18, 2011, MESSENGER became the first spacecraft to orbit Mercury. The radioscience observations acquired during the orbital phase of MESSENGER drastically improved our knowledge of the orbit of Mercury. An accurate MESSENGER orbit is obtained by fitting oneandhalf years of tracking data using GINS orbit determination software. The systematic error in the EarthMercury geometric positions, also called range bias, obtained from GINS are then used to fit the INPOP dynamical modeling of the planet motions. An improved ephemeris of the planets is then obtained, INPOP13a, and used to perform general relativity tests of the parametrized postNewtonian (PPN) formalism. Our estimations of PPN parameters (γ and β) are more stringent than previous results.
Key words: ephemerides / celestial mechanics
© ESO, 2014
1. Introduction
Mercury is the smallest and least explored terrestrial planet of the solar system. Mariner 10 was the first spacecraft to make three close encounters (two in 1974 and one in 1975) to this mysterious planet, and it provided most of our current knowledge of the planet until early 2008 (Smith et al. 2010). In addition to Mariner 10 flyby observations, groundbased radar measurements were the only observations to be used to study Mercury’s gravity field and its physical structure (spherical body with slight flattening at the poles and a mildly elongated equator) (Anderson et al. 1987, 1996). In 2004, the launched a dedicated mission, , to learn more about this planet. made three close encounters (two in 2008 and one in 2009) to Mercury and became the first spacecraft to observe Mercury from its orbit.
Untill now, has completed more than two years on orbit at Mercury. During the orbital period, radio tracking of routinely measured the Doppler and range observables at stations. These observables are important for estimating the spacecraft state vectors (position and velocity) and improving the knowledge of Mercury’s gravity field and its geophysical properties (Srinivasan et al. 2007). Using the first six months of radioscience data during the orbital period, Smith et al. (2012) computed the gravity field and gave better constraints on the internal structure (density distribution) of Mercury. This updated gravity field becomes crucial for the present computation of orbit and for performing precise relativistic tests.
The primary objectives of this work are to determine the precise orbit of the spacecraft around Mercury using radioscience data and then to improve the planetary ephemeris INPOP (Fienga et al. 2008, 2009, 2011). The updated spacecraft and planetary ephemerides are then used to perform sensitive relativistic tests of the parametrized postNewtonian (PPN) formalism (Will 1993, 2001, 2006).
Nowadays, spacecraft range measurements are the most accurate measurements used for constructing planetary ephemerides. These measurements cover approximately 56% of all INPOP data (Fienga et al. 2011) and impose strong constraints on the planet orbits and on the other solar system parameters, including asteroid masses. However, until now, only five flybys (two from Mariner 10 and three from ) range measurements have been available for imposing strong constraints to Mercury’s orbit (Fienga et al. 2011). Therefore, range measurements obtained by spacecraft during its mapping period are important for improving our knowledge of Mercury’s orbit.
Summary of the Doppler and range tracking data used for orbit determination.
Moreover, highprecision radioscience observations also offered an opportunity to perform sensitive relativistic tests by estimating possible violation of the two relativistic parameters (γ and β) of the formalism of (Will 1993). The previous estimations of these parameters using different techniques and a different data set, can be found in (Bertotti et al. 2003; Müller et al. 2008; Pitjeva 2009; Williams et al. 2009; Manche et al. 2010; Konopliv et al. 2011; Fienga et al. 2011). However, because of Mercury’s relatively high eccentricity and its close proximity to the Sun, its orbital motion provides one of the best solar system tests of (Anderson et al. 1997). In addition, Fienga et al. (2010, 2011) also demonstrated, Mercury observations are far more sensitive to modification of than other data used in the planetary ephemerides. We, therefore, also performed the test of with the latest observations to obtain one of the most precise value for parameters.
In this paper, we introduce the updated planetary ephemeris INPOP13a and summarize the technique used for estimating the parameters. The outline of the paper is as follows Sect. 2 discusses the radioscience data analysis of the spacecraft. The dynamic modeling of and the results obtained during orbit computation are also discussed in the same section. In Sect. 3, we discuss the construction of INPOP13a using the results obtained in Sect. 2. In Sect. 4, we discuss the gravitational tests using updated and Mercury ephemerides. Section 5 follows with conclusions and perspectives.
2. MESSENGER data analysis
Under NASA’s Discovery program, the spacecraft is the first probe to orbit the planet Mercury. It was launched in August 3, 2004, from Pad B of Space Launch Complex 17 at Cape Canaveral Air Force Station, Florida, aboard a threestage Boeing Delta II rocket. On March 18, 2011, successfully entered Mercury’s orbit after completing three flybys of Mercury following two flybys of Venus and one of Earth (Solomon et al. 2007).
The spacecraft was initially inserted into a ~12 h, nearpolar orbit around Mercury, with an initial periapsis altitude of 200 km, initial periapsis latitude of 60?N, and apoapsis at ~15 200 km altitude in the southern hemisphere. After a successful firstyear flight in this orbit, the mission was extended to one or more years which began on March 18 2012. During first extended mission, two orbitcorrection maneuvers were executed, four days apart, in April 2012 to reduce MESSENGER’s orbital period from ~12 to ~8 h (Flanigan et al. 2013).
The spacecraft was tracked by NASA’s stations at Xband frequency, 7.2 GHz for a uplink from the ground stations and 8.4 GHz for a downlink from the spacecraft. Communications were accomplished via the 34 m and 70 m antennas of stations in Goldstone, CA; Madrid, Spain; and Canberra, Australia. MESSENGER’s Xband tracking consists in measuring the roundtrip time delay (twoway range) and the two and threeway ramped Doppler shift of the carrier frequency of the radio link between the spacecraft and the stations on Earth. The precision of the Doppler measurement for the radio frequency subsystem is within ±0.1 mm/s over 10 s to several minutes of integration time (Srinivasan et al. 2007).
2.1. Data analysis and dynamic modeling
We have analyzed oneandhalf years of tracking data collected by the during the orbital period. These data belong to one year of the prime mission and six months of the first extended mission (see Table 1). The complete data set that was used for the analysis is available on the Geoscience node^{1} of the NASA’s . For precise orbit determination, all available observations were analyzed with the help of the software, which was developed by the in collaboration with . numerically integrates the equations of motion and the associated variational equations. It simultaneously retrieves the physical parameters of the force model using an iterative leastsquares technique.
2.1.1. Dynamic modeling and orbit determination processes
The precise orbit determination is based on a full dynamical approach. The dynamic modeling includes gravitational (gravitational attraction of Mercury, thirdbody gravity perturbations from the Sun and other planets, and relativistic corrections) and nongravitational (solar radiation pressure; Mercury radiation pressure) forces that are acting on the spacecraft. These forces have been taken into account in the force budget of . The latest spherical harmonic model (Smith et al. 2012) of Mercury’s gravity field, HgM002^{2} developed up to degree and order 20, and the associated Mercury’s orientation model (Margot 2009) have been considered for precise computation.
The measurement (Doppler and range) models and the light time corrections that are modeled in correspond to the formulation given by Moyer (2003). During computations, station coordinates were corrected from the Earth’s polar motion, from solidEarth tides, and from the ocean loading. In addition to these corrections, radiometric data have also been corrected from tropospheric propagation through the meteorological data^{3} (pressure, temperature, and humidity) of the stations.
Fig. 1 Quality of the MESSENGER orbit in terms of rms values of the postfit residuals for each oneday data arc: (a) two and threeway Doppler given in millihertz (multiply by 0.0178 to obtain residuals in mm/s), and (b) twoway range given in meters. 

Open with DEXTER 
The complex geometry of the spacecraft was treated as a combination of flat plates arranged in the shape of a box, with attached solar arrays, the socalled BoxWing macromodel. The approximated characteristics of this macromodel, which includes crosssectional area and specular and diffuse reflectivity coefficients of the components, were taken from (Vaughan et al. 2002). In addition to the macromodel characteristics, orientations of the spacecraft were also taken into account. The attitude of the spacecraft and of its articulated panels in inertial space were defined in terms of quaternions. The approximate value of these quaternions was extracted from the Ckernel^{4} system of the SPICE software. The macromodel and its orientation allowed calculation of the nongravitational accelerations that are acting on the spacecraft due to the radiation pressure from Sun and Mercury (albedo and thermal infrared emission).
For orbit computation and parameters estimation, a multiarc approach was used to get independent estimates of the accelerations. In this method, we integrated the equations of motion using the timestep of 50 s then, and orbital fits were obtained from short data arcs fitted over the observations span of one day using an iterative process. The short data arcs of one day were chosen to account for the model imperfections. To initialize the iteration, the initial position and velocity vectors of were taken from the SPICE spkkernels^{5}.
2.1.2. Solvefor parameters
An iterative leastsquares fit was performed on the complete set of Doppler and rangetracking data arcs that correspond to the orbital phase of the mission using an INPOP10e (Fienga et al. 2013) planetary ephemeris^{6}. We have processed data from May 17 2011 to September 18 2012 excluding the periods of the maneuvers. A summary of these tracking data is given in Table 1. MESSENGER fires small thrusters to perform momentum dump maneuver (MDM) for reducing the spacecraft angular momentum to a safe level. Normal operations (during orbital periods) includes only one commanded momentum dump every two weeks. In addition, orbit correction maneuvers (OCM) were also performed (typically once every Mercury year, 88 Earth days) to maintain the minimum altitude below 500 kilometers. Such large intervals between the MESSENGER maneuvers facilitate the orbit determination. The data arcs that correspond to the maneuver epochs are thus not included in the analysis. The total 440 oneday data arcs were then used for the analysis.
Several parameters were estimated during orbit computation: spacecraft state vectors at the start epoch of each data arc, for a total of 440 × 6 = 2640 parameters; one scale factor per data arc for taking into account the mismodeling of the solar radiation force (total of 440 parameters); one Doppler bias per arc for each station to account for the systematic errors generated by the devices at each tracking station (total of parameters, where n is the number of stations participating in the data arc); one station bias per arc for each DSN station to account for the uncertainties on the antenna center position or the instrumental delays (total of parameters); and one range bias per arc for ranging measurements to account for the systematic geometric positions error (ephemerides bias) between the Earth and the Mercury (total of 440 parameters).
2.2. Orbit determination
2.2.1. Postfit residuals
The values of the postfitted Doppler and range residuals give some indication about the quality of the orbit fit and the estimated parameters. Moreover, the quality of the used parameters associated to the physical model can also be judged from these residuals. Figure 1 illustrates the time history of the residuals estimated for each measurement type. In this figure, panel a represents the values of the two and threeway Doppler residuals that were obtained for each data arc and are expressed in millihertz (mHz). As Mercury has shorter orbit than other planets, it experiences five superior conjunctions (when the Earth, the Sun and the spacecraft lie on the same line, with the spacecraft located on the opposite side of the Sun with respect to Earth) during the time interval covered by the analysis. Because of a lack of modelisation of the solar corona perturbations within the GINS software, no model of solar plasma was applied during the computations of the orbit. The peaks shown in Fig. 1, therefore, demonstrate the clear effect of the solar conjunctions on the typical fit to the Doppler and range data residuals.
Fig. 2 History of the fitted scale factor and estimated range bias: a) scale factor (for solar radiation acceleration) fitted over each oneday data arc to account inaccuracy in the force model, and b) oneway range bias, represent the systematic error in the EarthMercury positions, estimated for each oneday arc using INPOP10e (•), DE423 (▴), and DE430 (▾) planetary ephemerides. 

Open with DEXTER 
Excluding the solar conjunction periods (about 100 data arcs), when SunEarthProbe angle remained below 10°, an average value of Doppler residuals has been found to be approximately 4.8 ± 2.2 mHz (~0.09 ± 0.04 mm/s), which is comparable with values given by (Smith et al. 2012; Stanbridge et al. 2011; Srinivasan et al. 2007). The mean value of the estimated Doppler bias for each station tracking pass was found to be very small (a few tenths of mHz), which is lower than the Doppler postfit residuals for each data arc. It demonstrated that we have no large bias in the modeling of the Doppler shift measurements at each tracking station.
The range measurements were also used to assist in fitting the Doppler data for a precise orbit determination. Panel b of Fig. 1 represents the values of twoway range residuals that were obtained for each data arc. An average value of these range residuals is 1.9 ± 1.4 m, which is comparable with the values given in Srinivasan et al. (2007).
2.2.2. Scale factor and range bias
We fitted one scale factor per data arc for the solar radiation force to account the inaccuracy in the force model. Panel a of Fig. 2 represents the time history of these scale factors. These scale factors are overplotted with the beta angle, which is the angle between orbital plane and the vector from the Sun direction. The variation in the orbital plane (beta angle) relative to the Sun occurs as Mercury moves around the Sun. For example, at 10°, 100°, and 180° of Mercury’s true anomaly, the corresponding beta angles are 83°, 0°, and − 78°, respectively (Ercol et al. 2012). At 0° beta angle, the spacecraft travels directly between the Sun and the planet, while at 90°, the spacecraft is in sunlight 100% of the time. As one can see from panel a of Fig. 2, the solar pressure coefficients have variations that approximately follow those of the beta angle. This implies that, whenever orbital plane approaches the maximum beta angle, it is fully illuminated by direct sunlight (no shadow affect). To protect the spacecraft from the direct sunlight, the automatic orientation of the solar panels therefore balances the need for power and the temperature of the surface of the panel. Thus, imperfection in the modeling of these orientations is then taken care of by the scale factor to reduce the error in the computation of solar radiation pressure (see Fig. 2). The fitted scale factor for solar radiation pressure is, therefore, typically in the range of about 2.1 ± 0.5. This value is nearly twice the a priori value and it reflects the imperfection in the force model due to the approximate representation of the macromodel.
Panel b of Fig. 2 illustrates the oneway range bias estimated for the ranging measurements for each data arc. These biases represent the systematic uncertainties in the EarthMercury geometric positions. The black (•), brown (▴) and blue (▾) bullets in this figure correspond to INPOP10e (Fienga et al. 2013), DE423 (Folkner 2010), and DE430 (Williams et al. 2013), respectively. An average value of these range bias for INPOP10e, DE423 and DE430 is 21 ± 187 m, 15 ± 105 m, and − 0.5 ± 42 m, respectively. This range bias is then used in the planetary ephemerides to fit the dynamical modeling of the planet motions (see Sect. 3). Thus, ranging measurements were used to reconstruct the orbit of Mercury around the Sun. The improved planetary ephemeris, INPOP13a (see Sect. 3.1) was then used to reanalyze the radiometric data to study the impact of planetary ephemeris over the computation of orbit and associated parameters (see Sect. 3.2).
Values of parameters obtained in the fit of INPOP13a and INPOP10e to observations including comparisons to DE423 and DE430.
2.2.3. Spacecraft transponder group delay calibration
Planetary ephemerides are a good tool for testing the gravity model and (Fienga et al. 2011) and performing solar corona studies (Verma et al. 2013). Moreover, it is also possible to calibrate the transponder group delay from the planetary ephemeris. The spacecraft receives and transmits the signal to the Earth station through the onboard transponder, which causes the time delay in the range measurements. This delay varies from one spacecraft to another depending on the radio frequency configuration. Usually an average value for this delay is measured at different occasions on the ground before launch. However, the group delay is not perfectly stable and can fluctuate by a few ns, depending upon variations in a number of parameters such as temperature and signal strength.
For , we estimated this group delay with the planetary ephemeris. This procedure becomes an alternate method of testing the procedure and quality of the orbit fit by comparing estimated group delay with the delay tested on the ground.Since the transponder delay does not affect the Doppler measurements we were therefore, able to compute the precise orbit of the spacecraft without considering the transponder delay in the range measurements. With this configuration, we then reanalyzed the entire radio tracking data (see Table 1). To check the precision on the knowledge of the spacecraft orbit, we compared the radial, alongtrack, and crosstrack components of the orbit for each data arc with the solution obtained in Sect. 2.2.1. An average rms value of radial, alongtrack, and crosstrack difference is 0.015 m, 0.16 m, and 0.19 m, respectively. Less than a meter level of differences in the orbit implies that the transponder delay has negligible impact on the orbit, since the spacecraft orbit is mostly constrained by the Doppler tracking data. However, there is a dramatic change in the estimation of range bias, which now includes ephemeris bias plus the bias due to the transponder delay. Using these range biases to fit the planetary ephemeris, we found a clear offset in the EarthMercury geocentric distances of about 410 ± 20 m (twoway) during the orbital period of the . This estimation of transponder delay is compatible with the one found during ground testing, which ranged from 1,356.89 ns (~407 m) to 1,383.74 ns (~415 m) (Srinivasan et al. 2007). Thus these results also suggested that there is not a large error in the spacecraft and the planetary orbit fit procedure.
Asteroid masses obtained with INPOP13a, significantly different from values found in INPOP10e, and inducing a change in the Earthplanets distances smaller than 5 meters over the fitting interval.
3. Improvement of planetary ephemeris, INPOP
Since 2003, INPOP planetary ephemerides have been built on a regular basis and provided to users thought the IMCCE website^{7}. The INPOP10e ephemeris was the latest release (Fienga et al. 2013) that was delivered as the official Gaia mission planetary ephemerides used for the navigation of the satellite as well as for the analysis of the data. Specific developments and analysis were done for the Gaia release such as the TCB timescale version or an accurate estimation of the INPOP link to ICRF. With the delivery of the MESSENGER radio science data, a new opportunity was offered to improve drastically our knowledge of the orbit of Mercury and to perform tests of gravity at a close distance from the Sun.
The use of the 1.5 year range measurements deduced from the previous analysis (see Sect. 2) is then a crucial chance to obtain better knowledge over the ~0.3 year Mercury orbit. The accuracy of a few meter for the MESSENGER range data will give big constraints over short period perturbations on Mercury’s orbit. The five flyby positions obtained with Mariner and the MESSENGER flybys will still be significant for the measurements of long period (10 or more years) perturbations (see Fig. 8). Only the addition of the BepiColombo range data will be able to disentangle such long period effects.
3.1. INPOP13a
The constants and dynamical modeling used for constructing the new ephemerides, INPOP13a, are similar to INPOP10e. A complete adjustment of the planet initial conditions (including Pluto and the Moon), the mass of the Sun, the oblateness of the Sun, the ratio between the mass of the Earth and the Moon, and 62 asteroid masses is operated. Values of the obtained parameters are given in Tables 2 and 3. Even if Mercury is not directly affected by the main belt asteroids, the use of the range measurements between MESSENGER and the Earth does have an impact on the Earth’s orbit and then could bring some information on asteroid masses perturbing the Earth orbit. On Table 3, we only gave the masses that are significantly different from those obtained with INPOP10e and inducing detectable signatures below five meters.These masses are also compared with the Konopliv et al. (2011) on the same table. The masses of the biggest objects differ from the two ephemerides inside their twosigma error bars, and one can notice the new determination of the mass of (51) Nemausa inducing slightly bigger perturbations on Mercury (7 m) and Venus (8 m) geocentric distances than on Mars (5 m).
Fig. 3 MESSENGER oneway range residuals obtained with INPOP13a, INPOP10a, DE423, and DE430. 

Open with DEXTER 
Table 4 gives the postfit residuals obtained with INPOP13a and compared with those obtained with INPOP10e. One can see a noticeable improvement in Mercury’s orbit over all the periods of the fit including direct radar observations. The result is of course more striking for MESSENGER range measurements that were deduced from Sect. 2, and not used for the fit of INPOP10e. In this particular case, the improvement reaches a factor of almost 16 on the estimation of the distance between Mercury and the Earth (see Fig. 3). The extrapolated residuals given in Table 4 are not really significant since INPOP10e was fitted over a very similar interval of time ending at about 2010.4 when INPOP13a was fitted up to 2011.4.
Figure 4 plots the differences between INPOP13a, INPOP10e and DE423 for Mercury geocentric right ascension, declination, and distance, and the EarthMoon barycenter longitudes, latitudes, and distances in the BCRS. These differences give estimations of the internal accuracy of INPOP13a. By comparison, the same differences between INPOP10a and DE421 are also plotted. They present the improvements reached since INPOP10a, clearly noticeable for the Mercury geocentric distances (a factor two between INPOP13aDE423 and INPOP10aDE421). They are less impressive for the EMB; however, one can notice that the clear systematic trend in the INPOP10aDE423 barycentric distances of the EMB is removed in INPOP13aDE423. The fact that the differences between INPOP13a and INPOP10e are smaller than the differences to DE ephemerides is mainly discussed in Fienga et al. (2011) and Fienga et al. (2013) by a different method of computing the orbit of the Sun relative to the solar system barycenter, as well as a different distribution of planetary and asteroid masses.
In conclusion, INPOP13a shows an important improvement in the Mercury orbit especially during the MESSENGER orbital and flyby phases of the mission. The improvement over the EMB orbit in the BCRS is less important but still a systematic trend noticeable in the EMB barycentric distance differences between INPOP10a and DE421 seems to be removed in the new comparisons.
Statistics of the residuals obtained after the INPOP13a fit.
Fig. 4 In Mercury panel, differences in geocentric Mercury right ascension (RA), declination (DE) and distances between INPOP13a, INPOP10e and DE423. In EMB panel, differences in BCRF longitudes, latitudes and distances of the EMB between INPOP13a, INPOP10e and DE423. Differences between INPOP10a and DE421 are also given. 

Open with DEXTER 
3.2. Reconstruction of MESSENGER orbit with INPOP13a
As given in Table 4, geometric distances between Earth and Mercury are ~16 times better in INPOP13a than the INPOP10e. To analyze the impact of the improvement of the planetary ephemeris on the spacecraft orbit, we reanalyzed the entire one and half years of radioscience data (see Table 1) using INPOP13a ephemeris.
Fig. 5 Comparison between INPOP13a and INPOP10e estimations of MESSENGER orbit: a) differences in the postfit Doppler residuals; b) differences in the postfit range residuals; c) differences in the periapsis altitude δp; d) differences in the apoapsis altitude δa. 

Open with DEXTER 
The dynamical modeling and orbit determination process for this analysis are the same as discussed in Sect. 2.1. To compare the results of this analysis with the one obtained from INPOP10e (see Sect. 2.2.1), the differences in the Doppler and range postfit residuals along with the changes that occurred in the periapsis and apoapsis altitudes of are plotted in Fig. 5.
An average value of these differences and its 1σ mean dispersion for Doppler, and range postfit residuals was estimated as 0.008 ± 0.04 mHz and 0.05 ± 0.3 m, respectively. These values are far below the estimated accuracy of 4.8 ± 2.2 mHz and 1.9 ± 1.4 m (see Sect. 2.2.1) for Doppler and range postfit residuals, respectively. In addition to these residuals, we also compared the orbit of computed with INPOP13a and INPOP10e ephemerides. The differences in the periapsis δp and apoapsis δa altitudes of due to the change in planetary ephemeris are plotted in panels c and d of Fig. 5. An average and 1σ dispersion of δp and δa was found as 0.05 ± 1.2 m and 0.03 ± 1.2 m, respectively. These values are also far below the required accuracy of 10 m (Srinivasan et al. 2007) for the orbit. This analysis is therefore consistent with the fact that change in the planetary ephemeris during the construction of the spacecraft orbit does not alter the radioscience analysis significantly.
4. Test of general relativity
4.1. General presentation
INPOP13a was built in the framework of in using the PPN formalism. A detailed description of the modeling used for the INPOP ephemerides is given in (Fienga et al. 2011). Specific relativistic timescales have been set up and integrated in INPOP (TCB, TDB) and mainly two parameters characterized the relativistic formalism in modern planetary ephemerides: the parameter β that measures the nonlinearity of gravity and γ, measuring the deflexion of light. In GR, both are supposed to be equal to 1 and were fixed to 1 for the INPOP13a construction. The GINS software used for the analysis of the radio science data and the reconstruction of the MESSENGER orbit is also coded in the PPN framework, including both β and γ PPN parameters.
Up to now, general relativity theory (GRT) has successfully described all available observations, and no clear observational evidence against has been identified. However, the discovery of Dark Energy which challenges GRT as a complete model for the macroscopic universe, and the continuing failure to merge GRT and quantum physics indicate that new physical ideas should be sought. To streamline this search it is indispensable to test GRT in all accessible regimes and to the highest possible accuracy.
Among all possibilities for testing GRT, the tests of the motion and light propagation in the solar system were historically the first ones, and they are still very important since they give highest accuracy since the dynamics of the solar system is well understood and supported by a long history of observational data. Concerning the Einstein field equations, the most important framework used for the tests in the solar system is the PPN formalism (such as Will 1993). The PPN formalism is a phenomenological scheme with ten dimensionless parameters covering certain class of metric theories of gravity, among them the β and γ parameters parts of the INPOP and GINS modelings. The tracking data of space missions give a good possibility to test GRT since the data is very sensitive to the GRTeffects in both dynamics of the spacecraft and signal propagation. However, some factors, such as navigation unknowns (AMDs, solar panel calibrations), planet unknowns (potential, rotation, etc.), effect of the solar plasma, or the correlation with planetary ephemerides limit this sort of gravity test. Dynamics of the solar system are, however, less affected by poorly modeled accelerations and technical unknowns. Up to now, the best constraints for β come from the planetary data in INPOP (Fienga et al. 2011). Constraints on other PPN parameters can be found in Will (2006). A number of theoretical models predict deviations of PPN parameters that are smaller than current constraints. Typical examples here are certain types of tensorscalar theories where cosmological evolution exhibits an attractor mechanism towards GRT (Damour & Nordtvedt 1993) or stringinspired scalartensor theories where the scalar field can decouple from matter (Damour & Polyakov 1994).
Another phenomenological test concerns the constancy of the Newtonian gravitational constant G in time. A variable G is produced say by alternative theories of gravity such tensorscalar theory (see e.g. Damour et al. 1990 and Uzan 2003) or some models of dark energy (Steinhardt & Wesley 2009; Alimi & Füzfa 2010). The ratio is now constrained at the level of 10^{13} with LLR analysis (Williams et al. 2004).
4.2. Estimation of PPN parameters, γ and β
In this section, we propose to use the improvement of Mercury’s orbit as an efficient tool for testing the consistency between planetary ephemerides built with MESSENGER radio science data and nonunity PPN parameters.
4.2.1. Method
A first estimation of PPN β and γ is possible by least square methods during the adjustment of the INPOP planetary ephemerides, and the results are given in Table 5. Figure 6 gives the correlations between the first 71 over the 343 parameters estimated in the adjustments. As one can see in Fig. 6 no correlation greater than 0.3 affects the determination of the PPN parameters β and γ, as well as the fit of the Sun oblateness, when the gravitational mass of the Sun is highly related to the Mercury and to the Earth orbits.
However to go further in the analysis of the uncertainties and the construction of acceptable intervals of violation of GR throught the PPN β and γ, we also considered the same method as the one that was used and described in Fienga et al. (2011) for determining acceptable intervals of violation of when the PPN formalism. Small variations in these two parameters near unity are imposed when constructing alternative planetary ephemerides that are fit over the whole data sets presented in Table 4 and with the same parameters and hypothesis as INPOP13a. A minimum of three iterations in the adjustment process is required for building new ephemerides, and comparisons between these ephemerides and INPOP13a are done to scale up what variations to are acceptable at the level of uncertainty of the present planetary ephemerides.
Fig. 6 Correlation between the first 71 (over 343) parameters estimated during the fit of the planetary ephemerides. The red frame frames the correlations related to the initial conditions of planet orbits and the blue rectangle frames the correlations related to the Sun J_{2} (JS), the PPN parameters β (BE) and γ (GA) and the gravitational mass of the Sun (GS). The magenta rectangle frames the correlations related to the gravitational masses of the first most perturbing asteroids including the gravitational mass of the asteroid ring (GR). m_{1},...m_{6} expresses the initial conditions of the Mercury orbit in equinoctial coordinates: semimajor axis, mean motion, k, h, q, and p respectively. The other planet initial conditions are indicated by the first letter of the planet (V for Venus, M for Mars etc...) and by the figures of the corresponding initial conditions as given for Mercury. 

Open with DEXTER 
Fig. 7 Variations in postfit residuals obtained for different values of PPN β (xaxis) and γ (yaxis). Panels a) and c) are obtained by considering the variations in the whole data sets when for panels b) and d), variations in the Mercury flyby data (from Mariner and MESSENGER missions) are excluded from the analysis. The dashed line indicates the limit in γ given by Bertotti et al. (2003). 

Open with DEXTER 
The improvement of Mercury’s orbit in INPOP13a justifies these new estimations. Indeed, Mercury played a historical role in testing gravity and in 1912 (Einstein 1912) and it is still the planet the most influenced by the gravitational potential of the Sun. Its orbit can then lead to the most efficient constraints on β, hence on γ in the PPN formalism. Before the recent input of MESSENGER flyby and orbital radio science data in the INPOP construction, Mars was the most constraining planet for the PPN parameters (Fienga et al. 2010). The reason was the long range of high accurate observations on Mars. The implementation of the first MESSENGER flyby data reduces the interval of violation of β to 50%. The first estimation of the γ interval of violation was made possible thanks to the gain in uncertainty on the Mercury orbit. With INPOP13a, even better improvement is achieved.
4.2.2. Results
The results obtained by direct least squares fit are presented in Table 5. As expected, the estimated uncertainties are very optimistic and a more detailed analysis is done based on the method proposed by Fienga et al. (2011).
Results obtained in terms of percentages of the variations in postfit residuals between a planetary ephemeris fitted and built up with PPN parameters different from one and INPOP13a are given in Fig. 7. Panel a in Fig. 7 gives the map of the variations in percent of the full dataset postfit residuals. Panel b in Fig. 7 gives the same map but without taking the variations of the Mercury flyby data into account. For panel b, the Mercury flyby data are indeed used in the ephemerides fit but not in the analysis of the postfit residuals for testing GR. The map of panels a and b is then dramatically different: where the limits for β and γ are stringent for the map including the Mariner data, the constraints are greatly enlarged for these two parameters. These phenomena were expected since the variations in PPN parameters induce longterm perturbations in the geocentric distances of Mercury as one can see in Fig. 8. Panels c and d are similar to a and b, but they are obtained with ephemeris INPOP10a. In this ephemeris, MESSENGER flyby data were included in its fits but not in the orbital data. By comparing panels a and c, one can see that the use of the MESSENGER orbital data significantly reduce the intervals of violation for both PPN parameters by a factor 10. The same manner, in the most pessimistic case and without considering the Mercury flybys in analysing of the variations in the postfit residual, one can see in panels b and d that the improvement of Mercury’s orbit is again crucial for reducting the violation intervals of PPN parameters.
Table 5 collects the acceptable violation intervals obtained from INPOP10a and INPOP13a. Values extracted from INPOP10a were obtained at 5% of postfit residual variations (Fienga et al. 2011).With INPOP13a, we extracted values from i) panel a of Fig. 7 obtained at 5%, but also at 10% and 25%; and ii) from panel b of Fig. 7 obtained at 5%, which is consistent with the 25% of intervals extracted from panel a.
All given intervals are compatible with GR with an uncertainty at least ten times smaller than our previous results with INPOP10a. In Table 5, comparisons to least squares estimations of other planetary ephemerides or Moon ephemerides like Pitjeva & Pitjev (2013), Konopliv et al. (2011), Müller et al. (2008), and Williams et al. (2009), as well as estimations deduced from VLBI observations Lambert & Le PoncinLafitte (2011), are also given. The most stringent published constraint for the PPN parameter γ has been obtained so far during a dedicated phase of the Cassini mission by Bertotti et al. (2003). This value is compatible with our 25% estimation when our 5% and 10% estimations give more restrictive intervals of GR violations.
Confirmations of the results presented in Table 5 will be obtained by the use of the radioscience data obtained during the future BepiColombo mission. In addition, the recovery of the Mariner flyby data would also be a great help for such confirmations. Unfortunately, the Mariner data seem to have been lost, and access to these data seems to be unrealistic. Indeed as one can see in Fig. 8, perturbations induced by a slight change in the parameters ((β − 1) = 1.5 × 10^{5} and (γ − 1) = − 2.2 × 10^{5}) inducing an effect of about six meters on the Mariner range data (12%) will induce a signature of about the same level at the BepiColombo epoch. With the improved BepiColombo radio science tracking, the expected accuracy in the range measurement is planned to be about 50 centimeters. With such accuracy, detecting the perturbations induced by the same modification of the PPN parameters should be done at 1200%! Two orders of magnitude are expected as a gain in the uncertainty for the β and γ estimations.
Fig. 8 Differences in geocentric Mercury distances between INPOP13a and a planetary ephemerides built with PPN β and γ different from 1. The indicated area shows intervals of time corresponding to Mariner observations, MESSENGER and the future BepiColombo. 

Open with DEXTER 
Fig. 9 Differences between solvefor parameters obtained from the solutions SOL_{Ref} and SOL_{βγ ≠ 1} (see text): panels a)–f) represent the changes in the initial state vectors, panel g) represents the changes in the scale factors estimated for solar radiation pressure, panel h) represents the changes in the estimated range bias. The error bars represent the 1σ uncertainties in the estimation of solvefor parameters obtained with the reference solution SOL_{Ref}. 

Open with DEXTER 
4.2.3. Impact on MESSENGER orbit
As stated previously, the software was modeled in the framework of the formalism which includes β and γ parameters. To analyze the combined impact of parameters over the orbit and the planetary ephemerides construction, we analyzed the entire one and half years of radioscience data again using the parameters that are different from unity. The same procedure as described in Sect. 2.1 has been used for reconstructing the orbit. Two sets of orbits were then built, one with β and γ equal to unity and the other with β and γ different from unity (in this case, β − 1 = γ − 1 = 1 × 10^{4}). Hereafter, the solution obtained from β and γ equal to unity is referred to as SOL_{Ref} (for the reference solution), and the solution corresponds to β and γ different from unity referred to as SOL_{βγ ≠ 1}.
To maintain consistency in constructing the orbit, we used corresponding planetary ephemerides that were built with the same configurations of parameters as used for SOL_{Ref} and SOL_{βγ ≠ 1}. In Fig. 9, we plotted the differences between solvefor parameters obtained for SOL_{Ref} and SOL_{βγ ≠ 1}. The error bars shown in the same figure represent the 1σ uncertainties in the estimation of solvefor parameters corresponding to SOL_{Ref}. From this figure, one can notice the differences in the parameters are always below the 1σ uncertainties. The estimated solvefor parameters for SOL_{βγ ≠ 1} are analogous to SOL_{Ref}, and there is no significant change in the orbit due to the change in parameters. In contrast, as shown in Fig. 7, this configuration of parameters (β − 1 = γ − 1 = 1 × 10^{4}) in the construction of planetary ephemerides led to ~65% of change in the postfit residuals, which shows that, the planetary ephemerides are more sensitive to GR effects. This can be explained from the fitting intervals of the data set. Usually planetary ephemerides are fitted over long intervals of times (see Table 4) to exhibit longterm effects, while a spacecraft orbit is usually constructed over much shorter intervals (usually one day to a few days) of data arcs to account for the model’s imperfections. The short fitting interval of the spacecraft orbit would absorb such effects.
Moreover, it is worth noticing that, unlike state vectors and scale factor F_{S} (see panels a–g of Fig. 9), the range bias differences between SOL_{Ref} and SOL_{βγ ≠ 1} solutions (see panel h) shows systematic behavior. This trend in the range bias can be explained from the contribution of the relativistic deflection of light by the Sun (a function of parameter γ, Shapiro 1964) in the light time computations. Explicitly this effect was not absorbed during the computation of range bias and it becomes important to examine this effect when constructing planetary ephemerides.
We, therefore, reconstruct the planetary ephemerides using the range bias obtained from SOL_{βγ ≠ 1} and the parameters β − 1 = γ − 1 = 1 × 10^{4}. The newly estimated postfit range bias is then compared with the range bias (prefit) corresponding to SOL_{βγ ≠ 1}. This investigation shows that, the postfit residuals are modified by ~6% for and ~1% for Mariner 10 with respect to prefit residuals. This modification in the residuals is negligible compared to a ~65% of change with respect to reference residuals obtained from INPOP13a. As a result, the supplementary contributions in the range bias due to the relativistic deflection between station and did not bring any significant change in the planetary ephemerides construction.
5. Conclusions
We analyzed one and half years of radioscience data of the spacecraft using orbit determination software GINS. An accurate orbit of was then constructed with the typical range of Doppler, and twoway range residuals of about 4.8 ± 2.2 mHz (~0.09 ± 0.04 mm/s), and 1.9 ± 1.4 m. Such accuracies are comparable to those in Smith et al. (2012); Stanbridge et al. (2011); Srinivasan et al. (2007). Range measurements obtained by the MESSENGER spacecraft during its mapping period were then used to construct improved planetary ephemerides called INPOP13a. This ephemeris showed an accuracy of about − 0.4 ± 8.4 m in the MercuryEarth geometric distances, which is two orders of improvement compared to DE423 and INPOP10e, and one order compared to the latest DE430.
Such high precision Mercury ephemeris allowed us to perform one of the best tests of formalism. To determine the acceptable intervals of the violation of through the parameters (β, γ), small variations of these two parameters near unity were imposed in the construction of alternative planetary ephemerides fitted over the whole data sets. The percentage difference between these ephemerides to INPOP13a are then used to defined the interval of parameters β and γ.
As expected, our estimations of PPN parameters are more stringent than previous results. We considered the 5%, 10% and 25% of changes in the postfit residuals. That the intervals correspond to these changes is compatible with with an uncertainty at least ten times smaller than our previous results with INPOP10a. Moreover, one of the best estimation of parameter γ has so far been estimated from the Cassini observations by Bertotti et al. (2003), which is compatible with our 25% estimation.
To further the accuracy of the parameters improve, and to confirm the results given in Table 5, one needs to analysis the radioscience data of the future BepiColombo mission.
Acknowledgments
We are very thankful to the CNES and Région FrancheComté, who gave us financial support. Part of this work useds GINS software, so we would like to acknowledge the CNES for providing access to this software. A. K. Verma is thankful to P. Rosenbatt and S. Le Maistre for fruitful discussions. We are also thankful to G.EspositoFarese for his constructive remarks and comments.
References
 Alimi, J.M., & Füzfa, A. 2010, in AIP Conf. Ser., 1241, 690 [Google Scholar]
 Anderson, J. D., Colombo, G., Espsitio, P. B., Lau, E. L., & Trager, G. B. 1987, Icarus, 71, 337 [NASA ADS] [CrossRef] [Google Scholar]
 Anderson, J. D., Jurgens, R. F., Lau, E. L., Slade, III, M. A., & Schubert, G. 1996, Icarus, 124, 690 [NASA ADS] [CrossRef] [Google Scholar]
 Anderson, J. D., Turyshev, S. G., Asmar, S. W., et al. 1997, Planet. Space Sci., 45, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Bertotti, B., Iess, L., & Tortora, P. 2003, Nature, 425, 374 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Damour, T., & Nordtvedt, K. 1993, Phys. Rev. D, 48, 3436 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Damour, T., & Polyakov, A. M. 1994, Gen. Relativ. Gravit., 26, 1171 [NASA ADS] [CrossRef] [Google Scholar]
 Damour, T., Gibbons, G. W., & Gundlach, C. 1990, Phys. Rev. Lett., 64, 123 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Ercol, C. J., M., B. S., & Holtzman, G. A. 2012, in Thermal and Fluids Analysis Workshop, TFAWS2012PT03 [Google Scholar]
 Fienga, A., Manche, H., Laskar, J., & Gastineau, M. 2008, A&A, 477, 315 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fienga, A., Laskar, J., Morley, T., et al. 2009, A&A, 507, 1675 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fienga, A., Laskar, J., Kuchynka, P., et al. 2010, in IAU Symp. 261, eds. S. A. Klioner, P. K. Seidelmann, & M. H. Soffel, 159 [Google Scholar]
 Fienga, A., Laskar, J., Kuchynka, P., et al. 2011, Celest. Mech. Dyn. Astron., 111, 363 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Fienga, A., Manche, H., Laskar, J., Gastineau, M., & Verma, A. 2013 [arXiv:1301.1510] [Google Scholar]
 Flanigan, S. H., OShaughnessy, D. J., Wilson, M. N., & Hill, T. A. 2013, in 23rd Space Flight Mechanics Meeting, AAS, 13 [Google Scholar]
 Folkner, W. M. 2010, JPL, interoffice memorandum, IOM 343R10001 [Google Scholar]
 Konopliv, A. S., Asmar, S. W., Folkner, W. M., et al. 2011, Icarus, 211, 401 [NASA ADS] [CrossRef] [Google Scholar]
 Lambert, S. B., & Le PoncinLafitte, C. 2011, A&A, 529, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Manche, H., Fienga, A., Laskar, J., et al. 2010, in Journées 2010 Systèmes de Référence Spatiotemporels [Google Scholar]
 Margot, J.L. 2009, Celest. Mech. Dyn. Astron., 105, 329 [NASA ADS] [CrossRef] [Google Scholar]
 Moyer, T. D. 2003, in Formulation for Observed and Computed Values of Deep Space Network Data Types for Navigation (John Wiley & Sons), Vol. 2 [Google Scholar]
 Müller, J., Soffel, M., & Klioner, S. A. 2008, J. Geodesy, 82, 133 [NASA ADS] [CrossRef] [Google Scholar]
 Pitjeva, E. V. 2009, in IAU Symp. (AAS), 261, 603 [Google Scholar]
 Pitjeva, E. V., & Pitjev, N. P. 2013, MNRAS, 432, 3431 [NASA ADS] [CrossRef] [Google Scholar]
 Shapiro, I. I. 1964, Phys. Rev. Lett., 13, 789 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Smith, D. E., Zuber, M. T., Phillips, R. J., et al. 2010, Icarus, 209, 88 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, D. E., Zuber, M. T., Phillips, R. J., et al. 2012, Science, 336, 214 [NASA ADS] [CrossRef] [Google Scholar]
 Solomon, S. C., McNutt, R. L., Gold, R. E., & Domingue, D. L. 2007, Space Sci. Rev., 131, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Srinivasan, D. K., Perry, M. E., Fielhauer, K. B., Smith, D. E., & Zuber, M. T. 2007, Space Sci. Rev., 131, 557 [NASA ADS] [CrossRef] [Google Scholar]
 Stanbridge, D. R., Williams, K. E., Taylor, A. H., et al. 2011, in AAS/AIAA Astrodynamics Specialist Conference, AAS 11548 [Google Scholar]
 Steinhardt, P. J., & Wesley, D. 2009, Phys. Rev. D, 79, 104026 [NASA ADS] [CrossRef] [Google Scholar]
 Uzan, J.P. 2003, Annales Henri Poincaré, 4, 347 [NASA ADS] [CrossRef] [Google Scholar]
 Vaughan, R. M., Haley, D. R., OShaughnessy, D. J., & Shapiro, H. S. 2002, Adv. Astronau. Sci., 109, 1139 [Google Scholar]
 Verma, A. K., Fienga, A., Laskar, J., et al. 2013, A&A, 550, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Will, C. 2001, Liv. Rev. Rel., 4, 4 [Google Scholar]
 Will, C. M. 1993, Theory and Experiment in Gravitational Physics (Cambridge University Press) [Google Scholar]
 Will, C. M. 2006, Liv. Rev. Rel., 9, 3 [Google Scholar]
 Williams, J. G., Turyshev, S. G., & Boggs, D. H. 2004, Phys. Rev. Lett., 93, 261101 [NASA ADS] [CrossRef] [Google Scholar]
 Williams, J. G., Turyshev, S. G., & Boggs, D. H. 2009, Int. J. Mod. Phys. D, 18, 1129 [NASA ADS] [CrossRef] [Google Scholar]
 Williams, J. G., Boggs, D. H., & Folkner, W. M. 2013, IOM 335JW, DB, WF20080314001 [Google Scholar]
All Tables
Values of parameters obtained in the fit of INPOP13a and INPOP10e to observations including comparisons to DE423 and DE430.
Asteroid masses obtained with INPOP13a, significantly different from values found in INPOP10e, and inducing a change in the Earthplanets distances smaller than 5 meters over the fitting interval.
All Figures
Fig. 1 Quality of the MESSENGER orbit in terms of rms values of the postfit residuals for each oneday data arc: (a) two and threeway Doppler given in millihertz (multiply by 0.0178 to obtain residuals in mm/s), and (b) twoway range given in meters. 

Open with DEXTER  
In the text 
Fig. 2 History of the fitted scale factor and estimated range bias: a) scale factor (for solar radiation acceleration) fitted over each oneday data arc to account inaccuracy in the force model, and b) oneway range bias, represent the systematic error in the EarthMercury positions, estimated for each oneday arc using INPOP10e (•), DE423 (▴), and DE430 (▾) planetary ephemerides. 

Open with DEXTER  
In the text 
Fig. 3 MESSENGER oneway range residuals obtained with INPOP13a, INPOP10a, DE423, and DE430. 

Open with DEXTER  
In the text 
Fig. 4 In Mercury panel, differences in geocentric Mercury right ascension (RA), declination (DE) and distances between INPOP13a, INPOP10e and DE423. In EMB panel, differences in BCRF longitudes, latitudes and distances of the EMB between INPOP13a, INPOP10e and DE423. Differences between INPOP10a and DE421 are also given. 

Open with DEXTER  
In the text 
Fig. 5 Comparison between INPOP13a and INPOP10e estimations of MESSENGER orbit: a) differences in the postfit Doppler residuals; b) differences in the postfit range residuals; c) differences in the periapsis altitude δp; d) differences in the apoapsis altitude δa. 

Open with DEXTER  
In the text 
Fig. 6 Correlation between the first 71 (over 343) parameters estimated during the fit of the planetary ephemerides. The red frame frames the correlations related to the initial conditions of planet orbits and the blue rectangle frames the correlations related to the Sun J_{2} (JS), the PPN parameters β (BE) and γ (GA) and the gravitational mass of the Sun (GS). The magenta rectangle frames the correlations related to the gravitational masses of the first most perturbing asteroids including the gravitational mass of the asteroid ring (GR). m_{1},...m_{6} expresses the initial conditions of the Mercury orbit in equinoctial coordinates: semimajor axis, mean motion, k, h, q, and p respectively. The other planet initial conditions are indicated by the first letter of the planet (V for Venus, M for Mars etc...) and by the figures of the corresponding initial conditions as given for Mercury. 

Open with DEXTER  
In the text 
Fig. 7 Variations in postfit residuals obtained for different values of PPN β (xaxis) and γ (yaxis). Panels a) and c) are obtained by considering the variations in the whole data sets when for panels b) and d), variations in the Mercury flyby data (from Mariner and MESSENGER missions) are excluded from the analysis. The dashed line indicates the limit in γ given by Bertotti et al. (2003). 

Open with DEXTER  
In the text 
Fig. 8 Differences in geocentric Mercury distances between INPOP13a and a planetary ephemerides built with PPN β and γ different from 1. The indicated area shows intervals of time corresponding to Mariner observations, MESSENGER and the future BepiColombo. 

Open with DEXTER  
In the text 
Fig. 9 Differences between solvefor parameters obtained from the solutions SOL_{Ref} and SOL_{βγ ≠ 1} (see text): panels a)–f) represent the changes in the initial state vectors, panel g) represents the changes in the scale factors estimated for solar radiation pressure, panel h) represents the changes in the estimated range bias. The error bars represent the 1σ uncertainties in the estimation of solvefor parameters obtained with the reference solution SOL_{Ref}. 

Open with DEXTER  
In the text 