Evidence for ice flow prior to trough formation in the martian north polar layered deposits

The relative importance of surface mass fluxes and ice flow in shaping the north polar layered deposits (NPLD), now or in the past, remains a fundamental and open question. Motivated by observation of an apparent ice divide on Gemina Lingula (also known as Titania Lobe), we propose a two-stage evolution leading to the present-day topography on that lobe of the NPLD. Ice flow approximately balances surface mass fluxes in the first stage, but in the second stage ice flow has minimal influence and topography is modified predominantly by the formation of troughs. We focus here on evidence for the first stage, by testing the fit of topography between troughs to an ice-flow model. We find that independent model fits on distinct flow paths closely match inter-trough topography, uniformly over a broad region on Gemina Lingula, with mutually consistent and physically reasonable fitting parameters. However, our model requires ice to occupy and flow in spaces where troughs currently incise the ice. We therefore infer that the troughs (and the distribution of mass balance that caused them) post-date deposition of the inter-trough material and its modification by flow. Because trough formation has apparently altered inter-trough topography very little, we infer that trough formation must have been rapid in comparison to the (still unknown) time-scale of flow since troughs began to form. We view the evidence for past flow as strong, but we do not think that topographic evidence alone can be conclusive. Observations of englacial stratigraphy using orbital sounding radars will yield conclusive tests of our inferred mechanism for the formation of inter-trough topography.     

A Terrestrial Analogy for Martian “Accublation Zones” Revealed by Airborne Ice-Penetrating Radar from the East Antarctic Ice Sheet

The surface of the northern martian ice cap is known to contain alternating white-dark scarp-trough features, and one method of their formation involves spatially coherent patterns of ice accumulation and ablation via sublimation. In this so-called “accublation” model, the internal ice sheet layers submerge into the ice sheet across the accumulation zone and emerge to the ice surface across the ablation region. Here we report on an Antarctic analogy to the martian accublation model, which demonstrates that surface mass balance conditions actually affect internal ice-sheet structure as predicted by D. A. Fisher (2000, Icarus144, 289-294). This analogy shows how it is possible for ancient ice to become exposed at the surface of an ice sheet, allowing the paleoclimate record stored within the ice to be measured without the need for deep ice coring.     

Martian polar and circum-polar sulfate-bearing deposits: Sublimation tills derived from the North Polar Cap

Previous spectroscopic studies have shown the presence of hydrated minerals in various kinds of sedimentary accumulations covering and encircling the martian North Polar Cap. More specifically, gypsum, a hydrated calcium sulfate, has been detected on Olympia Planum, a restricted part of the Circum-Polar Dune Field. To further constrain the geographical distribution and the process of formation and accumulation of these hydrated minerals, we performed an integrated morphological, structural and compositional analysis of a key area where hydrated minerals were detected and where the main polar landforms are present. By the development of a spectral processing method based on spectral derivation and by the acquisition of laboratory spectra of gypsum-ice mixtures we find that gypsum-bearing sediment is not restricted to the Olympia Planum dunes but is also present in all kinds of superficial sediment covering the surface of the North Polar Cap and the Circum-Polar Dune Field. Spectral signatures consistent with perchlorates are also detected on these deposits. The interpretation of landforms reveals that this gypsum-bearing sediment was released from the ice cap by sublimation. We thus infer that gypsum crystals that are now present in the Circum-Polar Dune Field derive from the interior of the North Polar Cap. Gypsum crystals that were initially trapped in the ice cap have been released by sublimation of the ice and have accumulated in the form of ablation tills at the surface of the ice cap. These gypsum-bearing sublimation tills are reworked by winds and are transported towards the Circum-Polar Dune Field. Comparison with sulfates found in terrestrial glaciers suggests that gypsum crystals in the martian North Polar Cap have formed by weathering of dust particles, either in the atmosphere prior to their deposition during the formation of the ice cap, and/or in the ice cap after their deposition.     

Morphologic and topographic analyses of debris aprons in the eastern Hellas region, Mars

