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.     

Solar wind plasma protrusion into the martian magnetosphere: ASPERA-3 observations

The ASPERA-3 experiment onboard the Mars Express spacecraft revealed, near the wake boundary of Mars, a spatially narrow, strip-like plasma structure composed of magnetosheath-like electrons and planetary ions. The peak electron energy often exceeds the peak energy at the bow shock that indicates a significant heating (acceleration) during the structure formation. It is shown that this structure is formed during efficient plasma penetration into the martian magnetosphere in the region near the terminator. The penetration of sheath electrons and their gradual heating (acceleration) is accompanied by a change of the ion composition from a solar wind plasma to a planetary plasma dominated by oxygen ions. A possible mechanism of plasma inflow to the magnetosphere is discussed.     

Asymmetry of plasma fluxes at Mars. ASPERA-3 observations and hybrid simulations

The asymmetry of fluxes of solar wind and planetary ions is studied by using the ASPERA-3 observations onboard the Mars Express spacecraft in February 2004 to March 2006. Due to the small scale of the Martian magnetosphere and its induced origin, the flow pattern near Mars is sensitive to the directions of the interplanetary magnetic and electric (-V×B) fields. Asymmetry of the magnetic field draping produces an asymmetry in plasma flows in the plane containing the IMF. The crustal magnetic fields on Mars also influence the flow pattern. Scavenging of planetary ions is less efficient in the regions of strong crustal magnetization and therefore the escape fluxes of planetary ions in the southern hemisphere are smaller. The results of the observations are compared to simulations based on a 3D hybrid model with several ion species.     

Observations of magnetic anomaly signatures in Mars Express ASPERA-3 ELS data

Mars Express (MEX) Analyser of Space Plasmas and Energetic Atoms (ASPERA-3) data is providing insights into atmospheric loss on Mars via the solar wind interaction. This process is influenced by both the interplanetary magnetic field (IMF) in the solar wind and by the magnetic ‘anomaly’ regions of the martian crust. We analyse observations from the ASPERA-3 Electron Spectrometer near to such crustal anomalies. We find that the electrons near remanent magnetic fields either increase in flux to form intensified signatures or significantly reduce in flux to form plasma voids. We suggest that cusps intervening neighbouring magnetic anomalies may provide a location for enhanced escape of planetary plasma. Initial statistical analysis shows that intensified signatures are mainly a dayside phenomenon whereas voids are a feature of the night hemisphere.     

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.     

Phasespace transport of a quasi-neutral multi-fluid plasma over the solar wind MHD termination shock

We study the physics of a multi-ion MHD shock, i.e. an MHD shock feature that forms when a supersonic flow of mixed ion populations is forced to adapt itself to a pressure obstacle further downstream. We shall describe this situation by using a multi-fluid approach for a mixture of ion populations with different specific masses and charges per ion species. First we calculate the effective electric potential that forces the plasma bulk to decelerate to the downstream bulk flow velocity which also then defines that system into which the downstream magnetic field is frozen-in. Then we calculate the unavoidable ion-specific overshoot velocities and gain from them, requesting energy conservation, the ion-specific contributions to the downstream thermal energies and pressures. The aim thereby is to find the solution for the MHD status of the downstream flow of the plasma mixture, specifically for a proton-electron plasma. We derive an implicit equation for the effective compression ratio and explicit relations for the different, downstream ion and electron temperatures as function of the multi-fluid compression ratio s. The resulting actual multi-fluid compression ratio s _eff is found by adding up all the partial downstream pressures and comparing it with the upstream ram pressure. As we can show, the electron pressure is the dominant contribution to the total downstream plasma pressure.     

Ion escape from Mars as a function of solar wind conditions: A statistical study

