The MESSENGER mission to Mercury: scientific objectives and implementation

Mercury holds answers to several critical questions regarding the formation and evolution of the terrestrial planets. These questions include the origin of Mercury's anomalously high ratio of metal to silicate and its implications for planetary accretion processes, the nature of Mercury's geological evolution and interior cooling history, the mechanism of global magnetic field generation, the state of Mercury's core, and the processes controlling volatile species in Mercury's polar deposits, exosphere, and magnetosphere. The MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) mission has been designed to fly by and orbit Mercury to address all of these key questions. After launch by a Delta 2925H-9.5, two flybys of Venus, and two flybys of Mercury, orbit insertion is accomplished at the third Mercury encounter. The instrument payload includes a dual imaging system for wide and narrow fields-of-view, monochrome and color imaging, and stereo; X-ray and combined gamma-ray and neutron spectrometers for surface chemical mapping; a magnetometer; a laser altimeter; a combined ultraviolet–visible and visible-near-infrared spectrometer to survey both exospheric species and surface mineralogy; and an energetic particle and plasma spectrometer to sample charged species in the magnetosphere. During the flybys of Mercury, regions unexplored by Mariner 10 will be seen for the first time, and new data will be gathered on Mercury's exosphere, magnetosphere, and surface composition. During the orbital phase of the mission, one Earth year in duration, MESSENGER will complete global mapping and the detailed characterization of the exosphere, magnetosphere, surface, and interior.     

Analysis of MESSENGER Gamma-Ray Spectrometer data from the Mercury flybys

During its three flybys of Mercury, the MESSENGER spacecraft made the first detection of gamma-ray emission from the planet's surface. With a closest approach distance of ∼200 km, the flybys provided an opportunity to measure elemental abundances of Mercury's near-equatorial regions, which will not be visited at low altitude during MESSENGER's orbital mission phase. Despite being limited by low planetary photon flux, sufficient counts were accumulated during the first two flybys to estimate bounds on abundances for some elements having relatively strong gamma-ray spectral peaks, including Si, Fe, Ti, K, and Th. Only for Si is the standard deviation σ sufficiently small to conclude that this element was detected with 99% confidence. Iron and potassium are detected at the 2−σ (95% confidence) level, whereas only upper bounds on Ti and Th can be determined. Relative to a Si abundance assumed to be 18 weight percent (wt%), 2−σ upper bounds have been estimated as 9.7 wt% for Fe, 7.0 wt% for Ti, 0.087 wt% for K, and 2.2 ppm for Th. The relatively low upper bound on K rules out some previously suggested models for surface composition for the regions sampled. Upper bounds on Fe/Si and Ti/Si ratios are generally consistent with Ti and Fe abundances estimated from the analysis of measurements by the MESSENGER Neutron Spectrometer during the flybys but are also permissive of much lower concentrations.     

Determination of the properties of Mercury's magnetic field by the MESSENGER mission

The MESSENGER mission to Mercury, to be launched in 2004, will provide an opportunity to characterize Mercury's internal magnetic field during an orbital phase lasting one Earth year. To test the ability to determine the planetary dipole and higher-order moments from measurements by the spacecraft's fluxgate magnetometer, we simulate the observations along the spacecraft trajectory and recover the internal field characteristics from the simulated observations. The magnetic field inside Mercury's magnetosphere is assumed to consist of an intrinsic multipole component and an external contribution due to magnetospheric current systems described by a modified Tsyganenko 96 model. Under the axis-centered-dipole approximation without correction for the external field the moment strength is overestimated by ∼4% for a simulated dipole moment of , and the error depends strongly on the magnitude of the simulated moment, rising as the moment decreases. Correcting for the external field contributions can reduce the error in the dipole term to a lower limit of ∼1-2% without a solar wind monitor. Dipole and quadrupole terms, although highly correlated, are then distinguishable at the level equivalent to an error in the position of an offset dipole of a few tens of kilometers. Knowledge of the external magnetic field is therefore the primary limiting factor in extracting reliable knowledge of the structure of Mercury's magnetic field from the MESSENGER observations.     