This paper presents new, detailed analyses of small-scale morphologic and topographic characteristics of martian debris aprons that support Viking-based interpretations of debris aprons as ice-rich flow features derived from local uplands. Fifty-four debris apron complexes in the eastern Hellas region of Mars were examined using Mars Global Surveyor data sets, including Mars Orbiter Camera images and Mars Orbiter Laser Altimeter topographic profiles. Consistent patterns in a suite of small-scale surface textures and geomorphic features observed throughout the population reflect a history of viscous flow and surface degradation through wind ablation and loss of contained ice. A wide variety of shapes seen in topographic profile reveal variations in distribution of contained ice and different stages of apron development and degradation. The results of this study provide new evidence consistent with multiple models of apron formation, including rock glacier, debris-covered glacier, and ice-rich landslide models. Typical eastern Hellas debris aprons formed from a series of large-scale events, emplacing debris that was enriched initially or later by ground ice, complemented by small-scale mass wasting of multiple styles and postemplacement flow of apron masses.     

The atmospheric circulation and dust activity in different orbital epochs on Mars

A general circulation model is used to evaluate changes to the circulation and dust transport in the martian atmosphere for a range of past orbital conditions. A dust transport scheme, including parameterized dust lifting, is incorporated within the model to enable passive or radiatively active dust transport. The focus is on changes which relate to surface features, as these may potentially be verified by observations. Obliquity variations have the largest impact, as they affect the latitudinal distribution of solar heating. At low obliquities permanent CO₂ ice caps form at both poles, lowering mean surface pressures. At higher obliquities, solar insolation peaks at higher summer latitudes near solstice, producing a stronger, broader meridional circulation and a larger seasonal CO₂ ice cap in winter. Near-surface winds associated with the main meridional circulation intensify and extend polewards, with changes in cap edge position also affecting the flow. Hence the model predicts significant changes in surface wind directions as well as magnitudes. Dust lifting by wind stress increases with obliquity as the meridional circulation and associated near-surface winds strengthen. If active dust transport is used, then lifting rates increase further in response to the larger atmospheric dust opacities (hence circulation) produced. Dust lifting by dust devils increases more gradually with obliquity, having a weaker link to the meridional circulation. The primary effect of varying eccentricity is to change the impact of varying the areocentric longitude of perihelion, l, which determines when the solar forcing is strongest. The atmospheric circulation is stronger when l aligns with solstice rather than equinox, and there is also a bias from the martian topography, resulting in the strongest circulations when perihelion is at northern winter solstice. Net dust accumulation depends on both lifting and deposition. Dust which has been well mixed within the atmosphere is deposited preferentially over high topography. For wind stress lifting, the combination produces peak net removal within western boundary currents and southern midlatitude bands, and net accumulation concentrated in Arabia and Tharsis. In active dust transport experiments, dust is also scoured from northern midlatitudes during winter, further confining peak accumulation to equatorial regions. As obliquity increases, polar accumulation rates increase for wind stress lifting and are largest for high eccentricities when perihelion occurs during northern winter. For dust devil lifting, polar accumulation rates increase (though less rapidly) with obliquity above o=25°, but increase with decreasing obliquity below this, thus polar dust accumulation at low obliquities may be increasingly due to dust lifted by dust devils. For all cases discussed, the pole receiving most dust shifts from north to south as obliquity is increased.     

Thermal contraction crack polygons on Mars: A synthesis from HiRISE, Phoenix, and terrestrial analog studies

Thermal contraction crack polygons are complex landforms that have begun to be deciphered on Earth and Mars by the combined investigative efforts of geomorphology, environmental monitoring, physical models, paleoclimate reconstruction, and geochemistry. Thermal contraction crack polygons are excellent indicators of the current or past presence of ground ice, ranging in ice content from weakly cemented soils to debris-covered massive ice. Relative to larger topographic features, polygons may form rapidly, and reflect climate conditions at the time of formation—preserving climate information as relict landforms in the geological record. Polygon morphology and internal textural characteristics can be used to distinguish surfaces modified by the seasonal presence of a wet active layer or dry active layer, and to delimit subsurface ice conditions. Analysis of martian polygon morphology and distribution indicates that geologically-recent thermal contraction crack polygons on Mars form predominantly in an ice-rich latitude-dependent mantle, more likely composed of massive ice deposited by precipitation than by cyclical vapor diffusion into regolith. Regional and local heterogeneities in polygon morphology can be used to distinguish variations in ice content, deposition and modification history, and to assess microclimate variation on timescales of ka to Ma. Analyses of martian polygon morphology, guided by investigations of terrestrial analog thermal contraction crack polygons, strongly suggest the importance of excess ice in the formation and development of many martian thermal contraction crack polygons—implying the presence of an ice-rich substrate that was fractured during and subsequent to obliquity-driven depositional periods and continually modified by ongoing vapor equilibration processes.     

