Long-lived lunar volcanism sustained by precession-driven core-mantle friction

Core-mantle friction induced by the precession of the Moon's spin axis is a strong heat source in the deep lunar mantle during the early phase of a satellite's evolution, but its influence on the long-term thermal evolution still remains poorly explored. Using a one-dimensional thermal evolution model, we show that core-mantle friction can sustain global-scale partial melting in the upper lunar mantle until ∼3.1 Ga, thus accounting for the intense volcanic activity on the Moon before ∼3.0 Ga. Besides, core-mantle friction tends to suppress the secular cooling of the lunar core and is unlikely to be an energy source for the long-lived lunar core dynamo. Our model also favours the transition of the Cassini state before the end of the lunar magma ocean phase (∼4.2 Ga), which implies a decreasing lunar obliquity over time after the solidification of the lunar magma ocean. Such a trend of lunar obliquity evolution may allow volcanically released water to be buried in the lunar regolith of the polar regions. As a consequence, local water ice could be more abundant than previously thought when considering only its accumulation caused by solar wind and comet spreading.

Gravity evidence for a heterogeneous crust of Mercury

We modeled gravity data to explore Mercury's internal structure and show the presence of crustal heterogeneities in density. We first evaluated the lithospheric flexure occurring in the spherical harmonic degree range 5-80, according to the flexural isostatic response curve. We thus estimated a mean elastic lithosphere thickness of about 30 [Formula: see text] 10 km and modeled the crust-mantle interface, which varies from 19 to 42 km depth, according to a flexural compensation model. The isostatic gravity anomalies were then obtained as the residual field with respect to the contributions from topography and lithospheric flexure. Isostatic anomalies are mainly related to density variations in the crust: gravity highs mostly correspond to large-impact basins suggesting intra-crustal magmatic intrusions as the main origin of these anomalies. Isostatic gravity lows prevail, instead, above intercrater plains and may represent the signature of a heavily fractured crust.

Remote detection of a lunar granitic batholith at Compton-Belkovich

Granites are nearly absent in the Solar System outside of Earth. Achieving granitic compositions in magmatic systems requires multi-stage melting and fractionation, which also increases the concentration of radiogenic elements1. Abundant water and plate tectonics facilitate these processes on Earth, aiding in remelting. Although these drivers are absent on the Moon, small granite samples have been found, but details of their origin and the scale of systems they represent are unknown2. Here we report microwave-wavelength measurements of an anomalously hot geothermal source that is best explained by the presence of an approximately 50-kilometre-diameter granitic system below the thorium-rich farside feature known as Compton-Belkovich. Passive microwave radiometry is sensitive to the integrated thermal gradient to several wavelengths depth. The 3-37-gigahertz antenna temperatures of the Chang'e-1 and Chang'e-2 microwave instruments allow us to measure a peak heat flux of about 180 milliwatts per square metre, which is about 20 times higher than that of the average lunar highlands3,4. The surprising magnitude and geographic extent of this feature imply an Earth-like, evolved granitic system larger than believed possible on the Moon, especially outside of the Procellarum region5. Furthermore, these methods are generalizable: similar uses of passive radiometric data could vastly expand our knowledge of geothermal processes on the Moon and other planetary bodies.

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.

Mercury's Crustal Thickness and Contractional Strain

The crust of Mercury has experienced contraction on a global scale. Contractional deformation is expressed by a broadly distributed network of lobate thrust fault scarps. The most likely principal source of stress is global contraction from cooling of Mercury's interior. Global contraction alone would be expected to result in a uniformly distributed population of thrust faults. Mercury's fault scarps, however, often occur in long, linear clusters or bands. An analysis of the contractional strain as a function of crustal thickness, estimated in two crustal thickness models (CT1 and CT2) derived from gravity and topography data obtained during the MESSENGER mission, indicates the greatest contractional strain occurs in crust 50-60 km thick. On Earth, mantle downwelling can thicken and compress overlying crust, regionally concentrating thrust faults. Clusters of lobate scarps collocated with regions of thick crust suggest downward mantle flow contributed to the localization of lithosphere-penetrating thrust faults.