On the oscillations in Mercury's obliquity

Mercury's capture into the 3:2 spin–orbit resonance can be explained as a result of its chaotic orbital dynamics. One major objective of MESSENGER and BepiColombo spatial missions is to accurately measure Mercury's rotation and its obliquity in order to obtain constraints on internal structure of the planet. Analytical approaches at the first-order level using the Cassini state assumptions give the obliquity constant or quasi-constant. Which is the obliquity's dynamical behavior deriving from a complete spin–orbit motion of Mercury simultaneously integrated with planetary interactions? We have used our SONYR model (acronym of Spin–Orbit N-bodY Relativistic model) integrating the spin–orbit N-body problem applied to the Solar System (Sun and planets). For lack of current accurate observations or ephemerides of Mercury's rotation, and therefore for lack of valid initial conditions for a numerical integration, we have built an original method for finding the libration center of the spin–orbit system and, as a consequence, for avoiding arbitrary amplitudes in librations of the spin–orbit motion as well as in Mercury's obliquity. The method has been carried out in two cases: (1) the spin–orbit motion of Mercury in the 2-body problem case (Sun–Mercury) where an uniform precession of the Keplerian orbital plane is kinematically added at a fixed inclination (S_2K case), (2) the spin–orbit motion of Mercury in the N-body problem case (Sun and planets) (S_n case). We find that the remaining amplitude of the oscillations in the S_n case is one order of magnitude larger than in the S_2K case, namely 4 versus 0.4 arcseconds (peak-to-peak). The mean obliquity is also larger, namely 1.98 versus 1.80 arcminutes, for a difference of 10.8 arcseconds. These theoretical results are in a good agreement with recent radar observations but it is not excluded that it should be possible to push farther the convergence process by drawing nearer still more precisely to the libration center. We note that the dynamically driven spin precession, which occurs when the planetary interactions are included, is more complex than the purely kinematic case. Nevertheless, in such a N-body problem, we find that the 3:2 spin–orbit resonance is really combined to a synchronism where the spin and orbit poles on average precess at the same rate while the orbit inclination and the spin axis orientation on average decrease at the same rate. As a consequence and whether it would turn out that there exists an irreducible minimum of the oscillation amplitude, quasi-periodic oscillations found in Mercury's obliquity should be to geometrically understood as librations related to these synchronisms that both follow a Cassini state. Whatever the open question on the minimal amplitude in the obliquity's oscillations and in spite of the planetary interactions indirectly acting by the solar torque on Mercury's rotation, Mercury remains therefore in a stable equilibrium state that proceeds from a 2-body Cassini state.     

The proximity of Mercury's spin to Cassini state 1 from adiabatic invariance

In determining Mercury's core structure from its rotational properties, the value of the normalized moment of inertia, C/MR², from the location of Cassini 1 is crucial. If Mercury's spin axis occupies Cassini state 1, its position defines the location of the state, where the axis is fixed in the frame precessing with the orbit. Although tidal and core-mantle dissipation drive the spin to the Cassini state with a time scale O(¹⁰⁵) years, the spin might still be displaced from the Cassini state if the variations in the orbital elements induced by planetary perturbations, which change the position of the Cassini state, cause the spin to lag behind as it attempts to follow the state. After being brought to the state by dissipative processes, the spin axis is expected to follow the Cassini state for orbit variations with time scales long compared to the 1000 year precession period of the spin about the Cassini state because the solid angle swept out by the spin axis as it precesses is an adiabatic invariant. Short period variations in the orbital elements of small amplitude should cause displacements that are commensurate with the amplitudes of the short period terms. The exception would be if there are forcing terms in the perturbations that are nearly resonant with the 1000 year precession period. The precision of the radar and eventual spacecraft measurements of the position of Mercury's spin axis warrants a check on the likely proximity of the spin axis to the Cassini state. How confident should we be that the spin axis position defines the Cassini state sufficiently well for a precise determination of C/MR²? By following simultaneously the spin position and the Cassini state position during long time scale orbital variations over past 3 million years [Quinn, T.R., Tremaine, S., Duncan, M., 1991. Astron. J. 101, 2287-2305] and short time scale variations for 20,000 years [JPL Ephemeris DE 408; Standish, E.M., private communication, 2005], we show that the spin axis will remain within one arcsec of the Cassini state after it is brought there by dissipative torques. In this process the spin is located in the orbit frame of reference, which in turn is referenced to the inertial ecliptic plane of J2000. There are no perturbations with periods resonant with the precession period that could cause large separations. We thus expect Mercury's spin to occupy Cassini state 1 well within the uncertainties for both radar and spacecraft measurements, with correspondingly tight constraints on C/MR² and the extent of Mercury's molten core. Two unlikely caveats for this conclusion are: (1) an excitation of a free spin precession by an unknown mechanism or (2) a displacement by a dissipative core mantle interaction that exceeds the measurement uncertainties.     