Influence of ice rheology and dust content on the dynamics of the north-polar cap of Mars

The evolution and dynamics of the north-polar cap (residual-ice-cap/layered-deposits complex) of Mars is simulated with a thermomechanical ice-sheet model. We consider a scenario with ice-free initial conditions at 5 Ma before present due to the large obliquities which prevailed prior to this time. The north-polar cap is then built up to its present shape, driven by a parameterized climate forcing (surface temperature, surface mass balance) based on the obliquity and eccentricity history. The effects of different ice rheologies and different dust contents are investigated. It is found that the build-up scenarios require an accumulation rate of approximately 0.15-0.2 mm a⁻¹ at present. The topography evolution is essentially independent of the ice dynamics due to the slow ice flow. Owing to the uncertainties associated with the ice rheology and the dust content, flow velocities can only be predicted within a range of two orders of magnitude. Likely present values are of the order of 0.1-1 mm a⁻¹, and a strong variation over the climatic cycles is found. For all cases, computed basal temperatures are far below pressure melting.     

Atmospheric variations and meteorite production on Mars

Through a combination of aerobraking (drag deceleration) and ablation, meteoroids which enter planetary atmospheres may be slowed sufficiently to soft-land as meteorites. Results of an earlier study suggest that the current 6 mbar atmosphere of Mars is sufficient to aerobrake significant numbers of small (<10 kg) asteroidal-type meteoroids into survivable, low-velocity (<500 m s⁻¹) impacts with the planet's surface. Since rates of meteorite production depend upon the density of Mars's atmosphere, they must also change as the martian climate changes. However, to date, martian meteorite production has received relatively little attention in the literature Here we expand upon our previous work to study martian meteorite production rates and how they depend upon variations of the martian atmosphere, and to estimate the ranges of mass, velocity and entry-angle that produce meteorites. We find that even the current atmosphere of Mars is sufficient to soft-land significant fractions of incident stony and iron objects, and that these fractions increase dramatically for denser martian atmospheres. Therefore, like impact cratering, meteorite populations may preserve evidence of past martian climates.     

Evidence for subsurface water ice in Korolev crater, Mars

Following the work of Kieffer and Titus (2001, Icarus 154, 162-180), we present results of thermal IR observations of Korolev crater, located at ∼73° latitude in the martian northern polar region. Similar to techniques employed by Titus et al. (2003, Science 299, 1048-1050), we use infrared images from the Thermal Emission Imaging System (THEMIS) aboard Mars Odyssey to identify several regions within the crater basin with distinct thermal properties that correlate with topography. The THEMIS results show these regions exhibit temperature variations, spatially within the crater and throughout the martian year. In addition to the variations identified in the THEMIS observations, Mars Global Surveyor Thermal Emission Spectrometer (TES) observations show differences in albedo and temperature of these regions on both daily and seasonal cycles. Modeling annual temperature variations of the surface, we use TES observations to examine the thermal properties of these regions. This analysis reveals the crater interior deposits are likely thick layers (several meters) of high thermal inertia material (water ice, or extremely ice-rich regolith). Spatial variations of the physical properties of these regions are likely due to topography and possibly variations in the subsurface material itself. The nature of these deposits may help constrain polar processes, as well as provide context for the polar lander mission, Phoenix.     

Radiative habitable zones in martian polar environments

