First ENA observations at Mars: Charge exchange ENAs produced in the magnetosheath

Measurements of energetic neutral atoms (ENA) generated in the magnetosheath at Mars are reported. These ENAs are the result of charge exchange collisions between solar wind protons and neutral oxygen and hydrogen in the exosphere of Mars. The peak of the observed ENA flux is . For the case studied here, i.e., the passage of Mars Express through the martian magnetosheath around 20:15 UT on 3 May 2004, the measurements agree with an analytical model of the ENA production at the planet. It is possible to find parameter values in the model such that the observed peak in the ENA count rate during the spacecraft passage through the magnetosheath is reproduced.     

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.     

Interferometric millimeter observations of water vapor on Mars and comparison with Mars Express measurements

We present interferometric mapping of the 225.9-GHz HDO and 203.4-GHz lines on Mars obtained with the IRAM Plateau de Bure facility (PdBI). The observations were performed during martian year 28 (MY28), at L_s=320.3° for the HDO line, and at L_s=324.3° for the line. The HDO line is detected at the eastern (morning) and western (evening) limbs of the northern hemisphere, corresponding to a water column density in the range 3-6 pr.-μm. The line is not detected, which is compatible with the column densities derived from the HDO line. Quasi-simultaneous far infrared measurements obtained by the Planetary Fourier Spectrometer (PFS) onboard the Mars Express spacecraft confirm our PdBI results, yielding a 5±1 pr.-μm meridionally constant water column abundance.Such a low water abundance during the southern mid-autumn of MY28 does not correspond to the standard martian climatology as observed during the previous years. It was however already retrieved from near-infrared observations performed by the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) onboard the Mars Reconnaissance Orbiter spacecraft [Smith, M.D., Wolff, M.J., Clancy, R.T., Murchie, S.L. 2009. CRISM observations of water vapor and carbon monoxide. J. Geophys. Res. 114, doi: 10.1029/2008JE003288]. Our observations thus confirm that the planet-encircling dust storm that occurred during MY28 significantly affected the martian water cycle. Our observations also demonstrate the usefulness of interferometric submillimeter observations to survey the martian water cycle from ground-based facilities.     

Rosetta-Alice observations of exospheric hydrogen and oxygen on Mars

The European Space Agency’s Rosetta spacecraft, en route to a 2014 encounter with comet 67P/Churyumov-Gerasimenko, made a gravity assist swing-by of Mars on 25 February 2007, closest approach being at 01:54 UT. The Alice instrument on board Rosetta, a lightweight far-ultraviolet imaging spectrograph optimized for in situ cometary spectroscopy in the 750-2000 Å spectral band, was used to study the daytime Mars upper atmosphere including emissions from exospheric hydrogen and oxygen. Offset pointing, obtained five hours before closest approach, enabled us to detect and map the H i Lyman-α and Lyman-β emissions from exospheric hydrogen out beyond 30,000 km from the planet’s center. These data are fit with a Chamberlain exospheric model from which we derive the hydrogen density at the 200 km exobase and the H escape flux. The results are comparable to those found from the Ultraviolet Spectrometer experiment on the Mariner 6 and 7 fly-bys of Mars in 1969. Atomic oxygen emission at 1304 Å is detected at altitudes of 400-1000 km above the limb during limb scans shortly after closest approach. However, the derived oxygen scale height is not consistent with recent models of oxygen escape based on the production of suprathermal oxygen atoms by the dissociative recombination of .     

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:

     

Oxygen ion escape at Mars in a hybrid model: High energy and low energy ions

We have studied the escape and energization of several O⁺ populations and an population at Mars by using a hybrid model. The quasi-neutral hybrid model, HYB-Mars model, included five oxygen ion populations making it possible to distinguish photoions from oxygen ions originating from charge exchange processes and from the ionosphere.We have identified two high-energy ion components and one low-energy ion component of oxygen. They have different spatial and energy distributions near Mars. The two high-energy oxygen ion components, consisting of a high-energy “beam” and a high-energy “halo”, have different origins. (1) The high-energy (>∼100 eV) “beam” of O⁺ and ions are originating from the ionosphere. These ions form a highly asymmetric spatial distribution of escaping oxygen ions with respect to the direction of the convective electric field in the solar wind. (2) The high-energy (>∼100 eV) “halo” component contains O⁺ ions which are formed from the oxygen neutral exosphere by extreme ultraviolet radiation (EUV) and by charge exchange processes. These energetic halo ions can be found all around Mars. (3) The low energy O⁺ and ions (<∼100 eV) form a relatively symmetric spatial distribution around the Mars-Sun line. They originate from the ionosphere and from charge exchange processes between protons and exospheric oxygen atoms.The existence of the low- and the high-energy oxygen components is in agreement with recent in situ plasma measurements made by the ASPERA-3 instrument on the Mars Express mission. The analysis of the escaping oxygen ions suggests that the global energization of escaping planetary ions in the martian tail is controlled by the convective electric field.     

