Issue 
A&A
Volume 522, November 2010



Article Number  L5  
Number of page(s)  3  
Section  Letters  
DOI  https://doi.org/10.1051/00046361/201015659  
Published online  29 October 2010 
Letter to the Editor
Lunar laser ranging test of the Nordtvedt parameter and a possible variation in the gravitational constant
^{1}
Institut für Erdmessung (IfE), Leibniz Universität Hannover,
Scheiderberg 50,
30167
Hannover,
Germany
email: hofmann@ife.unihannover.de
^{2}
QUEST – Centre for Quantum Engineering and Space Time Research,
Leibniz Universität Hannover, Welfengarten 1, 30167
Hannover,
Germany
Received:
30
August
2010
Accepted:
19
October
2010
Context. Forty years of lunar laser ranging (LLR) data provide an excellent basis to determine various parameters of the EarthMoon system as well as parameters related to gravitational physics.
Aims. We update the Institut für Erdmessung (IfE) LLR model taking the effect of a fluid lunar core into consideration. The temporal variation in the gravitational constant Ġ/G_{0} and the strong equivalence principle, parameterized by the Nordtvedt parameter η, are investigated.
Methods. A set of LLR observations from 1970 to 2009 was analysed and the parameters were determined by a least squares adjustment in two runs. After solving for classical Newtonian parameters (e.g. initial conditions for the lunar orbit and rotation) in the first run, relativistic parameters were determined in the second run.
Results. The upper limits to the gravitational constant and the Nordtvedt parameter were found to be Ġ/G_{0} = (−0.7 ± 3.8) × 10^{13} yr^{1} and η = (−0.6 ± 5.2) × 10^{4}.
Key words: Moon / gravitation / celestial mechanics / ephemerides
© ESO, 2010
1. Introduction
Since the first passive retroreflector on the lunar surface was deployed by the Apollo 11 mission, measurements of round trip travel times of laser pulses between stations on the Earth and retroreflectors on the Moon have been carried out for over 40 years. The weak signaltonoise ratio of the pulses is mainly caused by atmospheric beam divergence and corner cube diffraction, and makes lunar laser ranging (LLR) challenging. Therefore, only a few observatories worldwide are capable of tracking the Moon. The measured lunar returns over a period of 15 to 20 min are combined to socalled normal points (NP). Up to now, approximately 17 000 NPs were collected and serve as observations for the analysis. Various components of the EarthMoon system dynamics can be addressed, e.g., initial conditions for the lunar orbit and rotation, reflector and station coordinates, Earth orientation parameters (Biskupek & Müller 2009; Zerhouni & Capitaine 2009), nutation coefficients (Biskupek et al. 2009), and also studies of the lunar interior (Williams & Boggs 2009). In addition, LLR is used to test general relativity (Müller et al. 2008), e.g., the strong equivalence principle (EP) (Müller & Nordtvedt 1998; Williams et al. 2009; Nordtvedt et al. 1995), temporal variation in the gravitational constant G (Müller & Biskupek 2007; Williams & Folkner 2009), metric parameters, and preferred frame effects (Soffel et al. 2008). Here, we focus on the strong equivalence principle and a possible temporal variability of the gravitational constant after refining the analysis model for the lunar rotation taking into account the effect of a fluid core and an updated data series (March 1970 to December 2009). The knowledge of more tighter limits of these quantities may be helpful to verify Einstein’s theory or to distinguish between cosmological models and alternative theories of gravitation. Some theories predict for instance an EP violation at a level of m_{g}/m_{i} = 10^{18}–10^{12} with gravitational mass m_{g} and inertial mass m_{i} (Damour et al. 2002; Turyshev 2008), while some multidimensional theories expect a temporal variation in G of the order of Ġ/G = 10^{14}–10^{11} yr^{1} (Steinhardt & Wesley 2010; Sanders et al. 2010).
This Letter is organized as follows. Section 2 describes the analysis model and Sect. 3 shows the basis of the parameter estimation with a special focus on the investigated relativistic parameters. The results obtained are discussed in Sect. 4 and our conclusions are given in Sect. 5.
2. Analysis model
The analysis model used at Institut für Erdmessung (IfE) is based on Einstein’s theory of gravity. It is fully relativistic and complete up to the first postNewtonian (1/c^{2}) level (see Müller et al. 2008). The simplified observation equation for the stationreflector distance d (see Fig. 1) is (1)with the speed of light c, pulse travel time τ, vector between geocenter and selenocenter r_{EM}, the position vector of the observatories r_{obs}, the position vector of the reflectors r_{ref}, and a delay Δτ due to atmospheric and relativistic effects.
Fig. 1 LLR measurement setup with basis vectors. 
The inclusion of the lunar core refines the modelling of the lunar rotation and leads to an improved LLR analysis model. The root mean square of the postfit residuals of the whole data series yields a reduction of 4 mm or 1.1% with the core extension.
3. Parameter estimation
The model is adjusted to the observations by performing a weighted least squares fit with two parameter groups. The first group are the nonrelativistic Newtonian parameters, which are estimated in a first run of the analysis, with fixed relativistic parameters to their Einsteinian values, until a stable solution is obtained. These parameters include coordinates of LLR observatories and retroreflectors, GM_{E + M}, long periodic nutation coefficients, initial conditions for the lunar orbit and rotation, the lowest mass multipole moments of the Moon, lunar Love number, as well as a dissipation parameter and a lag angle indicating the lunar tidal acceleration. After convergence is achieved, single relativistic parameters are estimated together with the Newtonian parameters. All other relativistic parameters are fixed to their Einsteinian values. Tests relating to the gravitational constant and the EP are discussed in more detail below.
3.1. Gravitational constant
The gravitational constant G entering into the equations of motion is described as (8)where Δt = t − t_{0} is the time difference to the reference epoch t_{0} = J2000 and Ġ/G_{0} as linear trend of the gravitational constant. This yields Ġ/G_{0} = (−0.9 ± 4.0) × 10^{13} yr^{1} when only the updated data series is considered and Ġ/G_{0} = (−0.7 ± 3.8) × 10^{13} yr^{1} including the core extension. Owing to modelling inadequacies and remaining systematic errors, the formal error from the fit was scaled by a factor of 2 to get a more realistic uncertainty. The scaling factor was inferred by different tests from the determined uncertainty in Ġ/G_{0}, e.g., by varying values of some fixed Newtonian parameter, such as the EarthMoon orbit around the Sun, or the weighting of the individual measurements. For a comprehensive error analysis, we refer to Müller & Biskupek (2007).
3.2. Nordtvedt parameter
Compared to laboratory masses, Earth and Moon have a large amount of gravitational selfenergy. Therefore, they can be used as bodies for testing the strong EP with respect to the gravitational field of the Sun (Nordtvedt 1968a,b). If the EP is violated (Nordtvedt parameter η ≠ 0), the Earth and the Moon will be accelerated differently by the Sun, resulting in a shifted lunar orbit in the solar direction. In the IfE LLR analysis program, this can be handled by determining the difference of the ratio of gravitational to inertial mass Δ(m_{g}/m_{i}) between the Earth and Moon in the equations of motion or by synodic EarthMoon range oscillations towards the sun, parameterized by η, which leads to an equivalent result (Müller & Nordtvedt 1998). In the present study, only the second procedure with the range term Δr_{EM} was considered (9)where D is the synodic angle and the amplitude of A = 12.8 m is derived from theoretical calculations by Nordtvedt (1995). This yields η = (2.1 ± 5.3) × 10^{4} considering only the updated data series and η = (−0.6 ± 5.2) × 10^{4} with the core extension. When using the realistic uncertainty value of η = 5.2 × 10^{4}, Eq. (9) would lead to a maximum additional range of 7 mm along the EarthSun direction during the new and full Moon phases. For the mass ratio Δ(m_{g}/m_{i}), this yields Δ(m_{g}/m_{i}) = −4.45 × 10^{10}η = 0.3 ± 2.3 × 10^{13}, where we refer to Williams et al. (2009) for more details and Müller & Nordtvedt (1998) for a dedicated error analysis: this procedure is carried out in a similar way to that indicated in Sect. 3.1.
4. Discussion
We have succeeded in reducing the uncertainty in the current variation in G by a factor of 2, compared to previous IfE results (Müller et al. 2008). This is mainly due to the longer data series used and to a smaller extent a result of the refined modelling. The analysis shows a high correlation of Ġ/G_{0} with the lag angle, which indicates the secular tidal acceleration. At the moment, this prevents a more precise determination of Ġ/G_{0}. Together with Williams & Folkner (2009), this paper gives the actual highest accuracy of Ġ/G_{0} ≈ 4 × 10^{13} yr^{1} from pure LLR observations. Compared to measurements based on radar observations (Pitjeva 2005), the LLR results seem to be less accurate by a factor of about 8.
The Nordtvedt parameter has been improved by a factor of 1.5 and agrees with the value derived by Williams et al. (2009). This is again mainly a result of the longer data series and to a smaller extent the refined modelling. At the epochs of possible maximum amplitude, there are either no (new Moon) or fewer (full Moon) observations because of the too high background noise during the measurements. This means that, at the moment, the full potential of the EP test with LLR is not being exploited. New measurements, e.g., by the APOLLO station (Murphy et al. 2008), at full Moon phases will certainly lead to further improvement. In the parametrized postNewtonian (PPN) framework, η is related to the PPN parameter β and γ by η = 4β − γ −3. This equation is valid under the assumption of no compositioninduced EP violation and no preferred frame effects. By taking η from this paper and the parameter γ derived from the Cassini mission (Bertotti et al. 2003), the PPN parameter β has been determined to be β − 1 = (0.3 ± 1.3) × 10^{4}.
Both relativistic values discussed in this paper do not show any significant deviation from their Einsteinian values (Ġ/G_{0} = 0, η = 0) and confirm the upper limit to a possible violation of general relativity.
5. Conclusions
Owing to the long data series of observations, LLR has become a powerful tool for determining relativistic parameters with high accuracy. Here, the variations in the gravitational constant and the Nordtvedt parameter have been investigated by applying an analysis model that had been extended by incorporating the fluid lunar core influence in combination with an updated data series up to December 2009. We obtained Ġ/G_{0} = (−0.7 ± 3.8) × 10^{13} yr^{1} and η = (−0.6 ± 5.2) × 10^{4}. The results of these parameters have been improved by a factor of about 2
relative to previous IfE results, and confirm or even slightly improve the results of other analysis groups for an upper limit to a possible violation of general relativity.
Acknowledgments
Current LLR data are collected, archived and distributed under the auspices of the International Laser Ranging Service (ILRS). We acknowledge with thanks, that the more than 40 years of used LLR data has been obtained under the efforts of the personnel at the Observatoire de la Côte d’Azur in France (9067 NPs used in this analysis), the LURE Observatory in Maui, Hawaii (402 NPs), the McDonald Observatory in Texas (6551 NPs) as well as the Apache Point Observatory in New Mexico (600 NPs). We would also like to thank the International Space Science Institute (ISSI) in Bern for supporting this research. This research was funded by the Centre for Quantum Engineering and SpaceTime Research QUEST and the DFG, the German Research Foundation, within the research unit FOR584 “Earth rotation and global dynamic processes”.
References
 Bertotti, B., Iess, L., & Tortora, P. 2003, Nature, 425, 374 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Biskupek, L., & Müller, J. 2009, in Proceedings of the Journées 2008 Systèmes de référence spatiotemporels, ed. M. Soffel, & N. Capitaine, LohrmannObservatorium and Observatoire de Paris, 182–185, les Journées Systèmes de référence spatiotemporels and X. LohrmannKolloquium, Dresden, Germany, 22–24.09.2008 [Google Scholar]
 Biskupek, L., Müller, J., & Hofmann, F. 2009, in Proc. of IAG 2009 Scientific Assembly Geodesy for Planet Earth, International Association of Geodesy Symposia (Springer), accepted [Google Scholar]
 Damour, T., Piazza, F., & Veneziano, G. 2002, Phys. Rev. Lett., 89, 081601 [Google Scholar]
 Hinderer, J., Legros, H., & Amalvict, M. 1982, Geophys. J. Int., 71, 303 [NASA ADS] [CrossRef] [Google Scholar]
 McCarthy, D. D., & Petit, G. (ed.) 2004, IERS Conventions (2003), (Verlag des Bundesamtes für Kartographie und Geodäsie) [Google Scholar]
 Müller, J., & Biskupek, L. 2007, Classical and Quantum Gravity, 24, 4533 [Google Scholar]
 Müller, J., & Nordtvedt, K. 1998, Phys. Rev. D, 58, 062001 [NASA ADS] [CrossRef] [Google Scholar]
 Müller, J., Williams, J. G., & Turyshev, S. G. 2008, in Lasers, Clocks and DragFree Control: Exploration of Relativistic Gravity in Space, ed. H. Dittus, C. Lämmerzahl, & S. G. Turyshev, Astrophysics and Space Science Library, 349, 457 [Google Scholar]
 Murphy, T. W., Adelberger, E. G., Battat, J. B. R., et al. 2008, PASP, 120, 20 [Google Scholar]
 Nordtvedt, K. 1968a, Phys. Rev., 169, 1014 [NASA ADS] [CrossRef] [Google Scholar]
 Nordtvedt, K. 1968b, Phys. Rev., 169, 1017 [Google Scholar]
 Nordtvedt, K. 1995, Icarus, 114, 51 [NASA ADS] [CrossRef] [Google Scholar]
 Nordtvedt, K. L., Müller, J., & Soffel, M. 1995, A&A, 293, L73 [NASA ADS] [Google Scholar]
 Pitjeva, E. V. 2005, Astron. Lett., 31, 340 [NASA ADS] [CrossRef] [Google Scholar]
 Sanders, A. J., Gillies, G. T., & Schmutzer, E. 2010, [arXiv:1005.3226] [Google Scholar]
 Soffel, M., Klioner, S., Müller, J., & Biskupek, L. 2008, Phys. Rev. D, 78, 024033 [NASA ADS] [CrossRef] [Google Scholar]
 Steinhardt, P. J., & Wesley, D. 2010, [arXiv:1003.2815] [Google Scholar]
 Turyshev, S. G. 2008, Ann. Rev. Nucl. Part. Sci., 58, 207 [NASA ADS] [CrossRef] [Google Scholar]
 Williams, J. G., & Boggs, D. H. 2009, in Proc. of the 16th Int. Workshop on Laser Ranging, 13–17.08.2008, Poznań, Poland, ed. S. Schillak, 1, 101 [Google Scholar]
 Williams, J. G., & Folkner, W. M. 2009, IAU Symp. 261, BAAS, 41, 882 [Google Scholar]
 Williams, J. G., Boggs, D. H., Yoder, C. F., Ratcliff, J. T., & Dickey, J. O. 2001, J. Geophys. Res., 106, 27933 [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]
 Zerhouni, W., & Capitaine, N. 2009, A&A, 507, 1687 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Figures
Fig. 1 LLR measurement setup with basis vectors. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.