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.     

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.     

Ancient melting of mid-latitude snowpacks on Mars as a water source for gullies

We hypothesize that during past epochs of high obliquity seasonal snowfields at mid-latitudes melted to produce springtime sediment-rich surface flows resulting in gully formation. Significant seasonal mid-latitude snowfall does not occur on Mars today. General Circulation Model (GCM) results, however, suggest that under past climate conditions there may have been centimeters of seasonal mid-latitude snowfall [Mischna, M.A., Richardson, M.I., Wilson, R.J., McCleese, D.J., 2003. J. Geophys. Res. Planets 108, doi:10.1029/2003JE002051. 5062]. Gully locations have been tabulated by several researchers (e.g. [Heldmann, J.L., Mellon, M.T., 2004. Icarus 168, 285–304; Heldmann, J.L., Carlsson, E., Johansson, H., Mellon, M.T., Toon, O.B., 2007. Icarus 188, 324–344; Malin, M.C., Edgett, K.S., 2000. Science 288, 2330–2335]) and found to correspond to mid-latitude bands. A natural question is whether the latitudinal bands where the gullies are located correspond to areas where the ancient snowfalls may have melted, producing runoff which may have incised gullies. In this study we model thin snowpacks with thicknesses similar to those predicted by [Mischna, M.A., Richardson, M.I., Wilson, R.J., McCleese, D.J., 2003. J. Geophys. Res. Planets 108, doi:10.1029/2003JE002051. 5062]. We model these snowpacks under past climate regimes in order to determine whether snowmelt runoff could have occurred, and whether significant amounts of warm soil (T>273 K) existed on both poleward and equatorward slopes in the regions where gullies exist. Both warm soil and water amounts are modeled because soil and water may have mixed to form a sediment-rich flow. We begin by applying the snowpack model of Williams et al. [Williams, K.E., Toon, O.B., Heldmann, J.E., Mellon, M., 2008. Icarus 196, 565–577] to past climate regimes characterized by obliquities of 35° (600 ka before present) and 45° (5.5 ma before present), and to all latitudes between 70° N and 70° S. We find that the regions containing significant snowmelt runoff correspond to the regions identified by Heldmann and Mellon [Heldmann, J.L., Mellon, M.T., 2004. Icarus 168, 285–304], Heldmann et al. [Heldmann, J.L., Carlsson, E., Johansson, H., Mellon, M.T., Toon, O.B., 2007. Icarus 188, 324–344] and Malin and Edgett [Malin, M.C., Edgett, K.S., 2000. Science 288, 2330–2335] as containing large numbers of gullies. We find that the snowmelt runoff (>1 mm, with equivalent rainfall rates of 0.25 mm/h) and warm soil (>1 cm depth) would have occurred on slopes within the gullied latitudinal bands. The snowfall amounts modeled are predicted to be seasonal [Mischna, M.A., Richardson, M.I., Wilson, R.J., McCleese, D.J., 2003. J. Geophys. Res. Planets 108, doi:10.1029/2003JE002051. 5062], and our modeling finds that under the previous climate regimes there would have been meltwater present on the slopes in question for brief periods of time, on the order of days, each year. Our model provides a simple explanation for the latitudinal distribution of the gullies, and also suggests that the gullies date to times when water migrated away from the present poles to the mid-latitudes.     

Principal components analysis of Mars in the near-infrared

Principal components analysis and target transformation are applied to near-infrared image cubes of Mars in a study to disentangle the spectra into a small number of spectral endmembers and characterize the spectral information. The image cubes are ground-based telescopic data from the NASA Infrared Telescope Facility during the 1995 and 1999 near-aphelion oppositions when ice clouds were plentiful [ [Clancy et al., 1996] and [56]], and the 2003 near-perihelion opposition when ice clouds are generally limited to topographically high regions (volcano cap clouds) but airborne dust is more common [Martin, L.J., Zurek, R.W., 1993. J. Geophys. Res. 98 (E2), 3221-3246]. The heart of the technique is to transform the data into a vector space along the dimensions of greatest spectral variance and then choose endmembers based on these new “trait” dimensions. This is done through a target transformation technique, comparing linear combinations of the principal components to a mineral spectral library. In general Mars can be modeled, on the whole, with only three spectral endmembers which account for almost 99% of the data variance. This is similar to results in the thermal infrared with Mars Global Surveyor Thermal Emission Spectrometer data [Bandfield, J.L., Hamilton, V.E., Christensen, P.R., 2000. Science 287, 1626-1630]. The globally recovered surface endmembers can be used as inputs to radiative transfer modeling in order to measure ice abundance in martian clouds [Klassen, D.R., Bell III, J.F., 2002. Bull. Am. Astron. Soc. 34, 865] and a preliminary test of this technique is also presented.     

