Subsurface structure of Planum Boreum from Mars Reconnaissance Orbiter Shallow Radar soundings

We map the subsurface structure of Planum Boreum using sounding data from the Shallow Radar (SHARAD) instrument onboard the Mars Reconnaissance Orbiter. Radar coverage throughout the 1,000,000-km² area reveals widespread reflections from basal and internal interfaces of the north polar layered deposits (NPLD). A dome-shaped zone of diffuse reflectivity up to 12 μs (∼1-km thick) underlies two-thirds of the NPLD, predominantly in the main lobe but also extending into the Gemina Lingula lobe across Chasma Boreale. We equate this zone with a basal unit identified in image data as Amazonian sand-rich layered deposits [Byrne, S., Murray, B.C., 2002. J. Geophys. Res. 107, 5044, 12 pp. doi:10.1029/2001JE001615; Fishbaugh, K.E., Head, J.W., 2005. Icarus 174, 444-474; Tanaka, K.L., Rodriguez, J.A.P., Skinner, J.A., Bourke, M.C., Fortezzo, C.M., Herkenhoff, K.E., Kolb, E.J., Okubo, C.H., 2008. Icarus 196, 318-358]. Elsewhere, the NPLD base is remarkably flat-lying and co-planar with the exposed surface of the surrounding Vastitas Borealis materials. Within the NPLD, we delineate and map four units based on the radar-layer packets of Phillips et al. [Phillips, R.J., and 26 colleagues, 2008. Science 320, 1182-1185] that extend throughout the deposits and a fifth unit confined to eastern Gemina Lingula. We estimate the volume of each internal unit and of the entire NPLD stack (821,000 km³), exclusive of the basal unit. Correlation of these units to models of insolation cycles and polar deposition [Laskar, J., Levrard, B., Mustard, J.F., 2002. Nature 419, 375-377; Levrard, B., Forget, F., Montmessin, F., Laskar, J., 2007. J. Geophys. Res. 112, E06012, 18 pp. doi:10.1029/2006JE002772] is consistent with the 4.2-Ma age of the oldest preserved NPLD obtained by Levrard et al. [Levrard, B., Forget, F., Montmessin, F., Laskar, J., 2007. J. Geophys. Res. 112, E06012, 18 pp. doi:10.1029/2006JE002772]. We suggest a dominant layering mechanism of dust-content variation during accumulation rather than one of lag production during periods of sublimation.     

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

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

Reading the climate record of the martian polar layered deposits

The martian polar regions have layered deposits of ice and dust. The stratigraphy of these deposits is exposed within scarps and trough walls and is thought to have formed due to climate variations in the past. Insolation has varied significantly over time and caused dramatic changes in climate, but it has remained unclear whether insolation variations could be linked to the stratigraphic record. We present a model of layer formation based on physical processes that expresses polar deposition rates of ice and dust in terms of insolation. In this model, layer formation is controlled by the insolation record, and dust-rich layers form by two mechanisms: (1) increased summer sublimation during high obliquity, and (2) variations in the polar deposition of dust modulated by obliquity variations. The model is simple, yet physically plausible, and allows for investigations of the climate control of the polar layered deposits (PLD). We compare the model to a stratigraphic column obtained from the north polar layered deposits (NPLD) (Fishbaugh, K.E., Hvidberg, C.S., Byrne, S., Russel, P.S., Herkenhoff, K.E., Winstrup, M., Kirk, R. [2010a]. Geophys. Res. Lett., 37, L07201) and show that the model can be tuned to reproduce complex layer sequences. The comparison with observations cannot uniquely constrain the PLD chronology, and it is limited by our interpretation of the observed stratigraphic column as a proxy for NPLD composition. We identified, however, a set of parameters that provides a chronology of the NPLD tied to the insolation record and consistently explains layer formation in accordance with observations of NPLD stratigraphy. This model dates the top 500 m of the NPLD back to ～1 million years with an average net deposition rate of ice and dust of 0.55 mm a^?1. The model stratigraphy contains a quasi-periodic ～30 m cycle, similar to a previously suggested cycle in brightness profiles from the NPLD (Laskar, J., Levrard, B., Mustard, F. [2002]. Nature, 419, 375–377; Milkovich, S., Head, J.W. [2005]. J. Geophys. Res. 110), but here related to half of the obliquity cycles of 120 and 99 kyr and resulting from a combination of the two layer formation mechanisms. Further investigations of the non-linear insolation control of PLD formation should consider data from other geographical locations and include radar data and other stratigraphic datasets that can constrain the composition and stratigraphy of the NPLD layers.     