The global distribution of pyroclastic deposits on Mercury: The view from MESSENGER flybys 1–3

We present a global survey of candidate pyroclastic deposits on Mercury, derived from images obtained during MESSENGER flybys 1–3 that provided near-global coverage at resolutions between 5 and 0.5 km/pixel. Thirty-five deposits were identified and characterized and are located principally on the floors of craters, along rims of craters, and along the edge of the Caloris basin. Deposits are commonly centered on rimless, often irregularly shaped pits, mostly between 5 and 45 km in diameter. The deposits identified are generally similar in morphology and absolute reflectance to lunar pyroclastic deposits. Spectrally the deposits appear brighter and redder than background Mercury terrain. On the basis of the available coverage, the candidate pyroclastic deposits appear to be essentially globally distributed. The diameters of the deposits, when mapped to lunar gravity conditions, are larger than their lunar counterparts, implying that more abundant volatiles were present during the typical eruptive process than on the Moon. These observations indicate that if these deposits resulted from hawaiian-style eruptions, the volatile contents required would be between ～1600 and 16,000 ppm CO or an equivalent value of H₂O, CO₂, SO₂, or H₂S (for a more oxidizing interior), or N₂, S₂, CS₂, S₂Cl, Cl, Cl₂, or COS (for a more reducing interior). These abundances are much greater than those predicted by existing models for Mercury's formation. An apparent lack of small deposits, compared with the Moon, may be due to resolution effects, a topic that can be further assessed during the orbital phase of the MESSENGER mission. These results provide a framework within which orbital observations by MESSENGER and the future BepiColombo mission can be analyzed.     

Global MHD modeling of Mercury's magnetosphere with applications to the MESSENGER mission and dynamo theory

We use a global magnetohydrodynamic (MHD) model to simulate Mercury's space environment for several solar wind and interplanetary magnetic field (IMF) conditions in anticipation of the magnetic field measurements by the MESSENGER spacecraft. The main goal of our study is to assess what characteristics of the internally generated field of Mercury can be inferred from the MESSENGER observations, and to what extent they will be able to constrain various models of Mercury's magnetic field generation. Based on the results of our simulations, we argue that it should be possible to infer not only the dipole component, but also the quadrupole and possibly even higher harmonics of the Mercury's planetary magnetic field. We furthermore expect that some of the crucial measurements for specifying the Hermean internal field will be acquired during the initial fly-bys of the planet, before MESSENGER goes into orbit around Mercury.     

Measurement of the radius of Mercury by radio occultation during the MESSENGER flybys

The MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft completed three flybys of Mercury in 2008–2009. During the first and third of those flybys, MESSENGER passed behind the planet from the perspective of Earth, occulting the radio-frequency (RF) transmissions. The occultation start and end times, recovered with 0.1 s accuracy or better by fitting edge-diffraction patterns to the RF power history, are used to estimate Mercury's radius at the tangent point of the RF path. To relate the measured radius to the planet shape, we evaluate local topography using images to identify the high-elevation feature that defines the RF path or using altimeter data to quantify surface roughness. Radius measurements are accurate to 150 m, and uncertainty in the average radius of the surrounding terrain, after adjustments are made from the local high at the tangent point of the RF path, is 350 m. The results are consistent with Mercury's equatorial shape as inferred from observations by the Mercury Laser Altimeter and ground-based radar. The three independent estimates of radius from occultation events collectively yield a mean radius for Mercury of 2439.2±0.5 km.     

