Spatial inhomogeneity of the martian subsurface water distribution: implication from a global water cycle model

In an effort to test and to understand the global hydrogen distribution in the shallow subsurface of Mars retrieved by the Mars Odyssey gamma-ray spectrometer, the present state and movement of water are investigated by a coupled global subsurface-atmosphere water cycle model. It was found that the observed global subsurface hydrogen distribution is largely consistent with the modeled global water cycle, so a large fraction of hydrogen is likely to exist as water, at low and mid latitudes in the form of adsorbed water. Under the present climate the water content in the shallow subsurface becomes higher in the northern hemisphere than in the southern hemisphere as a result of global water cycle, regardless of the initial water distribution in the soil or adsorptive capacity. The higher annual maximum soil temperature in the south, stronger net northward transport of atmospheric water vapor, and the emission of vapor from the northern residual polar cap in northern summer contribute to this hemispheric asymmetry. The generally higher adsorptive capacity of clay minerals in the northern plains may further increase this bias. The longitudinal inhomogeneity is caused by several factors, such as thermal inertia, adsorptive capacity, and atmospheric surface pressure. The water abundance is locally high in low thermal inertia regions (e.g., Arabia Terra) and at deep places where the surface pressure is high (e.g., Hellas); it is low in soil with a low adsorptive capacity (e.g., Tharsis) and high thermal inertia regions (e.g., Solis Planum). Most of the soil humidity near the surface at low and mid latitudes may originate from the atmosphere. The model implies that the upper soil layer should be largely ice-free because otherwise an excessive sublimation and vapor emission into the atmosphere in warm seasons would violate the observational constraints. Moreover, the more uniform latitudinal variation of the observed hydrogen abundance near the surface compared to that of deeper layers is indicative of the presence of adsorbed water instead of ground ice because the adsorbed water content does not as steeply depend on latitude as the ground ice stability. Concerning the regolith mineralogy, montmorillonite can much better account for the observed water cycle than palagonite. While the presence of permanent ground ice appears likely in the polar region below a thin layer, large seasonal cycle of phase change between pore ice and adsorbed water may be possible. Regolith adsorption/desorption is neither negligible nor crucial for the seasonal atmospheric water cycle, but the surface-atmosphere coupling is a major prerequisite for the long-term evolution of subsurface water distribution.     

Recent low-latitude freeze-thaw on Mars