Formation of an Hesperian-aged sedimentary basin containing phyllosilicates in Coprates Catena,Mars

The extensive light-toned deposits in canyons and troughs in Valles Marineris provide evidence of formation through water-related processes. As such, these deposits offer a window to past conditions on Mars. We study a small outcrop of light-toned deposits in a closed trough in Coprates Catena, a chain of collapse pits to the south-east of the main Valles Marineris system. A well-exposed sequence of deposits on the base of the north wall of the trough offers a 220 m section for geochemical and morphologic analysis. Using CRISM data we identify the presence of both phyllosilicates and sulfates and/or opaline silica in the light toned deposits, which vary in relative strength with elevation. We observe a trend in the dominant mineralogical signal, with Al phyllosilicates occurring near the base of the deposits, both below and above a band of Fe/Mg phyllosilicates, before a transition to more sulfate- or opaline silica-rich material near the top of the section. This trend likely reflects a change in the chemistry of the water in which the deposits formed. Using a HiRISE Digital Elevation Model, we find that the layers in the light-toned deposits on both sides of the trough dip gently towards the center of the trough, with a dip direction that aligns with the strike of the trough, suggesting that the light-toned deposits formed after the trough. Our general morphologic and mineralogical observations fit well with significant amounts of water in the trough. The deposits are too small to be dated using crater counting techniques, however, our crater analysis suggests that the plains in which the trough formed are probably Late Hesperian in age. If the chemistry of the light-toned deposits reflects the primary depositional mineralogy, then this and other small troughs in Coprates Catena might provide evidence of limited phyllosilicate formation in this region towards the end of the Hesperian era on Mars.     

The effects of martian orbital variations upon the sublimation and relaxation of north polar troughs and scarps

The North Polar Layered Deposits (PLD) of Mars are climatologically significant because they represent the largest actively-exchanging reservoir of martian water. The kilometer-scale topography of the North PLD is dominated by troughs and scarps, which exhibit highly-correlated surface slopes and total depths. The most widespread theories of PLD evolution presume that the asymmetrical nature of North PLD troughs (characterized by equatorward-facing slopes that are generally steeper than poleward-facing slopes) is the result of preferential H₂O sublimation from equatorward-facing trough walls. However, our orbitally-modulated simulations indicate that there is no long-term sublimation advantage of equatorward-facing trough walls, because of the effects of obliquity upon the slope dependence of sublimation rate. In contrast, we propose that viscous relaxation of subsurface water is consistent with the slope and depth distributions of North PLD troughs and scarps. The results of our finite element simulations suggest that a mere 2 K difference in the subsurface temperatures of opposing trough walls is sufficient to produce the observed slope disparity, due to the slower rate of uplift of colder poleward-facing trough walls. Both our sublimation and relaxation simulations indicate that present-day North PLD troughs have formed since 5 Ma and are not sites of long-term deposition; additionally, our results imply that glacial flow may govern the large-scale evolution of the North PLD, especially at high obliquity.     

Chaos formation by sublimation of volatile-rich substrate: Evidence from Galaxias Chaos, Mars

Galaxias Chaos deviates significantly from other chaotic regions due to the lack of associated outflow channels, lack of big elevation differences between the chaos and the surrounding terrain and due to gradual trough formation.A sequence of troughs in different stages is observed, and examples of closed troughs within blocks suggest that the trough formation is governed by a local stress field rather than a regional stress field. Moreover, geomorphic evidence suggests that Galaxias Chaos is capped by Elysium lavas, which superpose an unstable subsurface layer that causes chaotic tilting of blocks and trough formation.Based on regional mapping we suggest a formation model, where Vastitas Borealis Formation embedded between Elysium lavas is the unstable subsurface material, because gradual volatile loss causes shrinkage and differential substrate movement. This process undermines the lava cap, depressions form and gradually troughs develop producing a jigsaw puzzle of blocks due to trough coalescence.Observations of chaos west of Elysium Rise indicate that this process might have been widespread along the contact between Vastitas Borealis Formation and Elysium lavas. However, the chaotic regions have probably been superposed by Elysium/Utopia flows to the NW of Elysium Rise, and partly submerged with younger lavas to the west.     

