Monte Carlo computations of the escape of atomic nitrogen from Mars

Monte Carlo simulations were carried out to compute the escape flux of atomic nitrogen for the low and high solar activity martian thermospheres. The total escape of atomic nitrogen at low and high solar activities was found to be 3.03×¹⁰⁵ and , respectively. The escape flux of atomic nitrogen at low and high solar activities from photodissociation of N₂ was found to be 2.75×¹⁰⁵ and , respectively. The remainder of the contribution is from dissociative recombination, which is only important at high solar activity were it comprises about 25% of the total escape. The relative contributions to the total N escape flux from thermal motion of the background atmosphere, winds and co-rotation, and photoionization and subsequent solar wind pickup are also considered here. We find that the total predicted escape fluxes are observed to increase by 20 and 25% at low and high solar activities owing to thermal motion of the background atmosphere. At low and high solar activities, we find that the co-rotation and wind velocities combined translate to a maximum transferable energy of ∼0.0103 and 0.0181 eV, respectively, and that the total escape flux contribution from winds and co-rotation is negligible. Photoionization was found to be a minor process only impacting those source atoms produced with energies close to the escape energy, between 1.5 and 2 eV. The contributions to the total escape fluxes at low and high solar activities from photoionization and subsequent solar wind pickup are found to be about 8 and 13%, respectively.     

Comment on the papers “Production of hot nitrogen atoms in the martian thermosphere” by F. Bakalian and “Monte Carlo computations of the escape of atomic nitrogen from Mars” by F. Bakalian and R.E. Hartle

In recent articles published in Icarus, Bakalian [2006. Icarus 183, 69-78] discusses and computes the production rates of hot nitrogen atoms in the martian thermosphere due to N₂ photodissociation and N⁺₂ dissociative recombination, and Bakalian and Hartle [2006. Icarus 183, 55-68] use a Monte Carlo code to compute the escape rates of nitrogen atoms from Mars due to photodissociation of N₂, dissociative recombination of N⁺₂, and pickup ion escape due to photoionization of N atoms above the ionopause. Bakalian concludes that “photodissociation of N₂ is the dominant escape mechanism in the martian atmosphere.” We will show that this conclusion is not supportable. In addition, both papers contain scientific errors, misrepresentations, inaccurate referencing, lack of proper attribution, and they fail to place these investigations into the existing extensive body of work on this subject.     

Estimation of the escape of photoelectrons from Mars in 2004 liberated by the ionization of carbon dioxide and atomic oxygen

Photoelectron peaks in the atmosphere of Mars caused by the ionization of carbon dioxide and atomic oxygen by solar 30.4 nm photons have been observed by the Electron Spectrometer (ELS), a component of the Mars Express (MEx) Analyzer of Space Plasmas and Energetic Atoms (ASPERA-3) experiment. Ionization mostly occurs at the Mars exobase with the majority of the photoionized electron flux trapped in the remanent and induced magnetic field, with a portion of that flux escaping the planet down its tail. Since Mars is overall charge neutral, the number of electrons must be identical to the number of ion charges which escape the planet. An estimate of the fraction of the total number of escaping electrons is obtained for the year 2004, specifically those produced by the ionization of carbon dioxide and atomic oxygen by solar 30.4 nm photons. In achieving this process, an illustrative example pass is used to show how the electron spectrum is adjusted for the potential on the spacecraft; then the region of the electron spectrum which shows photoelectron peaks is integrated over energy, yielding a flux of 5.74 × 10⁶ electrons/(cm² s sr). This technique is then applied to a subset of 22 sample averaged spectra from the 2004 data (5 January 2004 through 25 January 2005), yielding an average result of 4.15 × 10⁶ electrons/(cm² s sr) for the 22 cases. The observation cone of 33.75° is used to integrate over solid angle (assuming the flux is constant), giving 4.39 × 10⁶ electrons/(cm² s). This average value was taken as representative of the full data interval. Frequency of occurrence statistics showing about a 6.2% occurrence rate for the 2004 data is applied to give an average escape flux from Mars of 2.72 × 10⁵ electrons/(cm² s) during 2004. By estimating the outflow area as 1.16 × 10¹⁸ cm² at X = −1.5 R_Mars the electron escape rate of 3.14 × 10²³ electrons/s is obtained. Thus about 9.92 × 10³⁰ electrons or 16.5 Mmole of electrons escaped Mars during 2004 due to the ionization of carbon dioxide and atomic oxygen by the He 30.4 nm line. Due to the caveats of the analysis, these derived escape rates should be considered lower limits on the total electron escape rate from Mars.     