Stratigraphical evidence of late Amazonian periglaciation and glaciation in the Astapus Colles region of Mars

Recent modeling of the meteorological conditions during and following times of high obliquity suggests that an icy mantle could have been emplaced in western Utopia Planitia by atmospheric deposition during the late Amazonian period [Costard, F.M., Forget, F., Madeleine, J.B., Soare, R.J., Kargel, J.S., 2008. Lunar Planet. Sci. 39. Abstract 1274; Madeleine, B., Forget, F., Head, J.W., Levrard, B., Montmessin, F., 2007. Lunar Planet. Sci. 38. Abstract 1778]. Astapus Colles (ABa) is a late Amazonian geological unit — located in this hypothesized area of accumulation — that comprises an icy mantle tens of meters thick [Tanaka, K.L., Skinner, J.A., Hare, T.M., 2005. US Geol. Surv. Sci. Invest., Map 2888]. For the most part, this unit drapes the early Amazonian Vastitas Borealis interior unit (ABv_i); to a lesser degree it overlies the early Amazonian Vastitas Borealis marginal unit (ABv_m) and the early to late Hesperian UP plains unit HBu₂ [Tanaka, K.L., Skinner, J.A., Hare, T.M., 2005. US Geol. Surv. Sci. Invest., Map 2888]. Landscapes possibly modified by late-Amazonian periglacial processes [Costard, F.M., Kargel, J.S., 1995. Icarus 114, 93-112; McBride, S.A., Allen, C.C., Bell, M.S., 2005. Lunar Planet. Sci. 36. Abstract 1090; Morgenstern, A., Hauber, E., Reiss, D., van Gasselt, S., Grosse, G., Schirrmeister, L., 2007. J. Geophys. Res. 112, doi:10.1029/2006JE002869. E06010; Seibert, N.M., Kargel, J.S., 2001. Geophys. Res. Lett. 28, 899-902; Soare, R.J., Kargel, J.S., Osinski, G.R., Costard, F., 2007. Icarus 191, 95-112; Soare, R.J., Osinski, G.R., Roehm, C.L., 2008. Earth Planet. Sci. Lett. 272, 382-393] and glacial processes [Milliken, R.E., Mustard, J.F., Goldsby, D.L., 2003. J. Geophys. Res. 108 (E6), doi:10.1029/2002JE002005. 5057; Mustard, J.F., Cooper, C.D., Rifkin, M.K., 2001. Nature 412, 411-414; Tanaka, K.L., Skinner, J.A., Hare, T.M., 2005. US Geol. Surv. Sci. Invest., Map 2888] have been reported within the region. Researchers have assumed that the periglacial and glacial landscapes occur within the same geological unit, the ABa [i.e., Morgenstern, A., Hauber, E., Reiss, D., van Gasselt, S., Grosse, G., Schirrmeister, L., 2007. J. Geophys. Res. 112; doi:10.1029/2006JE002869. E06010; Tanaka, K.L., Skinner, J.A., Hare, T.M., 2005. US Geol. Surv. Sci. Invest., Map 2888]. In this study we use HiRISE (High Resolution Image Science Experiment, Mars Reconnaissance Orbiter) imagery to identify the stratigraphical separation of the two landscapes and show that periglacial landscape modification has occurred in the geological units that underlie the ABa, not in the ABa itself. Moreover, we suggest that the periglacial landscape extends well beyond the perimeter of the ABa and could be the product of “wet” cold-climate processes. These processes involve freeze-thaw cycles and intermittently stable liquid-water at or near the surface. By contrast, we propose that the ABa is a very recent late-Amazonian geological unit formed principally by “dry” cold-climate processes. These processes comprise accumulation (by atmospheric deposition) and ablation (by sublimation).     

