First ENA observations at Mars: Subsolar ENA jet

The Neutral Particle Detector (NPD), an Energetic Neutral Atom (ENA) sensor of the Analyzer of Space Plasmas and Energetic Atoms (ASPERA-3) on board Mars Express, detected intense fluxes of ENAs emitted from the subsolar region of Mars. The typical ENA fluxes are (4-7) × 10⁵ cm⁻² sr⁻¹ s⁻¹ in the energy range 0.3-3 keV. These ENAs are likely to be generated in the subsolar region of the martian exosphere. As the satellite moved away from Mars, the ENA flux decreased while the field of view of the NPD pointed toward the subsolar region. These decreases occurred very quickly with a time scale of a few tens of seconds in two thirds of the orbits. Such a behavior can be explained by the spacecraft crossing a spatially constrained ENA jet, i.e., a highly directional ENA emission from a compact region of the subsolar exosphere. This ENA jet is highly possible to be emitted conically from the subsolar region. Such directional ENAs can result from the anisotropic solar wind flow around the subsolar region, but this can not be explained in the frame of MHD models.     

First ENA observations at Mars: Solar-wind ENAs on the nightside

We present measurements with an Energetic Neutral Atom (ENA) imager on board Mars Express when the spacecraft moves into Mars eclipse. Solar wind ions charge exchange with the extended Mars exosphere to produce ENAs that can spread into the eclipse of Mars due to the ions' thermal spread. Our measurements show a lingering signal from the Sun direction for several minutes as the spacecraft moves into the eclipse. However, our ENA imager is also sensitive to UV photons and we compare the measurements to ENA simulations and a simplified model of UV scattering in the exosphere. Simulations and further comparisons with an electron spectrometer sensitive to photoelectrons generated when UV photons interact with the spacecraft suggest that what we are seeing in Mars' eclipse are ENAs from upstream of the bow shock produced in charge exchange with solar wind ions with a non-zero temperature. The measurements are a precursor to a new technique called ENA sounding to measure solar wind and planetary exosphere properties in the future.     

Propagation of tropospheric gravity waves into the upper atmosphere of Mars

We study the propagation of gravity waves in the martian atmosphere using a linearized one-dimensional full-wave model. Calculations are carried out for atmospheric parameters characteristic of Mars Orbiter Laser Altimeter (on Mars Global Surveyor MGS) observations of apparent gravity waves in high latitude clouds and MGS radio occultation measurements of temperature variations with height suggestive of gravity wave activity. Waves that reach the thermosphere produce fluctuations in density comparable in amplitude with the density variations detected in Mars Odyssey aerobraking data. Gravity waves of modest amplitude are found to deposit momentum and generate significant heating and cooling in the martian atmosphere. The largest heating and cooling effects occur in the thermosphere, at altitudes between about 130 and 150 km, with heating occurring at the lower altitudes and cooling taking place above.     

Traveling waves in the martian atmosphere from MGS TES Nadir data

