Distribution and state of H2O in the high-latitude shallow subsurface of mars

A quantitative model of the state, distribution, and migration of water in the shallow Martian regolith is presented. Reported results are confined to the region of the planet greater than 40° lat. The calculations take into account (1) expected thermal variations at all depths, latitudes, and times resulting from seasonal and astronomically induced insolation variations; (2) variations in atmospheric P_H2O and P_CO2 resulting from polar insolation variations and regolith adsorptive equilibria; (3) feedback effects related to latent heat and albedo variations resulting from condensation of atmospheric constituents; (4) two possible regolith mineralogies; (5) variable total H₂O content of the regolith; (6) kinetics of H₂O transport through the Martian atmosphere and regolith; and (7) equilibrium phase partitioning of H₂O between the condensed, adsorbed, and vapor phases. Results suggest that the adsorptive capacity of the regolith is important in controlling the state and distribution of high-latitude H₂O; unweathered mafic silicates favor the development of shallow ground ice at all temperate and polar latitudes, while heavily weathered clay-like regolith materials leads to a deeper ground ice interface and far more extensive quantities of adsorbed H₂O. The capacity of the high-latitude regolith for storage of H₂O and the total mass of H₂O exchanged between the atmosphere, polar cap, and subsurface over an obliquity cycle is found to be relatively independent of mineralogy. The maximum exchanged volume is found to be 3.0 × 10⁴ km³ of ice per cycle. Implications for the history of the polar caps and the origin of the layered terrain are discussed. Results also suggest that seasonal thermal waves act to force adsorbed H₂O into the solid phase over a wide variety of latitude/obliquity conditions. Seasonal phase cycling of regolith H₂O is most common at high latitudes and obliquities. Such phase behavior is highly dependent on regolith mineralogy. In a highly weathered regolith, adsorbed H₂O is annually forced into the solid phase at all latitudes ≥40° at obliquities greater than approximately 25°. Seasonal adsorption-freezing cycles which are predicted here may produce geomorphologic signatures not unlike those produced by terrestrial freeze-thaw cycles.     

The astronomical theory of climatic change on Mars

We examine the response of Martian climate to changes in solar energy deposition caused by variations of the Martian orbit and obliquity. We systematically investigate the seasonal cycles of carbon dioxide, water, and dust to provide a complete picture of the climate for various orbital configurations. We find that at low obliquity (15°) the atmospheric pressure will fall below 1 mbar; dust storms will cease; thick permanent CO₂ caps will form; the regolith will release CO₂; and H₂O polar ice sheets will develop as the permafrost boundaries move poleward. At high obliquity (35°) the annual average polar temperature will increase by about 10°K, slightly desorbing the polar regolith and causing the atmospheric pressure to increase by not more than 10 to 20 mbar. Summer polar ground temperatures as high as 273°K will occur. Water ice caps will be unstable and may disappear as the equilibrium permafrost boundary moves equatorward. However, at high eccentricity, polar ice sheets will be favored at one pole over the other. At high obliquity dust storms may occur during summers in both hemispheres, independent of the eccentricity cycle. Eccentricity and longitude of perihelion are most significant at modest obliquity (25°). At high eccentricity and when the longitude of perihelion is close to the location of solstice hemispherical asymmetry in dust-storm generation and in polar ice extent and albedo will occur.The systematic examination of the relation of climate and planetary orbit provides a new theory for the formation of the polar laminae. The terraced structure of the polar laminae originates when eccentricity and/or obliquity variations begin to drive water ice off the dusty permanent H₂O polar caps. Then a thin (meters) layer of consolidated dust forms on top of a dirty, slightly thicker (tens of meters) ice sheet and the composite is preserved as a layer of laminae composed predominately of water ice. Because of insolation variation on slopes, a series of poleward- and equatorward-facing scarps are formed where the edges of the laminae are exposed. Independently of orbital variations, these scarps propagate poleward both by erosion of the equatorward slopes and by deposition on the poleward slopes. Scarp propagation resurfaces and recycles the laminae forming the distinctive spiral bands of terraces observed and provides a supply of water to form new permanent ice caps. The polar laminae boundary marks the furthest eqautorward extension of the permanent H₂O caps as the orbit varies. The polar debris boundary marks the furthest equatorward extension of the annual CO₂ caps as the orbit varies.The Martian regolith is now a significant geochemical sink for carbon dioxide. CO₂ has been irreversibly removed from the atmosphere by carbonate formation. CO₂ has also benn removed by regolith adsorption. Polar temperature increases caused by orbital variations are not great enough     

