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 ).

The thickness of radar-bright deposits in Mercury's northern hemisphere from individual Mercury Laser Altimeter tracks.

The discovery of Mercury's radar-bright deposits has expanded our understanding of volatiles in the solar system. Key to deciphering the history and origin of the radar-bright deposits is an estimate of the volume of radar-bright material that in turn requires a measure of the average thickness of the deposits. In this study we investigate changes in topography across radar-bright deposits hosted in flat-floored, complex craters using individual edited Mercury Laser Altimeter (MLA) tracks. We compare the difference in heights of radar-bright regions and non-radar-bright regions of the crater floor and the difference of similarly sized and located regions in non-radar-bright craters and show that the two populations cannot be distinguished. The similarity of topography in these two sets of craters allows an upper limit of 15 m to be placed on the thickness of the radar-bright deposits.

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/m3, 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/m3. 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.

Age constraints of Mercury's polar deposits suggest recent delivery of ice.

Surface ice at the poles of Mercury appears as several-m-thick deposits that are composed of nearly pure water. We provide new age estimates of Mercury's polar deposits from combined analyses of Poisson statistics and direct observations of crater densities within permanently shadowed, radar-bright regions imaged by the MESSENGER spacecraft. These age estimates suggest that ice was delivered to Mercury within the last ~150 Myr. A single, recent impactor is one possible delivery mechanism that is consistent with our new age constraints, as well as the observed distinct reflectance boundaries of the polar deposits and the relative purity of the ice, as suggested by the Earth-based radar observations. In contrast to ice on Mercury, observations of the lunar poles are suggestive of a highly patchy distribution of surface frost. The patchiness of lunar polar deposits is consistent with long exposure times to the space weathering environment. Given enough time, the polar deposits on Mercury may age into a more heterogeneous spatial distribution, similar to that on the Moon.

Trilogy, a Planetary Geodesy Mission Concept for Measuring the Expansion of the Solar System.

The scale of the solar system is slowly changing, likely increasing as a result of solar mass loss, with additional change possible if there is a secular variation of the gravitational constant, G. The measurement of the change of scale could provide insight into the past and the future of the solar system, and in addition a better understanding of planetary motion and fundamental physics. Estimates for the expansion of the scale of the solar system are of order 1.5 cm year-1 AU-1, which over several years is an observable quantity with present-day laser ranging systems. This estimate suggests that laser measurements between planets could provide an accurate estimate of the solar system expansion rate. We examine distance measurements between three bodies in the inner solar system -- Earth's Moon, Mars and Venus -- and outline a mission concept for making the measurements. The concept involves placing spacecraft that carry laser ranging transponders in orbit around each body and measuring the distances between the three spacecraft over a period of several years. The analysis of these range measurements would allow the co-estimation of the spacecraft orbit, planetary ephemerides, other geophysical parameters related to the constitution and dynamics of the central bodies, and key geodetic parameters related to the solar system expansion, the Sun, and theoretical physics.

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.

Solar system expansion and strong equivalence principle as seen by the NASA MESSENGER mission.

The NASA MESSENGER mission explored the innermost planet of the solar system and obtained a rich data set of range measurements for the determination of Mercury's ephemeris. Here we use these precise data collected over 7 years to estimate parameters related to general relativity and the evolution of the Sun. These results confirm the validity of the strong equivalence principle with a significantly refined uncertainty of the Nordtvedt parameter η = (-6.6 ± 7.2) × 10-5. By assuming a metric theory of gravitation, we retrieved the post-Newtonian parameter ß = 1 + (-1.6 ± 1.8) × 10-5 and the Sun's gravitational oblateness, [Formula: see text] = (2.246 ± 0.022) × 10-7. Finally, we obtain an estimate of the time variation of the Sun gravitational parameter, [Formula: see text] = (-6.13 ± 1.47) × 10-14, which is consistent with the expected solar mass loss due to the solar wind and interior processes. This measurement allows us to constrain [Formula: see text] to be <4 × 10-14 per year.

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.

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.

ICESat/GLAS Altimetry Measurements: Received Signal Dynamic Range and Saturation Correction.

NASA's Ice, Cloud, and land Elevation Satellite (ICESat), which operated between 2003 and 2009, made the first satellite-based global lidar measurement of Earth's ice sheet elevations, sea-ice thickness and vegetation canopy structure. The primary instrument on ICESat was the Geoscience Laser Altimeter System (GLAS), which measured the distance from the spacecraft to Earth's surface via the roundtrip travel time of individual laser pulses. GLAS utilized pulsed lasers and a direct detection receiver consisting of a silicon avalanche photodiode (Si APD) and a waveform digitizer. Early in the mission, the peak power of the received signal from snow and ice surfaces was found to span a wider dynamic range than planned, often exceeding the linear dynamic range of the GLAS 1064-nm detector assembly. The resulting saturation of the receiver distorted the recorded signal and resulted in range biases as large as ~50 cm for ice and snow-covered surfaces. We developed a correction for this "saturation range bias" based on laboratory tests using a spare flight detector, and refined the correction by comparing GLAS elevation estimates to those derived from Global Positioning System (GPS) surveys over the calibration site at the salar de Uyuni, Bolivia. Applying the saturation correction largely eliminated the range bias due to receiver saturation for affected ICESat measurements over Uyuni and significantly reduced the discrepancies at orbit crossovers located on flat regions of the Antarctic ice sheet.

