Search

We model expected dynamo currents above, and resulting magnetic field profiles at, InSight's landing site on Mars, including for the first time the effect of electron-ion collisions. We calculate their diurnal and seasonal variability using inputs from global models of the Martian thermosphere, ionosphere, and magnetosphere. Modeled currents primarily depend on plasma densities and the strength of the neutral wind component perpendicular to the combined crustal and draped magnetic fields that thread the ionosphere. Negligible at night, currents are the strongest in the late morning and near solstices due to stronger winds and near perihelion due to both stronger winds and higher plasma densities from solar EUV photoionization. Resulting surface magnetic fields of tens of nanotesla and occasionally >100 nT may be expected at the InSight landing site. We expect currents and surface fields to vary significantly with changes in the draped magnetic field caused by Mars' dynamic solar wind environment. ©2019. American Geophysical Union. All Rights Reserved., R. J. L., C. M. F., and L. A. were funded by the NASA MAVEN mission prime contract. M. O. F. was funded by the NASA Insight mission prime contract. F. G. G. was funded by the Spanish National Research Council (CSIC) under intramural project CSIC 201450E022 and acknowledges financial support from the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award for the Instituto de Astrofisica de Andalucia-CSIC (SEV-2017-0709). F. F. was funded by the French Centre National De La Recherche Scientifique.