Titan's hydrodynamically escaping atmosphere

The upper atmosphere of Titan is currently losing mass at a rate , by hydrodynamic escape as a high density, slow outward expansion driven principally by solar UV heating by CH₄ absorption. The hydrodynamic mass loss is essentially CH₄ and H₂ escape. Their combined escape rates are restricted by power limitations from attaining their limiting rates (and limiting fluxes). Hence they must exhibit gravitational diffusive separation in the upper atmosphere with increasing mixing ratios to eventually become major constituents in the exosphere. A theoretical model with solar EUV heating by N₂ absorption balanced by HCN rotational line cooling in the upper thermosphere yields densities and temperatures consistent with the Huygens Atmospheric Science Investigation (HASI) data [Fulchignoni, M., and 42 colleagues, 2005. Nature 438, 785-791], with a peak temperature of ∼185-190 K between 3500-3550 km. This model implies hydrodynamic escape rates of and , or some other combination with a higher H₂ escape flux, much closer to its limiting value, at the expense of a slightly lower CH₄ escape rate. Nonthermal escape processes are not required to account for the loss rates of CH₄ and H₂, inferred by the Cassini Ion Neutral Mass Spectrometer (INMS) measurements [Yelle, R.V., Borggren, N., de la Haye, V., Kasprzak, W.T., Niemann, H.B., Müller-Wodarg, I., Waite Jr., J.H., 2006. Icarus 182, 567-576].     

Evaluating the meaning of “layer” in the martian north polar layered deposits and the impact on the climate connection

Using data from the High Resolution Imaging Science Experiment (HiRISE) aboard the Mars Reconnaissance Orbiter, we reassess the methods by which layers within the north polar layered deposits (NPLD) can be delineated and their thicknesses measured. Apparent brightness and morphology alone are insufficient for this task; high resolution topographic data are necessary. From these analyses, we find that the visible appearance of layers depends to a large degree on the distribution of younger, mantling deposits (which in turn is partially influenced by inherent layer properties) and on the shape and location of the particular outcrop. This younger mantle partially obscures layer morphology and brightness and is likely a cause of the gradational contacts between individual layers at this scale. High resolution images reveal that there are several layers similar in appearance to the well-known marker bed discovered by Malin, M., Edgett, K., 2001. J. Geophys. Res. 106, 23429-23570. The morphology, thicknesses , and separation distances of these marker beds, as gleaned from a high resolution stereo digital elevation model, lend insight into the connection between stratigraphy and climate.     

Arecibo radar observations of Phobos and Deimos

In November 2005, we observed the moons of Mars using the Arecibo 2380-MHz (13-cm) radar, obtaining a result for the OC radar albedo of Phobos (0.056±0.014) consistent with its previously reported radar albedo and implying an upper bound on its near-surface bulk density of . We detected Deimos by radar for the first time, finding its OC radar albedo to be 0.021±0.006, implying an upper bound on its near-surface density of , consistent with a high-porosity regolith. We briefly discuss reasons for these low radar albedos, Deimos' being possibly the lowest of any Solar System body yet observed by radar.     

Measurement of the meteoroid flux at Mars

In the fall of 2005, a dedicated meteor observing campaign was carried out by the Panoramic Camera (Pancam) onboard the Mars Exploration Rover (MER) Spirit to determine the viability of using MER cameras as meteor detectors and to obtain the first experimental estimate of the meteoroid flux at Mars. Our observing targets included both the sporadic meteoroid background and two predicted martian meteor showers: one associated with 1P/Halley and a potential stream associated with 2001/R1 LONEOS. A total of 353 images covering 2.7 h of net exposure time were analyzed with no conclusive meteor detections. From these data, an upper limit to the background meteoroid flux at Mars is estimated to be for meteoroids with mass larger than 4 g. For comparison, the estimated flux to this mass limit at the Earth is [Grün, E., Zook, H.A., Fechtig, H., Giese, R.H., 1985. Icarus 62, 244-272]. This result is qualitatively consistent, within error bounds, with theoretical models predicting martian fluxes of ∼50% that at Earth for meteoroids of mass ¹⁰⁻³-¹⁰¹ g [Adolfsson, L.G., Gustafson, B.A.S., Murray, C.D., 1996. Icarus 119, 144-152]. The MER cameras, even using the most sensitive mode of operation, should expect to see on average only one coincident meteor on of order 40-150 h of total exposure time based on these same theoretical martian flux estimates. To more meaningfully constrain these flux models, a longer total integrated exposure time or more sensitive camera is needed. Our analysis also suggests that the event reported as the first martian meteor [Selsis, F., Lemmon, M.T., Vaubaillon, J., Bell, J.F., 2005. Nature 435, 581] is more likely a grazing cosmic ray impact, which we show to be a major source of confusion with potential meteors in all Pancam images.     

