The early martian evolution—Constraints from basin formation ages

Impact basin formation ages give insight into the early evolution of a planet. The martian basins Hellas, Isidis and Argyre provide an important time-marker for the cessation of the magnetic dynamo and the crustal thickness distribution, both established before 4 Ga ago. No martian surfaces are older than 4.15 Ga based on crater count statistics, and all are younger than the oldest lunar ones. I show that the heavy bombardment period on the Moon and Mars evolved similarly, but endogenic processes have removed the oldest martian basin record. The basin-forming projectile population appears to be different from the impactor population observed today in the inner Solar System. It is yet uncertain whether the heavy bombardment period is cataclysmic or characterized by the decaying flux of planetary formation.     

A comparison of global models for the solar wind interaction with Mars

We present initial results from the first community-wide effort to compare global plasma interaction model results for Mars. Seven modeling groups participated in this activity, using MHD, multi-fluid, and hybrid assumptions in their simulations. Moderate solar wind and solar EUV conditions were chosen, and the conditions were implemented in the models and run to steady state. Model output was compared in three ways to determine how pressure was partitioned and conserved in each model, the location and asymmetry of plasma boundaries and pathways for planetary ion escape, and the total escape flux of planetary oxygen ions. The two participating MHD models provided similar results, while the five sets of multi-fluid and hybrid results were different in many ways. All hybrid results, however, showed two main channels for oxygen ion escape (a pickup ion ‘plume’ in the hemisphere toward which the solar wind convection electric field is directed, and a channel in the opposite hemisphere of the central magnetotail), while the MHD models showed one (a roughly symmetric channel in the central magnetotail). Most models showed a transition from an upstream region dominated by plasma dynamic pressure to a magnetosheath region dominated by thermal pressure to a low altitude region dominated by magnetic pressure. However, calculated escape rates for a single ion species varied by roughly an order of magnitude for similar input conditions, suggesting that the uncertainties in both the current and integrated escape over martian history as determined by models are large. These uncertainties are in addition to those associated with the evolution of the Sun, the martian dynamo, and the early atmosphere, highlighting the challenges we face in constructing Mars’ past using models.     

On the effect of the martian crustal magnetic field on atmospheric erosion

Without the shielding of a strong intrinsic magnetic field, the martian atmosphere directly interacts with the impacting solar wind. The neutral constituents of the atmospheric corona can be ionized, and then picked up and accelerated by the magnetic field and convection electric field in the solar wind. A significant fraction of pickup ions escape Mars’ gravitational pull and are lost to space. This non-thermal escape process of heavy species is an important mechanism responsible for atmospheric erosion. While there is a perception that the martian magnetic anomalies are significant for the ionospheric density distribution and the bow shock standoff location, little is known about the quantitative influence of the martian crustal magnetic field on the global distribution of escaping pickup ions. In this paper, we apply a newly developed Monte Carlo ion transport model to resolve the crustal field effect on the pickup oxygen ion distribution around Mars. The background magnetic and electric fields, in which test particles are followed, are calculated using an independent three-dimensional multispecies MHD model. The effects of the crustal magnetic field on particle escape are quantified by varying the crustal field orientation in the model setup and comparing the corresponding test particle simulation results. The comparison is made by turning on or off the crustal field or changing the local time of the strongest field from the dayside to the dawnside. It is found that without the protection of the crustal magnetic field, the total amount of atmospheric escape through the tail region would be enhanced by more than a factor of two. It is shown that the crustal magnetic field not only regionally deflects the solar wind around the martian atmosphere, but also has an important global effect on atmospheric erosion and thus on long-term atmospheric evolution.     

Comparison of small lunar landslides and martian gullies

Some lunar crater-wall landslides strongly resemble martian gullies, despite the lack of geologically active water on the Moon today or in the past. The lunar features indicate that alcove-channel-apron morphology, attributed on Mars to seepage of liquid water, can also form via a dry landslide mechanism. Therefore a more stringent test than just an alcove-channel-apron morphology is necessary to differentiate dry landslides from water carved gullies.     

Electric fields within the martian magnetosphere and ion extraction: ASPERA-3 observations