Stereo topographic models of Mercury after three MESSENGER flybys

From photogrammetric analysis of stereo images of Mercury obtained during three MESSENGER flybys, we have produced three digital terrain models (DTMs) that have a grid spacing of 1 km and together cover 30% of the planet's surface. The terrain models provide a rich source of information on the morphology of Mercury's surface, including details of tectonic scarp systems as well as impact craters and basins. More than 400 craters larger than 15 km in diameter are included in the models. Additionally, the models provide important test cases for the analysis of stereo image data to be collected during MESSENGER's orbital mission phase. Small lateral offsets and differences in trends between stereo DTMs and laser altimeter profiles may be due to remaining errors in spacecraft position, instrument pointing, or Mercury coordinate knowledge. Such errors should be resolved during the orbital mission phase, when more joint analyses of data and detailed orbit modeling will be possible.     

Gemini J ‐band observations of RX J0806.4–4123

The detection of near‐infrared (NIR) excess at the position of a star can indicate either a substellar companion or a disk around the respective star. In this work we probed whether a 2.5σ H ‐band flux enhancement at the position of the isolated neutron star RX J0806.4–4123 can be confirmed at another NIR wavelength. We observed RXJ0806.4–4123 in the J ‐band with Gemini South equipped with FLAMINGOS‐2. There was no significant detection of a J ‐band source at the neutron star position. However, similarly to the H ‐band we found a very faint (1.4σ) flux enhancement with a nominal magnitude of J = 24.8 ± 0.5. The overall NIR‐detection significance is 3.1σ. If real, this emission is too bright to come from the neutron star alone. Deeper near‐infrared observations are necessary to confirm or refute the potential NIR excess. The confirmation of such NIR excess could imply that there is a substellar companion or a disk around RXJ0806.4–4123. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Viking bistatic radar experiment: Summary of first-order results emphasizing north polar data

Bistatic radar observations of Mars' north polar region during 1977–1978 showed surface rms slope σ_β ranging from 1 to 6°; these values apply to horizontal scales of 1–100 m. Values of roughness tend to decrease with increasing latitude (especially over 65–80°N), but there are many exceptions. The smoothest surfaces (σ_β≤1°) appear to be inclusions within generally rougher (σ_β～3°) terrain, rather than broad expanses of very smooth material. The permanent north polar cap is relatively uniform with 2.5?σ_β?3.0°. Considerable structure has been found in echo spectra, indicating a heterogeneous and perhaps anisotropic scattering surface. Echo spectra obtained from the same region, but several months apart (1°<L_S<62°), show no significant differences in inferred roughness. Estimates of reflectivity and dielectric constant are systematically low in the polar region. This may indicate that surface material north of 65°N is less dense than that near the equator, but more study of these data is needed. Estimates of surface roughness and dielectric constant in the equatorial region are consistent with results from Earth-based measurements to the accuracy of our analysis.     

Long-period forcing of Mercury's libration in longitude

Planetary perturbations of Mercury's orbit lead to forced librations in longitude in addition to the 88-day forced libration induced by Mercury's orbital motion. The forced librations are a combination of many periods, but 5.93 and 5.66 years dominate. These two periods result from the perturbations by Jupiter and Venus respectively, and they lead to a 125-year modulation of the libration amplitude corresponding to the beat frequency. Other periods are also identified with Jupiter and Venus perturbations as well as with those of the Earth, and these and other periods in the perturbations cause several arc second fluctuations in the libration extremes. The maxima of these extremes are about 30″ above and the minima about 7″ above the superposed ∼^60″, 88-day libration during the 125-year modulation. Knowledge of the nature of the long-period forced librations is important for the interpretation of the details of Mercury's rotation state to be obtained from radar and spacecraft observations. We show that the measurement of the 88-day libration amplitude for the purposes of determining Mercury's core properties is not compromised by the additional librations, because of the latter's small amplitude and long period. If the free libration in longitude has an amplitude that is large compared with that of the forced libration, its ∼10-year period will dominate the libration spectrum with the 88-day forced libration and the long-period librations from the orbital perturbations superposed. If the free libration has an amplitude that is comparable to those of the long-period forced libration, it will be revealed by erratic amplitude, period and phase on the likely time span of a series of observations. However, a significant free libration component is not expected because of relatively rapid damping.     

