Radar measurements of Mercury’s north pole at 70 cm wavelength

We present radar imaging of Mercury using the Arecibo Observatory’s 70-cm wavelength radar system during the inferior conjunction of July 1999. At that time the sub-Earth latitude was ∼11°N and the highly reflective region at Mercury’s north pole that was first identified in radar images at the shorter wavelengths of 3.6 cm [Slade, M.A., Butler, B.J., Muhleman, D.O., 1992. Science 258, 635-640] and 13 cm [Harmon, J.K., Slade, M.A., 1992. Science 258, 640-643] was again clearly detected. The reflectivity averaged over a 75,000 km² region including the pole is similar to that measured at the other wavelengths over a comparable area, and the 70 cm circular polarization ratio of μ_C∼0.87 is possibly slightly lower. If this strong backscattering results from volume scattering in low absorption layers, the persistence of this effect over more than an order of magnitude change in wavelength scale has implications for the depth and thickness of the deposits responsible. The resolution of the radar maps at this wavelength is not sufficient to resolve individual craters, nor to discern features at other latitudes, but the planet’s total reflectivity is consistent with previous work and the scattering function suggests a surface roughness at this wavelength similar to the lunar highlands.     

Radar imagery of the southern Caloris region, Mercury

Radar imagery from July 2005 Arecibo observations has provided new information on surface relief over the southern portion of Caloris Basin and the smooth plains to the south of the basin. A lobe of smooth plains has been identified in the Mariner-unimaged region southwest of Mozart Crater that coincides precisely with topographically down-bowed terrain seen in earlier Arecibo radar altimetry. A 105-km-diameter crater has been found at 193.6° W, 25.6° N that appears to be the largest crater in the Caloris basin floor.     

The influence of Mercury’s inner core on its physical libration

As Mercury orbits the Sun, gravitational torques on its equatorial elliptical shape give rise to a planetary libration. The amplitude of Mercury’s libration, as determined from Earth-based radar speckle pattern observations, suggests that only the mantle participates in the motion. This indicates a decoupling between the core and the mantle, and therefore that the outermost part of the core must be fluid. If a solid inner core is present at the center of Mercury, the equatorial elliptical shape of the latter may become misaligned with that of Mercury’s mantle, leading to an internal gravitational torque between the two. If this torque is large, it may participate in the dynamics of Mercury’s libration. The goal of this work is to determine whether Mercury’s observed librations can be used to place constraints on the properties of its inner core. We present a comparison between predicted and observed librations for a range of interior models of Mercury, with various inner core sizes and fluid core densities. We show that a marginally better fit to observations can be achieved for interior models that have an inner core radius larger than 400 km. However, the improvement in fit is small, and it is not possible to draw robust conclusions on the size of Mercury’s inner core on the basis of existing libration data. Nevertheless, our study demonstrates that the influence of the inner core on the libration of Mercury could be detected with a decade worth of accurate observations.     

Observations of metallic species in Mercury’s exosphere

From observations of the metallic species sodium (Na), potassium (K), and magnesium (Mg) in Mercury’s exosphere, we derive implications for source and loss processes. All metallic species observed exhibit a distribution and/or line width characteristic of high to extreme temperature - tens of thousands of degrees K. The temperatures of refractory species, including magnesium and calcium, indicate that the source process for the atoms observed in the tail and near-planet exosphere are consistent with ion sputtering and/or impact vaporization of a molecule with subsequent dissociation into the atomic form. The extended Mg tail is consistent with a surface abundance of 5-8% Mg by number, if 30% of impact-vaporized Mg remains as MgO and half of the impact vapor condenses. Globally, ion sputtering is not a major source of Mg, but locally the sputtered source can be larger than the impact vapor source. We conclude that the Na and K in Mercury’s exosphere can be derived from a regolith composition similar to that of Luna 16 soil (or Apollo 17 orange glass), in which the abundance by number is 0.0027 (0.0028) for Na and 0.0006 (0.0045) for K.     

Beagle Rupes - Evidence for a basal decollement of regional extent in Mercury’s lithosphere

