# A Priori “Imprinting” of General Relativity Itself on Some Tests of It?

Authors: Lorenzo Iorio
Journal: Advances in Astronomy
 (2010)

---

## Abstract

We investigate the effect of possible a priori “imprinting” effects of general
relativity itself on satellite/spacecraft-based tests of it. We deal with some performed or proposed time-delay ranging experiments in the sun's gravitational field. It turns out that the “imprint” of general relativity on the Astronomical Unit and the solar gravitational constantGM⊙, not solved for in the so far performed spacecraft-based time-delay tests, induces an a priori bias of the order of
10-6 in typical solar system ranging experiments aimed to measure the space curvature PPN parameter γ. It is too small by one order of magnitude to be of concern for the performed Cassini experiment, but it would affect future planned or proposed tests aiming to reach a 10-7–10-9 accuracy in determining γ.


---

## Body

## 1. Introduction

Several space-based tests of the Einstein’s General Theory of Relativity (GTR) have been performed, or attempted, in the more or less recent past by following an “opportunistic” approach, that is, by suitably analyzing existing data sets of artificial satellites or interplanetary spacecrafts almost always built and launched for different original purposes (e.g., satellite geodesy, geodynamics, planetology, etc.). A cornerstone result was the Cassini radio science experiment [1] which led to constraining the deviation of the PPN parameter γ from its general relativistic value of unity at a 10-5 level [2, 3] through the measurement of the general relativistic time delay affecting the electromagnetic waves linking the Earth and the Cassini spacecraft during its journey to Saturn when they were in superior conjunction, that is, aligned and on opposite sides of the Sun. Many other high-precision space-based tests, aiming to reach accuracy levels as high as 10-7–10-9 in determining γ, have been proposed [4–8]. Indeed, some theoretical scenarios predict deviations from GTR at such a level [9–12]; for a recent overview see, for example, [13].In this paper we wish to critically discuss certain subtle issues pertaining the consistency of such data analyses interpreted as genuine tests of GTR. In Section2 we outline the problem in general terms. We will move to concrete examples in Section 3 in which we will deal with the Cassini-like ranging experiments (Section 3.1). Section 4 is devoted to the conclusions.


## 2. The Issue of the a Priori “Imprinting” of GTR Itself in Tests Dedicated to It: General Considerations

In general, in such tests huge data sets from man-made interplanetary probes are confronted with the corresponding predictions computed from given dynamical models of their motion in which the relativistic effect to be tested was explicitly included, with one or more solve-for parameters{P} accounting for it to be estimated in a least-square fashion along with many other ones {K}, not directly pertaining GTR.Of crucial importance for interpreting such data analyses as genuine tests of GTR is to clarify how the numerical values of the models' parameters{F} which have been kept fixed to certain reference values, that is, which have not been solved-for, have been originally obtained. From the point of view of testing GTR, it is not enough that the resulting postfit residuals of a certain directly observable quantity are statistically compatible with zero at a good level; the standard data reduction procedure used for the original goals of the missions exploited for the GTR test considered may not be valid, in principle, for performing a truly unbiased, genuine check of GTR which is not a “tautology”.Indeed, if the primary task of a space-based mission is, for example, to reach a certain astronomical target with a given accuracy, the only thing that is important to this aim is that the dynamical models adopted to predict the probe's motion are accurate enough; this is usually quantitatively judged by inspecting the postfit residuals of some directly measurable quantities like, for example, ranges or range-rates. How the parameters{F} entering the models have been obtained, that is, their a priori values, does not matter at all; the only important thing is that the resulting fit of an existing set of observations is good enough to minimize the observable residuals.Such an approach may, in principle, not be entirely adequate when the goal of the data analysis is testing a gravitational theory like GTR in an unambiguous, unbiased, and self-consistent way. In this case, the way the fixed parameters{F} of the models have been obtained does, in fact, matter. Indeed, if one or more of them {I} were previously obtained from different data of different bodies in such a way that they somehow retain a nonnegligible a priori “imprint” of the same effect we are now interested in, their use may bias the current test just towards the desired outcome yielding, for example, a very high accuracy corroboration. In this case, it would be more correct to use, if possible, values of such “imprinted” parameters {I} which have been obtained independently of the effect itself whose existence we are just testing in the present data analysis, even if the accuracy of such different values of the “suspect” parameters {I} was worse. Alternatively, if, for some reasons, such “unbiased” values are not available, {I} should be included, if possible, in the list of the solved-for parameters along with the one(s) {P} accounting for the effect to be tested, and the resulting covariance matrix should be checked to inspect the correlations among them. The price to be paid may be an overall accuracy of the test not so high as that previously obtained, but we would have more epistemologically consistent, reliable, and trustable tests (inprinciple, certain positive corroborations may have theirstatus changed).


