The formation and evolution of youthful gullies on Mars: Gullies as the late-stage phase of Mars’ most recent ice age

Gullies are extremely young erosional/depositional systems on Mars that have been carved by an agent that was likely to have been comprised in part by liquid water [Malin, M.C., Edgett, K.S., 2000. Evidence for recent groundwater seepage and surface runoff on Mars. Science 288, 2330-2335; McEwen, A.S. et al., 2007. A closer look at water-related geologic activity on Mars. Science 317, 1706-1709]. The strong latitude and orientation dependencies that have been documented for gullies require (1) a volatile near the surface, and (2) that insolation is an important factor for forming gullies. These constraints have led to two categories of interpretations for the source of the volatiles: (1) liquid water at depth beneath the melting isotherm that erupts suddenly (“groundwater”), and (2) ice at the surface or within the uppermost layer of soil that melts during optimal insolation conditions (“surface/near-surface melting”). In this contribution we synthesize global, hemispheric, regional and local studies of gullies across Mars and outline the criteria that must be met by any successful explanation for the formation of gullies. We further document trends in both hemispheres that emphasize the importance of top-down melting of recent ice-rich deposits and the cold-trapping of atmospherically-derived H₂O frost/snow as important components in the formation of gullies. This provides context for the incorporation of high-resolution multi-spectral and hyper-spectral data from the Mars Reconnaissance Orbiter that show that (1) cold-trapping of seasonal H₂O frost occurs at the alcove/channel-level on contemporary Mars; (2) gullies are episodically active systems; (3) gullies preferentially form in the presence of deposits plausibly interpreted as remnants of the Late Amazonian emplacement of ice-rich material; and (4) gully channels frequently emanate from the crest of alcoves instead of the base, showing that alcove generation is not necessarily a product of undermining and collapse at these locations, a prediction of the groundwater model. We interpret these various lines of evidence to mean that the majority of gullies on Mars are explained by the episodic melting of atmospherically emplaced snow/ice under spin-axis/orbital conditions characteristic of the last several Myr.     

Observations of martian gullies and constraints on potential formation mechanisms

The discovery of presumably geologically recent gully features on Mars (Malin and Edgett, 2000, Science 288, 2330-2335) has spawned a wide variety of proposed theories of their origin including hypotheses of the type of erosive material. To test the validity of gully formation mechanisms, data from the Mars Global Surveyor spacecraft has been analyzed to uncover trends in the dimensional and physical properties of the gullies and their surrounding terrain. We located 106 Mars Orbiter Camera (MOC) images that contain clear evidence of gully landforms, distributed in the southern mid and high latitudes, and analyzed these images in combination with Mars Orbiter Laser Altimeter (MOLA) and Thermal Emission Spectrometer (TES) data to provide quantitative measurements of numerous gully characteristics. Parameters we measured include apparent source depth and distribution, vertical and horizontal dimensions, slopes, orientations, and present-day characteristics that affect local ground temperatures. We find that the number of gully systems normalized to the number of MOC images steadily declines as one moves poleward of 30° S, reaches a minimum value between 60°-63° S, and then again rises poleward of 63° S. All gully alcove heads occur within the upper one-third of the slope encompassing the gully and the alcove bases occur within the upper two-thirds of the slope. Also, the gully alcove heads occur typically within the first 200 meters of the overlying ridge with the exception of gullies equatorward of 40° S where some alcove heads reach a maximum depth of 1000 meters. While gullies exhibit complex slope orientation trends, gullies are found on all slope orientations at all the latitudes studied. Assuming thermal conductivities derived from TES measurements as well as modeled surface temperatures, we find that 79% of the gully alcove bases lie at depths where subsurface temperatures are greater than 273 K and 21% of the alcove bases lie within the solid water regime. Most of the gully alcoves lie outside the temperature-pressure phase stability of liquid CO₂. Based on a comparison of measured gully features with predictions from the various models of gully formation, we find that models involving carbon dioxide, melting ground ice in the upper few meters of the soil, dry landslide, and surface snowmelt are the least likely to describe the formation of the martian gullies. Although some discrepancies still exist between prediction and observation, the shallow and deep aquifer models remain as the most plausible theories. Interior processes involving subsurface fluid sources are generally favored over exogenic processes such as wind and snowfall for explaining the origin of the martian gullies.     

Residual south polar cap of Mars: Stratigraphy, history, and implications of recent changes

