Exploring the Interior of Europa with the Europa Clipper.

The Galileo mission to Jupiter revealed that Europa is an ocean world. The Galileo magnetometer experiment in particular provided strong evidence for a salty subsurface ocean beneath the ice shell, likely in contact with the rocky core. Within the ice shell and ocean, a number of tectonic and geodynamic processes may operate today or have operated at some point in the past, including solid ice convection, diapirism, subsumption, and interstitial lake formation. The science objectives of the Europa Clipper mission include the characterization of Europa's interior; confirmation of the presence of a subsurface ocean; identification of constraints on the depth to this ocean, and on its salinity and thickness; and determination of processes of material exchange between the surface, ice shell, and ocean. Three broad categories of investigation are planned to interrogate different aspects of the subsurface structure and properties of the ice shell and ocean: magnetic induction, subsurface radar sounding, and tidal deformation. These investigations are supplemented by several auxiliary measurements. Alone, each of these investigations will reveal unique information. Together, the synergy between these investigations will expose the secrets of the Europan interior in unprecedented detail, an essential step in evaluating the habitability of this ocean world.

Small All-Range Lidar for Asteroid and Comet Core Missions.

We report the development of a new type of space lidar specifically designed for missions to small planetary bodies for both topographic mapping and support of sample collection or landing. The instrument is designed to have a wide dynamic range with several operation modes for different mission phases. The laser transmitter consists of a fiber laser that is intensity modulated with a return-to-zero pseudo-noise (RZPN) code. The receiver detects the coded pulse-train by correlating the detected signal with the RZPN kernel. Unlike regular pseudo noise (PN) lidars, the RZPN kernel is set to zero outside laser firing windows, which removes most of the background noise over the receiver integration time. This technique enables the use of low peak-power but high pulse-rate lasers, such as fiber lasers, for long-distance ranging without aliasing. The laser power and the internal gain of the detector can both be adjusted to give a wide measurement dynamic range. The laser modulation code pattern can also be reconfigured in orbit to optimize measurements to different measurement environments. The receiver uses a multi-pixel linear mode photon-counting HgCdTe avalanche photodiode (APD) array with near quantum limited sensitivity at near to mid infrared wavelengths where many fiber lasers and diode lasers operate. The instrument is modular and versatile and can be built mostly with components developed by the optical communication industry.

New Illumination and Temperature Constraints of Mercury's Volatile Polar Deposits.

Images from the Mercury Dual Imaging System (MDIS) aboard the MErcury Surface, Space ENvironment, GEochemistry, and Ranging mission reveal low-reflectance polar deposits that are interpreted to be lag deposits of organic-rich, volatile material. Interpretation of these highest-resolution images of Mercury's polar deposits has been limited by the available topography models, so local high-resolution (125 m pixel-1) digital elevation models (DEMs) were made using a combination of data from the Mercury Laser Altimeter (MLA) and from shape-from-shading techniques using MDIS images. Local DEMs were made for eight of Mercury's north polar craters; these DEMs were then used to create high-resolution simulated image, illumination, and thermal models. The simulated images reveal that the pixel brightness variations imaged within Mercury's low-reflectance deposits are consistent with scattered light reflecting off of topography and do not need to be explained by volatile compositional differences as previously suggested. The illumination and thermal models show that these low-reflectance polar deposits extend beyond the permanently shadowed region, more than 1.0 km in some locations, and correspond to a maximum surface temperature of greater than 250 K but less than 350 K. The low-reflectance boundaries of all eight polar deposits studied here show a close correspondence with the surface stability boundary of coronene (C24H12). While coronene should only be viewed as a proxy for the myriad volatile compounds that may exist in Mercury's polar deposits, coronene's surface stability boundary supports the idea that Mercury's low-reflectance polar deposits are composed of macromolecular organic compounds, consistent with the hypotheses of exogenous transport and in situ production.

Optical characterization of laser retroreflector arrays for lunar landers.

The Laser Retroreflector Array for Lunar Landers (LRALL) is a small optical instrument designed to provide a target for precision laser ranging from a spacecraft in lunar orbit, enabling geolocation of the lander and its instrument suite and establishing a fiducial maker on the lunar surface. Here we describe the optical performance of LRALL at visible and near-infrared wavelengths. Individual corner cube reflectors (CCRs) within LRALL were tested for surface flatness and dihedral angle values. We also imaged the far-field diffraction patterns of individual CCRs as well as the entire retroreflector array over the range of possible incident angles to extract the optical cross section as a function of viewing angle. We also measured the optical properties of one of the CCRs over the lunar temperature range (100-380 K) and found no significant temperature-dependent variance. The test results show LRALL meets the design criteria and can be ranged to elevation angles above 30° with respect to the instrument base from an orbital laser altimeter such as the Lunar Orbiter Laser Altimeter on the Lunar Reconnaissance Orbiter. This work summarizes the test data and serves as a guide for future laser ranging to these retroreflector arrays.

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.

Small and lightweight laser retro-reflector arrays for lunar landers.

A set of small and lightweight laser retro-reflector arrays (LRAs) was fabricated and tested for use on lunar landers under NASA's Commercial Lunar Payload Service program. Each array contains eight 1.27-cm-diameter corner cube retro-reflectors mounted on a dome-shaped aluminum structure. The arrays are 5.0 cm in diameter at the base, 1.6 cm in height, and 20 g in mass. They can be tracked by an orbiting laser altimeter, such as the Lunar Orbiter Laser Altimeter, from a distance of a few hundred kilometers or by a landing lidar on future lunar landers. The LRAs demonstrated a diffraction-limited optical performance. They were designed and tested to survive and function on the Moon for decades, well after the lander missions are completed.

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.

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.

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.

Advanced Illumination Modeling for Data Analysis and Calibration. Application to the Moon.

We present a new illumination modeling tool, called IllumNG, developed at NASA Goddard Space Flight Center (GSFC). We describe its capabilities to enhance the analysis and calibration of science data collected by planetary missions. We highlight these with examples making use of lunar data, particularly the topographic and radiometric measurements collected by the Lunar Orbiter Laser Altimeter (LOLA) instrument, with applications to radiometric measurements from other LRO instruments as well. The unique features of IllumNG are its accuracy and flexibility to handle multiple types of observers and light sources, and its ability to accurately model both singly- and doubly-scattered radiation to an observer.

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.

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.

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.