Photometric correction of Mercury's global color mosaic

During the third flyby of Mercury by the MESSENGER spacecraft, a dedicated disk-integrated photometric sequence was acquired with the wide-angle multispectral camera to observe Mercury's global photometric behavior in 11 spectral filters over as broad a range of phase angle as possible within the geometric constraints of the flyby. Extraction of disk-integrated measurements from images acquired during this sequence required careful accounting for scattered light and residual background effects. The photometric model fit to these measurements is shown to fit observed radiances at phase angles below 110°, possibly except where both solar incidence and emission angles are high (>70°). The complexity of the scattered light at wavelengths greater than 828 nm contributes to a less accurate photometric correction at these wavelengths. The model is used to correct the global imaging data set acquired at a variety of geometries to a common geometry of incidence angle=30°, emission angle=0°, and phase angle=30°, yielding a relatively seamless mosaic. The results here will be used to correct image mosaics of Mercury acquired in orbit.     

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.     

Seismic effects of the Caloris basin impact, Mercury

Striking geological features on Mercury's surface have been linked to tectonic disruption associated with the Caloris impact and have the potential to provide information on the interior structure of Mercury. The unusual disrupted terrain located directly at the antipode of the 1500-km-diameter Caloris basin could have plausibly formed as a consequence of focused seismic waves generated by the massive impact event. In this paper, we revisit the antipodal seismic focusing effects of the Caloris impact by developing physically consistent structure models for Mercury and parameterized seismic source models for the Caloris impact. If the focused seismic body waves caused the disrupted terrain, then the amplitudes of the waves and the areal extent of surface disruptions could be used for estimating the seismic energy imparted by the impact.In this study, we show that effects of direct body waves are small relative to those of focused guided waves. Two types of guided waves are produced by the Caloris impact. One is the conventional Rayleigh wave generated by the impact. The second is the mantle guided waves trapped between the core and the free surface. Mantle guided waves, not recognized in previous studies, may have played an important role in the creation of the disrupted terrain. We find that the early core state has only moderate effects on the antipodal response to the Caloris impact. The fact that the zone of predicted disruption for both fluid and solid core cases is smaller than the observed region of chaotic terrain suggests either that the antipodal response to the Caloris impact was modulated by the shallow structure of Mercury, or that the energy imparted by the impact was larger than those used in this study.     

On the chemistry of mantle and magmatic volatiles on Mercury

The surface of Mercury contains ancient volcanic features and signs of pyroclastic activity related to unknown magmatic volatiles. Here, the nature of possible magmatic volatiles (H, S, C, Cl, and N) is discussed in the contexts of formation and evolution of the planet, composition and redox state of its mantle, solubility in silicate melts, chemical mechanisms of magma degassing, and thermochemical equilibria in magma and volcanic gases. The low-FeO content in surface materials (<6 wt%) evaluated with remote sensing methods corresponds to less than 2.3 fO₂ log units below the iron-wüstite buffer. These redox conditions imply restricted involvement of hydrous species in nebular and accretion processes, and a limited loss of S, Cl, and N during formation and evolution of the planet. Reduced conditions correspond to high solubilities of these elements in magma and do not favor degassing. Major degassing and pyroclastic activity would require oxidation of melts in near-surface conditions. Low-pressure oxidation of graphite in moderately oxidized magmas causes formation of low-solubility CO. Decompression of reduced N-saturated melts involves oxidation of nitride melt complexes and could cause N₂ degassing. Putative assimilation of oxide (FeO, TiO₂, and SiO₂) rich rocks in magma chambers could have caused major degassing through oxidation of graphite and S-, Cl- and N-bearing melt complexes. However, crustal rocks may never have been oxidized, and sulfides, graphite, chlorides, and nitrides could remain in crystallized rocks. Chemical equilibrium models show that N₂, CO, S₂, CS₂, S₂Cl, Cl, Cl₂, and COS could be among the most abundant volcanic gases on Mercury. Though, these speciation models are restricted by uncertain redox conditions, unknown volatile content in magma, and the adopted chemical degassing mechanism.     

Simulations using terrestrial geological analogues to assess interpretability of potential geological features of the Hermean surface restituted by the STereo imaging Camera of the SIMBIOSYS package (BepiColombo mission)

