Mass and interior of Enceladus from Cassini data analysis

Gravity results are available from radio Doppler data acquired by the Deep Space Network during the encounter of the Cassini spacecraft with Enceladus in February 2005. We report the mass of Enceladus to be (1.0798±0.0016)×¹⁰²⁰ kg, which implies a density of . For a core made of hydrated silicates with a density of 2500 kg m⁻³ the core radius is ∼190 km and the quadrupole moment C₂₂∼1.4×¹⁰⁻³. If Enceladus is in hydrostatic equilibrium, the larger than previously anticipated density implies that the recently proposed secondary spin-orbit resonance cannot be present. Therefore, the source of endogenic activity of Enceladus remains unexplained.     

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.     

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.     

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.     

Observations of C2H6 and C2H2 in the stratosphere of Saturn

We have performed high-resolution spectral observations at mid-infrared wavelengths of C₂H₆ (12.16 μm), and C₂H₂ (13.45 μm) on Saturn. These emission features probe the stratosphere of the planet and provide information on the hydrocarbon photochemical processes taking place in that region of the atmosphere. The observations were performed using our cryogenic echelle spectrometer Celeste, in conjunction with the McMath-Pierce 1.5-m solar telescope in November and December 1994. We used Voyager IRIS CH₄ observations (7.67 μm) to derive a temperature profile on the saturnian atmosphere for the region of the stratosphere. This profile was then used in conjunction with height-dependent volume mixing ratios of each hydrocarbon to determine global abundances for ethane and acetylene. Our ground-based measurements indicate abundances of for C₂H₆ (1.0 mbar pressure level), and for C₂H₂ (1.6 mbar pressure level). We also derived new mixing ratios from the Voyager mid-latitude IRIS observations; 8.6±0.9×¹⁰⁻⁶ for C₂H₆ (0.1-3.0 mbar pressure level), and 1.6±0.2×¹⁰⁻⁷ for C₂H₂ (2.0 mbar pressure level).     

Meridional variations of temperature, C2H2 and C2H6 abundances in Saturn's stratosphere at southern summer solstice

Measurements of the vertical and latitudinal variations of temperature and C₂H₂ and C₂H₆ abundances in the stratosphere of Saturn can be used as stringent constraints on seasonal climate models, photochemical models, and dynamics. The summertime photochemical loss timescale for C₂H₆ in Saturn's middle and lower stratosphere (∼40-10,000 years, depending on altitude and latitude) is much greater than the atmospheric transport timescale; ethane observations may therefore be used to trace stratospheric dynamics. The shorter chemical lifetime for C₂H₂ (∼1-7 years depending on altitude and latitude) makes the acetylene abundance less sensitive to transport effects and more sensitive to insolation and seasonal effects. To obtain information on the temperature and hydrocarbon abundance distributions in Saturn's stratosphere, high-resolution spectral observations were obtained on September 13-14, 2002 UT at NASA's IRTF using the mid-infrared TEXES grating spectrograph. At the time of the observations, Saturn was at a LS≈270°, corresponding to Saturn's southern summer solstice. The observed spectra exhibit a strong increase in the strength of methane emission at 1230 cm⁻¹ with increasing southern latitude. Line-by-line radiative transfer calculations indicate that a temperature increase in the stratosphere of ≈10 K from the equator to the south pole between 10 and 0.01 mbar is implied. Similar observations of acetylene and ethane were also recorded. We find the 1.16 mbar mixing ratio of C₂H₂ at −1° and −83° planetocentric latitude to be and , respectively. The C₂H₂ mixing ratio at 0.12 mbar is found to be at −1° planetocentric latitude and at −83° planetocentric latitude. The 2.3 mbar mixing ratio of C₂H₆ inferred from the data is and at −1° and −83° planetocentric latitude, respectively. Further observations, creating a time baseline, will be required to completely resolve the question of how much the latitudinal variations of C₂H₂ and C₂H₆ are affected by seasonal forcing and/or stratospheric circulation.     

Gravitational tidal waves in Titan's upper atmosphere

Tidal waves driven by Titan's orbital eccentricity through the time-dependent component of Saturn's gravitational potential attain nonlinear, saturation amplitudes (|T^′|>10 K, , and ) in the upper atmosphere (?500 km) due to the approximate exponential growth as the inverse square root of pressure. The gravitational tides, with vertical wavelengths of ∼100-150 km above 500 km altitude, carry energy fluxes sufficient in magnitude to affect the energy balance of the upper atmosphere with heating rates in the altitude range of 500-900 km.     

