Abstract
The NASA MESSENGER mission explored the innermost planet of the solar system and obtained a rich data set of range measurements for the determination of Mercuryâs ephemeris. Here we use these precise data collected over 7 years to estimate parameters related to general relativity and the evolution of the Sun. These results confirm the validity of the strong equivalence principle with a significantly refined uncertainty of the Nordtvedt parameter ηâ=â(â6.6â±â7.2)âÃâ10â5. By assuming a metric theory of gravitation, we retrieved the post-Newtonian parameter βâ=â1â+â(â1.6â±â1.8)âÃâ10â5 and the Sunâs gravitational oblateness, \({{J}}_{2 \odot }\)â=â(2.246â±â0.022)âÃâ10â7. Finally, we obtain an estimate of the time variation of the Sun gravitational parameter, \(\dot{{G} {{M}}_ \odot} {\mathrm{/}}{{G}}{{M}}_ \odot\)â=â(â6.13â±â1.47)âÃâ10â14, which is consistent with the expected solar mass loss due to the solar wind and interior processes. This measurement allows us to constrain \(\left| {{\dot{ G}}} \right|{\mathrm{/}}{{G}}\) to be <4âÃâ10â14 per year.
Similar content being viewed by others
Introduction
Mercuryâs role in testing theories of gravitation has always been crucial because the strong gravitational mass of the Sun creates notable perturbations on its orbit. The precession of the closest distance of Mercury to the Sun (perihelion) first highlighted the limits of Newtonian physics and later validated the predictions of Einsteinâs theory of general relativity (GR)1. The precession of Mercuryâs perihelion is primarily due to third-body perturbations from other planets (~531.63â per Julian century2), and the relativity effect produces an additional perihelion shift of ~42.98â per Julian century3. The relativistic corrections to Mercuryâs heliocentric acceleration can be formulated based on the parameterized post-Newtonian (PPN) parameters β and γ, which respectively measure the nonlinearity in superposition of gravity and space-time curvature produced by a unit rest mass. Both parameters are zero in the Newtonian formulation and equal to 1 in GR (Methods).
The Sunâs interior structure and dynamics also affect Mercuryâs trajectory. The solar gravitational oblateness \({{J}}_{{\it{2}}{ \odot }}\) and the angular momentum Sâ are responsible for additional precession rates of ~0.029â per Julian century4 and ~0.002â per Julian century5, respectively. The latter perturbation, which is known as the gravitomagnetic EinsteinâLenseâThirring (ELT) effect, is related to the distortion of space-time induced by the rotation of the Sun.
In practice, strong correlations between γ, β, \({{J}}_{{\it{2}}{ \odot }}\), and Sâ limit the combined estimation of these parameters since they all primarily affect Mercuryâs orbit through the precession of its perihelion. For this reason, a priori assumptions are necessary to disentangle the effect of each parameter. The Nordtvedt parameter η, related to the equivalence principle (EP), can be used as a constraint between the PPN parameters γ and β6. The relationship between these coefficients is:
if we assume spatial isotropy, which implies that the PPN parameters α1 and α2 are equal to 0.
The orbit of Mercury is well-suited to test the EP, which describes the equality between gravitational and inertial masses. The EP has been partially demonstrated by laboratory experiments, to a precision of ~1âÃâ10â13 with recent torsion-balance tests7. However, these precise results only concern the weak EP, which is based on Galileoâs postulate that different objects fall with the same acceleration in a uniform gravitational field, independent of their composition and structure. Einstein extended this concept in his development of GR by introducing the strong EP (SEP). The SEP states that a uniform gravitational field is locally indistinguishable from an accelerated reference frame8. The contribution of the SEP to the gravitational-to-inertial mass ratio depends on the self-gravitational energy (ΩB) and the rest energy of the body (mIc2), as follows:
where mG and mI are the gravitational and inertial masses, respectively, ΩB is proportional to G(mG)2Râ1, G is the gravitational constant, c is the speed of light, and R is the radius of body B, respectively. The Nordtvedt parameter η must be zero to validate the SEP. To prove the SEP, the test mass used in the experiment needs to be sufficiently large so that the self-gravitational force is not negligible. For this reason, tests at the scale of the solar or planetary system are suitable to prove the SEP.
The most accurate estimations of η have been retrieved from lunar laser ranging (LLR) over the past 40 years9,10,11,12. The latest solution validates the SEP with an uncertainty of Ï Î· ~3.0âÃâ10â4 (Table 1). The coupling of the gravitational attraction of the Sun on the EarthâMoon system with the self-gravitational force of the Earth would provide a significant perturbation in the case of SEP violation. This effect would be measurable with LLR mm-precision data of the EarthâMoon distance12.
An equivalent dynamical effect on Mercuryâs orbit is due to the coupling between the Sunâs self-gravitational force and the gravitational attraction of other planets, mainly Jupiter. However, the main effect that a SEP violation has on the ephemeris of Mercury results from the implied redefinition of the solar system barycenter (SSB), which is negligible in the EarthâMoon case (Methods). A Nordtvedt parameter η of 1âÃâ10â5 results in discrepancies in the MercuryâEarth relative distance of ~3âm after 2 years13. Thus, the knowledge of Mercuryâs ephemeris to better than 1âm can yield better constraints on possible SEP violations than LLR. Furthermore, this dynamical perturbation of the Nordtvedt parameter is less correlated with other forces and thus separates the effects of \({ {J}}_{{\it{2}}{\it{ \odot }}}\) and β, given the constraint of Eq. 1.
The study of Mercuryâs orbit with a long-duration data set also gives a unique opportunity to detect the time variation of the gravitational constant G. The estimation of \(\dot{{G}}/{{G}}\) is not strongly correlated with other relativistic and solar parameters because its effect is quadratic in time. However, Mercuryâs orbit is perturbed by the combined effect of secular changes in G and Mâ as follows:
where \(\dot { {M}}_ \odot {\mathrm{/}}{ {M}}_ \odot\) is the Sunâs mass loss due to solar radiance and wind. A perturbation on Mercuryâs orbit induced by a \({ {G}}\dot { {M}}_ \odot /{ {G}}{ {M}}_ \odot\) of 5âÃâ10â14, which is ~10% of the Sunâs expected mass loss14, is on the order of ~2âm after 2 years, when projected on the EarthâMercury line-of-sight. An estimated time variation of \({ {G}}{ {M}}_ \odot\) combined with a \(\dot { {M}}_ \odot {\mathrm{/}}{ {M}}_ \odot\) value from heliophysics studies improves the knowledge of \(\dot{ {G}}/{ {G}}\). Such a study of heliophysics and relativity requires precise observations of Mercuryâs position and velocity over an extended period of time.
In this study, we focused on the radio science data of the NASA Mercury Surface, Space Environment, Geochemistry, and Ranging (MESSENGER) mission to investigate the interior structure and evolution of the Sun (GMâ, \({ {J}}_{{\it{2}} \odot }\), and \(\dot { {M}}_ \odot {\mathrm{/}}{ {M}}_ \odot\)) and theory of gravitation (β, η, and \(\dot{ {G}}/{ {G}}\)). Our results show improved estimates of the solar gravitational oblateness and the mass loss rate that are consistent with helioseismology and heliophysics theoretical studies, respectively. The accurate measurement of the time variation of the solar gravitational parameter enabled us to constrain \(\left| \dot { {G}} \right|/{ {G}}\) to be lower than 4.0âÃâ10â14 per year. Furthermore, we determined the Nordtvedt parameter η with a refined uncertainty that demonstrates that there are no violations of the SEP at the level of ~6â7âÃâ10â5.
Results
MESSENGER and Mercury-combined orbit determination
The MESSENGER mission collected spacecraft radio tracking data near Mercury between January 2008 and April 2015, which are well-suited to improve the ephemeris of the planet15. These data are range-rate (or Doppler) observables that measure the relative velocity in the line-of-sight between the spacecraft and a Deep Space Network (DSN) Earth station, and range observables of the relative distance between the spacecraft and the DSN station.
Doppler observables have been used extensively to determine the trajectory of spacecraft for navigation and geophysical parameter estimation, e.g., the gravitational field of Mercury16. On the other hand, range observables bear on the knowledge of Mercuryâs orbit. Range measurements have been analyzed by the Solar System Dynamics Group of the Jet Propulsion Laboratory (JPL), Institute of Applied Astronomy of the Russian Academy of Science, and the Institut de mécanique céleste et de calcul des éphémérides to determine the ephemeris of Mercury, estimating relativistic and heliophysics perturbative forces17,18,19,20.
Parallel and independent investigations so far have been conducted to exclusively determine either Mercuryâs geophysics or its ephemeris. The estimation of Mercuryâs gravity field relied on the assumption of planetâs orbits and GMs, including Mercury, from JPL development ephemeris (DE). On the other hand, the ephemeris work processed spacecraft range measurements only by using a pre-converged MESSENGER trajectory. Although both methods have successfully been used for interplanetary orbit determination, their piecemeal combination is not the best approach in the case of Mercury. Systematic errors in the MESSENGER orbits directly enter the range data, and are then absorbed into Mercuryâs estimated position, since the spacecraft trajectory is not adjusted in the ephemeris solution. Conversely, a mismodeling of Mercuryâs ephemeris leads to imperfect geophysical solutions.
We numerically integrate the spacecraft and planet orbits simultaneously in order to provide a comprehensive solution that includes geophysical, heliophysics, and relativity results together with their associated covariances. Here we focus on the results that provide new information on the interior of the Sun and on gravitational theories.
The accuracy of the heliophysics and relativity results largely depends on the precision of the range data. Measurement noise is incurred by the electronics of the radio frequency (RF) telecommunications subsystem onboard the spacecraft that relays spacecraft telemetry and performs as a radio science instrument21. The MESSENGER RF subsystem operated at X-band frequencies (7.2âGHz uplink, 8.4âGHz downlink) and its two opposite-viewing phased-array antennas (PAAs) were used to conduct the range data campaigns (Methods).
Figure 1 shows the level of noise of the range data collected over the entire MESSENGER mission. Each point represents the range data RMS (root mean square) during each full tracking pass, which usually provides one measurement every 5âmin for several hours, as a function of the SunâProbeâEarth (SPE) angle. The relative position of the Sun and Earth during MESSENGER radio observations strongly controls the tracking data quality. The solar plasma causes phase scintillations in the RF signal, increasing the noise of both range-rate and range measurements. At low SPE angles (near superior solar conjunctions), the Sun is located between Mercury and the Earth, and the mean level of noise increases from <0.5âm for 90°â<âSPEâ<â180° to ~1.5âm for 35°â<âSPEâ<â90° (Fig. 1). MESSENGER orbits close to superior solar conjunctions (SPEâ<â35°) are not included in the solution because of plasma-induced range errors higher than 3âm. The exclusion of these data does not degrade the phase sampling of Mercuryâs orbit over the full MESSENGER mission. Figure 2 shows a histogram of the number of the processed measurements vs. the Mercury orbital phase. A great portion of data was collected in proximity to Mercuryâs perihelion and aphelion, enabling precise measurements of the precession induced by solar and relativistic effects. The rest of Mercuryâs orbit was evenly sampled, and the level of noise is sufficiently uniform for all orbital phases.
Noise level of the MESSENGER range data. RMS of range measurements as a function of the SunâProbeâEarth angle, which illustrates the effect of the solar plasma on the data noise. Lower SPE angles produce higher noise since the signal passes through dense solar plasma closer to the Sun. The data collected near superior solar conjunction (SPEâ<â35°) were not included in the analysis. The figure also shows the antennas that were used to provide the downlink to the DSN station. The range data were always collected during tracking passes with fanbeam for uplink and PPAs for downlink reducing thermal noise effects
Data distribution throughout Mercuryâs orbit. Number of the analyzed measurements as function of the Mercury distance from the Sun in AU. Colors indicate the noise level distribution during each phase bin of Mercuryâs orbit. The greater part of the data was collected close to Mercuryâs perihelion and aphelion
This MESSENGER data set is used to determine at the same time the solar and relativistic parameters that provide fundamental information on the interior structure of the Sun (GMâ, \({ {J}}_{{\it{2}} \odot }\), Sâ, and \(\dot { {M}}_ \odot /{ {M}}_ \odot\)) and theory of gravitation (γ, β, η, and \(\dot{ {G}}/{ {G}}\)). Multiple separate experiments have previously established a comprehensive survey of these fundamental physics effects. Table 1 summarizes recent estimates of these parameters. Helioseismology studies22,23 enabled precise measurements of Sâ and \({ {J}}_{{\it{2}} \odot }\), which have been also recovered with ephemerides analysis17,18,19,20. Planetary ephemeris investigations, furthermore, provided the best estimates for β20,24,25 and \(\dot{ {G} { {M}}}_ \odot /{ {G}}{ {M}}_ \odot\)20,25. Several of these studies do not include the ELT effect, so their estimated \({ {J}}_{{\it{2}} \odot }\) must be scaled. A more recent determination of Mercuryâs ephemeris reported the estimation of \({ {J}}_{{\it{2}} \odot }\) by accounting for ELT accelerations24. LLR provided accurate estimates of η and \(\dot{ {G}}/{ {G}}\)9,10,11,12, which has also been determined by astrophysical studies26. The Cassini mission achieved the most precise measurement of the PPN parameter γ through the analysis of radio tracking data near superior solar conjunction27. Although the MESSENGER data are not strongly sensitive to γ and the Sunâs angular momentum, Sâ, this investigation provides a unique opportunity to simultaneously improve the knowledge of GMâ, \({ {J}}_{{\it{2}} \odot }\), β, η, and \({ {G}}\dot { {M}}_ \odot /{ {G}}{ {M}}_ \odot\).
Heliophysics and relativity solutions
The main effect of γ on the radio tracking data is the deflection and delay of photons by the curvature of space time produced by the Sun, and it is best measured at SPE angles lower than 10°. However, this geometry yields high plasma noise so the estimation of the PPN parameter γ with MESSENGER data is not feasible. Given the γ-effect on Mercuryâs dynamical equations cannot be separated from the dynamical perturbation due to the PPN parameter β nor ignored, γ was thus fixed to 1, while still considering the uncertainty obtained by Cassini27.
The strong correlation between Mercuryâs orbital perturbations due to the Sunâs gravitational oblateness and the ELT effect also does not allow the determination of the angular momentum of the Sun with the MESSENGER radio tracking data. The a priori value of Sâ adopted in this study is 190âÃâ1039âkgâm2âsâ1. A covariance analysis shows that the MESSENGER data sensitivity to the ELT effect yields \(\sigma _{{\it{S}}_{\it{ \odot }}}\)~40âÃâ1039âkgâm2âsâ1, which is ~30 times larger than the current best knowledge (Table 1). For this reason, the angular momentum of the Sun is not adjusted, but the ELT effect is of course included in the integration of both the Mercury and MESSENGER orbits. The ELT effect was not modeled in the JPL DE432 ephemerides18.
Our results are based on a global combined estimation of MESSENGER- and Mercury-related orbital dynamics (Methods). Table 2 shows the a priori and estimated values and uncertainties of the heliophysics and relativistic parameters. The Sunâs GM and J2 estimates are in good agreement with previous works based on Mercuryâs ephemeris analysis20,24,25. The solar gravitational flattening is notably improved and consistent with helioseismology results, which were based on solar internal rotation measurements22,23. By applying Eq. 1 as constraint, we assume a metric theory of gravitation. The Nordtvedt relation enables a highly accurate recovery of \({ {J}}_{{\it{2}} \odot }\) and β leading to a formal uncertainty of the gravitational flattening refined by, at least a factor of 3 compared to previous ephemeris studies20,24,25. However, the correlation between \({ {J}}_{{\it{2}} \odot }\) and β is still high (~0.9, Supplementary Table 1) because the estimation of η is limited by the accuracy of the range data (Methods). Four different cases were studied to assess the effects of a priori knowledge or constraints, if we do not assume a metric theory of gravitation or if we assume that βâ1, η, or both parameters are equal to 0. These tests generalize our results further and are shown in Supplementary Table 2. The Nordtvedt equation significantly benefits the estimation of β and \({ {J}}_{{\it{2}}_{\it{ \odot }}}\), but the η and \({ {G}}\dot { {M}}_ \odot /{ {G}}{ {M}}_ \odot\) estimates are always stable and near the values shown in Table 2. We note that an unconstrained solution yields a near-unity β-\({ {J}}_{{\it{2}}{\it{ \odot }}}\) correlation and values of βâ1â=â(â1.43â±â1.47)âÃâ10â4 and \({ {J}}_{{\it{2}}{\it{ \odot }}}\)â=â(2.10â±â0.15)âÃâ10â7. In case we do not adjust for β and η, the Sunâs gravitational oblateness converges to (2.271â±â0.003)âÃâ10â7 that is still within ~1-Ï of the constrained solution.
Both constrained and unconstrained solutions are consistent with Einsteinâs theory of GR. GR predictions of β and η values are within 1-Ï, as reported in Table 2. These results enable substantial enhancement of both βâ1 and η estimates, which are ~7 and ~5 times closer to 0 than LLR solutions, respectively. The knowledge of the PPN parameter β in this study is now comparable to the Cassini solution of the PPN parameter, γ.
Furthermore, Table 2 reports the estimation of \(\dot{ {G} { {M}}_ \odot} {\mathrm{/}}{ {G}}{ {M}}_ \odot\) that combines the temporal variations of both G and \({ {M}}_ \odot\). The retrieved negative rate is close to the theoretical computations of the Sunâs mass loss due to interior processes and solar wind. The fusion cycle that generates energy into the Sun relies on the conversion of hydrogen into helium, which is responsible for a solar mass reduction with a rate of ~â0.679âÃâ10â13 per year28. On the other hand, the solar wind contribution is more uncertain. The solar cycle significantly influences the solar mass loss rate due to solar wind. Estimates of the mass carried away with the solar wind showed rates between â(2â3)âÃâ10â14\({ {M}}_ \odot\) per year28, whereas numerical simulations of coupled corona and solar wind models provided rates between â(4.2â6.9)âÃâ10â14\({ {M}}_ \odot\) per year29. Therefore, a mean value of the total solar mass loss of â(0.9â1.1)âÃâ10â13\({ {M}}_ \odot\) per year would be expected, since the MESSENGER mission operated during ~2/3 of an entire solar cycle whose maximum occurred in proximity of the end of the mission.
Discussion
The estimated \(\dot{ {G} { {M}}_ \odot} /{ {G}}{ {M}}_ \odot\) represents one of the first experimental observations of the solar mass loss. Previous studies have demonstrated the feasibility of estimating this parameter by adjusting the planetary ephemerides20,25,28. Their results, which are consistent with our estimates, were limited by the data availability and possible mismodeling of spacecraft orbits. Our processing of the entire MESSENGER mission increased the solution sensitivity to a variation of \(\dot{ {G} { {M}}_ \odot} /{ {G}}{ {M}}_ \odot\), which has a quadratic dependence in time. Furthermore, our new technique, which consists of a double-integration and a combined estimation of both planet and spacecraft orbits, mitigates the systematic errors related to the spacecraft position and velocity.
The discrepancy between our solution and the computed \(\dot { {M}}_ \odot /{ {M}}_ \odot\) may be interpreted as an indirect measurement of the universal constant time variation. The reconstructed \(\dot G/G\) is lower than 4.0âÃâ10â14 per year with an uncertainty that is mainly limited by the knowledge of the solar interior evolution (2Ïâ=â5.0âÃâ10â14 per year). This result strengthens the hypothesis that \(\dot{ {G}}/{ {G}}\) is close to 0, improving the estimates of LLR studies by almost an order of magnitude.
To validate the accuracy of these results, we reintegrated the orbit of Mercury with our adjusted values. Figure 3 shows the required corrections to the MESSENGER range data to fit at the noise level shown in Fig. 1. The red, blue, and green dots are the measurement biases needed with the JPL DE430 and DE432 ephemerides, and our solution, respectively. This plot shows major improvements compared to previous JPL ephemerides. The DE430 Mercury trajectory is affected by 80-m amplitude errors that were corrected in the later DE432, with remaining 5â10âm errors. Our reintegrated ephemeris for Mercury, which is available on the NASA Goddard Space Flight Center (GSFC) Planetary Geodynamics Data Archive30, shows only 0.5â3âm biases over the full mission.
Temporal distribution of the range biases with three Mercuryâs ephemeris. The measurement biases are required to fit the MESSENGER range data at the noise level with the JPL DE430 (purple) and D432 (blue) ephemerides, and our integrated trajectory for Mercury (black). These biases were used to determine the quality of the ephemeris results. After convergence of the global solution, all the adjusted parameters (Methods) are applied in a final iteration, in which the range biases are adjusted instead of the Mercuryâs initial state, \({{G}}{{M}}_ \odot\), \({{J}}_{2 \odot }\), β, η, and \(\dot{{{ {G} { {M}}}}_ \odot} /{{G}}{{M}}_ \odot\). Large range biases suggest significant errors in the planetâs ephemeris
The stability of Mercuryâs orbit integration also depends on the ephemerides of the other bodies of the solar system that are provided by the JPL ephemerides. Therefore, we evaluated the changes in recovery of the heliophysics and relativity parameters when using DE430 or DE436 for planetary ephemerides and initial state of Mercury (instead of DE432 previously) and modeling the ephemerides of the asteroids with the JPL AST343DE43017. These solutions, which are reported in the Supplementary Table 3, show consistent results with the estimated values in Table 2. The formal uncertainties do not take into account probable errors in the planetsâ GMs and trajectories, as they rely only on the measurement accuracies and the correlation among the adjusted parameters. The fourth column of Table 2 shows the maximum estimation differences between the three cases, which rely on the JPL DE430, DE432, or DE436 for the other planet trajectories. These values may be interpreted as conservative uncertainty bounds that account for pessimistic errors due to mismodeling of GMs and orbits of the solar system bodies. The discrepancies between \(\dot{ {G} { {M}}_ \odot} {\mathrm{/}}{ {G}}{ {M}}_ \odot\) solutions are slightly larger than the formal uncertainties because of variations in Earthâs orbit, Jupiterâs GM and orbit, and SSB location between the different JPL DEs that our methodology cannot mitigate. Furthermore, the value of the solar gravitational constant is significantly affected by these discrepancies leading to estimates that are, however, still within 2- and 3-Ïs (Supplementary Table 3).
In conclusion, our analysis of Mercuryâs ephemeris with the MESSENGER data enhances the knowledge of the relativistic parameter η, confirming predictions of Einsteinâs theory of GR. We provide one of the first observations of the solar system expansion due to the solar mass loss. The negative rate of \(\dot{ {G} { {M}}_ \odot} /{ {G}}{ {M}}_ \odot\) is very close to theoretical computations of the Sunâs mass loss rate leading us to significantly constrain the universal constant time variation. These results are mainly limited by the uncertainty in planet and asteroid ephemerides that perturb Mercuryâs orbit. We demonstrate the potential of measuring the planetsâ relative distances over decadal timescales to provide a better understanding of the solar system and Sun evolution. To pursue these challenging scientific goals, future investigations employing precision ranging from a dedicated multi-spacecraft constellation at interplanetary scale may provide a leap in planetary science, heliophysics, and theoretical physics31.
Methods
Parametrized post-Newtonian formulation
The orbital dynamics of planets, satellites, and asteroids relies mainly on the gravitational attraction of the other bodies that are modeled as external point masses. Nevertheless, Newtonâs law of universal gravitation needs to be modified to include Einsteinâs relativistic corrections by means of the PPN n-body formalism17. The acceleration of Mercury due to the interaction with other point masses is, therefore, given by:
where r i , v i , a i , and μ i are the position, velocity, and acceleration vectors with respect to the SSB, and the gravitational parameter of the body i, respectively, c is the speed of light, and β and γ are the PPN parameters that measure the nonlinearity in superposition of gravity and space-time curvature produced by unit rest mass, respectively. This formulation has been applied to the orbital dynamics of both Mercury and MESSENGER. The bodies included in these integrations are all major bodies, including the Sun, Moon, planets, Pluto, and 343 asteroids in the main belt between Mars and Jupiter. The positions, velocities, accelerations, and gravitational parameters are obtained from the JPL DE432 ephemerides18 for planets, and the JPL AST343DE430 for asteroids17. However, to test the stability of the solution, the global estimation was repeated by using JPL DE430 and DE436 ephemerides as a priori. The results of Mercuryâs ephemeris, GMâ, \({ {J}}_{{\it{2}} \odot }\), β, η, and \(\dot{ {G} { {M}}}_ \odot /{ {G}}{ {M}}_ \odot\) are within the corrected Ï reported in Table 2.
LenseâThirring precession
The dynamical orbital equations of both Mercury and MESSENGER account for the LenseâThirring effect due to the Sunâs gravitomagnetic field that leads to a secular precession of the heliocentric longitude of the ascending node and argument of pericenter.
This precession is a prediction of GR, and, for this reason, it has been recently renamed as ELT effect. Einstein postulated the frame-dragging in the context of the general theory of relativity stating that non-static stationary distributions of mass-energy affect space-time. In 1918, Josef Lense and Hans Thirring derived the first frame-dragging effect predicting that the rotation of a massive body induces a distortion of space-time. The ELT effect has been measured with the LAGEOS satellites in orbit around the Earth32, and the gyroscopes of the Gravity Probe B33. A test on ELT effects was initially proposed with the NASA mission Juno in orbit about Jupiter34. The large angular momentum of the planet induces a significant ELT acceleration; however, it is also highly correlated with perturbations due to Jupiterâs orientation, which is currently not sufficiently defined35.
The ELT effect on Mercury due to the solar gravitomagnetic field is not negligible, and it may theoretically be used to measure the angular momentum of the Sun. Mercuryâs acceleration due to the ELT effect is:
where \({ {S}}_ \odot\) is the angular momentum of the Sun, G is the universal gravitational constant, \({\mathbf{r}}_{{\mathrm{M}} \odot }\) and \({\mathbf{v}}_{{\mathrm{M}} \odot }\) are the heliocentric position and velocity vectors of Mercury, respectively, and \(\widehat {p_ \odot }\) is the unit vector of the Sunâs pole direction, which relies on the right ascension \(\alpha _ \odot\)â=â286.13° and declination \(\delta _ \odot\)â=â63.87° of the pole defined in the International Celestial Reference Frame (ICRF)36.
The ELT effect on Mercuryâs orbit is mainly in the radial direction with a maximum acceleration of ~2âÃâ10â7âmâsâ2, assuming \({ {S}}_ \odot\)â=â190âÃâ10â39âkgâm2âsâ1. However, the perturbation induced by the ELT precession is strongly anti-correlated with the effect due to \({ {J}}_{2 \odot }\), and the recovery of \({ {S}}_ \odot\) is unachievable with the estimation of Mercuryâs ephemeris.
Strong equivalence principle
Milani et al.37 formulated, for the first time, a redefinition of the SSB due to violations of the SEP. This effect causes a significant indirect perturbation on Mercuryâs orbit that enables an accurate measurement of η by adjusting the planetâs ephemeris. These results provoked a scientific debate on the consequences of SEP violations for the modeling of planetary ephemerides. Ashby et al.38 presented an alternative approach that does not fully include the indirect perturbation presented by Milani et al.37, limiting the contribution of η on planetâs orbital dynamics. However, current planetary ephemerides studies are based on the hypothesis that the gravitational and inertial masses are equal to compute the SSB location17.
The SSB represents the origin of the ephemerides reference frame. The assumptions to compute its position are the conservation of mass/energy and the momentum of the solar system. The SSB is, then, approximated as follows13:
where \({\cal R}\) is equal to 0 if the SSB is the origin of the reference frame, r j is the relative distance of body j with respect to the SSB, and:
where \(GM_{\mathrm{j}}^{\mathrm{G}}\) is the gravitational mass parameter of body j, r jk is \(\left| {r_{k} - r_{j}} \right|\) and v j is the magnitude of the velocity of body j. This formulation is valid only when the SEP is not violated (ηâ=â0). The inertial masses should be used in the computation of the SSB, as follows, if we neglect terms of order 1/c2:
where r j is the position vector of body j with respect to the SSB, and \(M_{j}^{\mathrm{I}}\) is the inertial mass of body j. However, these inertial masses are unknown since the masses of the Sun, planets, and satellites are retrieved in space by means of their gravitational pull. A violation of the SEP may lead to an intrinsic mismodeling of the SSB position. To account for this effect, the position of the Sun should be redefined by:
where the sum includes planets and asteroids, \({{\Omega }}_{j}\) and M j are the self-gravitational energy and the mass of body j, respectively, and the symbol â stands for the Sun. The self-gravitational energy of the Sun, Earth, and Moon are â3.52âÃâ10â6, â4.64âÃâ10â10, and â1.88âÃâ10â11, respectively39,40. The self-gravitational energy of the other planets is computed by assuming uniform density \(\left( {\Omega _{j} = \frac{3}{5}\frac{{GM_{j}^2}}{R}} \right)\). We also tested other self-gravitational energy modeling for the other planets, but the η estimates only changed within 1-Ï since the Sunâs self-gravitational energy represents the dominant term.
The Sunâs position correction (Eq. 9) entails an indirect term in the heliocentric acceleration of Mercury. The partial derivative of Mercuryâs heliocentric acceleration (aM) with respect to η, which enables the estimation of this parameter by adjusting the planet ephemeris, is:
where the symbol M stands for Mercury, r kj is the position vector of the relative distance between bodies k and j, \({\mathbf{r}}_ \odot\) is the position vector of the Sun with respect to the SSB, and the last term is the indirect effect due to the correction of the SSB position neglecting terms of the order 1/η2. SEP violations would provide significant perturbations on Mercuryâs orbit that enable the measurement of η to high accuracy, and decorrelate the PPN parameters and \({ {J}}_{{\it{2}} \odot }\) if the Nordtvedt equation (Eq. 1) is applied as a priori constraint13. However, the correlation between \({ {J}}_{{\it{2}} \odot }\) and β is still ~0.9 even applying the Nordtvedt equation (Supplementary Table 1).
This a priori constraint approach was proposed for the first time in the simulations of the relativity experiment that will be conducted by the European Space Agency (ESA) mission BepiColombo13. One year of operations in orbit about Mercury will allow BepiColombo to collect 30-cm precision range data for the determination of Mercuryâs ephemeris. The results of those simulations showed lower correlation between \({ {J}}_{{\it{2}} \odot }\) and β (~â0.3) by using the Nordtvedt equation37. The stronger effect of this constraint on the BepiColombo solutions is mainly due to the more precise range data that will enable to determine a more accurate estimate of SEP violations. The accuracy of η estimation affects directly the correlation between \({{J}}_{{\it{2}} \odot }\) and β, if the constraint is applied. If we assume the Nordtvedt equation and to know η at the same level of BepiColombo results (~10â6), the correlation between \({ {J}}_{{\it{2}} \odot }\) and β drops to ~0.3 that is consistent with the simulation of the future ESA mission to Mercury41.
Time-variable gravitational constant
The time-varying gravitational parameter \(\dot{ {G} { {M}}_ \odot} /{ {G}}{ {M}}_ \odot\) is defined by the sum of the time-variations of the gravitational universal constant \(\dot G/G\) and the mass of the Sun \(\dot { {M}}_ \odot /{ {M}}_ \odot\). The additional term of Mercuryâs heliocentric acceleration due to \(\dot{ {G} { {M}}_ \odot} /{ {G}}{ {M}}_ \odot\) is:
where Ît is the difference between the current epoch and the reference epoch J2000 (1 January 2000 at 1200 UTC), and \({\mathbf{r}}_{{{M}} \odot }\) is the relative position vector between Mercury and the Sun.
Ephemeris and orbit determination
The results presented in this paper were obtained with the NASA GSFC orbit determination software GEODYN II, which has been used to determine geophysical parameters of, for example, the Earth, Moon, and Mars. We used GEODYN II to recover previous solutions of Mercuryâs gravity field, orientation, and tides assuming the JPL DE430 ephemeris of Mercury17. To estimate Mercuryâs ephemeris and the associated heliophysics and relativity parameters, we modified GEODYN II to numerically integrate the orbits of both MESSENGER and the planet Mercury simultaneously.
This software is based on a batch least-squares scheme that allows the combination of all observations within one batch (arc) for the estimation of the parameters of interest. The least-squares technique relies on an adjustment of model parameters to minimize the discrepancies between the computed observables and actual measurements (residuals). If the trajectory of the spacecraft alone is integrated, the only parameters that can be estimated are related to MESSENGERâs dynamics around Mercury (e.g., the gravity field of the planet). The simultaneous numerical integration of the planet ephemeris allows the adjustment of other model parameters, such as those from heliophysics and relativity that perturb the orbit of Mercury.
The MESSENGER orbital mission (2011â2015) was partitioned in 1499 1-day arcs. Three additional ~10-day arcs cover the three Mercury flybys in 2008â2009. The range data were weighted according to the contribution of the solar plasma that varied through the mission as expressed by the SunâEarthâProbe angle. For each arc, the Mercuryâs ephemeris is continuously integrated from the Flyby 1 initial epoch (7 January 2008 at 0000 UTC). We generated partial derivatives of the following MESSENGER-related parameters: spacecraft initial states, areas of the spacecraft sunshield, and solar panels, Mercuryâs gravity field up to degree and order 100 in spherical harmonics, and Mercuryâs Love number k2 and orientation (poleâs right ascension and declination). We also computed partial derivatives of the following Mercury-related parameters: planetâs initial state, \({ {G}}{ {M}}_ \odot\), \(\dot{ {G} { {M}}}_ \odot /{ {G}}{ {M}}_ \odot\), \({ {J}}_{{\it{2}}{\it{ \odot }}}\) and \({ {S}}_{\it{ \odot }}\), PPN parameters β and γ, and Nordtvedtâs parameter η.
The individual normal equations of all these arcs were combined and inverted to yield the final estimates of the geophysical, heliophysics, and relativity parameters. The orbit of the Earth is not integrated and adjusted in this study since the orbital accuracy of the Earth from the JPL DE432 ephemerides is comparable to the precision of the MESSENGER data.
MESSENGER data set
The data processed in this paper include the three MESSENGER flybys around Mercury, and the whole orbital mission. The three flybys occurred on 14 January 2008, 6 October 2008, and 29 September 2009, respectively. MESSENGER was inserted in a highly eccentric and near-polar orbit about Mercury on 18 March 2011. The initial period was ~12âh and the orbital periapsis was at ~200-km altitude and ~60°N latitude. Orbit-correction maneuvers (OCMs) were required to maintain the periapsis between 200 and 500âkm for the first year of operations. The third-body perturbation of the Sun combined with the high eccentricity of the orbit led to a significant drift of the periapsis altitude and latitude.
The mission was extended for a second year in March 2012. The OCMs became less frequent, and one of them was used to reduce the orbital period to ~8âh. A second extended mission (XM2) started in March 2013 and included a low-altitude campaign until Mercury impact on 28 April 2015. The fuel reserves enabled the spacecraft to maintain periapsis altitudes as low as 15â25âkm for several weeks. NASAâs DSN stations tracked the spacecraft during part of these passages from April to October 2014 leading to accurate measurements of Mercuryâs gravity at altitudes between 25 and 100âkm. In the last 6 months of the mission, the closest approaches of MESSENGER were occulted by Mercury and were thus not visible from the Earth. However, additional range-rate and range measurements were collected at low altitudes between 75 and 100âkm.
The data included in this study were collected over ~900 days. The greater part of the excluded data is because of high levels of plasma noise in proximity of superior solar conjunctions (SPEâ<â35°). Other arcs were also omitted because of the presence of OCMs or reaction wheel momentum desaturation maneuvers that imparted significant ÎVs leading to significant orbital errors.
Range-rate and range measurements
The analysis of the range data to estimate Mercuryâs ephemeris relies strongly on the accuracy of MESSENGER orbital reconstruction. The data collected during XM2, especially, are very sensitive to the quality of the spacecraft orbits. Uncompensated gravity anomalies of Mercuryâs gravity field may affect significantly the range residuals leading to inaccurate ephemeris solutions. To mitigate the effects of MESSENGER orbital errors in the determination of Mercuryâs ephemeris, both range-rate and range data have been analyzed in this study. This data set includes two-way and three-way coherent range-rate and two-way coherent range measurements. The difference between two- and three-way data is only related to the receiving station. The signal is transmitted by the DSN station and sent coherently back to the same (two-way) or a different (three-way) station by the spacecraft deep space transponder (DST). The two-way configuration guarantees highly accurate data thanks to the H-masers at the DSN ground stations. The three-way data require additional bias corrections due to the time delay between the oscillators at the transmitting and receiving stations. The biases of the three-way range-rate data are adjusted in the solution to mitigate this error source.
The Earth-spacecraft radio link was supported by diametrically opposite-facing PAAs for the high-gain downlink signal, and two fanbeam antennas to provide medium-gain uplink and downlink. Four low-gain antennas were also used to enable the instrumentation pointing towards the planet surface during tracking periods. However, the range data campaigns were always conducted with the front and back PPAs, as shown in Fig. 1. The gain level of the antennas influences significantly the level of noise of the range-rate data21. A major source of error for the range data is the internal spacecraft delay that was measured during ground testing with an uncertainty of ~12â14âns that leads to a range accuracy of <2âm21. Further tests in flight enabled to reconstruct a more precise delay time, which was necessary for science operations at Mercury. MESSENGER operated for ~11 years in space, and its instrumentations, including the transponder, coped with the effects of ageing. By interpolating the range data residuals, we were able to determine a linear trend of the time delay that is probably associated with the ageing of the spacecraft transponder. The rate of the mean time delay is ~0.45âns (~13.5âcm) per year, which provides a maximum offset of <1âm between January 2008 and April 2015. This effect is within the level of accuracy of the range data (1â2âm) that was retrieved during test laboratory results21.
Another source of range data error is given by station biases due to imperfect calibration. The accurate measurement of the ranging signal round-trip delay is made through digital signal processors at the DSN stations by correlating the uplink and downlink carriers that are coherently related. This calibration may lead to biases on the measured delay with a standard deviation of 1â3âm. To mitigate these calibration errors, range station biases for each tracking pass may be estimated in orbit determination. However, the estimation of the range station biases tends to absorb the uncompensated ephemeris mismodeling. For this reason, the range station biases are not estimated in this study of Mercuryâs ephemeris, heliophysics and GR.
These biases may instead be used to determine the quality of the ephemeris results. After convergence of the global solution, all the adjusted parameters (see Ephemeris and Orbit Determination) are applied in a final iteration, in which the range station biases are adjusted instead of the Mercuryâs initial state, GMâ, \({{J}}_{{\it{2}} \odot }\), β, η, and \(\dot{ {G} { {M}}_ \odot} /{ {G}}{ {M}}_ \odot\). Figure 3 shows the retrieved station biases that are within the expected range of calibration errors. To compare the quality of these results, Fig. 3 shows the range biases estimated by using JPL DE430 and DE432 original settings.
Data availability
The MESSENGER radio tracking data are available from the NASA Planetary Data System archive (http://pds-geosciences.wustl.edu/missions/messenger/rs.htm). The retrieved ephemeris of Mercury is available on the NASA GSFC Planetary Geodynamics Data Archive30.
References
Einstein, A. Relativity: The Special and General Theory (Princeton University Press, Princeton, 2015).
Clemence, G. M. The relativity effect in planetary motions. Rev. Mod. Phys. 19, 361â364 (1947).
Einstein, A. Erklarung der Perihelionbewegung der Merkur aus der allgemeinen Relativitatstheorie. Sitzubgsber. Preuss. Akad. Wiss. 47, 831â839 (1915).
Pireaux, S. & Rozelot, J.-P. Solar quadrupole moment and purely relativistic gravitation contributions to Mercuryâs perihelion advance. Astrophys. Space Sci. 284, 1159â1194 (2003).
Iorio, L. Constraining the angular momentum of the Sun with planetary orbital motions and general relativity. Sol. Phys. 281, 815â826 (2012).
Nordtvedt, K. Jr. Post-Newtonian metric for a general class of scalar-tensor gravitational theories and observational consequences. Astrophys. J. 161, 1059â1067 (1970).
Wagner, T. A., Schlamminger, S., Gundlach, J. H. & Adelberger, E. G. Torsion-balance tests of the weak equivalence principle. Class. Quant. Grav. 29, 184002 (2012).
Nordtvedt, K. Jr. Equivalence principle for massive bodies. II. Theory. Phys. Rev. 169, 1017â1025 (1968).
Williams, J. G. et al. New test of the equivalence principle from lunar laser ranging. Phys. Rev. Lett. 36, 551â554 (1976).
Williams, J. G., Turyshev, S. G. & Boggs, D. H. Progress in lunar laser ranging tests of relativistic gravity. Phys. Rev. Lett. 93, 261101 (2004).
Williams, J. G., Turyshev, S. G. & Boggs, D. H. Lunar laser ranging tests of the equivalence principle. Class. Quant. Grav. 29, 184004 (2012).
Müller, J., Hofmann, F., Fang, X. & Biskupek, L. Lunar laser ranging: recent results based on refined modelling. Earth Edge Sci. Sustain. Planet 139, 447â451 (2014).
Milani, A., Tommei, G., Vokrouhlický, D., Latorre, E. & Cicalò, S. Relativistic models for the BepiColombo radioscience experiment. Proc. Int. Astron. Union 5, 356â365 (2009).
Noerdlinger, P. D. Solar mass loss, the astronomical unit, and the scale of the solar system. Preprint at https://arxiv.org/pdf/0801.3807.pdf (2008).
Solomon, S. C. et al. The MESSENGER mission to Mercury: scientific objectives and implementation. Planet. Space Sci. 49, 1445â1465 (2001).
Mazarico, E. et al. The gravity field, orientation, and ephemeris of Mercury from MESSENGER observations after three years in orbit. J. Geophys. Res. Planets 119, 2417â2436 (2014).
Folkner, W. M., Williams, J. G., Boggs, D. H., Park, R. S., and Kuchynka, P., The planetary and lunar ephemerides DE430 and DE431. Interplanet. Netw. Prog. Rep., 196, 1-81, Jet Propulsion Laboratory, Pasadena, CA, USA. https://naif.jpl.nasa.gov/pub/naif/generic_kernels/spk/planets/de430_and_de431.pdf (2014).
Folkner, W. M. Planetary ephemeris DE432. Technical Report. Jet Propulsion Laboratory, Pasadena, CA, USA. https://naif.jpl.nasa.gov/pub/naif/generic_kernels/spk/planets/de432.pdf (2014).
Pitjeva, E. V. & Pitjev, N. P. Development of planetary ephemerides EPM and their applications. Celest. Mech. Dyn. Astron. 119, 237â256 (2014).
Fienga, A., Laskar, J., Exertier, P., Manche, H. & Gastineau, M. Numerical estimation of the sensitivity of INPOP planetary ephemerides to general relativity parameters. Celest. Mech. Dyn. Astron. 123, 325â349 (2015).
Srinivasan, D. K., Perry, M. E., Fielhauer, K. B., Smith, D. E. & Zuber, M. T. The radio frequency subsystem and radio science on the MESSENGER mission. Space Sci. Rev. 131, 557â571 (2007).
Mecheri, R., Abdelatif, T., Irbah, A., Provost, J. & Berthomieu, G. New values of gravitational moments J 2 and J 4 deduced from helioseismology. Sol. Phys. 222, 191â197 (2004).
Pijpers, F. P. Helioseismic determination of the solar gravitational quadrupole moment. Mon. Not. R. Astron. Soc. 297, L76âL80 (1998).
Park, R. S. et al. Precession of Mercuryâs Perihelion from Ranging to the MESSENGER Spacecraft. Astron. J. 153, 121 (2017).
Pitjeva, E. V. & Pitjev, N. P. Relativistic effects and dark matter in the Solar system from observations of planets and spacecraft. Mon. Not. R. Astron. Soc. 432, 3431â3437 (2013).
Zhu, W. W. et al. Testing theories of gravitation using 21-year timing of pulsar binary J1713+ 0747. Astrophys. J. 809, 41 (2015).
Bertotti, B., Iess, L. & Tortora, P. A test of general relativity using radio links with the Cassini spacecraft. Nature 425, 374â376 (2003).
Pitjeva, E. V. & Pitjev, N. P. Changes in the Sunâs mass and gravitational constant estimated using modern observations of planets and spacecraft. Sol. Syst. Res. 46, 78â87 (2012).
Pinto, R. F., Brun, A. S., Jouve, L. & Grappin, R. Coupling the solar dynamo and the corona: wind properties, mass, and momentum losses during an activity cycle. Astrophys. J. 737, 72 (2011).
NASAâs Planetary Geodynamics Laboratory, Goddard Space Flight Center, Greenbelt, MD, USA. https://pgda.gsfc.nasa.gov (2017).
Zuber, M.T. et al. From Copernicus to Newton to Einstein: toward a dynamical understanding of the Solar System. Planetary Science Vision 2050 Workshop. Preprint at https://www.hou.usra.edu/meetings/V2050/pdf/8074.pdf (2017).
Ciufolini, I. & Pavlis, E. C. A confirmation of the general relativistic prediction of the Lense-Thirring effect. Nature 431, 958â960 (2004).
Everitt, C. W. F. et al. Gravity probe B: final results of a space experiment to test general relativity. Phys. Rev. Lett. 106, 221101 (2011).
Iorio, L. A possible new test of general relativity with Juno. Class. Quant. Grav. 30, 195011 (2013).
Le Maistre, S., Folkner, W. M., Jacobson, R. A. & Serra, D. Jupiter spin-pole precession rate and moment of inertia from Juno radio-science observations. Planet. Space Sci. 126, 78â92 (2016).
Archinal, B. A. et al. Report of the IAU working group on cartographic coordinates and rotational elements: 2009. Celest. Mech. Dyn. Astron. 109, 101â135 (2011).
Milani, A., Vokrouhlický, D., Villani, D., Bonanno, C. & Rossi, A. Testing general relativity with the BepiColombo radio science experiment. Phys. Rev. D 66, 082001 (2002).
Ashby, N., Bender, P. L. & Wahr, J. M. Future gravitational physics tests from ranging to the BepiColombo Mercury planetary orbiter. Phys. Rev. D 75, 022001 (2007).
Anderson, J. D., Gross, M., Nordtvedt, K. L. & Turyshev, S. G. The solar test of the equivalence principle. Astrophys. J. 459, 365â370 (1996).
Congedo, G. & De Marchi, F. Testing the strong equivalence principle with spacecraft ranging towards the nearby Lagrangian points. Phys. Rev. D 93, 102003 (2016).
De Marchi, F., Tommei, G., Milani, A. & Schettino, G. Constraining the Nordtvedt parameter with the BepiColombo Radioscience experiment. Phys. Rev. D 93, 123014 (2016).
Acknowledgements
This work was conducted at the NASA Goddard Space Flight Center. We are grateful to J.B. Nicholas (EST, Inc.), D.E. Pavlis (SGT, Inc.), and D.D. Rowlands (NASA, GSFC) for their help with the GEODYN II software. A.G. thanks L. Iorio (MIUR), and A. Milani, (University of Pisa) for ideas and discussions. The data used in this paper are available at http://pds-geosciences.wustl.edu/missions/messenger/rs.htm.
Author information
Authors and Affiliations
Contributions
A.G., E.M. and S.G. performed radio tracking data processing and preliminary analysis of the MESSENGER orbits. A.G. developed updated relativistic and solar modeling in the NASA GSFC orbit determination software (GEODYN II). A.G., E.M., S.G., F.G.L., G.A.N., D.E.S. and M.T.Z. contributed to the interpretation of the results. A.G. wrote the manuscript with input from E.M., S.G., F.G.L., G.A.N., D.E.S. and M.T.Z.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Genova, A., Mazarico, E., Goossens, S. et al. Solar system expansion and strong equivalence principle as seen by the NASA MESSENGER mission. Nat Commun 9, 289 (2018). https://doi.org/10.1038/s41467-017-02558-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-017-02558-1