Lander radioscience for obtaining the rotation and orientation of Mars

The paper presents the concept, the objectives, the approach used, and the expected performances and accuracies of a radioscience experiment based on a radio link between the Earth and the surface of Mars. This experiment involves radioscience equipment installed on a lander at the surface of Mars. The experiment with the generic name lander radioscience (LaRa) consists of an X-band transponder that has been designed to obtain, over at least one Martian year, two-way Doppler measurements from the radio link between the ExoMars lander and the Earth (ExoMars is an ESA mission to Mars due to launch in 2013). These Doppler measurements will be used to obtain Mars’ orientation in space and rotation (precession and nutations, and length-of-day variations). More specifically, the relative position of the lander on the surface of Mars with respect to the Earth ground stations allows reconstructing Mars’ time varying orientation and rotation in space.Precession will be determined with an accuracy better by a factor of 4 (better than the 0.1% level) with respect to the present-day accuracy after only a few months at the Martian surface. This precession determination will, in turn, improve the determination of the moment of inertia of the whole planet (mantle plus core) and the radius of the core: for a specific interior composition or even for a range of possible compositions, the core radius is expected to be determined with a precision decreasing to a few tens of kilometers.A fairly precise measurement of variations in the orientation of Mars’ spin axis will enable, in addition to the determination of the moment of inertia of the core, an even better determination of the size of the core via the core resonance in the nutation amplitudes. When the core is liquid, the free core nutation (FCN) resonance induces a change in the nutation amplitudes, with respect to their values for a solid planet, at the percent level in the large semi-annual prograde nutation amplitude and even more (a few percent, a few tens of percent or more, depending on the FCN period) for the retrograde ter-annual nutation amplitude. The resonance amplification depends on the size, moment of inertia, and flattening of the core. For a large core, the amplification can be very large, ensuring the detection of the FCN, and determination of the core moment of inertia.The measurement of variations in Mars’ rotation also determines variations of the angular momentum due to seasonal mass transfer between the atmosphere and ice caps. Observations even for a short period of 180 days at the surface of Mars will decrease the uncertainty by a factor of two with respect to the present knowledge of these quantities (at the 10% level).The ultimate objectives of the proposed experiment are to obtain information on Mars’ interior and on the sublimation/condensation of CO₂ in Mars’ atmosphere. Improved knowledge of the interior will help us to better understand the formation and evolution of Mars. Improved knowledge of the CO₂ sublimation/condensation cycle will enable better understanding of the circulation and dynamics of Mars’ atmosphere.     

Mars geodesy, rotation and gravity

This review provides explanations of how geodesy, rotation and gravity can be addressed using radioscience data of an orbiter around a planet or of the lander on its surface.The planet Mars is the center of the discussion.The information one can get from orbitography and radioscience in general concerns the global static gravitational field, the time variation of the gravitational field induced by mass exchange between the atmosphere and the ice caps, the time variation of the gravitational field induced by...     

Mars geodesy,rotation and gravity

This review provides explanations of how geodesy, rotation and gravity can be addressed using radioscience data of an orbiter around a planet or of the lander on its surface.The planet Mars is the center of the discussion.The information one can get from orbitography and radioscience in general concerns the global static gravitational field, the time variation of the gravitational field induced by mass exchange between the atmosphere and the ice caps, the time variation of the gravitational field induced by the tides, the secular changes in the spacecraft's orbit induced by the little moons of Mars named Phobos and Deimos, the gravity induced by particular targets, the Martian ephemerides, and Mars' rotation and orientation.The paper addresses as well the determination of the geophysical parameters of Mars and, in particular, the state of Mars' core and its size, which is important for understanding the planet's evolution.Indeed, the state and dimension of the core determined from the moment of inertia and nutation depend in turn on the percentage of light elements in the core as well as on the core temperature, which is related to heat transport in the mantle.For example, the radius of the core has implications for possible mantle convection scenarios and, in particular, for the presence of a perovskite phase transition at the bottom of the mantle.This is also important for our understanding of the large volcanic province Tharsis on the surface of Mars.     

Numerical simulations of a Mars geodesy network experiment: Effect of orbiter angular momentum desaturation on Mars' rotation estimation