Multi-spacecraft missions after 2000 (Cluster II, THEMIS, Van Allen Probes, and MMS) have revolutionized our understanding of the causes, patterns and variability of a wide array of plasma phenomena in the terrestrial magnetospheric environment. ESCAPADE is a twin-spacecraft Mars mission concept that will similarly revolutionize our understanding of how solar wind momentum and energy flows throughout Mars’ magnetosphere to drive ion and sputtering escape, two processes which have helped shape Mars’ climate evolution over solar system history. ESCAPADE will measure magnetic field strength and topology, ion plasma distributions (separated into light and heavy masses), as well as suprathermal electron flows and thermal electron and ion densities, from coordinated elliptical, 200 km x ~7000 km orbits. ESCAPADE are small spacecraft ( ESCAPADE is due to complete its preliminary design review in August 2020, thereafter moving toward build, test, integration and launch two years later. We will report on science goals and objectives, mission design challenges, and provide a status update.

International audience; Many aspects of the Sun-Mars interaction have been investigated during solar transient events with measurements from multiple spacecrafts and also simulation efforts. Limited discussion has been paid to magnetic topology response to disturbed upstream conditions. The implications of topology changes include, but are not limited to, the pattern of energetic particle precipitation into the Martian atmosphere and the impact on cold ion escape during solar transient events as low-energy ion escape is dependent on magnetic topology. In this study, we investigate the magnetic topology response to the 2017 September interplanetary coronal mass ejection (ICME) event with measurements collected by the Mars Atmospheric and Volatile EvolutioN spacecraft. It is found that the interface between draped interplanetary magnetic field and closed field lines was moved from 800-1400 km in altitude during quiet conditions to 200-400 km after ICME arrived at Mars and then relaxed back to high altitudes again after the event. To gain insight into magnetic topology response on a global scale, we first validate magnetic topology from a time-dependent simulation with a single-fluid multispecies magnetohydrodynamic (MHD) model by comparing magnetic topology determined from Mars Atmospheric and Volatile EvolutioN data, which shows a good agreement. Then we present MHD predictions of global magnetic topology changes during this ICME event. In addition to a deeper interplanetary magnetic field penetration, MHD results suggest more open field lines in response to the ICME event.

Measurements provided by the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft are analyzed to investigate the Martian magnetotail configuration as a function of interplanetary magnetic field (IMF) BY. We find that the magnetotail lobes exhibit a ~45deg twist, either clockwise or counterclockwise from the ecliptic plane, up to a few Mars radii downstream. Moreover, the associated cross-tail current sheet is rotated away from the expected location for a Venus-like induced magnetotail based on nominal IMF draping. Data-model comparisons using magnetohydrodynamic simulations are in good agreement with the observed tail twist. Model field line tracings indicate that a majority of the twisted tail lobes are composed of open field lines, surrounded by draped IMF. We infer that dayside magnetic reconnection between the crustal fields and draped IMF creates these open fields and may be responsible for the twisted tail configuration, similar to what is observed at Earth.

We present results from a global Mars time-dependent MHD simulation under constant solar wind and solar radiation impact considering inherent magnetic field variations due to continuous planetary rotation. We calculate the 3-D shapes and locations of the bow shock (BS) and the induced magnetospheric boundary (IMB) and then examine their dynamic changes with time. We develop a physics-based, empirical algorithm to effectively summarize the multidimensional crustal field distribution. It is found that by organizing the model results using this new approach, the Mars crustal field shows a clear, significant influence on both the IMB and the BS. Specifically, quantitative relationships have been established between the field distribution and the mean boundary distances and the cross-section areas in the terminator plane for both of the boundaries. The model-predicted relationships are further verified by the observations from the NASA Mars Atmosphere and Volatile EvolutioN (MAVEN) mission. Our analysis shows that the boundaries are collectively affected by the global crustal field distribution, which, however, cannot be simply parameterized by a local parameter like the widely used subsolar longitude. Our calculations show that the variability of the intrinsic crustal field distribution in Mars-centered Solar Orbital itself may account for approx.60% of the variation in total atmospheric loss, when external drivers are static. It is found that the crustal field has not only a shielding effect for atmospheric loss but also an escape-fostering effect by positively affecting the transterminator ion flow cross-section area.

We report the first global, time-dependent simulation of the Mars upper atmospheric responses to a realistic solar flare event, an X8.2 eruption on 10 September 2017. The Mars Global Ionosphere-Thermosphere Model runs with realistically specified flare irradiance, giving results in reasonably good agreement with the Mars Atmosphere and Volatile EvolutioN spacecraft measurements. It is found that the ionized and neutral regimes of the upper atmosphere are significantly disturbed by the flare but react differently. The ionospheric electron density enhancement is concentrated below ~110-km altitude due to enhanced solar X-rays, closely following the time evolution of the flare. The neutral atmospheric perturbation increases with altitude and is important above ~150-km altitude, in association with atmospheric upwelling driven by solar extreme ultraviolet heating. It takes ~2.5 hr past the flare peak to reach the maximum disturbance and then additional ~10 hr to generally settle down to preflare levels.

In this Letter, we make use of sophisticated 3D numerical simulations to assess the extent of atmospheric ion and photochemical losses from Mars over time. We demonstrate that the atmospheric ion escape rates were significantly higher (by more than two orders of magnitude) in the past at $\sim 4$ Ga compared to the present-day value owing to the stronger solar wind and higher ultraviolet fluxes from the young Sun. We found that the photochemical loss of atomic hot oxygen dominates over the total ion loss at the current epoch whilst the atmospheric ion loss is likely much more important at ancient times. We briefly discuss the ensuing implications of high atmospheric ion escape rates in the context of ancient Mars, and exoplanets with similar atmospheric compositions around young solar-type stars and M-dwarfs., 6 pages, 4 figures, 1 table, accepted for publication in ApJ Letters

We study roles of the thermosphere and exosphere on the Martian ionospheric structure and ion escape rates in the process of the solar wind-Mars interaction. We employ a four-species multifluid MHD (MF-MHD) model to simulate the Martian ionosphere and magnetosphere. The $cold$ thermosphere background is taken from the Mars Global Ionosphere Thermosphere Model (M-GITM) and the $hot$ oxygen exosphere is adopted from the Mars exosphere Monte Carlo model - Adaptive Mesh Particle Simulator (AMPS). A total of four cases with the combination of 1D (globally averaged) and 3D thermospheres and exospheres are studied. The ion escape rates calculated by adopting 1D and 3D atmospheres are similar; however, the latter are required to adequately reproduce MAVEN ionospheric observations. In addition, our simulations show that the 3D hot oxygen corona plays an important role in preventing planetary molecular ions (O$_2^+$ and CO$_2^+$) escaping from Mars, mainly resulting from the mass loading of the high-altitude exospheric O$^+$ ions. The $cold$ thermospheric oxygen atom, however, is demonstrated to be the primary neutral source for O$^+$ ion escape during the relatively weak solar cycle 24., Comment: 21 pages, 10 figures, 4 tables, accepted for publication in Journal of Geophysical Research-Space Physics

International audience; Observations of the Mars upper atmosphere made from the Mars Atmosphere and Volatile Evolution (MAVEN) spacecraft have been used to determine the loss rates of gas from the upper atmosphere to space for a complete Mars year (16 Nov 2014 – 3 Oct 2016). Loss rates for H and O are sufficient to remove ∼2-3 kg/s to space. By itself, this loss would be significant over the history of the planet. In addition, loss rates would have been greater early in history due to the enhanced solar EUV and more-active Sun. Integrated loss, based on current processes whose escape rates in the past are adjusted according to expected solar evolution, would have been as much as 0.8 bar CO2 or 23 m global equivalent layer of H2O; these losses are likely to be lower limits due to the nature of the extrapolation of loss rates to the earliest times. Combined with the lack of surface or subsurface reservoirs for CO2 that could hold remnants of an early, thick atmosphere, these results suggest that loss of gas to space has been the dominant process responsible for changing the climate of Mars from an early, warmer environment to the cold, dry one that we see today.

We study the solar wind interaction with the Martian upper atmosphere during the 8 March 2015 interplanetary coronal mass ejection (ICME) by using a global multifluid MHD model. Comparison of the simulation results with observations from Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft shows good agreement. The total ion escape rate is increased by an order of magnitude, from 2.05 × 1024 s−1 (pre-ICME phase) to 2.25 × 1025 s−1 (ICME sheath phase), during this time period. Two major ion escape channels are illustrated: accelerated pickup ion loss through the dayside plume and ionospheric ion loss through the nightside plasma wake region. Interestingly, the tailward ion loss is significantly increased at the ejecta phase. Both bow shock and magnetic pileup boundary (BS and MPB) locations are decreased from (1.2RM, 1.57RM) at the pre-ICME phase to (1.16RM, 1.47RM), respectively, during the sheath phase along the dayside Mars-Sun line. Furthermore, both simulation and observational results indicate that there is no significant variation in the Martian ionosphere (at altitudes ≲ 200 km, i.e., the photochemical region) during this event.

An important mechanism in the generation of polar wind outflow is the ambipolar electric potential which assists ions in overcoming gravity and is a key mechanism for Terrestrial ionospheric escape. At Mars, open field lines are not confined to the poles, and outflow of ionospheric electrons is observed far into the tail. It has thus been hypothesized that a similar electric potential may be present at Mars, contributing to global ionospheric loss. However, no direct measurements of this potential have been made. In this pilot study, we examine photoelectron spectra measured by the Solar Wind Electron Analyzer instrument on the NASA Mars Atmosphere and Volatile EvolutioN (MAVEN) Mars Scout to put an initial upper bound on the total potential drop in the ionosphere of Mars of Φ♂≾⊥2V , with the possibility of a further ≾4.5 V potential drop above this in the magnetotail. If the total potential drop was close to the upper limit, then strong outflows of major ionospheric species (H+, O+, and O2+) would be expected. However, if most of the potential drop is confined below the spacecraft, as expected by current theory, then such a potential would not be sufficient on its own to accelerate O2+ to escape velocities, but would be sufficient for lighter ions. However, any potential would contribute to atmospheric loss through the enhancement of Jeans escape.

The nucleus of comet C/2013 A1 (Siding Spring) passed within 141,000 km of Mars on 19 October 2014. Thus, the cometary coma and the plasma it produces washed over Mars for several hours producing significant effects in the Martian magnetosphere and upper atmosphere. We present observations from Mars Atmosphere and Volatile EvolutioN's (MAVEN's) particles and field's instruments that show the Martian magnetosphere was severely distorted during the comet's passage. We note four specific major effects: (1) a variable induced magnetospheric boundary, (2) a strong rotation of the magnetic field as the comet approached, (3) severely distorted and disordered ionospheric magnetic fields during the comet's closest approach, and (4) unusually strong magnetosheath turbulence lasting hours after the comet left. We argue that the comet produced effects comparable to that of a large solar storm (in terms of incident energy) and that our results are therefore important for future studies of atmospheric escape, MAVEN's primary science objective.