Thanks to its location at low latitude and close to the terminator in the outbound view of Mercury obtained during MESSENGER’s first fly-by, the Beagle Rupes lobate scarp on Mercury has been particularly clearly imaged. This enables us to interpret it as a component of a linked fault system, consisting of a frontal scarp terminated by transpressive lateral ramps. The terrain bounded by these surface manifestations of faulting is the hanging-wall block of a thrust sheet and must be underlain by a basal decollement (a detachment horizon) constituting the fault zone at depth. The decollement must extend a minimum of 150 km eastwards from the frontal scarp, and at least 400 km if displacement is transferred to features interpreted as out-of-sequence thrusts and offset lateral ramps that appear to continue the linked fault system to the east. The depth of the basal decollement could be controlled by crustal stratigraphy or by rheological change within, or at the base of, the lithosphere. Previous interpretations of mercurian lobate scarps regard their thrusts as uniformly dipping and dying out at depth, lacking lateral ramps and any extensive detachment horizon. Anticipated improvements in image resolution and lighting geometry should make it possible to document what percentage of lobate scarps share the Beagle Rupes style of tectonics.     

A radar survey of main-belt asteroids: Arecibo observations of 55 objects during 1999-2003

We report Arecibo observations of 55 main-belt asteroids (MBAs) during 1999-2003. Most of our targets had not been detected previously with radar, so these observations more than double the number of radar-detected MBAs. Our bandwidth estimates constrain our targets' pole directions in a manner that is geometrically distinct from optically derived constraints. We present detailed statistical analyses of the disk-integrated properties (radar albedo and circular polarization ratio) of the 84 MBAs observed with radar through March 2003; all of these observations are summarized in the online supplementary information. Certain conclusions reached in previous studies are strengthened: M asteroids have higher mean radar albedos and a wider range of albedos than do other MBAs, suggesting that both metal-rich and metal-poor M-class objects exist; and C- and S-class MBAs have indistinguishable radar albedo distributions, suggesting that most S-class objects are chondritic. Also in accord with earlier results, there is evidence that primitive asteroids from outside the C taxon (F, G, P, and D) are not as radar-bright as C and S objects, but a convincing statistical test must await larger sample sizes. In contrast with earlier work, we find S-class MBAs to have higher circular polarization ratios than other MBAs, indicating greater near-surface structural complexity at decimeter scales, due to different mineralogy (material strength or loss tangent), a different impactor population, or both.     

Mercury: Radar images of the equatorial and midlatitude zones

Radar imaging results for Mercury's non-polar regions are presented. The dual-polarization, delay-Doppler images were obtained from several years of observations with the upgraded Arecibo S-band (λ12.6-cm) radar telescope. The images are dominated by radar-bright features associated with fresh impact craters. As was found from earlier Goldstone-VLA and pre-upgrade Arecibo imaging, three of the most prominent crater features are located in the Mariner-unimaged hemisphere. These are: “A,” an 85-km-diameter crater (348° W, 34° S) whose radar ray system may be the most spectacular in the Solar System; “B,” a 95-km-diameter crater (343° W, 58° N) with a very bright halo but less distinct ray system; and “C,” an irregular feature with bright ejecta and rays distributed asymmetrically about a 125-km source crater (246° W, 11° N). Due south of “C” lies a “ghost” feature (242° W, 27° S) that resembles “A” but is much fainter. An even fainter such feature is associated with Bartok Crater. These may be two of the best mercurian examples of large ejecta/ray systems observed in an intermediate state of degradation. Virtually all of the bright rayed craters in the Mariner 10 images show radar rays and/or bright rim rings, with radar rays being less common than optical rays. Radar-bright craters are particularly common in the H-7 quadrangle. Some diffuse radar albedo variations are seen that have no obvious association with impact ejecta. In particular, some smooth plains regions such as the circum-Caloris plains in Tir, Budh, and Sobkou Planitiae and the interiors of Tolstoj and “Skinakas” basins show high depolarized brightness relative to their surroundings, which is the reverse of the mare/highlands contrast seen in lunar radar images. Caloris Basin, on the other hand, appears dark and featureless in the images.     

Monte Carlo modeling of sodium in Mercury’s exosphere during the first two MESSENGER flybys