The residual south polar cap (RSPC) of Mars includes a group of different depositional units of CO₂ ice undergoing a variety of erosional processes. Complete summer coverage of the RSPC by ∼6-m/pixel data of the Context Imager (CTX) on Mars Reconnaissance Orbiter (MRO) has allowed mapping and inventory of the units in the RSPC. Unit maps and estimated thicknesses indicate the total volume of the RSPC is currently <380 km³, and represents less than 3% of the total mass of the current Mars atmosphere. Scarp retreat rates in the CO₂ ice derived from comparison of High Resolution Imaging Science Experiment (HiRISE) data with earlier images are comparable to those obtained for periods up to 3 Mars years earlier. These rates, combined with sizes of depressions suggest that the oldest materials were deposited more than 125 Mars years ago. Most current erosion is by backwasting of scarps 1-12 m in height. This backwasting is initiated by a series of scarp-parallel fractures. In the older, thicker unit these fractures form about every Mars year; in thinner, younger materials they form less frequently. Some areas of the older, thicker unit are lost by downwasting rather than by the scarp retreat. A surprising finding from the HiRISE data is the scarcity of visible layering of RSPC materials, a result quite distinct from previous interpretations of layers in lower resolution images. Layers ∼0.1 m thick are exposed on the upper surfaces of some areas, but their timescale of deposition is not known. Late summer albedo changes mapped by the CTX images indicate local recycling of ice, although the amounts may be morphologically insignificant. Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) data show that the primary material of all the different forms of the RSPC is CO₂ ice with only small admixtures of water ice and dust.     

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.     

Impacts onto H2O ice: Scaling laws for melting, vaporization, excavation, and final crater size

Shock-induced melting and vaporization of H₂O ice during planetary impact events are widespread phenomena. Here, we investigate the mass of shock-produced liquid water remaining within impact craters for the wide range of impact conditions and target properties encountered in the Solar System. Using the CTH shock physics code and the new 5-phase model equation of state for H₂O, we calculate the shock pressure field generated by an impact and fit scaling laws for melting and vaporization as a function of projectile mass, impact velocity, impact angle, initial temperature, and porosity. Melt production nearly scales with impact energy, and natural variations in impact parameters result in only a factor of two change in the predicted mass of melt. A fit to the π-scaling law for the transient cavity and transient-to-final crater diameter scaling are determined from recent simulations of the entire cratering process in ice. Combining melt production with π-scaling and the modified Maxwell Z-model for excavation, less than half of the melt is ejected during formation of the transient crater. For impact energies less than about 2 × 10²⁰ J and impact velocities less than about 5 km s⁻¹, the remaining melt lines the final crater floor. However, for larger impact energies and higher impact velocities, the phenomenon of discontinuous excavation in H₂O ice concentrates the impact melt into a small plug in the center of the crater floor.     

Recent rheologic processes on dark polar dunes of Mars: Driven by interfacial water?

In springtime on HiRISE images of the Southern polar terrain of Mars flow-like or rheologic features were observed. Their dark color is interpreted as partly defrosted surface where the temperature is too high for CO₂ but low enough for H₂O ice to be present there. These branching streaks grow in size and can move by an average velocity of up to about 1 m/day and could terminate in pond-like accumulation features. The phenomenon may be the result of interfacial water driven rheologic processes. Liquid interfacial water can in the presence of water ice exist well below the melting point of bulk water, by melting in course of interfacial attractive pressure by intermolecular forces (van der Waals forces e.g.), curvature of water film surfaces, and e.g. by macroscopic weight, acting upon ice. This melting phenomenon can be described in terms of “premelting of ice”. It is a challenging consequence, that liquid interfacial water unavoidably must in form of nanometric layers be present in water ice containing soil in the subsurface of Mars. It is the aim of this paper to study possible rheologic consequences in relation to observations, which seem to happen at sites of dark polar dunes on Mars at present. The model in this work assumes that interfacial water accumulates at the bottom of a translucent water-ice layer above a dark and insolated ground. This is warmed up towards the melting point of water. The evolving layer of liquid interfacial water between the covering ice sheet and the heated ground is assumed to drive downward directed flow-like features on slopes, and it can, at least partially, infiltrate (seep) into a porous ground. There, in at least temporarily cooler subsurface layers, the infiltrated liquid water refreezes and forms ice. The related stress built-up is shown to be sufficient to cause destructive erosive processes. The above-mentioned processes may cause change in the structure and thickness of the covering ice and/or may cause the movement of dune grains. All these processes may explain the observed springtime growing and downward extension of the slope streaks analyzed here.     

Possible role of brines in the darkening and flow-like features on the Martian polar dunes based on HiRISE images

