Ionospheric corrections of MARSIS subsurface sounding signals with filters including collision frequency

Ionospheric corrections of Mars advanced radar for subsurface and ionosphere sounding (MARSIS) subsurface sounding signals are necessary before they can be further analyzed. Usually the ionosphere correction only considers the phase dispersion owing to the electron densities. In this paper we show that if the electron-neutral collision frequency is included in the correction filter, the signal-to-noise ratio of the processed signal can be further maximized and the spatial resolution of the signal be improved. Three different models of the ionosphere profile have been studied, and it is shown that a uniform slab model, of both densities and collisions, is feasible and probably the best choice when implementing the correction filter including collisions. The physical significance of the parameters obtained with the uniform model (equivalent parameters) is discussed, and it is shown that they are useful for the estimation of the real physical parameters of the ionosphere. We developed a recursive, random search algorithm to facilitate the realization of the correction process with the filter including collision frequency.     

MARSIS surface reflectivity of the south residual cap of Mars

The south residual cap of Mars is commonly described as a thin and bright layer of CO₂-ice. The Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) is a low-frequency radar on board Mars Express operating at the wavelength between 55 and 230 m in vacuum. The reflection of the radar wave on a stratified medium like the residual cap can generate interferences, causing weaker surface reflections compared to reflections from a pure water ice surface. In order to understand this anomalous low reflectivity, we propose a stratified medium model, which allows us to estimate both the thickness and the dielectric constant of the optically thin slab. First, we consider the residual cap as single unit and show that the decrease in the reflected echo strength is well explained by a mean thickness of 11 m and a mean dielectric constant of 2.2. This value of dielectric constant is close to the experimental value 2.12 for pure CO₂-ice. Second, we study the spatial variability of the radar surface reflectivity. We observe that the reflectivity is not homogeneous over the residual cap. This heterogeneity can be modeled either by variable thickness or variable dielectric constant. The surface reflectivity shows that two different units comprise the residual cap, one central unit with high reflectivity and surrounding, less reflective units.     

Diagnostic of the dayside ionosphere of Mars using the Total Electron Content measurement by the NEIGE/Netlander experiment: An assessment study

The NEtlander Ionosphere and Geodesy Experiment (NEIGE) of the Netlander mission to Mars will measure Doppler shifts affecting the radio links between ground stations and an orbiter. The experiment has two complementary scientific objectives which are the monitoring of the structure and dynamics of the ionosphere of Mars and the precise determination of Mars orientation parameters. The horizontal variation of the Total Electron Content (TEC) of the ionosphere will be derived from the so-called “geometric-free” combination of the Doppler shifts which affect radio links at two frequencies (in the UHF and S bands) between the Netlander microstations on the Mars surface and the data-relay orbiter. We describe a new method for retrieving the horizontal profile of the absolute value of the TEC. Simulations have allowed to evaluate the precision in the determination of the TEC using this method. We show that the daytime TEC can be retrieved with a precision of for a nominal accuracy of on the orbital pseudo-velocity, which represents a relative precision of a few percent. A preliminary analysis of the sensitivity of the TEC to the physical parameters which control the ionosphere has been performed. For this purpose, we have used a new one-dimensional ionospheric model based on the solution of coupled kinetic, fluid and MHD equations. This model describes the suprathermal electron component, the thermal plasma component as well as the induced horizontal magnetic field. The code which provides the vertical electron density profile has been used to study the variation of the TEC with the solar zenith angle and with the induced magnetic field at the top of the ionosphere. In particular, we show that NEIGE will allow to diagnose the penetration into the daytime ionosphere of an induced magnetic field.     

Top layers charaterization of the Martian surface: Permittivity estimation based on geomorphology analysis

Mars Express spacecraft inserted successfully Martian orbit at the end of 2003. On board this probe, a radar instrument called MARSIS (for Mars Advanced Radar for Surface and Ionosphere Sounding) is looking for water inside the first kilometers of Martian crust. To support MARSIS planning and data inversion, Laboratoire de Planétologie de Grenoble developed a MARSIS signal simulator.We show in this paper that MARSIS can also characterize some surface features, in addition to subsurface water and ionosphere sounding. We study a Martian surface region of special interest: Nilokeras Mensae, inside Acidalia Planitia. We discuss the previous geological studies of this region, and show the geomorphologies analyze of this surface area could lead to a simple terrain model. Then, we present a possible data inversion scheme and applying the MARSIS simulator, we test a first radar data inversion.Finally, we will show that complete dielectric characteristics of surface top layers can be retrieved, at least as often Mars Express flies over some layered terrain (at wavelength scale).     

MARSIS radar sounder observations in the vicinity of Ma'adim Vallis, Mars