Thickness and structure of the martian crust from InSight seismic data

A planet's crust bears witness to the history of planetary formation and evolution, but for Mars, no absolute measurement of crustal thickness has been available. Here, we determine the structure of the crust beneath the InSight landing site on Mars using both marsquake recordings and the ambient wavefield. By analyzing seismic phases that are reflected and converted at subsurface interfaces, we find that the observations are consistent with models with at least two and possibly three interfaces. If the second interface is the boundary of the crust, the thickness is 20 ± 5 kilometers, whereas if the third interface is the boundary, the thickness is 39 ± 8 kilometers. Global maps of gravity and topography allow extrapolation of this point measurement to the whole planet, showing that the average thickness of the martian crust lies between 24 and 72 kilometers. Independent bulk composition and geodynamic constraints show that the thicker model is consistent with the abundances of crustal heat-producing elements observed for the shallow surface, whereas the thinner model requires greater concentration at depth.

Groundwater production from geothermal heating on early Mars and implication for early martian habitability

In explaining extensive evidence for past liquid water, the debate on whether Mars was primarily warm and wet or cold and arid 4 billion years (Ga) ago has continued for decades. The Sun's luminosity was ~30% lower 4 Ga ago; thus, most martian climate models struggle to elevate the mean surface temperature past the melting point of water. Basal melting of ice sheets may help resolve that paradox. We modeled the thermophysical evolution of ice and estimate the geothermal heat flux required to produce meltwater on a cold, arid Mars. We then analyzed geophysical and geochemical data, showing that basal melting would have been feasible on Mars 4 Ga ago. If Mars were warm and wet 4 Ga ago, then the geothermal flux would have even sustained hydrothermal activity. Regardless of the actual nature of the ancient martian climate, the subsurface would have been the most habitable region on Mars.

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.

Geodetic evidence that Mercury has a solid inner core

Geodetic analysis of radio tracking measurements of the MESSENGER spacecraft while in orbit about Mercury has yielded new estimates for the planet's gravity field, tidal Love number, and pole coordinates. The derived right ascension (α = 281.0082° ± 0.0009°; all uncertainties are 3 standard deviations) and declination (δ =61.4164° ± 0.0003°) of the spin pole place Mercury in the Cassini state. Confirmation of the equilibrium state with an estimated mean (whole-planet) obliquity ϵ of 1.968 ± 0.027 arcmin enables the confident determination of the planet's normalized polar moment of inertia (0.333 ± 0.005), which indicates a high degree of internal differentiation. Internal structure models generated by a Markov-Chain Monte Carlo process and consistent with the geodetic constraints possess a solid inner core with a radius (r _ic ) between 0.3 and 0.7 that of the outer core (r _oc ).

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.

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.

Present-day heat flow model of Mars

Until the acquisition of in-situ measurements, the study of the present-day heat flow of Mars must rely on indirect methods, mainly based on the relation between the thermal state of the lithosphere and its mechanical strength, or on theoretical models of internal evolution. Here, we present a first-order global model for the present-day surface heat flow for Mars, based on the radiogenic heat production of the crust and mantle, on scaling of heat flow variations arising from crustal thickness and topography variations, and on the heat flow derived from the effective elastic thickness of the lithosphere beneath the North Polar Region. Our preferred model finds heat flows varying between 14 and 25 mW m^-2, with an average value of 19 mW m^-2. Similar results (although about ten percent higher) are obtained if we use heat flow based on the lithospheric strength of the South Polar Region. Moreover, expressing our results in terms of the Urey ratio (the ratio between total internal heat production and total heat loss through the surface), we estimate values close to 0.7-0.75, which indicates a moderate contribution of secular cooling to the heat flow of Mars (consistent with the low heat flow values deduced from lithosphere strength), unless heat-producing elements abundances for Mars are subchondritic.

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.

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.

Lunar exploration: opening a window into the history and evolution of the inner Solar System