Regolith-atmosphere exchange of water and carbon dioxide on Mars: Effects on atmospheric history and climate change

Exchange of CO₂ and H₂O between the Mars regolith and the atmosphere-cap system plays an important role in governing the evolution of the martian atmosphere and the martian climate. Most of the exchangeable CO₂ (perhaps one or two orders of magnitude more than the atmospheric inventory) is currently adsorbed on the deep regolith, and can be “cryopumped” to a large quasipermanent CO₂ cap (not now present) during lowest Mars obliquity (θ). During the obliquity driven regolith-cap CO₂ exchange cycle, the atmospheric pressure varies harmonically between ~0.1 mb (lowest Θ) and ? 20 mb (highest Θ). The regolith buffer plays only a small or negligible role in the seasonal CO₂ pressure variations caused by atmosphere-cap exchange because adsorption greatly inhibits diffusion of the seasonal “pressure wave” into the regolith. In contrast, thermally driven H₂O seasonal exchange between the atmosphere and regolith appears to be in large part responsible for observed seasonal variations in the small atmospheric H₂O inventory. Long term exchange of H₂O may be dominated by transfer between the polar caps and ice in the regolith. Available and potential tests of regolith-atmospheric-cap volatile exchange models using ground-based and spacecraft-based techniques are discussed.     

A conceptual model for explanation of Albedo changes in Martian craters

We propose a conceptual model to interpret AM/PM high albedo events (HAEs) in crater interiors at the Martian seasonal polar caps. This model consists of two components: (1) a relatively permanent high-albedo water–ice body exposed in a crater interior and (2) a variable crater albedo in response to aerosol optical depth, dust contamination, and H₂O/CO₂ frost deposits or sublimes in four phases, based on temperature and solar longitude changes. Two craters (Korolev crater of fully exposed water–ice layer and ‘Louth’ crater of partially exposed water–ice layer) are used to demonstrate the model. This model explains the HAEs and their seasonal changes and suggests that many crater-like features formed in the last episodic advance of the polar ice cap in the last high obliquity period should have water–ice exposed or covered. For the AM-only HAEs craters, there seems no need of a water–ice layer to be fully exposed, but a subsurface water–ice layer (or ice-rich regolith) is a necessary condition.     

PFS-MEX observation of ices in the residual south polar cap of Mars

The Mars Express spacecraft has a highly inclined orbit around Mars and so has been able to observe the south pole of Mars in illuminated conditions at the end of the southern summer (L_s=330). Spectra from the planetary Fourier spectrometer (PFS) short wavelength (SW) channel were recorded over the permanent ice cap to study its composition in terms of CO₂ ice and H₂O ice. Models are fitted to the observed data, which include a spatial mixture of soil (not covered by ice) and CO₂ frost (with a specific grain size and a small amount of included dust and H₂O ice). Two different kinds of spectra were observed: those over the permanent polar cap with almost pure CO₂ ice, negligible water ice, no soil fraction required, and bright; and those over mixed terrain (at the edge of the cap or near troughs) containing a significant soil spatial fraction, more water ice and smaller CO₂ grain size. The amount of water ice given by fits to scaled albedo models is less than 10 ppm by weight. When using multi-stream reflectance models with the appropriate lighting geometry, the water amount must be 2-5 times greater than the albedo fit (less than 50 ppm). At the periphery of the residual polar cap, we found a region almost completely covered by water frost, modeled as a mixture of micron-sized and sub-mm sized grains. Our result using a granular mixture of micron-sized grains of water ice and dust with the CO₂ grains is different from the modeling of OMEGA polar cap observations using molecular mixtures.     

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.     

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.     

Explaining the redox imbalance between the H and O escape fluxes at Mars by the oxidation of methane