The MARSIS radar experiment aboard the ESA Mars Express satellite has recorded several unusual reflections in the Ma'adim Vallis region of Mars. These reflections display a wide variety of morphologies which are very different from those of reflections seen beneath the Polar Layered Deposits, Medusae Fossae Formation and Dorsa Argentea Formation. Their morphologies are sometimes very laterally extensive, parabolic or hyperbolic, and apparently deep, but they can also appear horizontal and shallow. Aided by a geological map of the Ma'adim Vallis region, the morphological, locational and temporal characteristics of the reflections have been studied individually in an attempt to constrain their origin. While some may be subsurface reflections based on their shallow morphologies and correlation with the Eridania Planitia basin network, all of the reflections are ambiguous to some degree, displaying characteristics that do not allow a definite subsurface- or possibly ionospheric-sourced mechanism to be proposed for their creation. Those with more exaggerated morphologies are regarded as being much more likely to result from ionospheric distortion rather than subsurface inhomogeneity.     

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.     

Dayside induced magnetic field in the ionosphere of Mars

The Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) onboard the Mars Express spacecraft has occasionally displayed surprising features. One such feature is the occurrence of a series of broadband, low-frequency echoes at equally spaced delay times after the sounder transmitter pulse. The interval between the echoes has been shown to be at the cyclotron period of electrons orbiting in the local magnetic field. The electrons are believed to be accelerated by the large voltages applied to the antenna by the sounder transmitter. Measurements of the period of these “electron cyclotron echoes” provide a simple technique for determining the magnitude of the magnetic field near the spacecraft. These measurements are particularly useful because Mars Express carries no magnetometer, so this is the only method available for measuring the magnetic field magnitude. Using this technique, results are presented showing the large scale structure of the draped field inside the magnetic pile-up boundary. The magnitude of the draped field is shown to vary from about 40 nT at a solar zenith angle of about 25°, to about 25 nT at a solar zenith angle of 90°. The results compare favorably with similar results from the Mars Global Surveyor spacecraft. A fitting technique is developed to derive the vector direction and magnitude of the draped magnetic field in cases where the spacecraft passes through regions with significant variation in the crustal field. The magnetic field directions are consistent with current knowledge of the draping geometry of the magnetic field around Mars.     

Evidence for topographic effects on the Martian ionosphere

The Radio occultation experiment on board Mariner 9 has been used to demonstrate that the altitude of the main electron density peak in the Martian ionosphere is closely related to the height of Mars’ surface at the occultation point. This is direct evidence for topographic effects on the Martian ionosphere. Modeling indicates that topographic-induced modulations of the neutral density in the upper atmosphere can account for the observed ionospheric effects. The neutral density modulation is likely to be caused by nonmigrating tides in the Martian thermosphere.     

Shallow radar (SHARAD) sounding observations of the Medusae Fossae Formation, Mars

The SHARAD (shallow radar) sounding radar on the Mars Reconnaissance Orbiter detects subsurface reflections in the eastern and western parts of the Medusae Fossae Formation (MFF). The radar waves penetrate up to 580 m of the MFF and detect clear subsurface interfaces in two locations: west MFF between 150 and 155° E and east MFF between 209 and 213° E. Analysis of SHARAD radargrams suggests that the real part of the permittivity is ∼3.0, which falls within the range of permittivity values inferred from MARSIS data for thicker parts of the MFF. The SHARAD data cannot uniquely determine the composition of the MFF material, but the low permittivity implies that the upper few hundred meters of the MFF material has a high porosity. One possibility is that the MFF is comprised of low-density welded or interlocked pyroclastic deposits that are capable of sustaining the steep-sided yardangs and ridges seen in imagery. The SHARAD surface echo power across the MFF is low relative to typical martian plains, and completely disappears in parts of the east MFF that correspond to the radar-dark Stealth region. These areas are extremely rough at centimeter to meter scales, and the lack of echo power is most likely due to a combination of surface roughness and a low near-surface permittivity that reduces the echo strength from any locally flat regions. There is also no radar evidence for internal layering in any of the SHARAD data for the MFF, despite the fact that tens-of-meters scale layering is apparent in infrared and visible wavelength images of nearby areas. These interfaces may not be detected in SHARAD data if their permittivity contrasts are low, or if the layers are discontinuous. The lack of closely spaced internal radar reflectors suggests that the MFF is not an equatorial analog to the current martian polar deposits, which show clear evidence of multiple internal layers in SHARAD data.     

On the sublimation of ice particles on the surface of Mars; with applications to the 2007/8 Phoenix Scout mission