A comprehensive study of the solar wind interaction with the Martian upper atmosphere is presented. Three global models: the 3-D Mars multifluid Block Adaptive Tree Solar-wind Roe Upwind Scheme MHD code (MF-MHD), the 3-D Mars Global Ionosphere Thermosphere Model (M-GITM), and the Mars exosphere Monte Carlo model Adaptive Mesh Particle Simulator (M-AMPS) were used in this study. These models are one-way coupled; i.e., the MF-MHD model uses the 3-D neutral inputs from M-GITM and the 3-D hot oxygen corona distribution from M-AMPS. By adopting this one-way coupling approach, the Martian upper atmosphere ion escape rates are investigated in detail with the combined variations of crustal field orientation, solar cycle, and Martian seasonal conditions. The calculated ion escape rates are compared with Mars Express observational data and show reasonable agreement. The variations in solar cycles and seasons can affect the ion loss by a factor of ∼3.3 and ∼1.3, respectively. The crustal magnetic field has a shielding effect to protect Mars from solar wind interaction, and this effect is the strongest for perihelion conditions, with the crustal field facing the Sun. Furthermore, the fraction of cold escaping heavy ionospheric molecular ions [( 2+ and/or 2+)/Total] are inversely proportional to the fraction of the escaping (ionospheric and corona) atomic ion [O+/Total], whereas 2+ and 2+ ion escape fractions show a positive linear correlation since both ion species are ionospheric ions that follow the same escaping path.