From a comparison between the different observations of Martian methane existing today, including the new TES methane maps (Fonti and Marzo, 2010), we show that all sets of data are globally consistent with each other, and that a well definite seasonal cycle of methane has been at work for at least 10 yr. With a simple model of the balance between the loss fluxes of H and O, using up-to-date values of the escape fluxes, we show that the long-standing enigma of the imbalance between H and O escape fluxes may be solved by assuming that the missing sink of oxygen is the oxidation of methane. If no H₂ is released together with CH₄, a good agreement is found between the present CH₄ flux and the value imposed by the balance between H and O escape fluxes, an average over the last ≈10³ yr. If H₂ is released together with CH₄, as expected if CH₄ originates in serpentinization, the average level of CH₄ during the last 10³ yr should have been at least ten times lower than the present one. The lack of present H₂ release could suggest a long-term storage of methane in the subsurface under the form of clathrates, whereas H₂ has been lost to the atmosphere shortly after being produced. We suggest that the thin layer of CO₂ ice covering the permanent southern polar cap could result from the release of methane since the end of the last obliquity transition (time scale: 1 Myr), at an average rate of 0.1 Mt yr^?1, consistent with the values derived from: (i) the present observations of methane (time scale: 10 yr), (ii) the estimate from the observed imbalance between the H and O escape fluxes (time scale: 1 kyr). If so, the present release of methane from subsurface clathrates would have acted at a similar rate since at least 3 Myr.     

Remnant Massifs of Layered Deposits at High Northern Latitudes of Mars

Ice deposits of the permanent northern polar cap of Mars exhibit a clearly expressed layered structure connected with consecutive accumulation of a solid phase of H₂O (ice) with different impurity fractions of dust. The observed structure of the northern polar cap deposits is direct evidence of the uneven accumulation of ice associated with climate variations on Mars over a long period of time. Outside the boundaries of the present northern polar cap, there are isolated remnant massifs of bright deposits with a clearly distinguishable layered structure. In this paper, we analyze the morphology, spatial distribution, and material composition of remnant massifs within the latitude belt of 60°–85°. The similarity in the structure and material composition found in the remnant massifs and the layered deposits of the northern polar cap of Mars apparently suggests that these formations were genetically connected in the past climatic epochs.     

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.     

Thick and thin models of the evolution of carbon dioxide on Mars

We present results from a new simulation code that accounts for the evolution of the reservoirs of carbon dioxide on Mars, from its early years to the present. We establish a baseline model parameter set that produces results compatible with the present (i.e., Patm?6.5 mbar with permanent CO₂ ice cap) for a wide range of initial inventories. We find that the initial inventory of CO₂ broadly determines the evolutionary course of the reservoirs of CO₂. The reservoirs include the atmosphere, ice cap, adsorbed CO₂ in the regolith, and carbonate rocks. We track the evolution of the free inventory: the atmosphere, ice cap and regolith. Simulations begin at 4.53 Gyr before present with a rapid loss of free inventory to space in the early Noachian. Models that assume a relatively small initial inventory (?5 bar) have pronounced minima in the free inventory of CO₂ toward the end of the Noachian. Under baseline parameters, initial inventories below ∼4.5 bar result in a catastrophic loss of the free inventory to space. The current free inventory would be then determined by the balance between outgassing, sputtering losses and chemical weathering following the end of the late bombardment. We call these “thin” models. They generically predict small current free inventories in line with expectations of a small present CO₂ ice cap. For “thick” models, with initial inventories ?5 bar, a surplus of 300-700 mbar of free CO₂ remains during the late-Noachian. The histories of free inventory in time for thick models tend to converge within the last 3.5 Gyr toward a present with an ice cap plus atmospheric inventory of about 100 mbar. For thick models, the convergence is largely due to the effects of chemical weathering, which draws down higher free inventories more rapidly than the low. Thus, thick models have ?450 mbar carbonate reservoirs, while thin models have ?200 mbar. Though both thick and thin scenarios can reproduce the current atmospheric pressure, the thick models imply a relatively large current CO₂ ice cap and thin models, little or none. While the sublimation of a massive cap at a high obliquity would create a climate swing of greenhouse warming for thick models, under the thin model, mean temperatures and pressures would be essentially unaffected by increases in obliquity.     