The scientific objectives of a geodetic experiment based on a network of landers, such as NEIGE (NEtlander Ionosphere and Geodesy Experiment) are to improve the current knowledge of Mars' interior and atmosphere dynamics. Such a network science experiment allows monitoring the motions of the Martian rotation axis with a precision of a few centimeters (or milli-arc-seconds (mas)) over annual and sub-annual periods. Thereto, besides radio tracking of a Mars orbiter from the Earth, radio Doppler shifts between this orbiter and several landers at the planet's surface will be performed. From the analysis of these radio Doppler data, it is possible to reconstruct the orbiter motion and Mars' orientation in space. The errors on the orbit determination (position and velocity of the orbiter) have an impact on the geodetic parameters determination from the Doppler shifts and must be removed from the signal in order to achieve a high enough accuracy. In this paper, we perform numerical simulations of the two Doppler signals involved in such an experiment to estimate the impact of the spacecraft angular momentum desaturations on the determination of Mars' orientation variations. The attitude control of the orbiter needs such desaturation maneuvers regularly repeated. They produce velocity variations that must be taken into account when determining the orbit. For our simulations, we use a priori models of the Martian rotation and introduce the spacecraft velocity variations induced by each desaturation event. By a least-squares adjustment of the simulated Doppler signals, we then estimate the orbiter velocity variations and the parameters of the Mars' rotation model. We show that these velocity variations are ill resolved when the spacecraft is not tracked, therefore requiring a near-continuous tracking from the Earth to accurately determine the orbit. In such conditions we show that only 15- of lander-orbiter tracking per week allows recovering Mars' orientation parameters with a precision of a few mas over a period of 1 Martian year.     

Tidally induced surface displacements, external potential variations, and gravity variations on Mars

We have used and extended Roosbeek’s tidal potential for Mars to calculate tidal displacements, gravity variations, and external gravitational potential variations. The tides on Mars are caused by the Sun, and to a lesser degree by the natural satellites Phobos (8%, relative to the Sun) and Deimos (0.08%, relative to the Sun). To determine the reaction of Mars to the tidal forcing, the Love numbers h, l, and k and the gravimetric factor δ were calculated for interior models of Mars with different state, density, and radius of the core and for models which include mantle anelasticity. The latitude dependence and frequency dependence of the Love numbers have been taken explicitly into account. The Love numbers are about three times smaller than those for the Earth and are very sensitive to core changes; e.g., a difference of about 30% is found between a model with a liquid core and an otherwise similar model with a solid core. Tidal displacements on Mars are much smaller than on Earth due to the smaller tidal potential, but also due to the smaller reaction of Mars (smaller Love numbers). For both the tidal diplacement and the tidal external potential perturbations, the tidal signal is at the limit of detection and is too small to permit properties of Mars’s interior to be inferred. On the other hand, the Phobos tidally induced gravity changes, which are subdiurnal with typical periods shorter than 12 h, can be measured very precisely by the very broad band seismometer with thermal control of the seismological experiment SEIS of the upcoming NetLander mission. It is shown that the Phobos-induced gravity tides could be used to study the Martian core.     

Geodesy constraints on the interior structure and composition of Mars