Isotope fractionation in the photochemical escape of O from Mars

The production of energetic and escaping neutral O atoms at the current epoch in the martian thermosphere is thought to be dominated by the dissociative recombination process:

     

Ion escape at Mars: Comparison of a 3-D hybrid simulation with Mars Express IMA/ASPERA-3 measurements

We have analysed ion escape at Mars by comparing ASPERA-3/Mars Express ion measurements and a 3-D quasi-neutral hybrid model. As Mars Express does not have a magnetometer onboard, the analysed IMA data are from an orbit when the IMF clock angle was possible to determine from the magnetic field measurements of Mars Global Surveyor. We found that fast escaping planetary ions were observed at the place which, according to the 3-D model, is anticipated to contain accelerated heavy ions originating from the martian ionosphere. The direction of the interplanetary magnetic field was found to affect noticeably which regions can be magnetically connected to Mars Express and to the overall 3-D Mars-solar wind interaction.     

First ENA observations at Mars: ENA emissions from the martian upper atmosphere

The neutral particle detector (NPD) on board Mars Express has observed energetic neutral atoms (ENAs) from a broad region on the dayside of the martian upper atmosphere. We show one such example for which the observation was conducted at an altitude of 570 km, just above the induced magnetosphere boundary (IMB). The time of flight spectra of these ENAs show that they had energies of 0.2-2 keV/amu, with an average energy of ∼1.1 keV/amu. Both the spatial distribution and the energy of these ENAs are consistent with the backscattered ENAs, produced by an ENA albedo process. This is the first observation of backscattered ENAs from the martian upper atmosphere. The origin of these ENAs is considered to be the solar wind ENAs that are scattered back by collision processes in the martian upper atmosphere. The particle flux and energy flux of the backscattered ENAs are and , respectively.     

A comparison of global models for the solar wind interaction with Mars

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.     

Mass composition of the escaping plasma at Mars

Data from the Ion Mass Analyzer (IMA) sensor of the ASPERA-3 instrument suite on Mars Express have been analyzed to determine the mass composition of the escaping ion species at Mars. We have examined 77 different ion-beam events and we present the results in terms of flux ratios between the following ion species: CO⁺₂/O⁺ and O⁺₂/O⁺. The following ratios averaged over all events and energies were identified: CO⁺₂/O⁺ = 0.2 and O⁺₂/O⁺ = 0.9. The values measured are significantly higher, by a factor of 10 for O⁺₂/O⁺, than a contemporary modeled ratio for the maximum fluxes which the martian ionosphere can supply. The most abundant ion species was found to be O⁺, followed by O⁺₂ and CO⁺₂. We estimate the loss of CO⁺₂ to be by using the previous measurements of Phobos-2 in our calculations. The dependence of the ion ratios in relation to their energy ranges we studied, 0.3-3.0 keV, indicated that no clear correlation was found.     

Rebuttal to “Comment on the papers “Production of hot nitrogen atoms in the martian thermosphere” by F. Bakalian (2006), and “Monte Carlo computations of the escape of atomic nitrogen from Mars” by F. Bakalian and R.E. Hartle (2006)"

This response is to address some of the comments made by Dr. Fox and also to clarify some points in the Bakalian [2006. Icarus 183, 69-73] and Bakalian and Hartle [2006. Icarus 183, 55-68] papers. It is regrettable that some of the statements in the Bakalian and Bakalian and Hartle papers were misinterpreted.     

