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.     

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.     

The Beagle 2 environmental sensors: science goals and instrument description

A suite of instruments on the Beagle 2 Mars lander was designed and built in order to investigate the environmental conditions at the landing site. The sensor suite was capable of measuring air temperature at two heights, surface level pressure, wind speed and direction, saltated particle momentum, UV flux (diffuse and direct at five wavelengths), the total accumulated radiation dose and investigating the nature of the oxidising environment. The scientific goals of the instruments are discussed within the context of current understanding of the environmental conditions on Mars, and the instruments themselves are described in detail. Beagle 2 landed on Mars in late 2003, as part of the ESA Mars Express mission. The expected lifetime of the lander on the surface was 180 sols, with a landing site in Isidis Planitia, but has not responded to attempts to contact it, and has now been declared lost. The Environmental Sensor Suite (ESS) was intended to monitor and characterise the current local meteorological parameters, investigating specific areas of scientific interest raised from previous missions, most notably dust transport and transient phenomena, and additionally to add context to the conditions that any possible martian micro-organisms would have to face. The design of the instrument suite was strongly influenced by mass limitations, with eight sensor subsystems having a total mass of approximately 100 g. Although Beagle 2 has been now declared lost, the scientific goals of an Environmental Sensors Suite still remain a valid target for any future astrobiology orientated missions.     

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.     

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.     

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.     

Mars Exploration Rover candidate landing sites as viewed by THEMIS

The analysis of six landing sites that were candidates for the two NASA Mars Exploration Rovers (MER) benefited from recently available image data from the Thermal Emission Imaging Spectrometer (THEMIS) onboard the 2001 Mars Odyssey spacecraft. The combination of daytime and nighttime thermal infrared images from THEMIS supplemented by additional data sets has lead to new or expanded insights into the nature of each landing site. In Meridiani Planum, a layer of lighter-toned, higher thermal inertia material is observable just below the hematite-bearing layer. Gusev Crater displays a more complex stratigraphy than previously observed, including an upper layer with lobate margins. The highest inertia unit of southern Isidis Planitia is confined to topographic lows in the rim/basin margin and does not appear to be due to highland material transported onto the basin floor. The enigmatic, ovoid, blocky terrain on the floor of Melas Chasma displays higher thermal inertia than its surroundings, an indication that it contains coarser or more indurated material than the adjacent aeolian bedforms. The myriad channels of Athabasca Valles display distinctive thermal signatures despite the presence of a bright layer of dust covering the region. The presence of alluvial fans produced from spur-and-gulley erosion of the walls of Eos Chasma demonstrates that mass movements have occurred following the canyon scouring floods.     

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

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.     

Lithospheric structure and tectonics at Isidis Planitia, Mars

We characterize the lithospheric structure of Isidis Planitia on Mars by analyzing Mars Global Surveyor and Mars Odyssey gravity and topography data using a flexural model of a thin elastic shell including bending and membrane stresses. Isidis Planitia is a circular, relatively flat plain formed near the end of the Early Noachian, at the edge of the highlands-lowlands boundary and the site of a large free-air gravity anomaly, features consistent with modification and filling of an impact basin. Our results suggest that the bulk density of the fill material inside Isidis must be more than 2600 kg m⁻³ and higher densities are probable. A comparison of the faulting observed at Nili Fossae to the predicted zone of extensional strain northwest of Isidis constrains the thickness of the elastic lithosphere to be 100-180 km thick beneath the basin at the time of loading. We also find that loads outside of the basin play a significant role in the interpretation of the tectonics; simplified models tend to overestimate the lithospheric thickness. We place relatively narrow bounds on the thermal gradient (3.4-6.5 K km⁻¹) and heat flux (13.6-26 mW m⁻²) at Isidis at the time of loading.     

Major episodes of geologic history of Isidis Planitia on Mars