A global solution for the Mars static and seasonal gravity, Mars orientation, Phobos and Deimos masses, and Mars ephemeris

With the collection of six years of MGS tracking data and three years of Mars Odyssey tracking data, there has been a continual improvement in the JPL Mars gravity field determination. This includes the measurement of the seasonal changes in the gravity coefficients (e.g., , , , , , ) caused by the mass exchange between the polar ice caps and atmosphere. This paper describes the latest gravity field MGS95J to degree and order 95. The improvement comes from additional tracking data and the adoption of a more complete Mars orientation model with nutation, instead of the IAU 2000 model. Free wobble of the Mars' spin axis, i.e. polar motion, has been constrained to be less than 10 mas by looking at the temporal history of and . A strong annual signature is observed in , and this is a mixture of polar motion and ice mass redistribution. The Love number solution with a subset of Odyssey tracking data is consistent with the previous liquid outer core determination from MGS tracking data [Yoder et al., 2003. Science 300, 299-303], giving a combined solution of k₂=0.152±0.009 using MGS and Odyssey tracking data. The solutions for the masses of the Mars' moons show consistency between MGS, Odyssey, and Viking data sets; Phobos GM=(7.16±0.005)×¹⁰⁻⁴ km³/s² and Deimos GM=(0.98±0.07)×¹⁰⁻⁴ km³/s². Average MGS orbit errors, determined from differences in the overlaps of orbit solutions, have been reduced to 10-cm in the radial direction and 1.5 m along the spacecraft velocity and normal to the orbit plane. Hence, the ranging to the MGS and Odyssey spacecraft has resulted in position measurements of the Mars system center-of-mass relative to the Earth to an accuracy of one meter, greatly reducing the Mars ephemeris errors by several orders of magnitude, and providing mass estimates for Asteroids 1 Ceres, 2 Pallas, 3 Juno, 4 Vesta, and 324 Bamberga.     

Long term wind erosion on Mars

Wind erosion and deposition are powerful agents of surface change on Mars. Erosion is sensitive to the atmospheric density, so feedback between orbit variations and atmospheric density can enhance the sensitivity of erosion rates to orbital parameters. We use statistics derived from a 1 Gyr integration of the spin axis of Mars, coupled with runs of the NASA Ames 3-D general circulation model (GCM) at a variety of orbital conditions and pressures, to explore this feedback. We find that wind erosion in the GCM is associated with two factors: baroclinic zone winds, and strong cross-equator solstice flows. Both of these factors are influenced by topography, producing an asymmetry in the erosion pattern between the north and the south. The erosion model, averaged over 1 Gyr, produces potential erosion rates of in the north and in the south, which increase by an order of magnitude in an early 28 mbar atmosphere. The stability of these potential erosion patterns over geological time indicates that the lowland regions of Mars are continuously eroded, and that wind is capable of revealing ancient sedimentary deposits, or eroding substantial deposits that may have otherwise been preserved.     

The density of the upper martian atmosphere measured by Lyman-α absorption with Mars Express SPICAM

Absorption of interplanetary Lyman-α emission by Mars’ nightside lower thermosphere was observed by Mars Express Spectrometer for Investigation of Characteristics of the Atmosphere of Mars (SPICAM), and is analyzed to derive the CO₂ density at 110 km during a martian year. The observed density seasonal variability is consistent with recent observations obtained by stellar occultations, proving that this method, though not as accurate as stellar occultations could be used complementary to them to characterize large variations of thermospheric density on Mars and provide a better spatial coverage by Lyman-α imagery.     

Continuous monitoring of nightside upper thermospheric mass densities in the martian southern hemisphere over 4 martian years using electron reflectometry

