Walter S. Kiefer

The Crust of the Moon as Seen by GRAIL

The Holy GRAIL? The gravity field of a planet provides a view of its interior and thermal history by revealing areas of different density. GRAIL, a pair of satellites that act as a highly sensitive gravimeter, began mapping the Moons gravity in early 2012. Three papers highlight some of the results from the primary mission. Zuber et al. (p. 668, published online 6 December) discuss the overall gravity field, which reveals several new tectonic and geologic features of the Moon. Impacts have worked to homogenize the density structure of the Moons upper crust while fracturing it extensively. Wieczorek et al. (p. 671, published online 6 December) show that the upper crust is 35 to 40 kilometers thick and less dense—and thus more porous—than previously thought. Finally, Andrews-Hanna et al. (p. 675, published online 6 December) show that the crust is cut by widespread magmatic dikes that may reflect a period of expansion early in the Moons history. The Moons gravity field shows that the lunar crust is less dense and more porous than was thought. High-resolution gravity data obtained from the dual Gravity Recovery and Interior Laboratory (GRAIL) spacecraft show that the bulk density of the Moons 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 Holy GRAIL? The gravity field of a planet provides a view of its interior and thermal history by revealing areas of different density. GRAIL, a pair of satellites that act as a highly sensitive gravimeter, began mapping the Moons gravity in early 2012. Three papers highlight some of the results from the primary mission. Zuber et al. (p. 668, published online 6 December) discuss the overall gravity field, which reveals several new tectonic and geologic features of the Moon. Impacts have worked to homogenize the density structure of the Moons upper crust while fracturing it extensively. Wieczorek et al. (p. 671, published online 6 December) show that the upper crust is 35 to 40 kilometers thick and less dense—and thus more porous—than previously thought. Finally, Andrews-Hanna et al. (p. 675, published online 6 December) show that the crust is cut by widespread magmatic dikes that may reflect a period of expansion early in the Moons history. The Moons gravity map shows that the crust is cut by extensive magmatic dikes, perhaps implying a period of early expansion. 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 Moons radius by 0.6 to 4.9 kilometers early in lunar history, consistent with predictions of thermal models.

A mantle plume model for the equatorial highlands of Venus

The Equatorial Highlands of Venus consist of four main structures, Atla, Beta, Ovda, and Thetis regiones. Each has a circular to oval-shaped planform and rises 4–5 km above the mean planetary radius. These highlands are associated with long-wavelength geoid highs, with amplitudes ranging from 35 m at Ovda to 120 m at Atla. They also contain topographic valleys, interpreted as extensional rift zones, and Beta is known to contain shield volcanoes. These characteristics are all consistent with the Equatorial Highlands being formed by mantle plumes. An alternative model, in which Ovda and Thetis are interpreted as spreading centers analogous to terrestrial mid-ocean ridges, fails to explain most of the observed geoid anomalies and topography in these regions. Some smaller highlands, such as Bell Regio, Eistla Regio, and the Hathor/Innini/Ushas region, may also be plume related, but most coronae are unlikely to be the direct result of plume activity. We have modeled plumes using a cylindrical, axisymmetric finite element code and a depth-dependent, Newtonian rheology. We compare our model results with profiles of geoid and topography across Atla, Beta, Ovda, and Thetis; our best model fits are for Beta and Atla. Assuming whole mantle convection and that Earth and Venus have similar mantle heat flows, Venus must lack an Earth-like low-viscosity zone in its upper mantle in order satisfy the observed geoid and topography for these features. This conclusion is consistent with the long-wavelength admittance spectrum of Venus and with the observed differences in the slopes of the geoid spectra for the two planets. One explanation for the different viscosity structures of the two planets could be that the mantle of Venus is drier than Earths mantle.

Lunar interior properties from the GRAIL mission