This work identifies and describes features of the changing seasonal frost-covered surface of Mars based on HiRISE images, and analyses the possibility that ephemeral liquid brine formation produces them. Because some of these dark features show flow-like appearance, and salts on Mars are present, liquid brines might be also present, possibly accounting for the changing droplet-like features on the Phoenix lander.We observed in-situ darkening and movement of dark features (or movement of the darkening front) on seasonal frost-covered polar dunes. Darkening and brightening may happen within several meters from each other during local spring. Darkening always starts from the bottom and moves up, while brightening progresses from top and moves toward the bottom between the small dune ripples. Brightening occurs during the springtime warming on time scales of several days close to the sites of darkening; therefore, dark material falling from the air, and refreezing of bright ice on it, does not adequately explain the observations. Interpreting the observations as brine-related melting or refreezing also poses problems, but because brine may engulf salt grains or ice blocks, phase changes here could be influenced by factors other than temperature values, and could produce the observations.Analysis of absolute albedo changes indicates that the flow-like features are the darkest at their lower frontal end, sometimes darker than the dark spot from which they originate. A bright halo (white collar) also forms around these spots, possibly due to refreezing. Inside the observed larger spots an outer gray area surrounds the central darkest cores, which is about 10 cm lower than the surrounding bright CO₂ ice. At those places, most or all of the CO₂ ice deposited earlier has disappeared, and H₂O ice is present. Observations of dark flow features moving on the top of this H₂O rich layer suggest even if the flow features start as dry dune avalanches of rolling grains, their dark material heated by solar insolation is in contact with H₂O ice and may produce brines.     

Groundwater discharge and gully formation on martian slopes

Young gullies and gully deposits on walls of martian craters have been cited as evidence that liquid water flowed on the surface of Mars relatively recently. Effects of variable environmental conditions at the surface of Mars are modeled and applied to the case of groundwater emergence from shallow aquifers to investigate whether groundwater is a viable source to enable the erosion of these gullies. The model includes detailed treatment of ice growth in the aquifer. Model results indicate that groundwater discharge can be maintained under the current environmental conditions if the aquifer permeability is like that of terrestrial gravel or higher, if the aquifer is 350 K or warmer, or if the aquifer is a brine with a freezing point depressed to 250 K or below. Groundwater discharge cannot be maintained for the conservative case of a cold, pure water, semi-pervious aquifer. Cold (275 K) pure water pervious (gravel) aquifers, warm (350 K) pure water semi-pervious aquifers, and cold (275 K) CaCl₂ brine semi-pervious aquifers all exhibit a dependence of discharge on season, latitude and slope orientation in our modeling. Seasonal, latitudinal and azimuthal discharge variations are strongest for cold CaCl₂ brine semi-pervious aquifers, with discharges from this aquifer type favoring equator-facing slopes at mid and high southern latitudes. At all latitudes and slope azimuths under our nominal conditions, the cold pure water pervious aquifer, the cold pure water semi-pervious aquifer and the cold CaCl₂ brine semi-pervious aquifer all freeze completely shortly after the simulations are started. Discharge restarts in the summer for the cold pure water pervious aquifer and the cold brine aquifer, but discharge does not restart for the cold pure water semi-pervious aquifer. The warm pure water semi-pervious aquifer maintains daily seeps throughout the year at all but high latitudes. In the case of the cold pure water pervious aquifer, approximately 500,000 m³ of water could be discharged from a mid-latitude, 150-m thick aquifer with a 20-m wide seepage face orientated towards the equator or the pole after a single undermining-induced event before ice growth seals the seepage face. For a brine semi-pervious aquifer with the same dimensions, 200-300 m³ of water could be released from a mid-latitude 20-m wide equator-facing seepage face before the fresh exposure is sealed for the fall and winter seasons. Our results do not rule out groundwater emergence as a means of creating some recent gullies, but they indicate that rather special and perhaps unusual conditions would be required.     

Erosional modification and gully formation at Meteor Crater, Arizona: Insights into crater degradation processes on Mars