The influence of solar EUV and solar wind conditions on ion escape at Mars is investigated using ion data from the Aspera-3 instrument on Mars Express, combined with solar wind proxy data obtained from the Mars Global Surveyor (MGS) spacecraft. A solar EUV flux proxy based on data from the Earth position, scaled and shifted in time for Mars, is used to study relatively long time scale changes related to solar EUV variability. Data from May 2004 until November 2005 has been used. A clear dependence on the strength of the subsolar magnetic field as inferred from MGS measurements is seen in the ion data. The region of significant heavy ion flows is compressed and the heavy ion flux density is higher for high subsolar magnetic field strength. Because of the difference in outflow area, the difference in estimated total outflow is somewhat less than the difference in average flux density. We confirm previous findings that escaping planetary ions are mainly seen in the hemisphere into which the solar wind electric field is pointed. The effect is more pronounced for the high subsolar magnetic field case.The average ion motion has a consistent bias towards the direction of the solar wind electric field, but the main motion is in the antisunward direction. The antisunward flow velocity increases with tailward distance, reaching above at 2 to 3 martian radii downtail from Mars for O⁺ ions. Different ion species reach approximately the same bulk flow energy. We did not find any clear correlation between the solar EUV flux and the ion escape distribution or rate, probably because the variation of the solar EUV flux over our study interval was too small. The results indicate that the solar wind and its magnetic field directly interacts with the ionosphere of Mars, removing more ions for high subsolar magnetic field strength. The interaction region and the tail heavy ion flow region are not perfectly shielded from the solar wind electric field, which accelerates particles over relatively large tail distances.     

Ionospheric plasma acceleration at Mars: ASPERA-3 results

The Analyzer of Space Plasma and Energetic Atoms (ASPERA) on-board the Mars Express spacecraft (MEX) measured penetrating solar wind plasma and escaping/accelerated ionospheric plasma at very low altitudes (250 km) in the dayside subsolar region. This implies a direct exposure of the martian topside atmosphere to solar wind plasma forcing leading to energization of ionospheric plasma. The ion and electron energization and the ion outflow from Mars is surprisingly similar to that over the magnetized Earth. Narrow “monoenergetic” cold ion beams, ion beams with broad energy distributions, sharply peaked electron energy spectra, and bidirectional streaming electrons are particle features also observed near Mars. Energized martian ionospheric ions (O⁺, O⁺₂, CO⁺₂, etc.) flow in essentially the same direction as the external sheath flow. This suggests that the planetary ion energization couples directly to processes in the magnetosheath/solar wind. On the other hand, the beam-like distribution of the energized plasma implies more indirect energization processes like those near the Earth, i.e., energization in a magnetized environment by waves and/or parallel (to B) electric fields. The general conditions for martian plasma energization are, however, different from those in the Earth's magnetosphere. Mars has a weak intrinsic magnetic field and solar wind plasma may therefore penetrate deep into the dense ionospheric plasma. Local crustal magnetization, discovered by Acuña et al. [Acuña, M.J., Connerey, J., Ness, N., Lin, R., Mitchell, D., Carlsson, C., McFadden, J., Anderson, K., Rème, H., Mazelle, C., Vignes, D., Wasilewski, P., Cloutier, P., 1999. Science 284, 790-793], provide some dayside shielding against the solar wind. On the other hand, multiple magnetic anomalies may also lead to “hot spots” facilitating ionospheric plasma energization. We discuss the ASPERA-3 findings of martian ionospheric ion energization and present evidences for two types of plasma energization processes responsible for the low- and mid-altitude plasma energization near Mars: magnetic field-aligned acceleration by parallel electric fields and plasma energization by low frequency waves.     

Three-dimensional computer modeling of dynamic reconnection in the magnetotail

Magnetospheric substorm in the magnetotail region is studied numerically by means of a three-dimensional MHD code. The analytic solution for the quiet magnetotail is emloyed as an initial configuration. The localized solar wind is modeled to enter the simulation domain through the boundaries located in the magnetotail lobe region. As a result of the interaction between the solar wind and the magnetosphere, the magnetic field lines are stretched, and the plasma sheet becomes thinner and thinner. When the current-driven resistivity is generated, magnetic reconnection is triggered by this resistivity. The resulting plasma jetting is found to be super-magnetosonic. Although the plasmoid formation and its tailward motion is not quite clear as in the two-dimensional simulation, which is mainly because of the numerical model chosen for the present simulation, the rarification of the plasmas near thex-point is observed. Field-aligned currents are observed in the late expansive stage of the magnetospheric substorm. These field-aligned currents flow from the tail toward the ionosphere on the dawn side and from the ionosphere toward the tail on the dusk side, namely in the same sence of the region 1 current. As the field-aligned currents develop, it is found that the cross-tail current in the Earthside midnight section of the magneticx-point is reduced.     