Outside polar latitudes, features corresponding to surface thaw have yet to be identified on Mars. The youthful gully landforms observed at mid-high latitude [Malin, M., Edgett, K., 2000. Science 288, 2330-2335] are the nearest candidate, but the source (and nature) of the gully carving agent remains controversial [e.g., Musselwhite, D.S., Swindle, T.D., Lunine, J.I., 2001. Geophys. Res. Lett. 28, 1283-1285; Mellon, M.T., Phillips, R.J., 2001. J. Geophys. Res. 106, 1-15; Knauth, L.P., Burt, D.M., 2002. Icarus 158, 267-271; Costard, F., Forget, F., Mangold, N., Peulvast, J.P., 2002. Science 295, 110-113; Christensen, P.R., 2003. Nature 422, 45-48; Treiman, A.H., 2003. J. Geophys. Res. 108]. At higher obliquity than the present epoch, near-surface ground ice should be present globally [Mellon, M.T., Jakosky, B.M., 1995. J. Geophys. Res. 100 (E6), 11781-11799], populated by condensation of atmospheric water vapour in the top few metres of the regolith, or emplaced as dusty ice sheets reaching down towards the equator. The latitudinal restriction of these gullies to regions poleward of ±30° appears to argue against a thaw component to their formation—since ground ice is present and stable at all latitudes at high obliquity, the current (low) obliquity regime should result in ground ice thaw at low latitudes, where insolation and daytime temperatures are currently greatest, and this is not observed. A previously undescribed meltwater sequence in the Cerberus plains, at 20° N/187° E, shows that comparable, but much more continuous, and mappable melting and surface runoff have occurred in the geologically recent past at near-equatorial latitudes on Mars. Polygonal ground in the Cerberus plains is seen by the Mars Global Surveyor Mars Orbiter Camera (MOC) to suffer sequential, regional-scale volatile-loss consistent with thaw of near-surface ground ice under periglacial conditions. This degradation is continuously sampled by a single MOC strip, showing an icy landscape undergoing thaw modification and collapse, and may form the first evidence of equatorial wet-based glaciation during late Amazonian time, with indications of melting within the last million years. The dissolution and re-formation of polygonal ground links this landform to freeze-thaw processes, providing the conclusion to a question that has been the subject of debate for three decades—whether Mars' polygonal grounds require ice to form—and a consistent explanation for the fate of the water that carved the great outflow channels, much of which may still reside as ground ice in the regolith. This thaw occurs in the Cerberus Formation; deposits that are considered to be magmatic in origin, and the type formation for late-stage, “plains-style” volcanism on Mars [Keszthelyi, L., McEwen, A.S., Thordarson T., 2000. J. Geophys. Res. 105, 15027-15049]. By superposing large numbers of small impact craters, polygonal ground in the Cerberus plains sustains previous suggestions of a non-magmatic origin for this and other landforms in the region [Page, D.P., Murray, J.B., 2006. Icarus 183, 46-54]. Together, these periglacial landforms document evidence of climate change much younger than is currently recognised by crater counts, with important implications for age constraints on young surfaces and absolute age determinations by this method. It is tentatively suggested that this melting may be occurring today, with a striking correspondence between permafrost thaw in the Cerberus plains, the high atmospheric methane flux currently observed over this region [Mumma, M.J., Novak, R.E., DiSanti, M.A., Bonev, B.P., Dello Russo, N., 2004. Bull. Am. Astron. Soc. 36, 1127; Krasnopolsky, V.A., Maillard, J.P., Owen, T.C., 2004. Icarus 172, 537-547; Formisano, V., Atreya, S., Encrenaz, T., Ignatiev, N., Giuranna, M., 2004. Science 306, 1758-1761], and the only latitude zone on Mars—equatorward of 30° N—where melting of ground ice is thought possible in the current climate [Haberle, R.M., McKay, C.P., Schaeffer, J., Cabrol, N.A., Grin, E.A., Zent, A.P., Quinn, R., 2001. J. Geophys. Res. 106 (E10), 23317-23326; Lobitz, B., Wood, B.L., Averner, M.M., McKay, C.P., 2001. Proc. Natl. Acad. Sci. 98, 2132-2137]. Low-latitude polygonal ground as transient, and hydrologically active over wide areas transforms our understanding of the recent climatic evolution of Mars, supporting models of atmospheric water-ice migration [Mischna, M., Richardson, M.I., Wilson, R.J., McCleese, D.J., 2003. J. Geophys. Res. 108 (E6). 5062], complex, volatile stratigraphies [Clifford, S.M., Parker, T.J., 2001. Icarus 154, 40-79], and hypothesised, geologically recent ‘ice ages’ [Head, J.W., Mustard, J.F., Kreslavsky, M.A., Milliken, R.E., Marchant, D.R., 2003. Nature 426, 797-802]. The temporal coincidence of glacial epochs on the Earth and Mars during the Quaternary and latest Amazonian would suggest a coupled system linking both [Sagan, C., Young, A.T., 1973. Nature 243, 459].     

A 1 Gyr climate model for Mars: new orbital statistics and the importance of seasonally resolved polar processes