Production of hot nitrogen atoms in the martian thermosphere

The energy and altitude distributions for nitrogen atoms produced from photodissociation of N₂ and dissociative recombination of N⁺₂ have been computed for the low and high solar activity martian thermospheres. We find that the fraction of nitrogen atoms with E>Eesc at the exobase is ∼2.5% for photodissociation as compared to ∼22.6% for dissociative recombination. However, the peak magnitudes of the production rate profiles for dissociative recombination are a factor of ∼1000 and ∼200 times smaller than those of photodissociation at low and high solar activities, respectively. Thus, our findings suggest that photodissociation of N₂ is the dominant escape mechanism for atomic nitrogen for both the low and high solar activity martian thermospheres. The general consensus in the field prior to these research results was that dissociative recombination of N⁺₂ was the main escape mechanism. We find that it is the dissociation of N₂ along repulsive states that results in the production of very energetic nitrogen atoms.     

Model of photoelectron impact ionization within the high latitude ionosphere at Mars: Comparison of calculated and measured electron density

The production rate, ion density and electron density are calculated between longitudes 0° and 360° E due to incident radiation of wavelength range 1-102.57 nm in the dayside atmosphere of Mars. These calculations are made by using global analytical yield spectrum (AYS) model at solar zenith angle 80° between latitudes 50° and 70° N for spring equinox and medium solar activity condition. These conditions are appropriate for Mars Global Surveyor (MGS) Phase 2 aerobraking period during which both the accelerometer and the radio occultation data are used. The calculated results are compared with MGS radio occultation measurements carried out at different latitudes (64.7°-67.3° N) and longitudes (0°-360° E) in December 1998 between solar zenith angle 78° and 81°. This measurement shows primary and secondary ionization peaks, which are varying with longitudes. Our calculation suggests that first peak is produced by photoionization and photoelectron impact ionization processes due to absorption of solar EUV radiation (9-102.57 nm). The second peak is produced by photoelectron impact ionization of soft X-ray photon (1-9 nm). There is a good agreement between our calculation and measurement as far as the maximum and the minimum values of primary peak altitude/peak density of electrons are concerned. However, the calculated values of secondary peak density and peak altitude are higher than the measured values by a factor of 1.5-2.0 and 1.1, respectively. The secondary peak is brought into agreement with the measurement using low X-ray flux by a factor of 2 to 3 below 9 nm. The longitudinal distribution of calculated and measured peak density and peak altitude are fitted by least-square method with 0.95 confidence limits.     

The nitrogen cycle on Mars: Impact decomposition of near-surface nitrates as a source for a nitrogen steady state

Nitrogen isotopes appear to be escaping from Mars at approximately the primordial ratio ¹⁴N/¹⁵N ≈ 275 and to have an atmospheric nitrogen depletion time scale of about 800 Myr. For the standard model of a progressive decline of an initial inventory of atmospheric nitrogen, having no source of N, the agreement of the isotopic ratio of escaping N with primitive nitrogen would be coincidental. Here we propose a steady state model in which nitrates, produced early in Mars' history, are later decomposed by the current impact flux. The detection of near-surface nitrates can discriminate between the standard and the steady state models. Based on current estimates of N loss to space, we predict a quantity of nitrates equivalent to 60 ± 30 mbars for a steady state, or a global layer of about 3 m of pure NaNO₃.     

Ionospheric plasma acceleration at Mars: ASPERA-3 results

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

Atmospheric variations and meteorite production on Mars

Through a combination of aerobraking (drag deceleration) and ablation, meteoroids which enter planetary atmospheres may be slowed sufficiently to soft-land as meteorites. Results of an earlier study suggest that the current 6 mbar atmosphere of Mars is sufficient to aerobrake significant numbers of small (<10 kg) asteroidal-type meteoroids into survivable, low-velocity (<500 m s⁻¹) impacts with the planet's surface. Since rates of meteorite production depend upon the density of Mars's atmosphere, they must also change as the martian climate changes. However, to date, martian meteorite production has received relatively little attention in the literature Here we expand upon our previous work to study martian meteorite production rates and how they depend upon variations of the martian atmosphere, and to estimate the ranges of mass, velocity and entry-angle that produce meteorites. We find that even the current atmosphere of Mars is sufficient to soft-land significant fractions of incident stony and iron objects, and that these fractions increase dramatically for denser martian atmospheres. Therefore, like impact cratering, meteorite populations may preserve evidence of past martian climates.     