Observations made by the ASPERA-3 experiment onboard the Mars Express spacecraft found within the martian magnetosphere beams of planetary ions. In the energy (E/q)-time spectrograms these beams are often displayed as dispersive-like, ascending or descending (whether the spacecraft moves away or approach the planet) structures. A linear dependence between energy gained by the beam ions and the altitude from the planet suggests their acceleration in the electric field. The values of the electric field evaluated from ion energization occur close to the typical values of the interplanetary motional electric field. This suggests an effective penetration of the solar wind electric field deep into the martian magnetosphere or generation of large fields within the magnetosphere. Two different classes of events are found. At the nominal solar wind conditions, a ‘penetration’ occurs near the terminator. At the extreme solar wind conditions, the boundary of the induced magnetosphere moves to a more dense upper atmosphere that leads to a strong scavenging of planetary ions from the dayside regions.     

Radar absorption due to a corotating interaction region encounter with Mars detected by MARSIS

Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) is a subsurface and topside ionosphere radar sounder aboard the European Space Agency spacecraft Mars Express, in orbit at Mars since 25 December 2003, and in operation since 17 June 2005. The ionospheric sounding mode of MARSIS is capable of detecting the reflection of the sounding wave from the martian surface. This ability has been used in previous work to show that the surface reflection is absorbed and disappears during periods when high fluxes of energetic particles are incident on the ionosphere of Mars. These absorption events are believed to be the result of increased collisional damping of the sounding wave, caused by increased electron density below the spacecraft, in turn caused by impact ionization from the impinging particles. In this work we identify two absorption events that were isolated during periods when the surface reflection is consistently visible and when Mars is nearly at opposition. The visibility of the surface reflection is viewed in conjunction with particle and photon measurements taken at both Mars and Earth. Both absorption events are found to coincide with Earth passing through solar wind speed and ion flux signatures indicative of a corotating interaction region (CIR). The two events are separated by an interval of approximately 27 days, corresponding to one solar rotation. The first of the two events coincides with abruptly enhanced particle fluxes seen in situ at Mars. Simultaneous with the particle enhancement there are an abrupt decrease in the intensity of electron oscillations, typically seen by the Mars Express particle instrument ASPERA-3 between the magnetic pileup boundary and the martian bow shock, and a sharp drop in the solar wind pressure, seen in the proxy quantity based on MGS magnetometer observations. The decrease in oscillation intensity is therefore the probable effect of a relaxation of the martian bow shock. The second absorption event does not show a particle enhancement and complete ASPERA-3 data during that time are unavailable. Other absorption events are the apparent result of solar X-ray and XUV enhancements. We conclude that surface reflection absorption events are sometimes caused by enhanced ionospheric ionization from high energy particles accelerated by the shocks associated with a CIR. A full statistical analysis of CIRs in relation to observed absorption events in conjunction with a quantitative analysis of the deposition of ionization during space weather events is needed for a complete understanding of this phenomenon. If such analyses can be carried out, radar sensing of the martian ionosphere might be useful as a space weather probe.     

Simulating OH/H2O formation by solar wind at the lunar surface

We simulate the OH/H₂O production from the action of keV protons on the lunar regolith using a vacuum chamber and a mass analyzer to examine the molecular products released from olivine and SiO₂ powders during their irradiation by deuterium ions. The measured mass spectra, showing the OD/D₂O signature, confirm the possibility of OH/H₂O formation on the lunar surface by solar-wind hydrogen.     

First ENA observations at Mars: ENA emissions from the martian upper atmosphere

The neutral particle detector (NPD) on board Mars Express has observed energetic neutral atoms (ENAs) from a broad region on the dayside of the martian upper atmosphere. We show one such example for which the observation was conducted at an altitude of 570 km, just above the induced magnetosphere boundary (IMB). The time of flight spectra of these ENAs show that they had energies of 0.2-2 keV/amu, with an average energy of ∼1.1 keV/amu. Both the spatial distribution and the energy of these ENAs are consistent with the backscattered ENAs, produced by an ENA albedo process. This is the first observation of backscattered ENAs from the martian upper atmosphere. The origin of these ENAs is considered to be the solar wind ENAs that are scattered back by collision processes in the martian upper atmosphere. The particle flux and energy flux of the backscattered ENAs are and , respectively.     