The BepiColombo space mission is one of the European Space Agency's cornerstone projects; it is planned for launch in 2013 to study the planet Mercury. One of the imaging instruments of BepiColombo is a STereo Camera (STC), whose main scientific objective is the global stereo mapping of the entire surface of Mercury. STC will permit the generation of a Digital Terrain Model (DTM) of Mercury's surface, improving the interpretation of morphological features at different scales and clarifying the stratigraphic relationships between different geological units.To evaluate the effectiveness of the STC-derived DTM for geological purposes, a series of simulations has been performed to find out to what extent the errors expected in the DTM may prevent the correct classification and interpretation of geological features. To meet this objective, Earth analogues (a crater, a lava cone and an endogenous dome) of likely components of the Hermean surface, small enough to be near the detection limit of the STC, were selected and a photorealistic three-dimensional (3D) model of each feature was generated using Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) stereo images. Stereoscopic pairs of synthetic images of each feature were then generated from the 3D model at different locations along the BepiColombo orbit. For each stereo pair, the corresponding Hermean DTM was computed using image correlation and compared to the reference data to assess the loss of detail and interpretability. Results show that interpretation and quantitative analysis of simple craters morphologies and small volcanic features should be possible all along the periherm orbit arc. At the apoherm only the larger features can be unequivocally distinguished, but they will be reconstructed to a poor approximation.     

Exposure of spectrally distinct material by impact craters on Mercury: Implications for global stratigraphy

MESSENGER’s Mercury Dual Imaging System (MDIS) obtained multispectral images for more than 80% of the surface of Mercury during its first two flybys. Those images have confirmed that the surface of Mercury exhibits subtle color variations, some of which can be attributed to compositional differences. In many areas, impact craters are associated with material that is spectrally distinct from the surrounding surface. These deposits can be located on the crater floor, rim, wall, or central peak or in the ejecta deposit, and represent material that originally resided at depth and was subsequently excavated during the cratering process. The resulting craters make it possible to investigate the stratigraphy of Mercury’s upper crust. Studies of laboratory, terrestrial, and lunar craters provide a means to bound the depth of origin of spectrally distinct ejecta and central peak structures. Excavated red material (RM), with comparatively steep (red) spectral slope, and low-reflectance material (LRM) stand out prominently from the surrounding terrain in enhanced-color images because they are spectral end-members in Mercury’s compositional continuum. Newly imaged examples of RM were found to be spectrally similar to the relatively red, high-reflectance plains (HRP), suggesting that they may represent deposits of HRP-like material that were subsequently covered by a thin layer (∼1 km thick) of intermediate plains. In one area, craters with diameters ranging from 30 km to 130 km have excavated and incorporated RM into their rims, suggesting that the underlying RM layer may be several kilometers thick. LRM deposits are useful as stratigraphic markers, due to their unique spectral properties. Some RM and LRM were excavated by pre-Tolstojan basins, indicating a relatively old age (>4.0 Ga) for the original emplacement of these deposits. Detailed examination of several small areas on Mercury reveals the complex nature of the local stratigraphy, including the possible presence of buried volcanic plains, and supports sequential buildup of most of the upper ∼5 km of crust by volcanic flows with compositions spanning the range of material now visible on the surface, distributed heterogeneously across the planet. This emerging picture strongly suggests that the crust of Mercury is characterized by a much more substantial component of early volcanism than represented by the phase of mare emplacement on Earth’s Moon.     

The sodium exosphere of Mercury: Comparison between observations during Mercury's transit and model results

In this study we compare the sodium exosphere observations made by Schleicher et al. [Schleicher, H., and 4 colleagues, 2004. Astron. Astrophys. 425, 1119-1124] with the result of a detailed numerical simulation. The observations, made during the transit of Mercury across the solar disk on 7 May 2003, show a maximum of sodium emission near the polar regions, with north prevalence, and the presence of a dawn-dusk asymmetry. We interpret this distribution as the resulting effect of two combined processes: the solar wind proton precipitation causing chemical alteration of the surface, freeing the sodium atoms from their bounds in the crystalline structure on the surface, and the subsequent photon-stimulated and thermal desorption of the sodium atoms. While we find that the velocity distribution of photon desorbed sodium can explain the observed exosphere population, thermal desorption seems to play a minor role only causing a smearing at the locations where Na atoms are released on the dayside. The observed and simulated distributions agree very well with this hypothesis and indicate that the combination of the proposed processes is able to explain the observed features.     