Hydrogeological modification of Meteor Crater produced a spectacular set of gullies throughout the interior wall in response to rainwater precipitation, snow melting, and possible groundwater discharge. The crater wall has an exceptionally well-developed centripetal drainage pattern consisting of individual alcoves, channels, and fans. Some of the gullies originate from the rim crest and others from the middle crater wall where a lithologic transition occurs; broad gullies occur along the crater corner radial faults. Deeply incised alcoves are well developed on the soft Coconino Sandstone exposed on the middle crater wall, beneath overlying dolomite. In general, the gully locations are along crater wall radial fractures and faults, which are favorable locales of erosion due to preferential rock breakup from faulting, and groundwater flow/discharge; these structural discontinuities are also the locales where the surface runoff from rain precipitation and snow melting can preferentially flow, causing erosion and crater degradation. Channels are well developed on the talus deposits and alluvial fans on the periphery of the crater floor. Caves exposed on the lower crater level point to percolation of surface runoff and selective discharge through fractures on the crater wall. In addition, lake sediments on the crater floor provide significant evidence of a past pluvial climate, when the water table was higher, and groundwater may have seeped from springs on the crater wall. Although these hydrological processes continue at Meteor Crater today, conditions at the crater are much more arid than they were soon after impact, reflecting a climatic shift. This climate shift and the hydrological modifications observed at Meteor Crater provide insights for landscape sculpturing on Mars during various parts of its history.     

Micrometer-sized ice particles for planetary-science experiments - I. Preparation, critical rolling friction force, and specific surface energy

Coagulation models assume a higher sticking threshold for micrometer-sized ice particles than for micrometer-sized silicate particles. However, in contrast to silicates, laboratory investigations of the collision properties of micrometer-sized ice particles (in particular, of the most abundant H₂O-ice) have not been conducted yet. Thus, we used two different experimental methods to produce micrometer-sized H₂O-ice particles, i.e. by spraying H₂O droplets into liquid nitrogen and by spraying H₂O droplets into a cold nitrogen atmosphere. The mean particle radii of the ice particles produced with these experimental methods are (1.49 ± 0.79) μm and (1.45 ± 0.65) μm. Ice aggregates composed of the micrometer-sized ice particles are highly porous (volume filling factor: ? = 0.11 ± 0.01) or rather compact (volume filling factor: ? = 0.72 ± 0.04), depending on the method of production. Furthermore, the critical rolling friction force of F_Roll,ice = (114.8 ± 23.8) × 10⁻¹⁰ N was measured for micrometer-sized ice particles, which exceeds the critical rolling friction force of micrometer-sized SiO₂ particles . This result implies that the adhesive bonding between micrometer-sized ice particles is stronger than the bonding strength between SiO₂ particles. An estimation of the specific surface energy of micrometer-sized ice particles, derived from the measured critical rolling friction forces and the surface energy of micrometer-sized SiO₂ particles, results in γ_ice = 0.190 J m⁻².     

Gully formation on Mars: Two recent phases of formation suggested by links between morphology, slope orientation and insolation history

The unusual 80 km diameter Noachian-aged Asimov crater in Noachis Terra (46°S, 5°E) is characterized by extensive Noachian-Hesperian crater fill and a younger superposed annulus of valleys encircling the margins of the crater floor. These valleys provide an opportunity to study the relationships of gully geomorphology as a function of changing slope orientation relative to solar insolation. We found that the level of development of gullies was highly correlated with slope orientation and solar insolation. The largest and most complex gully systems, with the most well-developed fluvial landforms, are restricted to pole-facing slopes. In contrast, gullies on equator-facing slopes are smaller, more poorly developed and integrated, more highly degraded, and contain more impact craters. We used a 1D version of the Laboratoire de Météorologie Dynamique GCM, and slope geometries (orientation and angle), driven by predicted spin-axis/orbital parameter history, to assess the distribution and history of surface temperatures in these valleys during recent geological history. Surface temperatures on pole-facing slopes preferential for water ice accumulation and subsequent melting are predicted to occur as recently as 0.5-2.1 Ma, which is consistent with age estimates of gully activity elsewhere on Mars. In contrast, the 1D model predicts that water ice cannot accumulate on equator-facing slopes until obliquities exceed 45°, suggesting they are unlikely to have been active over the last 5 Ma. The correlation of the temperature predictions and the geological evidence for age differences suggests that there were two phases of gully formation in the last few million years: an older phase in which top-down melting occurred on equator-facing slopes and a younger more robust phase on pole-facing slopes. The similarities of small-scale fluvial erosion features seen in the gullies on Mars and those observed in gullies cut by seasonal and perennial snowmelt in the Antarctic Dry Valleys supports a top-down melting origin for these gullies on Mars.     

Dependence of the effective rate constants for the hydrogenation of CO on the temperature and composition of the surface