Accretion, loss, and fractionation of martian water

A recently published model for accretion of water by Mars from asteroids and comets beyond 2.5 AU calls for Mars to have accreted 6-27% of Earth's ocean (equivalent to 600-2700 meters depth on Mars' surface) with a D/H ratio 1.2-1.6 times Standard Mean Ocean Water. This model is shown to provide enough outgassed water to account for the subsequent evolution of D/H because of loss of hydrogen and the D/H ratio in water of hydration found in martian meteoric minerals, but enhanced rate loss of hydrogen from early Mars is necessary to permit adequate fractionation to occur during the lifetime of the planet.     

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.     

A 3D hybrid simulation study of the electromagnetic field distributions in the lunar wake

A three-dimensional hybrid code is used to study the electromagnetic disturbances in the solar wind that arise due to the absorption effect of the Moon. Due to the nearly insulating nature of the Moon, interplanetary magnetic fields (IMFs) can move through the interior without hindrance. However, the near-vacuum created in the wake region due to the lunar absorption effect will lead to enhancement of the strength of the magnetic field by a factor of about 1.4 in the middle of the lunar wake and lead to depletions at two sides. The situations arising from different orientations of the interplanetary magnetic fields relative to the radial direction are compared. Asymmetries of the inward diffusions both along and perpendicular to the field lines are also observed. The electric field formed from the plasma convection could reach a magnitude of 0.2–0.8 mV/m at the border of the wake. The role of the electric field on the inward accelerations is important to the geometry of the lunar wake.     

A plasma flow velocity boundary at Mars from the disappearance of electron plasma oscillations

The Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) on the Mars Express (MEX) spacecraft is capable of measuring ionospheric electron density by the use of two main methods: remote radar sounding and from the excitation of local plasma oscillations. The frequency of the locally excited electron plasma oscillations is used to measure the local electron density. However, plasma oscillations are not observed when the plasma flow velocity is higher than about 160 km/s, which occurs mainly in the solar wind and magnetosheath. As a consequence, in many passes, there is a sudden disappearance of the plasma oscillations as the spacecraft enters into the magnetosheath. This fact allows us to identify a flow velocity boundary on the dayside, between the ionosphere of Mars and the shocked solar wind. This paper summarizes the results of the measurement of 552 orbits mostly over a period from August 4, 2005 to August 17, 2007. The boundary points found using MARSIS have been verified by measurements from the Analyzer of Space Plasma and Energetic Atoms (ASPERA-3) Electron Spectrometer (ELS) instrument on Mars Express. The average position of the flow velocity boundary is compared to flow velocity simulations computed using hybrid model and other boundaries. The boundary altitude is slightly lower than the magnetic pile-up boundary determined using Phobos 2 and Mars Global Surveyor (MGS) crossings, but it is in good agreement with the induced magnetospheric boundary determined by ASPERA-3. Investigation of the effect of the crustal magnetic field revealed that the flow velocity boundary is raised at the locations with strong crustal magnetic fields.     

Lofted charged dust distribution above the Moon surface

We developed kinetic theory for the charging processes of small dust grains near the lunar surface due to interaction with the anisotropic solar wind plasma. Once charged, these dust grains, which are exposed to the electric field in the sheath region near the lunar surface, could loft and distribute around such heights off the surface where they reach equilibrium with the local gravitational force. Analytical solutions were derived for the charging time, grain floating potential, and grain charge, characterizing the charging processes of small dust grains in a two-component and in a multi-component solar wind plasma, and further highlighting the unique features presented by the high streaming plasma velocity. We have also formulated a novel kinetic theory of sheath formation around an absorbing planar surface immersed in the anisotropic solar wind plasma in the case of a negligible photoelectric effect and presented solutions for the sheath structure. In this study we combined the results from these analyses and provided estimates for the size distribution function of dust that is expected to be lofted in regions dominated by the solar wind plasma, such as near the terminator and in nearby shadowed craters. Corresponding to the two dominant streaming velocity peaks of 300 and 800 km/s, mean dust diameters of 500 and 350 nm, respectively, are expected to be found at equilibrium at heights of relevance to exploration operations, e.g., around 1.5 m height off the lunar surface. In shadowed craters near the terminator region, where isotropic plasma should be dominating, we estimate mean lofted dust diameter of 800 nm around the same 1.5 m height off the lunar surface. The generally applicable solutions could be used to readily calculate the expected lofted size distribution near the lunar surface as a function of plasma parameters, dust grain composition, and other parameters of interest.     