We present a Monte Carlo model of the distribution of neutral sodium in Mercury’s exosphere and tail using data from the Mercury Atmospheric and Surface Composition Spectrometer (MASCS) on the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft during the first two flybys of the planet in January and September 2008. We show that the dominant source mechanism for ejecting sodium from the surface is photon-stimulated desorption (PSD) and that the desorption rate is limited by the diffusion rate of sodium from the interior of grains in the regolith to the topmost few monolayers where PSD is effective. In the absence of ion precipitation, we find that the sodium source rate is limited to ∼10⁶-10⁷ cm⁻² s⁻¹, depending on the sticking efficiency of exospheric sodium that returns to the surface. The diffusion rate must be at least a factor of 5 higher in regions of ion precipitation to explain the MASCS observations during the second MESSENGER flyby. We estimate that impact vaporization of micrometeoroids may provide up to 15% of the total sodium source rate in the regions observed. Although sputtering by precipitating ions was found not to be a significant source of sodium during the MESSENGER flybys, ion precipitation is responsible for increasing the source rate at high latitudes through ion-enhanced diffusion.     

Transient surface liquid in Titan’s polar regions from Cassini

Cassini RADAR images of Titan’s south polar region acquired during southern summer contain lake features which disappear between observations. These features show a tenfold increases in backscatter cross-section between images acquired one year apart, which is inconsistent with common scattering models without invoking temporal variability. The morphologic boundaries are transient, further supporting changes in lake level. These observations are consistent with the exposure of diffusely scattering lakebeds that were previously hidden by an attenuating liquid medium. We use a two-layer model to explain backscatter variations and estimate a drop in liquid depth of approximately 1-m-per-year. On larger scales, we observe shoreline recession between ISS and RADAR images of Ontario Lacus, the largest lake in Titan’s south polar region. The recession, occurring between June 2005 and July 2009, is inversely proportional to slopes estimated from altimetric profiles and the exponential decay of near-shore backscatter, consistent with a uniform reduction of 4 ± 1.3 m in lake depth.Of the potential explanations for observed surface changes, we favor evaporation and infiltration. The disappearance of dark features and the recession of Ontario’s shoreline represents volatile transport in an active methane-based hydrologic cycle. Observed loss rates are compared and shown to be consistent with available global circulation models. To date, no unambiguous changes in lake level have been observed between repeat images in the north polar region, although further investigation is warranted. These observations constrain volatile flux rates in Titan’s hydrologic system and demonstrate that the surface plays an active role in its evolution. Constraining these seasonal changes represents the first step toward our understanding of longer climate cycles that may determine liquid distribution on Titan over orbital time periods.     

Constraints on Mercury’s surface composition from MESSENGER and ground-based spectroscopy

The composition and chemistry of Mercury’s regolith has been calculated from MESSENGER MASCS 0.3-1.3 μm spectra from the first flyby, using an implementation of Hapke’s radiative transfer-based photometric model for light scattering in semi-transparent porous media, and a linear spectral mixing algorithm. We combine this investigation with linear spectral fitting results from mid-infrared spectra and compare derived oxide abundances with mercurian formation models and lunar samples. Hapke modeling results indicate a regolith that is optically dominated by finely comminuted particles with average area weighted grain size near 20 μm. Mercury shows lunar-style space weathering, with maturation-produced microphase iron present at ∼0.065 wt.% abundance, with only small variations between mature and immature sites, the amount of which is unable to explain Mercury’s low brightness relative to the Moon. The average modal mineralogies for the flyby 1 spectra derived from Hapke modeling are 35-70% Na-rich plagioclase or orthoclase, up to 30% Mg-rich clinopyroxene, <5% Mg-rich orthopyroxene, minute olivine, ∼20-45% low-Fe, low-Ti agglutinitic glass, and <10% of one or more lunar-like opaque minerals. Mercurian average oxide abundances derived from Hapke models and mid-infrared linear fitting include 40-50 wt.% SiO₂, 10-35 wt.% Al₂O₃, 1-8 wt.% FeO, and <25 wt.% TiO₂; the inferred rock type is basalt. Lunar-like opaques or glasses with high Fe and/or Ti abundances cannot on their own, or in combination, explain Mercury’s low brightness. The linear mixing results indicate the presence of clinopyroxenes that contain up to 21 wt.% MnO and the presence of a Mn-rich hedenbergite. Mn in M1 crystalline lattice sites of hedenbergite suppresses the strong 1 and 2 μm crystal field absorption bands and may thus act as a strong darkening agent on Mercury. Also, one or more of thermally darkened silicates, Fe-poor opaques and matured glasses, or Mercury-unique Ostwald-ripened microphase iron nickel may lower the albedo. A major part of the total microphase iron present in Mercury’s regolith is likely derived from FeO that is not intrinsic to the crust but has been subsequently delivered by exogenic sources.     

