High‒degree gravity models from GRAIL primary mission data

• Published in: Journal of geophysical research. Planets Vol. 118; no. 8; pp. 1676 - 1698
• Main Authors: Lemoine, Frank G., Goossens, Sander, Sabaka, Terence J., Nicholas, Joseph B., Mazarico, Erwan, Rowlands, David D., Loomis, Bryant D., Chinn, Douglas S., Caprette, Douglas S., Neumann, Gregory A., Smith, David E., Zuber, Maria T.
• Format: Journal Article
• Language: English
• Published: 01.08.2013
• Subjects: KBRR data; lunar gravity; orbit determination; planetary geodesy
• DOI: 10.1002/jgre.20118
• ISSN: 2169-9097

We have analyzed Ka‒band range rate (KBRR) and Deep Space Network (DSN) data from the Gravity Recovery and Interior Laboratory (GRAIL) primary mission (1 March to 29 May 2012) to derive gravity models of the Moon to degree 420, 540, and 660 in spherical harmonics. For these models, GRGM420A, GRGM540A, and GRGM660PRIM, a Kaula constraint was applied only beyond degree 330. Variance‒component estimation (VCE) was used to adjust the a priori weights and obtain a calibrated error covariance. The global root‒mean‒square error in the gravity anomalies computed from the error covariance to 320×320 is 0.77 mGal, compared to 29.0 mGal with the pre‒GRAIL model derived with the SELENE mission data, SGM150J, only to 140×140. The global correlations with the Lunar Orbiter Laser Altimeter‒derived topography are larger than 0.985 between ℓ=120 and 330. The free‒air gravity anomalies, especially over the lunar farside, display a dramatic increase in detail compared to the pre‒GRAIL models (SGM150J and LP150Q) and, through degree 320, are free of the orbit‒track‒related artifacts present in the earlier models. For GRAIL, we obtain an a posteriori fit to the S‒band DSN data of 0.13 mm/s. The a posteriori fits to the KBRR data range from 0.08 to 1.5 μm/s for GRGM420A and from 0.03 to 0.06 μm/s for GRGM660PRIM. Using the GRAIL data, we obtain solutions for the degree 2 Love numbers, k20=0.024615±0.0000914, k21=0.023915±0.0000132, and k22=0.024852±0.0000167, and a preliminary solution for the k30 Love number of k30=0.00734±0.0015, where the Love number error sigmas are those obtained with VCE. Key Points Analyze satellite‐to‐satellite and DSN data to the GRAIL spacecraft Determine gravity field solutions to 660x660 in spherical harmonics Use variance component estimation to calibrate solution errors