Seasonal melting of surface water ice condensing in martian gullies

In this work we consider when and how much liquid water during present climate is possible within the gullies observed on the surface of Mars. These features are usually found on poleward directed slopes. We analyse the conditions for melting of H₂O ice, which seasonally condenses within the gullies. We follow full annual cycle of condensation and sublimation of atmospheric CO₂ and H₂O, accounting for the heat and mass transport in the soil. During the summer, once the facets of the gullies are exposed to the Sun the water ice can melt and evaporate. Two mid latitude locations in both hemispheres are considered. The model includes both the rough geometry of the gullies as well as the slope of the surface where the gullies appear. It is an extension of the model developed to calculate condensation of CO₂ ice in troughs of different sizes, including polygonal features on Mars (Kossacki and Markiewicz, 2002, Icarus 160, 73; Kossacki et al., 2003, Planet. Space Sci. 51, 569). We have found, that water ice accumulated during winter can undergo transition to the liquid phase after complete sublimation of CO₂ ice. The amount of liquid water depends on water content in the atmosphere and on the local wind speed. It is probably not enough to destabilise the slope and cause flow of the surface material. However, even the small amounts of liquid water predicted, can play an important role in surface chemistry, in increasing the cohesive strength of the soil's surface layer and possibly may have some exobiological implications.     

Structure of the martian wake

We present the first results from the ion mass analyzer IMA of the ASPERA-3 instrument on-board of Mars Express. More than 200 orbits for May 2004-September 2004 time interval have been selected for the statistical study of the distribution of the atmospheric origin ions in the planetary wake. This study shows that the martian magnetotail consists of two different ion regimes. Planetary origin ions of the first regime form the layer adjacent to the magnetic pile-up boundary. These ions are accelerated to energy greater than 2000 eV and exhibit a gradual decreasing of energy down to the planetary tail. The second plasma regime is observed in the planetary shadow. The heavy ions (considered as planetary ones) are accelerated to the energy of the solar wind protons. Obviously the acceleration mechanism is different for the different plasma regimes. Study of two plasma regimes in the frame referred to the interplanetary magnetic field (IMF) direction (we used MGS magnetometer data to obtain the IMF clock angle) clearly shows their spatial anisotropy. The monoenergetic plasma in the planetary shadow is observed only in the narrow angular sector around the positive direction of the interplanetary electric field.     

Monte-Carlo simulations of atmospheric muon production: Implication of the past martian environment

In order to evaluate the obliquity-driven atmospheric-density path length effect on nuclide production rate on Mars, we performed a Monte-Carlo simulation to produce the number of secondary particles such as muons, neutrons and protons in the martian atmosphere and to simulate that production of ¹⁰Be and ³⁶Cl in the martian regolith by muons and neutrons depends on how much atmosphere had been present for the past 10 million years. The vertical profile of the present martian atmosphere to generate secondary particles has been determined based on the data provided by the Viking missions. For other thickness profiles, we scaled Linsley's atmospheric model. Atmospheric shower has been generated with the SIBYLL 2.1 for high-energy hadronic interactions and EHSA for low energy photonuclear interactions. With increasing atmospheric thickness, more primary interactions occur in the atmosphere. Consequently the proton flux is reduced and the secondary cosmic ray flux, such as muons or energetic neutrons increases at surface. The result indicates that the muon production is more sensitive to obliquity-driven atmospheric variations than proton reduction. A thicker atmosphere would result in enhanced nuclide production at a place deeper than 5 m below the surface and the nuclides present in detectable concentrations. Application to the polar deposit is described.     

First ENA observations at Mars: Solar-wind ENAs on the nightside