Experiments on the hydrogenation of CO on pure CO and CO-H₂O mixed ice have been performed at temperatures between 8 and 20 K. We obtained temperature and compositional dependence of the effective reaction rate constants. Results indicate that hydrogenation proceeds efficiently on pure solid CO and CO-H₂O mixed ice at temperatures below 10 and 20 K, respectively. Rate constants for pure CO decreased significantly at 12 K compared to those obtained with CO-H₂O mixed ice. Hydrogenation of CO at temperatures greater than 12 K were catalyzed by the H₂O adjacent to the CO. The importance of the experimental results for some relevant astrophysical environments has also been outlined.     

On liquid phases in cometary nuclei

In this paper we review the relevant literature and investigate conditions likely to lead to melting of H₂O ice, methanol (CH₃OH) ice, ethane (C₂H₆) ice and other volatile ices in cometary nuclei. On the basis of a heat balance model which takes account of volatiles loss, we predict the formation of occasional aqueous and hydrocarbon liquid phases in subsurface regions at heliocentric distances, r_h of 1–3 AU, and 5–12 AU, respectively. Low triple-point temperatures and low vapour pressures of C₂H₆, C₃H₈, and some higher-order alkanes and alkenes, favour liquid phase formation in cometary bodies at high r_h. Microporosity and the formation of a stabilization crust occluding the escape of volatiles facilitate liquid-phase formation. Characteristics of the near-surface which favour subsurface melting include; low effective surface emissivity (at low r_h), high amorphous carbon content, average pore sizes of ～10 μm or less, presence of solutes (e.g. CH₃OH), mixtures of C₂–C₆ hydrocarbons (for melting at high r_h), diurnal thermal cycling, and slow rotation rate. Applying the principles of soil mechanics, capillary forces are shown to initiate pre-melting phenomena and subsequent melting, which is expected to impart considerable strength of ～10⁴ Pa in partially saturated layers, reducing porosity and permeability, enhancing thermal conductivity and heat transfer. Diurnal thermal cycling is expected to have a marked effect on the composition and distribution of H₂O ice in the near-surface leading to frost heave-type phenomena even where little if any true melting occurs. Where melting does take place, capillary suction in the wetted zone has the potential to enhance heat transfer via capillary wetting in a low-gravity environment, and to modify surface topography creating relatively smooth flat-bottomed features, which have a tendency to be located within small depressions. An important aspect of the “wetted layer” model is the prediction that diurnal melt–freeze cycles alter the mixing ratio vs. depth of solutes present, or of other miscible components, largely through a process of fractional crystallization, but also potentially involving frost heave. Wetted layers are potentially durable and can involve significant mass transport of volatile materials in the near-surface, increasing in extent over many rotations of the nucleus prior to and just after perihelion passage, and causing stratification and trapping of the lowest-melting mixtures at depths of several metres. A possible mechanism for cometary outbursts is proposed involving a heat pulse reaching the liquid phase in the deepest wetted zone, leading to supersaturation and triggering the sudden release under pressure of dissolved gases, in particular CO₂, CO, CH₄ or N₂, contained beneath a consolidated near-surface layer. This study indicates that liquid water can persist for long periods of time in the near-surface of some intermediate-sized bodies (10²–10³ km radius) within protoplanetary discs.     

Investigating gully flow emplacement mechanisms using apex slopes

The origin of the martian gullies has been much debated since their discovery by the Mars Orbiter Camera (MOC, Malin, M.C., Edgett, K.S. [2000]. Science 288, 2330-2335). Several previous studies have looked at slope gradients in and around gullies, but none have used Digital Elevation Models (DEMs) from the High Resolution Imaging Science Experiment (HiRISE, McEwen, A.S., and 14 colleagues [2007]. J. Geophys. Res. 112 (E05), E0505S02), which has a pixel scale down to 25 cm/pixel. We use five 1 m/post HiRISE DEMs to measure gully apex slopes, the local channel gradient at the upslope extent of the gully debris apron, which marks a shift from erosion to deposition. The apex slope provides information about whether a flow was likely a typical dry granular flow (begins depositing on slopes ∼21°) or fluidized by some extra mechanism (depositing on shallower slopes). We find that 72% of the 75 gully fans studied were likely emplaced by fluidized flows. Relatively old gullies appear more likely to have hosted fluidized flows than relatively fresh gullies. This suggests a time and location dependent fluidizing agent, possibly liquid water produced in a different climate as previously proposed. Our results do not provide evidence for water-rich flows in gullies today.     

Possible fossil H2O liquid-ice interfaces in the Martian crust