The biologically damaging solar ultraviolet (UV) radiation (quantified by the DNA-weighted dose) reaches the martian surface in extremely high levels. Searching for potentially habitable UV-protected environments on Mars, we considered the polar ice caps that consist of a seasonally varying CO2 ice cover and a permanent H2O ice layer. It was found that, though the CO2 ice is insufficient by itself to screen the UV radiation, at approximately 1 m depth within the perennial H2O ice the DNA-weighted dose is reduced to terrestrial levels. This depth depends strongly on the optical properties of the H2O ice layers (for instance snow-like layers). The Earth-like DNA-weighted dose and Photosynthetically Active Radiation (PAR) requirements were used to define the upper and lower limits of the northern and southern polar Radiative Habitable Zone (RHZ) for which a temporal and spatial mapping was performed. Based on these studies we conclude that photosynthetic life might be possible within the ice layers of the polar regions. The thickness varies along each martian polar spring and summer between approximately 1.5 and 2.4 m for H2O ice-like layers, and a few centimeters for snow-like covers. These martian Earth-like radiative habitable environments may be primary targets for future martian astrobiological missions. Special attention should be paid to planetary protection, since the polar RHZ may also be subject to terrestrial contamination by probes.     

Possible remnants of a frozen mud lake in southern Elysium, Mars

In this work we estimate the minimum persistence time of subsurface ice in water rich sediment layers remaining after sublimation of a martian lake. We simulate sublimation of ice from layers of different granulations and thicknesses. Presented results assume insolation and atmospheric conditions characteristic for the present day southern Elysium, where data from Mars Express have identified surface features possibly indicating the very recent presence of a frozen body of water [Murray et al., 2005. Nature 434, 352-356]. The age of these features is estimated to be several million years. On this time scale, we find that most of the water ice must have sublimated away, however remnant ice at a few percent level cannot be excluded. This amount of water ice is sufficient for chemical cementation of the observed features and explains their relatively pristine appearance, without significant signs of erosion.     

Ice grain size and the rheology of the martian polar deposits

The history and dynamics of the martian polar deposits (MPD), the largest known water reservoirs on Mars, are of great interest, but estimates of ice grain size are required before detailed modeling can be performed. We clarify the microphysical processes that may control grain size in the MPD. If the MPD are ∼2% dust by mass, the maximum ice grain size is ∼1 mm due to grain boundary pinning by silicate microparticles. Relatively dusty layers in the MPD will have smaller grain sizes. If MPD ice has a very low impurity content and has experienced a significant amount of strain, grains may reach a steady state size of ∼1.5 to 3 mm due to dynamic recrystallization, wherein a steady state grain size is maintained due to the balance of grain growth and destruction during flow. If the near-bed ice in the MPD is warmed close to its melting point and has been extensively sheared, grain sizes at its base may be between 10 and 40 mm, by analogy with warm, dirty, near-bed ice in terrestrial ice sheets.     

Effects of soil heterogeneity on martian ground-ice stability and orbital estimates of ice table depth

Data from the Mars Odyssey Gamma-Ray Spectrometer (GRS) instrument suite and results from numerical simulations of subsurface ground-ice stability have been used to estimate the depth of martian ground-ice. Geographic correlation between these estimates is remarkable; the relative ice table depth distributions also agree well. However, GRS-based estimates of ice table depth are generally deeper than predictions based on ground-ice stability simulations. This discrepancy may be related to heterogeneities in the martian surface such as rocks, dust, and albedo variations. We develop a multi-dimensional numerical model of ground-ice stability in a heterogeneous martian subsurface and use it to place the first quantitative constraints on the response of the ice table to meter-scale heterogeneities. We find that heterogeneities produce significant undulations/topography in the ice table at horizontal length scales of a few meters. Decimeter scale rocks create localized areas of deep ice, producing a vertical depression of 10-60 cm in the ice table over a horizontal range of 1-2 rock radii. Decimeter scale dust lenses produce locally shallow ice; however the magnitude of the vertical deflection of the ice table is small (1-4 cm). The effects of decimeter scale albedo variations on the ice table are nearly negligible, although albedo very weakly enhances the effects of dark rocks and bright dust on the ice table. Additionally, we investigate the role played by rocks in estimates of ice table depth based on orbital data. Surface rocks can account for more than half of the discrepancy between ice table depths inferred from GRS data and those predicted by theoretical ice-stability simulations that utilize thermophysical observations. Our results have considerable relevance to the up-coming Mars Scout Mission, Phoenix, because they indicate that the uncertainty in the ice table depth of a given region is greater than differences between current depth estimates. Likewise, small-scale depth variability due to heterogeneities at the eventual landing site is potentially greater than differences between current depth estimates.     