New results from a 1 Gyr integration of the martian orbit are presented along with a seasonally resolved energy balance climate model employed to illuminate the gross characteristics of the long-term atmospheric pressure evolution. We present a new analysis of the statistical variation of the martian obliquity and precession prior to and subsequent to the formation of the Tharsis uplift, and explore the long term effects on the martian climate. We find that seasonal polar cycles have a critical influence on the ability for the regolith to release CO₂ at high obliquities, and find that the atmospheric CO₂ actually decreases at high obliquities due to the cooling effect of polar deposits at latitudes where seasonal caps form. At low obliquity, the formation of massive, permanent polar caps depends critically on the values of the frost albedo, A_frost, and frost emissivity, ?_frost. Using our model with values of A_frost=0.67 and ?_frost=0.55, matched to the NASA Ames General Circulation Model (GCM) results (Haberle et al., 1993, J. Geophys. Res. 98, 3093-3123, and Haberle et al., 2003, Icarus 161, 66-89), we find that permanent caps only form at low obliquities (<13°), suggesting that any permanent deposits on the surface of Mars today may be residuals left over from a period of very low obliquity, or are the result of mechanisms not represented by this model. Thus, contrary to expectations, the martian atmospheric pressure is remarkable static over time, and decreases both at high and low obliquity. Also, from our one billion year orbital model, we present new results on the fraction of time Mars is expected to experience periods of low obliquity and high obliquity.     

The evolution of exposed ice in a fresh mid-latitude crater on Mars

Recent observations of the surface of Mars have shown several fresh mid-latitude craters. Some of these craters show exposed ice (Byrne, S. et al. [2009]. Science 325, 1674-1676.). In some craters, albedo of ice slowly decreases, while in others, it remains nearly constant. We attempt to determine influence of the regolith structure on the rate of sublimation of ice. For this purpose we performed numerical simulations describing evolution of the exposed ice in model craters located at middle latitudes.We consider a new model for the structure and evolution of the material at- and beneath the crater floors. In contrast to the previous study by Dundas and Byrne (Dundas, C.M., Byrne, S. [2010]. Icarus 206, 716-728.) we do not investigate sublimation of dirty ice, and the related formation of a sublimation lag. Instead, we consider sublimation of a pure ice layer on top of layered regolith. In our model the observed reflectivity decreases due to the sublimation-driven changes of the optical properties of thinning clean ice. This offers an alternative to the deposition of the dust embedded in ice (sublimation lag).We have shown that in our model among many parameters affecting ice sublimation rate, volumetric fraction of water ice in the subsurface beneath the crater has the strongest influence. Hence observed darkening of the ice patch on the crater floor might be sufficient to determine the content of water ice in the subsurface. Our calculations show that an albedo decrease of fresh ice patches in mid-latitude craters can be explained by either strong dust sedimentation or, if this is excluded, by sublimation of a thin layer of water ice from the regolith with large thermal inertia. This is consistent with a large volumetric fraction of water ice beneath the crater floor and contributes to evidence for an extended subsurface water reservoir on Mars.The overall conclusion of our work is that a thin post-impact surface ice coating over ice-rich ground beneath the crater floors is consistent with the observations.     

THEMIS observations of Mars aerosol optical depth from 2002-2008

We use infrared images obtained by the Thermal Emission Imaging System (THEMIS) instrument on-board Mars Odyssey to retrieve the optical depth of dust and water ice aerosols over more than 3.5 martian years between February 2002 (MY 25, Ls=330°) and December 2008 (MY 29, Ls=183°). These data provide an important bridge between earlier TES observations and recent observations from Mars Express and Mars Reconnaissance Orbiter. An improvement to our earlier retrieval [Smith, M.D., Bandfield, J.L., Christensen, P.R., Richardson, M.I., 2003. J. Geophys. Res. 108, doi:10.1029/2003JE002114] to include atmospheric temperature information from THEMIS Band 10 observations leads to much improved retrievals during the largest dust storms. The new retrievals show moderate dust storm activity during Mars Years 26 and 27, although details of the strength and timing of dust storms is different from year to year. A planet-encircling dust storm event was observed during Mars Year 28 near Southern Hemisphere Summer solstice. A belt of low-latitude water ice clouds was observed during the aphelion season during each year, Mars Years 26 through 29. The optical depth of water ice clouds is somewhat higher in the THEMIS retrievals at ∼5:00 PM local time than in the TES retrievals at ∼2:00 PM, suggestive of possible local time variation of clouds.     

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.     

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.     

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.     

Ozone abundance on Mars from infrared heterodyne spectra: II. Validating photochemical models