Interferometric millimeter observations of water vapor on Mars and comparison with Mars Express measurements

We present interferometric mapping of the 225.9-GHz HDO and 203.4-GHz lines on Mars obtained with the IRAM Plateau de Bure facility (PdBI). The observations were performed during martian year 28 (MY28), at L_s=320.3° for the HDO line, and at L_s=324.3° for the line. The HDO line is detected at the eastern (morning) and western (evening) limbs of the northern hemisphere, corresponding to a water column density in the range 3-6 pr.-μm. The line is not detected, which is compatible with the column densities derived from the HDO line. Quasi-simultaneous far infrared measurements obtained by the Planetary Fourier Spectrometer (PFS) onboard the Mars Express spacecraft confirm our PdBI results, yielding a 5±1 pr.-μm meridionally constant water column abundance.Such a low water abundance during the southern mid-autumn of MY28 does not correspond to the standard martian climatology as observed during the previous years. It was however already retrieved from near-infrared observations performed by the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) onboard the Mars Reconnaissance Orbiter spacecraft [Smith, M.D., Wolff, M.J., Clancy, R.T., Murchie, S.L. 2009. CRISM observations of water vapor and carbon monoxide. J. Geophys. Res. 114, doi: 10.1029/2008JE003288]. Our observations thus confirm that the planet-encircling dust storm that occurred during MY28 significantly affected the martian water cycle. Our observations also demonstrate the usefulness of interferometric submillimeter observations to survey the martian water cycle from ground-based facilities.     

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.     

Retrieval algorithm for atmospheric dust properties from Mars Global Surveyor Thermal Emission Spectrometer data during global dust storm 2001A

We present a new parameter retrieval algorithm for Mars Global Surveyor Thermal Emission Spectrometer data. The algorithm uses Newtonian first-order sensitivity functions of the infrared spectrum in response to variations in physical parameters to fit a model spectrum to the data at 499, 1099, and 1301 cm⁻¹. The algorithm iteratively fits the model spectrum to data to simultaneously retrieve dust extinction optical depth, effective radius, and surface temperature. There are several sources of uncertainty in the results. The assumed dust vertical distribution can introduce errors in retrieved optical depth of a few tens of percent. The assumed dust optical constants can introduce errors in both optical depth and effective radius, although the systematic nature of these errors will not affect retrieval of trends in these parameters. The algorithm does not include the spectral signature of water ice, and hence data needs to be filtered against this parameter before the algorithm is applied. The algorithm also needs sufficient dust spectral signature, and hence surface-to-atmosphere temperature contrast, to successfully retrieve the parameters. After the application of data filters the algorithm is both relatively accurate and very fast, successfully retrieving parameters, as well as meaningful parameter variability and trends from tens of thousands of individual spectra on a global scale (Elteto, A., Toon, O.B. [2010]. Icarus, this issue). Our results for optical depth compare well with TES archive values when corrected by the single scattering albedo. Our results are on average 1–4 K higher in surface temperatures from the TES archive values, with greater differences at higher optical depths. Our retrieval of dust effective radii compare well with the retrievals of Wolff and Clancy (Wolff, M.J., Clancy, R.T. [2003]. J. Geophys. Res. 108 (E9), 5097) for the corresponding data selections from the same orbits.     

Quantitative compositional analysis of martian mafic regions using the MEx/OMEGA reflectance data 1. Methodology, uncertainties and examples of application