HiRISE observations of gas sublimation-driven activity in Mars’ southern polar regions: III. Models of processes involving translucent ice

Enigmatic surface features, known as ‘spiders’, found at high southern martian latitudes, are probably caused by sublimation-driven erosion under the seasonal carbon dioxide ice cap. The Mars Reconnaissance Orbiter (MRO) High Resolution Imaging Science Experiment (HiRISE) has imaged this terrain in unprecedented details throughout southern spring. It has been postulated [Kieffer, H.H., Titus, T.N., Mullins, K.F., Christensen, P.R., 2000. J. Geophys. Res. 105, 9653-9700] that translucent CO₂ slab ice traps gas sublimating at the ice surface boundary. Wherever the pressure is released the escaping gas jet entrains loose surface material and carries it to the top of the ice where it is carried downslope and/or downwind and deposited in a fan shape. Here we model two stages of this scenario: first, the cleaning of CO₂ slab ice from dust, and then, the breaking of the slab ice plate under the pressure built below it by subliming ice. Our modeling results and analysis of HiRISE images support the gas jet hypothesis and show that outbursts happen very early in spring.     

Orbital effects on Mercury’s escaping sodium exosphere

We present results from coronagraphic imaging of Mercury’s sodium tail over a 7° field of view. Several sets of observations made at the McDonald Observatory since May 2007 show a tail of neutral sodium atoms stretching more than 1000 Mercury radii (R_m) in length, or a full degree of sky. However, no tail was observed extending beyond 120 R_m during the January 2008 MESSENGER fly-by period, or during a similar orbital phase of Mercury in July 2008. Large changes in Mercury’s heliocentric radial velocity cause Doppler shifts about the Fraunhofer absorption features; the resultant change in solar flux and radiation pressure is the primary cause of the observed variation in tail brightness. Smaller fluctuations in brightness may exist due to changing source rates at the surface, but we have no explicit evidence for such changes in this data set. The effects of radiation pressure on Mercury’s escaping atmosphere are investigated using seven observations spanning different orbital phases. Total escape rates of atmospheric sodium are estimated to be between 5 and 13 × 10²³ atoms/s and show a correlation to radiation pressure. Candidate sources of Mercury’s sodium exosphere include desorption by UV sunlight, thermal desorption, solar wind channeled along Mercury’s magnetic field lines, and micro-meteor impacts. Wide-angle observations of the full extent of Mercury’s sodium tail offer opportunities to enhance our understanding of the time histories of these source rates.     

New Cassini RADAR results for Saturn’s icy satellites