The effects of the martian regolith on GCM water cycle simulations

This paper describes General Circulation Model (GCM) simulations of the martian water cycle focusing on the effects of an adsorbing regolith. We describe the 10-layer regolith model used in this study which has been adapted from the 1-D model developed by Zent, A.P., Haberle, R.M., Houben, H.C., Jakosky, B.M. [1993. A coupled subsurface-boundary layer model of water on Mars. J. Geophys. Res. 98 (E2), 3319-3337, February]. Even with a 30-min timestep and taking into account the effect of surface water ice, our fully implicit scheme compares well with the results obtained by Zent, A.P., Haberle, R.M., Houben, H.C., Jakosky, B.M. [1993. A coupled subsurface-boundary layer model of water on Mars. J. Geophys. Res. 98 (E2), 3319-3337, February]. This means, however, that the regolith is not able to reproduce the diurnal variations in column water vapour abundance of up to a factor of 2-3 as seen in some observations, with only about 10% of the atmospheric water vapour column exchanging with the subsurface on a daily basis. In 3-D simulations we find that the regolith adsorbs water preferentially in high latitudes. This is especially true in the northern hemisphere, where perennial subsurface water ice builds up poleward of 60° N at depths which are comparable to the Odyssey observations. Much less ice forms in the southern high latitudes, which suggests that the water ice currently present in the martian subsurface is not stable under present conditions and is slowly subliming and being deposited in the northern hemisphere. When initialising the model with an Odyssey-like subsurface water ice distribution the model is capable of forcing the simulated water cycle from an arbitrary state close to the Mars Global Surveyor Thermal Emission Spectrometer observations. Without the actions of the adsorbing regolith the equilibrated water cycle is found to be a factor of 2-4 too wet. The process by which this occurs is by adsorption of water during northern hemisphere summer in northern mid and high latitudes where it remains locked in until northern spring when the seasonal CO₂ ice cap retreats. At this time the water diffuses out of the regolith in response to increased temperature and is returned to the residual water ice cap by eddie transport.     

Near-infrared spectral monitoring of Triton with IRTF/SpeX I: establishing a baseline for rotational variability

We present eight new 0.8 to 2.4 μm spectral observations of Neptune's satellite Triton, obtained at IRTF/SpeX during 2002 July 15-22 UT. Our objective was to determine how Triton's near-infrared spectrum varies as Triton rotates, and to establish an accurate baseline for comparison with past and future observations. The most striking spectral change detected was in Triton's nitrogen ice absorption band at 2.15 μm; its strength varies by about a factor of two as Triton rotates. Maximum N₂ absorption approximately coincides with Triton's Neptune-facing hemisphere, which is also the longitude where the polar cap extends nearest Triton's equator. More subtle rotational variations are reported for Triton's CH₄ and H₂O ice absorption bands. Unlike the other ices, Triton's CO₂ ice absorption bands remain nearly constant as Triton rotates. Triton's H₂O ice is shown to be crystalline, rather than amorphous. Triton's N₂ ice is confirmed to be the warmer, hexagonal, β N₂ phase, and its CH₄ is confirmed to be highly diluted in N₂ ice.     

Martian volatiles: Their degassing history and geochemical fate