Planetary ENA imaging: Effects of different interaction models for Mars

We present and compare energetic neutral atom (ENA) images that are calculated from plasma parameters given by three different simulation models of the interaction between the solar wind and Mars. The images are calculated by combining a model for the ion flow with a model of the neutral atmosphere using the cross-sections for the charge exchange collisions. The three ion models are: an empirical model that is based on Phobos 2 measurements; a three-dimensional hybrid simulation; and a three-dimensional MHD simulation. For the empirical and MHD models the images are obtained by integration of the ENA emission along lines of sight to a virtual ENA instrument. In the case of the hybrid model images are obtained by summing the contributions from all ions, whose positions, velocities, and weights are saved in files at regular intervals.Differences between the models can be detected in the images, for example the hybrid model produces ENA emissions from a larger region than the MHD model does. An asymmetry in the oxygen ion density develops in the hybrid model and can be seen in the oxygen ENA images. The images are influenced by finite gyro radius effects, which are included in the hybrid model but not in the other two. The total production rates of hydrogen ENAs are , , and for the empirical, hybrid and MHD models respectively.This study shows the importance of considering both the type of simulation model used and the proper inclusion of relevant physical phenomena and boundary conditions, when modelling the interaction between planets and the solar wind. Although the different models agree fairly well in terms of macroscopic plasma parameters they produce ENA images that differ substantially.     

Solar wind-magnetosphere interaction as simulated by a 3-D em particle code

We have studied the solar wind-magnetosphere interaction using a 3-D electromagnetic particle code. The results for an unmagnetized solar wind plasma streaming past a dipole magnetic field show the formation of a magnetopause and a magnetotail, the penetration of energetic particles into cusps and radiation belt and dawn-dusk asymmetries. The effects of interplanetary magnetic field (IMF) have been investigated in a similar way as done by MHD simulations. The simulation results with a southward IMF show the shrunk magnetosphere with great particle entry into the cusps and nightside magnetosphere. This is a signature of a magnetic reconnection at the dayside magnetopause. After a quasi-stable state is established with an unmagnetized solar wind we switched on a solar wind with an northward IMF. In this case the significant changes take place in the magnetotail. The waving motion was seen in the magnetotail and its length was shortened. This phenomena are consistent with the reconnections which occur at the high latitude magnetopause. In our simulations kinetic effects will determine the self-consistent anomalous resistivity in the magnetopause that causes reconnections.Deceased January 24, 1993; R. Bunemanet al. 1993.     

The magnetic field draping direction at Mars from April 1999 through August 2004

Using more than five years of data from the magnetometer and electron reflectometer (MAG/ER) on Mars Global Surveyor (MGS), we derive the draping direction of the magnetic field above a given latitude band in the northern hemisphere. The draping direction varies on timescales associated with the orbital period of Mars and with the solar rotation period. We find that there is a strongly preferred draping direction when Mars is in one solar wind sector, but the opposite direction is not preferred as strongly for the other solar wind sector. This asymmetry occurs at or below the magnetic pileup boundary (MPB), is observed preferentially on field lines that connect to the collisional ionosphere, and is independent of planetary longitude. The observations could be explained by a hemispherical asymmetry in the access of field lines to the low-altitude ionosphere, or possibly from global modification of the low-altitude solar wind interaction by crustal magnetic fields. We show that the draping direction affects both the penetration of sheath plasma to 400 km altitudes on the martian dayside and the radial component of the magnetic field on the planetary night side.     

A petrogenetic model for the origin and compositional variation of the martian basaltic meteorites

