Xiaohua Fang

The spatial distribution of planetary ion fluxes near Mars observed by MAVEN

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

Pickup oxygen ion velocity space and spatial distribution around Mars

[1] We report a newly created highly parallelized global test particle model for resolving the pickup oxygen ion distribution around Mars. The background magnetic and convection electric fields are calculated using a three-dimensional multispecies magnetohydrodynamic model, which includes the effect of the Martian crustal magnetic field. In addition to photo-ionization, charge exchange collisions and solar wind electron impact ionization are included for the pickup ion generation. The most novel feature of our model is that more than one billion test particles are launched in the simulation domain in total. This corresponds to a profound enhancement by at least 3 orders of magnitude in the total number, compared to all existing test particle models. This substantial improvement enables an unprecedented examination of the pickup ion flux distribution in velocity space, which is not achievable in previous simulation studies due to the insufficient statistics arising from the limited number of test particles. Using the velocity space distribution of pickup O + ions as a tool, the Mars-solar wind interaction can be investigated in a unique way. It is shown that the velocity space distribution is highly non-Maxwellian, exhibiting non-gyrotropic and non-symmetric distributions, including many beam-like features. In the tail region, pickup ions have a prominent outflowing component in the whole energy range. The energy examination of particles traveling across the tail region shows that the acceleration highly depends on the source region where the particles originate. The strong convection electric field in the magnetosheath region is favorable to the pickup ion acceleration.

Effects of crustal field rotation on the solar wind plasma interaction with Mars

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.

MHD model results of solar wind interaction with Mars and comparison with MAVEN plasma observations

The Mars Atmosphere and Volatile EvolutioN mission (MAVEN), launched on 18 November 2013, is now in its primary science phase, orbiting Mars with a 4.5 h period. In this study, we use a time-dependent MHD model to interpret plasma observations made by MAVEN particle and field instruments. Detailed comparisons between the model and the relevant plasma observations from MAVEN are presented for an entire Mars rotation under relatively quiet solar wind conditions. Through comparison along MAVEN orbits, we find that the time-dependent multispecies single-fluid MHD model is able to reproduce the main features of the plasma environment around Mars. Using the model results, we find that photoionization beyond the terminator is the dominant ion source as compared with day-night transport in maintaining the nightside ionosphere. Model results also show that both the time-varying solar wind conditions and the continuously rotating crustal field work together to control the ion escape variation with time.

Response of Mars O+ pickup ions to the 8 March 2015 ICME: Inferences from MAVEN data‐based models

We simulate and compare three phases of the Mars-solar wind interaction with the 8 March interplanetary coronal mass ejection (ICME) event using Mars Atmosphere and Volatile EvolutioN (MAVEN) mission observations in order to derive heavy ion precipitation and escape rates. The MAVEN observations provide the initial conditions for three steady state MHD model cases, which reproduce the observed features in the solar wind density, velocity, and magnetic field seen along the MAVEN orbit. Applying the MHD results to a kinetic test particle model, we simulate global precipitation and escape maps of O+ during the (1) pre-ICME phase, (2) sheath phase, and (3) ejecta phase. We find that the Case 1 had the lowest precipitation and escape rates of 9.5 × 1025 and 4.1 × 1025 s−1, Case 2 had the highest rates of 9.5 × 1025 and 4.1 × 1025 s−1, and Case 3 had rates of 3.2 × 1025 and 1.3 × 1025 s−1, respectively. Additionally, Case 2 produced a high-energy escaping plume >10 keV, which mirrored corresponding STATIC observations.

Modeling of the O+ pickup ion sputtering efficiency dependence on solar wind conditions for the Martian atmosphere

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 Mars Atmosphere and Volatile Evolution (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 2x10(24)s(-1) at nominal solar wind condition and can be enhanced about 50 times when both the interplanetary magnetic field (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 nightside 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.

Test particle comparison of heavy atomic and molecular ion distributions at Mars

This study uses the Mars Test Particle simulation to create virtual detections of O+, O2+, and CO2+ in an orbital configuration in the Mars space environment. These atomic and molecular planetary pickup ions are formed when the solar wind directly interacts with the neutral atmosphere, causing the ions to be accelerated by the background convective electric field. The subsequent ion escape is the subject of great interest, specifically with respect to which species dominates ion loss from Mars. O+ is found to be the dominant escaping ion because of the large sources of transported ions in the low-energy ( 1 keV) range. O2+ and CO2+ are observed at these energy ranges but with much lower fluxes and are generally only found in the tail between 10 eV and 1 keV. Using individual particle traces, we reveal the origin and trajectories of the low-energy downtail O+ populations and high-energy polar O+ populations that contribute to the total escape. Comparing them against O2+ and CO2+ reveals that the extended hot oxygen corona contributes to source regions of high- and low-energy escaping ions. Additionally, we present results for solar minimum and maximum conditions with respect to ion fluxes and energies in order to robustly describe the physical processes controlling planetary ion distributions and atmospheric escape.

Martian low‐altitude magnetic topology deduced from MAVEN/SWEA observations

The Mars Atmosphere and Volatile Evolution (MAVEN) mission has obtained comprehensive particle and magnetic field measurements. The Solar Wind Electron Analyzer (SWEA) 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 (< 400 km) in sunlight, the northern hemisphere is found to be dominated by closed field lines (both ends intersecting the collisional atmosphere), with more day-night connections through cross-terminator closed field lines than in the south. Although draped field lines with ~100-km-amplitude vertical fluctuations that intersect the electron exobase (~160-220 km) in two locations could appear to be closed at the spacecraft, a more likely explanation is provided by crustal magnetic fields, which naturally have the required geometry. Around 30% of the time, we observe open field lines from 200-400 km, which implies three distinct topological layers over the northern hemisphere: closed field lines below 200 km, open field lines with footpoints at lower latitudes that pass over the northern hemisphere from 200-400 km, and draped IMF above 400 km. This study also identifies open field lines with one end attached to the dayside ionosphere and the other end connected with the solar wind, providing a path for ion outflow.

Statistical studies on Mars atmospheric sputtering by precipitating pickup O+: Preparation for the MAVEN mission

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.

Ionization due to electron and proton precipitation during the August 2011 storm

The parameterizations of monoenergetic particle impact ionization in Fang et al. (2010) (Fang2010) and Fang et al. (2013) (Fang2013) are applied to the complex energy spectra measured by DMSP F16 satellite to calculate the ionization rates from electron and ion precipitations for a Northern Hemisphere pass from 0030 UT to 0106 UT on 6 August 2011. Clear enhancement of electron flux is found in the polar cap. The mean electron energy in the polar cap is mostly above 100 eV, while the mean energy in the auroral zone is typically above 1 keV. At the same time, F16 captures a strong Poynting flux enhancement in the polar cap, which is comparable to those in the auroral zone. The particle impact ionization rates using Fang2010 and Fang2013 parameterizations show clear enhancement at F region altitudes mainly due to the low-energy precipitating electrons, peaking probably in the cusp but also showing enhanced levels throughout most of the polar cap region. The general circulation models (GCMs), National Center for Atmospheric Research Thermosphere-Ionosphere-Electrodynamics General Circulation Model, and Global Ionosphere-Thermosphere Model, using their default empirical formulations of particle impact ionization, do not capture the observed features shown in the total particle ionization rate applying the Fang2010 and Fang2013 parameterizations to DMSP measurements. The difference between GCM simulations and Fang2010 and Fang2013 applied to DMSP data is due to the difference of both the inputs to the models and the parameterization of the ionization rates.