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.     

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.     

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.     

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.     

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.     

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.     

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.     

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).     

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.     

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.     

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.     

Interfacial liquid water on Mars and its potential role in formation of hill and dune gullies

Gullies are among the most intriguing structures identified on the surface of Mars. Most common are gullies located on the slopes of craters which are probably formed by liquid water transported by shallow aquifers (Heldmann, J.L., Carlsson, E., Johansson, H., Mellon, M.T., Toon, O.B. [2007]. Icarus 188, 324-344). Two particular types of gullies are found on slopes of isolated hills and dunes. The hill-slope gullies are located mostly at 50°S, which is at the high end of latitudes of bulk of the gullies found so far. The dune gullies are found in several locations up to 65°S (Reiss, D., Jaumann, R., Kereszturi, A., Sik, A., Neukum, G. [2007]. Lunar Planet. Sci. XXXVIII. Abstract 1993), but the best known are those in Russel crater at 54°S. The hill and dune gullies are longer than others making the aquifers explanation for their formation unlikely (Balme, M., Mangold, N., Baratoux, D., Costard, F., Gosselin, M., Masson, P., Pnet, P., Neukum, G. [2006]. J. Geophys. Res. 111. doi:10.1029/2005JE002607). Recently it has been noted that thin liquid films of interfacial water can play a role in rheological processes on the surface of Mars (Moehlmann, D. [2008]. Icarus 195, 131-139. Kereszturi, A., Moehlmann, D., Berczi, Sz., Ganti, T., Kuti, A., Sik, A., Horvath, A. [2009]. Icarus 201, 492-503.). Here we try to answer the question whether interfacial liquid water may occur on Mars in quantities large enough to play a role in formation of gullies. To verify this hypothesis we have calculated thermal models for hills and dunes of various steepness, orientation and physical properties. We find that within a range of average expected values of parameters it is not possible to have more than a few monolayers of liquid water at depths greater than a centimeter. To create subsurface interfacial water film significantly thicker and hence to produce conditions for the slope instability, parameters have to be chosen to have their extreme realistic values or an additional source of surface heating is needed. One possibility for additional heating is a change of atmospheric conditions due to a local dust storm. We conclude that if interfacial water is responsible for the formation of the hill-slope gullies, our results may explain why the hill gullies are rare.     

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.     

Could giant impacts cripple core dynamos of small terrestrial planets?

Large impacts not only create giant basins on terrestrial planets but also heat their interior by shock waves. We investigate the impacts that have created the largest basins existing on the planets: Utopia on Mars, Caloris on Mercury, Aitken on Moon, all formed at ∼4 Ga. We determine the impact-induced temperature increases in the interior of a planet using the “foundering” shock heating model of Watters et al. (Watters, W.A., Zuber, M.T., Hager, B.H. [2009]. J. Geophys. Res. 114, E02001. doi:10.1029/2007JE002964). The post-impact thermal evolution of the planet is investigated using 2D axi-symmetric convection in a spherical shell of temperature-dependent viscosity and thermal conductivity, and pressure-dependent thermal expansion. The impact heating creates a superheated giant plume in the upper mantle which ascends rapidly and develops a strong convection in the mantle of the sub-impact hemisphere. The upwelling of the plume rapidly sweeps up the impact-heated base of the mantle away from the core-mantle boundary and replaces it with the colder surrounding material, thus reducing the effects of the impact-heated base of the mantle on the heat flux out of core. However, direct shock heating of the core stratifies the core, suppresses the pre-existing thermal convection, and cripples a pre-existing thermally-driven core dynamo. It takes about 17, 4, and 5 Myr for the stratified cores of Mars, Mercury, and Moon to exhaust impact heat and resume global convection, possibly regenerating core dynamos.     

Variability of CO concentrations in the Venus troposphere from Venus Express/VIRTIS using a Band Ratio Technique

A fast method is presented for deriving the tropospheric CO concentrations in the Venus atmosphere from near-infrared spectra using the night side 2.3 μm window. This is validated using the spectral fitting techniques of Tsang et al. [Tsang, C.C.C., Irwin, P.G.J., Taylor, F.W., Wilson, C.F., Drossart, P., Piccioni, G., de Kok, R., Lee, C., Calcutt, S.B., and the Venus Express/VIRTIS Team, 2008a. Tropospheric carbon monoxide concentrations and variability on Venus with Venus Express/VIRTIS-M observations. J. Geophys. Res. 113, doi: 10.1029/2008JE003089. E00B08] to show that monitoring CO in the deep atmosphere can be done quickly using large numbers of observations, with minimal effect from cloud and temperature variations. The new method is applied to produce some 1450 zonal mean CO profiles using data from the first eighteen months of operation from the Visible and Infrared Thermal Imaging Spectrometer infrared mapping subsystem (VIRTIS-M-IR) on Venus Express. These results show many significant long- and short-term variations from the mean equator-to-pole increasing trend previously found from earlier Earth- and space-based observations, including a possible North-South dichotomy, with interesting implications for the dynamics and chemistry of the lower atmosphere of Venus.     

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.     