Impact processing of nitrogen on early Mars

An intense impact flux upon a planet having a CO₂ + N₂ atmosphere, such as Mars, provides energy to synthesize nitric oxide, NO, which is likely converted into nitrate minerals. The same impact flux can decompose nitrate minerals if present in the crust. We build a numerical model to study the effects of early impact processes on the evolution of nitrogen in a dominantly CO₂ atmosphere. We model the period of intense post-accretionary bombardment, the roughly 500 Myr period after crustal stabilization that locks in previously accreted volatiles. A best-guess, “fiducial” set of parameters is chosen, with a fixed “veneer” of post-accretionary impactors (δR=950 m thick), assumed to contain carbon at 1 wt% (fg=0.01), with a molar C/N ratio of 18, an initial atmospheric pressure of 1 bar (with CO₂/N₂ = 36), and a power law impactor mass distribution slope b=0.75. This model produces a nitrate reservoir RNO₃?0.5×¹⁰¹⁹ moles, equivalent to ∼30 mbars of N₂, during the intense impact phase. Starting with 1 bar, the atmosphere grows to 2.75 bars. Results of models with variations of parameter values show that RNO₃ responds sluggishly to changes in parameter values. To significantly limit the size of this reservoir, one is required to limit the initial total atmospheric pressure be less than about 0.5 bars, and the impactor volatile content fg to be less than 0.003. The value of fg substantially determines whether the atmosphere grows or not; when fg=0.01, the atmosphere gains about 1.7 bars, while for fg=0.003, the atmosphere gains less than 200 mbars, and for fg=0.001, it loses about 400 mbars. Impact erosion is a minor sink of N, constituting generally less than 10% of the total supply. The loss of impactor volatile plumes can take almost 50% of incoming N and C under fiducial parameters, when atmospheric pressures are low. This nitrogen does not significantly interact with Mars, and hence is not properly delivered. When the initial N is greater than the delivered N, most of the nitrogen ends up as nitrates; when delivered N is larger, most nitrogen ends up in the atmosphere. The reason for this dichotomy seems to be that initial nitrogen is present during the whole bombardment, while delivered N, on average, only experiences half the bombardment. The operating caveat here is that the above results are all conditioned on the assumption that impact processes dominate this period of Mars atmospheric evolution.     

THEMIS observations of Mars aerosol optical depth from 2002-2008

We use infrared images obtained by the Thermal Emission Imaging System (THEMIS) instrument on-board Mars Odyssey to retrieve the optical depth of dust and water ice aerosols over more than 3.5 martian years between February 2002 (MY 25, Ls=330°) and December 2008 (MY 29, Ls=183°). These data provide an important bridge between earlier TES observations and recent observations from Mars Express and Mars Reconnaissance Orbiter. An improvement to our earlier retrieval [Smith, M.D., Bandfield, J.L., Christensen, P.R., Richardson, M.I., 2003. J. Geophys. Res. 108, doi:10.1029/2003JE002114] to include atmospheric temperature information from THEMIS Band 10 observations leads to much improved retrievals during the largest dust storms. The new retrievals show moderate dust storm activity during Mars Years 26 and 27, although details of the strength and timing of dust storms is different from year to year. A planet-encircling dust storm event was observed during Mars Year 28 near Southern Hemisphere Summer solstice. A belt of low-latitude water ice clouds was observed during the aphelion season during each year, Mars Years 26 through 29. The optical depth of water ice clouds is somewhat higher in the THEMIS retrievals at ∼5:00 PM local time than in the TES retrievals at ∼2:00 PM, suggestive of possible local time variation of clouds.     

