Search

The Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight) mission was launched on 5 May 2018 and successfully landed at Elysium Planitia (4.5oN, 135.9oE)on Mars on 26 November 2018. The InSight Lander carries a magnetometer to measure disturbances from the Martian ionosphere. In order to understand the daily variations in the magnet field measurements on Martian surface, in this study, we use the time-dependent MHD model to study how plasma conditions vary with local time above insight landing site using solar wind condition from MAVEN observation. Significant diurnal variations can be seen in all plasma quantities due to solar wind interactions and planetary rotation. The induced magnetic field is mainly in the same direction as the upstream IMF. However, it seems that the variations seen by the Insight magnetometer cannot be only due to the interaction of the solar wind. We also add a neutral wind effect in our simulations to further investigate possible causes of surface field changes.

Mars Atmosphere and Volatile EvolutioN (MAVEN) mission observations show clear evidence of the occurrence of the magnetic reconnection process in the Martian plasma tail. In this study, we use soph ...

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.

International audience; Global models of the interaction of the solar wind with the Martian upper atmosphere have proved to be valuable tools for investigating both the escape to space of the Martian atmosphere and the physical processes controlling this complex interaction. The many models currently in use employ different physical assumptions, but it can be difficult to directly compare the effectiveness of the models since they are rarely run for the same input conditions. Here we present the results of a model comparison activity, where five global models (single‐fluid MHD, multi‐fluid MHD, multi‐fluid electron pressure MHD, and two hybrid models) were run for identical conditions corresponding to a single orbit of observations from the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft.We find that low altitude ion densities are very similar across all models, and are comparable to MAVEN ion density measurements from periapsis. Plasma boundaries appear generally symmetric in all models and vary only slightly in extent. Despite these similarities there are clear morphological differences in ion behavior in other regions such as the tail and southern hemisphere. These differences are observable in ion escape loss maps, and are necessary to understand in order to accurately use models in aiding our understanding of the martian plasma environment.

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

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.

We develop a parametric fit to the results of a detailed magnetohydrodynamic (MHD) study of the response of ion escape rates (O+, O 2 + and CO 2 + ) to strongly varied solar forcing factors, as a way to efficiently extend the MHD results to different conditions. We then use this to develop a second, evolutionary model of solar forced ion escape. We treat the escape fluxes of ion species at Mars as proportional to the product of power laws of four factors – that of the EUV flux Reuv, the solar wind particle density Rρ, its velocity (squared) R v 2 , and the interplanetary magnetic field pressure R B 2 , where forcing factors are expressed in units of the current epoch-averaged values. Our parametric model is: ϕ ( i ) = ϕ 0 ( i ) R euv α ( i ) R ρ β ( i ) R v 2 γ ( i ) R B 2 δ ( i ) , where ϕ(i) is the escape flux of ion i. We base our study on the results of just six provided MHD model runs employing large forcing factor variations, and thus construct a successful, first-order parametric model of the MHD program. We perform a five-dimensional least squares fit of this power law model to the MHD results to derive the flux normalizations and the indices of the solar forcing factors. For O+, we obtain the values, 1.73 × 1024 s−1, 0.782, 0.251, 0.382, and 0.214, for ϕ0, α, β, γ, and δ, respectively. For O 2 + , the corresponding values are 1.68 × 1024 s−1, −0.393, 0.798, 0.967, and 0.533. For CO 2 + , they are 8.66 × 1022 s−1, −0.427, 1.083, 1.214, and 0.690. The fit reproduces the MHD results to an average error of about 5%, suggesting that the power laws are broadly representative of the MHD model results. Our analysis of the MHD model shows that by itself an increase in REUV enhances O+ loss, but suppresses the escape of O 2 + and CO 2 + , whereas increases in solar wind (i.e., in R ρ , R v 2 , and R B 2 , with Reuv constant) favors the escape of heavier ions more than light ions. The ratios of escaping ions detectable at Mars today can be predicted by this parametric fit as a function of the solar forcing factors. We also use the parametric model to compute escape rates over martian history. This second parametric model expresses ion escape functions of one variable (per ion), ϕ(i) = ϕ0(i)(t/t0)−ξ(i). The ξ(i) are linear combinations of the epoch-averaged ion escape sensitivities, which are seen to increase with ion mass. We integrate the CO 2 + and oxygen ion escape rates over time, and find that in the last 3.85 Gyr, Mars would have lost about 25 - 0.19 + 85 mbars of CO 2 + , and 0.64 - 0.34 + 0.62 m of water (from O+ and O 2 + ) from ion escape.

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.

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

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.