Observations of Mars and cosmochemical considerations imply that the total inventory of degassed volatiles on Mars is 10² to 10³ times that present in Mars' atmosphere and polar caps. The degassed volatiles have been physically and chemically incorporated into a layer of unconsolidated surface rubble (a “megaregolith”) up to 2km thick. Tentative lines of evidence suggest a high concentration (～5g/cm²) of ⁴⁰ Ar in the atmosphere of Mars. If correct, this would be consistent with a degassing model for Mars in which the Martian “surface” volatile inventory is presumed identical to that of Earth but scaled to Mars' smaller mass and surface area. The implied inventory would be: (⁴⁰Ar) = 4g/cm², (H₂O) = 1 × 10⁵g/cm², (CO₂) = 7 × 10³g/cm², (N₂) = 3 × 10²g/cm², (Cl) = 2 × 10³g/cm², and (S) = 2 × 10²g/cm². Such a model is useful for testing, but differences in composition and planetary energy history may be anticipated between Mars and Earth on theoretical grounds. Also, the model demands huge regolith sinks for the volatiles listed.If the regolith were in physical equilibrium with the atmosphere, as much as 2 × 10⁴g/cm² of H₂O could be stored in it as hard-frozen permafrost, or 5 × 10⁴g/cm² if equilibrium with the atmosphere were inhibited. Spectral measurements of Martian regolith material and laboratory measurement of weathering kinetics on simulated regolith material suggest large amounts of hydrated iron oxides and clay minerals exist in the regolith; the amount of chemically bound H₂O could be from 1 × 10⁴ to 4 × 10⁴g/cm². In an Earth-analogous model, a 2 km mixed regolith must contain the following concentrations of other volatile-containing compounds by weight: carbonates = 1.5%, nitrates = 0·3%, chlorides = 0.6%, and sulfates = 0.1%. Such concentrations would be undetectable by current Earth-based spectral reflectance measurements, and (except the nitrates) formation of the “required” amounts of these compounds could result from interaction of adsorbed H₂O and ice with primary silicates expected on Mars. Most of the CO₂ could be physically adsorbed on the regolith.Thus, maximum amounts of H₂O and other volatiles which could be stored in the Mars regolith are marginally compatible with those required by an Earth-analogous model, although a lower atmospheric ⁴⁰Ar concentration and regolith volatile inventory would be easier to reconcile with observational constraints. Differences in the ratios of H₂O and other volatiles to ⁴⁰Ar between surface volatiles on the real Mars and on an Earth-analogous Mars could result from and reflect differences in bulk composition and time history of degassing between Mars and Earth. Models relating Viking-observable parameters, e.g., (⁴⁰Ar) and (³⁶Ar), to the time history and overall intensity of Mars degassing are given.     

Injection of dust into the martian atmosphere: Evidence from the viking gas exchange experiment

It is proposed that dust storms on Mars that develop during predawn hours may be triggered by a freeze/thaw dust injection process. The model is based on a phenomenon that was observed during the Viking Gas Exchange experiments on Mars, in which adsorbed gas was catastrophically desorbed from soil samples when exposed to humidification at ∼5°C. Similar conditions may develop at midlatitude locations on Mars near perihelion, and a similar humidification-driven desorption process might occur in the soil column. If soils are dampened during humikification, desorbed gases in confined pore spaces could possibly reach 8.6 bar. Diurnal freezing may possibly cause H₂O to crystallize within the pores, possibly producing cohesive soil failure, release of the trapped gas, and explosive injection of freeze-dried powdery overburden dust into the atmospheric column. The process could potentially occur at 5–20 cm depth, and the freeze/thaw dust injection event may initiate after 10:00 PM local time (20°S lat). Dust would be injected at velocities approaching 450 m sec⁻¹ and it would remain in the atmosphere for several hours before settling out. The plumes could potentially regenerate diurnally until the growing atmospheric dust load produced sufficient dampening of the diurnal thermal wave to prevent freeze/thaw. Seasonal replenishment of H₂O could potentially occur by upward migration from depth during the period between 150 and 475 sols after perihelion. The model was experimentally tested and the results were in good agreement with predictions, although a factor of 14 times more gas evolved from the laboratory samples than from the Viking samples. Most of the characteristics of the predawn storms could possibly be adequately explained by the freeze/thaw injection model, including (1) predawn onsets, (2) postperihelion seasonal occurence, (3) daily recurrence during the initial phases of the storms, and (4) generation of blue clouds (H₂O ice) at the storm sites. The process may possibly occur over widespread locations at midlatitudes during seasonal retreat of “tempofrost” from these latitude belts. Permanent low albedo features in these latitude belts may possibly be regions of preferential humidification-induced dust entrainment and net dust removal. The H₂O injected into the atmosphere may potentially be a major source of H₂SO₄ and HCl aerosols, which may possibly chemically react with the regolith to form soluble sulfate and chloride salts. Mg²⁺ may be preferentially depleted from the dust.     

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.     

Carbon dioxide: Adsorption on palagonite and partitioning in the Martian regolith