Experimental studies related to the sublimation of ice, in bulk or as small particles, alone or mixed with dust similar to that expected on the surface of Mars, are reported. The experiments, a cloud physics particle sublimation model, and a convection model presented by Ingersoll, all indicate a strong dependence of sublimation rate on temperature, and this appears to be the dominant factor, assuming that the relative humidity of the air is fairly low. In addition the rate of loss of water vapour appears to depend primarily on exposed surface area and less on particle size and the total mass of the sample, or the mass of ice in the sample. The 2007/8 Phoenix Scout mission plans to obtain and analyse samples of sub-surface ice from about 70° N on Mars. A concern is that these samples, in the form of ice chips of size about 1 mm diameter, could be prone to sublimation when exposed for prolonged periods (many hours) to a relatively warm and dry atmosphere. Our laboratory simulations confirm that this could be a problem if particles are simply left lying on the surface, but also indicate that samples kept suitably cold and collected together in confined piles will survive long enough for the collection and delivery (to the analysis instruments) procedure to be completed.     

Apparent thermal inertia and the surface heterogeneity of Mars

Thermal inertia derivation techniques generally assume that surface properties are uniform at horizontal scales below the footprint of the observing instrument and to depths of several decimeters. Consequently, surfaces with horizontal or vertical heterogeneity may yield apparent thermal inertia which varies with time of day and season. To investigate these temporal variations, we processed three Mars years of Mars Global Surveyor Thermal Emission Spectrometer observations and produced global nightside and dayside seasonal maps of apparent thermal inertia. These maps show broad regions with diurnal and seasonal differences up to 200 J m⁻² K⁻¹s^−1/2 at mid-latitudes (60° S to 60° N) and 600 J m⁻² K⁻¹s^−1/2 or greater in the polar regions. We compared the seasonal mapping results with modeled apparent thermal inertia and created new maps of surface heterogeneity at 5° resolution, delineating regions that have thermal characteristics consistent with horizontal mixtures or layers of two materials. The thermal behavior of most regions on Mars appears to be dominated by layering, with upper layers of higher thermal inertia (e.g., duricrusts or desert pavements over fines) prevailing in mid-latitudes and upper layers of lower thermal inertia (e.g., dust-covered rock, soils with an ice table at shallow depths) prevailing in polar regions. Less common are regions dominated by horizontal mixtures, such as those containing differing proportions of rocks, sand, dust, and duricrust or surfaces with divergent local slopes. Other regions show thermal behavior that is more complex and not well-represented by two-component surface models. These results have important implications for Mars surface geology, climate modeling, landing-site selection, and other endeavors that employ thermal inertia as a tool for characterizing surface properties.     

Stability of hydrous minerals on the martian surface

The presence of water-bearing minerals on Mars has long been discussed, but little or no data exist showing that minerals such as smectites and zeolites may be present on the surface in a hydrated state (i.e., that they could contain H₂O molecules in their interlayer or extra-framework sites, respectively). We have analyzed experimental thermodynamic and X-ray powder diffraction data for smectite and the most common terrestrial zeolite, clinoptilolite, to evaluate the state of hydration of these minerals under martian surface conditions. Thermodynamic data for clinoptilolite show that water molecules in its extra-framework sites are held very strongly, with enthalpies of dehydration for Ca-clinoptilolite up to three times greater than that for liquid water. Using these data, we calculated the Gibbs free energy of hydration of clinoptilolite and smectite as a function of temperature and pressure. The calculations demonstrate that these minerals would indeed be hydrated under the very low-P (H₂O) conditions existing on Mars, a reflection of their high affinities for H₂O. These calculations assuming the partial pressure of H₂O and the temperature range expected on Mars suggest that, if present on the surface, zeolites and Ca-smectites could also play a role in affecting the diurnal variations in martian atmospheric H₂O because their calculated water contents vary considerably over daily martian temperature ranges. The open crystal structure of clinoptilolite and existing hydration and kinetic data suggest that hydration/dehydration are not kinetically limited. Based on these calculations, it is possible that hydrated zeolites and clay minerals may explain some of the recent observations of significant amounts of hydrogen not attributable to water ice at martian mid-latitudes.     

MARCI and MOC observations of the atmosphere and surface cap in the north polar region of Mars