Ozone is an important observable tracer of martian photochemistry, including odd hydrogen (HO_x) species important to the chemistry and stability of the martian atmosphere. Infrared heterodyne spectroscopy with spectral resolution ?¹⁰⁶ provides the only ground-based direct access to ozone absorption features in the martian atmosphere. Ozone abundances were measured with the Goddard Infrared Heterodyne Spectrometer and the Heterodyne Instrument for Planetary Wind and Composition at the NASA Infrared Telescope Facility on Mauna Kea, Hawai'i. Retrieved total ozone column abundances from various latitudes and orbital positions (LS=40°, 74°, 102°, 115°, 202°, 208°, 291°) are compared to those predicted by the first three-dimensional gas phase photochemical model of the martian atmosphere [Lefèvre, F., Lebonnois, S., Montmessin, F., Forget, F., 2004. J. Geophys. Res. 109, doi:10.1029/2004JE002268. E07004]. Observed and modeled ozone abundances show good agreement at all latitudes at perihelion orbital positions (LS=202°, 208°, 291°). Observed low-latitude ozone abundances are significantly higher than those predicted by the model at aphelion orbital positions (LS=40°, 74°, 115°). Heterogeneous loss of odd hydrogen onto water ice cloud particles would explain the discrepancy, as clouds are observed at low latitudes around aphelion on Mars.     

Observations of martian gullies and constraints on potential formation mechanisms: II. The northern hemisphere

The formation process(es) responsible for creating the observed geologically recent gully features on Mars has remained the subject of intense debate since their discovery. We present new data and analysis of northern hemisphere gullies from Mars Global Surveyor data which is used to test the various proposed mechanisms of gully formation. We located 137 Mars Orbiter Camera (MOC) images in the northern hemisphere that contain clear evidence of gully landforms and analyzed these images in combination with Mars Orbiter Laser Altimeter (MOLA) and Thermal Emission Spectrometer (TES) data to provide quantitative measurements of numerous gully characteristics. Parameters we measured include apparent source depth and distribution, vertical and horizontal dimensions, slopes, orientations, and present-day characteristics that affect local ground temperatures. Northern hemisphere gullies are clustered in Arcadia Planitia, Tempe Terra, Acidalia Planitia, and Utopia Planitia. These gullies form in craters (84%), knobby terrain (4%), valleys (3%), other/unknown terrains (9%) and are found on all slope orientations although the majority of gullies are equator-facing. Most gullies (63%) are associated with competent rock strata, 26% are not associated with strata, and 11% are ambiguous. Assuming thermal conductivities derived from TES measurements as well as modeled surface temperatures, we find that 95% of the gully alcove bases with adequate data coverage lie at depths where subsurface temperatures are greater than 273 K and 5% of the alcove bases lie within the solid water regime. The average alcove length is 470 m and the average channel length is 690 m. Based on a comparison of measured gully features with predictions from the various models of gully formation, we find that models involving carbon dioxide, melting ground ice in the upper few meters of the soil, dry landslide, and surface snowmelt are the least likely to describe the formation of the martian gullies. Although some discrepancies still exist between prediction and observation, the shallow and deep aquifer models remain as the most plausible theories. Interior processes involving subsurface fluid sources are generally favored over exogenic processes such as wind and snowfall for explaining the origin of the martian gullies. These findings gleaned from the northern hemisphere data are in general agreement with analyses of gullies in the southern hemisphere [Heldmann, J.L., Mellon, M.T., 2004. Icarus 168, 285-304].     

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.     

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.     

Water distribution in Martian permafrost regions from joint analysis of HEND (Mars Odyssey) and MOLA (Mars Global Surveyor) data

