Mars high resolution gravity fields from MRO, Mars seasonal gravity, and other dynamical parameters

With 2 years of tracking data collection from the MRO spacecraft, there is noticeable improvement in the high frequency portion of the spherical harmonic Mars gravity field. The new JPL Mars gravity fields, MRO110B and MRO110B2, show resolution near degree 90. Additional years of MGS and Mars Odyssey tracking data result in improvement for the seasonal gravity changes which compares well to global circulation models and Odyssey neutron data and Mars rotation and precession (). Once atmospheric dust is accounted for in the spacecraft solar pressure model, solutions for Mars solar tide are consistent between data sets and show slightly larger values (k₂ = 0.164 ± 0.009, after correction for atmospheric tide) compared to previous results, further constraining core models. An additional 4 years of Mars range data improves the Mars ephemeris, determines 21 asteroid masses and bounds solar mass loss (dGM_Sun/dt < 1.6 × 10⁻¹³ GM_Sun year⁻¹).     

Extension of atmospheric dust loading to high altitudes during the 2001 Mars dust storm: MGS TES limb observations

Mars Global Surveyor (MGS) visible (solarband bolometer) and thermal infrared (IR) spectral limb observations from the Thermal Emission Spectrometer (TES) support quantitative profile retrievals for dust opacity and particle sizes during the 2001 global dust event on Mars. The current analysis considers the behavior of dust lifted to altitudes above 30 km during the course of this storm; in terms of dust vertical mixing, particle sizes, and global distribution. TES global maps of visible (solarband) limb brightness at 60 km altitude indicate a global-scale, seasonally evolving (over 190-240° solar longitudes, L_S) longitudinal corridor of vertically extended dust loading (which may be associated with a retrograde propagating, wavenumber 1 Rossby wave). Spherical radiative transfer analysis of selected limb profiles for TES visible and thermal IR radiances provide quantitative vertical profiles of dust opacity, indicating regional conditions of altitude-increasing dust mixing ratios. Observed infrared spectral dependences and visible-to-infrared opacity ratios of dust scattering over 30-60 km altitudes indicate particle sizes characteristic of lower altitudes (cross-section weighted effective radius, ), during conditions of significant dust transport to these altitudes. Conditions of reduced dust loading at 30-60 km altitudes present smaller dust particle sizes . These observations suggest rapid meridional transport at 30-80 km altitudes, with substantial longitudinal variation, of dust lifted to these altitudes over southern hemisphere atmospheric regions characterized by extraordinary (m/s) vertical advection velocities. By L_S=230° dust loading above 50 km altitudes decreased markedly at southern latitudes, with a high altitude (60-80 km) haze of fine (likely) water ice particles appearing over 10°S-40°N latitudes.     

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.     

The Planetary Fourier Spectrometer (PFS) onboard the European Mars Express mission

The Planetary Fourier Spectrometer (PFS) for the Mars Express mission is an infrared spectrometer optimised for atmospheric studies. This instrument has a short wave (SW) channel that covers the spectral range from 1700 to (1.2-) and a long-wave (LW) channel that covers 250- (5.5-). Both channels have a uniform spectral resolution of . The instrument field of view FOV is about 1.6^° (FWHM) for the Short Wavelength channel (SW) and 2.8^° (FWHM) for the Long Wavelength channel (LW) which corresponds to a spatial resolution of 7 and 12 km when Mars is observed from an height of 250  km. PFS can provide unique data necessary to improve our knowledge not only of the atmosphere properties but also about mineralogical composition of the surface and the surface-atmosphere interaction.The SW channel uses a PbSe detector cooled to 200-220 K while the LW channel is based on a pyroelectric (L_iT_aO₃) detector working at room temperature. The intensity of the interferogram is measured every 150 nm of physical mirrors displacement, corresponding to 600 nm optical path difference, by using a laser diode monochromatic light interferogram (a sine wave), whose zero crossings control the double pendulum motion. PFS works primarily around the pericentre of the orbit, only occasionally observing Mars from large distances. Each measurements take 4 s, with a repetition time of 8.5 s. By working roughly 0.6 h around pericentre, a total of 330 measurements per orbit will be acquired 270 looking at Mars and 60 for calibrations. PFS is able to take measurements at all local times, facilitating the retrieval of surface temperatures and atmospheric vertical temperature profiles on both the day and the night side.     

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.     

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.     

The dust trail of Comet 67P/Churyumov-Gerasimenko between 2004 and 2006

We report on observations of the dust trail of Comet 67P/Churyumov-Gerasimenko (CG) in visible light with the Wide Field Imager at the ESO/MPG 2.2 m telescope at 4.7 AU before aphelion, and at with the MIPS instrument on board the Spitzer Space Telescope at 5.7 AU both before and after aphelion. The comet did not appear to be active during our observations. Our images probe large dust grains emitted from the comet that have a radiation pressure parameter β<0.01. We compare our observations with simulated images generated with a dynamical model of the cometary dust environment and constrain the emission speeds, size distribution, production rate and geometric albedo of the dust. We achieve the best fit to our data with a differential size distribution exponent of −4.1, and emission speeds for a β=0.01 particle of 25 m/s at perihelion and 2 m/s at 3 AU. The dust production rate in our model is on the order of 1000 kg/s at perihelion and 1 kg/s at 3 AU, and we require a dust geometric albedo between 0.022 and 0.044. The production rates of large (>) particles required to reproduce the brightness of the trail are sufficient to also account for the coma brightness observed while the comet was inside 3 AU, and we infer that the cross-section in the coma of CG may be dominated by grains of the order of .     

