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.     

On the population, physical decay and orbital distribution of Jupiter family comets: Numerical simulations

We study the Jupiter family comet (JFC) population assumed to come from the Scattered Disk and transferred to the Jupiter’s zone through gravitational interactions with the Jovian planets. We shall define as JFCs those with orbital periods and Tisserand parameters in the range 2<T?3.1, while those comets coming from the same source, but that do not fulfill the previous criteria (mainly because they have periods ) will be called ‘non-JFCs’. We performed a series of numerical simulations of fictitious comets with a purely dynamical model and also with a more complete dynamical-physical model that includes besides nongravitational forces, sublimation and splitting mechanisms. With the dynamical model, we obtain a poor match between the computed distributions of orbital elements and the observed ones. However with the inclusion of physical effects in the complete model we are able to obtain good fits to observations. The best fits are attained with four splitting models with a relative weak dependence on q, and a mass loss in every splitting event that is less when the frequency is high and vice versa. The mean lifetime of JFCs with radii and is found to be of about 150-200 revolutions (∼. The total population of JFCs with radii within Jupiter’s zone is found to be of 450±50. Yet, the population of non-JFCs with radii in Jupiter-crossing orbits may be ∼4 times greater, thus leading to a whole population of JFCs + non-JFCs of ∼2250±250. Most of these comets have perihelia close to Jupiter’s orbit. On the other hand, very few non-JFCs reach the Earth’s vicinity (perihelion distances ) which gives additional support to the idea that JFCs and Halley-type comets have different dynamical origins. Our model allows us to define the zones of the orbital element space in which we would expect to find a large number of JFCs. This is the first time, to our knowledge, that a physico-dynamical model is presented that includes sublimation and different splitting laws. Our work helps to understand the role played by these erosion effects in the distribution of the orbital elements and lifetimes of JFCs.     

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.     

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.     

Titan haze distribution and optical properties retrieved from recent observations

The VIMS instrument onboard Cassini observed the north polar region of Titan at 113° phase angle, 28 December 2006. On this spectral image, a vast polar cloud can be seen northward to 62°N, and elsewhere, the haze appears as the dominant source of scattering. Because the surface does not appear in the wavelength range between 0.3 and , this spectro-image is ideal to study airborn scatterers both in methane bands and windows. In this work, we study this image, along with another image taken at 13° phase angle. This image probe both the atmosphere and the surface from pole to pole. First, we characterise the spatial distribution of the haze layer above 100 km between 80°S and 70°N. We find a north south asymmetry with a haze opacity increasing by a factor 3 from the south pole to the equator, then a constant value up to about 30°N and a decrease of a factor 2 between 30°N and about 60°N. Beyond 60°N, we can see the influence of the north polar cloud, even in the band, but no polar haze accumulation. The fact that the north polar region is still in the polar night is a possible explanation. No haze accumulation is observed in the southern polar region. Secondly, we partly identify the origin of spectral features in the 2.8-μm methane window, which are found to be due to deuterated methane (CH₃D). This allows the analyse of this window and to retrieve the opacity of scatterers layer below 80 km (presumably made of aerosols and condensate droplets) between 35°N and 60°N. Finally, we constrained the values and the spectral behaviour of the imaginary part of the aerosol refractive index in the range between 0.3 and . To do so, we studied the 2.8-μm window with the image taken at 113° phase angle. To complete the analysis, we studied the transmission through the haze layer in the 3.4-μm band observed in solar occultation mode with VIMS, and we analysed the single scattering albedo retrieved with DISR instrument between 0.4 and . The imaginary part of the refractive index that we find for Titan aerosols follows Khare et al. (Khare, B.N. et al. [1984]. Icarus 60, 127-137) optical constant up to and becomes constant beyond this wavelength at least up to . It also has a prominent peak at and a secondary peak at , which indicates material rich in C-H bonds, with much less N-H bonds than in Khare et al. (1984) tholins.     

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:

     

The E ring in the vicinity of Enceladus: I. Spatial distribution and properties of the ring particles

Saturn's diffuse E ring is the largest ring of the Solar System and extends from about (Saturn radius RS=60,330 km) to at least encompassing the icy moons Mimas, Enceladus, Tethys, Dione, and Rhea. After Cassini's insertion into her saturnian orbit in July 2004, the spacecraft performed a number of equatorial as well as steep traversals through the E ring inside the orbit of the icy moon Dione. Here, we report about dust impact data we obtained during 2 shallow and 6 steep crossings of the orbit of the dominant ring source—the ice moon Enceladus. Based on impact data of grains exceeding 0.9 μm we conclude that Enceladus feeds a torus populated by grains of at least this size along its orbit. The vertical ring structure at agrees well with a Gaussian with a full-width-half-maximum (FWHM) of ∼4200 km. We show that the FWHM at is due to three-body interactions of dust grains ejected by Enceladus' recently discovered ice volcanoes with the moon during their first orbit. We find that particles with initial speeds between 225 and 235 m s⁻¹ relative to the moon's surface dominate the vertical distribution of dust. Particles with initial velocities exceeding the moon's escape speed of 207 m s⁻¹ but slower than 225 m s⁻¹ re-collide with Enceladus and do not contribute to the ring particle population. We find the peak number density to range between 16×¹⁰⁻² m⁻³ and 21×¹⁰⁻² m⁻³ for grains larger 0.9 μm, and 2.1×¹⁰⁻² m⁻³ and 7.6×¹⁰⁻² m⁻³ for grains larger than 1.6 μm. Our data imply that the densest point is displaced outwards by at least with respect of the Enceladus orbit. This finding provides direct evidence for plume particles dragged outwards by the ambient plasma. The differential size distribution for grains >0.9 μm is described best by a power law with slopes between 4 and 5. We also obtained dust data during ring plane crossings in the vicinity of the orbits of Mimas and Tethys. The vertical distribution of grains >0.8 μm at Mimas orbit is also well described by Gaussian with a FWHM of ∼5400 km and displaced southwards by ∼1200 km with respect to the geometrical equator. The vertical distribution of ring particles in the vicinity of Tethys, however, does not match a Gaussian. We use the FWHM values obtained from the vertical crossings to establish a 2-dimensional model for the ring particle distribution which matches our observations during vertical and equatorial traversals through the E ring.     

Orbits and masses of Saturn's coorbital and F-ring shepherding satellites

We have obtained numerically integrated orbits for Saturn's coorbital satellites, Janus and Epimetheus, together with Saturn's F-ring shepherding satellites, Prometheus and Pandora. The orbits are fit to astrometric observations acquired with the Hubble Space Telescope and from Earth-based observatories and to imaging data acquired from the Voyager spacecraft. The observations cover the 38 year period from the 1966 Saturn ring plane crossing to the spring of 2004. In the process of determining the orbits we have found masses for all four satellites. The densities derived from the masses for Janus, Epimetheus, Prometheus, and Pandora in units of g cm⁻³ are , , , and , respectively.     

Microgravity experiments on the collisional behavior of saturnian ring particles

In this paper we present results of two novel experimental methods to investigate the collisional behavior of individual macroscopic icy bodies. The experiments reported here were conducted in the microgravity environments of parabolic flights and the Bremen drop tower facility. Using a cryogenic parabolic-flight setup, we were able to capture 41 near-central collisions of 1.5-cm-sized ice spheres at relative velocities between 6 and . The analysis of the image sequences provides a uniform distribution of coefficients of restitution with a mean value of and values ranging from ε=0.06 to 0.84. Additionally, we designed a prototype drop-tower experiment for collisions within an ensemble of up to one hundred cm-sized projectiles and performed the first experiments with solid glass beads. We were able to statistically analyze the development of the kinetic energy of the entire system, which can be well explained by assuming a granular ‘fluid’ following Haff’s law with a constant coefficient of restitution of ε=0.64. We could also show that the setup is suitable for studying collisions at velocities of <5 mm s⁻¹ appropriate for collisions between particles in Saturn’s dense main rings.     

Eclipsing binary Trojan asteroid Patroclus: Thermal inertia from Spitzer observations

We present mid-infrared observations of the binary L5-Trojan system (617) Patroclus-Menoetius before, during, and after two shadowing events, using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. For the first time, we effectively observe changes in asteroid surface temperature in real time, allowing the thermal inertia to be determined very directly. A new detailed binary thermophysical model is presented which accounts for the system’s known mutual orbit, arbitrary component shapes, and thermal conduction in the presence of eclipses.We obtain two local thermal-inertia values, representative of the respective shadowed areas: and . The average thermal inertia is estimated to be , potentially with significant surface heterogeneity. This first thermal-inertia measurement for a Trojan asteroid indicates a surface covered in fine regolith. Independently, we establish the presence of fine-grained (<a few μm) silicates on the surface, based on emissivity features near 10 and similar to those previously found on other Trojans.We also report V-band observations and report a lightcurve with complete rotational coverage. The lightcurve has a low amplitude of peak-to-peak, implying a roughly spherical shape for both components, and is single-periodic with a period equal to the period of the mutual orbit, indicating that the system is fully synchronized.The diameters of Patroclus and Menoetius are 106±11 and , respectively, in agreement with previous findings. Taken together with the system’s known total mass, this implies a bulk mass density of , significantly below the mass density of L4-Trojan asteroid (624) Hektor and suggesting a bulk composition dominated by water ice.All known physical properties of Patroclus, arguably the best studied Trojan asteroid, are consistent with those expected in icy objects with devolatilized surface (extinct comets), consistent with what might be implied by recent dynamical modeling in the framework of the Nice Model.     

Numerical models of Saturn's long-lived anticyclones

New measurements of the dynamical properties of the long-lived Saturn's anticyclonic vortex known as “Brown Spot” (BS), discovered during the Voyager 1 and 2 flybys in 1980-1981 at latitude 43.1° N, and model simulations using the EPIC code, have allowed us to constrain the vertical wind shear and static stability in Saturn's atmosphere (vertically from pressure levels from 10 mbar to 10 bars) at this latitude. BS dynamical parameters from Voyager images include its size as derived from cloud albedo gradient (6100 km East-West times 4300 km North-South), mean tangential velocity ( at 2400 km from center) and mean vorticity (4.0±1.5×¹⁰⁻⁵ s⁻¹), lifetime >1 year, drift velocity relative to Voyager's System III rotation rate, mean meridional atmospheric wind profile at cloud level at its latitude and interactions with nearby vortices (pair orbiting and merging). An extensive set of numerical experiments have been performed to try to reproduce this single vortex properties and its observed mergers with smaller anticyclones by varying the vertical structure of the zonal wind and adjusting the static stability of the lower stratosphere and upper troposphere. Within the context of the EPIC model atmosphere, our simulations indicate that BS's drift velocity, longevity and merging behavior are very sensitive to these two atmospheric properties. The best results at the BS latitude occur for static stability conditions that use a Brunt-Väisäla frequency constant in the upper troposphere (from 0.5 to 10 bar) above 3.2×¹⁰⁻³ s⁻¹ and suggest that the wind speed slightly decays below the visible cloud deck from ∼0.5 to 10 bar at a rate per scale height. Changing the vortex latitude within the band domain introduces latitude oscillations in the vortex but not a significant meridional migration. Simulated mergers always showed orbiting movements with a typical merging time of about three days, very close to the time-span observed in the interaction of real vortices. Although these results are not unique in view of the unknowns of Saturn's deep atmosphere, they serve to constrain realistically its structure for ongoing Cassini observations.     

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.     

Characteristics of Saturn's polar atmosphere and auroral electrons derived from HST/STIS, FUSE and Cassini/UVIS spectra

Ultraviolet (UV) spectra of Saturn's aurora obtained with the Hubble Space Telescope Imaging Spectrograph (STIS), the Cassini Ultraviolet Imaging Spectrograph (UVIS) and the Far Ultraviolet Spectroscopic Explorer (FUSE) have been analyzed. Comparisons between the observed spectra and synthetic models of electron-excited H₂ have been used to determine various auroral characteristics. Far ultraviolet (FUV: 1200-1700 Å) STIS and UVIS spectra exhibit, below 1400 Å, weak absorption due to methane, with a vertical column ranging between 1.4×¹⁰¹⁵ and . Using the low-latitude Moses et al. [Moses, J.I., Bézard, B., Lellouch, E., Feuchtgruber, H., Gladstone, G.R., Allen, M., 2000. Icarus, 143, 244-298] atmospheric model of Saturn and an electron energy-H₂ column relationship, these methane columns are converted into the mean energy of the primary precipitating electrons, estimated to lie in the range 10-18 keV. This result is confirmed by the study of self-absorption with UVIS and FUSE extreme ultraviolet (EUV: 900-1200 Å) spectra. Below 1200 Å, it is seen that transitions connecting to the v^″<2 vibrational levels of the H₂ electronic ground state are partially self-absorbed by H₂ molecules overlying the auroral emission. Because of its low spectral resolution (∼5.5 Å), the UVIS EUV spectrum we analyzed does not allow us to unequivocally determine reasonable ranges of temperatures and H₂ columns. On the other hand, the high spectral resolution (∼0.2 Å) of the FUSE LiF1a and LiF2a EUV spectra we examined resolve the H₂ rotational lines and makes it possible to determine the H₂ temperature. The modeled spectrum best fitting the FUSE LiF1a observation reveals a temperature of 500 K and self-absorption by a H₂ vertical column of . When converted to energy of precipitating electrons, this H₂ column corresponds to primary electrons of ∼10 keV. The model that best fits the LiF2a spectrum is characterized by a temperature of 400 K and is not self-absorbed, making this segment ideal to determine the H₂ temperature at the altitude of the auroral emission. The latter value is in agreement with temperatures obtained from infrared polar spectra. Self-absorption is detectable in the LiF2a segment for H₂ columns exceeding , which sets the maximum mean energy determined from the FUSE observations to ∼15 keV. The total electron energy range of 10-18 keV deduced from FUV and EUV observations places the auroral emission peak between the 0.1 and 0.3 μbar pressure levels. These values should be seen as an upper limit, since most of the Voyager UVS spectra of Saturn's aurora examined by Sandel et al. [Sandel, B.R., Shemansky, D.E., Broadfoot, A.L., Holberg, J.B., Smith, G.R., 1982. Science 215, 548] do not exhibit methane absorption. The auroral H₂ emission is thus likely located above but close to the methane homopause. The H₂ auroral brightness in the 800-1700 Å bandwidth varies from 2.9 kR to 139 kR, comparable to values derived from FUV Faint Object Camera (FOC) and STIS images.