The Mars Express Observatoire pour la Minéralogie, l'Eau, les Glaces et l'Activité (OMEGA) collected an unprecedented visible and near-infrared hyperspectral dataset covering the low albedo regions of Mars. We investigate the ability to infer modal abundance of surfaces of these regions from a radiative transfer model developed by Shkuratov et al. [Shkuratov, Y., Starukhina, L., Hoffmann, H., Arnold, G., 1999. Icarus 137, 235-246] and adapted to basaltic surfaces by Poulet and Erard [Poulet F., Erard, S., 2004. J. Geophys. Res. 109 (E2), doi:10.1029/2003JE002179]. From OMEGA measurements of mafic surfaces, we develop several sensitivity tests to assess the extent to which the model can be applied to predict pyroxene composition (high-calcium phase and low-calcium phase), abundance of almost neutral components (plagioclase) in the near-infrared wavelength as well as grain sizes, by using a library of selected end-members. Results of the sensitivity tests indicate that the scattering model can estimate both abundances and grain sizes of major basaltic materials of low albedo regions within uncertainties (±5 to 15 vol%). The model is then applied to data from dissected cratered terrains located in Terra Meridiani. The derived grain size including uncertainties is in the 50-500 μm range. This is consistent with the thermal inertia and albedo of this region, which indicates a fine sand-sized surface with little dust. The abundances of plagioclase (43-57%) and pyroxenes (35-45±10%, including 11±5% of low-calcium phase) are in good agreement with previous basalt-like compositions of low albedo regions from thermal infrared spectral measurements. The method presented in this paper will provide a valuable tool for evaluating the modal mineralogy of other mafic regions of Mars observed in the near-infrared wavelength range.     

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.     

Valley network-fed, open-basin lakes on Mars: Distribution and implications for Noachian surface and subsurface hydrology

A new catalog of 210 open-basin lakes (lakes with outlet valleys) fed by valley networks shows that they are widely distributed in the Noachian uplands of Mars. In order for an outlet valley to form, water must have ponded in the basin to at least the level of the outlet. We use this relationship and the present topography to directly estimate the minimum amount of water necessary to flood these basins in the past. The volumes derived for the largest lakes (∼3×¹⁰⁴ to ∼2×¹⁰⁵ km³) are comparable to the largest lakes and small seas on modern Earth, such as the Caspian Sea, Black Sea, and Lake Baikal. We determine a variety of other morphometric properties of these lakes and their catchments (lake area, mean depth, volume, shoreline development, outlet elevation, and watershed area). Most candidate lakes have volumes proportional to and commensurate with their watershed area, consistent with precipitation as their primary source. However, other lakes have volumes that are anomalously large relative to their watershed areas, implying that groundwater may have been important in their filling. Candidate groundwater-sourced lakes are generally concentrated in the Arabia Terra region but also include the Eridania basin [Irwin, R.P., Howard, A.D., Maxwell, T.A., 2004a. J. Geophys. Res. 109, doi: 10.1029/2004JE002287. E12009; Irwin, R.P., Watters, T.R., Howard, A.D. Zimbelman, J.R., 2004b. J. Geophys. Res. 109, doi: 10.1029/2004JE002248. E09011] and several lakes near the dichotomy boundary. This areal distribution is broadly consistent with where groundwater should have reached the surface as predicted by current models. Both surface runoff and groundwater flow appear to have been important sources for lakes and lake chains, suggesting a vertically integrated hydrological system, the absence of a global cryosphere, and direct communication between the surface and subsurface hydrosphere of early Mars.     

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.     

Shocked plagioclase signatures in Thermal Emission Spectrometer data of Mars

The extensive impact cratering record on Mars combined with evidence from SNC meteorites suggests that a significant fraction of the surface is composed of materials subjected to variable shock pressures. Pressure-induced structural changes in minerals during high-pressure shock events alter their thermal infrared spectral emission features, particularly for feldspars, in a predictable fashion. To understand the degree to which the distribution and magnitude of shock effects influence martian surface mineralogy, we used standard spectral mineral libraries supplemented by laboratory spectra of experimentally shocked bytownite feldspar [Johnson, J.R., Hörz, F., Christensen, P., Lucey, P.G., 2002b. J. Geophys. Res. 107 (E10), doi:10.1029/2001JE001517] to deconvolve Thermal Emission Spectrometer (TES) data from six relatively large (>50 km) impact craters on Mars. We used both TES orbital data and TES mosaics (emission phase function sequences) to study local and regional areas near the craters, and compared the differences between models using single TES detector data and 3×2 detector-averaged data. Inclusion of shocked feldspar spectra in the deconvolution models consistently improved the rms errors compared to models in which the spectra were not used, and resulted in modeled shocked feldspar abundances of >15% in some regions. However, the magnitudes of model rms error improvements were within the noise equivalent rms errors for the TES instrument [Hamilton V., personal communication]. This suggests that while shocked feldspars may be a component of the regions studied, their presence cannot be conclusively demonstrated in the TES data analyzed here. If the distributions of shocked feldspars suggested by the models are real, the lack of spatial correlation to crater materials may reflect extensive aeolian mixing of martian regolith materials composed of variably shocked impact ejecta from both local and distant sources.     

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.     