Knowledge of the interior structure of Mars is of fundamental importance to the understanding of its past and present state as well as its future evolution. The most prominent interior structure properties are the state of the core, solid or liquid, its radius, and its composition in terms of light elements, the thickness of the mantle, its composition, the presence of a lower mantle, and the density of the crust. In the absence of seismic sounding only geodesy data allow reliably constraining the deep interior of Mars. Those data are the mass, moment of inertia, and tides. They are related to Mars’ composition, to its internal mass distribution, and to its deformational response to principally the tidal forcing of the Sun. Here we use the most recent estimates of the moment of inertia and tidal Love number k₂ in order to infer knowledge about the interior structure of the Mars.We have built precise models of the interior structure of Mars that are parameterized by the crust density and thickness, the volume fractions of upper mantle mineral phases, the bulk mantle iron concentration, and the size and the sulfur concentration of the core. From the bulk mantle iron concentration and from the volume fractions of the upper mantle mineral phases, the depth dependent mineralogy is deduced by using experimentally determined phase diagrams. The thermoelastic properties at each depth inside the mantle are calculated by using equations of state. Since it is difficult to determine the temperature inside the mantle of Mars we here use two end-member temperature profiles that have been deduced from studies dedicated to the thermal evolution of Mars. We calculate the pressure and temperature dependent thermoelastic properties of the core constituents by using equations state and recent data about reference thermoelastic properties of liquid iron, liquid iron-sulfur, and solid iron. To determine the size of a possible inner core we use recent data on the melting temperature of iron-sulfur.Within our model assumptions the geodesy data imply that Mars has no solid inner core and that the liquid core contains a large fraction of sulfur. The absence of a solid inner is in agreement with the absence of a global magnetic field. We estimate the radius of the core to be 1794 ± 65 km and its core sulfur concentration to be 16 ± 2 wt%. We also show that it is possible for Mars to have a thin layer of perovskite at the bottom of the mantle if it has a hot mantle temperature. Moreover a chondritic Fe/Si ratio is shown to be consistent with the geodesy data, although significantly different value are also possible. Our results demonstrate that geodesy data alone, even if a mantle temperature is assumed, can almost not constrain the mineralogy of the mantle and the crust. In order to obtain stronger constraints on the mantle mineralogy bulk properties, like a fixed Fe/Si ratio, have to be assumed.     

Atmospheric angular momentum variations of Earth, Mars and Venus at seasonal time scales

Atmospheric angular momentum variations of a planet are associated with the global atmospheric mass redistribution and the wind variability. The exchange of angular momentum between the fluid layers and the solid planet is the main cause for the variations of the planetary rotation at seasonal time scales. In the present study, we investigate the angular momentum variations of the Earth, Mars and Venus, using geodetic observations, output of state-of-the-art global circulation models as well as assimilated data. We discuss the similarities and differences in angular momentum variations, planetary rotation and angular momentum exchange for the three terrestrial planets. We show that the atmospheric angular momentum variations for Mars and Earth are mainly annual and semi-annual whereas they are expected to be “diurnal” on Venus. The wind terms have the largest contributions to the LOD changes of the Earth and Venus whereas the matter term is dominant on Mars due to the CO₂ sublimation/condensation. The corresponding LOD variations (ΔLOD) have similar amplitudes on Mars and Earth but are much larger on Venus, though more difficult to observe.     

Some problems related to the origin of methane on Mars

The following problems related to the origin of methane on Mars have been considered. (1) Laboratory simulations of the impact phenomena confirm effective heterogeneous chemistry between the products of the fireball. This chemistry lowers the fireball freezing temperature from 2000 to 750 K for methane and to 1100 K for CO/CO₂. Production of methane on Mars by cometary impacts is 0.8% of the total production. A probability that the observed methane on Mars came from impact of a single comet is 0.0011. (2) The PFS observations of variations of methane on Mars require a very effective heterogeneous loss of methane. Heterogeneous effect of dust is half that of the surface rocks. Thermochemical equilibrium requires production, not loss, of methane. Existing kinetic data show a very low efficiency of heterogeneous reactions of methane. Highly reactive superoxide ions generated by the solar UV photons on the martian rocks cannot remove methane. The required efficiency of heterogeneous loss of methane on Mars is higher than that on Earth by a factor of ?1000, although the expected efficiency on Earth is stronger than that on Mars because of the liquid ocean and the abundant oxygen. All these inconsistencies may be removed if variations of the rock reflectivity contribute to the PFS observations of methane on Mars. The PFS data on H₂CO, HCl, HF, and HBr also raise doubts. (3) Although geologic sources of methane are possible, the lack of current volcanism, hydrothermal activity, hot spots, and very low seepage of gases from the interior are not favorable for geologic methane. Any proposed geological source of methane on Mars should address these problems. Some weak points in the suggested geologic sources are discussed. (4) Measurements of ¹³C/¹²C and D/H in methane would be difficult because of the low methane abundance. These ratios are mostly sensitive to a temperature of methane formation and cannot distinguish between biogenic and low-temperature geologic sources. Their analysis requires the carbon isotope ratio in CO₂ on Mars, which is known with the insufficient accuracy, and D/H in water, which is different in the atmosphere, polar caps, regolith and interior. Therefore, the stable isotope ratios may not give a unique answer on the origin of methane. (5) Ethane and propane react with OH much faster than methane. If their production relative to methane is similar to that on Earth, then their expected abundances on Mars are of a few parts per trillion. (6) Loss of SO₂ in the reaction with peroxide on ice is smaller than its gas-phase loss by an order of magnitude. The overall results strengthen the biogenic origin of martian methane and its low variability.     