Without the shielding of a substantial intrinsic dipole magnetic field, the atmospheres of Mars and Venus are particularly susceptible to similar atmospheric ion energization and scavenging processes. However, each planet has different attributes and external conditions controlling its high altitude planetary ion spatial and energy distributions. This paper describes analogous test particle simulations in background MHD fields that allow us to compare the properties and fates, precipitation or escape, of the mainly O+ atmospheric pick-up ions at Mars and Venus. The goal is to illustrate how atmospheric and planetary scales affect the upper atmospheres and space environments of our terrestrial planet neighbors. The results show the expected convection electric field-related hemispheric asymmetries in both precipitation and escape, where the degree of asymmetry at each planet is determined by the planetary scale and local interplanetary field strength. At Venus, the kinetic treatment of O+ reveals a strong nightside source of precipitation while Mars' crustal fields complicate the simple asymmetry in ion precipitation and drive a dayside source of precipitation. The pickup O+ escape pattern at both Venus and Mars exhibits low energy tailward escape, but Mars exhibits a prominent, high energy ‘polar plume’ feature in the hemisphere of the upward convection electric field while the Venus ion wake shows only a modest poleward concentration. The overall escape is larger at Venus than Mars ( 2.1 × 10 25 and 4.3 × 10 24 at solar maximum, respectively), but the efficiency (likelihood) of O+ escaping is 2–3 times higher at Mars. The consequences of these comparisons for pickup ion related atmospheric energy deposition, loss rates, and detection on spacecraft including PVO, VEX, MEX and MAVEN are considered. In particular, both O+ precipitation and escape show electric field controlled asymmetries that grow with energy, while the O+ fluxes and energy spectra at selected spatial locations show characteristic signatures of the pickup related acceleration and precipitation.

International audience; With electron and magnetic field data obtained by the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft, we have identified closed magnetic field lines, with both foot points embedded in the dayside ionosphere, extending up to 6,200 km altitude into the Martian tail. This topology is deduced from photoelectrons produced in the dayside ionosphere being observed traveling both parallel and antiparallel to the magnetic field. At trapped-zone pitch angles (within a range centered on 90° where electrons magnetically reflect before interacting with the atmosphere), cases with either solar wind electrons or photoelectrons have been found, indicating different formation mechanisms for these closed loops. These large closed loops are present in MHD simulations. The case with field-aligned photoelectrons mixed with solar wind electrons having trapped-zone pitch angles is likely to be associated with reconnection, while the case with photoelectrons at all pitch angles is probably due to closed field lines being pulled tailward by the surrounding plasma flow. By utilizing an algorithm for distinguishing photoelectrons from solar wind electrons in pitch angle-resolved energy spectra, we systematically map the spatial distribution and occurrence rate of these closed magnetic loops over the region sampled by the MAVEN orbit. We find that the occurrence rate ranges from a few percent to a few tens of percent outside of the optical shadow and less than one percent within the shadow. These observations can be used to investigate the general magnetic topology in the tail, which is relevant to cold ion escape, reconnection, and flux ropes.