Digging deep for ice in Isidis Planitia—New constraints on subsurface volatile concentrations from thermal modeling

The Isidis Planitia region on Mars usually is regarded as a comparably attractive site for landing missions based on engineering constraints such as elevation and smooth regional topography. The Mars Express landed element Beagle 2 was deployed to this area, and the southern margin of the basin was selected as one of the backup landing sites for the NASA Mars Exploration Rovers.Especially in the context of the Beagle 2 mission, Isidis Planitia has been discussed as a place which might have experienced a volatile-rich history with associated potential for biological activity [e.g. Bridges et al., 2003. Selection of the landing site in Isidis Planitia of Mars Probe Beagle 2. J. Geophys. Res. 108(E1), 5001, doi: 10.1029/2001JE001820]. However the measurements of by the GRS instrument on Mars Odyssey indicate a maximum inferred water abundance of only 3 wt% in the upper few meters of the surface [Feldman et al., 2004. Global distribution of near-surface hydrogen on Mars. J. Geophys. Res. 109, E09006, doi: 10.1029/2003JE002160]. Based on these measurements this area seems to be one of the driest spots in the equatorial region of Mars.To support future landing site selections we took a more detailed look at the minimum burial depth of stable ice deposits in this area, focusing as an example on the planned Beagle 2 landing site. We are especially interested in the likelihood of ground ice deposits within the range of proposed subsurface sampling tools as drills or ‘mole’-like devices [Richter et al., 2002. Development and testing of subsurface sampling devices for the Beagle 2 Lander. Planet. Space Sci. 50, 903-913] given reasonable physical constraints for the surface and near surface material.For a mission like ExoMars [Kminek, G., Vago, J.L., 2005. The Aurora Exploration Program—The ExoMars Mission. In: Proceedings of the 35th Lunar and Planetary Science Conference, abstract no. 1111, 15-19 March 2004, League City, TX] with a focus on finding traces of fossil life the area might be of potential interest, because these traces would be better conserved in the dry soil. Modeling and measurement indicate that Isidis Planitia is indeed a dry place and any hypothetical ground ice deposits in this region are out of range of currently proposed sampling devices.     

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.     

Suprathermal electron fluxes on the nightside of Mars: ASPERA-3 observations

Recently aurora-type UV emissions were discovered on the nightside of Mars [Bertaux, J.-L., Leblanc, F., Witasse, O., et al., 2005. Discovery of an aurora on Mars. Nature 439, doi:10.1038/nature03603]. It was suggested that these emissions are produced by suprathermal electrons with energies of tens of eV, rather than by the electrons with spectra peaked above 100 eV [Leblanc, F., Witasse, O., Winningham J., et al., 2006. Origin of the martian aurora observed by spectroscopy for investigation of characteristics of the atmosphere of Mars (SPICAM) onboard Mars Express. J. Geophys. Res. 111, A09313, doi:10.1029/2006JA011763]. In this paper we present observations of fluxes of suprathermal electrons (E_e≈30-100 eV) on the Martian nightside by the ASPERA-3 experiment onboard the Mars Express spacecraft. Narrow spikes of suprathermal electrons are often observed in energy-time spectrograms of electron fluxes at altitudes between 250 and 600 km. These spikes are spatially organized and form narrow strips in regions with strong upward or downward crustal magnetic field. The values of electron fluxes in such events generally could explain the observed auroral UV emissions although a question of their origin (transport from the dayside or local precipitation) remains open.     

Observation of the hydrogen corona with SPICAM on Mars Express