Characterization of terrestrial exoplanets based on the phase curves and albedos of Mercury, Venus and Mars

The empirically derived phase curves of terrestrial planets strongly distinguish between airless Mercury, cloud-covered Venus, and the intermediate case of Mars. The function for Mercury is steeply peaked near phase angle zero due to powerful backscattering from its surface, while that for Venus has 100 times less contrast and exhibits a brightness excess near 170° due to Mie scattering from droplets in the atmosphere. The phase curve of Mars falls between those of Mercury and Venus, and there are variations in luminosity due to the planet’s rotation, seasons, and atmospheric states. The phase function and geometric albedo of the Earth are estimated from published albedos values. The curves for Mercury, Venus and Mars are compared to that of the Earth as well as theoretical phase functions for giant planets. The parameters of these different phase functions can be used to characterize exoplanets.     

Numerical assessment of the effects of topography and crustal thickness on martian seismograms using a coupled modal solution-spectral element method

The past 4 decades of Mars exploration have provided much information about the Mars surface, when its interior structure remains relatively poorly constrained. Today available data are compatible with a large range of model parameters. Seismology is able to provide valuable additional data but the number of seismographs will likely be quite limited, specially in the early-stage of future Mars seismic networks. It is thus of importance to be able to correctly isolate effects induced by the crust structure. Mars topography is characterized by spectacular reliefs like the Tharsis bulge or the Hellas basin and by the so-called “Mars dichotomy”: the north hemisphere is made up of low-altitude plains above a relatively thin crust when the south hemisphere is characterized by a thick crust sustaining high reliefs. The aim of this paper is to study the effects induced on seismograms by the topography of the surface and crust-mantle discontinuities. Synthetic seismograms were computed using the coupled spectral element-modal solution method, which reduces the numerical cost by limiting the use of the spectral element method to the regions where lateral variations, like the presence of a topography, are considered. Due to numerical cost, this study is limited to long period and thus focuses on surface waves, mainly on long period Rayleigh waves. We show that reliefs like the Tharsis bulge or the Hellas basin can induce an apparent velocity anomaly up to 0.5% when only the surface topography is introduced. Apparent anomalies can raise up to 1.0% when the surface topography is fully compensated by a mirror-image topography of the crust-mantle discontinuity. Travel-time of surface wave are systematically increased for seismometers in the north hemisphere of Mars and decreased in the south hemisphere. When comparing effects on seismograms by the Earth and Mars topography, we found them to be larger for the Earth. It is due to the fact that we work with a seismic velocity model of Mars with a mean crust thickness of 110 km when the crust thickness has a mean value of 50 km for the Earth. When changing the Mars model for a thinner crust with a mean thickness of 50 km, effects by the topography on Mars seismograms becomes of the same order when not larger than what is observed on the Earth.     

Possible liquid water origin for Atacama Desert mudflow and recent gully deposits on Mars

Evidence of recent gully activity on Mars has been reported based on the formation of new light toned deposits within the past decade, the origin of which remains controversial. Analogous recent light toned gully features have formed by liquid water activity in the Atacama Desert on Earth. These terrestrial deposits leave no mineralogical trace of water activity but rather show an albedo difference due to particle size sorting within a fine-grained mudflow. Therefore, spectral differences indicating varying mineralogy between a recent gully deposit and the surrounding terrain may not be the most relevant criteria for detecting water flow in arid environments. Instead, variation in particle size between the deposit and surrounding terrain is a possible discriminator to identify a water-based flow. We show that the Atacama deposit is similar to the observed Mars gully deposits, and both are consistent with liquid water activity. The light-toned Mars gully deposits could have formed from dry debris flows, but a liquid water origin cannot be ruled out because not all liquid water flows leave hydrated minerals behind on the surface. Therefore, the Mars deposits could be remnant mudflows that formed on Mars within the last decade.     

