Gravitation of the Moon

The gravitational field of the Moon has been measured by tracking the radio signals emitted by orbiting spacecraft. The principle used depends on the Doppler effect, whereby the line-of-sight spacecraft acceleration can be measured by small shifts in frequency of the radio signal, and the measurement of the distance from the spacecraft to a station on Earth. Since the gravitational field of the Moon affects the orbit of a spacecraft, one can use this tracking data to detect gravity anomalies.

Most low lunar orbits are unstable. Detailed data collected has shown that for low lunar orbit the only "stable" orbits are at inclinations near 27°, 50°, 76°, and 86°.^[2] Because of the Moon's synchronous rotation it is not possible to track spacecraft from Earth much beyond the limbs of the Moon, so until the recent Gravity Recovery and Interior Laboratory (GRAIL) mission the far-side gravity field was not well mapped.

The missions with accurate Doppler tracking that have been used for deriving gravity fields are in the accompanying table. The table gives the mission spacecraft name, a brief designation, the number of mission spacecraft with accurate tracking, the country of origin, and the time span of the Doppler data. Apollos 15 and 16 released subsatellites. The Kaguya/SELENE mission had tracking between 3 satellites to get far-side tracking. GRAIL had very accurate tracking between 2 spacecraft and tracking from Earth.

The accompanying table below lists lunar gravity fields. The table lists the designation of the gravity field, the highest degree and order, a list of mission IDs that were analyzed together, and a citation. Mission ID LO includes all 5 Lunar Orbiter missions. The GRAIL fields are very accurate; other missions are not combined with GRAIL.

A major feature of the Moon's gravitational field is the presence of mascons, which are large positive gravity anomalies associated with some of the giant impact basins. These anomalies significantly influence the orbit of spacecraft around the Moon, and an accurate gravitational model is necessary in the planning of both crewed and uncrewed missions. They were initially discovered by the analysis of Lunar Orbiter tracking data:^[16] navigation tests prior to the Apollo program showed positioning errors much larger than mission specifications.

Mascons are in part due to the presence of dense mare basaltic lava flows that fill some of the impact basins.^[17] However, lava flows by themselves cannot fully explain the gravitational variations, and uplift of the crust-mantle interface is required as well. Based on Lunar Prospector gravitational models, it has been suggested that some mascons exist that do not show evidence for mare basaltic volcanism.^[3] The huge expanse of mare basaltic volcanism associated with Oceanus Procellarum does not cause a positive gravity anomaly. The center of gravity of the Moon does not coincide exactly with its geometric center, but is displaced toward the Earth by about 2 kilometers.^[18]

The gravitational constant G is less accurate than the product of G and masses for Earth and Moon. Consequently, it is conventional to express the lunar mass M multiplied by the gravitational constant G. The lunar GM = 4902.8001 km³/s² from GRAIL analyses.^[12][11][19] The mass of the Moon is M = 7.3458 × 10²² kg and the mean density is 3346 kg/m³. The lunar GM is 1/81.30057 of the Earth's GM.^[20]

For the lunar gravity field, it is conventional to use an equatorial radius of R = 1738.0 km. The gravity potential is written with a series of spherical harmonic functions P_nm. The gravitational potential V at an external point is conventionally expressed as positive in astronomy and geophysics, but negative in physics. Then, with the former sign,

$${\displaystyle V=\left({\frac {GM}{r}}\right)-\left({\frac {GM}{r}}\right)\sum \left({\frac {R}{r}}\right)^{n}J_{n}P_{n,0}(\sin \phi )+\left({\frac {GM}{r}}\right)\sum \left({\frac {R}{r}}\right)^{n}[C_{n,m}P_{n,m}(\sin \phi )\cos(m\lambda )+S_{n,m}P_{n,m}(\sin \phi )\sin(m\lambda )]}$$ 

where r is the radius to an external point with r ≥ R, φ is the latitude of the external point, and λ is the east longitude of the external point. Note that the spherical harmonic functions P_nm can be normalized or unnormalized affecting the gravity coefficients J_n, C_nm, and S_nm. Here we will use unnormalized functions and compatible coefficients. The P_n0 are called Legendre polynomials and the P_nm with m≠0 are called the Associated Legendre polynomials, where subscript n is the degree, m is the order, and m ≤ n. The sums start at n = 2. The unnormalized degree-2 functions are

$${\displaystyle {\begin{aligned}P_{2,0}&={\frac {3}{2}}\sin ^{2}\!\phi -{\frac {1}{2}}\\[1ex]P_{2,1}&=3\sin \phi \cos \phi \\[1ex]P_{2,2}&=3\cos ^{2}\phi \end{aligned}}}$$ 

Note that of the three functions, only P₂₀(±1)=1 is finite at the poles. More generally, only P_n0(±1)=1 are finite at the poles.

The gravitational acceleration of vector position r is

$${\displaystyle {\begin{aligned}{\frac {d^{2}r}{dt^{2}}}&=\nabla V\\[1ex]&={\partial V \over \partial r}e_{r}+{\frac {1}{r}}{\partial V \over \partial \phi }e_{\phi }+{\frac {1}{r\cos \phi }}{\partial V \over \partial \lambda }{e_{\lambda }}\end{aligned}}}$$ 

where e_r, e_φ, and e_λ are unit vectors in the three directions.

The unnormalized gravity coefficients of degree 2 and 3 that were determined by the GRAIL mission are given in Table 1.^[12][11][19] The zero values of C₂₁, S₂₁, and S₂₂ are because a principal axis frame is being used. There are no degree-1 coefficients when the three axes are centered on the center of mass.

The J₂ coefficient for an oblate shape to the gravity field is affected by rotation and solid-body tides whereas C₂₂ is affected by solid-body tides. Both are larger than their equilibrium values showing that the upper layers of the Moon are strong enough to support elastic stress. The C₃₁ coefficient is large.

In January 2022 China was reported by the South China Morning Post to have built a small (60 centimeters in diameter) research facility to simulate low lunar gravity with the help of magnets.^[21][22] The facility was reportedly partly inspired by the work of Andre Geim (who later shared the 2010 Nobel Prize in Physics for his research on graphene) and Michael Berry, who both shared the Ig Nobel Prize in Physics in 2000 for the magnetic levitation of a frog.^[21][22]

• Magnetic field of the Moon
• Micro-g environment