Gravity field of the Moon from the Gravity Recovery and Interior Laboratory (GRAIL) mission

Thermal asymmetry in the Moon's mantle inferred from monthly tidal response

The Moon undergoes periodic tidal forcing due to its eccentric and oblique orbit around the Earth1. The response to this tidal interaction drives temporal changes in the lunar gravity field and is sensitive to the satellite's internal structure2-4. We use data from the NASA GRAIL spacecraft5-9 to recover the time-varying lunar gravity field, including a degree-3 gravitational tidal Love number, k3. Here, we report our estimated value of k3 = 0.0163 ± 0.0007, which is about 72% higher than that expected for a spherically symmetric moon10. Such a large k3 can be explained if the elastic shear modulus of the mantle varies by about 2-3% between the nearside and farside4, providing an observational demonstration of lateral heterogeneities in the deep lunar interior. This asymmetric structure suggests preservation of a predominantly thermal anomaly of roughly 100-200 K in the nearside mantle that formed surface mare regions 3-4 billion years ago11 and could influence the spatial distribution of deep moonquakes12.

Io's tidal response precludes a shallow magma ocean

Io experiences tidal deformation as a result of its eccentric orbit around Jupiter, which provides a primary energy source for Io's continuing volcanic activity and infrared emission1. The amount of tidal energy dissipated within Io is enormous and has been suggested to support the large-scale melting of its interior and the formation of a global subsurface magma ocean. If Io has a shallow global magma ocean, its tidal deformation would be much larger than in the case of a more rigid, mostly solid interior2. Here we report the measurement of Io's tidal deformation, quantified by the gravitational tidal Love number k2, enabled by two recent flybys of the Juno spacecraft. By combining Juno3,4 and Galileo5-7 Doppler data from the NASA Deep Space Network and astrometric observations, we recover Re(k2) of 0.125 ± 0.047 (1σ) and the tidal dissipation parameter Q of 11.4 ± 3.6 (1σ). These measurements confirm that a shallow global magma ocean in Io does not exist and are consistent with Io having a mostly solid mantle2. Our results indicate that tidal forces do not universally create global magma oceans, which may be prevented from forming owing to rapid melt ascent, intrusion and eruption8,9, so even strong tidal heating-such as that expected on several known exoplanets and super-Earths10-may not guarantee the formation of magma oceans on moons or planetary bodies.

Unveiling magmatic structures and connectivity beneath the lunar Oceanus Procellarum region from GRAIL gravity data

The Oceanus Procellarum region, characterized by its vast basaltic plains and pronounced volcanic activity, serves as a focal point for understanding the volcanic history of the Moon. Here we present density models of the magmatic structures beneath Oceanus Procellarum, derived from Gravity Recovery and Interior Laboratory (GRAIL) mission data. The models uncover pronounced linear magmatic structures along the Procellarum's western border and significant intrusions within the northern and southern Marius Hills. Crucially, they reveal three narrow near-horizontal sheeted magmatic conduits, 80-150 km long, extending from near-surface to 6-7 km depth. These magmatic conduits connect the Marius Hills' northern and southern intrusions and bridge them with the Procellarum western border structures, suggesting that they likely served as central pathways facilitating magma transport across various volcanic systems. These discoveries reveal widespread magmatic connectivity beneath Oceanus Procellarum and underscore the critical role of lateral magma transport processes in shaping the Moon's volcanic evolution.

Geological context of the Chang'e-6 landing area and implications for sample analysis

Research on returned samples can provide ground truth for the study of the geological evolution history of the Moon. However, previous missions all collected samples from the near side of the Moon, which is significantly different from the far side of the Moon in terms of the thickness of the lunar crust, magma activity, and composition. Therefore, the samples from the far side of the Moon are of great significance for a comprehensive understanding of the history of the Moon. China's Chang'e-6 (CE-6) probe has successfully landed on the lunar far side and will return samples in the coming days. With the precise location of the CE-6 landing site, a detailed analysis of the geological background is conducted in this research. The landing site of CE-6 is within the Apollo crater, which is inside the largest impact basin on the Moon, i.e., the South Pole-Aitken (SPA) basin. According to the numerical simulation of the formation process of the SPA basin, CE-6 landed at the edge of the SPA impact melting zone, which is presumably composed of impact melt of the lunar mantle. The Apollo crater subsequently excavated deep material again, which constitutes the basement of the CE-6 landing area. Later, erupted basalt covered these basement rocks, and they also constitute the main source of the CE-6 samples. Based on the dating method of crater size-frequency distribution, we find that the basalt is ∼2.50 Ga. The CE-6 samples also possibly contain basement rocks as excavated and ejected by craters, and they can provide crucial information for our understanding of lunar geological history along with the basalt samples.

Evidence for structural control of mare volcanism in lunar compressional tectonic settings