We present measurements with an Energetic Neutral Atom (ENA) imager on board Mars Express when the spacecraft moves into Mars eclipse. Solar wind ions charge exchange with the extended Mars exosphere to produce ENAs that can spread into the eclipse of Mars due to the ions' thermal spread. Our measurements show a lingering signal from the Sun direction for several minutes as the spacecraft moves into the eclipse. However, our ENA imager is also sensitive to UV photons and we compare the measurements to ENA simulations and a simplified model of UV scattering in the exosphere. Simulations and further comparisons with an electron spectrometer sensitive to photoelectrons generated when UV photons interact with the spacecraft suggest that what we are seeing in Mars' eclipse are ENAs from upstream of the bow shock produced in charge exchange with solar wind ions with a non-zero temperature. The measurements are a precursor to a new technique called ENA sounding to measure solar wind and planetary exosphere properties in the future.     

First ENA observations at Mars: Charge exchange ENAs produced in the magnetosheath

Measurements of energetic neutral atoms (ENA) generated in the magnetosheath at Mars are reported. These ENAs are the result of charge exchange collisions between solar wind protons and neutral oxygen and hydrogen in the exosphere of Mars. The peak of the observed ENA flux is . For the case studied here, i.e., the passage of Mars Express through the martian magnetosheath around 20:15 UT on 3 May 2004, the measurements agree with an analytical model of the ENA production at the planet. It is possible to find parameter values in the model such that the observed peak in the ENA count rate during the spacecraft passage through the magnetosheath is reproduced.     

The magnetic field near Mars: A comparison between a hybrid model, Mars Global Surveyor and Mars Express observations

Mars Express (MEX) does not carry its own magnetometer which complicates interpretation of ASPERA-3/MEX ion measurements. The direction of the interplanetary magnetic field (IMF) is especially important because it, among other things, determines the direction of the convective electric field and orientation of the cross tail current sheet and tail lobes. In this paper we present a case study to show the properties of the magnetic field near Mars in a quasi-neutral hybrid (QNH) model at the orbits where the Mars Global Surveyor (MGS) has made measurements, present a method to derive the IMF clock angle by comparing fields in a hybrid model and the direction of the magnetic field measured by MGS by deriving the IMF clock angle. We also use H⁺ ring velocity distribution observations upstream of the bow shock measured by the IMA/ASPERA-3 instrument on board MEX spacecraft. These observations are used to indirectly provide the orientation of the IMF. We use a QNH model (HYB-Mars) where ions are modeled as particles while electrons form a mass-less charge neutralizing fluid. We found that the direct MGS and non-direct IMA observations of the orientation magnetic field vectors in non-crustal magnetic field regions are consistent with the global magnetic field draping pattern predicted by the global model.     

Electron reflectometry in the martian atmosphere

The technique of electron reflectometry, a method for remote estimation of planetary magnetic fields, is expanded from its original use of mapping crustal magnetic fields at the Moon to achieving the same purpose at Mars, where the presence of a substantial atmosphere complicates matters considerably. The motion of solar wind electrons, incident on the martian atmosphere, is considered in detail, taking account of the following effects: the electrons' helical paths around the magnetic field lines to which they are bound, the magnetic mirror force they experience due to converging field lines in the vicinity of crustal magnetic anomalies, their acceleration/deceleration by electrostatic potentials, their interactions with thermal plasma, their drifts due to magnetic field line curvature and perpendicular electric fields and their scattering off, and loss of energy through a number of different processes to, atmospheric neutrals. A theoretical framework is thus developed for modeling electron pitch angle distributions expected when a spacecraft is on a magnetic field line which is connected to both the martian crust and the interplanetary magnetic field. This framework, along with measured pitch angle distributions from the Mars Global Surveyor (MGS) Magnetometer/Electron Reflectometer (MAG/ER) experiment, can be used to remotely measure crustal magnetic field magnitudes and atmospheric neutral densities at ∼180 km above the martian datum, as well as estimate average parallel electric fields between 200 and 400 km altitude. Detailed analysis and full results, concerning the crustal magnetic field and upper thermospheric density of Mars, are left to two companion papers.