Impact of atmospheric gravity waves on the jovian ionosphere

We investigate the effects of atmospheric gravity waves on the vertical and horizontal structure of the ionosphere of Jupiter. The presented non-linear, two-dimensional model of the jovian ionosphere allows for spatially and temporally varying neutral wind and temperature fields and tracks the time evolution of six ionospheric species, , and . An analytical approach is used to validate the model results for linear, small-amplitude waves and to elucidate the mechanisms that leads to perturbations in the density of the main ion species, H⁺ and . We demonstrate that the long-lived H⁺ ions are perturbed directly by wave dynamics whereas short-lived ions such as are perturbed by chemical interactions with other perturbed ion species. The model is then applied using larger gravity wave amplitudes consistent with observations. Atmospheric gravity waves propagating at high altitudes create layers of enhanced electron density similar to the system of layers observed during the J0-ingress radio occultation of the Galileo spacecraft. Our best fit to the J0-ingress observation is achieved using an 82 min period forcing wave with horizontal and vertical wavelengths of 500 km and 60 km respectively, and peaks at 510 km above the 1 bar pressure level. We further investigate the effects of the wave-induced ion flux on the background ionospheric structure and demonstrate that in the presence of a gravity wave the background density profiles of the H⁺ and ions are significantly modified. We also find that the column density of has variations that can exceed 10% as the wave propagates.     

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

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.     

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.     

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.     

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.     

A solar spectrum for PFS data analysis

A measured calibrated solar radiance in the range 1.2-, with the spectral sampling of does not exist. When studying the measured Planetary Fourier Spectrometer (PFS) spectra of the Earth's or Mars's atmosphere we discover that the most used solar spectrum contains several important errors. Here we present a “calibrated” solar radiance in the wavelength range 1.2-, with the spectral resolution of PFS , which we are going to use for studying Martian spectra. This spectrum has been assembled using measurements from Kitt Peak and from ATMOS Spacelab experiment (uncalibrated high resolution) and theoretical results, together with low resolution calibrated continuum. This is the best we can have in this moment to be used with PFS, while waiting to have good solar calibrated radiances. Examples of solar lines at Mars are given.     

The structure and time variability of the ring atmosphere and ionosphere

The saturnian system is subject to constant bombardment by interplanetary meteoroids and irradiation by solar UV photons. Both effects release neutral molecules from the icy ring particles either in the form of impact water vapor or gas emission in the form of H₂O, O₂ and H₂. The observations of the Cassini spacecraft during its orbit insertion have shown the existence of molecular and atomic oxygen ions. Subsequent modeling efforts have led to the picture that an exospheric population of neutral oxygen molecules is probably maintained in the vicinity of the rings via photolytic-decomposition of ice and surface reactions. At the same time, ionized products O⁺ and ions move along the magnetic field lines and, depending on the optical local thickness rings, can thread through the ring plane or impact a ring particle, the ion principal sink. In addition, collisional interactions between the ions and neutrals will change the scale height of the ions and produce a scattered component of O₂ molecules and O atoms which can be injected into Saturn’s upper atmosphere or the inner magnetosphere. The ring atmosphere, therefore, serves as a source of ions throughout Saturn’s magnetosphere. If photolysis of ice is the dominant source of O₂, then the complex structure of the ring atmosphere/ionosphere and the injection rate of neutral O₂ will be subject to modulation by the seasonal variation of Saturn along its orbit. In this work, we show how the physical properties of the ring oxygen atmosphere, the scattered component, and the magnetospheric ion source rate vary as the ring system goes through the cycle of solar insolation. In particular, it is shown that the magnetopheric ions should be nearly depleted at Saturn’s equinox if O₂ is produced mainly by photolysis of the ring material.     

Calculation of the enrichment of the giant planet envelopes during the “late heavy bombardment”

The giant planets of our solar system possess envelopes consisting mainly of hydrogen and helium but are also significantly enriched in heavier elements relatively to our Sun. In order to better constrain how these heavy elements have been delivered, we quantify the amount accreted during the so-called “late heavy bombardment”, at a time when planets were fully formed and planetesimals could not sink deep into the planets. On the basis of the “Nice model”, we obtain accreted masses (in terrestrial units) equal to for Jupiter, and for Saturn. For the two other giant planets, the results are found to depend mostly on whether they switched position during the instability phase. For Uranus, the accreted mass is with an inversion and without an inversion. Neptune accretes in models in which it is initially closer to the Sun than Uranus, and otherwise. With well-mixed envelopes, this corresponds to an increase in the enrichment over the solar value of 0.033±0.001 and 0.074±0.007 for Jupiter and Saturn, respectively. For the two other planets, we find the enrichments to be 2.1±1.4 (w/ inversion) or 1.2±0.7 (w/o inversion) for Uranus, and 2.0±1.2 (w/ inversion) or 2.7±1.6 (w/o inversion) for Neptune. This is clearly insufficient to explain the inferred enrichments of ∼4 for Jupiter, ∼7 for Saturn and ∼45 for Uranus and Neptune.     