Throughout the northern equatorial region of Mars, extensive areas have been uniformly stripped, roughly to a constant depth. These terrains vary widely in their relative ages. A model is described here to explain this phenomenon as reflecting the vertical distribution of H₂O liquid and ice in the crust. Under present conditions the Martian equatorial regions are stratified in terms of the stability of water ice and liquid water. This arises because the temperature of the upper 1 or 2 km is below the melting point of ice and liquid is stable only at greater depth. It is suggested here that during planetary outgassing earlier in Martian history H₂O was injected into the upper few kilometers of the crust by subsurface and surface volcanic eruption and lateral migration of the liquid and vapor. As a result, a discontinuity in the physical state of materials developed in the Martian crust coincident with the depth of H₂O liquid-ice phase boundary. Material above the boundary remained pristine; material below underwent diagenetic alteration and cementation. Subsequently, sections of the ice-laden zone were erosionally stripped by processes including eolian deflation, gravitational slump and collapse, and fluvial transport due to geothermal heating and melting of the ice. The youngest plains which display this uniform stripping may provide a minimum stratigraphic age for the major period of outgassing of the planet. Viking results suggest that the total amount of H₂O outgassed is less than half that required to fill the ice layer, hence any residual liquid eventually found itself in the upper permafrost zone or stored in the polar regions. Erosion stopped at the old liquid-ice interface due to increased resistance of subjacent material and/or because melting of ice was required to mobilize the debris. Water ice may remain in uneroded regions, the overburden of debris preventing its escape to the atmosphere. Numerous morphological examples shown in Viking and Mariner 9 images suggest interaction of impact, volcanic, and gravitational processes with the ice-laden layer. Finally, volcanic eruptions into ice produces a highly oxidized friable amorphous rock, palagonite. Based on spectral reflectance properties, these materials may provide the best analog to Martian surface materials. They are easily eroded, providing vast amounts of eolian debris, and have been suggested (Toulmin et al., 1977) as possible source rocks for the materials observed at the Viking landing sites.     

Seasonal melting of surface water ice condensing in martian gullies

In this work we consider when and how much liquid water during present climate is possible within the gullies observed on the surface of Mars. These features are usually found on poleward directed slopes. We analyse the conditions for melting of H₂O ice, which seasonally condenses within the gullies. We follow full annual cycle of condensation and sublimation of atmospheric CO₂ and H₂O, accounting for the heat and mass transport in the soil. During the summer, once the facets of the gullies are exposed to the Sun the water ice can melt and evaporate. Two mid latitude locations in both hemispheres are considered. The model includes both the rough geometry of the gullies as well as the slope of the surface where the gullies appear. It is an extension of the model developed to calculate condensation of CO₂ ice in troughs of different sizes, including polygonal features on Mars (Kossacki and Markiewicz, 2002, Icarus 160, 73; Kossacki et al., 2003, Planet. Space Sci. 51, 569). We have found, that water ice accumulated during winter can undergo transition to the liquid phase after complete sublimation of CO₂ ice. The amount of liquid water depends on water content in the atmosphere and on the local wind speed. It is probably not enough to destabilise the slope and cause flow of the surface material. However, even the small amounts of liquid water predicted, can play an important role in surface chemistry, in increasing the cohesive strength of the soil's surface layer and possibly may have some exobiological implications.     

Origin,Bulk Chemical Composition And Physical Structure Of The Galilean Satellites Of Jupiter: A Post-Galileo Analysis