Induced magnetic field effects at planet Mercury

At Mercury's surface external magnetic field contributions caused by magnetospheric current systems play a much more important role than at Earth. They are subjected to temporal variations and therefore will induce currents in the large conductive iron core. These currents give rise to an additional magnetic field superposing the planetary field. We present a model to estimate the size of the induced fields using a magnetospheric magnetic field model with time-varying magnetopause position. For the Hermean interior we assume a two-layer conductivity distribution. We found out that about half of the surface magnetic field is due to magnetospheric or induced currents. The induced fields achieve 7-12% of the mean surface magnetic intensity of the internal planetary field, depending on the core size. The magnetic field was also modeled for a satellite moving along a polar orbit in the Hermean magnetosphere, showing the importance of a careful separation of the magnetic field measurements.     

First measurements with the Physikalisches Institut Radiometric Experiment (PHIRE)

We have constructed an experiment to perform bidirectional reflectance distribution function (BRDF) measurements of laboratory samples, and have used the experiment to characterize a sample of JSC-1 lunar regolith simulant. Characterizations relied on in-plane BRDF measurements in visible and near-infrared (NIR) bandpasses. The optical properties of the simulant sample were found to be similar to those observed for bright, lunar highland regions. Reflectance models (Hapke 1981. Bidirectional reflectance spectroscopy 1. Theory. J. Geophys. Res. 86(B4), 3,039−3,054; 1984. Bidirectional reflectance spectroscopy 3. Correction for macroscopic roughness. Icarus 59, 41−59; 1986. Bidirectional reflectance spectroscopy 4. The extinction coefficient and the opposition effect. Icarus 67, 264−280; 2002. Bidirectional reflectance spectroscopy 5. The coherent backscatter opposition effect and anisotropic scattering. Icarus 157, 523−534) made excellent fits to fixed incidence angle, variable emission angle data sets. However, the models were not found to extrapolate well to fixed, near-zero phase angle data at varying incidence angles, and no solutions were found that provided simultaneous, high quality fits to the two types of data sets. Except for the single-scattering albedo, the best-fit parameters of the fixed incidence angle data were statistically the same in the visible and NIR. Correlations between the reflectance model parameters were systematically examined, and strong correlations were found between single-scattering albedo and the two two-stream Henyey-Greenstein scattering parameters and, to a lesser extent, the small-scale mean surface roughness.     

Evolution of Mercury's obliquity

Mercury has a near-zero obliquity, i.e. its spin axis is nearly perpendicular to its orbital plane. The value of the obliquity must be known precisely in order to constrain the size of the planet's core within the framework suggested by Peale [Peale, S.J., 1976. Nature 262, 765-766]. Rambaux and Bois [Rambaux, N., Bois, E., 2004. Astron. Astrophys. 413, 381-393] have suggested that Mercury's obliquity varies on thousand-year timescales due to planetary perturbations, potentially ruining the feasibility of Peale's experiment. We use a Hamiltonian approach (free of energy dissipation) to study the spin-orbit evolution of Mercury subject to secular planetary perturbations. We can reproduce an obliquity evolution similar to that of Rambaux and Bois [Rambaux, N., Bois, E., 2004. Astron. Astrophys. 413, 381-393] if we integrate the system with a set of initial conditions that differs from the Cassini state. However the thousand-year oscillations in the obliquity disappear if we use initial conditions corresponding to the equilibrium position of the Cassini state. This result indicates that planetary perturbations do not force short-period, large amplitude oscillations in the obliquity of Mercury. In the absence of excitation processes on short timescales, Mercury's obliquity will remain quasi-constant, suggesting that one of the important conditions for the success of Peale's experiment is realized. We show that interpretation of data obtained in support of this experiment will require a precise knowledge of the spin-orbit configuration, and we provide estimates for two of the critical parameters, the instantaneous Laplace plane orientation and the orbital precession rate from numerical fits to ephemeris data. Finally we provide geometrical relationships and a scheme for identifying the correct initial conditions required in numerical integrations involving a Cassini state configuration subject to planetary perturbations.