Environmental predictions for the Beagle 2 lander, based on GCM climate simulations

The Mars climate database (MCD) is a database of statistics based on output from physically consistent numerical model simulations which describe the climate and surface environment of Mars. It is used here to predict the meteorological environment of the Beagle 2 lander site. The database was constructed directly on the basis of output from multiannual integrations of two general circulation models, developed jointly at Laboratoire de Météorologie Dynamique du Centre National de la Recherche Scientifique, France, and the University of Oxford, UK. In an atmosphere with dust opacities similar to that observed by Mars Global Surveyor, the predicted surface temperature at the time of landing (L_s=322°, 13:00 local time), is , and varying between ∼186 and over the Martian day. The predicted air temperature at above the surface, the height of the fully extended Beagle 2 robot arm, is at the time of landing. The expected mean wind near the surface on landing is north-westerly in direction, becoming more southerly over the mission. An increase in mean surface pressure is expected during the mission. Heavy global dust storm predictions are discussed; conditions which may only occur in the extreme as the expected time of landing is around the end of the main dust storm period. Past observations show approximately a one in five chance of a large-scale dust storm in a whole Mars year over the landing region, Isidis Planitia. This statistic results from observations of global, encircling, regional and local dust storms but does not include any small-scale dust “events” such as dust devils.     

Observations of the inner coma of C/1995O1 (Hale-Bopp)—gas and dust production

Hale-Bopp (C/1995 O1) was the most productive recent comet observed in terms of gas and dust output. Since its discovery in 1995 at a distance of 7.14 AU from the Sun, the comet has been well observed, revealing the dynamics of a rare and large comet. Hale-Bopp showed strong emissions of the principle cometary gases CN, C₃, and C₂, as well as an abundance of dust. The production rates of these gases were found to be 1.45×10²⁸, 1.71×10²⁸, and , respectively, with dust production, in terms of Afρ, , as measured in the green continuum (5260 Å). The observations for this paper are presented in two groups spanning 10 days each, one group centered near 32 days prior to and the other 21 days after perihelion. The averages of dust and gas production rates show a slightly higher value for each prior to perihelion than after perihelion, consistent with a possible peak in production a few weeks prior to perihelion passage.     

Low-speed impacts between rubble piles modeled as collections of polyhedra, 2

We present the results of additional calculations involving the collisions of km-scale rubble piles. In new work, we used the Open Dynamics Engine (ODE), an open-source library for the simulation of rigid-body dynamics that incorporates a sophisticated collision-detection and resolution routine. We found that using ODE resulted in a speed-up of approximately a factor of 30 compared with previous code. In this paper we report on the results of almost 1200 separate runs, the bulk of which were carried out with 1000-2000 elements. We carried out calculations with three different combinations of the coefficients of friction η and (normal) restitution ?: low (η=0,?=0.8), medium (η=0,?=0.5), and high (η=0.5,?=0.5) dissipation.For target objects of ∼1 km in radius, we found reduced critical disruption energy values in head-on collisions from 2 to 100 J kg⁻¹ depending on dissipation and impactor/target mass ratio. Monodisperse objects disrupted somewhat more easily than power-law objects in general. For oblique collisions of equal-mass objects, mildly off-center collisions (b/b₀=0.5) seemed to be as efficient or possibly more efficient at collisional disruption as head-on collisions. More oblique collisions were less efficient and the most oblique collisions we tried (b/b₀=0.866) required up to ∼200 J kg⁻¹ for high-dissipation power-law objects. For calculations with smaller numbers of elements (total impactor or 200 elements) we found that collisions were more efficient for smaller numbers of more massive elements, with values as low as for low-dissipation cases. We also analyzed our results in terms of the relations proposed by Stewart and Leinhardt [Stewart, S.T., Leinhardt, Z.M., 2009. Astrophys. J. 691, L133-L137] where where Q_R is the impact kinetic energy per unit total mass m_i+m_T. Although there is a significant amount of scatter, our results generally bear out the suggested relation.     

Airglow on Mars: Some model expectations for the OH Meinel bands and the O2 IR atmospheric band

This work presents model calculations of the diurnal airglow emissions from the OH Meinel bands and the O₂ IR atmospheric band in the neutral atmosphere of Mars. A time-dependent photochemical model of the lower atmosphere below 80 km has been developed for this purpose. Special emphasis is placed on the nightglow emissions because of their potential to characterize the atomic oxygen profile in the 50-80 km region. Unlike on Earth, the OH Meinel emission rates are very sensitive to the details of the vibrational relaxation pathway. In the sudden death and collisional cascade limits, the maximum OH Meinel column intensities for emissions originating from a fixed upper vibrational level are calculated to be about 300 R, for transitions v^′=9→v^′?8, and 15,000 R, for transitions v^′=1→v^′=0, respectively. During the daytime the 1.27 μm emission from O₂(), primarily formed from ozone photodissociation, is of the order of MegaRayleighs (MR). Due to the long radiative lifetime of O₂(), a luminescent remnant of the dayglow extends to the dark side for about two hours. At night, excited molecular oxygen is expected to be produced through the three body reaction O + O + CO₂. The column emission of this nighttime component of the airglow is estimated to amount to 25 kR. Both nightglow emissions, from the OH Meinel bands and the O₂ IR atmospheric band, overlap in the 50-80 km region. Photodissociation of CO₂ in the upper atmosphere and the subsequent transport of the atomic oxygen produced to the emitting layer are revealed as key factors in the nightglow emissions from these systems. The Mars 5 upper constraint for the product [H][O₃] is revised on the basis of more recent values for the emission probabilities and collisional deactivation coefficients.