Patched Local Lunar Gravity Solutions Using GRAIL Data

We present a method to determine local gravity fields for the Moon using Gravity Recovery and Interior Laboratory (GRAIL) data. We express gravity as gridded gravity anomalies on a sphere, and we estimate adjustments to a background global start model expressed in spherical harmonics. We processed GRAIL Ka‐band range‐rate data with a short‐arc approach, using only data over the area of interest. We determine our gravity solutions using neighbor smoothing constraints. We divided the entire Moon into 12 regions and 2 polar caps, with a resolution of [Formula: see text] (which is equivalent to degree and order 1199 in spherical harmonics), and determined the optimal smoothing parameter for each area by comparing localized correlations between gravity and topography for each solution set. Our selected areas share nodes with surrounding areas and they are overlapping. To mitigate boundary effects, we patch the solutions together by symmetrically omitting the boundary parts of overlapping solutions. Our new solution has been iterated, and it has improved correlations with topography when compared to a fully iterated global model. Our method requires fewer resources, and can easily handle regionally varying resolution or constraints. The smooth model describes small‐scale features clearly, and can be used in local studies of the structure of the lunar crust.