We jointly analyze data from the High-Energy Neutron Detector (HEND) onboard the NASA Mars Odyssey spacecraft and data from the Mars Orbiter Laser Altimeter (MOLA) onboard the Mars Global Surveyor spacecraft. The former instrument measures the content of hydrogen (in the form of H₂O or OH) in the subsurface layer of soil and the latter instrument measures the surface albedo with respect to the flux of solar energy. We have checked the presence of a correlation between these two data sets in various Martian latitude bands. A significant correlation has been found between these data at latitudes poleward of 40° in the northern hemisphere and at latitudes 40°–60° in the southern hemisphere. This correlation is interpreted as evidence for the presence of stable water ice in these regions under a dry layer of soil whose thickness is determined by the condition for equilibrium between the condensation of water from the atmosphere and its sublimation when heated by solar radiation. For these regions, we have derived an empirical relation between the flux of absorbed solar radiation and the thickness of the top dry layer. It allows the burial depth of the water ice table to be predicted with a sub-kilometer resolution based on near-infrared albedo measurements. We have found no correlation in the southern hemisphere at latitudes >60°, although neutron data also suggest that water ice is present in this region under a layer of dry soil. We conclude that the thickness of the dry layer in this region does not correspond to the equilibrium condition between the water ice table and the atmosphere.     

Sorted clastic stripes, lobes and associated gullies in high-latitude craters on Mars: Landforms indicative of very recent, polycyclic ground-ice thaw and liquid flows

Self-organised patterns of stone stripes, polygons, circles and clastic solifluction lobes form by the sorting of clasts from fine-grained sediments in freeze-thaw cycles. We present new High Resolution Imaging Science Experiment (HiRISE) images of Mars which demonstrate that the slopes of high-latitude craters, including Heimdal crater - just 25 km east of the Phoenix Landing Site - are patterned by all of these landforms. The order of magnitude improvement in imaging data resolution afforded by HiRISE over previous datasets allows not only the reliable identification of these periglacial landforms but also shows that high-latitude fluviatile gullies both pre- and post-date periglacial patterned ground in several high-latitude settings on Mars. Because thaw is inherent to the sorting processes that create these periglacial landforms, and from the association of this landform assemblage with fluviatile gullies, we infer the action of liquid water in a fluvio-periglacial context. We conclude that these observations are evidence of the protracted, widespread action of thaw liquids on and within the martian regolith. Moreover, the size frequency statistics of superposed impact craters demonstrate that this freeze-thaw environment is, at least in Heimdal crater, less than a few million years old. Although the current martian climate does not favour prolonged thaw of water ice, observations of possible liquid droplets on the strut of the Phoenix Lander may imply significant freezing point depression of liquids sourced in the regolith, probably driven by the presence of perchlorates in the soil. Because perchlorates have eutectic temperatures below 240 K and can remain liquid at temperatures far below the freezing point of water we speculate that freeze-thaw involving perchlorate brines provides an alternative “low-temperature” hypothesis to the freeze-thaw of more pure water ice and might drive significant geomorphological work in some areas of Mars. Considering the proximity of Heimdal crater to the Phoenix Landing Site, the presence of such hydrated minerals might therefore explain the landforms described here. If this is the case then the geographical distribution of martian freeze-thaw landforms might reflect relatively high temperatures (but still below 273 K) and the locally elevated concentration of salts in the regolith.     

O2 1.27 μm emission maps as derived from OMEGA/MEx data

We present a method to derive the 2D maps of the O₂ (a¹Δ_g) airglow emission at 1.27 μm from the OMEGA/MEx nadir observations. The OMEGA imaging capabilities allow monitoring the 2D distribution, daily and seasonal variation of the O₂ emission intensities with a detection limit of 4 MR. The highest values, of the order of ∼31 MR, are found on the south pole for 11 h < LT < 13 h, during the early spring (186° < Ls < 192°) of martian year (MY) 27, according to the Mars Year numbering scheme of Clancy et al. [Clancy, R.T., Wolff, M.J., Christensen, P.R., 2003. Mars aerosol studies with the MGS TES emission phase function observations: Optical depths, particle sizes, and ice cloud types versus latitude and solar longitude. J. Geophys. Res. 108. doi: 10.1029/2003JE002058]. In the polar regions the day-by-day variability, associated with polar vortex turbulences, is obtained of the order of 30-50% as predicted by the model [Lefévre, F., Lebonnois, S., Montmessin, F., Forget, F., 2004. Three-dimensional modeling of ozone on Mars. J. Geophys. Res. 109, E07004. doi: 10.1029/2004JE002268] and found by SPICAM [Perrier, S., Bertaux, J.-L., Lebonnois, S., Korablev, O., Fedorova, A., 2006. Global distribution of total ozone on Mars from SPICAM/MEX UV measurements. J. Geophys. Res. 111, E09S06. doi: 10.1029/2006JE002681]. In the considered set of data a maximum of the O₂ emission is observed between 11 h and 15 h LT in the latitude range 70-85° during early spring on both hemispheres, while for the southern autumn-winter season a maximum is found between 50° and 60° in the southern hemisphere for MY28. Increase of intensity of the O₂ emission observed from Ls 130° to 160° at southern high latitudes may be explained by increase of solar illumination conditions in the maps acquired during the considered period.Atmospheric waves crossing the terminator on the southern polar regions are observed for the first time during the MY28 early spring. The spatial scale of the waves ranges from 100 to 130 km, and the intensity fluctuations are of the order of 4MR.This study confirms the high potentiality of O₂ (a¹Δ_g) day glow as a passive tracer of the martian atmosphere dynamics at high latitudes.     

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.     

