Mercury surface composition: Integrating petrologic modeling and remote sensing data to place constraints on FeO abundance

A significant opaque component in Mercury’s crust is inferred based on albedo and spectral observations. Previous workers have favored iron-titanium bearing oxide minerals as the spectrally neutral opaque. A consequence of this hypothesis is that Mercury’s surface would have a high FeO content. An array of remote sensing techniques have not provided definitive constraints on the FeO content of Mercury’s surface. However, spectral observations have not detected a diagnostic 1 μm absorption band and have thus limited the FeO in coexisting silicates to <2 wt.% FeO. In this paper, we assess equilibrium among oxide and silicate minerals to constrain the distribution of iron between opaque oxides and silicates under a variety of environmental conditions. Equilibrium modeling is favored here because the geologic process that produced Mercury’s low-albedo intermediate terrain must have occurred globally, which favors a common widespread igneous process. Based on our modeling, we find that iron-rich ilmenite cannot occur with silicates that do not display a 1 μm absorption feature unless plagioclase abundances are high. However, such high plagioclase abundances are precluded by Mercury’s low albedo. Incorporating equilibrium crystallization modeling with spectral and albedo constraints we find the iron abundance of Mercury’s intermediate terrain is ?10 wt.% FeO. This intermediate iron composition matches constraints provided by visible albedo and total neutron absorption observed by MESSENGER. In fact, the total neutron absorption of mixtures of oxide, plagioclase, olivine and pyroxene for the oxide abundances estimated for Mercury, favor Mg-rich members of the ilmenite-geikielite solid-solution series. This work offers compositional constraints for Fe, Ti, and Mg that will be testable by various MESSENGER instrument data sets after it begins its orbital mission.     

Possible impact melt and debris flows at Tooting Crater, Mars

Candidate examples of impact melt flows and debris flows have been identified at Tooting crater, an extremely young (<2 Myr), 29 km diameter impact crater in Amazonis Planitia, Mars. Using HiRISE and CTX images, and stereo-derived digital elevation models derived from these images, we have studied the rim and interior wall of Tooting crater to document the morphology and topography of several flow features in order to constrain the potential flow formation mechanisms. Four flow types have been identified; including possible impact melt sheets and three types of debris flows. The flow features are all located within 2 km of the rim crest on the southern rim or lie on the southern interior wall of the crater ∼1500 m below the rim crest. Extensive structural failure has modified the northern half of the crater inner wall and we interpret this to have resulted in the destruction of any impact melt emplaced, as well as volatile-rich wall rock. The impact melt flows are fractured on the meter to decameter scale, have ridged, leveed lobes and flow fronts, and cover an area >6 km × 5 km on the southern rim. The debris flows are found on both the inner wall and rim of the crater, are ∼1-2 km in length, and vary from a few tens of meters to >300 m in width. These flows exhibit varying morphologies, from a channelized, leveed flow with arcuate ridges in the channel, to a rubbly flow with a central channel but no obvious levees. The flows indicate that water existed within the target rocks at the time of crater formation, and that both melt and fluidized sediment was generated during this event.     

Super-resolution of THEMIS thermal infrared data: Compositional relationships of surface units below the 100 meter scale on Mars

The Thermal Emission Imaging System (THEMIS) has provided the highest spatial resolution (100 m/pixel) thermal infrared (TIR) data of the surface of Mars to date. These data have enabled the discovery of many small-scale compositional units and helped to better constrain surface processes operating at these scales. However, with higher-resolution visible (VIS) instruments revealing smaller-scale surficial differences, there exists a need to detect and map compositional variability using TIR data at scales below 100 m. Because it is unlikely there will be a higher-resolution TIR instrument sent to Mars in the near future, creative image processing techniques commonly classified under the umbrella of “super-resolution” can be employed to improve or enhance the spatial resolution of the THEMIS TIR data. These approaches typically integrate another higher-resolution dataset and can either be qualitative for visual appeal, quantitative for data accuracy, or some combination of both. The super-resolution approach presented here produces enhanced TIR images that are radiometrically accurate, but also visually appealing. For the technique to be successfully applied, multi- to hyper-spectral data from two different spectral regions are required (e.g., the THEMIS TIR and VIS datasets). The focus here is to introduce this new super-resolution methodology and demonstrate its ability using existing THEMIS IR and VIS data. The quartzofeldspathic deposit in northern Syrtis Major was selected because of the spectral variability detected using the original IR resolution data and to better constrain the relationship between the small-scale surface morphology and areal extent of the deposit as well as its formation process. Despite being associated with the central peaks of two craters, the results here show no positive correlation between the small rocky outcrops and the quartzofeldspathic unit. A gradational contact exists between the unit and basaltic sands within the intercrater eolian material. The super-resolution approach offers an alternative approach to traditional sub-pixel deconvolution identification and provides a higher-resolution IR dataset for thermophysical and spectral analysis on Mars.     