Numerical simulation of the winter polar wave clouds observed by Mars Global Surveyor Mars Orbiter Laser Altimeter

In 1998, the Mars Orbiter Laser Altimeter revealed the presence of isolated or quasi-periodic thick clouds during the martian polar night. They are believed to be composed of CO₂ ice particles and to be tilted against the wind direction, a feature characteristic of vertically propagating orographic gravity waves. To support that interpretation, we present here numerical simulations with a two-dimensional anelastic model of stratified shear flow that includes simple CO₂ ice microphysics. In some of the simulations presented, the orography is an idealized trough, with dimensions characteristic of the many troughs that shape the Mars polar cap. In others, it is near the real orography. In the polar night conditions, our model shows that gravity waves over the north polar cap are strong enough to induce adiabatic cooling below the CO₂ frost point. From this cooling, airborne heterogeneous nucleation of CO₂ ice particles occurs from the ground up to the altitude of the polar thermal inversion. Although the model predicts that clouds can be present above 15 km, only low altitude clouds can backscatter the Laser beams of MOLA at a detectable level. Accordingly, the shape of the Laser echoes is related to the shape of the clouds at low level, but do not necessarily coincide with the top of the clouds. The model helps to interpret the cloud patterns observed by MOLA. Above an isolated orographic trough, an isolated extended sloping cloud tilted against the wind is obtained. The model shows that the observed quasi-periodic clouds are due to the succession of small-scale topographic features, rather than to the presence of resonant trapped lee waves. Indeed, the CO₂ condensation greatly damps the buoyancy force, essential for the maintenance of gravity waves far from their sources. Simulations with realistic topography profiles show the cloud response is sensitive to the wind direction. When the wind is directed upslope of the polar cap, on the one hand, a large scale cloud, modulated by small-scale waves, forms just above the ground. On the other hand, when the wind is directed downslope, air is globally warmed, and periodic ice clouds induced by small-scale orography form at altitudes higher than 3-5 km above the ground. In both cases, a good agreement between the simulated echoes and the observed one is obtained. According to our model, we conclude that the observed clouds are quasi-stationary clouds made of moving ice particles that successively grow and sublimate by crossing cold and warm phases of orographic gravity waves generated by the successive polar troughs. We also find that the rate of ice precipitation is relatively weak, except when there is a large scale air dynamical cooling.     

Stratigraphic analysis of the northern polar layered deposits of Mars: Implications for recent climate history

We examine the stratigraphy of the polar layered deposits (PLD) within the north polar cap of Mars to assess its layer continuity, correlations, cyclicity and structure and implications for the recent climate record. PLD sequences characterized using Fourier analysis and curve shape matching algorithms show that layers correlate throughout the upper part of the PLD. We tested for cyclicity and found that the uppermost ∼300 m contain a dominant wavelength layer packet of ∼30 m, interpreted to be a climate signal related to the 51 kyr precession cycle. Directly below this region we document a section of polar layered deposits ∼100 m thick without a dominant periodic signal; this is interpreted to represent a phase of low net accumulation and lag deposits formed during the last ice age, about 0.5-2 Ma ago. We further analyzed layer structure by combining these results with three-dimensional determinations of layer orientation (strike and dip) to assess the internal stratigraphy of the PLD and its implications for polar history. We show that individual layers within the PLD stratigraphy are not horizontal (no dip) but rather show broad variation in elevation with distance. Correlations suggest that the layer strikes and dips broadly follow present cap surface topography. Local variations in layer orientations in the vicinity of the troughs suggest that (1) trough structures were present at the time of layer accumulation and (2) dips may have been influenced by ice flow and/or static ice accumulation in the presence of preexisting troughs. This new information favors models in which the troughs are long-term structures of the PLD rather than (1) recently eroded into the PLD, or (2) very active and laterally migrating around the PLD. Our results strongly support the hypothesis that significant volumes of polar volatiles are mobilized and transported equator-ward during periods of increased obliquity. Our results predict that the upper ∼300 m of the north polar PLD accumulated in the last 500 ka, yielding net accumulation rates of ∼0.06 cm/yr. The presence and albedo of the no periodic signal zone suggest that polar net accumulation rates are very low and that dust rich lag deposits form during periods of sustained high obliquity. Layer sequences in the south polar and equatorial regions are examined and compared to those in the north; rhythmic sequences are observed in both regions but no direct correlations to the dominant signals of the north polar deposits have yet been found. These new techniques and observations provide a paradigm for further analysis of recent polar history (the upper kilometer of the record) and a basis for extending assessments to the lower part of the polar deposits and to other cyclic deposits in the geological record of Mars.     