International audience; Sputtering of the Martian atmosphere by O+ pickup ions has been proposed as a potentially important process in the early evolution of the Martian atmosphere. In preparation for the MAVEN mission, we performed a study using a Monte Carlo model coupled to a molecular dynamic calculation to investigate the cascade sputtering effects in the region of the Martian exobase. Pickup ion fluxes based on test particle simulations in an MHD model for three different solar wind conditions are used to examine the local and global sputtering efficiencies. The resultant sputtering escape rate is 2 × 1024 s− 1 at nominal solar wind condition, and can be enhanced about 50 times when both the IMF strength and the solar wind pressure increase. It is found that when the IMF strength becomes stronger, both the pickup ion precipitation energies and the resultant sputtering efficiencies increase. The related escape flux, hot component, and atmospheric energy deposition deduced from the MAVEN measurements may reveal clues about the prominent enhanced sputtering effects. Significant hemispheric asymmetries can be observed related to the solar wind electric fields. The shielding by the crustal fields and the recycling onto the night-side due to different magnetic field draping features can also lead to regional variations of sputtering efficiencies. The results suggest that disturbed or enhanced solar wind conditions provide the best prospects for detecting sputtering effects for MAVEN mission.

[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.

The magnetosphere contains a significant amount of ionospheric O{+}, particularly during geomagnetically active times. The presence of ionospheric plasma in the magnetosphere has a notable impact on magnetospheric composition and processes. We present a new multifluid MHD version of the BATS-R-US model of the magnetosphere to track the fate and consequences of ionospheric outflow. The multi-fluid MHD equations are presented as are the novel techniques for overcoming the formidable challenges associated with solving them. Our new model is then applied to the May 4, 1998 and March 31, 2001 geomagnetic storms. The results are juxtaposed with traditional single- fluid MHD and multispecies MHD simulations from a previous study, thereby allowing us to assess the benefits of using a more complex model with additional physics. We find that our multi-fluid MHD model (with outflow) gives comparable results to the multi-species MHD model (with outflow), including a more strongly negative Dst, reduced CPCP, and a drastically improved magnetic field at geosynchronous orbit, as compared to single-fluid MHD with no outflow. Significant differences in composition and magnetic field are found between the multi-species and multi-fluid approach further away from the Earth. We further demonstrate the ability to explore pressure and bulk velocity differences between H{+} and O(+}, which is not possible when utilizing the other techniques considered.

[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.

[1] The global features of the interaction between Saturn’s magnetospheric plasma flow and Titan’s atmosphere/ ionosphere are simulated by using a 3-D, multi-species, high spatial resolution, global MHD model. Our model uses a spherical grid structure leading to very good (36 km) altitude resolution in the ionospheric region of Titan. The model also provides good resolution and meaningful results in the upstream and wake regions. Titan’s atmosphere and ionosphere are approximated by 10 neutral and 7 ion species. Calculations, which make predictions for the anticipated results to be obtained by the Cassini spacecraft during its first close flyby (TA) of Titan, are presented. INDEX TERMS: 2427 Ionosphere: Ionosphere/ atmosphere interactions (0335); 2732 Magnetospheric Physics: Magnetosphere interactions with satellites and rings; 2753 Magnetospheric Physics: Numerical modeling. Citation: Ma, Y.-J., A. F. Nagy, T. E. Cravens, I. V. Sokolov, J. Clark, and K. C. Hansen (2004), 3-D global MHD model prediction for the first close flyby of Titan by Cassini, Geophys. Res. Lett., 31, L22803

[1] We present the results of model calculations, using our new, four-species, spherical MHD model. Our results are compared with the relevant and limited available data. The resulting comparisons help us to increase our understanding of the interaction processes between the solar wind and the Martian atmosphere/ionosphere. This new model with a spherical grid structure allowed us to use small (� 10 km) radial grid spacing in the ionospheric region. We found that the calculated bow shock positions agree reasonably well with the observed values. The calculated results vary with interplanetary magnetic field orientation, solar cycle conditions, and subsolar location. We found that our calculated ion densities, with parameters corresponding to solar cycle minimum conditions, reproduced the Viking 1 observed ion densities well. The calculated solar cycle maximum densities, above � 140 km, are also consistent with the appropriate Mars Global Surveyor radio occultation electron densities. Both the calculated solar cycle maximum and solar cycle minimum total transterminator and escape fluxes are significantly smaller than our previously published values. This decrease is due to the improved temperature values used for the recombination rates in this new model, which in turn results in lower ion densities and lower fluxes. INDEX TERMS: 2780 Magnetospheric Physics: Solar wind interactions with unmagnetized bodies; 6026 Planetology: Comets and Small Bodies: Ionospheres— composition and chemistry; 6028 Planetology: Comets and Small Bodies: Ionospheres—structure and dynamics; 2728 Magnetospheric Physics: Magnetosheath; KEYWORDS: Mars, MHD, bow shock, ionosphere, solar wind interaction

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.