Abstract— The major element, trace element, and isotopic compositional ranges of the martian basaltic meteorite source regions have been modeled assuming that planetary differentiation resulted from crystallization of a magma ocean. The models are based on low to high pressure phase relationships estimated from experimental runs and estimates of the composition of silicate Mars from the literature. These models attempt to constrain the mechanisms by which the martian meteorites obtained their superchondritic CaO/Al₂O₃ ratios and their source regions obtained their parent/daughter (⁸⁷Rb/⁸⁶Sr, ¹⁴⁷Sm/¹⁴⁴Nd, and ¹⁷⁶Lu/¹⁷⁷Hf) ratios calculated from the initial Sr, Nd, and Hf isotopic compositions of the meteorites. High pressure experiments suggest that majoritic garnet is the liquidus phase for Mars relevant compositions at or above 12 GPa. Early crystallization of this phase from a martian magma ocean yields a liquid characterized by an elevated CaO/Al₂O₃ ratio and a high Mg#. Olivine‐pyroxene‐garnet‐dominated cumulates that crystallize subsequently will also be characterized by superchondritic CaO/Al₂O₃ ratios. Melting of these cumulates yields liquids with major element compositions that are similar to calculated parental melts of the martian meteorites. Furthermore, crystallization models demonstrate that some of these cumulates have parent/daughter ratios that are similar to those calculated for the most incompatible‐element‐depleted source region (i.e., that of the meteorite Queen Alexandra [QUE] 94201). The incompatible‐element abundances of the most depleted (QUE 94201‐like) source region have also been calculated and provide an estimate of the composition of depleted martian mantle. The incompatible‐element pattern of depleted martian mantle calculated here is very similar to the pattern estimated for depleted Earth's mantle. Melting the depleted martian mantle composition reproduces the abundances of many incompatible elements in the parental melt of QUE 94201 (e.g., Ba, Th, K, P, Hf, Zr, and heavy rare earth elements) fairly well but does not reproduce the abundances of Rb, U, Ta and light rare earth elements. The source regions for meteorites such as Shergotty are successfully modeled as mixtures of depleted martian mantle and a late stage liquid trapped in the magma ocean cumulate pile. Melting of this hybrid source yields liquids with major element abundances and incompatible‐element patterns that are very similar to the Shergotty bulk rock.     

Martian xenon components in Shergotty mineral separates: Locations,sources, and trapping mechanisms

Abstract— Isotopic signatures and concentrations of xenon have been measured in Shergotty mineral separates by laser step heating. Martian atmosphere and ‘martian interior’ xenon are present, as is a spallation component. Martian atmospheric xenon is 5–10 times more concentrated in opaque minerals (magnetite, ilmenite, and pyrrhotite) and maskelynite than in pyroxenes, perhaps reflecting grain size variation. This is shown to be consistent with shock incorporation. A component consisting of solar xenon with a fission contribution, similar to components previously identified in martian meteorites and associated with the martian interior, is best defined in the pyroxene‐dominated separates. This component exhibits a consistent ¹²⁹Xe (¹²⁹Xe/¹³²Xe ?1.2) excess over solar/planetary (¹²⁹Xe/¹³²Xe ?1.04). We suggest that gas present in the melt, perhaps a mixture of interior xenon and martian atmosphere, was incorporated into the pyroxenes in Shergotty as the minerals crystallized.     

Continuous monitoring of nightside upper thermospheric mass densities in the martian southern hemisphere over 4 martian years using electron reflectometry