Thermal and topographic tests of Europa chaos formation models from Galileo E15 observations

Stereo topography of an area near Tyre impact crater, Europa, reveals chaos regions characterised by marginal cliffs and domical topography, rising to 100-200 m above the background plains. The regions contain blocks which have both rotated and tilted. We tested two models of chaos formation: a hybrid diapir model, in which chaos topography is caused by thermal or compositional buoyancy, and block motion occurs due to the presence of near-surface (1-3 km) melt; and a melt-through model, in which chaos regions are caused by melting and refreezing of the ice shell. None of the hybrid diapir models tested generate any melt within 1-3 km of the surface, owing to the low surface temperature. A model of ocean refreezing following melt-through gives effective elastic thicknesses and ice shell thicknesses of 0.1-0.3 and 0.5-2 km, respectively. However, for such low shell thicknesses the refreezing model requires implausibly large lateral density contrasts (50-100 kg m⁻³) to explain the elevation of the centres of the chaos regions. Although a global equilibrium ice shell thickness of ≈2 km is possible if Europa's mantle resembles that of Io, it is unclear whether local melt-through events are energetically possible. Thus, neither of the models tested here gives a completely satisfactory explanation for the formation of chaos regions. We suggest that surface extrusion of warm ice may be an important component of chaos terrain formation, and demonstrate that such extrusion is possible for likely ice parameters.     

Formation of a terraced fan deposit in Coprates Catena, Mars

We have studied a terraced fan deposit with unique characteristics located within a trough of Coprates Catena. The fan has an average length of 6.8 km, and is approximately 44 km² in area and 18 km³ in volume. The fan's broad contributing valley is approximately 35 km long and it noticeably increases in depth about 12.8 km before it intersects the trough, where a rounded knickpoint marks the transition between flat-floored upstream and V-shaped downstream cross-sections. A 14-km-long channel with no apparent source enters the contributing valley from the south. A much smaller sinuous channel has incised along a smaller V-shaped valley in the uppermost eastern portion of the fan deposit. We explored several possible origins for the terraced fan, including mass wasting, volcanic flow, alluvial fan, and delta. We propose that water sourced from volcanic melting of ice eroded and transported material along the contributing valley. This material was then deposited as a delta in a lake within the trough. The concentric terraces are most likely the result of shoreline or ice cover erosion during drops in lake level. A light-toned layered deposit to the east of the fan deposit along the floor of the trough may represent a sedimentary unit formed during the terminal stages of the lake. Although other terraced fans have been identified on Mars, the Coprates Catena fan is unique because it has many more terraces and its surface was incised by a channel and associated valley. The identification of several other valleys to the east suggests that volcanic melting of volatiles during the Hesperian Period created favorable conditions for water flow along the plains in this region.     

Martian Seasonal CO2 Ice in Polygonal Troughs in Southern Polar Region: Role of the Distribution of Subsurface H2O Ice

The Mars Orbiter Camera onboard the Mars Global Surveyor has obtained several images of polygonal features in the southern polar region. In images taken during the end of the southern spring, when the surrounding surface is free of the seasonal frost, CO₂ ice still appears to be present within the polygonal troughs. In Earth's polar regions, polygons such as these are indicative of water ice in the ground below. We analyzed the seasonal evolution of the thermal state and the CO₂ content of these features. Our 2-D model includes condensation and sublimation of the CO₂ ice, a self consistent treatment of the variations of the thermal properties of the regolith, and the seasonal variations of the local atmospheric pressure which we take from the results of a general circulation model. We find that the residence time of seasonal CO₂ ice in troughs depends not only on atmospheric opacity and albedo of the CO₂ ice, but also and most significantly on the distribution of water ice in the regolith. Optical properties of the atmosphere and surface CO₂ ice can be independently obtained from observations. To date this is not true about the distribution of water ice below the surface. Our analysis quantifies the dependence of the seasonal cycle of the CO₂ ice within the troughs on the assumed distribution of the water ice below the surface. We show that presence of water ice in the ground at a depth smaller than the depth of the troughs reduces winter condensation rate of CO₂ ice. This is due to higher heat flux conducted from the water ice rich regolith toward the facets of the troughs.     