Deep versus shallow origin of gravity anomalies, topography and volcanism on Earth, Venus and Mars

The relation between gravity anomalies, topography and volcanism can yield important insights about the internal dynamics of planets. From the power spectra of gravity and topography on Earth, Venus and Mars we infer that gravity anomalies have likely predominantly sources below the lithosphere up to about spherical harmonic degree l=30 for Earth, 40 for Venus and 5 for Mars. To interpret the low-degree part of the gravity spectrum in terms of possible sublithospheric density anomalies we derive radial mantle viscosity profiles consistent with mineral physics. For these viscosity profiles we then compute gravity and topography kernels, which indicate how much gravity anomaly and how much topography is caused by a density anomaly at a given depth. With these kernels, we firstly compute an expected gravity-topography ratio. Good agreement with the observed ratio indicates that for Venus, in contrast to Earth and Mars, long-wavelength topography is largely dynamically supported from the sublithospheric mantle. Secondly, we combine an empirical power spectrum of density anomalies inferred from seismic tomography in Earth’s mantle with gravity kernels to model the gravity power spectrum. We find a good match between modeled and observed gravity power spectrum for all three planets, except for 2?l?4 on Venus. Density anomalies in the Venusian mantle for these low degrees thus appear to be very small. We combine gravity kernels and the gravity field to derive radially averaged density anomaly models for the Martian and Venusian mantles. Gravity kernels for l?5 are very small on Venus below ≈800 km depth. Thus our inferences on Venusian mantle density are basically restricted to the upper 800 km. On Mars, gravity anomalies for 2?l?5 may originate from density anomalies anywhere within its mantle. For Mars as for Earth, inferred density anomalies are dominated by l=2 structure, but we cannot infer whether there are features in the lowermost mantle of Mars that correspond to Earth’s Large Low Shear Velocity Provinces (LLSVPs). We find that volcanism on Mars tends to occur primarily in regions above inferred low mantle density, but our model cannot distinguish whether or not there is a Martian analog for the finding that Earth’s Large Igneous Provinces mainly originate above the margins of LLSVPs.     

Mars

Mars is the fourth planet out from the sun. It is a terrestrial planet with a density suggesting a composition roughly similar to that of the Earth. Its orbital period is 687 days, its orbital eccentricity is 0.093 and its rotational period is about 24 hours. Mars has two small moons of asteroidal shapes and sizes (about 11 and 6 km mean radius), the bigger of which, Phobos, orbits with decreasing semimajor orbit axis. The decrease of the orbit is caused by the dissipation of tidal energy in the Martian mantle. The other satellite, Deimos, orbits close to the synchronous position where the rotation period of a planet equals the orbital period of its satellite and has hardly evolved with time. Mars has a tenous atmosphere composed mostly of CO with strong winds and with large scale aeolian transport of surface material during dust storms and in sublimation-condensation cycles between the polar caps. The planet has a small magnetic field, probably not generated by dynamo action in the core but possibly due to remnant magnetization of crustal rock acquired earlier from a stronger magnetic field generated by a now dead core dynamo. A dynamo powered by thermal power alone would have ceased a few billions of years ago as the core cooled to an extent that it became stably stratified. Mars' topography and its gravity field are dominated by the Tharsis bulge, a huge dome of volcanic origin. Tharsis was the major center of volcanic activity, a second center is Elysium about 100° in longitude away. The Tharsis bulge is a major contributor to the non-hydrostaticity of the planet's figure. The moment of inertia factor together with the mass and the radius presently is the most useful constraint for geophysical models of the Martian interior. It has recently been determined by Doppler range measurements to the Mars Pathfinder Lander to be (Folkner et al. 1997). In addition, models of the interior structure use the chemistry of the SNC meteorites which are widely believed to have originated on Mars. According to the models, Mars is a differentiated planet with a 100 to 200 km thick basaltic crust, a metallic core with a radius of approximately half the planetary radius, and a silicate mantle. Mantle dynamics is essential in forming the elements of the surface tectonics. Models of mantle convection find that the pressure-induced phase transformations of -olivine to -spinel, -spinel to -spinel, and -spinel to perovskite play major roles in the evolution of mantle flow fields and mantle temperature. It is not very likely that the -spinel to perovskite transition is present in Mars today, but a few 100 km thick layer of perovskite may have been present in the lower mantle immediately above the core-mantle boundary early in the Martian history when mantle temperatures were hotter than today. The phase transitions act to reduce the number of upwellings to a few major plumes which is consistent with the bipolar distribution of volcanic centers of Mars. The phase transitions also cause a partial layering of the lower mantle which keeps the lower mantle and the core from extensive cooling over the past aeons. A relatively hot, fluid core is the most widely accepted explanation for the present lack of a self-generated magnetic field. Growth of an inner core which requires sub-liquidus temperatures in the core would have provided an efficient mechanism to power a dynamo up to the present day. Received 10 May 1997     

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.     