Wavenumber spectra of planetary-scale disturbances in the Mars atmosphere

Wavenumber spectra of the martian atmosphere covering zonal wavenumbers s=1-6 were obtained as a function of latitude and season for the first time from the temperatures measured by the Thermal Emission Spectrometer onboard the Mars Global Surveyor. The stationary component tends to peak at s=2, where the martian topography has large amplitude, and drops rapidly at higher wavenumbers. The transient component in the middle and high latitudes tends to peak at s=1, which is lower than the most unstable wavenumber based on linear theories, and exhibits spectral slopes much flatter than the stationary component. In the equatorial region, the spectra of the transient component are almost flat, indicating that the organization of large-scale structures is less efficient in this region. The spectral shapes are similar between the 0.5 and 2.2 hPa surfaces, except that the slopes are slightly steeper at 0.5 than at 2.2 hPa, probably due to selective vertical transmission at low wavenumbers. The seasonal variation is relatively large in the middle and high latitudes, where the maximum power occurs in winter and the minimum occurs in summer, with an exception that the transient component is maximum in spring in the southern hemisphere. Intensification of s=1 transient waves is observed around the period of the initiation of global dust storms.     

The depth of the convective boundary layer on Mars

We are using observations obtained with Mars Express to explore the structure and dynamics of the martian lower atmosphere. We consider a series of radio occultation experiments conducted in May-August 2004, when the season on Mars was midspring of the northern hemisphere. The measurements are widely distributed in latitude and longitude, but the local time remained within a narrow range, 17.0-17.2 h. Most of the atmospheric profiles retrieved from these data contain a distinct, well-mixed convective boundary layer (CBL). We have accurately determined the depth of the CBL and its spatial variations at fixed local time through analysis of these profiles. The CBL extends to a height of 3-10 km above the surface at the season and locations of these measurements. Its depth at fixed local time is clearly correlated with variations in surface elevation on planetary scales, with a weaker dependence on spatial variations in surface temperature. In general, the CBL is deep (8-10 km) where the surface elevation is high, as in Tharsis Montes and Syrtis Major, and shallow (4-6 km) where the surface elevation is low, as in Amazonis and Utopia. This variability results from the combined effects of conditions near the surface and in the atmosphere above the CBL. Convection arises from solar heating of the ground, and the impact of this heat source on thermal structure is largest where the surface pressure and atmospheric density are smallest, at high surface elevations. The vertical extent of the CBL is in turn constrained by the static stability of the overlying atmosphere. These results greatly reduce the long-standing uncertainty concerning the depth of the CBL.     

Heavy ion escape from Mars, influence from solar wind conditions and crustal magnetic fields

We have used more than 4 years of Mars Express ion data to estimate the escape of heavy ions ( and ) from Mars. To take the limited field of view of the instrument into account, the data has been binned into spatial bins and angular bins to create average distribution functions for different positions in the near Mars space. The net escape flux for the studied low solar activity period, between May 2007 and May 2011, is 2.0 ± 0.2 × 10²⁴ s⁻¹. The escape has been calculated independently for four different quadrants in the Y_MSO − Z_MSO plane, south, dusk, north and dawn. Escape is highest from the northern and dusk quadrants, 0.6 ± 0.1 × 10²⁴ s⁻¹, and smallest from the south and dawn quadrants, 0.4 ± 0.1 × 10²⁴ s⁻¹. The flux ratio of molecular ( and ) to O⁺ ions is 0.9 ± 0.1, averaged over all quadrants. The flux difference between the north and south quadrants is statistically significant, and is presumed to be due to the presence of significant crustal magnetic fields in the southern hemisphere, reducing the outflow. The difference between the dawn and dusk quadrants is likely due to the magnetic tension associated with the nominal Parker angle spiral, which should lead to higher average magnetic tension on the dusk side. The escape increases during periods of high solar wind flux and during times when co-rotating interaction regions (CIR) affect Mars. In the latter case the increase is a factor 2.4-2.9 as compared to average conditions.