The initial responses of hot liquid water released under low atmospheric pressures: Experimental insights

Experiments have been performed to simulate the shallow ascent and surface release of water and brines under low atmospheric pressure. Atmospheric pressure was treated as an independent variable and water temperature and vapor pressure were examined as a function of total pressure variation down to low pressures. The physical and thermal responses of water to reducing pressure were monitored with pressure transducers, temperature sensors and visible imaging. Data were obtained for pure water and for solutions with dissolved NaCl or CO₂. The experiments showed the pressure conditions under which the water remained liquid, underwent a rapid phase change to the gas state by boiling, and then solidified because of removal of latent heat. Liquid water is removed from phase equilibrium by decompression. Solid, liquid and gaseous water are present simultaneously, and not at the 611 Pa triple point, because dynamic interactions between the phases maintain unstable temperature gradients. After phase changes stop, the system reverts to equilibrium with its surroundings. Surface and shallow subsurface pressure conditions were simulated for Mars and the icy satellites of the outer Solar System. Freezing by evaporation in the absence of wind on Mars is shown to be unlikely for pure water at pressures greater than c. 670 Pa, and for saline solutions at pressures greater than c. 610 Pa. The physical nature of ice that forms depends on the salt content. Ice formed from saline water at pressures less than c. 610 Pa could be similar to terrestrial sea ice. Ice formed from pure water at pressures less than c. 100 Pa develops a low thermal conductivity and a ‘honeycomb’ structure created by sublimation. This ice could have a density as low as c. 450 kg m⁻³ and a thermal conductivity as low as 1.6 W m⁻¹ K⁻¹, and is highly reflective, more akin to snow than the clear ice from which it grew. The physical properties of ice formed from either pure or saline water at low pressures will act to reduce the surface temperature, and hence rate of sublimation, thereby prolonging the lifespan of any liquid water beneath.     

Evaporation effects on the formation of martian gullies

In order to investigate the formation of martian gullies and the stability of fluids on Mars, we examined about 120 gully images. Twelve HiRISE images contained a sufficient number of Transverse Aeolian Ridges (TARs) associated with the gullies to make the following measurements: overall gully length, length of the alcove, channel and apron, and we also measured the frequency of nearby TARs. Six of the 12 images examined showed a statistically significant negative correlation between overall gully length (alcove, channel and apron length) and TAR frequency. Previous experimental work from our group has shown that at temperatures below ∼200 K, evaporation rate increases by about an order of magnitude as wind speed increases from 0 to ∼15 m/s. Thus the negative correlations we observe between gully length and dune frequency can be explained by formation at temperatures below ∼200 K where wind speed/evaporation is a factor governing gully length. In these cases evaporation of the fluid carving the gully was a constraint on their dimensions. Cases where there is no correlation between gully length and TAR frequency, can be explained by formation at temperatures >200 K. The temperatures are consistent with Global Circulation Model and Thermal Emission Spectrometer (TES) data for these latitudes. The temperatures suggested by these trends are consistent with the fluid responsible for gully formation being a strong brine, such as Fe₂(SO₄)₃ which has a eutectic temperature of ∼200 K. We also find that formation timescales for gullies are 10⁵-10⁶ years.     

Flow patterns of lobate debris aprons and lineated valley fill north of Ismeniae Fossae, Mars: Evidence for extensive mid-latitude glaciation in the Late Amazonian