Diagnostic of the dayside ionosphere of Mars using the Total Electron Content measurement by the NEIGE/Netlander experiment: An assessment study

The NEtlander Ionosphere and Geodesy Experiment (NEIGE) of the Netlander mission to Mars will measure Doppler shifts affecting the radio links between ground stations and an orbiter. The experiment has two complementary scientific objectives which are the monitoring of the structure and dynamics of the ionosphere of Mars and the precise determination of Mars orientation parameters. The horizontal variation of the Total Electron Content (TEC) of the ionosphere will be derived from the so-called “geometric-free” combination of the Doppler shifts which affect radio links at two frequencies (in the UHF and S bands) between the Netlander microstations on the Mars surface and the data-relay orbiter. We describe a new method for retrieving the horizontal profile of the absolute value of the TEC. Simulations have allowed to evaluate the precision in the determination of the TEC using this method. We show that the daytime TEC can be retrieved with a precision of for a nominal accuracy of on the orbital pseudo-velocity, which represents a relative precision of a few percent. A preliminary analysis of the sensitivity of the TEC to the physical parameters which control the ionosphere has been performed. For this purpose, we have used a new one-dimensional ionospheric model based on the solution of coupled kinetic, fluid and MHD equations. This model describes the suprathermal electron component, the thermal plasma component as well as the induced horizontal magnetic field. The code which provides the vertical electron density profile has been used to study the variation of the TEC with the solar zenith angle and with the induced magnetic field at the top of the ionosphere. In particular, we show that NEIGE will allow to diagnose the penetration into the daytime ionosphere of an induced magnetic field.     

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 :

     

Thermal disturbances caused by lander shadowing and the measurability of the martian planetary heat flow

A measurement of the martian planetary heat flow requires the determination of the subsurface temperature gradient, which is affected by surface insolation. I investigate the propagation of thermal disturbances caused by lander shadowing and derive measurement requirements for in situ heat flow experiments. I find that for short term measurements spanning 180 sol, a measurement depth of at least 2 m is needed to guarantee a stable thermal environment directly underneath the lander for Moon-like thermal conductivities of . For extremely large conductivities of , this depth needs to be increased to 4 m, but if the probe can be deployed outside the lander structure, the respective depths can be decreased by 1 m. For long term measurements spanning at least a full martian year heat flow perturbations are smaller than 5% below a depth of 3 m directly underneath the lander. Outside the lander structure, essentially unperturbed measurements may be conducted at depths of 0.5 and 1.5 m for thermal conductivities of 0.02 and , respectively.     

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.     

On the mean oxygen isotope composition of the Solar System

Since the first discovery of extraordinary oxygen isotope compositions in carbonaceous meteorites by Clayton et al. [Clayton, R.N., Grossman, L., Mayeda, T.K., 1973. Science 182, 485-488], numerous studies have been done to explain the unusual mass-independent isotope fractionation, but the problem is still unresolved to this day. Clayton's latest interpretation [Clayton, R.N., 2002. Nature 415, 860-861] sheds new light on the problem, and possible hypotheses now seem to be fairly well defined. A key issue is to resolve whether the oxygen isotopes in the Solar System represented by the Sun (solar oxygen) are the same as oxygen isotopes in planetary objects such as bulk meteorites, Mars, Earth, and Moon, or whether the solar oxygen is more similar to the lightest oxygen isotopes observed in CAIs (Calcium Aluminum-rich Inclusions) in primitive meteorites. Here, we examined the problem using oxygen isotope analytical data of about 400 bulk meteorite samples of various classes or types (data compiled by K. Lodders). We used in our discussion exclusively the parameter , a direct measure of the degree of mass-independent isotope fractionation of oxygen isotopes. When is arranged according to a characteristic size of their host planetary object, it shows a systematic trend: (1) values scatter around zero; (2) the scatter from the mean () decreases with increasing representative size of the respective host planetary object. This systematic trend is easily understood on the basis of a hierarchical scenario of planetary formation, that is, larger planetary objects have formed by progressive accretion of planetesimals by random sampling over a wide spectrum of proto-solar materials. If this progressive random sampling of planetesimals were the essential process of planetary formation, the isotopic composition of planetary oxygen should approach that of the solar oxygen. To test this random sampling hypothesis, we applied a multiscale, multistep bootstrap statistical method [Shimodaira, H., 2004. Ann. Statist. 32, 2616-2641] to the meteorite oxygen isotope data, and deduced a σ-N relation, where σ is the standard deviation of , and N is the representative size of a host planetary object. If we assign 200 and 500 km as a representative sizes of the chondrite and achondrite parent bodies, the observed σ of agree well with the values predicted by the σ-N relation. A common mean value of for all planetary objects also agrees with the progressive random sampling process. Therefore, we conclude that the solar oxygen is the same as planetary oxygen, but differs from CAI oxygen. The conclusion implies that a massive enrichment in ¹⁷O and ¹⁸O resulting from CO self-shielding, a current influential interpretation of CAI-O, did not occur.     

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:

     

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.