Further observations of regional dust storms and baroclinic eddies in the northern hemisphere of Mars

We have investigated the near-surface meteorology in the northern hemisphere of Mars through detailed analysis of data obtained with Mars Global Surveyor in January-August 2005. The season in the northern hemisphere ranged from midsummer through winter solstice of Mars Year (MY) 27. We examined composite, wide-angle images from the Mars Orbiter Camera and compiled a catalog of the dust storms that occurred in this interval. As in previous martian years, activity in the northern hemisphere was dominated by regional “flushing” dust storms that sweep southward through the major topographic basins, most frequently in Acidalia Planitia. We also used atmospheric profiles retrieved from radio occultation experiments to characterize eddy activity near the surface at high northern latitudes. There are strong correlations between the two sets of observations, which allowed us to identify three factors that influence the timing and location of the regional dust storms: (1) transitions among baroclinic wave modes, which strongly modulate the intensity of meridional winds near the surface, (2) storms zones, which impose strong zonal variations on the amplitude of some baroclinic eddies, and (3) stationary waves, which further modulate the wind field near the surface. The flushing dust storms ceased abruptly in midautumn, possibly in response to source depletion, CO₂ condensation, a shift in the period of the baroclinic eddies, and changes in the tidal wind field near the surface. Our results extend the meteorological record of the northern hemisphere, substantiate the findings of previous investigations, and further illuminate the climatic impact of baroclinic eddies.     

Investigations of the variability of dust particle sizes in the martian atmosphere using the NASA Ames General Circulation Model

We present a Mars General Circulation Model (GCM) numerical investigation of the physical processes (i.e., wind stress and dust devil dust lifting and atmospheric transport) responsible for temporal and spatial variability of suspended dust particle sizes. Measurements of spatial and temporal variations in airborne dust particles sizes in the martian atmosphere have been derived from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) spectral and emission phase function data [Wolff, M.J., Clancy, R.T., 2003. J. Geophys. Res. (Planets) 108 (E9), doi:10.1029/2003JE002057. 1-1; Clancy, R.T., Wolff, M.J., Christensen, P.R., 2003. J. Geophys. Res. (Planets) 108 (E9), doi:10.1029/2003JE002058. 2-1]. The range of dust particle sizes simulated by the NASA Ames GCM is qualitatively consistent with TES-derived observations of effective dust particle size variability. Model results suggest that the wind stress dust lifting scheme (which produces regionally confined dust lifting) is the process responsible for the majority of the dust particle size variability in the martian atmosphere. Additionally, model results suggest that atmospheric transport processes play an important role in the evolution of atmospheric dust particles sizes during substantial dust storms on Mars. Finally, we show that including the radiative effects of a spatially variable particle size distribution significantly influences thermal and dynamical fields during the dissipation phase of the simulated global dust storm.     

Methane adsorption on a martian soil analog: An abiogenic explanation for methane variability in the martian atmosphere