A variety of Late Amazonian landforms on Mars have been attributed to the dynamics of ice-related processes. Evidence for large-scale, mid-latitude glacial episodes existing within the last 100 million to 1 billion years on Mars has been presented from analyses of lobate debris aprons (LDA) and lineated valley fill (LVF) in the northern and southern mid-latitudes. We test the glacial hypothesis for LDA and LVF along the dichotomy boundary in the northern mid-latitudes by examining the morphological characteristics of LDA and LVF surrounding two large plateaus, proximal massifs, and the dichotomy boundary escarpment north of Ismeniae Fossae (centered at 45.3°N and 39.2°E). Lineations and flow directions within LDA and LVF were mapped using images from the Context (CTX) camera, the Thermal Emission Imaging Spectrometer (THEMIS), and the High Resolution Stereo Camera (HRSC). Flow directions were then compared to topographic contours derived from the Mars Orbiter Laser Altimeter (MOLA) to determine the down-gradient components of LDA and LVF flow. Observations indicate that flow patterns emerge from numerous alcoves within the plateau walls, are integrated over distances of up to tens of kilometers, and have down-gradient flow directions. Smaller lobes confined within alcoves and superposed on the main LDA and LVF represent a later, less extensive glacial phase. Crater size-frequency distributions of LDA and LVF suggest a minimum (youngest) age of 100 Ma. The presence of ring-mold crater morphologies is suggestive that LDA and LVF are formed of near-surface ice-rich bodies. From these observations, we interpret LDA and LVF within our study region to result from formerly active debris-covered glacial flow, consistent with similar observations in the northern mid-latitudes of Mars. Glacial flow was likely initiated from the accumulation and compaction of snow and ice on plateaus and in alcoves within the plateau walls as volatiles were mobilized to the mid-latitudes during higher obliquity excursions. Together with similar analyses elsewhere along the dichotomy boundary, these observations suggest that multiple glacial episodes occurred in the Late Amazonian and that LDA and LVF represent significant reservoirs of non-polar ice sequestered below a surface lag for hundreds of millions of years.     

Geologic mapping of the Hi’iaka and Shamshu regions of Io

We produced regional geologic maps of the Hi’iaka and Shamshu regions of Io’s antijovian hemisphere using Galileo mission data to assess the geologic processes that are involved in the formation of Io’s mountains and volcanic centers. Observations reveal that these regions are characterized by several types of volcanic activity and features whose orientation and texture indicate tectonic activity. Among the volcanic features are multiple hotspots and volcanic vents detected by Galileo, one at each of the major paterae: Hi’iaka, Shamshu, and Tawhaki. We mapped four primary types of geologic units: flows, paterae floors, plains, and mountains. The flows and patera floors are similar, but are subdivided based upon emplacement environments and mechanisms. The floors of Hi’iaka and Shamshu Paterae have been partially resurfaced by dark lava flows, although portions of the paterae floors appear bright and unchanged during the Galileo mission; this suggests that the floors did not undergo complete resurfacing as flooding lava lakes. However, the paterae do contain compound lava flow fields and show the greatest activity near the paterae walls, a characteristic of Pele type lava lakes. Mountain materials are tilted crustal blocks that exhibit varied degrees of degradation. Lineated mountains have characteristic en echelon grooves that likely formed as a result of gravitational sliding. Undivided mountains are partially grooved but exhibit evidence of slumping and are generally lower elevation than the lineated units. Debris lobes and aprons are representative of mottled mountain materials. We have explored the possibility that north and south Hi’iaka Mons were originally one structure. We propose that strike-slip faulting and subsequent rifting separated the mountain units and created a depression which, by further extension during the rifting event, became Hi’iaka Patera. This type of rifting and depression formation is similar to the mechanism of formation of terrestrial pull-apart basins. With comparison to other regional maps of Io and global studies of paterae and mountains, this work provides insight into the general geologic evolution of Io.     

Volatile retention from cometary impacts on the Moon