## 3. Application to Some Concrete Cases

### 3.1. The Cassini Radio Science Test and Other Proposed High-Accuracy Space-Based Measurements ofγ

To be more definite, let us look at the Cassini radio science test. In that case, the radiotechnical data of the spacecraft traveling to Saturn were contrasted with a set of dynamical models by JPL of its motion and electromagnetic waves propagation in such a way that a correctionΔγ to the GTR-predicted value of the PPN parameter γ was solved for, among other parameters, obtaining [2]
(1)Δγ≡|γ-1|=(2.1±2.3)×10-5;
other authors got [3]
(2)Δγ≡|γ-1|=(-1.3±5.2)×10-5.Now, a physical parameter which is, of course, crucial in such a test is the gravitational parameterGM⊙ of the Sun, which is the source of the relativistic time delay put on the test. It was not estimated [2, 3], so that its numerical value was kept fixed to the standard reference figure of the JPL DE ephemerides used. It does, in principle, contain an a priori “imprinting” by GTR itself through the same effect itself that was just tested with Cassini, in particular by γ itself. Indeed, the numerical value of GM⊙ comes from the fixed value of the defining Gaussian constant (see on the WEB http://ssd.jpl.nasa.gov/?constants.) (3)k=0.01720209895au3/2d-1,
and from the value of the Astronomical Unit, not estimated in the Cassini tests,
(4)AU=1.49597870691×1011(±3)m
through
(5)GM⊙=k2AU3d-2=1.32712440018×1020(±8×109)m3s-2.(Here we will use au for the symbol of the Astronomical Unit, like m for the meter, while AU will denote its numerical value in m.)AU was, in fact, obtained just through a combination of radar ranging of Mercury, Venus, and Mars, laser ranging of the moon (making use of light reflectors left on the lunar surface by Apollo astronauts), and timing of signals returned from spacecrafts as they orbit or make close passes and objects in the solar system [14]; thus, it is affected in a nonnegligible way, given the level of accuracy of the techniques adopted, by GTR itself and, in particular, by γ which enters the PPN expressions for the time delay and the bending of traveling electromagnetic waves. Thus, there exists, in principle, the possibility that the high-accuracy results of the Cassini radio science tests may retain an a priori “imprint” of GTR itself through GM⊙ (and the Astronomical Unit as well).Let us put our hypothesis on the test by making some concrete calculations; for the sake of clarity, we will refer to the Cassini radio science tests, but the conclusions may be considered valid also for any of the many proposedγ-dedicated missions.The GTR time delay experienced by electromagnetic waves propagating from point 1 to point 2 is(6)Δt=2Rgcln(r1+r2+r12+Rgr1+r2-r12+Rg),
where Rg=2GM⊙/c2 is the Sun's Schwarzschild radius in which G is the Newtonian gravitational constant, M⊙ is the solar mass, c is the speed of light in vacuum, r1 is the heliocentric coordinate distance to point 1, r2 is the heliocentric coordinate distance to point 2, and r12 is the distance between the points 1 and 2. Equation (6) is the expression actually used in the JPL's Orbit Determination Program (ODP) to analyze interplanetary ranging to planets and probes. In order to quantitatively evaluate the level of “imprinting” by GTR itself in the used value of the Astronomical Unit, let us assume that r1 equal to the Earth-Sun distance and let us vary r2 within 0.38 au and 1.5 au to account for the ranging to inner planets; the maximum effect occurs at the superior conjunction, that is, when n1≈-n2 and r12≈r1+r2.(Here n≡r/r.) It turns out that Δtranging≈4×10-4 s, which is certainly not negligible with respect to the accuracy of the order of 10-8 s with which the light time for 1 au τA is actually measured (http://ssd.jpl.nasa.gov/?constants).(The value in km of the Astronomical Unit is obtained by measuring at a given epoch the distance between the Earth and a target body (a planet or a probe orbiting it) by multiplying c times the round trip travel time τ of electromagnetic waves sent from the Earth and reflected back by the target body, and confronting it with the distance, expressed in AU, between the Earth and the target body at the same epoch as predicted by some accurate dynamical ephemeris [14].) As a consequence, the quantitative impact of the interplanetary ranging in the inner solar system on the determination of the Astronomical Unit is of the order
(7)dAU=cΔtranging=1.14291×105m,
not negligible with respect to the meter-level accuracy in measuring the Astronomical Unit; thus, dAU/AU=8×10-7. Differentiating (5) with respect to au and (7) yields
(8)dGMGM⊙=2×10-6.
Thus, we conclude that the technique adopted to determine the numerical values of the Astronomical Unit and of the Sun's GM induced an a priori “imprint” of GTR on them of 8×10-7 and 2×10-6, respectively.Let us apply this result to a typical radio science experiment in the solar system withr1 fixed to the Earth-Sun distance. By writing r1/2=x1/2 au, with x1/2 expressing distances in Astronomical Units, differentiation of (6) with respect to au and GM⊙, and (7)-(8) yields an “imprinting” effect of the order of
(9)δ(Δt)Δt|GTR=2×10-6
for r2 up to tens of AU; it turns out that the largest contribution comes from dGM⊙. (The Cassini test was performed with r2=7.43 au [2].)It is too small by one order of magnitude with respect to the performed Cassini radio science tests, but it should be taken into account in the future, more accurate experiments whose expected accuracy is of the order of 10-7–10-9, in the sense that the a priori bias of GTR in the future determinations of deviations of γ from unity will be as large as, or even larger than the effects one will test, unless either GM⊙ will be estimated as well along with γ itself or a value obtained independently of it will be adopted.About the first point, we mention thatγ and β were never estimated in those solutions in which, among other things, also AU was fitted; instead, they were kept fixed [15–19]. Moreover, Fienga et al. [19] used a modified version of their latest planetary ephemerides, named INPOP08d, in which they kept fixed AU and simultaneously fitted GM⊙ to the observations along with other parameters, but they did not simultaneously estimate the PPN parameters as well, which were kept fixed. Instead, for the purposes outlined in this paper it would be important to simultaneously fit γ and AU (or GM⊙) and inspect the covariance matrix for the correlation among AU (or GM⊙) and γ. (Recently, W. M. Folkner (private communication to the author, March 2010) produced a global solution aimed to specifically test general relativity with the latest version of the JPL DE ephemerides. In it β, γ and AU were simultaneously fitted obtaining a correlation between γ and AU of  -1.34×10-2.)Concerning the last point, a possible choice may consist, in principle, of using a figure forGM⊙ obtained from measurements of the solar gravitational redshift
(10)Z≡νrec-νemνrec≈GM⊙R⊙c2
which, at first order, is independent of GTR itself. (Indeed, it depends on the coefficient g00 of the spacetime metric tensor; only the PPN parameter β enters g00 with the term of order 𝒪(c-4), at present, undetectable.)Latest measurements of the IR oxygen triplet 7772–7775, extrapolated to the Sun's limb [20], yield an accuracy of the order of a few percent levels. Projects to improve it with future missions have been proposed [21, 22]; for example, the use of the Magneto-Optical Filter technique, developed by Cacciani et al. [21], would allow reaching a relative accuracy of 10-6. Anyway, it must be noted that extracting GM⊙ from the measured gravitational redshift also requires the knowledge of the Sun's radius R⊙, which is uncertain at a 10-4 level. Indeed, the commonly accepted value for the solar radius was for a long time [23]
(11)R⊙=695.99Mm(1Mm=106m), although it was not clear how such a figure was obtained; moreover, no error bar was released. Subsequent observations of solar f-mode frequencies were induced by Schou et al. [24] and Antia [25], who used data from the Michelson Doppler Imager (MDI) mounted on the Solar and Heliospheric Observatory (SOHO) satellite and from the Global Oscillation Network Group (GONG) network; to conclude that the actual solar radius was 0.3-0.2 Mm smaller, indeed, the estimate by Schou et al. [24] is
(12)R⊙=695.68±0.02Mm.
In the work of Brown and Christensen-Dalsgaard [26] the value
(13)R⊙=695.509±0.026Mm
is reported from the Earth-based High Altitude Observatory's Solar Diameter Monitor campaign; it differs by about 0.5 Mm from the standard value by Allen [23]. In the work of Takata and Gough [27] we find, from SOHO/MDI frequency data for p-mode frequencies,
(14)R⊙=695.69±0.14Mm.
Note that the values by Brown and Christensen-Dalsgaard [26] and Takata and Gough [27] are mutually inconsistent, although at a 1-σ level only, the estimates by Schou et al. [24] and Brown and Christensen-Dalsgaard [26] are not consistent at 3-σ level. The f-mode and p-mode measurements [24, 27] are, instead, consistent with each other. By the way, it can be noted that the discrepancies among the different best estimates are larger than the associated errors; thus, as a conservative evaluation of the accuracy in knowing the solar radius we will assume such differences.


## 3.1. The Cassini Radio Science Test and Other Proposed High-Accuracy Space-Based Measurements ofγ

To be more definite, let us look at the Cassini radio science test. In that case, the radiotechnical data of the spacecraft traveling to Saturn were contrasted with a set of dynamical models by JPL of its motion and electromagnetic waves propagation in such a way that a correctionΔγ to the GTR-predicted value of the PPN parameter γ was solved for, among other parameters, obtaining [2]
(1)Δγ≡|γ-1|=(2.1±2.3)×10-5;
other authors got [3]
(2)Δγ≡|γ-1|=(-1.3±5.2)×10-5.Now, a physical parameter which is, of course, crucial in such a test is the gravitational parameterGM⊙ of the Sun, which is the source of the relativistic time delay put on the test. It was not estimated [2, 3], so that its numerical value was kept fixed to the standard reference figure of the JPL DE ephemerides used. It does, in principle, contain an a priori “imprinting” by GTR itself through the same effect itself that was just tested with Cassini, in particular by γ itself. Indeed, the numerical value of GM⊙ comes from the fixed value of the defining Gaussian constant (see on the WEB http://ssd.jpl.nasa.gov/?constants.) (3)k=0.01720209895au3/2d-1,
and from the value of the Astronomical Unit, not estimated in the Cassini tests,
(4)AU=1.49597870691×1011(±3)m
through
(5)GM⊙=k2AU3d-2=1.32712440018×1020(±8×109)m3s-2.(Here we will use au for the symbol of the Astronomical Unit, like m for the meter, while AU will denote its numerical value in m.)AU was, in fact, obtained just through a combination of radar ranging of Mercury, Venus, and Mars, laser ranging of the moon (making use of light reflectors left on the lunar surface by Apollo astronauts), and timing of signals returned from spacecrafts as they orbit or make close passes and objects in the solar system [14]; thus, it is affected in a nonnegligible way, given the level of accuracy of the techniques adopted, by GTR itself and, in particular, by γ which enters the PPN expressions for the time delay and the bending of traveling electromagnetic waves. Thus, there exists, in principle, the possibility that the high-accuracy results of the Cassini radio science tests may retain an a priori “imprint” of GTR itself through GM⊙ (and the Astronomical Unit as well).Let us put our hypothesis on the test by making some concrete calculations; for the sake of clarity, we will refer to the Cassini radio science tests, but the conclusions may be considered valid also for any of the many proposedγ-dedicated missions.The GTR time delay experienced by electromagnetic waves propagating from point 1 to point 2 is(6)Δt=2Rgcln(r1+r2+r12+Rgr1+r2-r12+Rg),
where Rg=2GM⊙/c2 is the Sun's Schwarzschild radius in which G is the Newtonian gravitational constant, M⊙ is the solar mass, c is the speed of light in vacuum, r1 is the heliocentric coordinate distance to point 1, r2 is the heliocentric coordinate distance to point 2, and r12 is the distance between the points 1 and 2. Equation (6) is the expression actually used in the JPL's Orbit Determination Program (ODP) to analyze interplanetary ranging to planets and probes. In order to quantitatively evaluate the level of “imprinting” by GTR itself in the used value of the Astronomical Unit, let us assume that r1 equal to the Earth-Sun distance and let us vary r2 within 0.38 au and 1.5 au to account for the ranging to inner planets; the maximum effect occurs at the superior conjunction, that is, when n1≈-n2 and r12≈r1+r2.(Here n≡r/r.) It turns out that Δtranging≈4×10-4 s, which is certainly not negligible with respect to the accuracy of the order of 10-8 s with which the light time for 1 au τA is actually measured (http://ssd.jpl.nasa.gov/?constants).(The value in km of the Astronomical Unit is obtained by measuring at a given epoch the distance between the Earth and a target body (a planet or a probe orbiting it) by multiplying c times the round trip travel time τ of electromagnetic waves sent from the Earth and reflected back by the target body, and confronting it with the distance, expressed in AU, between the Earth and the target body at the same epoch as predicted by some accurate dynamical ephemeris [14].) As a consequence, the quantitative impact of the interplanetary ranging in the inner solar system on the determination of the Astronomical Unit is of the order
(7)dAU=cΔtranging=1.14291×105m,
not negligible with respect to the meter-level accuracy in measuring the Astronomical Unit; thus, dAU/AU=8×10-7. Differentiating (5) with respect to au and (7) yields
(8)dGMGM⊙=2×10-6.
Thus, we conclude that the technique adopted to determine the numerical values of the Astronomical Unit and of the Sun's GM induced an a priori “imprint” of GTR on them of 8×10-7 and 2×10-6, respectively.Let us apply this result to a typical radio science experiment in the solar system withr1 fixed to the Earth-Sun distance. By writing r1/2=x1/2 au, with x1/2 expressing distances in Astronomical Units, differentiation of (6) with respect to au and GM⊙, and (7)-(8) yields an “imprinting” effect of the order of
(9)δ(Δt)Δt|GTR=2×10-6
for r2 up to tens of AU; it turns out that the largest contribution comes from dGM⊙. (The Cassini test was performed with r2=7.43 au [2].)It is too small by one order of magnitude with respect to the performed Cassini radio science tests, but it should be taken into account in the future, more accurate experiments whose expected accuracy is of the order of 10-7–10-9, in the sense that the a priori bias of GTR in the future determinations of deviations of γ from unity will be as large as, or even larger than the effects one will test, unless either GM⊙ will be estimated as well along with γ itself or a value obtained independently of it will be adopted.About the first point, we mention thatγ and β were never estimated in those solutions in which, among other things, also AU was fitted; instead, they were kept fixed [15–19]. Moreover, Fienga et al. [19] used a modified version of their latest planetary ephemerides, named INPOP08d, in which they kept fixed AU and simultaneously fitted GM⊙ to the observations along with other parameters, but they did not simultaneously estimate the PPN parameters as well, which were kept fixed. Instead, for the purposes outlined in this paper it would be important to simultaneously fit γ and AU (or GM⊙) and inspect the covariance matrix for the correlation among AU (or GM⊙) and γ. (Recently, W. M. Folkner (private communication to the author, March 2010) produced a global solution aimed to specifically test general relativity with the latest version of the JPL DE ephemerides. In it β, γ and AU were simultaneously fitted obtaining a correlation between γ and AU of  -1.34×10-2.)Concerning the last point, a possible choice may consist, in principle, of using a figure forGM⊙ obtained from measurements of the solar gravitational redshift
(10)Z≡νrec-νemνrec≈GM⊙R⊙c2
which, at first order, is independent of GTR itself. (Indeed, it depends on the coefficient g00 of the spacetime metric tensor; only the PPN parameter β enters g00 with the term of order 𝒪(c-4), at present, undetectable.)Latest measurements of the IR oxygen triplet 7772–7775, extrapolated to the Sun's limb [20], yield an accuracy of the order of a few percent levels. Projects to improve it with future missions have been proposed [21, 22]; for example, the use of the Magneto-Optical Filter technique, developed by Cacciani et al. [21], would allow reaching a relative accuracy of 10-6. Anyway, it must be noted that extracting GM⊙ from the measured gravitational redshift also requires the knowledge of the Sun's radius R⊙, which is uncertain at a 10-4 level. Indeed, the commonly accepted value for the solar radius was for a long time [23]
(11)R⊙=695.99Mm(1Mm=106m), although it was not clear how such a figure was obtained; moreover, no error bar was released. Subsequent observations of solar f-mode frequencies were induced by Schou et al. [24] and Antia [25], who used data from the Michelson Doppler Imager (MDI) mounted on the Solar and Heliospheric Observatory (SOHO) satellite and from the Global Oscillation Network Group (GONG) network; to conclude that the actual solar radius was 0.3-0.2 Mm smaller, indeed, the estimate by Schou et al. [24] is
(12)R⊙=695.68±0.02Mm.
In the work of Brown and Christensen-Dalsgaard [26] the value
(13)R⊙=695.509±0.026Mm
is reported from the Earth-based High Altitude Observatory's Solar Diameter Monitor campaign; it differs by about 0.5 Mm from the standard value by Allen [23]. In the work of Takata and Gough [27] we find, from SOHO/MDI frequency data for p-mode frequencies,
(14)R⊙=695.69±0.14Mm.
Note that the values by Brown and Christensen-Dalsgaard [26] and Takata and Gough [27] are mutually inconsistent, although at a 1-σ level only, the estimates by Schou et al. [24] and Brown and Christensen-Dalsgaard [26] are not consistent at 3-σ level. The f-mode and p-mode measurements [24, 27] are, instead, consistent with each other. By the way, it can be noted that the discrepancies among the different best estimates are larger than the associated errors; thus, as a conservative evaluation of the accuracy in knowing the solar radius we will assume such differences.


## 4. Conclusions

We have investigated the impact of possible a priori “imprinting” effects of GTR itself on satellite/spacecraft data analyses specifically designed to test some general relativistic predictions. In particular, we considered the time-delay experiments conducted or proposed in the Sun's field with ranging to interplanetary spacecraft.Concerning the time-delay tests, the numerical values for the Astronomical Unit and the solar gravitational constant currently adopted in the ephemerides used so far to analyze spacecrafts' data retain an a priori “imprint” of GTR itself of the order of8×10-7 and 2×10-6, respectively. As a consequence, the bias in typical solar system radioscience experiments is of the order of 10-6, which is one order of magnitude smaller than the accuracy level reached in the performed Cassini experiment, but it would be of concern for future planned tests aiming to measure deviations of the PPN parameter γ from its general relativistic value at 10-7–10-9 level.In order to have genuine, unambiguous, and unbiased tests of GTR which are not “tautologic”, it would be necessary to either estimate the suspect parameters as well along with those accounting for the relativistic effect of interest or use values for them obtained independently from GTR itself.



---

*Source: 735487-2010-03-18.xml*