Cassini radar tracks on Saturn’s icy satellites through the end of the Prime Mission in 2008 have increased the number of radar albedo estimates from 10 (Ostro et al., 2006) to 73. The measurements sample diverse subradar locations (and for Dione, Rhea, and Iapetus almost always use beamwidths less than half the target angular diameters), thereby constraining the satellites’ global radar albedo distributions. The echoes result predominantly from volume scattering, and their strength is thus strongly sensitive to ice purity and regolith maturity. The combination of the Cassini data set and Arecibo 13-cm observations of Enceladus, Tethys, Dione, Rhea (Black et al., 2007), and Iapetus (Black et al., 2004) discloses an unexpectedly complex pattern of 13-to-2-cm wavelength dependence. The 13-cm albedos are generally smaller than 2-cm albedos and lack the correlation seen between 2-cm and optical geometric albedos. Enceladus and Iapetus are the most interesting cases. We infer from hemispheric albedo variations that the E-ring has a prominent effect on the 13-cm radar “lightcurve”. The uppermost trailing-side regolith is too fresh for meteoroid bombardment to have developed larger-scale heterogeneities that would be necessary to elevate the 13-cm radar albedo, whereas all of Enceladus is clean and mature enough for the 2-cm albedo to be uniformly high. For, Iapetus, the 2-cm albedo is strongly correlated with optical albedo: low for the optically dark, leading-side material and high for the optically bright, trailing-side material. However, Iapetus’ 13-cm albedo values show no significant albedo dichotomy and are several times lower than 2-cm values, being indistinguishable from the weighted mean of 13-cm albedos for main-belt asteroids, 0.15 ± 0.10. The leading side’s optically dark contaminant must be present to depths of at least one to several decimeters, so 2-cm albedos can mimic the optical dichotomy; however, it does not have to extend any deeper than that. The fact that both hemispheres of Iapetus look Asteroid-like at 13 cm means that coherent backscattering itself is not nearly as effective as it is at 2 cm. Since Iapetus’ entire surface is mature regolith, the wavelength dependence must involve composition, not structure. Either the composition is a function of depth everywhere (with electrical loss much greater at depths greater than a decimeter or two), or the intrinsic electrical loss of some pervasive constituent is much higher at 13 cm than at 2 cm. Ammonia is a candidate for such a contaminant. If ammonia’s electrical properties do not depend on frequency, and if ammonia is globally much less abundant within the upper one or two decimeters than at greater depths, then coherent backscattering would effectively be shut down at 13 cm, explaining the Asteroid-like 13-cm albedo.     

Arecibo radar observations of Phobos and Deimos

In November 2005, we observed the moons of Mars using the Arecibo 2380-MHz (13-cm) radar, obtaining a result for the OC radar albedo of Phobos (0.056±0.014) consistent with its previously reported radar albedo and implying an upper bound on its near-surface bulk density of . We detected Deimos by radar for the first time, finding its OC radar albedo to be 0.021±0.006, implying an upper bound on its near-surface density of , consistent with a high-porosity regolith. We briefly discuss reasons for these low radar albedos, Deimos' being possibly the lowest of any Solar System body yet observed by radar.     

A method to infer past surface mass balance and topography from internal layers in martian polar layered deposits

Internal layers in ice masses can be detected with ice-penetrating radar. In a flowing ice mass, each horizon represents a past surface that has been subsequently buried by accumulation, and strained by ice flow. These layers retain information about relative spatial patterns of accumulation and ablation (mass balance). Internal layers are necessary to accurately infer mass-balance patterns because the ice-surface shape only weakly reflects spatial variations in mass balance. Additional rate-controlling information, such as the layer age, the ice temperature, or the ice-grain sizes and ice-crystal fabric, can be used to infer the absolute rate of mass balance. To infer mass balance from the shapes of internal layers, we solve an inverse problem. The solution to the inverse problem is the best set or sets of unknown boundary conditions or initial conditions that, when used in our calculation of ice-surface elevation and internal-layer shape, generate appropriate predictions of observations that are available. We also show that internal layers can be used to infer martian paleo-surface topography from a past era of ice flow, even though the topography may have been largely altered by subsequent erosion. We have successfully inferred accumulation rates and surface topography from internal layers in Antarctica. Using synthetic data, we demonstrate the ability of this method to solve the corresponding inverse problem to infer accumulation and ablation rates, as well as the surface topography, for martian ice. If past ice flow has affected the shapes of martian internal layers, this method is necessary to infer the spatial pattern and rate of mass balance.     

Arecibo radar observations of Rhea, Dione, Tethys, and Enceladus

We have measured the bulk radar reflectance properties of the mid-size saturnian satellites Rhea, Dione, Tethys, and Enceladus with the Arecibo Observatory's 13 cm wavelength radar system during the 2004 through 2007 oppositions of the Saturn system. Comparing to the better studied icy Galilean satellites, we find that the total reflectivities of Rhea and Tethys are most similar to Ganymede while Dione is most similar to Callisto. Enceladus' reflectivity falls between those of Ganymede and Europa. The mean circular polarization ratios of the saturnian satellites range from ∼0.8 to 1.2, and are on average lower than those of the icy Galilean satellites at this wavelength although still larger than expected for single reflections off the surface. The ratio for the trailing hemisphere of Enceladus may be the exception with a value ?0.56. The 13 cm wavelength radar albedos and polarization ratios may be systematically lower than similar results from the Cassini orbiter's RADAR instrument at 2.2 cm wavelength [Ostro, S.J., and 19 colleagues, 2006. Icarus 183, 479-490]. Overall, these reflectivities and polarization properties, together with the shapes of the echo spectra, suggest subsurface multiple scattering to be the dominant reflection mechanism although operating less efficiently than on the large icy moons of Jupiter. All these saturnian moons and icy jovian moons are atmosphere-less, low temperature water ice surfaces, and any differences in radar properties may be indicative of differences in composition or the effects of various processes that modify the regolith structure. The degree of variation in radar properties with wavelength on each satellite may constrain the thickness and efficiency of the scattering layer.     