We used MGS-MOC and MRO-MARCI daily mapping images of the North Polar Region of Mars from 16 August 2005 (L_s = 270°) to 21 May 2009 (L_s = 270°), covering portions of three consecutive martian years (MY 27-MY 29), to observe the seasonal behavior of the polar ice cap and atmospheric phenomena. The rate of cap regression was similar in MY 28 and MY 29, but was advanced by 3.5° of L_s (∼7-8 sols) in MY 29. The spatial and temporal behaviors of dust and condensate clouds were similar in the two years and generally in accord with prior years. Dust storms (>100 km²) were observed in all seasons, with peak activity occurring at L_s = 10-20° from 50°N to 70°N and at L_s = 135-140° from 70°N to 90°N. The most active quadrant was 0-90°W in MY 28, shifting to 180-270°W in MY 29. The majority of regional storms in both years developed in longitudes from 10°W to 60°W. During late summer the larger storms obscure the North Polar Region in a cloud of dust that transitions to north polar hood condensate clouds around autumnal equinox.Changes in the distribution of perennial ice deposits, especially in Olympia Planum, were observed between the 2 years, with the MY 29 ice distribution being the most extensive observed to date. Modeling suggests that the small, bright ice patches on the residual cap are not the result of slope or elevation effects. Rather we suggest that they are the result of local meteorological effects on ice deposition. The annual darkening and brightening of peripheral areas of the residual cap around summer solstice can be explained by the sublimation of a brighter frost layer revealing an underlying darker, ice rich layer that itself either sublimes to reveal brighter material below or acts as a cold trap, attracting condensation of water vapor that brightens the surface. An alternative explanation invokes transport and deposition of dust on the surface from the cap interior, and later removal of that dust. The decrease in cap albedo and accompanying increase in near surface atmospheric stability may be related to the annual minimum of polar storm activity near northern summer solstice.     

Aram Chaos and its constraints on the surface heat flux of Mars

The surface heat flux of a planet is an important parameter to characterize its internal activity and to determine its thermal evolution. Here we report on a new method to constrain the surface heat flux of Mars during the Hesperian. For this, we explore the consequences for the martian surface heat flux from a recently presented new hypothesis for the formation of Aram Chaos (Zegers, T.E., Oosthoek, J.H.P., Rossi, A.P., Blom, J.K., Schumacher, S. [2010]. Earth Planet. Sci. Lett. 297, 496-504. doi:10.1016/j.epsl.2010.06.049.). In this hypothesis the chaotic terrain is thought to have formed by melting of a buried ice sheet. The slow sedimentation and burial of the ice sheet led to an increased thermal insulation of the ice and subsequently to a temperature increase high enough to trigger melting and the formation of the subsurface lake. As these processes highly depend on the thermal properties of the subsurface and especially on the surface heat flux, it is possible to constrain the latter by using numerical simulations. Based on the hypothesis for the formation of Aram Chaos, we conducted an extensive parameter study to determine the parameter settings leading to sufficient melting of the buried ice sheet. We find that the surface heat flux in the Aram Chaos region during the Hesperian was most likely between 20 and 45 mW m⁻² with a possible maximum value of up to 60 mW m⁻².     

Energetic Neutral Atoms (ENA) at Mars: Properties of the hydrogen atoms produced upstream of the martian bow shock and implications for ENA sounding technique around non-magnetized planets

We have studied the interaction of fast solar wind hydrogen atoms with the martian atmosphere by a three-dimensional Monte Carlo simulation. These energetic neutral hydrogen atoms, H-ENAs, are formed upstream of the martian bow shock. Both H-ENAs scattered and non-scattered from the martian atmosphere/exosphere were studied. The colliding H-ENAs were found to scatter both to the dayside and nightside. On the dayside they contribute to the so-called H-ENA albedo. On the nightside the heated and scattered hydrogen atoms were found also in the martian wake. The density, the energy distribution function and the direction of the velocity of H-ENAs on the nightside are presented. The present study describes a novel “ENA sounding” technique in which energetic neutral atoms are used to derive information of the properties of planetary exosphere and atmosphere in a similar manner as the solar wind photons are used to derive atmospheric densities by measuring the scattered UV light. A detailed study of the direction and energy of the scattered and non-scattered H-ENAs suggest that the ENA sounding is a method to study the interaction between the planetary atmosphere and the solar wind and to monitor the density, and likely also the magnetization, of the planetary upper atmosphere. Already present-day ENA instrument should be capable to detect the analyzed particle fluxes.     

Sensitivity of a Titan ionospheric model to the ion-molecule reaction parameters

Various aspects of ion-molecule reactions for Titan ionospheric chemistry modeling are reviewed in this work: temperature/collision energy effects on reaction rates and, more importantly, on products distributions; differential reactivity of isomers of ions; reactivity of excited states of ions; pathways to the building of complex ions. We evaluate here the extent to which these points affect model predictions. We find that the present limiting factors to model predictivity are the model incompleteness for heavy ion production pathways; the differential reactivity of isomers; and, to a lesser degree, the temperature effects on the branching ratios of ion-molecule reactions. Extensive experimental studies are required to fill these gaps in ion-molecule reactivity knowledge.