Modeling of Venus, Mars, and Titan

Increased computer capacity has made it possible to model the global plasma and neutral dynamics near Venus, Mars and Saturn??s moon Titan. The plasma interactions at Venus, Mars, and Titan are similar because each possess a substantial atmosphere but lacks a global internally generated magnetic field. In this article three self-consistent plasma models are described: the magnetohydrodynamic (MHD) model, the hybrid model and the fully kinetic plasma model. Chamberlain and Monte Carlo models of the Martian exosphere are also described. In particular, we describe the pros and cons of each model approach. Results from simulations are presented to demonstrate the ability of the models to capture the known plasma and neutral dynamics near the three objects.     

Calculation of the Initial Magnetic Field for Mercury’s Magnetosphere Hybrid Model

Several types of numerical models are used to analyze the interactions of the solar wind flow with Mercury’s magnetosphere, including kinetic models that determine magnetic and electric fields based on the spatial distribution of charges and currents, magnetohydrodynamic models that describe plasma as a conductive liquid, and hybrid models that describe ions kinetically in collisionless mode and represent electrons as a massless neutralizing liquid. The structure of resulting solutions is determined not only by the chosen set of equations that govern the behavior of plasma, but also by the initial and boundary conditions; i.e., their effects are not limited to the amount of computational work required to achieve a quasi-stationary solution. In this work, we have proposed using the magnetic field computed by the paraboloid model of Mercury’s magnetosphere as the initial condition for subsequent hybrid modeling. The results of the model have been compared to measurements performed by the Messenger spacecraft during a single crossing of the magnetosheath and the magnetosphere. The selected orbit lies in the terminator plane, which allows us to observe two crossings of the bow shock and the magnetopause. In our calculations, we have defined the initial parameters of the global magnetospheric current systems in a way that allows us to minimize paraboloid magnetic field deviation along the trajectory of the Messenger from the experimental data. We have shown that the optimal initial field parameters include setting the penetration of a partial interplanetary magnetic field into the magnetosphere with a penetration coefficient of 0.2.     

The response of the Hermean magnetosphere to the interplanetary magnetic field

The interaction between the solar wind and Mercury is anticipated to be unique because of Mercury’s relatively weak intrinsic magnetic field and tenuous neutral exosphere. In this paper the role of the IMF in determining the structure of the Hermean magnetosphere is studied using a new self-consistent three-dimensional quasi-neutral hybrid model. A comparison between a pure northward and southward IMF shows that the general morphology of the magnetic field, the position and shape of the bow shock and the magnetopause as well as the density and velocity of the solar wind in the magnetosheath and in the magnetosphere are quite similar in these two IMF situations. A Parker spiral IMF case, instead, produces a magnetosphere with a substantial north–south asymmetric plasma and magnetic field configuration. In general, this study illustrates quantitatively the role of IMF when the solar wind interacts with a weakly magnetised planetary body.     

The plasma Environment of Mars

Space Science Reviews -     

RPC-ICA: The Ion Composition Analyzer of the Rosetta Plasma Consortium

The Ion Composition Analyzer (ICA) is part of the Rosetta Plasma Consortium (RPC). ICA is designed to measure the three-dimensional distribution function of positive ions in order to study the interaction between the solar wind and cometary particles. The instrument has a mass resolution high enough to resolve the major species such as protons, helium, oxygen, molecular ions, and heavy ions characteristic of dusty plasma regions. ICA consists of an electrostatic acceptance angle filter, an electrostatic energy filter, and a magnetic momentum filter. Particles are detected using large diameter (100 mm) microchannel plates and a two-dimensional anode system. ICA has its own processor for data reduction/compression and formatting. The energy range of the instrument is from 25 eV to 40 keV and an angular field-of-view of 360° × 90° is achieved through electrostatic deflection of incoming particles.     

IMF Direction Derived from Cycloid-Like Ion Distributions Observed by Mars Express

Although the Mars Express (MEX) does not carry a magnetometer, it is in principle possible to derive the interplanetary magnetic field (IMF) orientation from the three dimensional velocity distribution of pick-up ions measured by the Ion Mass Analyser (IMA) on board MEX because pick-up ions' orbits, in velocity phase space, are expected to gyrate around the IMF when the IMF is relatively uniform on a scale larger than the proton gyroradius. During bow shock outbound crossings, MEX often observed cycloid distributions (two dimensional partial ring distributions in velocity phase space) of protons in a narrow channel of the IMA detector (only one azimuth for many polar angles). We show two such examples. Three different methods are used to derive the IMF orientation from the observed cycloid distributions. One method is intuitive (intuitive method), while the others derive the minimum variance direction of the velocity vectors for the observed ring ions. These velocity vectors are selected either manually (manual method) or automatically using simple filters (automatic method). While the intuitive method and the manual method provide similar IMF orientations by which the observed cycloid distribution is well arranged into a partial circle (representing gyration) and constant parallel velocity, the automatic method failed to arrange the data to the degree of the manual method, yielding about a 30° offset in the estimated IMF direction. The uncertainty of the derived IMF orientation is strongly affected by the instrument resolution. The source population for these ring distributions is most likely newly ionized hydrogen atoms, which are picked up by the solar wind.     