Do ice caves exist on Mars?

We have developed a numerical model for assessing the lifetime of ice deposits in martian caves that are open to the atmosphere. Our model results and sensitivity tests indicate that cave ice would be stable over significant portions of the surface of Mars. Ice caves on Earth commonly occur in lava tubes, and Mars has been significantly resurfaced by volcanic activity during its history, including the two main volcanic provinces, the Tharsis and Elysium rises. These areas, known or suspected of having subsurface caves and related voids are among the most favorable regions for the occurrence of ice stability. The martian ice cave model predicts regions which, if caves occur, would potentially be areas of astrobiological importance as well as possible water sources for future human missions to Mars.     

Modeling sublimation of ice exposed by new impacts in the martian mid-latitudes

New impacts in the martian mid-latitudes have exposed near-surface ice. This ice is observed to slowly fade over timescales of months. In the present martian climate, exposed surface ice is unstable during summer months in the mid-latitudes and will sublimate. We model the sublimation of ice at five new impact sites and examine the implications of its persistence. Even with generally conservative assumptions, for most reasonable choices of parameters it is likely that over a millimeter of sublimation occurred in the period during which the ice was observed to fade. The persistence of visible ice through such sublimation suggests that the ice is relatively pure rather than pore-filling. Such ice could be analogous to the nearly pure ice observed by the Phoenix Lander in the “Dodo-Goldilocks” trench and suggests that the high ice contents reported by the Mars Odyssey Gamma Ray Spectrometer at high latitudes extend to the mid-latitudes. Our observations are consistent with a model of the martian ice table in which a layer with high volumetric ice content overlies pore-filling ice, although other structures are possible.     

Albedo control of seasonal South Polar cap recession on Mars

Over the last few decades, General Circulation Models (GCM) have been used to simulate the current martian climate. The calibration of these GCMs with the current seasonal cycle is a crucial step in understanding the climate history of Mars. One of the main climatic signals currently used to validate GCMs is the annual atmospheric pressure cycle. It is difficult to use changes in seasonal deposits on the surface of Mars to calibrate the GCMs given the spectral ambiguities between CO₂ and H₂O ice in the visible range. With the OMEGA imaging spectrometer covering the near infra-red range, it is now possible to monitor both types of ice at a spatial resolution of about 1 km. At global scale, we determine the change with time of the Seasonal South Polar Cap (SSPC) crocus line, defining the edge of CO₂ deposits. This crocus line is not symmetric around the geographic South Pole. At local scale, we introduce the snowdrop distance, describing the local structure of the SSPC edge. Crocus line and snowdrop distance changes can now be used to calibrate GCMs. The albedo of the seasonal deposits is usually assumed to be a uniform and constant parameter of the GCMs. In this study, albedo is found to be the main parameter controlling the SSPC recession at both global and local scale. Using a defrost mass balance model (referred to as D-frost) that incorporates the effect of shadowing induced by topography, we show that the global SSPC asymmetry in the crocus line is controlled by albedo variations. At local scale, we show that the snowdrop distance is correlated with the albedo variability. Further GCM improvements should take into account these two results. We propose several possibilities for the origin of the asymmetric albedo control. The next step will be to identify and model the physical processes that create the albedo differences.     

A historical search for habitable ice at the Phoenix landing site