We have characterized the annual behavior of martian atmospheric traveling waves in the MGS TES data set from the first two martian years of mapping. There is a high degree of repeatability between the two years. They are dominated by strong low zonal wavenumber waves with high amplitudes near the polar jets, strongest in late northern fall and early northern winter. The m=1 waves have amplitudes up to about 20 K, are vertically extended, and occasionally extend even into the tropics. Periods for m=1 range from 2.5 to 30 sols. Much weaker waves were identified in the south, with amplitudes less than about 3.5 K. Traveling waves with m=2 and m=3 are also seen, but their amplitudes are typically limited to less than 4 K, and are generally more confined near the surface. In the north, they are more evident in fall and spring rather than winter solstice, which is clearly dominated by m=1 waves. Some evidence of storm tracks has been identified in the data, with accentuated weather-related temperature perturbations near longitudes 200° to 320° E for both the southern and northern hemispheres near latitude ±65° at the surface. Some evidence was also found for a sharpening of longitudinal gradients into what may be frontal systems. EP flux divergences show the waves extracting energy from the zonal mean winds. When the m=1 waves were strongest, decelerations of the zonal jet of order 30 m/(s sol) were measured. Above 1 scale height, the waves extract energy from the jet predominately through barotropic processes, but their character is overall mixed barotropic/baroclinic. Inertial instabilities may exist at altitude on the equatorward flanks of the polar jets, and marginal stability extends through to the tropics. This may explain the coordination of the tropical behavior of the waves with that centered along the polar jet, consistent with the ideas expressed in Wilson et al. (2002, Geophys. Res. Lett. 29, #1684) and similar to those in Barnes et al. (1993, J. Geophys. Res. 98, 3125-3148). Throughout the year, there exist large regions with the meridional gradient of PV less than zero, but they are strongest near winter solstice. Poleward of the winter jet, the regions of instability reach the surface, equatorward they do not. These regions, satisfying a necessary criterion for instability, likely explain the genesis of the waves, and perhaps also their bimodal character between surface (faster waves) and altitude (slow m=1 waves).     

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.     

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.     

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.     

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.     

The effects of topographically-controlled thermal tides in the martian upper atmosphere as seen by the MGS accelerometer

Mars Global Surveyor accelerometer observations of the martian upper atmosphere revealed large variations in density with longitude during northern hemisphere spring at altitudes of 130-160 km, all latitudes, and mid-afternoon local solar times (LSTs). This zonal structure is due to tides from the surface. The zonal structure is stable on timescales of weeks, decays with increasing altitude above 130 km, and is dominated by wave-3 (average amplitude 22% of mean density) and wave-2 (18%) harmonics. The phases of these harmonics are constant with both altitude and latitude, though their amplitudes change significantly with latitude. Near the South Pole, the phase of the wave-2 harmonic changes by 90° with a change of half a martian solar day while the wave-3 phase stays constant, suggesting diurnal and semidiurnal behaviour, respectively. We use a simple application of classical tidal theory to identify the dominant tidal modes and obtain results consistent with those of General Circulation Models. Our method is less rigorous, but simpler, than the General Circulation Models and hence complements them. Topography has a strong influence on the zonal structure.     

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.     

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 :

     

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.     

Investigations of the variability of dust particle sizes in the martian atmosphere using the NASA Ames General Circulation Model

We present a Mars General Circulation Model (GCM) numerical investigation of the physical processes (i.e., wind stress and dust devil dust lifting and atmospheric transport) responsible for temporal and spatial variability of suspended dust particle sizes. Measurements of spatial and temporal variations in airborne dust particles sizes in the martian atmosphere have been derived from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) spectral and emission phase function data [Wolff, M.J., Clancy, R.T., 2003. J. Geophys. Res. (Planets) 108 (E9), doi:10.1029/2003JE002057. 1-1; Clancy, R.T., Wolff, M.J., Christensen, P.R., 2003. J. Geophys. Res. (Planets) 108 (E9), doi:10.1029/2003JE002058. 2-1]. The range of dust particle sizes simulated by the NASA Ames GCM is qualitatively consistent with TES-derived observations of effective dust particle size variability. Model results suggest that the wind stress dust lifting scheme (which produces regionally confined dust lifting) is the process responsible for the majority of the dust particle size variability in the martian atmosphere. Additionally, model results suggest that atmospheric transport processes play an important role in the evolution of atmospheric dust particles sizes during substantial dust storms on Mars. Finally, we show that including the radiative effects of a spatially variable particle size distribution significantly influences thermal and dynamical fields during the dissipation phase of the simulated global dust storm.     

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.     

A plasma flow velocity boundary at Mars from the disappearance of electron plasma oscillations

The Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) on the Mars Express (MEX) spacecraft is capable of measuring ionospheric electron density by the use of two main methods: remote radar sounding and from the excitation of local plasma oscillations. The frequency of the locally excited electron plasma oscillations is used to measure the local electron density. However, plasma oscillations are not observed when the plasma flow velocity is higher than about 160 km/s, which occurs mainly in the solar wind and magnetosheath. As a consequence, in many passes, there is a sudden disappearance of the plasma oscillations as the spacecraft enters into the magnetosheath. This fact allows us to identify a flow velocity boundary on the dayside, between the ionosphere of Mars and the shocked solar wind. This paper summarizes the results of the measurement of 552 orbits mostly over a period from August 4, 2005 to August 17, 2007. The boundary points found using MARSIS have been verified by measurements from the Analyzer of Space Plasma and Energetic Atoms (ASPERA-3) Electron Spectrometer (ELS) instrument on Mars Express. The average position of the flow velocity boundary is compared to flow velocity simulations computed using hybrid model and other boundaries. The boundary altitude is slightly lower than the magnetic pile-up boundary determined using Phobos 2 and Mars Global Surveyor (MGS) crossings, but it is in good agreement with the induced magnetospheric boundary determined by ASPERA-3. Investigation of the effect of the crustal magnetic field revealed that the flow velocity boundary is raised at the locations with strong crustal magnetic fields.     

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.     

Oxygen and carbon isotope ratios in the martian atmosphere

Oxygen and carbon isotope ratios in the martian CO₂ are key values to study evolution of volatiles on Mars. The major problems in spectroscopic determinations of these ratios on Mars are uncertainties associated with: (1) equivalent widths of the observed absorption lines, (2) line strengths in spectroscopic databases, and (3) thermal structure of the martian atmosphere during the observation. We have made special efforts to reduce all these uncertainties. We observed Mars using the Fourier Transform Spectrometer at the Canada–France–Hawaii Telescope. While the oxygen and carbon isotope ratios on Mars were byproducts in the previous observations, our observation was specifically aimed at these isotope ratios. We covered a range of 6022 to 6308 cm⁻¹ with the highest resolving power of ν/δν=3.5×10⁵ and a signal-to-noise ratio of 180 in the middle of the spectrum. The chosen spectral range involves 475 lines of the main isotope, 184 lines of ¹³CO₂, 181 lines of CO¹⁸O, and 119 lines of CO¹⁷O. (Lines with strengths exceeding 10⁻²⁷ cm at 218 K are considered here.) Due to the high spectral resolution, most of the lines are not blended. Uncertainties of retrieved isotope abundances are in inverse proportion to resolving power, signal-to-noise ratio, and square root of the number of lines. Laboratory studies of the CO₂ isotope spectra in the range of our observation achieved an accuracy of 1% in the line strengths. Detailed observations of temperature profiles using MGS/TES and data on temperature variations with local time from two GCMs are used to simulate each absorption line at various heights in each part of the instrument field of view and then sum up the results. Thermal radiation of Mars' surface and atmosphere is negligible in the chosen spectral range, and this reduces errors associated with uncertainties in the thermal structure on Mars. Using a combination of all these factors, the highest accuracy has been achieved in measuring the CO₂ isotope ratios: ¹³C/¹²C = 0.978 ± 0.020 and ¹⁸O/¹⁶O = 1.018 ± 0.018 times the terrestrial standards. Heavy isotopes in the atmosphere are enriched by nonthermal escape and sputtering, and depleted by fractionation with solid-phase reservoirs. The retrieved ratios show that isotope fractionation between CO₂ and oxygen and carbon reservoirs in the solid phase is almost balanced by nonthermal escape and sputtering of O and C from Mars.     

The loss of water from Mars: Numerical results and challenges

The simulation of Mars is a very challenging effort. However, simulations are a major method of addressing the issues of the solar wind interaction with Mars. Further, it is via simulations that issues such as water loss from Mars via solar wind pick up of ionospheric ions will be addressed. This paper discusses some of the issues raised during the Chapman Conference on Solar Wind Interactions with Mars, SWIM. It also addresses numerical issues and the authors attempts to address them, coupled with results of preliminary simulations of Mars.