Details are presented of an improved technique to use atmospheric absorption of magnetically reflecting solar wind electrons to constrain neutral mass densities in the nightside martian upper thermosphere. The helical motion of electrons on converging magnetic field lines, through an extended neutral atmosphere, is modeled to enable prediction of loss cone pitch angle distributions measured by the Magnetometer/Electron Reflectometer (MAG/ER) experiment on Mars Global Surveyor at 400 km altitude. Over the small fraction of Mars' southern hemisphere (∼2.5%) where the permanent crustal magnetic fields are both open to the solar wind and sufficiently strong as to dominate the variable induced martian magnetotail field, spherical harmonic expansions of the crustal fields are used to prescribe the magnetic field along the electron's path, allowing least-squares fitting of measured loss cones, in order to solve for parameters describing the vertical neutral atmospheric mass density profile from 160 to 230 km. Results are presented of mass densities in the southern hemisphere at 2 a.m. LST at the mean altitude of greatest sensitivity, 180 km, continuously over four martian years. Seasonal variability in densities is largely explained by orbital and latitudinal changes in dayside insolation that impacts the nightside through the resulting thermospheric circulation. However, the physical processes behind repeatable rapid, late autumnal cooling at mid-latitudes and near-aphelion warming at equatorial latitudes is not fully clear. Southern winter polar warming is generally weak or nonexistent over several Mars years, in basic agreement with MGS and MRO accelerometer observations. The puzzling response of mid-latitude densities from 160° to 200° E to the 2001 global dust storm suggests unanticipated localized nightside upper thermospheric lateral and vertical circulation patterns may accompany such storms. The downturn of the 11-year cycle of solar EUV flux is likely responsible for lower aphelion densities in 2004 and 2006 (Mars years 27 and 28).     

Parametric analysis of modeled ion escape from Mars

We develop a parametric fit to the results of a detailed magnetohydrodynamic (MHD) study of the response of ion escape rates (O⁺, and ) 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 R_euv, the solar wind particle density R_ρ, its velocity (squared) R_v², and the interplanetary magnetic field pressure R_B², where forcing factors are expressed in units of the current epoch-averaged values. Our parametric model is: , 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 × 10²⁴ s⁻¹, 0.782, 0.251, 0.382, and 0.214, for ?₀, α, β, γ, and δ, respectively. For , the corresponding values are 1.68 × 10²⁴ s⁻¹, −0.393, 0.798, 0.967, and 0.533. For , they are 8.66 × 10²² s⁻¹, −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 R_EUV enhances O⁺ loss, but suppresses the escape of and , whereas increases in solar wind (i.e., in , and R_B², with R_euv 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) = ?₀(i)(t/t₀)^−ξ(i). The ξ(i) are linear combinations of the epoch-averaged ion escape sensitivities, which are seen to increase with ion mass. We integrate the and oxygen ion escape rates over time, and find that in the last 3.85 Gyr, Mars would have lost about mbars of , and of water (from O⁺ and ) from ion escape.     

Photochemical escape of oxygen from Mars: A comparison of the exobase approximation to a Monte Carlo method

The non-thermal escape of neutral O atoms from Mars at the current epoch is largely due to dissociative recombination of :

     

Ion escape from Mars as a function of solar wind conditions: A statistical study

The influence of solar EUV and solar wind conditions on ion escape at Mars is investigated using ion data from the Aspera-3 instrument on Mars Express, combined with solar wind proxy data obtained from the Mars Global Surveyor (MGS) spacecraft. A solar EUV flux proxy based on data from the Earth position, scaled and shifted in time for Mars, is used to study relatively long time scale changes related to solar EUV variability. Data from May 2004 until November 2005 has been used. A clear dependence on the strength of the subsolar magnetic field as inferred from MGS measurements is seen in the ion data. The region of significant heavy ion flows is compressed and the heavy ion flux density is higher for high subsolar magnetic field strength. Because of the difference in outflow area, the difference in estimated total outflow is somewhat less than the difference in average flux density. We confirm previous findings that escaping planetary ions are mainly seen in the hemisphere into which the solar wind electric field is pointed. The effect is more pronounced for the high subsolar magnetic field case.The average ion motion has a consistent bias towards the direction of the solar wind electric field, but the main motion is in the antisunward direction. The antisunward flow velocity increases with tailward distance, reaching above at 2 to 3 martian radii downtail from Mars for O⁺ ions. Different ion species reach approximately the same bulk flow energy. We did not find any clear correlation between the solar EUV flux and the ion escape distribution or rate, probably because the variation of the solar EUV flux over our study interval was too small. The results indicate that the solar wind and its magnetic field directly interacts with the ionosphere of Mars, removing more ions for high subsolar magnetic field strength. The interaction region and the tail heavy ion flow region are not perfectly shielded from the solar wind electric field, which accelerates particles over relatively large tail distances.