Stratigraphical and morphological evidence for pingo genesis in the Cerberus plains

‘Rootless’ debris cones (or pseudocraters) occur in platy, patterned ground throughout the Cerberus plains of Mars and are thought to represent the products of explosive magma-ice interaction [Lanagan et al., 2001. Geophys. Res. Lett. 28, 2365-2368; Fagents et al., 2002. In: Smellie, J.L., Chapman, M.G. (Eds.), Volcano-Ice Interaction on Earth and Mars. In: Geol. Soc. Spec. Publ., vol. 202, pp. 295-317]. Requiring lava and water interspersed, they are central to theories of multiple magmatic and aqueous flood events [Burr et al., 2002. Icarus 159, 53-73; Berman, D.C., Hartmann, W.K., 2002. Icarus 159, 1-17] and widespread sheet volcanism [Keszthelyi et al., 2000. J. Geophys. Res. 105, 15027-15049] in the region during the late Amazonian (a region reported to have been occupied by water bodies ranging from lakes to oceans [Scott et al., 1995. Map of Mars showing channels and possible paleolake basins. USGS Miscellaneous Investigations Series, Map I-2461 (1:30,000,000)]). The nature of the platy substrate is the subject of debate, with evidence given for lava [Keszthelyi et al., 2000. J. Geophys. Res. 105, 15027-15049; Plescia, J.B., 2003. Icarus 164, 79-95] and ice [Brakenridge, G.R., 1993. Lunar Planet. Sci. XXIV (Part 1), 175-176; Rice et al., 2002. Lunar Planet. Sci. XXXIII. Abstract #2026; Murray et al., 2005. Nature 434, 352-355]. The superposition relationships of cones and platy deposits in the channels of the Athabasca Valles precludes a magmatic origin, indicating later formation as permafrost mounds (or ‘pingos’), with implications for geologically recent flood volcanism, age constraints on young surfaces and recent climate change on Mars.     

Numerical modeling of protocore destabilization during planetary accretion: Methodology and results

We developed and tested an efficient 2D numerical methodology for modeling gravitational redistribution processes in a quasi spherical planetary body based on a simple Cartesian grid. This methodology allows one to implement large viscosity contrasts and to handle properly a free surface and self-gravitation. With this novel method we investigated in a simplified way the evolution of gravitationally unstable global three-layer structures in the interiors of large metal-silicate planetary bodies like those suggested by previous models of cold accretion [Sasaki, S., Nakazawa, K., 1986. J. Geophys. Res. 91, 9231-9238; Karato, S., Murthy, V.R., 1997. Phys. Earth Planet Interios 100, 61-79; Senshu, H., Kuramoto, K., Matsui, T., 2002. J. Geophys. Res. 107 (E12), 5118. 10.1029/2001JE001819]: an innermost solid protocore (either undifferentiated or partly differentiated), an intermediate metal-rich layer (either continuous or disrupted), and an outermost silicate-rich layer. Long-wavelength (degree-one) instability of this three-layer structure may strongly contribute to core formation dynamics by triggering planetary-scale gravitational redistribution processes. We studied possible geometrical modes of the resulting planetary reshaping using scaled 2D numerical experiments for self-gravitating planetary bodies with Mercury-, Mars- and Earth-size. In our simplified model the viscosity of each material remains constant during the experiment and rheological effects of gravitational energy dissipation are not taken into account. However, in contrast to a previously conducted numerical study [Honda, R., Mizutani, H., Yamamoto, T., 1993. J. Geophys. Res. 98, 2075-2089] we explored a freely deformable planetary surface and a broad range of viscosity ratios between the metallic layer and the protocore (0.001-1000) as well as between the silicate layer and the protocore (0.001-1000). An important new prediction from our study is that realistic modes of planetary reshaping characterized by a high viscosity protocore and low viscosity molten silicate and metal [Senshu, H., Kuramoto, K., Matsui, T., 2002. J. Geophys. Res. 107 (E12), 5118. 10.1029/2001JE001819] may result in the transient exposure of the protocore to the planetary surface and a strongly (up to 8% of the planetary diameter) aspherical deviation of the planetary shape during the early stages of core formation. Exposure of the protocore might happen in the early stages of iron core formation. This process may conceivably convert a large amount of potential energy into temperature increase and a transient strongly non-uniform depth of the magma ocean around the protoplanet. Our simplified model also predicts that the time for metallic core formation out of the metal-rich layer depends mainly on the dynamics of the deformation of the solid strong protocore. In nature this dynamics will be strongly dependent on the effective viscosity of the protocore, which should generally have non-Newtonian pressure-, temperature-, and stress-dependent rheology with strong thermomechanical feedbacks from gravitational energy dissipation.     

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].     

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.