Hybrid simulations of the O ion escape from Venus: Influence of the solar wind density and the IMF x component

As an initial effort to study the evolution of the Venus atmosphere, the influence of the solar wind density and the interplanetary magnetic field (IMF) x component (the x-axis points from Venus towards the Sun) on the O⁺ ion escape rate from Venus is investigated using a three-dimensional quasi-neutral hybrid (HYB-Venus) model. The HYB-Venus model is first applied to a case of the high-density (100 cm⁻³) solar wind interaction with Venus selected from the Pioneer Venus Orbiter observations to demonstrate its capability for the study. Two sets of simulations with a wide range of solar wind densities and different IMF x components are then performed. It is found that the O⁺ ion escape rate increases with increasing solar wind density. The O⁺ ion escape rate saturates when the solar wind density becomes high (above 100 cm⁻³). The results also suggest that the IMF x component enhances the O⁺ ion escape rate, given a fixed IMF component perpendicular to the x-axis. Finally, the results imply a higher ion loss rate for early-Venus, when solar conditions were dramatically different.     

Plasma Flow and Related Phenomena in Planetary Aeronomy

Understanding the processes involved in the interaction of solar system bodies with plasma flows is fundamental to the entire field of space physics. The features of the interaction can be very different, depending upon the properties of the incident plasma as well as the nature of the obstacle. The properties of the atmosphere/ionosphere associated with the obstacle are of particular importance into understanding the plasma interaction process, especially for non-magnetized obstacle. This paper discusses in detail the roles of the atmosphere and ionosphere systems of plasma interaction around Venus, Mars, comets and some particular satellites. The coupling between magnetosphere and ionosphere is also discussed for Earth and Giant planets.     

Energisation of O+ and O+ 2 Ions at Mars: An Analysis of a 3-D Quasi-Neutral Hybrid Model Simulation

We have studied the loss of O⁺ and O⁺ ₂ ions at Mars with a numerical model. In our quasi-neutral hybrid model ions (H⁺, He⁺⁺, O⁺, O⁺ ₂) are treated as particles while electrons form a massless charge-neutralising fluid. The employed model version does not include the Martian magnetic field resulting from the crustal magnetic anomalies. In this study we focus the Martian nightside where the ASPERA instrument on the Phobos-2 spacecraft and recently the ASPERA-3 instruments on the Mars Express spacecraft have measured the proprieties of escaping atomic and molecular ions, in particular O⁺ and O⁺ ₂ ions. We study the ion velocity distribution and how the escaping planetary ions are distributed in the tail. We also create similar types of energy-spectrograms from the simulation as were obtained from ASPERA-3 ion measurements. We found that the properties of the simulated escaping planetary ions have many qualitative and quantitative similarities with the observations made by ASPERA instruments. The general agreement with the observations suggest that acceleration of the planetary ions by the convective electric field associated with the flowing plasma is the key acceleration mechanism for the escaping ions observed at Mars.     

The Analyzer of Space Plasmas and Energetic Atoms (ASPERA-3) for the Mars Express Mission

The general scientific objective of the ASPERA-3 experiment is to study the solar wind – atmosphere interaction and to characterize the plasma and neutral gas environment with within the space near Mars through the use of energetic neutral atom (ENA) imaging and measuring local ion and electron plasma. The ASPERA-3 instrument comprises four sensors: two ENA sensors, one electron spectrometer, and one ion spectrometer. The Neutral Particle Imager (NPI) provides measurements of the integral ENA flux (0.1–60 keV) with no mass and energy resolution, but high angular resolution. The measurement principle is based on registering products (secondary ions, sputtered neutrals, reflected neutrals) of the ENA interaction with a graphite-coated surface. The Neutral Particle Detector (NPD) provides measurements of the ENA flux, resolving velocity (the hydrogen energy range is 0.1–10 keV) and mass (H and O) with a coarse angular resolution. The measurement principle is based on the surface reflection technique. The Electron Spectrometer (ELS) is a standard top-hat electrostatic analyzer in a very compact design which covers the energy range 0.01–20 keV. These three sensors are located on a scanning platform which provides scanning through 180^∘ of rotation. The instrument also contains an ion mass analyzer (IMA). Mechanically IMA is a separate unit connected by a cable to the ASPERA-3 main unit. IMA provides ion measurements in the energy range 0.01–36 keV/charge for the main ion components H⁺, He⁺⁺, He⁺, O⁺, and the group of molecular ions 20–80 amu/q. ASPERA-3 also includes its own DC/DC converters and digital processing unit (DPU).