We report new CO₂ adsorption measurements on palagonites. These results are used together with earlier results on basalt and nontronite adsorption to derive a “generic” relationship which is valid to within a factor of 3 for likely mixtures of basalt and weathering products of basalt. The relationship involves only t, P_CO2, and the specific surface area, and is relatively insensitive to mineralogy. It is used to predict the distribution and exchange of CO₂ on Mars. We conclude: (1) One to two orders of magnitude more CO₂ is adsorbed on the regolith than is present in the atmosphere and cap. (2) Nonetheless, most of the initially degassed CO₂ must have been lost to space or must be present as carbonates, especially if there was enough degassed CO₂ to provide a significant early greenhouse effect. (3) Given the derived relationship, the CO₂ vapor pressure curve, and the constraint that the system exhibits the current P_CO2 at the current obliquity, it is possible to predict approximately the atmospheric pressure at any obliquity (with or without a cap) without knowing the total available CO₂ inventory, the regolith mass, the regolith distribution, or its mineralogy, any better than those parameters are currently known.     

High latitude patterned grounds on Mars: Classification, distribution and climatic control

Patterned grounds such as polygonal features and slope stripes are the signature of the presence of ground ice and of temperature variations in cold regions on Earth. Identifying similar features on Mars is important to understand its past climate as well as the ground ice distribution. In this study, young patterned grounds are classed and mapped from the systematical analysis of Mars Observer Camera high resolution images. These features are located poleward of 55° latitude which fits the distribution of ground ice found by the Neutron Spectrometer onboard Mars Odyssey. Thermal contraction due to seasonal temperature variations is the predominant process of formation of polygons formed by cracks which sizes vary from 15 to 300 m. The small (<40 m) widespread polygons are very recent and degraded by the desiccation of ground ice from the cracks which enhances the effect of ice sublimation. The large polygons (50 to 300 m) located only around the south CO₂ polar cap indicate the presence of ground ice and thus outline the limit of the CO₂ ice cap. They could be due to the blanketing of water ice deposits by the advances and retreats of the residual CO₂ ice cap during the last thousand years. Large (50-250 m) and homogeneous polygons similar to ice wedge polygons, hillslope stripes and solifluction lobes may indicate that specific environments such as crater floors and hillslopes could have been submitted to freeze-thaw cycles, possibly related to higher summer temperatures in periods of obliquity higher than 35°. These interpretations must be strengthened by higher resolution images such as those of the HiRise mission of the Mars Reconnaissance Orbiter because locations with past seasonal thaw could be of major interest as potential landing sites for the Phoenix mission.     

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.     

Constraints on the composition of the martian south polar cap from gravity and topography

The polar caps of Mars have long been acknowledged to be composed of unknown proportions of water ice, solid CO₂ (dry ice), and dust. Gravity and topography data are here analyzed over the southern cap to place constraints on its density, and hence composition. Using a localized spectral analysis combined with a lithospheric flexure model of ice cap loading, the best fit density of the volatile-rich south polar layered deposits is found to be 1271 kg m⁻³ with 1-σ limits of 1166 and 1391 kg m⁻³. The best fit elastic thickness of this geologically young deposit is 140 km, though any value greater than 102 km can fit the observations. The best fit density implies that about 55% dry ice by volume could be sequestered in these deposits if they were completely dust free. Alternatively, if these deposits were completely free of solid CO₂, the dust content would be constrained to lie between about 14 and 28% by volume. The bulk thermal conductivity of the polar cap is not significantly affected by these maximum allowable concentrations of dust. However, even if a moderate quantity of solid CO₂ were present as horizontal layers, the bulk thermal conductivity of the polar cap would be significantly reduced. Reasonable estimates of the present day heat flow of Mars predict that dry ice beneath the thicker portions of the south polar cap would have melted. Depending on the quantity of solid CO₂ in these deposits today, it is even possible that water ice could melt where the cap is thickest. If independent estimates for either the dust or CO₂ content of the south polar cap could be obtained, and if radar sounding data could determine whether this polar cap is presently experiencing basal melting or not, it would be possible to use these observations to place tight constraints on the present day heat flow of Mars.