We have mapped the area of Isidis Planitia (1–27°N, 75–103°E) in order to assess the geologic history of this region using modern data sets such as MOLA topography and the high-resolution images provided by the HRSC, CTX, and HiRISE cameras. Results of our mapping show that the geologic history of Isidis Planitia consists of three principal episodes. (1) Impact dominated episode (Noachian, until ～3.8 Ga): During this time, the oldest materials in the study area were formed mostly by impact reworking and mass-wasting. Other processes (e.g., volcanism and fluvial/glacial activity) likely operated at this time but played a subordinate role. (2) An episode related to volcanic and fluvial/glacial activities (late Noachian–early Amazonian, ～3.8–2.8 Ga): Volcanism appears as the most important process at the beginning of this episode (～3.8–3.5 Ga) and was responsible for the formation of a large circum-Isidis volcanic province by the early Hesperian epoch. Volcanic materials covered large portions of the Isidis rim, almost completely buried the previous crater record on the floor of the Isidis basin, and probably were the major contributors to the filling of the basin. Fluvial/glacial processes prevailed closer to the end of the episode (early Hesperian–early Amazonian, ～3.5–2.8 Ga) and were responsible for widespread resurfacing in the Isidis Planitia region, mostly at ～3.1–3.4 Ga. Glaciers and/or ice sheets probably resulted in a massive glaciation of the rim and the floor of the Isidis basin. The total volume of material eroded from the Isidis rim by glacial and fluvial activity is estimated to be about 35,000–50,000 km³, which is equivalent to a composite layer about 40–60 m thick on entire floor of the basin. More important, however, is that the eroded materials were likely saturated with ice/water and could form wet deposits on the floor. (3) Wind-dominated episode (since early Amazonian, ～2.8 Ga): Wind activity dominated the later geologic history of Isidis Planitia but resulted only in minor modification of the surface.     

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.     

Boulder Analysis on the Oxia Planum ExoMars 2022 Rover Landing Site: Scientific and Engineering Perspectives

Solar System Research - The Oxia Planum region has been chosen as the landing site for the future ESA ExoMars 2022 rover for both scientific value and engineering safety (Ivanov et al., 2020). The...     

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.     

Weathering of Fe-bearing minerals under Martian conditions, investigated by Mössbauer spectroscopy

The surface of Mars is covered by weathered material. Mars' rusty red colour in particular is commonly ascribed to ferric iron-bearing minerals. The planet's surface is generally iron rich. Mössbauer spectroscopy is a powerful tool for quantitative mineralogical analysis of Fe-bearing minerals. Consequently, the miniaturized Mössbauer spectrometer MIMOS II is part of the payload of NASA's twin Mars Exploration Rovers “Spirit” and “Opportunity”, and ESA's ill-fated Mars Express lander “Beagle 2”. Both Mars Exploration Rovers are currently conducting successful surface operations on Mars. In this paper, we give a brief insight into mission operations with respect to the reconstruction of local weathering scenarios at the landing sites, which in turn will help to illuminate the climatic history of the planet. Mössbauer spectra obtained in preparation of the mission from the SNC meteorites Nakhla, Dar al Gani 476, and Sayh al Uhaymir, show weathering and other alteration features. Preliminary results of laboratory weathering experiments on Fe-bearing minerals (olivine and pyroxene) show the importance of analysing individual minerals to understand weathering of more complex mineral assemblages like, e.g., basalt.     

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.     

Beagle-2 position determination from the returned camera panoramas using MOLA data

After the Beagle-2 lander of the Mars Express mission comes to rest on the surface of Isidis Planitia in late December 2003 to carry out a range of geochemistry and exobiology experiments, there will be considerable interest in determining its exact location. This work considers the feasibility of identifying topographic features seen in the Mars Global Surveyor MOLA dataset in images of the horizon returned by the lander, and the probability of observing lesser features identifiable in orbital imagery. By taking bearings from such features, and attempting to match the configuration back to the available data, it may be possible to determine the spacecraft's position with high precision. Since the MOLA data is fairly coarse compared to the area of the landing ellipse, the range of visibility and likelihood of observation of each of the resolved craters in the area is considered. For the more numerous smaller craters and many small knobs a probabilistic view is taken.