(42355) Typhon-Echidna: Scheduling observations for binary orbit determination

We describe a strategy for scheduling astrometric observations to minimize the number required to determine the mutual orbits of binary transneptunian systems. The method is illustrated by application to Hubble Space Telescope observations of (42355) Typhon-Echidna, revealing that Typhon and Echidna orbit one another with a period of 18.971±0.006 days and a semimajor axis of 1628±29 km, implying a system mass of (9.49±0.52)×¹⁰¹⁷ kg. The eccentricity of the orbit is 0.526±0.015. Combined with a radiometric size determined from Spitzer Space Telescope data and the assumption that Typhon and Echidna both have the same albedo, we estimate that their radii are and , respectively. These numbers give an average bulk density of only , consistent with very low bulk densities recently reported for two other small transneptunian binaries.     

Estimating the masses of Saturn’s A and B rings from high-optical depth N-body simulations and stellar occultations

We have completed a series of local N-body simulations of Saturn’s B and A rings in order to identify systematic differences in the degree of particle clumping into self-gravity wakes as a function of orbital distance from Saturn and dynamical optical depth (a function of surface density). These simulations revealed that the normal optical depth of the final configuration can be substantially lower than one would infer from a uniform distribution of particles. Adding more particles to the simulation simply piles more particles onto the self-gravity wakes while leaving relatively clear gaps between the wakes. Estimating the mass from the observed optical depth is therefore a non-linear problem. These simulations may explain why the Cassini UVIS instrument has detected starlight at low incidence angles through regions of the B ring that have average normal optical depths substantially greater than unity at some observation geometries [Colwell, J.E., Esposito, L.W., Srem?evi?, M., Stewart, G.R., McClintock, W.E., 2007. Icarus 190, 127-144]. We provide a plausible internal density of the particles in the A and B rings based upon fitting the results of our simulations with Cassini UVIS stellar occultation data. We simulated Cassini-like occultations through our simulation cells, calculated optical depths, and attempted to extrapolate to the values that Cassini observes. We needed to extrapolate because even initial optical depths of >4 (σ > 240 g cm⁻²) only yielded final optical depths no greater than 2.8, smaller than the largest measured B ring optical depths. This extrapolation introduces a significant amount of uncertainty, and we chose to be conservative in our overall mass estimates. From our simulations, we infer the surface density of the A ring to be , which corresponds to a mass of . We infer a minimum surface density of for Saturn’s B ring, which corresponds to a minimum mass estimate of . The A ring mass estimate agrees well with previous analyses, while the B ring is at least 40% larger. In sum, our lower limit estimate is that the total mass of Saturn’s ring system is 120-200% the mass of the moon Mimas, but significantly larger values would be plausible given the limitations of our simulations. A significantly larger mass for Saturn’s rings favors a primordial origin for the rings because the disruption of a former satellite of the required mass would be unlikely after the decay of the late heavy bombardment of planetary surfaces.     

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.     

Saturn A ring surface mass densities from spiral density wave dispersion behavior

We have undertaken an analysis of the Voyager photopolarimeter (PPS) stellar occultation data of Saturn's A ring. The Voyager PPS observed the bright star δ Scorpii as it was occulted by Saturn's main rings during the spacecraft flyby of the Saturn system in 1981. The occultation measurement produced a ring profile with radial resolution of approximately 100 m, and radial structure is evident in the profile down to the resolution limit. We have applied an autoregressive technique to the data for estimating the power spectrum as a function of radius, which has allowed us to identify 40 spiral density waves in Saturn's A ring, associated with the strongest torques due to forcing from the moons. The majority of the detected waves are observed to disperse linearly in regions beginning a few kilometers from the resonance location. We have used the dispersion behavior for those waves to calculate local surface mass densities in the vicinity of each wave. We find that the inner three-quarters of the A ring (up to the beginning of the Encke gap) has an average surface mass density of , while the outer region has an average surface mass density of . The two regions have different mean surface mass densities with a significance of approximately 0.999993, as estimated with a T-statistic, which corresponds to about 4.5σ. While the mean optical depth of the A ring increases slightly with increasing distance from Saturn, we find that it is not significantly correlated with the surface mass density; the two quantities having a linear Pearson's correlation coefficient of r_corr≈−0.03. The variation of mass density, independent of PPS optical depth, is consistent with previous conjectures that the particle size distribution and composition are not constant across the entire A ring, particularly in the very outer portion. We estimate the mass of Saturn's A ring from our surface mass density estimates as 4.9×10²¹ gm, or 8.61×10⁻⁹ of the mass of Saturn, roughly equivalent to the mass of a 110-km diameter icy satellite. This mass is almost 25% smaller than estimates from previous studies, but is well within the expected errors of the derived mass densities. We also identified three previously unstudied features which exhibit linear dispersion. The strongest of these features is tentatively identified as the Janus 13:11 density wave. The other two features do not fall near any known satellite resonances and may represent density waves created by previously undetected satellites.     