Radar observations of Comet P/2005 JQ5 (Catalina)

Arecibo radar observations of Comet P/2005 JQ5 (Catalina) have produced the first delay-Doppler images of a comet nucleus and the first radar detection of large-grain ejection from a Jupiter-family comet. The nucleus is small (1.4 km diameter), rough, and rapidly rotating. The large (>cm) grains have low velocities (∼1 m/s) and a low production rate.     

Identification and measurement of neutron-absorbing elements on Mercury’s surface

MESSENGER Neutron Spectrometer (NS) observations of cosmic-ray-generated thermal neutrons provide the first direct measurements of Mercury’s surface elemental composition. Specifically, we show that Mercury’s surface is enriched in neutron-absorbing elements and has a measured macroscopic neutron-absorption cross section of 45-81 × 10⁻⁴ cm²/g, a range similar to the neutron absorption of lunar basalts from Mare Crisium. The expected neutron-absorbing elements are Fe and Ti, with possible trace amounts of Gd and Sm. Fe and Ti, in particular, are important for understanding Mercury’s formation and how its surface may have changed over time through magmatic processes. With neutron Doppler filtering - a neutron energy separation technique based on spacecraft velocity - we demonstrate that Mercury’s surface composition cannot be matched by prior models, which have characteristically low abundances of Fe, Ti, Gd, and Sm. While neutron spectroscopy alone cannot separate the relative contributions of individual neutron-absorbing elements, these results provide strong new constraints on the nature of Mercury’s surface materials. For example, if all the measured neutron absorption were due to the presence of an Fe-Ti oxide and that oxide were ilmenite, then Mercury’s surface would have an ilmenite content of 7-18 wt.%. This result is in general agreement with the inference from color imaging and visible-near-infrared spectroscopy that Mercury’s overall low reflectance is consistent with a surface composition that is enriched in Fe-Ti oxides. The incorporation of substantial Fe and Ti in oxides would imply that the oxygen fugacity of basalts on Mercury is at the upper range of oxygen fugacities inferred for basalts on the Moon.     

Constraints on Mercury’s Na exosphere: Combined MESSENGER and ground-based data

We have used observations of sodium emission obtained with the McMath-Pierce solar telescope and MESSENGER’s Mercury Atmospheric and Surface Composition Spectrometer (MASCS) to constrain models of Mercury’s sodium exosphere. The distribution of sodium in Mercury’s exosphere during the period January 12-15, 2008, was mapped using the McMath-Pierce solar telescope with the 5″ × 5″ image slicer to observe the D-line emission. On January 14, 2008, the Ultraviolet and Visible Spectrometer (UVVS) channel on MASCS sampled the sodium in Mercury’s anti-sunward tail region. We find that the bound exosphere has an equivalent temperature of 900-1200 K, and that this temperature can be achieved if the sodium is ejected either by photon-stimulated desorption (PSD) with a 1200 K Maxwellian velocity distribution, or by thermal accommodation of a hotter source. We were not able to discriminate between the two assumed velocity distributions of the ejected particles for the PSD, but the velocity distributions require different values of the thermal accommodation coefficient and result in different upper limits on impact vaporization. We were able to place a strong constraint on the impact vaporization rate that results in the release of neutral Na atoms with an upper limit of 2.1 × 10⁶ cm⁻² s⁻¹. The variability of the week-long ground-based observations can be explained by variations in the sources, including both PSD and ion-enhanced PSD, as well as possible temporal enhancements in meteoroid vaporization. Knowledge of both dayside and anti-sunward tail morphologies and radiances are necessary to correctly deduce the exospheric source rates, processes, velocity distribution, and surface interaction.