Mapping iron abundances on the surface of Mercury: Predicted spatial resolution of the MESSENGER Gamma-Ray Spectrometer

To illustrate the spatial resolution of measurements of Mercury's surface elemental composition by the Gamma-Ray Spectrometer on the MESSENGER spacecraft after one year of orbital observations, we have simulated a global coverage map of the 846-keV iron gamma-ray count rate. The simulated map suggests that distinct geologic units larger than 800 km in horizontal dimension will be discernable when the difference in Fe abundance between adjacent geologic units exceeds 4 wt%. These results imply that the MESSENGER Gamma-Ray Spectrometer dataset will provide useful information for regional geological studies of the surface of Mercury.     

Radar imagery of Mercury’s putative polar ice: 1999-2005 Arecibo results

We present an updated survey of Mercury’s putative polar ice deposits, based on high-resolution (1.5-km) imaging with the upgraded Arecibo S-band radar during 1999-2005. The north pole has now been imaged over a full range of longitude aspects, making it possible to distinguish ice-free areas from radar-shadowed areas and thus better map the distribution of radar-bright ice. The new imagery of the south pole, though derived from only a single pair of dates in 2005, improves on the pre-upgrade Arecibo imagery and reveals many additional ice features. Some medium-size craters located within 3° of the north pole show near-complete ice coverage over their floors, central peaks, and southern interior rim walls and little or no ice on their northern rim walls, while one large (90 km) crater at 85°N shows a sharp ice-cutoff line running across its central floor. All of this is consistent with the estimated polar extent of permanent shading from direct sunlight. Some craters show ice in regions that, though permanently shaded, should be too warm to maintain unprotected surface ice owing to indirect heating by reflected and reradiated sunlight. However, the ice distribution in these craters is in good agreement with models invoking insulation by a thin dust mantle. Comparisons with Goldstone X-band radar imagery indicate a wavelength dependence that could be consistent with such a dust mantle. More than a dozen small ice features have been found at latitudes between 67° and 75°. All of this low-latitude ice is probably sheltered in or under steep pole-facing crater rim walls, although, since most is located in the Mariner-unimaged hemisphere, confirmation must await imaging by the MESSENGER orbiter. These low-latitude features are concentrated toward the “cold longitudes,” possibly indicating a thermal segregation effect governed by indirect heating. The radar imagery places the corrected locations of the north and south poles at 7°W, 88.35°N and 90°W, 88.7°S, respectively, on the original Mariner-based maps.     

Librations and obliquity of Mercury from the BepiColombo radio-science and camera experiments

A major goal of the BepiColombo mission to Mercury is the determination of the structure and state of Mercury's interior. Here the BepiColombo rotation experiment has been simulated in order to assess the ability to attain the mission goals and to help lay out a series of constraints on the experiment's possible progress. In the rotation experiment pairs of images of identical surface regions taken at different epochs are used to retrieve information on Mercury's rotation and orientation. The idea is that from observations of the same patch of Mercury's surface at two different solar longitudes of Mercury the orientation of Mercury can be determined, and therefore also the obliquity and rotation variations with respect to the uniform rotation.The estimation of the libration amplitude and obliquity through pattern matching of observed surface landmarks is challenging. The main problem arises from the difficulty to observe the same landmark on the planetary surface repeatedly over the MPO mission lifetime, due to the combination of Mercury's 3:2 spin-orbit resonance, the absence of a drift of the MPO polar orbital plane and the need to combine data from different instruments with their own measurement restrictions.By assuming that Mercury occupies a Cassini state and that the spacecraft operates nominally we show that under worst case assumptions the annual libration amplitude and obliquity can be measured with a precision of, respectively, 1.4 arcseconds (as) and 1.0 as over the nominal BepiColombo MPO lifetime with about 25 landmarks for rather stringent illumination restrictions. The outcome of the experiment cannot be easily improved by simply relaxing the observational constraints, or increasing the data volume.     