The evolution of dust deposits in the Martian north polar region

The origin and evolution of two major eolian deposits of the Martian north polar region, the layered deposits and the debris mantle, are examined. Both apparently result from deposition of dust along with the seasonal CO₂ frost cap. Dust deposited onto the perennial ice is incorporated into the layered deposits, while dust deposited directly onto the surface becomes part of the debris mantle. Climatically induced fluctuation of the perennial ice margin has influenced the evolution of both units. Periodic exposure to the atmosphere has allowed erosion of curvillinear troughs in the surface of the layered deposits. Intervening periods of deposition may have resulted in gradual poleward migration of the trough forms, leaving behind sets of low-amplitude surface undulations in former trough locations. Advance and retreat of the perennial ice margin has also probably resulted in a fine interfingering of the layered deposits-debris mantle contract. Limited post-depositional stripping of the debris mantle has been accomplished by intense winds blowing outward from the pole.     

Putative ice flows on Europa: Geometric patterns and relation to topography collectively constrain material properties and effusion rates

Europa's surface exhibits numerous small dome-like and lobate features, some of which have been attributed to fluid emplacement of ice or slush on the surface. We perform numerical simulations of non-Newtonian flows to assess the physical conditions required for these features to result from viscous flows. Our simulations indicate that the morphology of an ice flow on Europa will be, at least partially, influenced by pre-existing topography unless the thickness of the flow exceeds that of the underlying topography by at least an order of magnitude. Three classes of features can be identified on Europa. First, some (possibly most) putative flow-like features exhibit no influence from the pre-existing topography such as ridges, although their thicknesses are generally on the same order as those of ridges. Therefore, flow processes probably cannot explain the formation of these features. Second, some observed features show modest influence from the underlying topography. These might be explained by ice flows with wide ranges of parameters (ice temperatures >230 K, effusion rates >10⁷ m³ year⁻¹, and a wide range of grain sizes), although surface uplift (e.g., by diapirism) and in situ disaggregation provide an equally compelling explanation. Third, several observed features are completely confined by pre-existing topographic structures on at least one side; these are the best known candidates for flow features on Europa. If these features resulted from solid-ice flows, then temperatures >260 K and grain sizes <2 μm are required. Such small grain sizes seem unlikely; low-viscosity flows such as ice slurries or brines provide a better explanation for these features. Our results provide theoretical support for the view that many of Europa's lobate features have not resulted from solid-ice flows.     

Viscous relaxation on comets

Three Jupiter family comets have now been observed by spacecraft with the surprising result that these comets lack unambiguous impact craters. Large-scale topography generally appears to be softened on these comets, although sharp topography is preserved at small scales. We find that viscous relaxation of water ice may explain these observations, given reasonable assumptions about ice grain size and temperatures attained in the interiors. We suggest that both the shapes and the cratering records of Jupiter family comets may be substantially modified and no longer reflect cometary formation processes or collisional evolution in the Kuiper Belt.     

Amazonian northern mid-latitude glaciation on Mars: A proposed climate scenario

Recent geological observations in the northern mid-latitudes of Mars show evidence for past glacial activity during the late Amazonian, similar to the integrated glacial landsystems in the Dry Valleys of Antarctica. The large accumulation of ice (many hundreds of meters) required to create the observed glacial deposits points to significant atmospheric precipitation, snow and ice accumulation, and glacial flow. In order to understand the climate scenario required for these conditions, we used the LMD (Laboratoire de Météorologie Dynamique) Mars GCM (General Circulation Model), which is able to reproduce the present-day water cycle, and to predict past deposition of ice consistent with geological observations in many cases. Prior to this analysis, however, significant mid-latitude glaciation had not been simulated by the model, run under a range of parameters.In this analysis, we studied the response of the GCM to a wider range of orbital configurations and water ice reservoirs, and show that during periods of moderate obliquity (? = 25-35°) and high dust opacity (τ_dust = 1.5-2.5), broad-scale glaciation in the northern mid-latitudes occurs if water ice deposited on the flanks of the Tharsis volcanoes at higher obliquity is available for sublimation. We find that high dust contents of the atmosphere increase its water vapor holding capacity, thereby moving the saturation region to the northern mid-latitudes. Precipitation events are then controlled by topographic forcing of stationary planetary waves and transient weather systems, producing surface ice distribution and amounts that are consistent with the geological record. Ice accumulation rates of ∼10 mm yr⁻¹ lead to the formation of a 500-1000 m thick regional ice sheet that will produce glacial flow patterns consistent with the geological observations.     