Strength of Faults on Mars from MOLA Topography

The stresses which must be maintained on faults bounding the rift topography at Tempe Fossae—the “North Tempe Rift” (NTR)—and Valles Marineris (VM) on Mars are estimated, using a simple elastic model and topographic data from the Mars Orbiter Laser Altimeter (MOLA). The absence of rift-flank uplift at the NTR is consistent with an elastic thickness, T_e, of 20 km or greater at the time of rift formation. The maximum resolved shear stresses on bounding faults due to this topography do not therefore exceed 20 MPa, similar to the inferred strength of terrestrial faults. Elastic thickness estimates at VM are mostly around 50 km or greater. Therefore, for canyon widths of ∼400 km, the bounding faults of VM, if present, must be able to withstand stresses of up to approximately 100 MPa. However, if the fault-controlled sections of the canyons do not exceed 150 km in width, as suggested by geomorphological analysis, the fault strength required is only 20 MPa. Although the maximum resolved shear stresses required to support the topography at VM may need to be greater than the stresses which terrestrial faults can support, at least some faults on Mars are no stronger than similar features on Earth. This observation is consistent with the existence of liquid water in the shallow subsurface of Mars at the time the faults were active. On Venus, plate tectonics is probably prevented by the frictional resistance to motion across strong faults. On Mars, it is more likely that the large thickness of the elastic layer of the lithosphere and the possible positive buoyancy of the crust are responsible for the observed lack of plate tectonics.     

Do ice caves exist on Mars?

We have developed a numerical model for assessing the lifetime of ice deposits in martian caves that are open to the atmosphere. Our model results and sensitivity tests indicate that cave ice would be stable over significant portions of the surface of Mars. Ice caves on Earth commonly occur in lava tubes, and Mars has been significantly resurfaced by volcanic activity during its history, including the two main volcanic provinces, the Tharsis and Elysium rises. These areas, known or suspected of having subsurface caves and related voids are among the most favorable regions for the occurrence of ice stability. The martian ice cave model predicts regions which, if caves occur, would potentially be areas of astrobiological importance as well as possible water sources for future human missions to Mars.     

Variations of the magnetic field near Mars caused by magnetic crustal anomalies

The possible avenues for photoelectron transport were determined during southern hemisphere winter at Mars by using a mapping analysis of the theoretical magnetic field. Magnetic field line tracing was performed by superposing two magnetic field models: (1) magnetic field derived from a three-dimensional (3D) self-consistent quasi-neutral hybrid model which does not contain the Martian crustal magnetic anomalies and (2) a 3D map of the magnetic field associated with the magnetic anomalies based on Mars Global Surveyor magnetic field measurements. It was found that magnetic field lines connected to the nightside of the planet are mainly channeled within the optical shadow of the magnetotail whereas magnetic field lines connected to the dayside of the planet are observed to form the remainder of the magnetosphere. The simulation suggests that the crustal anomalies create “a magnetic shield” by decreasing the region near Mars which is magnetically connected to the Martian magnetosphere. The rotation of Mars causes periodic changes in magnetic connectivity, but not to qualitative changes in the overall magnetic field draping around Mars.     

Ionospheric photoelectrons: Comparing Venus, Earth, Mars and Titan