Details are presented of an improved technique to use atmospheric absorption of magnetically reflecting solar wind electrons to constrain neutral mass densities in the nightside martian upper thermosphere. The helical motion of electrons on converging magnetic field lines, through an extended neutral atmosphere, is modeled to enable prediction of loss cone pitch angle distributions measured by the Magnetometer/Electron Reflectometer (MAG/ER) experiment on Mars Global Surveyor at 400 km altitude. Over the small fraction of Mars' southern hemisphere (∼2.5%) where the permanent crustal magnetic fields are both open to the solar wind and sufficiently strong as to dominate the variable induced martian magnetotail field, spherical harmonic expansions of the crustal fields are used to prescribe the magnetic field along the electron's path, allowing least-squares fitting of measured loss cones, in order to solve for parameters describing the vertical neutral atmospheric mass density profile from 160 to 230 km. Results are presented of mass densities in the southern hemisphere at 2 a.m. LST at the mean altitude of greatest sensitivity, 180 km, continuously over four martian years. Seasonal variability in densities is largely explained by orbital and latitudinal changes in dayside insolation that impacts the nightside through the resulting thermospheric circulation. However, the physical processes behind repeatable rapid, late autumnal cooling at mid-latitudes and near-aphelion warming at equatorial latitudes is not fully clear. Southern winter polar warming is generally weak or nonexistent over several Mars years, in basic agreement with MGS and MRO accelerometer observations. The puzzling response of mid-latitude densities from 160° to 200° E to the 2001 global dust storm suggests unanticipated localized nightside upper thermospheric lateral and vertical circulation patterns may accompany such storms. The downturn of the 11-year cycle of solar EUV flux is likely responsible for lower aphelion densities in 2004 and 2006 (Mars years 27 and 28).     

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.     

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.     

Search for Phobos and Deimos gas/dust tori using in situ observations from Mars Global Surveyor MAG/ER

More than 490 elliptical aerobraking and science phasing orbits made by Mars Global Surveyor (MGS) in 1997 and 1998 provide unprecedented coverage of the solar wind in the vicinity of the orbits of the martian moons Phobos and Deimos. We have performed a comprehensive survey of magnetic field perturbations in the solar wind to search for possible signatures of solar wind interaction with dust or gas escaping from the moons. A total of 1246 solar wind disturbance events were identified and their distribution was examined relative to Phobos, the Phobos orbit, and the Deimos orbit. We find that the spatial distribution of solar wind perturbations does not increase near or downstream of Phobos, Phobos’ orbit, or Deimos’ orbit, which would have been expected if there is significant outgassing or dust escape from the martian moons. Of the 1246 magnetic field perturbation events found in the MGS data set, 11 events were found within 2000 km of the Phobos orbit, while three events were found within 2000 km of the Deimos orbit. These events were analyzed in detail and found to likely have other causes than outgassing/dust escape from the martian moons. Thus we conclude that the amount of gas/dust escaping the martian moons is not significant enough to induce detectable magnetic field perturbations in the solar wind. In essence we have not found any clear evidence in the MGS magnetic field data for outgassing or dust escape from the martian moons.     

On the convective mechanism for formation of the plasma sheet in the magnetospheric tail

Calculation of stationary distributions of the most important plasma parameters (particle energy, density, field-aligned and transversal pressure) is performed for a model magnetotail plasma sheet which is formed by convecting plasma mantle particles injected into the closed geomagnetic field line tubes. Computations have been done for two convection models: (i) a model of completely adiabatic particle motion with conservation of the first two invariants and (ii) a model with a strong pitch-angle diffusion which maintains isotropy. It is found that in both cases the heating and compression of the plasma are somewhat more effective than is necessary to account for the observed gradients of magnetic field in the magnetospheric tail. A leakage of accelerated particles through the dawn and dusk edges of the plasma sheet is proposed as a possible mechanism for maintenance of stationary convection in the magnetotail. The question of the dependence of the stationary magnetotail parameters on the solar wind state is discussed briefly.     

Features of the Martian Magnetic Field Structure

Based on the single-fluid MHD model of Mars space simulation, this paper has studied the magnetic field structure in the near-Mars space and investigated the influence of Martian crustal magnetic anomalies on the magnetic field structure. In the process of the solar wind interaction with Mars, the bow shock and magnetic pile-up region are produced. The interplanetary magnetic lines are curved and deformed while they are towed toward the two poles by the solar wind. The majority of magnetic lines bypass the two poles, then leave behind a ‘V-shaped’ structure in the magnetotail behind Mars. In the crust of Mars, the local magnetic anomalies have a noticeable influence on the magnetic field structure. The magnetic anomalies at different positions and in different intensities interact with the solar wind to form the mini-magnetospheres of different structures and morphologies, such as the towed mini-magnetosphere and the mini-magnetosphere with open magnetic lines. The local magnetic anomalies have changed the near-Mars magnetic field structure, and probably changed the plasma distribution as well.