The origin of Jupiter and the Galilean satellite system is examinedin the light of the new data that has been obtained by the NASA Galileo Project. In particular, special attention is given to a theory of satellite origin which was put forward at the start of the Galileo Mission and on the basis of which several predictions have now been proven successful (Prentice, 1996a–c). These predictions concern the chemical composition of Jupiter's atmosphere and the physical structure of the satellites. According to the proposed theory of satellite origin, each of the Galilean satellites formed by chemical condensation and gravitational accumulation of solid grains within a concentricfamily of orbiting gas rings. These rings were cast off equatorially by the rotating proto-Jovian cloud (PJC) which contracted gravitationally to form Jupiter some 4 billion years ago. The PJC formed from the gas and grains left over from the gas ring that had been shed at Jupiter's orbit by the contracting proto-solar cloud (PSC). Supersonic turbulentconvection provides the means for shedding discrete gas rings.The temperatures T_n of the system of gas rings shed by the PSCand PJC vary with their respective mean orbital radii R_n (n = 0, 1, 2, Ϊ ) according as T_n ∝ R_n ^-0.9. If the planet Mercury condenses at 1640 K, so accounting for the high density ofthat planet via a process of chemical fractionation between iron and silicates, then T_n at Jupiter's orbit is 158 K. Only 35% of the water vapour condenses out. Thus fractionation between rock and ice, together with an enhancement in the abundance of solids relative to gas which takes place through gravitational sedimentation of solids onto the mean orbit of the gas ring, ensures nearly equal proportions of rock and ice in each of Ganymede and Callisto. Io and Europa condense above the H₂O ice point and consist solely of hydrated rock (h-rock). The Ganymedan condensate consists of h-rock and H₂O ice. For Callisto, NH₃ ice makes up ∼5% of the condensate mass next to h-rock (∼50%) and H₂O ice (∼45%). Detailed thermal and structural models for each of Europa, Ganymedeand Callisto are constructed on the basis of the above initial bulk chemicalcompositions. For Europa (E), a predicted 2-zone model consisting of a dehydrated rock core of mass 0.912M_E and a 150 km thick frozen mantle of salty H₂O yields a moment-of-inertiacoefficient which matches the Galileo Orbiter gravity measurement. For Ganymede (G), a 3-zone model possessing an inner core of solid FeS and mass ∼0.116M_G, and an outer H₂O ice mantle of mass ∼0.502M_G is needed to explain the gravity data.Ganymede's native magnetic field was formed by thermoremanent magnetization of Fe₃O₄. A new Callisto (C) model is proposed consisting of a core of mass 0.826M_C containing a uniform mixture of h-rock (60% by mass) and H₂O and NH₃ ices, and capped by a mantle of pure ice. This model may have the capacity to yield a thin layer of liquid NH₃ċ2H₂O at the core boundary, in line with Galileo's discovery of an induced magnetic field This revised version was published online in July 2006 with corrections to the Cover Date.     

Spatial variability,composition and thickness of the seasonal north polar cap of Mars in mid-spring

The planetary fourier spectrometer (PFS) for the Mars express mission (MEX) is an infrared spectrometer operating in the wavelength range from 1.2 to 45 μm by means of two spectral channels, called SWC (short wavelength channel) and LWC (long wavelength channel), covering, respectively, 1.2–5.5 and 5.5–45 μm.The middle-spring Martian north polar cap (Ls∼40°) has been observed by PFS/MEX in illuminated conditions during orbit 452. The SWC spectra are here used to study the cap composition in terms of CO₂ ice, H₂O ice and dust content. Significant spectral variation is noted in the cap interior, and regions of varying CO₂ ice grain sizes, water frost abundance, CO₂ ice cover and dust contamination can be distinguished. In addition, we correlate the infrared spectra with an image acquired during the same orbit by the OMEGA imaging spectrometer and with the altimetry from MOLA data. Many of the spectra variations correlate with heterogeneities noted in the image, although significant spectral variations are not discernible in the visible. The data have been divided into five regions with different latitude ranges and strong similarities in the spectra, and then averaged. Bi-directional reflectance models have been run with the appropriate lighting geometry and used to fit the observed data, allowing for CO₂ ice and H₂O ice grain sizes, dust and H₂O ice contaminations in the form of intimate granular mixtures and spatial mixtures.A wide annulus of dusty water ice surrounds the recessing CO₂ seasonal cap. The inner cap exhibits a layered structure with a thin CO₂ layer with varying concentrations of dark dust, on top of an H₂O ice underneath ground. In the best-fits, the ices beneath the top layer have been considered as spatial mixtures. The results are still very good everywhere in the spectral range, except where the CO₂ ice absorption coefficients are such that even a thin layer is enough to totally absorb the incoming radiation (i.e. the band is saturated). This only happens around 3800 cm⁻¹, inside the strong 2.7-μm CO₂ ice absorption band. The effect of finite snow depth has been investigated through a layered albedo model. The thickness of the CO₂ ice deposits increases with latitude, ranging from 0.5–1 g cm⁻² within region II to 60–80 g cm⁻² within the highest-latitude (up to 84°N) region V.Region I is at the cap edge and extends from 65°N to 72°N latitude. No CO₂ ice is present in this region, which consists of relatively large grains of water ice (20 μm), highly contaminated by dust (0.15 wt%). The adjacent region II is a narrow region [76–79°N] right at the edge of the north residual polar cap. This region is very distinct in the OMEGA image, where it appears to surround the whole residual cap. The CO₂ ice features are barely visible in these spectra, except for the strong saturated 2.7 μm band. It basically consists of a thin layer of 5-mm CO₂ ice on top of an H₂O ice layer with the same composition as region I. A third interesting region III is found all along the shoulder of the residual cap [79–81°N]. It extends over 1.5 km in altitude and over only 2° of latitude and consists of CO₂ ice with a large dust content. It is an admixture of CO₂ ice (3–4 mm), with several tens of ppm by mass of water ice and more than 2 ppt by mass of dust. The surface temperatures have been retrieved from the LWC spectra for each observation. We found an increase in the surface temperature in this region, indicating a spatial mixture of cold CO₂ ice and warmer dust/H₂O ice. Region IV is close to the top of the residual cap [81–84°N]; it is much brighter than region III, with a dust content 10 times lower than the latter. The CO₂ grain size is 3 mm and strong CO₂ ice features are present in the data, indicating a thicker CO₂ ice layer than in region II (1–2 g cm⁻²). The final region V is right at the top of the residual cap (⩾84°N). It is “pure” CO₂ ice (no dust) of 5 mm grain sizes, with 30 ppm by weight of water ice. The CO₂ ice features are very pronounced and the 2.7 μm band is saturated. The optical thickness is close to the semi-infinite limit (30–40 g cm⁻²). Assuming a snowpack density of 0.5 g cm⁻³, we get a minimum thickness of 1–2 cm for the top-layer of regions II and III, 4–10 cm for region IV, and ⩾60–80 cm thickness for region V. These values are in close agreement with several recent results for the south seasonal polar cap.These results should provide new, useful constraints in models of the Martian climate system and volatile cycles.     