The lunar geological record contains a rich archive of the history of the inner Solar System, including information relevant to understanding the origin and evolution of the Earth-Moon system, the geological evolution of rocky planets, and our local cosmic environment. This paper provides a brief review of lunar exploration to-date and describes how future exploration initiatives will further advance our understanding of the origin and evolution of the Moon, the Earth-Moon system and of the Solar System more generally. It is concluded that further advances will require the placing of new scientific instruments on, and the return of additional samples from, the lunar surface. Some of these scientific objectives can be achieved robotically, for example by in situ geochemical and geophysical measurements and through carefully targeted sample return missions. However, in the longer term, we argue that lunar science would greatly benefit from renewed human operations on the surface of the Moon, such as would be facilitated by implementing the recently proposed Global Exploration Roadmap.

Lunar bulk chemical composition: a post-Gravity Recovery and Interior Laboratory reassessment

New estimates of the thickness of the lunar highlands crust based on data from the Gravity Recovery and Interior Laboratory mission, allow us to reassess the abundances of refractory elements in the Moon. Previous estimates of the Moon fall into two distinct groups: earthlike and a 50% enrichment in the Moon compared with the Earth. Revised crustal thicknesses and compositional information from remote sensing and lunar samples indicate that the crust contributes 1.13-1.85 wt% Al2O3 to the bulk Moon abundance. Mare basalt Al2O3 concentrations (8-10 wt%) and Al2O3 partitioning behaviour between melt and pyroxene during partial melting indicate mantle Al2O3 concentration in the range 1.3-3.1 wt%, depending on the relative amounts of pyroxene and olivine. Using crustal and mantle mass fractions, we show that that the Moon and the Earth most likely have the same (within 20%) concentrations of refractory elements. This allows us to use correlations between pairs of refractory and volatile elements to confirm that lunar abundances of moderately volatile elements such as K, Rb and Cs are depleted by 75% in the Moon compared with the Earth and that highly volatile elements, such as Tl and Cd, are depleted by 99%. The earthlike refractory abundances and depleted volatile abundances are strong constraints on lunar formation processes.

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.

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.

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

Spacecraft-to-spacecraft tracking observations from the Gravity Recovery and Interior Laboratory (GRAIL) have been used to construct a gravitational field of the Moon to spherical harmonic degree and order 420. The GRAIL field reveals features not previously resolved, including tectonic structures, volcanic landforms, basin rings, crater central peaks, and numerous simple craters. From degrees 80 through 300, over 98% of the gravitational signature is associated with topography, a result that reflects the preservation of crater relief in highly fractured crust. The remaining 2% represents fine details of subsurface structure not previously resolved. GRAIL elucidates the role of impact bombardment in homogenizing the distribution of shallow density anomalies on terrestrial planetary bodies.

Internal structure and early thermal evolution of Mars from Mars Global Surveyor topography and gravity

Topography and gravity measured by the Mars Global Surveyor have enabled determination of the global crust and upper mantle structure of Mars. The planet displays two distinct crustal zones that do not correlate globally with the geologic dichotomy: a region of crust that thins progressively from south to north and encompasses much of the southern highlands and Tharsis province and a region of approximately uniform crustal thickness that includes the northern lowlands and Arabia Terra. The strength of the lithosphere beneath the ancient southern highlands suggests that the northern hemisphere was a locus of high heat flow early in martian history. The thickness of the elastic lithosphere increases with time of loading in the northern plains and Tharsis. The northern lowlands contain structures interpreted as large buried channels that are consistent with northward transport of water and sediment to the lowlands before the end of northern hemisphere resurfacing.

Improved gravity field of the moon from lunar prospector

An improved gravity model from Doppler tracking of the Lunar Prospector (LP) spacecraft reveals three new large mass concentrations (mascons) on the nearside of the moon beneath the impact basins Mare Humboltianum, Mendel-Ryberg, and Schiller-Zucchius, where the latter basin has no visible mare fill. Although there is no direct measurement of the lunar farside gravity, LP partially resolves four mascons in the large farside basins of Hertzsprung, Coulomb-Sarton, Freundlich-Sharonov, and Mare Moscoviense. The center of each of these basins contains a gravity maximum relative to the surrounding basin. The improved normalized polar moment of inertia (0.3932 +/- 0.0002) is consistent with an iron core with a radius of 220 to 450 kilometers.