We address the important question of whether the newly discovered exoplanet, Proxima Centauri b (PCb), is capable of retaining an atmosphere over long periods of time. This is done by adapting a sophisticated multi-species MHD model originally developed for Venus and Mars, and computing the ion escape losses from PCb. The results suggest that the ion escape rates are about two orders of magnitude higher than the terrestrial planets of our Solar system if PCb is unmagnetized. In contrast, if the planet does have an intrinsic dipole magnetic field, the rates are lowered for certain values of the stellar wind dynamic pressure, but they are still higher than the observed values for our Solar system's terrestrial planets. These results must be interpreted with due caution, since most of the relevant parameters for PCb remain partly or wholly unknown., 7 pages, 2 figures, accepted for publication in ApJL

International audience; The Mars Atmosphere and Volatile Evolution mission has obtained comprehensive particle and magnetic field measurements. The Solar Wind Electron Analyzer provides electron energy-pitch angle distributions along the spacecraft trajectory that can be used to infer magnetic topology. This study presents pitch angle-resolved electron energy shape parameters that can distinguish photoelectrons from solar wind electrons, which we use to deduce the Martian magnetic topology and connectivity to the dayside ionosphere. Magnetic topology in the Mars environment is mapped in three dimensions for the first time. At low altitudes (

The crustal remnant field on Mars rotates with the planet at a period of 24 h 37 min, constantly varying the magnetic field configuration interacting with the solar wind. Until now, there has been no self-consistent modeling investigation on how this varying magnetic field affects the solar wind plasma interaction. Here we include the rotation of this localized crustal field in a multispecies single-fluid MHD model of Mars and simulate an entire day of solar wind interaction under normal solar wind conditions. The MHD model results are compared with Mars Global Surveyor (MGS) magnetic field observations and show very close agreement, especially for the field strength along almost all of the 12 orbits on the day simulated. Model results also show that the ion escape rates slowly vary with rotation, generally anticorrelating with the strength of subsolar magnetic crustal sources, with some time delay. In addition, it is found that in the intense crustal field regions, the densities of heavy ion components enhance significantly along the MGS orbit, implying strong influence of the crustal field on the ionospheric structures.

The 3-D multifluid Block Adaptive Tree Solar-wind Roe Upwind Scheme (BATS-R-US) MHD code (MF-MHD) is coupled with the 3-D Mars Thermospheric general circulation model (MTGCM). The ion escape rate from the Martian upper atmosphere is investigated by using a one-way coupling approach, i.e., the MF-MHD model incorporates the effects of 3-D neutral atmosphere profiles from the MTGCM model. The calculations are carried out for two cases with different solar cycle conditions. The calculated total ion escape flux (the sum of three major ionospheric species, O+, O2+, and CO2+) for solar cycle maximum conditions (6.6×1024 s−1) is about 2.6 times larger than that of solar cycle minimum conditions (2.5×1024 s−1). Our simulation results show good agreement with recent observations of 2–3×1024 s−1 (O+, O2+, and CO2+) measured near solar cycle minimum conditions by Mars Express. An extremely high solar wind condition is also simulated which may mimic the condition of coronal mass ejections or corotating interaction regions passing Mars. Simulation results show that it can lead to a significant value of the escape flux as large as 4.3×1025s−1.

Magnetic fields inconsistent with draped interplanetary magnetic fields and crustal fields have been observed on Mars. Considering the discovery of a global looping magnetic field around the Venusian magnetotail and the similarities in the solar wind interactions between Mars and Venus, we use Mars Atmosphere and Volatile Evolution observations to investigate the global looping field on Mars and its formation mechanism. It is found that a global looping field also exists on Mars; therefore, this type of global looping field is a common feature of unmagnetized planetary bodies with ionospheres, and therefore should also exist on Titan and near-Sun comets. The comparison of the looping fields on Mars and Venus shows that the looping field is stronger on Mars. Solar wind azimuthal flows around the magnetotail toward the − magnetotail polar region () are observed. We illustrate that the looping field can be formed by bending the draped field lines with these azimuthal flows, and that these azimuthal flows are associated with heavy ion plumes along the + direction that are expected to be stronger on Mars than Venus. The current system associated with the looping field and its possible connection with the nightside ionosphere formations and ion escapes on Mars and Venus are discussed.

[1] The Mars Test Particle model is used (with background parameters from a magnetohydrodynamic code) to simulate the transport of O+ ions in the near-Mars space environment to study the source processes responsible for ion escape. The MHD values at this altitude are used to inject an ionospheric outflow source of ions for the Mars Test Particle (MTP). The resulting loss distributions (in both real and velocity space) from this ionospheric source term are compared against those from high-altitude ionization mechanisms, in particular photoionization, charge exchange, and electron impact ionization, each of which has its own source regions, albeit overlapping. For the nominal MHD settings, this ionospheric outflow source contributes only 10% to the total O+ loss rate at solar maximum, predominantly via the central tail region. This percentage has very little dependence on the initial temperature, but a change in the initial ion density or bulk velocity directly alters this loss through the central tail. A density or bulk velocity increase of a factor of 10 makes the ionospheric outflow loss comparable in magnitude to the loss from the combined high-altitude sources. The spatial and velocity space distributions of escaping O+ are examined and compared for the various source terms to identify features specific to each ion source mechanism. For solar minimum conditions, the nominal MHD ionospheric outflow settings yield a 27% contribution to the total O+ loss rate, i.e., roughly equal to any one of the three high-altitude source terms with respect to escape.

[1] While the Mars upper atmosphere is continuously bombarded by charged particles of solar and planetary origins, the energy flux carried is often not sufficient to significantly affect the neutral atmosphere. However, we show that this is not the case during major space weather events. By applying two Mars global models—a Monte Carlo model for simulating pickup O+ precipitation at the exobase and a thermosphere-ionosphere model for assessing its global impact, we find that the thermospheric effects of reentering ions can change from negligible to very important when upstream solar wind conditions vary from normal to extreme. The atmospheric response under the most extreme conditions includes dramatic neutral temperature enhancement, significant neutral composition and wind changes, and increased importance of sputtering loss and possibly even thermal escape of heavy species.

We present a time-dependent MHD study of the controlling effects of the Mars crustal field on atmospheric escape. We calculate globally integrated planetary ion loss rates under quiet solar conditions considering the continuous rotation of crustal anomalies with the planet. It is found that the rotating crustal field plays an important role in controlling atmospheric escape. Significant time variation of ∼20% and ∼50% is observed during the entire rotation period for O+ and for O2+ and CO2+, respectively. The control is exerted mainly through two processes. First, the crustal magnetic pressure over the subsolar regime controls solar wind penetration and mass loading and therefore the escaping planetary ion source. There is a strong negative correlation between the magnetic pressure and ion loss, with a time lag of

International audience; Coupling between the lower and upper atmosphere, combined with loss of gas from the upper atmosphere to space, likely contributed to the thin, cold, dry atmosphere of modern Mars. To help understand ongoing ion loss to space, the Mars Atmosphere and Volatile Evolution (MAVEN) spacecraft made comprehensive measurements of the Mars upper atmosphere, ionosphere, and interactions with the Sun and solar wind during an interplanetary coronal mass ejection impact in March 2015. Responses include changes in the bow shock and magnetosheath, formation of widespread diffuse aurora, and enhancement of pick-up ions. Observations and models both show an enhancement in escape rate of ions to space during the event. Ion loss during solar events early in Mars history may have been a major contributor to the long-term evolution of the Mars atmosphere.

International audience; Two of the primary goals of the MAVEN mission are to determine how the rate of escape of Martian atmospheric gas to space at the current epoch depends upon solar influences and planetary parameters and to estimate the total mass of atmosphere lost to space over the history of the planet. Along with MAVEN’s suite of nine science instruments, a collection of complementary models of the neutral and plasma environments of Mars’ upper atmosphere and near-space environment are an indispensable part of the MAVEN toolkit, for three primary reasons. First, escaping neutrals will not be directly measured by MAVEN and so neutral escape rates must be derived, via models, from in situ measurements of plasma temperatures and neutral and plasma densities and by remote measurements of the extended exosphere. Second, although escaping ions will be directly measured, all MAVEN measurements are limited in spatial coverage, so global models are needed for intelligent interpolation over spherical surfaces to calculate global escape rates. Third, MAVEN measurements will lead to multidimensional parameterizations of global escape rates for a range of solar and planetary parameters, but further global models informed by MAVEN data will be required to extend these parameterizations to the more extreme conditions that likely prevailed in the early solar system, which is essential for determining total integrated atmospheric loss. We describe these modeling tools and the strategies for using them in concert with MAVEN measurements to greater constrain the history of atmospheric loss on Mars.

International audience; We present the results of an initial effort to statistically map the fluxes of planetary ions on a closed surface around Mars. Choosing a spherical shell ~1000 km above the planet, we map both outgoing and incoming ion fluxes (with energies >25 eV) over a 4 month period. The results show net escape of planetary ions behind Mars and strong fluxes of escaping ions from the northern hemisphere with respect to the solar wind convection electric field. Planetary ions also travel toward the planet, and return fluxes are particularly strong in the southern electric field hemisphere. We obtain a lower bound estimate for planetary ion escape of ~3 × 1024 s−1, accounting for the ~10% of ions that return toward the planet and assuming that the ~70% of the surface covered so far is representative of the regions not yet visited by Mars Atmosphere and Volatile EvolutioN (MAVEN).

International audience; With the upcoming MAVEN mission, the role of escape in the evolution of the Martian atmosphere is investigated in more detail. This work builds on our previous modeling of the atmospheric impact of the the pickup O+ sputtering effects for varioussolar wind parameters, solar EUV intensities, and the surface crustal field distributions. Relationships between the incident ions properties and the ejected hot neutral components, often referred to as atmospheric sputtering, are derived for application to proposed MAVEN ion spectrometer measurements of precipitating O+. We show how our simulation results can be used to constrain the sputtering effects under present conditions and to interpolate toward estimates of sputtering efficiencies occurring in earlier epochs. Present-day sputtering under typical circumstance is estimated to be weak, but possibly detectable as an exospheric enhancement. The ultimate goal of estimating the importance of atmospheric sputtering effects on the evolution ofthe Martian atmosphere can be better deduced by the combining MAVEN measurements with models, and the sputtering response relations derived here.

Atmospheric photoelectrons have been observed well above the ionosphere of Mars by the ASPERA-3 ELS instrument on Mars Express. To systematically interpret these observations, field lines from two global MHD simulations were analyzed for Connectivity to the dayside ionosphere (allowing photoelectron escape). It is found that there is a hollow cylinder behind the planet from 1–2 R M away from the Mars-Sun line that has a high probability of containing magnetic field lines with Connectivity to the dayside ionosphere. These results are in complete agreement with the ELS statistics. It is concluded that the high-altitude photoelectrons are the result of direct magnetic Connectivity to the dayside at the moment of the measurement, and no extra trapping or bouncing mechanisms are needed to explain the data.

This paper presents recent progress in understanding the Solar Wind-Mars interaction using a sophisticated global magnetohydrodynamic (MHD) model. Mars has only localized crustal magnetic fields, so the solar wind plasma flow interacts directly with the Mars atmosphere/ionosphere system. Such an interaction generates an induced current in the ionosphere, modifies the magnetic field environment around Mars, and more importantly, causes the erosion of the Mars atmosphere. The non-uniformly distributed crustal magnetic field also plays an important role in the interaction process, which is modulated by planetary rotation. Recent advances in computing power allow the inclusion of the continuous crustal field rotation in the simulation with a time-dependent MHD model. Model results have been validated with observations from previous and ongoing Mars missions. The validated time-dependent MHD model is useful in quantifying the variation of ion loss rates with planet rotation and the internal response time scale of the Martian ionosphere.

Without a global intrinsic magnetic field, oxygen ions in the Martian atmospheric corona can be picked up by the solar wind directly. The majority of the pickup ions escape, while some can precipitate into the atmosphere, producing sputtering of atmospheric constituents, which may play a significant role in the loss of neutral atmosphere. With its widely distributed crustal magnetic fields, Mars is unique. While it has been shown that Mars's crustal fields can alter the global distribution of escaping ions, little is known about the influence of the crustal fields on the spatial distribution and energy deposition of precipitating pickup ions. In this paper, the global distribution and energy spectrum of the precipitating pickup oxygen ions are calculated using the test particle method in the electric and magnetic fields set up by MHD simulation. Cases of different crustal field conditions are compared to show the influence of the crustal fields on the oxygen ion precipitation. We find, using a test particle code, that the crustal fields likely change the spatial distribution of the precipitation, resulting in irregularly distributed low-energy ion precipitation parches and wide high-energy ion precipitation belts on both the dayside and nightside. The crustal fields increase O(+) ion precipitation substantially, enhancing energy deposition especially on the nightside. Our results imply that atmospheric sputtering, and accompanying neutral escape, might happen globally and may be enhanced by a factor of almost 2 with the presence of the crustal magnetic fields.

[1] Our newly developed 3‐D, multifluid MHD model is used to study the interaction of the solar wind with Mars. This model is based on the BATS‐R‐US code, using a spherical grid and a radial resolution equal to 10 km in the ionospheric regions. We solve separate continuity, momentum, and energy equations for each ion fluid and run our model for both solar minimum and maximum conditions. We obtain asymmetric densities, velocities, and magnetic pileup in the plane containing both the direction of the solar wind and the convective electric field. These asymmetries are the result of the decoupling of the individual ions; therefore, our model is able to account for the respective dynamics of the ions and to show new physical processes that could not be observed by the single‐fluid model. Our results are consistent with the measured bow shock and magnetic pileup locations and with the Viking‐observed ion densities. We also compute the escape fluxes for both solar minimum and solar maximum conditions and compare them to the single‐fluid results and the observed values from Mars Express.

International audience; 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.

[1] A 3D, multi-species, non-ideal MHD numerical code was used to calculate the ion escape fluxes from Mars. The calculations were carried out for six cases with different nominal solar wind, solar cycle and crustal field orientation conditions and the total escape fluxes (the sum of the three major ionospheric species, O+, O2+, and CO2+) varied by about an order of magnitude from 2.7 × 1023 to 2.4 × 1024 sec−1. These results were compared to the recently measured Mars Express results of 3.2 × 1023 sec−1 (O+, O2+, and CO2+), which were obtained near solar cycle minimum conditions, indicating a good agreement between measured and calculated fluxes. We also calculated the escape flux for “extremely” high solar wind conditions which leads to a flux of 3 × 1025 sec−1.

[ 1] The Cassini spacecraft passed by Titan on 26 October 2004 (Ta flyby) and 13 December 2004 (Tb flyby). In both cases the Cassini Orbiter entered Titan's ionosphere and flew through Titan's dynamic wake region. In this paper, we present our simulation results of these two flybys using our three-dimensional multispecies MHD model. This model is a slightly updated version of the model outlined by Ma et al. (2004a); the main difference is the inclusion of magnetic diffusion into the equations. The calculations used the best available upstream plasma and magnetic field parameters obtained by the Cassini instrument complement. The calculated parameters, corresponding to the Cassini flybys near the closest approach, are compared with the relevant observed values. There is a reasonably good but clearly not perfect agreement between the measured and calculated values. Some of these differences are believed to be due to the uncertainties and time variability associated with the upstream parameters and some differences must definitely be the result of the uncertainties in the parameters selected for the model, as well as the limitations associated with the MHD approximations.

interaction of the solar wind with Mars. The three ions considered are H + ,O 2 , and O + , representing the solar wind and the two major ionospheric ion species, respectively. The calculations indicate that the presence of a hot oxygen corona does not, within the resolution and accuracy of the model, lead to any significant effect on the dayside bow shock and ionopause positions. Next the trans-terminator fluxes and escape fluxes down the tail were calculated neglecting the effects of the crustal magnetic field. The calculated flux values are consistent with the measured escape fluxes and the calculated limiting fluxes from the dayside ionosphere. Finally, a 60-order harmonic expansion model of the measured magnetic field was incorporated into the model. The crustal magnetic field did not cause major distortions in the bow shock but certainly had an important effect within the magnetosheath and on the apparent altitude of the ionopause. The model results also indicated the presence of ‘‘minimagnetocylinders,’’ consistent with the MGS observations. We also recalculated the trans-terminator and escape fluxes, for the nominal solar wind case, in the presence of the crustal magnetic field and found, as expected, that there is a decrease in the calculated escape flux; however, it is still reasonably close to the value estimated from the Phobos-2 observations. INDEX TERMS: 2780 Magnetospheric Physics: Solar wind interactions with unmagnetized bodies; 2459 Ionosphere: Planetary ionospheres (5435, 5729, 6026, 6027, 6028); 5440 Planetology: Solid Surface Planets: Magnetic fields and magnetism; 2728 Magnetospheric Physics: Magnetosheath; KEYWORDS: Mars, MHD, bow shock, escape flux, solar wind interaction, crustal magnetic field Citation: Ma, Y., A. F. Nagy, K. C. Hansen, D. L. DeZeeuw, T. I. Gombosi, and K. G. Powell, Three-dimensional multispecies MHD studies of the solar wind interaction with Mars in the presence of crustal fields, J. Geophys. Res., 107(A10), 1282, doi:10.1029/2002JA009293, 2002.