Radius and limb topography of Mercury obtained from images acquired during the MESSENGER flybys

Analysis of images obtained by the MESSENGER spacecraft during its three flybys of Mercury yields a new estimate for the planet's mean radius of 2439.25±0.69 km, in agreement with results from Mariner 10 and Earth-based observations, as well as with MESSENGER altimeter and occultation data. The mean equatorial radius and polar radius are identical to within error, suggesting that rotational oblateness is negligible when compared with other sources of topography. This result is consistent with the small gravitational oblateness of the planet. Minor differences in radius obtained at different locations reflect regional variations in topography. Residual topography along three limb profiles has a dynamic range of 7.4 km and a root-mean-square roughness of 0.8 km over hemispherical scales. Following MESSENGER's entry into orbit about Mercury in March 2011, we expect considerable additional improvements to our knowledge of Mercury's size and shape.     

Mercury crater statistics from MESSENGER flybys: Implications for stratigraphy and resurfacing history

The primary crater population on Mercury has been modified by volcanism and secondary craters. Two phases of volcanism are recognized. One volcanic episode that produced widespread intercrater plains occurred during the period of the Late Heavy Bombardment and markedly altered the surface in many areas. The second episode is typified by the smooth plains interior and exterior to the Caloris basin, both of which have a different crater size-frequency distribution than the intercrater plains, consistent with a cratering record dominated by a younger population of impactors. These two phases may have overlapped as parts of a continuous period of volcanism during which the volcanic flux tended to decrease with time. The youngest age of smooth plains volcanism cannot yet be determined, but at least small expanses of plains are substantially younger than the plains associated with the Caloris basin. The spatial and temporal variations of volcanic resurfacing events can be used to reconstruct Mercury's geologic history from images and compositional and topographic data to be acquired during the orbital phase of the MESSENGER mission.     

Resonant forcing of Mercury's libration in longitude

The period of free libration of Mercury's longitude about the position it would have had if it were rotating uniformly at 1.5 times its orbital mean motion is close to resonance with Jupiter's orbital period. The Jupiter perturbations of Mercury's orbit thereby lead to amplitudes of libration at the 11.86 year period that may exceed the amplitude of the 88 day forced libration determined by radar. Mercury's libration in longitude may be thus dominated by only two periods of 88 days and 11.86 years, where other periods from the planetary perturbations of the orbit have much smaller amplitudes.     

Constraints on the composition of the martian south polar cap from gravity and topography

The polar caps of Mars have long been acknowledged to be composed of unknown proportions of water ice, solid CO₂ (dry ice), and dust. Gravity and topography data are here analyzed over the southern cap to place constraints on its density, and hence composition. Using a localized spectral analysis combined with a lithospheric flexure model of ice cap loading, the best fit density of the volatile-rich south polar layered deposits is found to be 1271 kg m⁻³ with 1-σ limits of 1166 and 1391 kg m⁻³. The best fit elastic thickness of this geologically young deposit is 140 km, though any value greater than 102 km can fit the observations. The best fit density implies that about 55% dry ice by volume could be sequestered in these deposits if they were completely dust free. Alternatively, if these deposits were completely free of solid CO₂, the dust content would be constrained to lie between about 14 and 28% by volume. The bulk thermal conductivity of the polar cap is not significantly affected by these maximum allowable concentrations of dust. However, even if a moderate quantity of solid CO₂ were present as horizontal layers, the bulk thermal conductivity of the polar cap would be significantly reduced. Reasonable estimates of the present day heat flow of Mars predict that dry ice beneath the thicker portions of the south polar cap would have melted. Depending on the quantity of solid CO₂ in these deposits today, it is even possible that water ice could melt where the cap is thickest. If independent estimates for either the dust or CO₂ content of the south polar cap could be obtained, and if radar sounding data could determine whether this polar cap is presently experiencing basal melting or not, it would be possible to use these observations to place tight constraints on the present day heat flow of Mars.