A time-resolved energy balance model in the latitude range targeted by Phoenix, and extending back in time over the past 10 Ma, has been developed and used to predict the time-varying temperature field in ground ice over scales ranging from minutes to millions of years. The temperature history is compared to the population doubling times of terrestrial psychrophiles as a function of temperature, and the lifetime of analog microbe spores against de-activation by galactic cosmic rays (GCR), in order to assess the habitability of ground ice and surrounding materials that may be sampled by Phoenix. Metrics are derived to quantify “habitability” and compare different model configurations, including total and maximum continuous time, per year, that ground ice temperatures exceed various thresholds, maximum and average dormancy periods, and maximum and average consecutive growing seasons. The key unknowns in assessing the position, and hence the temperature, of the ground ice table at high northern latitude is the fate of the perennial north polar cap at high obliquity. If enough H₂O ice can persist at polar latitudes to buffer at least the high-latitude atmosphere at all orbital configurations, ground ice is found to be relatively shallow over much of the past 10 Ma, and regularly achieves temperatures in excess of those required for the growth of terrestrial psychrophiles. The dry overburden expected at the landing site can easily be sampled by Phoenix, and includes the “sweet spot” that is characterized by the optimal habitability metrics over the past 10 Ma. If the atmosphere is buffered only by low-latitude ice deposits at obliquities greater than about 30°, the frequency and duration of habitable ice is considerably diminished, and the intervening dormancy periods, during which cosmic ray damage accumulates, are correspondingly longer. In all cases, the maximum dormancy period that must be survived by putative martian psychrophiles is at least an order of magnitude greater than the amount of time required to reduce terrestrial psychrophile spore viability by 10⁻⁶ (∼7×¹⁰⁴ years). Depending on the fate of high-obliquity polar ice, the maximum dormancy period can exceed 4×¹⁰⁶ years, a factor of 60 longer than terrestrial psychrophile spore lifetimes. Habitability of martian ground ice is therefore dependent on putative martian psychrophiles developing robustness against GCR deactivation at least an order of magnitude greater than their terrestrial counterparts. Simulations of ground ice throughout the 65° N-72° N latitude range accessible to Phoenix suggest that higher-latitude ground ice has better habitability metrics, although the discrepancy is less than an order of magnitude for all metrics and across the entire latitude range.     

Preservation of Late Amazonian Mars ice and water-related deposits in a unique crater environment in Noachis Terra: Age relationships between lobate debris tongues and gullies

The Amazonian period of Mars has been described as static, cold, and dry. Recent analysis of high-resolution imagery of equatorial and mid-latitude regions has revealed an array of young landforms produced in association with ice and liquid water; because near-surface ice in these regions is currently unstable, these ice-and-water-related landforms suggest one or more episodes of martian climate change during the Amazonian. Here we report on the origin and evolution of valley systems within a degraded crater in Noachis Terra, Asimov Crater. The valleys have produced a unique environment in which to study the geomorphic signals of Amazonian climate change. New high-resolution images reveal Hesperian-aged layered basalt with distinctive columnar jointing capping interior crater fill and providing debris, via mass wasting, for the surrounding annular valleys. The occurrence of steep slopes (>20°), relatively narrow (sheltered) valleys, and a source of debris have provided favorable conditions for the preservation of shallow-ice deposits. Detailed mapping reveals morphological evidence for viscous ice flow, in the form of several lobate debris tongues (LDT). Superimposed on LDT are a series of fresh-appearing gullies, with typical alcove, channel, and fan morphologies. The shift from ice-rich viscous-flow formation to gully erosion is best explained as a shift in martian climate, from one compatible with excess snowfall and flow of ice-rich deposits, to one consistent with minor snow and gully formation. Available dating suggests that the climate transition occurred >8 Ma, prior to the formation of other small-scale ice-rich flow features identified elsewhere on Mars that have been interpreted to have formed during the most recent phases of high obliquity. Taken together, these older deposits suggest that multiple climatic shifts have occurred over the last tens of millions of years of martian history.     

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.     

Martian hillside gullies and icelandic analogs

We report observations of Icelandic hillside gully systems that are near duplicates of gullies observed on high-latitude martian hillsides. The best Icelandic analogs involve basaltic talus slopes at the angle of repose, with gully formation by debris flows initiated by ground water saturation, and/or by drainage of water from upslope cliffs. We report not only the existence of Mars analog gullies, but also an erosional sequence of morphologic forms, found both on Mars and in Iceland. The observations support hypotheses calling for creation of martian gullies by aqueous processes. Issues remain whether the water in each case comes only from surficial sources, such as melting of ground ice or snow, or from underground sources such as aquifers that gain surface access in hillsides. Iceland has many examples of the former, but the latter mechanism is not ruled out. Our observations are consistent with the martian debris flow mechanism of F. Costard et al. (2001c, Science295, 110-113), except that classic debris flows begin at midslope more frequently than on Mars. From morphologic observations, we suggest that some martian hillside gully systems not only involve significant evolution by extended erosive activity, but gully formation may occur in episodes, and the time interval since the last episode is considerably less than the time interval needed to erase the gully through normal martian obliteration processes.