Impacts of comets and asteroids play an important role in volatile delivery on the Moon. We use a novel method for tracking vapor masses that reach escape velocity in hydrocode simulations of cometary impacts to explore the effects of volatile retention. We model impacts on the Moon to find the mass of vapor plume gravitationally trapped on the Moon as a function of impact velocity. We apply this result to the impactor velocity distribution and find that the total impactor mass retained on the Moon is approximately 6.5% of the impactor mass flux. Making reasonable assumptions about water content of comets and the comet size-frequency distribution, we derive a water flux for the Moon. After accounting for migration and stability of water ice at the poles, we estimate a total 1.3×10⁸-4.3×10⁹ metric tons of water is delivered to the Moon and remains stable at the poles over 1 Ga. A factor of 30 uncertainty in the estimated cometary impact flux is primarily responsible for this large range of values. The calculated mass of water is sufficient to account for the neutron fluxes poleward of 75° observed by Lunar Prospector. A similar analysis for water delivery to the Moon via asteroid impacts shows that asteroids provide six times more water mass via impacts than comets.     

High outflow channels on Mars indicate Hesperian recharge at low latitudes and the presence of Canyon Lakes

A number of martian outflow channels were carved by discharges from large dilational fault zones. These channels were sourced by groundwater, not surface water, and when observed on high-standing plateaus they provide indicators of elevated paleo-groundwater levels. We identify three outflow channels of Hesperian age that issued from a 750-km-long fault zone extending from Candor Chasma to Ganges Chasma. Two of these channels, Allegheny Vallis and Walla Walla Vallis, have sources >2500 m above the topographic datum, too high to be explained by discharge from a global aquifer that was recharged solely in the south polar region. The indicated groundwater levels likely required regional sources of recharge at low latitudes. The floodwaters that erupted from Ophir Cavus to form Allegheny Vallis encountered two ridges that restricted the flow, forming temporary lakes. The flow probably breached or overtopped these obstructions quickly, catastrophically draining the lakes and carving several scablands. After the last obstacle had been breached, a single main channel formed that captured all subsequent flow. We performed hydrologic analyses of this intermediate phase of the flooding, prior to incision of the channel to its present depth. Using floodwater depths of 30-60 m, we calculated flow velocities of 6-15 m s⁻¹ and discharges in the range of . Locally higher flow velocities and discharges likely occurred when the transient lakes were drained. Variable erosion at the channel and scabland crossing of MOLA pass 10644 suggests that the upper 25-30 m may consist of poorly consolidated surface materials underlain by more cohesive bedrock. We infer that an ice-covered lake with a surface elevation >2500 m probably existed in eastern Candor Chasma because this canyon is intersected by the Ophir Catenae fault system from which Allegheny Vallis and Walla Walla Vallis originated. We introduce a new hydrology concept for Mars in which the groundwater system was augmented by recharge from canyon lakes that were formed when water was released by catastrophic melting of former ice sheets in Tharsis by effusions of flood basalts. This model could help to reconcile the expected presence of a thick cryosphere during the Hesperian with the abundant evidence for groundwater as a source for some of the circum-Chryse outflow channels.     

Perched lobe formation in the Gulf of Cadiz: Interactions between gravity processes and contour currents (Algarve Margin, Southern Portugal)

The Gulf of Cadiz is swept by the strong saline Mediterranean Outflow Water (MOW). On the Algarve Margin (South Portugal), this current has constructed fine-grained contourite drifts. This margin is dissected by the Portimao Canyon and three short channels that only incise the upper slope, and are absent on a terrace located at mid-slope depths along the Algarve Margin. High-resolution seismic profiles and sediment cores highlight the original architecture of the sedimentary deposits on this terrace. Coarse-grained lenticular chaotic bodies formed during major relative sea-level lowstands are intercalated within the drift. The lobate shape and sandy nature of the lenticular chaotic bodies and their location at the mouths of the three short channels suggest they are gravity-generated deposits that are perched on the middle continental slope.In the Gulf of Cadiz, the interaction between contour current and gravity processes is strongly controlled by climatic variations and relative sea-level changes during the late Quaternary. During cold periods when sea-level was low, erosion intensified on the continental shelf and the deepest part of MOW was active. Sediment was transported downslope through the channels and deposited on sedimentary lobes perched on the mid-slope terrace. During warm periods when relative sea-level was high, the supply of sediment from the shelf was shut off and the shallowest part of MOW was more active. Contourite drifts fill the channels and bury the sandy lobes.