Recent observations suggest methane in the martian atmosphere is variable on short spatial and temporal scales. However, to explain the variability by loss reactions requires production rates much larger than expected. Here, we report results of laboratory studies of methane adsorption onto JSC-Mars-1, a martian soil simulant, and suggest that this process could explain the observations. Uptake coefficient (γ) values were measured as a function of temperature using a high-vacuum Knudsen cell able to simulate martian temperature and pressure conditions. Values of γ were measured from 115 to 135 K, and the data were extrapolated to higher temperatures with more relevance to Mars. Adsorptive uptake was found to increase at lower temperatures and larger methane partial pressures. Although only sub-monolayer methane surface coverage is likely to exist under martian conditions, a very large mineral surface area is available for adsorption as atmospheric methane can diffuse meters into the regolith. As a result, significant methane may be temporarily lost to the regolith on a seasonal time scale. As this weak adsorption is fully reversible, methane will be re-released into the atmosphere when surface and subsurface temperatures rise and so no net loss of methane occurs. Heterogeneous interaction of methane with martian soil grains is the only process proposed thus far which contains both rapid methane loss and rapid methane production mechanisms and is thus fully consistent with the reported variability of methane on Mars.     

Transverse Aeolian Ridges (TARs) on Mars II: Distributions, orientations, and ages

Transverse Aeolian Ridges (TARs), 10 m scale, ripple-like aeolian bedforms with simple morphology, are widespread on Mars but it is unknown what role they play in Mars’ wider sediment cycle. We present the results of a survey of all Mars Global Surveyor Narrow angle images in a pole-to-pole study area, 45° longitude wide.Following on from the classification scheme and preliminary surveys of Balme et al. (Balme, M.R., Berman, D.C., Bourke, M.C., Zimbelman, J.R. [2008a]. Geomorphology 101, 703-720) and Wilson and Zimbelman (Wilson, S.A., Zimbelman, J.R. [2004]. J. Geophys. Res. 109 (E10). doi:10.1029/2004JE002247) we searched more than 10,000 images, and found that over 2000 reveal at least 5% areal cover by TARs. The mean TAR areal cover in the study area is about 7% (3% in the northern hemisphere and 11% in the southern hemisphere) but TARs are not homogenously distributed - they are concentrated in the mid-low latitudes and almost absent poleward of 35°N and 55°S. We found no clear correlation between TAR distribution and any of thermal inertia, kilometer-scale roughness, or elevation. We did find that TARs are less common at extremes of elevation.We found that TARs are most common near the equator (especially in the vicinity of Meridiani Planum, in which area they have a distinctive “barchan-like” morphology) and in large southern-hemisphere impact craters. TARs in the equatorial band are usually associated with outcrops of layered terrain or steep slopes, hence their relative absence in the northern hemisphere. TARs in the southern hemisphere are most commonly associated with low albedo, intercrater dune fields. We speculate that the mid-latitude mantling terrain (e.g., Mustard, J.F., Cooper, C.D., Rifkin, M.K. [2001]. Nature 412, 411-414; Kreslavsky, M.A., Head, J.W. [2002]. J. Geophys. Res. 29 (15). doi:10.1029/2002GL015392) could also play a role in covering TARs or inhibiting saltation.We compared TAR distribution with general circulation model (GCM) climate data for both surface wind shear stress and wind direction. We performed GCM runs at various obliquity values to simulate the effects of changing obliquity on recent Mars climate. We found good general agreement between TAR orientation and GCM wind directions from present day obliquity conditions in many cases, but found no good correlation between wind shear stress and TAR distribution.We performed preliminary high resolution crater count studies of TARs in both equatorial and southern intracrater dunefield settings and compared these to superposition relationships between TARs and large dark dunes. Our results show that TARs near dunefield appear to be younger than TARs in the equatorial regions. We infer that active saltation from the large dunes keeps TARs active, but that TARs are not active under present day condition when distal to large dunes - perhaps supporting the interpretation that TARs are granule ripples.We conclude that local geology, rather than wind strength, controls TAR distribution, but that their orientation matches present-day regional wind patterns in most cases. We suggest that TARs are likely most (perhaps only) active today when they are proximal to large dark dune fields.