The Gravity Recovery and Interior Laboratory (GRAIL) mission has sampled lunar gravity with unprecedented accuracy and resolution. The lunar GM, the product of the gravitational constant G and the mass M, is very well determined. However, uncertainties in the mass and mean density, 3345.56 ± 0.40 kg/m3, are limited by the accuracy of G. Values of the spherical harmonic degree-2 gravity coefficients J2 and C22, as well as the Love number k2 describing lunar degree-2 elastic response to tidal forces, come from two independent analyses of the 3 month GRAIL Primary Mission data at the Jet Propulsion Laboratory and the Goddard Space Flight Center. The two k2 determinations, with uncertainties of ~1%, differ by 1%; the average value is 0.02416 ± 0.00022 at a 1 month period with reference radius R = 1738 km. Lunar laser ranging (LLR) data analysis determines (C − A)/B and (B − A)/C, where A < B < C are the principal moments of inertia; the flattening of the fluid outer core; the dissipation at its solid boundaries; and the monthly tidal dissipation Q = 37.5 ± 4. The moment of inertia computation combines the GRAIL-determined J2 and C22 with LLR-derived (C − A)/B and (B − A)/C. The normalized mean moment of inertia of the solid Moon is Is/MR2 = 0.392728 ± 0.000012. Matching the density, moment, and Love number, calculated models have a fluid outer core with radius of 200–380 km, a solid inner core with radius of 0–280 km and mass fraction of 0–1%, and a deep mantle zone of low seismic shear velocity. The mass fraction of the combined inner and outer core is ≤1.5%.

GRAIL gravity constraints on the vertical and lateral density structure of the lunar crust

We analyzed data from the Gravity Recovery and Interior Laboratory (GRAIL) mission using a localized admittance approach to map out spatial variations in the vertical density structure of the lunar crust. Mare regions are characterized by a distinct decrease in density with depth, while the farside is characterized by an increase in density with depth at an average gradient of ∼35 kg m −3 km −1 and typical surface porosities of at least 20%. The Apollo 12 and 14 landing site region has a similar density structure to the farside, permitting a comparison with seismic velocity profiles. The interior of the South Pole-Aitken (SP-A) impact basin appears distinct with a near-surface low-density (porous) layer 2-3 times thinner than the rest of the farside. This result suggests that redistribution of material during the large SP-A impact likely played a major role in sculpting the lunar crust.

Lunar impact basins revealed by Gravity Recovery and Interior Laboratory measurements

New gravity measurements greatly improve the Moon’s preserved impact basin inventory. 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.

An inversion of gravity and topography for mantle and crustal structure on Mars

Analysis of the gravity and topography of Mars presently provides our primary quantitative constraints on the internal structure of Mars. We present an inversion of the long-wavelength (harmonic degree ≤ 10) gravity and topography of Mars for lateral variations of mantle temperature and crustal thickness. Our formulation incorporates both viscous mantle flow (which most prior studies have neglected) and isostatically compensated density anomalies in the crust and lithosphere. Our nominal model has a 150-km-thick high-viscosity surface layer over an isoviscous mantle, with a core radius of 1840 km. It predicts lateral temperature variations of up to a few hundred degrees Kelvin relative to the mean mantle temperature, with high temperature under Tharsis and to a lesser extent under Elysium and cool temperatures elsewhere. Surprisingly, the model predicts crustal thinning beneath Tharsis. If correct, this implies that thinning of the crust by mantle shear stresses dominates over thickening of the crust by volcanism. The major impact basins (Hellas, Argyre, Isidis, Chryse, and Utopia) are regions of crustal thinning, as expected. Utopia is also predicted to be a region of hot mantle, which is hard to reconcile with the surface geology. An alternative model for Utopia treats it as a mascon basin. The Utopia gravity anomaly is consistent with the presence of a 1.2 to 1.6 km thick layer of uncompensated basalt, in good agreement with geologic arguments about the amount of volcanic fill in this area. The mantle thermal structure is the dominant contributor to the observed geoid in our inversion. The mantle also dominates the topography at the longest wavelengths, but shorter wavelengths (harmonic degrees ≥4) are dominated by the crustal structure. Because of the uncertainty about the appropriate numerical values for some of the models input parameters, we have examined the sensitivity of the model results to the planetary structural model (core radius and core and mantle densities), the mantles viscosity stratification, and the mean crustal thickness. The model results are insensitive to the specific thickness or viscosity contrast of the high-viscosity surface layer and to the mean crustal thickness in the range 25 to 100 km. Models with a large core radius or with an upper mantle low-viscosity zone require implausibly large lateral variations in mantle temperature.