CCD photometry of 23 minor planets

We present CCD photometric observations of 23 main-belt asteroids, of which 8 have never been observed before; thus, the data of these objects are the first in the literature. The majority showed well-detectable light variations, exceeding 0^m1. We have determined synodic periods for 756 Lilliana (936), 1270 Datura (34), 1400 Tirela (1336), 1503 Kuopio (998), 3682 Welther (359), 7505 Furushu (414) and 11436 1969 QR (123), while uncertain period estimates were possible for 469 Argentina (123), 546 Herodias (104) and 1026 Ingrid (53). The shape of the lightcurves of 3682 Welther changed on a short time-scale and showed dimmings that might be attributed to eclipses in a binary system. For the remaining objects, only lower limits of the periods and amplitudes were concluded.     

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.     

Estimating Mercury's 88-day libration amplitude from orbit

In this paper it is derived that the libration of Mercury can be described by where Φ₀ is the unknown libration amplitude, M is Mercury's mean anomaly and K=−9.483. Φ₀ can be determined by comparing pairs of images of the same landmarks taken by an orbiter at different positions of Mercury. If the angle between the orbit plane of a polar orbiter and Mercury's line of periapsis is between −60° and 60° and if one landmark at the equator is imaged per day with a relative precision of , then the libration amplitude can be determined in two Mercury years (176 days) with an accuracy of or better, which is sufficient to answer the question whether Mercury has a solid or fluid core.     

Production, ionization and redistribution of O2 in Saturn's ring atmosphere

Molecular oxygen produced by the decomposition of icy surfaces is ubiquitous in Saturn's magnetosphere. A model is described for the toroidal O₂ atmosphere indicated by the detection of and O⁺ over the main rings. The O₂ ring atmosphere is produced primarily by UV photon-induced decomposition of ice on the sunlit side of the ring. Because O₂ has a long lifetime and interacts frequently with the ring particles, equivalent columns of O₂ exist above and below the ring plane with the scale height determined by the local ring temperature. Energetic particles also decompose ice, but estimates of their contribution over the main rings appear to be very low. In steady state, the O₂ column density over the rings also depends on the relative efficiency of hydrogen to oxygen loss from the ring/atmosphere system with oxygen being recycled on the grain surfaces. Unlike the neutral density, the ion densities can differ on the sunlit and shaded sides due to differences in the ionization rate, the quenching of ions by the interaction with the ring particles, and the northward shift of the magnetic equator relative to the ring plane. Although O⁺ is produced with a significant excess energy, is not. Therefore, should mirror well below those altitudes at which ions were detected. However, scattering by ion-molecule collisions results in much larger mirror altitudes, in ion temperatures that go through a minimum over the B-ring, and in the redistribution of both molecular hydrogen and oxygen throughout the magnetosphere. The proposed model is used to describe the measured oxygen ion densities in Saturn's toroidal ring atmosphere and its hydrogen content. The oxygen ion densities over the B-ring appear to require either significant levels of UV light scattering or ion transmission through the ring plane.     

The orbit, mass, size, albedo, and density of (65489) Ceto/Phorcys: A tidally-evolved binary Centaur

Hubble Space Telescope observations of Uranus- and Neptune-crossing object (65489) Ceto/Phorcys (provisionally designated 2003 FX₁₂₈) reveal it to be a close binary system. The mutual orbit has a period of 9.554±0.011 days and a semimajor axis of 1840±48 km. These values enable computation of a system mass of (5.41±0.42)×¹⁰¹⁸ kg. Spitzer Space Telescope observations of thermal emission at 24 and 70 μm are combined with visible photometry to constrain the system's effective radius and geometric albedo . We estimate the average bulk density to be , consistent with ice plus rocky and/or carbonaceous materials. This density contrasts with lower densities recently measured with the same technique for three other comparably-sized outer Solar System binaries (617) Patroclus, (26308) 1998 SM₁₆₅, and (47171) 1999 TC₃₆, and is closer to the density of the saturnian irregular satellite Phoebe. The mutual orbit of Ceto and Phorcys is nearly circular, with an eccentricity ?0.015. This observation is consistent with calculations suggesting that the system should tidally evolve on a timescale shorter than the age of the Solar System.