Unstable extension of Enceladus' lithosphere

Regions near Enceladus' equator, Sarandib and Diyar Planitia, contain extensive sets of parallel ridges and troughs that may be diagnostic of the region's formation conditions. We present photoclinometry profiles across these ridges and troughs, which indicate that they are periodic, low-slope features with dominant wavelengths of 3 to 4 km and amplitudes between 100 and 400 m. The morphology of these terrains is consistent with formation via unstable extension of the lithosphere. Our numerical modeling demonstrates that unstable extension can generate large-scale topography under Enceladus-like conditions. Comparison of our photoclinometry profiles with the dominant wavelengths produced by our numerical model permits estimation of the background heat flow at the time the Sarandib-Diyar province formed. We estimate heat flows of 110 to , suggesting that resurfacing of the planitiae was accompanied by strong, localized heating. The extension necessary to produce the ridges and troughs may have been caused by now-inactive diapirs, internal phase changes, or other mechanisms. Our heat flux estimates imply elastic thickness at the time of resurfacing of 0.4 to 1.4 km, which are sufficient to have allowed satellite reorientation if the province was underlain by a low-density region. It is therefore plausible that Enceladus has experienced multiple heating events, each leading to localized resurfacing and global reorientation.     

Small-scale trench in the north polar region of Mars: Evolution of surface frost and ground ice concentration

In this paper we attempt to answer the question, how formation of a small-scale trench in the martian regolith affects local distribution of the subsurface ice. We are especially interested in the consequences of digging a trench to search for buried ice, as has been done during the Phoenix Mars Lander mission. However, the results may be also applicable for natural troughs, or cracks. We present results of simulations of diurnal exchange of water between the regolith and the atmosphere. Our model includes the heat and vapor migration in the regolith surrounding the trench, as well as formation of diurnal frost. We take into account scattering of light in the atmosphere and on the trench facets, as well as changes of atmospheric humidity on diurnal and seasonal time scales. Our calculations show, that the measurements of ice content in a sample obtained within one, or two days from the beginning of digging should not be affected. However, on somewhat longer time scale at the south facing site of the trench the regolith can be significantly depleted from ice. This effect should be taken into account if the excavation and taking samples from different depths will be performed in stages separated in time by a month, or more.     

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.     

Volcano/ground ice interactions in Elysium Planitia,Mars

The occurence within Elysium Planitia of meltwater deposits, possible pseudocraters, collapse features within troughs, and outflow channels indicates that a layer of subsurface volatiles existed at the time of volcanic activity within this area. The pseudocraters are interpreted to be indicators of near-surface volatiles, while meltwater deposits and the degree of preservation of trough walls and floors are thought to signify greater volatile depths. A latitudinal variation in the distribution of these features indicates either that the depth to the volatile layer increased from less than about 50 m at 35°N to greater than 600 m at 24°N, or that an ice wedge that existed at 35°N thinned to nonexistence at 24°N. Braided distributary channel systems within the chaotic terrain north of Elysium Planitia show that ephemeral lakes were repeatedly created and drained at this locality. The existence of volatiles contemporaneous with volcanic activity permits a search to be made for explosively generated landforms predicted to exist by previous theoretical models. Morphological evidence for strombolian, vulcanian and plinian eruptions is lacking within western Elysium Planitia; there are no identifiable cinder cones, pyroclastic flow deposits, or mantled areas indicative of large airfall deposits at an image resolution of 50–150 m/pixel. However, the pseudocraters indicate that small-scale phreatomagmatic activity may have taken place.