Mantle downwelling and crustal convergence: A model for Ishtar Terra, Venus

The Ishtar Terra region contains the highest topography known on Venus, over 10 km above mean planetary radius, as well as abundant tectonic features, many of apparently compressional origin. These characteristics suggest that Ishtar is a crustal convergence zone overlying downwelling mantle. In order to explore quantitatively the implications of this hypothesis for Ishtars origin, we present models of viscous crustal flow driven by gradients in topography. Assuming a free-slip surface boundary condition, we find that if the crustal convergence hypothesis is correct, then the crustal thickness in the plains surrounding Ishtar can be no more than about 25 km thick. This upper bound assumes a cold (10 K km−1) geotherm and the stiffest available diabase flow law. If the geothermal gradient is larger or the rheology is weaker, the crust must be even thinner for net crustal convergence to be possible. This upper bound is in good agreement with several independent estimates of crustal thickness of 15–30 km in the plains of Venus based on modeling of the spacing of tectonic features and of impact crater relaxation. If the surface layer of Venus provides a no-slip boundary, then our models allow the crustal thickness in the plains to be up to 50 km, but the likely existence of faults that cut through the crust makes a no-slip surface layer unlikely. Our upper bound on crustal thickness is much less than that derived from an Airy isostasy model of Ishtars gravity anomaly. Much of the observed gravity anomaly must be due to density anomalies in the mantle beneath Ishtar. Although we treat Ishtar as a crustal convergence zone, our crustal flow model shows that under some circumstances near-surface material may actually flow away from Ishtar, providing a possible explanation for grabenlike structures in Fortuna Tessera.

High pressure, near‐liquidus phase equilibria of the Home Plate basalt Fastball and melting in the Martian mantle

Near-liquidus phase equilibria experiments have been conducted on a synthetic Fastball basalt composition, as analyzed at Home Plate plateau of Mars (Gusev Crater), to test if it represents a primitive mantle derived melt and place constraints on the temperature of the ancient mantle and on the lithosphere-asthenosphere boundary of Mars. The Fastball basalt is multiply saturated with olivine and orthopyroxene at ∼1.2 GPa and 1430°C. Based on melting models, we predict that the Fastball composition could be produced by 13–23% equilibrium melting of the Martian mantle, with extraction of melt from the base of ∼105 km thick lithosphere. The multiple saturation for Fastball also constrains the potential temperature of the Martian mantle to be approximately 1480–1530°C with an initial melting pressure of 4.0–4.7 GPa. This potential temperature is much lower than that of the terrestrial mantle derived from similarly ancient magmas, i.e., komatiites.

The formation of Mercury's smooth plains

There has been extensive debate about whether Mercurys smooth plains are volcanic features or impact ejecta deposits. We present new indirect evidence which supports a volcanic origin for two different smooth plains units. In Borealis Planitia, stratigraphic relations indicate at least two distinct stages of smooth plains formation. At least one of these stages must have had a volcanic origin. In the Hilly and Lineated Terrain, Petrarch and several other anomalously shallow craters apparently have been volcanically filled. Areally extensive smooth plains volcanism evidently occurred at these two widely separated areas on Mercury. These results, combined with work by other researchers on the circum-Caloris plains and the Tolstoi basin, show that smooth plains volcanism was a global process on Mercury. Present data suggest to us that the smooth and intercrater plains may represent two distinct episodes of volcanic activity on Mercury and that smooth plains volcanism may have been triggered by the Caloris impact. High-resolution and multispectral imaging from a future Mercury spacecraft could resolve many of the present uncertainties in our understanding of plains formation on Mercury.