Adsorptive fractionation of HDO on JSC MARS-1 during sublimation with implications for the regolith of Mars

A chamber was constructed to simulate the boundary between the ice table, regolith and atmosphere of Mars and to examine fractionation between H₂O and HDO during sublimation under realistic martian conditions of temperature and pressure. Thirteen experimental runs were conducted with regolith overlying the ice. The thickness and characteristic grain size of the regolith layer as well as the temperature of the underlying ice was varied. From these runs, values for the effective diffusivity, taking into account the effects of adsorption, of the regolith were derived. These effective diffusivities ranged from 1.8 × 10⁻⁴ m² s⁻¹ to 2.2 × 10⁻³ m² s⁻¹ for bare ice and from 2.4 × 10⁻¹¹ m² s⁻¹ to 2.0 × 10⁻⁹ m² s⁻¹ with an adsorptive layer present. From these, latent heats of adsorption of 8.6 ± 2.6 kJ mol⁻¹ and 9.3 ± 2.8 kJ mol⁻¹ were derived at ice-surface temperatures above 223 ± 8 K and 96 ± 28 kJ mol⁻¹ and 104 ± 31 kJ mol⁻¹ respectively for H₂O and HDO were derived at colder temperatures. For temperatures below 223 K, the effective diffusivity of HDO was found to be lower than the diffusivity of H₂O by 40% on average, suggesting that the regolith was adsorptively fractionating the sublimating gas with a fractionation factor of 1.96 ± 0.74. Applying these values to Mars predicts that adsorbed water on the regolith is enriched in HDO compared to the atmosphere, particularly where the regolith is colder. Based on current observations, the D/H ratio of the regolith may be as high as 21 ± 8 times VSMOW at 12°S and L_S = 357° if the regolith is hydrated primarily by the atmosphere, neglecting any hydration from subsurface ice.     

A model for the interior structure, evolution, and differentiation of Callisto

The recently measured dimensionless moment of inertia (MoI) factor for Callisto of 0.3549±0.0042 (Anderson et al., 2001, Icarus, 153, 157-161) poses a problem: its value cannot be explained by a model in which Callisto is completely differentiated into an ice shell above a rock shell and an iron core such as its neighboring satellite Ganymede nor can it be explained by a model of a homogeneous, undifferentiated ice-rock satellite. We show that Callisto may be incompletely differentiated into an outer ice-rock shell in which the volumetric rock concentration is close to the primordial one at the surface and decreases approximately linearly with depth, an ice mantle mostly depleted of rock, and an about 1800 km rock-ice core in which the rock concentration is close to the close-packing limit. The ice-rock shell thickness depends on uncertain rheology parameters and the heat flow and can be roughly 50 to 150 km thick. We show that if Callisto accreted from a mix of metal bearing rock and ice and if the average size of the rocks was of the order of meters to tens of meters, then Callisto may have experienced a gradual, but still incomplete unmixing of the two components. An ocean in Callisto at a depth of 100-200 km is difficult to obtain if the ice is pure H₂O and if the ice-rock lithosphere is 100 km or more thick; a water ocean is more plausible for ice contaminated by ammonia, methane or salts; or for pure H₂O at a depth of 400-600 km.