One of the long-standing enigmas for lunar tectonic-thermal evolution is the spatiotemporal association of contractional wrinkle ridges and basaltic volcanism in a compressional regime. Here, we show that most of the 30 investigated volcanic (eruptive) centers are linked to contractional wrinkle ridges developed above preexisting basin basement-involved ring/rim normal faults. Based on the tectonic patterns associated with the basin formation and mass loading and considering that during the subsequent compression the stress was not purely isotropic, we hypothesize that tectonic inversion produced not only thrust faults but also reactivated structures with strike-slip and even extensional components, thus providing a valid mechanism for magma transport through fault planes during ridge faulting and folding of basaltic layers. Our findings suggest that lunar syn-tectonic mare emplacement along reactivated inherited faults provides important records of basin-scale structure-involved volcanism, which is more complex than previously considered.

Widespread impact-generated porosity in early planetary crusts

NASA’s Gravity Recovery and Interior Laboratory (GRAIL) spacecraft revealed the crust of the Moon is highly porous, with ~4% porosity at 20 km deep. The deep lying porosity discovered by GRAIL has been difficult to explain, with most current models only able to explain high porosity near the lunar surface (first few kilometers) or inside complex craters. Using hydrocode routines we simulated fracturing and generation of porosity by large impacts in lunar, martian, and Earth crust. Our simulations indicate impacts that produce 100–1000 km scale basins alone are capable of producing all observed porosity within the lunar crust. Simulations under the higher surface gravity of Mars and Earth suggest basin forming impacts can be a primary source of porosity and fracturing of ancient planetary crusts. Thus, we show that impacts could have supported widespread crustal fluid circulation, with important implications for subsurface habitable environments on early Earth and Mars.

Large impacts can create deep lying porosity far away from the crater. This result explains GRAIL’s findings and suggests impacts could support widespread fluid circulation, which has implications for habitable environments on early Earth and Mars.

The Identification of Impact Craters from GRAIL-Acquired Gravity Data by U-Net Architecture

The identification of impact craters on the Moon and other planetary bodies is of great significance to studying and constraining the dynamical process and evolution of the Solar System. Traditionally, this has been performed through the visual examination of images. Due to the effect of overburden, some structural features cannot be effectively identified from optical images, resulting in limitations in the scope, efficiency and accuracy of identification. In this paper, we investigate the viability of convolutional neural networks (CNNs) to perform the detection of impact craters from GRAIL-acquired gravity data. The ideal values of each hyperparameter in U-net architecture are determined after dozens of iterations of model training, testing and evaluation. The final model was evaluated by the Loss function with the low value of 0.04, indicating that the predicted output of the model reached a relatively high fitting degree with the prior labelled output. The comparative results with different methods show that the proposed method has a clear detection of the target features, with an accuracy of more than 80%. In addition, the detection results of the whole image account for 83% of the number of manually delineated gravity anomalies. The proposed method can still achieve the same quality for the identification of the gravity anomalies caused by impact craters under the condition that the resolution of GRAIL gravity data are not superior. Our results demonstrate that the U-net architecture can be a very effective tool for the rapid and automatic identification of impact craters from gravity map on the Moon, as well as other Solar System bodies.

Distinguishing the Origin of Asteroid (16) Psyche

The asteroid (16) Psyche may be the metal-rich remnant of a differentiated planetesimal, or it may be a highly reduced, metal-rich asteroidal material that never differentiated. The NASA Psyche mission aims to determine Psyche's provenance. Here we describe the possible solar system regions of origin for Psyche, prior to its likely implantation into the asteroid belt, the physical and chemical processes that can enrich metal in an asteroid, and possible meteoritic analogs. The spacecraft payload is designed to be able to discriminate among possible formation theories. The project will determine Psyche's origin and formation by measuring any strong remanent magnetic fields, which would imply it was the core of a differentiated body; the scale of metal to silicate mixing will be determined by both the neutron spectrometers and the filtered images; the degree of disruption between metal and rock may be determined by the correlation of gravity with composition; some mineralogy (e.g., modeled silicate/metal ratio, and inferred existence of low-calcium pyroxene or olivine, for example) will be detected using filtered images; and the nickel content of Psyche's metal phase will be measured using the GRNS.

Lunar Megaregolith Structure Revealed by GRAIL Gravity Data