Evidence for surface water ice in the lunar polar regions using reflectance measurements from the Lunar Orbiter Laser Altimeter and temperature measurements from the Diviner Lunar Radiometer Experiment.

We find that the reflectance of the lunar surface within 5 ° of latitude of the South Pole increases rapidly with decreasing temperature, near ~110K, behavior consistent with the presence of surface water iceThe North polar region does not show this behavior, nor do South polar surfaces at latitudes more than 5° from the pole. This South pole reflectance anomaly persists when analysis is limited to surfaces with slopes less than 10° to eliminate false detection due to the brightening effect of mass wasting, and also when the very bright south polar crater Shackleton is excluded from the analysis. We also find that south polar regions of permanent shadow that have been reported to be generally brighter at 1064 nm do not show anomalous reflectance when their annual maximum surface temperatures are too high to preserve water ice. This distinction is not observed at the North Pole. The reflectance excursion on surfaces with maximum temperatures below 110K is superimposed on a general trend of increasing reflectance with decreasing maximum temperature that is present throughout the polar regions in the north and south; we attribute this trend to a temperature or illumination-dependent space weathering effect (e.g. Hemingway et al. 2015). We also find a sudden increase in reflectance with decreasing temperature superimposed on the general trend at 200K and possibly at 300K. This may indicate the presence of other volatiles such as sulfur or organics. We identified and mapped surfaces with reflectances so high as to be unlikely to be part of an ice-free population. In this south we find a similar distribution found by Hayne et al. 2015 based on UV properties. In the north a cluster of pixels near that pole may represent a limited frost exposure.

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.

Comparison of areas in shadow from imaging and altimetry in the north polar region of Mercury and implications for polar ice deposits.

Earth-based radar observations and results from the MESSENGER mission have provided strong evidence that permanently shadowed regions near Mercury's poles host deposits of water ice. MESSENGER's complete orbital image and topographic datasets enable Mercury's surface to be observed and modeled under an extensive range of illumination conditions. The shadowed regions of Mercury's north polar region from 65°N to 90°N were mapped by analyzing Mercury Dual Imaging System (MDIS) images and by modeling illumination with Mercury Laser Altimeter (MLA) topographic data. The two independent methods produced strong agreement in identifying shadowed areas. All large radar-bright deposits, those hosted within impact craters ≥6 km in diameter, collocate with regions of shadow identified by both methods. However, only ∼46% of the persistently shadowed areas determined from images and ∼43% of the permanently shadowed areas derived from altimetry host radar-bright materials. Some sizable regions of shadow that do not host radar-bright deposits experience thermal conditions similar to those that do. The shadowed craters that lack radar-bright materials show a relation with longitude that is not related to the thermal environment, suggesting that the Earth-based radar observations of these locations may have been limited by viewing geometry, but it is also possible that water ice in these locations is insulated by anomalously thick lag deposits or that these shadowed regions do not host water ice.

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.

Planetary science. Low-altitude magnetic field measurements by MESSENGER reveal Mercury's ancient crustal field.

Magnetized rocks can record the history of the magnetic field of a planet, a key constraint for understanding its evolution. From orbital vector magnetic field measurements of Mercury taken by the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft at altitudes below 150 kilometers, we have detected remanent magnetization in Mercury's crust. We infer a lower bound on the average age of magnetization of 3.7 to 3.9 billion years. Our findings indicate that a global magnetic field driven by dynamo processes in the fluid outer core operated early in Mercury's history. Ancient field strengths that range from those similar to Mercury's present dipole field to Earth-like values are consistent with the magnetic field observations and with the low iron content of Mercury's crust inferred from MESSENGER elemental composition data.

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.

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.

Asymmetric distribution of lunar impact basins caused by variations in target properties.

Maps of crustal thickness derived from NASA's Gravity Recovery and Interior Laboratory (GRAIL) mission revealed more large impact basins on the nearside hemisphere of the Moon than on its farside. The enrichment in heat-producing elements and prolonged volcanic activity on the lunar nearside hemisphere indicate that the temperature of the nearside crust and upper mantle was hotter than that of the farside at the time of basin formation. Using the iSALE-2D hydrocode to model impact basin formation, we found that impacts on the hotter nearside would have formed basins with up to twice the diameter of similar impacts on the cooler farside hemisphere. The size distribution of lunar impact basins is thus not representative of the earliest inner solar system impact bombardment.