*Koji MATSUMOTO1, Ryuhei YAMADA1, Fuyuhiko KIKUCHI1, Shunichi KAMATA2, Takahiro IWATA3, Yoshiaki ISHIHARA3, Hideo HANADA1, Sho SASAKI4
(1.RISE Project Office, NAOJ, 2.University of California Santa Cruz, 3.JAXA, 4.Osaka University)
Keywords:Moon, internal structure, gravity field, tidal Love number, GRAIL
Internal structure and composition of the Moon provide important clue and constraints on theories for how the Moon formed and evolved. The Apollo seismic network has contributed to the internal structure modeling. Efforts have been made to detect the lunar core from the noisy Apollo data (e.g., [1], [2]), but there is scant information about the structure below the deepest moonquakes at about 1000 km depth. On the other hand, there have been geodetic studies to infer the deep structure of the Moon. For example, LLR (Lunar Laser Ranging) data analyses detected a displacement of the lunar pole of rotation, indicating that dissipation is acting on the rotation arising from a fluid core [3]. Bayesian inversion using geodetic data weakly suggests a fluid core and partial melt in the lower mantle region [4]. Further improvements in determining the second-degree gravity coefficients and the Love numbers will help us to better constrain the lunar internal structure. Recent analyses of GRAIL data have achieved the improved k2 accuracy; JPL solution is 0.02405±0.00018 [5], and GSFC solution is 0.02427±0.00026 [6]. The two solutions are consistent with each other within their error bounds, and the accuracy of k2 is now about 1 %. By introducing the improved gravity coefficients and k2 from GRAIL mission, the updated LLR data analysis has also resulted in a better h2 determination. Such accurately-determined Love numbers will contribute to constrain the structure of the lunar deep interior, such as the radius of the possible liquid core. It is difficult, however, to tightly constrain the internal structure from the geodetic data only because there are trade-offs among the structures of crust, mantle, and core. The combination of the Apollo seismic data and the geodetic data therefore afford the key to better determination of the lunar interior structure. We included geodetic data of the mass, the mean moment of inertia, the Love numbers h2 and k2, and 262 P and S travel time data in the analysis.Markov Chain Monte Carlo (MCMC) method is used to infer the model parameters. When we used a five-layer model consisting of crust, upper-mantle, mid-mantle, lower-mantle, and core, the core radius is estimated to be 483±22 km, and the core density values tend to be sampled around the assumed lower limit of 3600 kg/m3. However, the inferred core radius is significantly larger than the magnetic constraint from SELENE data [7] which predicts the upper bound of the core radius to be 400 km. This discrepancy might be attributed to a possible low velocity layer above the core-mantle boundary which was not included in the five-layer model. We will discuss the results when such a low velocity layer is taken into account. [1] Weber et al. (2011), Science, 331, 309-312, doi:10.1126/science.1199375[2] Garcia et al. (2011), PEPI, doi:10.1016/j.pepi.2011.06.015[3] Williams et al. (2001), JGR, 106, E11, 27,933-27,968[4] Khan and Mosegaard (2005), GRL, 32, L22203, doi:10.1029/2005GL023985[5] Konopliv et al (2013), JGR, 118, doi:10.1002/jgre.20097[6] Taken from the PDS label of GRAIL Derived Data Products[7] Simizu et al. (2013), Icarus, 222, doi:10.1016/j.icarus.2012.10.029