We use gravity data from NASA's GRAIL mission to characterize the porosity structure of the upper lunar crust. We analyze the gravitational anomalies produced by the porosity of craters with diameters D between 10 and 30 km. We find that the gravitational signature of craters changes significantly atD = 16 . 4 - 0.6 + 1.4 km, which is related to a discrete change in porosity at a depth ∼3-5 km. We propose that this discrete porosity change reveals the location of the boundary between large-scale basin ejecta and the deeper less porous portion of the megaregolith, known as the structurally disturbed crust. The ejecta thickness can help constrain models of material transport and mixing on the Moon and, because the ejecta layer acts as an insulating blanket, models of heat flow and magmatism.

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 0.15 ° × 0.15 ° (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.

Large impact cratering during lunar magma ocean solidification

The lunar cratering record is used to constrain the bombardment history of both the Earth and the Moon. However, it is suggested from different perspectives, including impact crater dating, asteroid dynamics, lunar samples, impact basin-forming simulations, and lunar evolution modelling, that the Moon could be missing evidence of its earliest cratering record. Here we report that impact basins formed during the lunar magma ocean solidification should have produced different crater morphologies in comparison to later epochs. A low viscosity layer, mimicking a melt layer, between the crust and mantle could cause the entire impact basin size range to be susceptible to immediate and extreme crustal relaxation forming almost unidentifiable topographic and crustal thickness signatures. Lunar basins formed while the lunar magma ocean was still solidifying may escape detection, which is agreeing with studies that suggest a higher impact flux than previously thought in the earliest epoch of Earth-Moon evolution.

Observations, Meteorites, and Models: A Preflight Assessment of the Composition and Formation of (16) Psyche

Some years ago, the consensus was that asteroid (16) Psyche was almost entirely metal. New data on density, radar properties, and spectral signatures indicate that the asteroid is something perhaps even more enigmatic: a mixed metal and silicate world. Here we combine observations of Psyche with data from meteorites and models for planetesimal formation to produce the best current hypotheses for Psyche's properties and provenance. Psyche's bulk density appears to be between 3,400 and 4,100 kg m-3. Psyche is thus predicted to have between ~30 and ~60 vol% metal, with the remainder likely low-iron silicate rock and not more than ~20% porosity. Though their density is similar, mesosiderites are an unlikely analog to bulk Psyche because mesosiderites have far more iron-rich silicates than Psyche appears to have. CB chondrites match both Psyche's density and spectral properties, as can some pallasites, although typical pallasitic olivine contains too much iron to be consistent with the reflectance spectra. Final answers, as well as resolution of contradictions in the data set of Psyche physical properties, for example, the thermal inertia measurements, may not be resolved until the NASA Psyche mission arrives in orbit at the asteroid. Despite the range of compositions and formation processes for Psyche allowed by the current data, the science payload of the Psyche mission (magnetometers, multispectral imagers, neutron spectrometer, and a gamma-ray spectrometer) will produce data sets that distinguish among the models.

Improving the geometry of Kaguya extended mission data through refined orbit determination using laser altimetry

The Japan Aerospace Exploration Agency's (JAXA) Kaguya spacecraft carried a suite of instruments to map the Moon and its environment globally. During its extended mission, the average altitude was 50 km or lower, and Kaguya science products using these data hence have an increased spatial resolution. However, the geodetic position quality of these products is much worse than that of those acquired during the primary mission (at an altitude of 100 km) because of reduced radiometric tracking and frequent thrusting to maintain spacecraft attitude after the loss of momentum wheels. We have analyzed the Kaguya tracking data using gravity models based on the Gravity Recovery and Interior Laboratory (GRAIL) mission, and by making use of a new data type based on laser altimeter data collected by Kaguya: we adjust the spacecraft orbit such that the altimetry tracks fit a precise topographic basemap based on the Lunar Reconnaissance Orbiter's (LRO) Lunar Orbiter Laser Altimeter (LOLA) data. This results in geodetically accurate orbits tied to the precise LOLA/LRO frame. Whereas previously archived orbits show errors at the level of several a level of several tens of meters. When altimetry data are not available, the combination of GRAIL gravity and radio tracking results in an orbit precision of around several hundreds of meters for the low-altitude phase of the extended mission. Our greatly improved orbits result in better geolocation of the Kaguya extended mission data set.

Impact Melt Facies in the Moon's Crisium Basin: Identifying, Characterizing, and Future Radiogenic Dating

Both Earth and the Moon share a common history regarding the epoch of large basin formation, though only the lunar geologic record preserves any appreciable record of this Late Heavy Bombardment. The emergence of Earth's first life is approximately contemporaneous with the Late Heavy Bombardment; understanding the latter informs the environmental conditions of the former, which are likely necessary to constrain the mechanisms of abiogenesis. While the relative formation time of most of the Moon's large basins is known, the absolute timing is not. The timing of Crisium Basin's formation is one of many important events that must be constrained and would require identifying and dating impact melt formed in the Crisium event. To inform a future lunar sample dating mission, we thus characterized possible outcrops of impact melt. We determined that several mare lava-embayed kipukas could contain impact melt, though the rim and central peaks of the partially lava-flooded Yerkes Crater likely contain the most pure and intact Crisium impact melt. It is here where future robotic and/or human missions could confidently add a key missing piece to the puzzle of the combined issues of early Earth-Moon bombardment and the emergence of life.

Density Structure of the Von Kármán Crater in the Northwestern South Pole-Aitken Basin: Initial Subsurface Interpretation of the Chang'E-4 Landing Site Region

The Von Kármán Crater, within the South Pole-Aitken (SPA) Basin, is the landing site of China's Chang'E-4 mission. To complement the in situ exploration mission and provide initial subsurface interpretation, we applied a 3D density inversion using the Gravity Recovery and Interior Laboratory (GRAIL) gravity data. We constrain our inversion method using known geological and geophysical lunar parameters to reduce the non-uniqueness associated with gravity inversion. The 3D density models reveal vertical and lateral density variations, 2600-3200 kg/m3, assigned to the changing porosity beneath the Von Kármán Crater. We also identify two mass excess anomalies in the crust with a steep density contrast of 150 kg/m3, which were suggested to have been caused by multiple impact cratering. The anomalies from recovered near surface density models, together with the gravity derivative maps extending to the lower crust, are consistent with surface geological manifestation of excavated mantle materials from remote sensing studies. Therefore, we suggest that the density distribution of the Von Kármán Crater indicates multiple episodes of impact cratering that resulted in formation and destruction of ancient craters, with crustal reworking and excavation of mantle materials.

GRAIL-identified gravity anomalies in Oceanus Procellarum: Insight into subsurface impact and magmatic structures on the Moon

Four, quasi-circular, positive Bouguer gravity anomalies (PBGAs) that are similar in diameter (~90-190 km) and gravitational amplitude (>140 mGal contrast) are identified within the central Oceanus Procellarum region of the Moon. These spatially associated PBGAs are located south of Aristarchus Plateau, north of Flamsteed crater, and two are within the Marius Hills volcanic complex (north and south). Each is characterized by distinct surface geologic features suggestive of ancient impact craters and/or volcanic/plutonic activity. Here, we combine geologic analyses with forward modeling of high-resolution gravity data from the Gravity Recovery and Interior Laboratory (GRAIL) mission in order to constrain the subsurface structures that contribute to these four PBGAs. The GRAIL data presented here, at spherical harmonic degrees 6-660, permit higher resolution analyses of these anomalies than previously reported, and reveal new information about subsurface structures. Specifically, we find that the amplitudes of the four PBGAs cannot be explained solely by mare-flooded craters, as suggested in previous work; an additional density contrast is required to explain the high-amplitude of the PBGAs. For Northern Flamsteed (190 km diameter), the additional density contrast may be provided by impact-related mantle uplift. If the local crust has a density ~2800 kg/m³, then ~7 km of uplift is required for this anomaly, although less uplift is required if the local crust has a lower mean density of ~2500 kg/m³. For the Northern and Southern Marius Hills anomalies, the additional density contrast is consistent with the presence of a crustal complex of vertical dikes that occupies up to ~37% of the regionally thin crust. The structure of Southern Aristarchus Plateau (90 km diameter), an anomaly with crater-related topographic structures, remains ambiguous. Based on the relatively small size of the anomaly, we do not favor mantle uplift, however understanding mantle response in a region of especially thin crust needs to be better resolved. It is more likely that this anomaly is due to subsurface magmatic material given the abundance of volcanic material in the surrounding region. Overall, the four PBGAs analyzed here are important in understanding the impact and volcanic/plutonic history of the Moon, specifically in a region of thin crust and elevated temperatures characteristic of the Procellarum KREEP Terrane.

History of the Terminal Cataclysm Paradigm: Epistemology of a Planetary Bombardment That Never (?) Happened

This study examines the history of the paradigm concerning a lunar (or solar-systemwide)terminal cataclysm (also called “Late Heavy Bombardment” or LHB), a putative, brief spikein impacts at ~3.9 Ga ago, preceded by low impact rates. We examine origin of the ideas, why theywere accepted, and why the ideas are currently being seriously revised, if not abandoned. Thepaper is divided into the following sections:1. Overview of paradigm.2. Pre-Apollo views (1949-1969).3. Initial suggestions of cataclysm (ca. 1974).4. Ironies.5. Alternative suggestions, megaregolith evolution (1970s).6. Impact melt rocks “establish” cataclysm (1990).7. Imbrium redux (ca. 1998).8. Impact melt clasts (early 2000s).9. Dating of front-side lunar basins?10. Dynamical models “explain” the cataclysm (c. 2000s).11. Asteroids as a test case.12. Impact melts predating 4.0 Ga ago (ca. 2008-present.).13. Biological issues.14. Growing doubts (ca. 1994-2014).15. Evolving Dynamical Models (ca. 2001-present).16. Connections to lunar origin.17. Dismantling the paradigm (2015-2018).18. “Megaregolith Evolution Model” for explaining the data.19. Conclusions and new directions for future work.The author hopes that this open-access discussion may prove useful for classroom discussionsof how science moves forward through self-correction of hypotheses.

A surface gravity traverse on Mars indicates low bedrock density at Gale crater

Gravimetry, the precise measurement of gravitational fields, can be used to probe the internal structure of Earth and other planets. The Curiosity rover on Mars carries accelerometers normally used for navigation and attitude determination. We have recalibrated them to isolate the signature of the changing gravitational acceleration as the rover climbs through Gale crater. The subsurface rock density is inferred from the measured decrease in gravitational field strength with elevation. The density of the sedimentary rocks in Gale crater is 1680 ± 180 kilograms per cubic meter. This value is lower than expected, indicating a high porosity and constraining maximum burial depths of the rocks over their history.

Orbit determination of the Lunar Reconnaissance Orbiter: Status after Seven Years

The Lunar Reconnaissance Orbiter (LRO) has been orbiting the Moon since 2009, obtaining unique and foundational datasets important to understanding the evolution of the Moon and the Solar System. The high-resolution data acquired by LRO benefit from precise orbit determination (OD), limiting the need for geolocation and co-registration tasks. The initial position knowledge requirement (50 meters) was met with radio tracking from ground stations, after combination with LOLA altimetric crossovers. LRO-specific gravity field solutions were determined and allowed radio-only OD to perform at the level of 20 meters, although secular inclination changes required frequent updates. The high-accuracy gravity fields from GRAIL, with <10 km spatial resolution, further improved the radio-only orbit reconstruction quality (<10 meters). However, orbit reconstruction is in part limited by the 0.3-0.5 mm/s measurement noise level in S-band tracking. One-way tracking through Laser Ranging can supplement the tracking available for OD with 28-Hz ranges with 20-cm single-shot precision, but is available only on the nearside (the lunar hemisphere facing the Earth due to tidal locking). Here, we report on the status of the OD effort since the beginning of the mission, a period spanning more than seven years. We describe modeling improvements and the use of new measurements. In particular, the LOLA altimetric data give accurate, uniform, and independent information about LRO's orbit, with a different sensitivity and geometry which includes coverage over the lunar farside and is not tied to ground-based assets. With SLDEM2015 (a combination of the LOLA topographic profiles and the Kaguya Terrain Camera stereo images), another use of altimetry is possible for OD. We extend the 'direct altimetry' technique developed for the ICESat mission to perform OD and adjust spacecraft position to minimize discrepancies between LOLA tracks and SLDEM2015. Comparisons with the radio-only orbits are used to evaluate this new tracking type, of interest for the OD of future lunar orbiters carrying a laser altimeter. LROC NAC images also provide independent accuracy estimation, through the repeated views taken of anthropogenic features for instance.

In-flight characterization of the lunar orbiter laser altimeter instrument pointing and far-field pattern

The Lunar Orbiter Laser Altimeter (LOLA) aboard the Lunar Reconnaissance Orbiter (LRO) has collected nearly seven billion measurements of surface height on the Moon with an absolute accuracy of ∼1  m and a precision of ∼10  cm. Converting time-of-flight laser altimeter measurements to topographic elevations requires accurate knowledge of the laser pointing with respect to the spacecraft body-fixed coordinate system. To that end, we have utilized altimetric crossovers from LOLA, as well as bidirectional observations of the LOLA laser and receiver boresight via an Earth-based laser tracking ground station. Based on a sample of ∼780,000 globally distributed crossovers from the circular-orbit phase of LRO's mission (∼27 months), we derive corrections to the LOLA laser boresight. These corrections improve the cross-track and along-track agreement of the crossovers by 24% and 33%, respectively, yielding RMS residuals of ∼10  m. Since early in the LRO mission, the bidirectional laser tracking experiments have confirmed a pointing anomaly when the LOLA instrument is facing toward deep space or the night side of the Moon and have allowed the reconstruction of the laser far-field pattern and receiver telescope pointing. By conducting such experiments shortly after launch and nearly eight years later, we have directly measured changes in the laser characteristics and obtained critical data to understand the laser behavior and refine the instrument pointing model. The methods and results presented here are also relevant to the design, fabrication, and operation of future planetary laser altimeters and their long-term behavior in the space environment.

Ring faults and ring dikes around the Orientale basin on the Moon

The Orientale basin is the youngest and best-preserved multiring impact basin on the Moon, having experienced only modest modification by subsequent impacts and volcanism. Orientale is often treated as the type example of a multiring basin, with three prominent rings outside of the inner depression: the Inner Rook Montes, the Outer Rook Montes, and the Cordillera. Here we use gravity data from NASA's Gravity Recovery and Interior Laboratory (GRAIL) mission to reveal the subsurface structure of Orientale and its ring system. Gradients of the gravity data reveal a continuous ring dike intruded into the Outer Rook along the plane of the fault associated with the ring scarp. The volume of this ring dike is ~18 times greater than the volume of all extrusive mare deposits associated with the basin. The gravity gradient signature of the Cordillera ring indicates an offset along the fault across a shallow density interface, interpreted to be the base of the low-density ejecta blanket. Both gravity gradients and crustal thickness models indicate that the edge of the central cavity is shifted inward relative to the equivalent Inner Rook ring at the surface. Models of the deep basin structure show inflections along the crust-mantle interface at both the Outer Rook and Cordillera rings, indicating that the basin ring faults extend from the surface to at least the base of the crust. Fault dips range from 13-22° for the Cordillera fault in the northeastern quadrant, to 90° for the Outer Rook in the northwestern quadrant. The fault dips for both outer rings are lowest in the northeast, possibly due to the effects of either the direction of projectile motion or regional gradients in pre-impact crustal thickness. Similar ring dikes and ring faults are observed around the majority of lunar basins.

Research and Development of Electrostatic Accelerometers for Space Science Missions at HUST

High-precision electrostatic accelerometers have achieved remarkable success in satellite Earth gravity field recovery missions. Ultralow-noise inertial sensors play important roles in space gravitational wave detection missions such as the Laser Interferometer Space Antenna (LISA) mission, and key technologies have been verified in the LISA Pathfinder mission. Meanwhile, at Huazhong University of Science and Technology (HUST, China), a space accelerometer and inertial sensor based on capacitive sensors and the electrostatic control technique have also been studied and developed independently for more than 16 years. In this paper, we review the operational principle, application, and requirements of the electrostatic accelerometer and inertial sensor in different space missions. The development and progress of a space electrostatic accelerometer at HUST, including ground investigation and space verification are presented.

Evidence for a Low Bulk Crustal Density for Mars from Gravity and Topography

Knowledge of the average density of the crust of a planet is important in determining its interior structure. The combination of high-resolution gravity and topography data has yielded a low density for the Moon's crust, yet for other terrestrial planets the resolution of the gravity field models has hampered reasonable estimates. By using well-chosen constraints derived from topography during gravity field model determination using satellite tracking data, we show that we can robustly and independently determine the average bulk crustal density directly from the tracking data, using the admittance between topography and imperfect gravity. We find a low average bulk crustal density for Mars, 2582 ± 209 kg m^-3. This bulk crustal density is lower than that assumed until now. Densities for volcanic complexes are higher, consistent with earlier estimates, implying large lateral variations in crustal density. In addition, we find indications that the crustal density increases with depth.

Gravity field of the Orientale basin from the Gravity Recovery and Interior Laboratory Mission

The Orientale basin is the youngest and best-preserved major impact structure on the Moon. We used the Gravity Recovery and Interior Laboratory (GRAIL) spacecraft to investigate the gravitational field of Orientale at 3- to 5-kilometer (km) horizontal resolution. A volume of at least (3.4 ± 0.2) × 106 km3 of crustal material was removed and redistributed during basin formation. There is no preserved evidence of the transient crater that would reveal the basin's maximum volume, but its diameter may now be inferred to be between 320 and 460 km. The gravity field resolves distinctive structures of Orientale's three rings and suggests the presence of faults associated with the outer two that penetrate to the mantle. The crustal structure of Orientale provides constraints on the formation of multiring basins.

Formation of the Orientale lunar multiring basin

Multiring basins, large impact craters characterized by multiple concentric topographic rings, dominate the stratigraphy, tectonics, and crustal structure of the Moon. Using a hydrocode, we simulated the formation of the Orientale multiring basin, producing a subsurface structure consistent with high-resolution gravity data from the Gravity Recovery and Interior Laboratory (GRAIL) spacecraft. The simulated impact produced a transient crater, ~390 kilometers in diameter, that was not maintained because of subsequent gravitational collapse. Our simulations indicate that the flow of warm weak material at depth was crucial to the formation of the basin's outer rings, which are large normal faults that formed at different times during the collapse stage. The key parameters controlling ring location and spacing are impactor diameter and lunar thermal gradients.

Small-scale density variations in the lunar crust revealed by GRAIL

Data from the Gravity Recovery and Interior Laboratory (GRAIL) mission have revealed that ~98% of the power of the gravity signal of the Moon at high spherical harmonic degrees correlates with the topography. The remaining 2% of the signal, which cannot be explained by topography, contains information about density variations within the crust. These high-degree Bouguer gravity anomalies are likely caused by small-scale (10's of km) shallow density variations. Here we use gravity inversions to model the small-scale three-dimensional variations in the density of the lunar crust. Inversion results from three non-descript areas yield shallow density variations in the range of 100-200 kg/m3. Three end-member scenarios of variations in porosity, intrusions into the crust, and variations in bulk crustal composition were tested as possible sources of the density variations. We find that the density anomalies can be caused entirely by changes in porosity. Characteristics of density anomalies in the South Pole-Aitken basin also support porosity as a primary source of these variations. Mafic intrusions into the crust could explain many, but not all of the anomalies. Additionally, variations in crustal composition revealed by spectral data could only explain a small fraction of the density anomalies. Nevertheless, all three sources of density variations likely contribute. Collectively, results from this study of GRAIL gravity data, combined with other studies of remote sensing data and lunar samples, show that the lunar crust exhibits variations in density by ±10% over scales ranging from centimeters to 100's of kilometers.

Asteroid bombardment and the core of Theia as possible sources for the Earth's late veneer component

The silicate Earth contains Pt-group elements in roughly chondritic relative ratios, but with absolute concentrations <1% chondrite. This veneer implies addition of chondrite-like material with 0.3-0.7% mass of the Earth's mantle or an equivalent planet-wide thickness of 5-20 km. The veneer thickness, 200-300 m, within the lunar crust and mantle is much less. One hypothesis is that the terrestrial veneer arrived after the moon-forming impact within a few large asteroids that happened to miss the smaller Moon. Alternatively, most of terrestrial veneer came from the core of the moon-forming impactor, Theia. The Moon then likely contains iron from Theia's core. Mass balances lend plausibility. The lunar core mass is ~1.6 × 10²¹ kg and the excess FeO component in the lunar mantle is 1.3-3.5 × 10²¹ kg as Fe, totaling 3-5 × 10²¹ kg or a few percent of Theia's core. This mass is comparable to the excess Fe of 2.3-10 × 10²¹ kg in the Earth's mantle inferred from the veneer component. Chemically in this hypothesis, Fe metal from Theia's core entered the Moon-forming disk. H₂O and Fe₂O₃ in the disk oxidized part of the Fe, leaving the lunar mantle near a Fe-FeO buffer. The remaining iron metal condensed, gathered Pt-group elements eventually into the lunar core. The silicate Moon is strongly depleted in Pt-group elements. In contrast, the Earth's mantle contained excess oxidants, H₂O and Fe₂O₃, which quantitatively oxidized the admixed Fe from Theia's core, retaining Pt-group elements. In this hypothesis, asteroid impacts were relatively benign with ~1 terrestrial event that left only thermophile survivors.

Lunar impact basins revealed by Gravity Recovery and Interior Laboratory measurements

Observations from the Gravity Recovery and Interior Laboratory (GRAIL) mission indicate a marked change in the gravitational signature of lunar impact structures at the morphological transition, with increasing diameter, from complex craters to peak-ring basins. At crater diameters larger than ~200 km, a central positive Bouguer anomaly is seen within the innermost peak ring, and an annular negative Bouguer anomaly extends outward from this ring to the outer topographic rim crest. These observations demonstrate that basin-forming impacts remove crustal materials from within the peak ring and thicken the crust between the peak ring and the outer rim crest. A correlation between the diameter of the central Bouguer gravity high and the outer topographic ring diameter for well-preserved basins enables the identification and characterization of basins for which topographic signatures have been obscured by superposed cratering and volcanism. The GRAIL inventory of lunar basins improves upon earlier lists that differed in their totals by more than a factor of 2. The size-frequency distributions of basins on the nearside and farside hemispheres of the Moon differ substantially; the nearside hosts more basins larger than 350 km in diameter, whereas the farside has more smaller basins. Hemispherical differences in target properties, including temperature and porosity, are likely to have contributed to these different distributions. Better understanding of the factors that control basin size will help to constrain models of the original impactor population.

Structure and evolution of the lunar Procellarum region as revealed by GRAIL gravity data

The Procellarum region is a broad area on the nearside of the Moon that is characterized by low elevations, thin crust, and high surface concentrations of the heat-producing elements uranium, thorium, and potassium. The region has been interpreted as an ancient impact basin approximately 3,200 kilometres in diameter, although supporting evidence at the surface would have been largely obscured as a result of the great antiquity and poor preservation of any diagnostic features. Here we use data from the Gravity Recovery and Interior Laboratory (GRAIL) mission to examine the subsurface structure of Procellarum. The Bouguer gravity anomalies and gravity gradients reveal a pattern of narrow linear anomalies that border Procellarum and are interpreted to be the frozen remnants of lava-filled rifts and the underlying feeder dykes that served as the magma plumbing system for much of the nearside mare volcanism. The discontinuous surface structures that were earlier interpreted as remnants of an impact basin rim are shown in GRAIL data to be a part of this continuous set of border structures in a quasi-rectangular pattern with angular intersections, contrary to the expected circular or elliptical shape of an impact basin. The spatial pattern of magmatic-tectonic structures bounding Procellarum is consistent with their formation in response to thermal stresses produced by the differential cooling of the province relative to its surroundings, coupled with magmatic activity driven by the greater-than-average heat flux in the region.

The tidal-rotational shape of the Moon and evidence for polar wander

The origin of the Moon's large-scale topography is important for understanding lunar geology, lunar orbital evolution and the Moon's orientation in the sky. Previous hypotheses for its origin have included late accretion events, large impacts, tidal effects and convection processes. However, testing these hypotheses and quantifying the Moon's topography is complicated by the large basins that have formed since the crust crystallized. Here we estimate the large-scale lunar topography and gravity spherical harmonics outside these basins and show that the bulk of the spherical harmonic degree-2 topography is consistent with a crust-building process controlled by early tidal heating throughout the Moon. The remainder of the degree-2 topography is consistent with a frozen tidal-rotational bulge that formed later, at a semi-major axis of about 32 Earth radii. The probability of the degree-2 shape having both tidal-heating and frozen shape characteristics by chance is less than 1%. We also infer that internal density contrasts eventually reoriented the Moon's polar axis by 36 ± 4°, to the configuration we observe today. Together, these results link the geology of the near and far sides, and resolve long-standing questions about the Moon's large-scale shape, gravity and history of polar wander.

GRGM900C: A degree 900 lunar gravity model from GRAIL primary and extended mission data

We have derived a gravity field solution in spherical harmonics to degree and order 900, GRGM900C, from the tracking data of the Gravity Recovery and Interior Laboratory (GRAIL) Primary (1 March to 29 May 2012) and Extended Missions (30 August to 14 December 2012). A power law constraint of 3.6 ×10-4/ℓ2 was applied only for degree ℓ greater than 600. The model produces global correlations of gravity, and gravity predicted from lunar topography of ≥ 0.98 through degree 638. The model's degree strength varies from a minimum of 575-675 over the central nearside and farside to 900 over the polar regions. The model fits the Extended Mission Ka-Band Range Rate data through 17 November 2012 at 0.13 μm/s RMS, whereas the last month of Ka-Band Range-Rate data obtained from altitudes of 2-10 km fit at 0.98 μm/s RMS, indicating that there is still signal inherent in the tracking data beyond degree 900.

High-resolution local gravity model of the south pole of the Moon from GRAIL extended mission data

We estimated a high-resolution local gravity field model over the south pole of the Moon using data from the Gravity Recovery and Interior Laboratory's extended mission. Our solution consists of adjustments with respect to a global model expressed in spherical harmonics. The adjustments are expressed as gridded gravity anomalies with a resolution of 1/6° by 1/6° (equivalent to that of a degree and order 1080 model in spherical harmonics), covering a cap over the south pole with a radius of 40°. The gravity anomalies have been estimated from a short-arc analysis using only Ka-band range-rate (KBRR) data over the area of interest. We apply a neighbor-smoothing constraint to our solution. Our local model removes striping present in the global model; it reduces the misfit to the KBRR data and improves correlations with topography to higher degrees than current global models.

KEY POINTS

We present a high-resolution gravity model of the south pole of the Moon Improved correlations with topography to higher degrees than global models Improved fits to the data and reduced striping that is present in global models.

Detection of the lunar body tide by the Lunar Orbiter Laser Altimeter

The Lunar Orbiter Laser Altimeter instrument onboard the Lunar Reconnaissance Orbiter spacecraft collected more than 5 billion measurements in the nominal 50 km orbit over ∼10,000 orbits. The data precision, geodetic accuracy, and spatial distribution enable two-dimensional crossovers to be used to infer relative radial position corrections between tracks to better than ∼1 m. We use nearly 500,000 altimetric crossovers to separate remaining high-frequency spacecraft trajectory errors from the periodic radial surface tidal deformation. The unusual sampling of the lunar body tide from polar lunar orbit limits the size of the typical differential signal expected at ground track intersections to ∼10 cm. Nevertheless, we reliably detect the topographic tidal signal and estimate the associated Love number h2 to be 0.0371 ± 0.0033, which is consistent with but lower than recent results from lunar laser ranging.

KEY POINTS

Altimetric data are used to create radial constraints on the tidal deformationThe body tide amplitude is estimated from the crossover dataThe estimated Love number is consistent with previous estimates but more precise.

The origin of lunar mascon basins

High-resolution gravity data from the Gravity Recovery and Interior Laboratory spacecraft have clarified the origin of lunar mass concentrations (mascons). Free-air gravity anomalies over lunar impact basins display bull's-eye patterns consisting of a central positive (mascon) anomaly, a surrounding negative collar, and a positive outer annulus. We show that this pattern results from impact basin excavation and collapse followed by isostatic adjustment and cooling and contraction of a voluminous melt pool. We used a hydrocode to simulate the impact and a self-consistent finite-element model to simulate the subsequent viscoelastic relaxation and cooling. The primary parameters controlling the modeled gravity signatures of mascon basins are the impactor energy, the lunar thermal gradient at the time of impact, the crustal thickness, and the extent of volcanic fill.

The crust of the Moon as seen by GRAIL

High-resolution gravity data obtained from the dual Gravity Recovery and Interior Laboratory (GRAIL) spacecraft show that the bulk density of the Moon's highlands crust is 2550 kilograms per cubic meter, substantially lower than generally assumed. When combined with remote sensing and sample data, this density implies an average crustal porosity of 12% to depths of at least a few kilometers. Lateral variations in crustal porosity correlate with the largest impact basins, whereas lateral variations in crustal density correlate with crustal composition. The low-bulk crustal density allows construction of a global crustal thickness model that satisfies the Apollo seismic constraints, and with an average crustal thickness between 34 and 43 kilometers, the bulk refractory element composition of the Moon is not required to be enriched with respect to that of Earth.

Ancient igneous intrusions and early expansion of the Moon revealed by GRAIL gravity gradiometry

The earliest history of the Moon is poorly preserved in the surface geologic record due to the high flux of impactors, but aspects of that history may be preserved in subsurface structures. Application of gravity gradiometry to observations by the Gravity Recovery and Interior Laboratory (GRAIL) mission results in the identification of a population of linear gravity anomalies with lengths of hundreds of kilometers. Inversion of the gravity anomalies indicates elongated positive-density anomalies that are interpreted to be ancient vertical tabular intrusions or dikes formed by magmatism in combination with extension of the lithosphere. Crosscutting relationships support a pre-Nectarian to Nectarian age, preceding the end of the heavy bombardment of the Moon. The distribution, orientation, and dimensions of the intrusions indicate a globally isotropic extensional stress state arising from an increase in the Moon's radius by 0.6 to 4.9 kilometers early in lunar history, consistent with predictions of thermal models.