The sunlit portion of planetary ionospheres is sustained by photoionization. This was first confirmed using measurements and modelling at Earth, but recently the Mars Express, Venus Express and Cassini-Huygens missions have revealed the importance of this process at Mars, Venus and Titan, respectively. The primary neutral atmospheric constituents involved (O and CO₂ in the case of Venus and Mars, O and N₂ in the case of Earth and N₂ in the case of Titan) are ionized at each object by EUV solar photons. This process produces photoelectrons with particular spectral characteristics. The electron spectrometers on Venus Express and Mars Express (part of ASPERA-3 and 4, respectively) were designed with excellent energy resolution (ΔE/E=8%) specifically in order to examine the photoelectron spectrum. In addition, the Cassini CAPS electron spectrometer at Saturn also has adequate resolution (ΔE/E=16.7%) to study this population at Titan. At Earth, photoelectrons are well established by in situ measurements, and are even seen in the magnetosphere at up to 7R_E. At Mars, photoelectrons are seen in situ in the ionosphere, but also in the tail at distances out to the Mars Express apoapsis (∼3R_M). At both Venus and Titan, photoelectrons are seen in situ in the ionosphere and in the tail (at up to 1.45R_V and 6.8R_T, respectively). Here, we compare photoelectron measurements at Earth, Venus, Mars and Titan, and in particular show examples of their observation at remote locations from their production point in the dayside ionosphere. This process is found to be common between magnetized and unmagnetized objects. We discuss the role of photoelectrons as tracers of the magnetic connection to the dayside ionosphere, and their possible role in enhancing ion escape.     

Equilibrium rotational stability and figure of Mars

Studies extending over three decades have concluded that the current orientation of the martian rotation pole is unstable. Specifically, the gravitational figure of the planet, after correction for a hydrostatic form, has been interpreted to indicate that the rotation pole should move easily between the present position and a site on the current equator, 90° from the location of the massive Tharsis volcanic province. We demonstrate, using general physical arguments supported by a fluid Love number analysis, that the so-called non-hydrostatic theory is an inaccurate framework for analyzing the rotational stability of planets, such as Mars, that are characterized by long-term elastic strength within the lithosphere. In this case, the appropriate correction to the gravitational figure is the equilibrium rotating form achieved when the elastic lithospheric shell (of some thickness LT) is accounted for. Moreover, the current rotation vector of Mars is shown to be stable when the correct non-equilibrium theory is adopted using values consistent with recent, independent estimates of LT. Finally, we compare observational constraints on the figure of Mars with non-equilibrium predictions based on a large suite of possible Tharsis-driven true polar wander (TPW) scenarios. We conclude, in contrast to recent comparisons of this type based on a non-hydrostatic theory, that the reorientation of the pole associated with the development of Tharsis was likely less than 15° and that the thickness of the elastic lithosphere at the time of Tharsis formation was at least ∼50 km. Larger Tharsis-driven TPW is possible if the present-day gravitational form of the planet at degree 2 has significant contributions from non-Tharsis loads; in this case, the most plausible source would be internal heterogeneities linked to convection.     

Inversion of relief on Mars

Relief inversion has been invoked to explain a number of geomorphic features of the martian surface. Terrestrial relief inversion occurs when former depressions become elevated because their fill is more resistant to erosion than the surrounding terrain. It is a common product of long-term landscape evolution on Earth, especially in relatively stable intra-cratonic settings and flat, or near flat lying successions. The inverted relief will preserve relicts of former land surfaces and is therefore older than the surrounding terrain. Relief inversion can occur by a range of processes, including infill of depressions by intrinsically resistant material, selective secondary cementation via diagenesis and weathering, or surface armouring. We examine a number of possible cases of inverted relief on Mars that appear to have formed by these three processes. We suggest that the most likely cementing agents for surface induration are iron oxides, silica, and sulfates. Possible cementation mechanisms include fluid mixing during regional groundwater flow, cooling of hydrothermal or basinal fluids as they near the surface, and evaporation and sublimation of near surface water. Wind action appears the most common erosive process on Mars capable of the regional landscape lowering necessary for relief inversion to occur, unlike on Earth where both deflation and runoff are important. Preliminary crater densities of selected features show that the tops of the proposed inverted relief have considerably more craters than the surrounding plains, as is predicted by the inversion hypothesis. More accurate dating of inverted surfaces and the adjacent areas may provide a simple way of measuring the degree of erosion over time in at least some areas of Mars.