A series of seven dedicated Lyman-α observations of exospheric atomic hydrogen in the martian corona were performed in March 2005 with the ultraviolet spectrometer SPICAM on board Mars Express. Two types of observations are analyzed, observations at high illumination (for a solar zenith angle SZA equal to 30°) and observations at low illumination (for a solar zenith angle equal to 90° (terminator), and near the south pole). The measured Lyman-α emission is interpreted as purely resonant scattering of solar photons. Because the Lyman-α emission is optically thick, we use a forward model to analyze this data. Below the exobase, the hydrogen density is described by a diffusive model and above the exobase, it follows Chamberlain's approach without satellite particles. For different hydrogen density profiles between 80 and 50,000 km, the volume emission rates are computed by solving the radiative transfer equation. Such an approach has been used to analyze the Mariner 6, 7 exospheric Lyman-α data during the late 1960s. A reasonable fit of the set of observations is obtained assuming an exobase temperature between 200 and 250 K and an exobase density of ∼1-4 × 10⁵ cm⁻³ in good agreement with photochemical models. However, when considering the average exospheric temperature of 200 K measured by other methods [Leblanc, F., Chaufray, J.Y., Witasse, O., Lilensten, J., Bertaux, J.-L., 2006a. J. Geophys. Res. 111 (E9), doi:10.1029/2005JE002664. E09S11; Leblanc, F., Chaufray, J.-Y., Bertaux, J.-L., 2007. Geophys. Res. Lett. 34, doi:10.1029/2006GL028437. L02206; Bougher, S.W., Engel, S., Roble, R.G., Foster, B., 2000. J. Geophys. Res. 105, 17669-17692; Mazarico, E., Zuber, M.T., Lemoine, F.G., Smith, D.E., 2007. J. Geophys. Res. 112, doi:10.1029/2006JE002734. E05014] a supplementary hot population is needed above the exobase to reconcile Lyman-α measurements with these previous measurements, particularly for the observations at low SZA. In this case, the densities and temperatures at the exobase for the two populations are 1.0±0.2×¹⁰⁵ cm⁻³ and T=200 K and 1.9±0.5×¹⁰⁴ cm⁻³ and T>500 K for the cold and hot populations respectively at low SZA. At high SZA, the densities and temperatures are equal to 2±0.4×¹⁰⁵ cm⁻³ and T=200 K and n=1.2±0.5×¹⁰⁴ cm⁻³ and T>500 K. These high values of the hot hydrogen component are not presently supported by the theory. Moreover, it is important to underline that the two population model remains relatively poorly constrained by the limited range of altitude covered by the present set of SPICAM measurements and cannot be unambiguously identified because of the global uncertainty of our model and of SPICAM measurements. For a one population solution, an average water escape rate of 7.5 × 10⁻⁴ precipitable μm/yr is estimated, yielding a lifetime of 13,000 years for the average present water vapor content of 10 precipitable microns.     

HiRISE observations of gas sublimation-driven activity in Mars’ southern polar regions: IV. Fluid dynamics models of CO2 jets

The High Resolution Imaging Science Experiment (HiRISE) onboard Mars Reconnaissance Orbiter (MRO) has been used to monitor the seasonal evolution of several regions at high southern latitudes on Mars and, in particular, the jet-like activity which may result from the process described by Kieffer (Kieffer, H.H. [2007]. J. Geophys. Res. (Planets) 112, E08005. doi:10.1029/2006JE002816) involving translucent CO₂ ice. In this work, we concentrate on attempting to model the dusty CO₂ gas jets using a computational fluid dynamics code. Models that included surface slopes of up to 20° (as an analogy to the jet activity seen in “Inca City”, 81°S, 296°E), wind (from 0 to 6 m s⁻¹), variable vent cross-section and length, particles (including a particle size distribution) and mass loading (with dust to gas ratios exceeding 1) were investigated. The structure of the resulting gas jets, the particle distribution within the jets, the deposition patterns (including their dependence on particle size), and the appearance of jets when viewed from different orientations (including from a nadir-pointing camera) have been investigated for a range of input parameters. The results provide predictions for the size-dependency of altitudes of particles within a plume and the distribution of particle sizes in the deposition fans. Where slopes are a dominant influence, larger particles are expected to be seen furthest from the vent. Where wind is dominant, smaller particles should travel to larger distances. Models producing deposition patterns consistent in length (∼80 m) and form with fans observed by HiRISE on MRO have been demonstrated. The models also suggest that downward flow of gas produced by drag effects from particles falling from the jet under gravity could provide a mechanism for the production of bright haloes which are observed to surround dark fan deposits in MOC, HiRISE and CRISM.